孫雨薇，李艷軍*，劉永昌，薛 飛，朱守鴻，王翔飛，孫 杰

(1 石河子大學(xué) 農(nóng)學(xué)院，新疆生產(chǎn)建設(shè)兵團(tuán)綠洲生態(tài)農(nóng)業(yè)重點(diǎn)實(shí)驗(yàn)室，新疆石河子832003；2 石河子大學(xué) 藥學(xué)院，新疆石河子 832003)

?

陸地棉PIP亞家族基因的克隆及表達(dá)特征分析

孫雨薇1，李艷軍1*，劉永昌1，薛 飛1，朱守鴻1，王翔飛2，孫 杰1

(1 石河子大學(xué) 農(nóng)學(xué)院，新疆生產(chǎn)建設(shè)兵團(tuán)綠洲生態(tài)農(nóng)業(yè)重點(diǎn)實(shí)驗(yàn)室，新疆石河子832003；2 石河子大學(xué) 藥學(xué)院，新疆石河子 832003)

該研究以陸地棉苯基香豆?jié)M芐基醚還原酶( phenylcoumaran benzylic ether reductase，PCBER)氨基酸序列為探針，利用Blastp從陸地棉基因組數(shù)據(jù)庫(kù)中發(fā)現(xiàn)了6個(gè)同源性較高的基因。根據(jù)6個(gè)基因序列設(shè)計(jì)引物，利用 RT-PCR技術(shù)從陸地棉纖維細(xì)胞中克隆出了這6個(gè)基因的全長(zhǎng)cDNA序列，分別命名為GhPCBER1、GhPCBER2、GhPCBER3、GhIFR、GhPLR1和GhPLR2。多重序列比對(duì)和進(jìn)化樹分析發(fā)現(xiàn)，6個(gè)蛋白均含有PIP類型蛋白的所有保守性基序和活性殘基，屬于PIP亞家族。實(shí)時(shí)熒光定量PCR結(jié)果顯示，除GhPLR1之外其他5個(gè)PIP亞家族基因均在纖維細(xì)胞中優(yōu)勢(shì)或特異表達(dá)；在纖維發(fā)育過(guò)程中，GhPCBER1、GhPCBER2、GhPCBER3和GhIFR的表達(dá)均表現(xiàn)為先上升后下降，GhPCBER1和GhPCBER2在花后21 d表達(dá)量達(dá)到最高，GhPCBER3和GhIFR在花后18 d達(dá)到最高，GhPLR1和GhPLR2在纖維中的表達(dá)量呈持續(xù)上升趨勢(shì)。根據(jù)基因的表達(dá)特征，推測(cè)PIP亞家族可能在棉纖維的發(fā)育過(guò)程中發(fā)揮著重要作用。

棉花；纖維；PIP；表達(dá)分析

棉花是全球性的重要經(jīng)濟(jì)作物，是紡織工業(yè)最主要的天然纖維來(lái)源。成熟棉纖維含有90%以上的纖維素和少量的非纖維素碳水化合物，如木糖葡聚糖、葡聚糖和果膠多糖等[1-2]。纖維素的合成是棉纖維發(fā)育過(guò)程中主要的代謝途徑，近期的研究表明棉纖維中存在苯丙烷代謝途徑且在纖維次生壁增厚期優(yōu)勢(shì)表達(dá)[3-5]，是纖維發(fā)育過(guò)程中僅次于纖維素代謝的第二大代謝途徑[4]，該途徑產(chǎn)生的次生代謝產(chǎn)物可能對(duì)棉纖維的發(fā)育與品質(zhì)形成具有重要影響[5-6]。苯丙烷代謝途徑是植物最重要的次生代謝途徑之一[7]，能夠產(chǎn)生木質(zhì)素、木脂素、黃酮、異黃酮、花青素、酚類化合物等次生代謝產(chǎn)物[8-9]。木質(zhì)素和木脂素是兩種不同的物質(zhì)，木質(zhì)素是植物體內(nèi)一種大分子有機(jī)物，是具有三維結(jié)構(gòu)的芳香族高聚物，與纖維素、半纖維素合稱為“三素”是形成植物骨架的主要成分。木脂素是由苯丙素雙分子聚合而成的天然化合物，是一種植物雌激素。在苯丙烷代謝過(guò)程中肉桂醛經(jīng)肉桂醇脫氫酶(cinnamyl alcohol dehydrogenase，CAD)還原為香豆醇(p-Coumaryl，alcohol)、松柏醇(coniferyl alcohol)和芥子醇(sinapyl alcohol)這3種肉桂醇即是合成木質(zhì)素的主要單體。其中松柏醇能進(jìn)一步形成木脂素。棉纖維中苯丙烷類物質(zhì)的發(fā)現(xiàn)為棉纖維的品質(zhì)改良提供了新的思路。利用基因工程手段從棉纖維中克隆苯丙烷物質(zhì)生物合成相關(guān)基因并研究其功能，有助于弄清棉纖維細(xì)胞中這類物質(zhì)的存在情況及其在棉纖維發(fā)育及纖維品質(zhì)形成中的作用，為利用基因工程手段改良棉纖維的品質(zhì)奠定理論基礎(chǔ)。

PIP亞家族是同源性較高的SDR(short-chain dehydrogenase/reductase)家族成員[10]，是木脂素和異黃酮等植物次生代謝產(chǎn)物的關(guān)鍵酶。包括苯基香豆?jié)M芐基醚還原酶(phenylcoumaran benzylic ether reductase，PCBER)，異黃酮還原酶(isoflavone reductase，IFR)和松脂醇-落葉松脂醇還原酶(pinoresinol-lariciresinol reductases，PLR)[11]，3種蛋白不僅有很高的氨基酸序列相似性并且都具有一個(gè)連續(xù)的α/β NADPH結(jié)合結(jié)構(gòu)域和較小的底物結(jié)合域[12-13]。PCBER是苯丙烷代謝產(chǎn)物木脂素合成過(guò)程中的關(guān)鍵酶，在棉花的子葉、莖、葉片、根和發(fā)育的纖維中均發(fā)現(xiàn)了PCBER的轉(zhuǎn)錄產(chǎn)物[14-17]，研究發(fā)現(xiàn)陸地棉種皮中有高水平的PCBER轉(zhuǎn)錄[12]，棉花無(wú)絮突變體和野生型發(fā)育起始期的棉纖維比較蛋白質(zhì)組研究發(fā)現(xiàn)PCBER相應(yīng)的蛋白點(diǎn)僅存在于野生型棉纖維中[12,18]。PLR是松脂醇-落葉松酯醇還原酶，在8-8’連接結(jié)構(gòu)的木脂素合成途徑中，它可以催化松脂醇和落葉松脂醇還原為開環(huán)異落葉松樹脂酚。近期的研究中發(fā)現(xiàn)，擬南芥中的PCBER具有弱PLR酶活性，可催化相同的底物松脂醇和落葉松脂醇還原為開環(huán)異落葉松樹脂酚[10]。IFR是異黃酮類化合物生物合成途徑中的重要酶，參與抵御各種生物和非生物脅迫，銀杏中的IFR具有PCBER酶活性[19-20]。目前已經(jīng)從多種維管植物種中分離出PCBER、PLR和IFR，PIP亞家族基因轉(zhuǎn)化植物后，對(duì)于植物中多種代謝產(chǎn)物的含量產(chǎn)生影響[19-23]。目前棉花中尚未見PIP亞家族基因的報(bào)道，本研究利用陸地棉PCBER蛋白序列為探針序列同源比對(duì)出多個(gè)相似性序列，其中包括PCBER、IFR和PLR，為進(jìn)一步研究PIP亞家族基因在棉纖維發(fā)育及品質(zhì)形成中的功能奠定了基礎(chǔ)。

1 材料和方法

1.1 材 料

植物材料為新疆陸地棉(GossypiumhirsutumL.)品種‘新陸早33號(hào)’,由石河子大學(xué)棉花研究所提供。將‘新陸早33號(hào)’播種于石河子大學(xué)實(shí)驗(yàn)場(chǎng)。在棉田盛花期對(duì)當(dāng)日花掛牌標(biāo)記，以3 d為一時(shí)間節(jié)點(diǎn)，摘取3～24 DPA(day post anthesis，開花后天數(shù))棉鈴，室內(nèi)剝?nèi)∨咧榧袄w維；將‘新陸早33號(hào)’種子用濃硫酸脫絨后種于花盆中(營(yíng)養(yǎng)土︰蛭石=2︰1)，置于人工氣候室，待長(zhǎng)至2片真葉時(shí)將其小心拔出，水培3 d，取其根、莖、葉，從田間摘取當(dāng)日花。將上述植物材料液氮速凍后保存于-80℃，用于提取RNA。

1.2 方 法

1.2.1 不同組織RNA提取及cDNA第一鏈合成 利用CTAB/酸酚法提取棉花組織總RNA[24]。使用ND1000紫外分光光度計(jì)定量，根據(jù)A260/280和A260/230比值對(duì)樣品純度進(jìn)行鑒定，用1.2%瓊脂糖凝膠電泳檢測(cè)所提RNA的完整性。對(duì)質(zhì)量合格的RNA按照Reverse Transcriptase M-MLV(RNase H-)說(shuō)明書完成cDNA第一鏈的合成。

表1 PIP亞家族基因ORF擴(kuò)增及qPCR引物序列

1.2.2 棉花PIP亞家族基因的克隆 以陸地棉PCBER蛋白序列(GenBank登錄號(hào)ABN12322.1)作為查找序列，在陸地棉的基因組數(shù)據(jù)庫(kù)(https://www.cottongen.org)中進(jìn)行Blastp比對(duì)，獲得6個(gè)具有較高同源性的基因，利用Blastp比對(duì)和進(jìn)化樹分析，發(fā)現(xiàn)6個(gè)基因包含3個(gè)PCBER，1個(gè)IFR和2個(gè)PLR基因，分別對(duì)6個(gè)基因進(jìn)行命名(表1)。每個(gè)基因在陸地棉基因組中均具有2個(gè)拷貝，分別位于A染色體組和D染色體組上，一一對(duì)應(yīng)(表1)。設(shè)計(jì)擴(kuò)增6個(gè)基因ORF的引物(表1)，以棉纖維cDNA為模板對(duì)6個(gè)基因進(jìn)行PCR擴(kuò)增，反應(yīng)體系為：10×buffer 2 μL、2.5 mmol/L dNTPs 1.6 μL、5 U/ μLEX-taqDNA 聚合酶0.25 μL、10 μmol/L上下游引物各1 μL，補(bǔ) ddH2O至20 μL。反應(yīng)程序?yàn)椋?4 ℃預(yù)變性3 min；94 ℃變性30 s, 52～57 ℃退火30 s，72 ℃延伸1 min；重復(fù)31個(gè)循環(huán)；72 ℃延伸10 min。將PCR擴(kuò)增產(chǎn)物回收后連接至pGEM-T Easy載體，轉(zhuǎn)化大腸桿菌XL1-Blue感受態(tài)細(xì)胞，挑取單菌落，提取質(zhì)粒進(jìn)行PCR和酶切鑒定，正確的質(zhì)粒送至生工生物工程(上海)股份有限公司測(cè)序。

1.2.3 棉花PIP亞家族基因生物信息學(xué)分析 利用Expasy中的Protparam軟件進(jìn)行氨基酸基本理化性質(zhì)分析；用SOPMA軟件對(duì)蛋白序列進(jìn)行二級(jí)結(jié)構(gòu)預(yù)測(cè)；用TMHMM2.0在線工具對(duì)氨基酸序列跨膜結(jié)構(gòu)進(jìn)行預(yù)測(cè)分析，用SignalP 4.1和TargetP 1.1 Server軟件進(jìn)行信號(hào)肽預(yù)測(cè)和亞細(xì)胞定位分析；氨基酸序列的同源比對(duì)及多序列比對(duì)采用Clustal W和DNAMAN進(jìn)行分析。利用MAGE 6.06軟件采用鄰近相接法(neighbor-joining method)構(gòu)建系統(tǒng)進(jìn)化樹，Bootstrap參數(shù)設(shè)為重復(fù)檢測(cè)1 000次。

1.2.4 實(shí)時(shí)熒光定量PCR分析 利用qPCR技術(shù)分析棉花PIP亞家族基因在棉花不同組織及不同發(fā)育時(shí)期纖維中的時(shí)空表達(dá)模式，根據(jù)6個(gè)基因的cDNA序列，分別設(shè)計(jì)特異性qPCR引物(表1)，以陸地棉His3為內(nèi)參基因，以棉花幼苗的根、莖、葉、花和不同發(fā)育時(shí)期(花后3、6、9、12、15、18、21和24 d)的棉纖維cDNA為模板進(jìn)行qPCR擴(kuò)增。反應(yīng)程序?yàn)椋?4 ℃預(yù)變性1 min；95 ℃變性15 s，56～57 ℃ 退火20 s，72 ℃延伸30 s，45個(gè)循環(huán)。每個(gè)樣設(shè)置3個(gè)實(shí)驗(yàn)重復(fù)，將數(shù)據(jù)輸入Excel中按照2-ΔΔCt法分析。

2 結(jié)果與分析

2.1 PIP亞家族基因全長(zhǎng)cDNA序列克隆

以陸地棉PCBER蛋白序列作為查詢序列，在陸地棉基因組數(shù)據(jù)庫(kù)中進(jìn)行Blastp比對(duì)，發(fā)現(xiàn)了6個(gè)同源基因。根據(jù)Blastp比對(duì)結(jié)果，將6個(gè)基因分別命名為GhPCBER1、GhPCBER2、GhPCBER3、GhIFR、GhPLR1和GhPLR2。陸地棉基因組中6個(gè)PIP亞家族基因的基因組序列全長(zhǎng)為3 756～4 779 bp，CDS(coding sequence)長(zhǎng)度為921～942 bp，編碼306～313個(gè)氨基酸(表2)。根據(jù)基因的序列設(shè)計(jì)擴(kuò)增ORF的引物，以陸地棉纖維cDNA為模板對(duì)6個(gè)基因進(jìn)行PCR擴(kuò)增，結(jié)果如圖1所示，擴(kuò)增條帶與預(yù)期結(jié)果相一致。

2.2 PIP亞家族基因編碼蛋白的生物信息學(xué)分析

根據(jù)6個(gè)基因的測(cè)序結(jié)果，運(yùn)用ProtParam軟件預(yù)測(cè)相應(yīng)蛋白的基本理化性質(zhì)(表3)，6個(gè)蛋白分子量為33 854.9～35 199.6 Da，理論等電點(diǎn)為5.52～7.77，負(fù)電荷的氨基酸殘基數(shù)(Asp+Glu)是35～41，正電荷的氨基酸殘基數(shù)(Arg+Lys)是32～36；脂肪系數(shù)為91.75～101.96，氨基酸殘基總的親疏水性平均系數(shù)均為負(fù)值，分布在-0.131～ -0.056之間，說(shuō)明蛋白的親水性較強(qiáng)。不穩(wěn)定系數(shù)是22.55～43.71，其中GhPLR1和GhPLR2不穩(wěn)定系數(shù)大于40，為不穩(wěn)定蛋白，其余為穩(wěn)定蛋白。利用SOPMA在線預(yù)測(cè)分析蛋白的二級(jí)結(jié)構(gòu)(表4)，結(jié)果表明該蛋白主要由α螺旋、β轉(zhuǎn)角、延伸直鏈和無(wú)規(guī)則卷曲4種結(jié)構(gòu)形式構(gòu)成。其中，GhPCBER1和GhPCBER2蛋白二級(jí)結(jié)構(gòu)主要由無(wú)規(guī)則卷曲組成，GhPCBER3、GhPLR2主要由α螺旋組成，GhIFR和GhPLR1以α螺旋和無(wú)規(guī)則卷曲為主要組成結(jié)構(gòu)。TMHMM2.0預(yù)測(cè)表明6個(gè)蛋白均未形成跨膜結(jié)構(gòu)域。用SignalP 4.1和TargetP 1.1軟件對(duì)氨基酸序列跨膜預(yù)測(cè)結(jié)果顯示PIP亞家族的6個(gè)蛋白均無(wú)信號(hào)肽，是非分泌蛋白。對(duì)6個(gè)蛋白的信號(hào)肽及亞細(xì)胞定位預(yù)測(cè)顯示PIP亞家族蛋白均定位于細(xì)胞質(zhì)中。

表2 PIP亞家族基因在陸地棉基因組數(shù)據(jù)庫(kù)中的序列信息

2.3 PIP蛋白多重序列比對(duì)和蛋白進(jìn)化樹分析

為了進(jìn)一步明確棉花PIP亞家族6個(gè)基因的典型結(jié)構(gòu)域及與其他植物PIP亞家族基因的同源性，選取已報(bào)道的部分植物PIP亞家族基因，利用Clustal X和DNAMAN軟件對(duì)不同植物中PIP亞家族基因進(jìn)行多重序列比對(duì)，發(fā)現(xiàn)均有與NADPH結(jié)合的“GXXGXXG”保守結(jié)合域(圖2)。利用MEGA6.06軟件對(duì)不同植物PIP亞家族序列進(jìn)行進(jìn)化樹分析，結(jié)果顯示所有基因聚為三大類，大致可分為PCBER、PLR和IFR，而PCBER與IFR的親緣關(guān)系比PLR近。GhPCBER1、GhPCBER2和GhPCBER3與連翹和火炬松的PCBER聚為一類，GhPLR1、GhPLR2與北美紅崖柏和亞麻的PLR聚為一類，GhIFR與楊樹和琴葉擬南芥的IFR聚為一類(圖3)。

M. DNA marker Ⅲ; 1～6. GhPCBER1、GhPCBER2、GhPCBER3、GhIFR、GhPLR1和GhPLR2的 PCR擴(kuò)增結(jié)果圖1 陸地棉PIP亞家族基因 PCR 擴(kuò)增電泳M. DNA marker Ⅲ; 1-6 represent the PCR results of GhPCBER1, GhPCBER2, GhPCBER3, GhIFR, GhPLR1 and GhPLR2, respectivelyFig.1 PCR results of PIP subfamily genes

基因名稱Nameα螺旋Alphahelixβ轉(zhuǎn)角Betaturn延伸直鏈Extendedstrand無(wú)規(guī)則卷曲RandomcoilGhPCBER127.929.4225.0037.66GhPCBER231.1710.3923.7034.74GhPCBER338.249.8021.5730.39GhIFR32.3511.1123.8632.68GhPLR130.4511.8625.9631.73GhPLR235.4610.5426.2027.80

表3 蛋白的基本理化性質(zhì)

黑色方框?yàn)镻IP亞家族與NADPH結(jié)合的GXXGXXG保守結(jié)合域；Fi1PCBER.連翹(AF242491)；GhPCBER1、GhPCBER2、GhPCBER3、GhIFR、GhPLR1、GhPLR2.陸地棉(Gh_D01G0977、Gh_A03G0704、Gh_D12G0604、Gh_D02G2151、Gh_A08G1368、Gh_D08G1661)；PtIFR.楊樹(XP_002313789.1)；LuPLR.亞麻(CAH60858)圖2 陸地棉PIP與其它PIP蛋白的多重序列比對(duì)The black box is a conserved NADPH binding domain of GXXGXXG which is combined with PIP subfamily；Fi1PCBER. Forsythia suspensa(AF242491); GhPCBER1，GhPCBER2，GhPCBER3，GhIFR，GhPLR1，GhPLR2. G. hirsutum L.(Gh_D01G0977，Gh_A03G0704，Gh_D12G0604，Gh_D02G2151，Gh_A08G1368，Gh_D08G1661)；PtIFR. Populus trichocarpa (XP_002313789.1); LuPLR. Linum usitatissimum (CAH60858)Fig.2 Multiple sequence alignment of PIP and other PIP proteins

系統(tǒng)樹各分支上數(shù)字是Bootstrap 1 000次循環(huán)檢驗(yàn)的置信度；標(biāo)尺代表遺傳距離圖3 陸地棉PIP與其他 PIP蛋白的進(jìn)化樹分析The numbers on the tree branches represent bootstrap confidence values as bootstrap is 1 000；The scale bar represents genetic distanceFig.3 Phylogenic relationship of PIP and other PIP proteins

R.根；S.莖；L.葉；F.花；3～24. 花后3、6、9、12、15、18、21、24 d的纖維細(xì)胞圖4 陸地棉PIP亞家族基因的實(shí)時(shí)熒光定量 PCR 分析R. Root; S. Stem; L. Leaf; F. Flower; 3-24. Fiber cells at 3, 6, 9, 12, 15, 18, 21 and 24 day post anthesis, respectivelyFig.4 Relative expression of PIP subfamily genes

2.4 PIP亞家族基因在棉花不同組織中的表達(dá)特征分析

為研究棉花PIP亞家族6個(gè)基因在陸地棉中的表達(dá)模式，根據(jù)每個(gè)基因的cDNA序列設(shè)計(jì)特異性引物，對(duì)其在棉花不同組織和不同發(fā)育時(shí)期纖維中的表達(dá)模式進(jìn)行實(shí)時(shí)熒光定量分析。結(jié)果(圖4)表明，在6個(gè)PIP亞家族基因中，GhPLR1在根和莖中表達(dá)量明顯較高，在纖維中表達(dá)無(wú)優(yōu)勢(shì)，其它5個(gè)基因均在纖維中特異或優(yōu)勢(shì)表達(dá)。GhPCBER3在纖維中特異表達(dá)，GhPCBER1、GhPCBER2、GhIFR和GhPLR2在纖維中優(yōu)勢(shì)表達(dá)，在根和莖中有少量表達(dá)，在葉和花中幾乎不表達(dá)。GhPLR1和GhPLR2在纖維中的表達(dá)量呈上升趨勢(shì)；GhPCBER1、GhPCBER2和GhIFR在纖維中的表達(dá)量均呈先升高后降低的趨勢(shì)。6個(gè)基因在纖維發(fā)育起始期的表達(dá)量均較低，在次生壁合成期表達(dá)量較高，推測(cè)它們?cè)诶w維次生壁合成期可能具有一定的功能。

3 討 論

前期研究認(rèn)為棉纖維次生壁中不存在木質(zhì)素[25]。Fan等[5]利用熒光觀察、傅里葉轉(zhuǎn)換紅外線光譜(FTER)分析及高效液相色譜(HPLC)分析發(fā)現(xiàn)成熟棉纖維中含有木質(zhì)素，還發(fā)現(xiàn)參與單體木質(zhì)醇生物合成的CAD基因在纖維次生壁合成時(shí)期高度表達(dá)。Han等[6]克隆了WLIM1轉(zhuǎn)錄因子，該基因可以調(diào)控苯丙烷代謝相關(guān)基因的表達(dá)，該基因在棉花中過(guò)量表達(dá)致使棉纖維木質(zhì)素含量增多，長(zhǎng)度變長(zhǎng)，強(qiáng)力增加，推測(cè)木脂素對(duì)棉纖維的品質(zhì)具有影響。植物中的次生代謝產(chǎn)物包括木脂素、黃酮、異黃酮、花青素、酚類化合物等[8-9]，關(guān)于它們?cè)诿蘩w維中的存在情況及與棉纖維品質(zhì)的相關(guān)性未見報(bào)道。PCBER、IFR和PLR是PIP亞家族基因，是合成木脂素和異黃酮過(guò)程中重要的酶基因，本研究對(duì)陸地棉中PIP亞家族基因進(jìn)行克隆與分析，對(duì)于明確棉纖維中次生代謝產(chǎn)物的存在情況及它們對(duì)纖維發(fā)育和品質(zhì)形成的影響具有重要意義。

本研究以陸地棉PCBER蛋白為探針序列，從陸地棉基因組數(shù)據(jù)庫(kù)中發(fā)現(xiàn)了12個(gè)同源基因，根據(jù)cDNA序列相似性劃分為6對(duì)，每對(duì)中的兩個(gè)成員均分別位于A和D染色體組上，且具有較高的序列相似性。由于每對(duì)基因中的2個(gè)成員序列相似性高(大于97%)，難以進(jìn)行區(qū)分，因此將它們視為基因的兩個(gè)拷貝。本研究共設(shè)計(jì)6個(gè)擴(kuò)增ORF的引物和6個(gè)qPCR引物，克隆了6個(gè)基因。

由生物信息學(xué)分析發(fā)現(xiàn)，6個(gè)蛋白的三級(jí)結(jié)構(gòu)均具有PIP的典型結(jié)構(gòu)域，即靠近N端區(qū)域都存在一段特異性的保守序列“GXXGXXG”即NADPH結(jié)合位點(diǎn)，說(shuō)明6個(gè)基因?qū)儆赑IP亞家族基因。qPCR分析發(fā)現(xiàn)，6個(gè)PIP亞家族基因均在纖維發(fā)育次生壁合成期表達(dá)量較高，暗示PIP亞家族在纖維次生壁合成期具有一定的功能。6個(gè)基因中，其中GhPCBER3在棉纖維中呈現(xiàn)特異表達(dá)，GhPCBER1、GhPCBER2、GhIFR和GhPLR2呈優(yōu)勢(shì)表達(dá)，PIP亞家族基因編碼植物次生代謝產(chǎn)物合成過(guò)程中的關(guān)鍵酶，暗示次生代謝產(chǎn)物在棉纖維的發(fā)育過(guò)程中發(fā)揮著關(guān)鍵作用。

[1] MEINERT M C, DELMER D P. Changes in biochemical composition of the cell wall of the cotton fiber during development[J].PlantPhysiology, 1977, 59(6):1 088-1 097.

[2] HUWYLER H R, FRANZ G, MEIER H. Changes in the composition of cotton fiber cell walls during development[J].Planta, 1979, 146: 635-642.

[3] GOU J Y, WANG L J, CHEN S P,etal. Gene expression and metabolite profiles of cotton fiber during cell elongation and secondary cell wall synthesis[J].CellResearch, 2007, 17: 422-434.

[4] ALGHAZI Y, BOUROT S, ARIOLI T,etal. Transcript profiling during fiber development identifies pathways in secondary metabolism and cell wall structure that may contribute to cotton fiber quality[J].PlantCellPhysiology, 2009, 50: 1 364-1 381.

[5] FAN L, SHI W J, HU W R,etal. Molecular and biochemical evidence for phenylpropanoid synthesis and presence of wall-linked phenolics in cotton fibers[J].JournalofIntegrativePlantBiology, 2009, 51: 626-637.

[6] HAN L B, LI Y B, WANG H Y,etal. The dual functions of WLIM1 a in cell elongation and secondary wall formation in developing cotton fibers[J].ThePlantCell, 2013, 25: 4 421-4 438.

[7] BOERJAN W, RALPH J, BAUCHER M. Lignin biosynthesis[J].AnnualReviewPlantBiology, 2003, 54: 519-546.

[8] GRAND C, BOUDET A, BOUDET A M. Isoenzymes of hydroxycinnamate：CoA ligase from poplar stems, properties and tissue distribution[J].Planta, 1983, 158: 225-229.

[9] WALLIS P J, RHODES M J. Mutiple forms of hydroxycinnamate: CoA ligase in etiolated pea seedlings[J].Phytochemistry, 1977, 16: 1 891-1 894.

[10] GANG D R, DINKOVA-KOSTOVA A T, DAVIN L B,etal. Phylogenetic links in plant defense systems: lignans, isoflavonoids, and their reductases[C]// HEDIN P A, HOLLINGWORTH R M, MASTER E P,etal. Phytochemicals for Pest Control. Washington DC: ACS Symposium Series, 1997: 58-89.

[12] TURLEY R. Expression of a phenylcoumaran benzylic ether reductase-like protein in the ovules ofGossypiumhirsutum[J].BiologiaPlantarum, 2008, 52(4): 759-762.

[13] MIN T, KASAHARA H, BEDGAR D L,etal. Crystal structures of pinoresinol-lariciresinol and phenylcoumaran benzylic ether reductases and their relationship to isoflavone reductases[J].JournalofBiologicalChemistry, 2003, 278(50): 50 714-50 723.

[14] DOWD C, WILSON I, MCFADDEN H. Gene expression profile changes in cotton root and hypocotyl tissues in response to infection withFusariumoxysporumf. sp.vasinfectum[J].MolecularPlant-MicrobeInteracttions, 2004, 17: 654-667.

[15] HAIGLER CH, ZHANG D, WILKERSON C G. Biotechnological improvement of cotton fibre maturity[J].PhysiologiaPlantarum, 2005, 124: 285-294.

[16] TALIERCIO E, ALLEN R D, ESSENBERG M,etal. Analysis of ESTs from multipleGossypiumhirsutumtissues and identification of SSRs[J].Genome, 2006, 49: 306-319.

[17] UDALL J A, SWANSON J M, HALLER K,etal. A global assembly of cotton ESTs[J].GenomeResearch, 2006, 16: 441-450.

[18] TURLEY R B, FERGUSON D L. Changes of ovule proteins during early fiber development and a normal and a fiberless line of cotton (GossypiumhirsutumL.)[J].JournalofPlantPhysiology, 1996, 149: 695-702.

[19] CHENG Q, LI N H, DONG L D,etal. Overexpression of soybean isoflavone reductase (GmIFR) enhances resistance toPhytophthorasojaein soybean[J].FrontPlantSci., 2015, 6: 1 024.

[20] CHENG H, LI L L, XU F,etal. Expression patterns of an isoflavone reductase-like gene and its possibleroles in secondary metabolism inGinkgobiloba[J].PlantCellRep, 2013, 32:637-650.

[21] NUOENDAGULA, NAOFUMI K, TETSUYA M,etal. Expression and functional analyses of a putative phenylcoumaran benzylic ether reductase inArabidopsisthaliana[J].PlantCellRep, 2016, 35: 513-526.

[22] NICULAES C, MORREEL K, KIM H,etal. Phenylcoumaran benzylic ether reductase prevents accumulation of compounds formed under oxidative conditions in poplar xylem[J].PlantCell, 2014, 26(9):3 775-3 791.

[23] BANERJEE A, CHATTOPADHYAY S. Effect of over-expression ofLinumusitatissimumpinoresinol lariciresinol reductase. (LuPLR) gene in transgenicPhyllanthusamarus[J].PlantCellTissOrganCult, 2010, 103:315-323.

[24] 蔣建雄,張?zhí)煺?利用CTAB酸酚法提取棉花組織總RNA[J].棉花學(xué)報(bào),2003,15:166-167.

JIANG J X , ZHANG T Z. Extraction of total RNA in cotton tissues with CTAB acidic phenolic method[J].CottonSci., 2003, 15:166-167.

[25] KIM H J, TRIPLETT B A. Cotton fiber growth in plant andinvitro. Models for plant cell elongation and cell wall biogenesis[J].PlantPhysiology, 2001, 127: 1 361-1 366.

(編輯:宋亞珍)

Cloning and Expression Analysis of PIP Subfamily Genes in Gossypium hirsutum

SUN Yuwei1, LI Yanjun1*, LIU Yongchang1, XUE Fei1, ZHU Shouhong1, WANG Xiangfei2, SUN Jie1

(1 College of Agronomy，Shihezi University，Key Oasis Eco-agriculture Laboratory of Production and Construction Group，Shihezi，Xinjiang 832003，China；2 College of Pharmacy，Shehezi University, Shihezi，Xinjiang 832003，China)

Using the protein sequence ofGhPCBERas a probe, six of genes showing high sequence homology were obtained from a genome database ofGossypiumhirsutumL.tetraploid cotton with Blastp search. These sequences were used to design primers, and then six genes were isolated from cotton fiber ofG.hirsutumL. by using RT-PCR technique and designated asGhPCBER1,GhPCBER2,GhPCBER3,GhIFR,GhPLR1 andGhPLR2, respectively. Multiple sequence alignment and phylogenetic tree analysis revealed that all of the 6 proteins belong to PIP subfamily since the conserved motifs and active residues in their sequences. Real time fluorescent quantitative PCR results showed that in addition toGhPLR1, all the other 5 PIP subfamily genes were specifically or preferentially expressed in fiber cells. Of 6 genes, the expression pattern ofGhPCBER1,GhPCBER2,GhPCBER3 andGhIFRshowed increasing firstly and then decreasing trend during cotton fiber developing stage. The expression level ofGhPCBER1 andGhPCBER2 reached a maximum at 21 day post anthesis, andGhIFRandGhPCBER3 reached a maximum at 18 day post anthesis. The expression pattern ofGhPLR1 andGhPLR2 showed an ascendant trend with cotton fiber developing. According to the characteristics of gene expression, it is speculated that the PIP subfamily may play an important role in cotton fiber development.

cotton; fiber; PIP; expression analysis

1000-4025(2016)10-1948-07

10.7606/j.issn.1000-4025.2016.10.1948

2016-06-21;修改稿收到日期:2016-09-07

國(guó)家自然科學(xué)基金(31460360)

孫雨薇(1990-)，女，在讀碩士研究生，主要從事棉花分子育種研究。E-mail：sunyuwei2014@aliyun.com

*通信作者:李艷軍,博士，副教授，主要從事棉花分子育種研究。E-mail：lyj20022002@sina.com.cn

Q785;Q789

A