摘要：在不同濃度（0、20%、30%、40%、50%）PEG處理下考察木薯（Manihot esculenta Crantz）葉片中脯氨酸含量的動態(tài)變化。從木薯中分別克隆了MeP5CS和MeP5CR基因，對其序列進(jìn)行生物信息學(xué)分析，并分析了MeP5CS和MeP5CR基因在20% PEG脅迫處理下各組織中的表達(dá)。結(jié)果表明，脯氨酸含量受滲透脅迫快速誘導(dǎo)，在一定范圍內(nèi)隨PEG濃度升高而明顯升高，表現(xiàn)出較好的相關(guān)性。序列分析表明，MeP5CS基因cDNA全長2 217 bp，編碼738個氨基酸，含有P5CS保守結(jié)構(gòu)域；MeP5CR基因cDNA全長828 bp，編碼275個氨基酸，含有P5CR保守結(jié)構(gòu)域，它們分別與麻風(fēng)樹和蓖麻中的P5CS、P5CR親緣關(guān)系較近。MeP5CS和MeP5CR基因在轉(zhuǎn)錄水平受到滲透脅迫的調(diào)控。

關(guān)鍵詞：木薯（Manihot esculenta Crantz）；MeP5CS；MeP5CR；基因克隆；干旱；表達(dá)分析

中圖分類號：Q78 文獻(xiàn)標(biāo)識碼：A 文章編號：0439-8114（2016）15-4024-05

DOI：10.14088/j.cnki.issn0439-8114.2016.15.057

Abstract： Dynamic changes of proline content were investigated in cassava（Manihot esculenta Crantz） leaves in response to different concentration of PEG（polyethylene glycol） treatments（0，20%，30%，40% and 50%）. MeP5CS and MeP5CR were cloned in cassava and their sequences were analysed by bioinformatics method. The expression of MeP5CS and MeP5CR genes in different organs was analysed under 20% PEG treatment. The results showed that proline content were rapidly induced by osmotic stress，and at a certain extent，the increase of proline content was linearly correlated with PEG concentrations. Sequence analysis revealed that full length cDNA sequence of MeP5CS was 2 217 bp，which encoded 738 amino acids and contained P5CS conserved domain. The full length cDNA sequence of MeP5CR was 828 bp，which encoded 275 amino acids and contained P5CR conserved domain. They both had close genetic relationship with P5CS and P5CR，respectively，in Jatropha curcas and Ricinus communis. qRT-PCR results showed that the expression of both MeP5CS and MeP5CR were regulated by osmotic stress. Key words： cassava（Manihot esculenta Crantz）； MeP5CS； MeP5CR； gene clone； drought； expression analysis

在干旱、鹽漬等脅迫條件下，植物體內(nèi)脯氨酸含量顯著增加[1]。一定程度上，脯氨酸含量反映了植物的抗逆性，具有防止細(xì)胞脫水的作用，同時抗旱性強(qiáng)的品種往往能積累較多的脯氨酸[2]。因此，測定脯氨酸含量可以作為抗旱育種的生理指標(biāo)之一。木薯（Manihot esculenta Crantz）是重要的糧食作物和經(jīng)濟(jì)作物，具有抗旱的特性[3]。然而，長時間或較為嚴(yán)重的干旱脅迫嚴(yán)重地影響了木薯的生長和發(fā)育（如光合作用下降等），導(dǎo)致木薯塊根產(chǎn)量減少[3]。目前，對木薯抗旱分子機(jī)制的了解還十分有限。

近年來，探討脯氨酸合成相關(guān)酶基因在轉(zhuǎn)基因植物中的作用機(jī)制成了脯氨酸基因工程研究中的熱點。將脯氨酸生物合成關(guān)鍵酶基因?qū)胫仓?，以期獲得轉(zhuǎn)基因株系提高其對滲透脅迫的抵御或耐受力，進(jìn)而培育出高抗逆性的新品種。脯氨酸合成途徑可以分為兩種：谷氨酸途徑和鳥氨酸途徑[4]。大部分研究認(rèn)為，在受干旱脅迫時，谷氨酸途徑是脯氨酸生物合成的主要途徑，而在該途徑中存在兩個關(guān)鍵酶，即吡咯啉-5-羧酸合成酶（P5CS）和吡咯琳-5-羧酸還原酶（P5CR）[5]，研究者從甘蔗[6]、柑橘[7]和煙草[8]等植物中克隆了這兩個酶基因，并通過轉(zhuǎn)基因的方法研究其功能，但在木薯中鮮見相關(guān)報道。

本研究在不同濃度（0、20%、30%、40%和50%）PEG脅迫處理下考察了木薯葉片中脯氨酸含量的動態(tài)變化，從木薯中分別克隆了MeP5CS和MeP5CR基因，對其序列進(jìn)行初步的生物信息學(xué)分析，并分析了木薯MeP5CS和MeP5CR基因在20% PEG脅迫處理下各組織中的表達(dá)，旨在闡述干旱脅迫對木薯葉片脯氨酸含量和介入脯氨酸代謝基因表達(dá)的影響，為后續(xù)轉(zhuǎn)基因工作奠定基礎(chǔ)。

1 材料與方法

1.1 材料及干旱處理

試驗材料木薯為華南124（SC124），具有較強(qiáng)的抗旱性[1]，試驗在溫室大棚中進(jìn)行。將粗細(xì)均勻一致、長約15 cm的木薯種莖扦插于塑料盆（高18.8 cm，上直徑18.5 cm，下直徑14.8 cm）種植，栽培基質(zhì)選用蛭石和沙子同體積混合，種植10 d后進(jìn)行間苗，保留1苗/盆。種植60 d后（植株約0.5 m高，6～8片葉），在上午8∶00（記為0 h）用不同濃度（0、20%、30%、40%、50%）的PEG 6 000溶液進(jìn)行模擬干旱處理，以澆灌自來水為對照，在處理后0、3、6、12、24、48 h分別收集木薯葉片，采用磺基水楊酸法[2]測定脯氨酸含量。每個處理3次重復(fù)。

1.2 基因克隆

用RNA Plant Plus Reagent試劑盒[天根生化科技（北京）有限公司]提取木薯SC124葉片總RNA，用Revert Aid First Strand cDNA Synthesis Kit 試劑盒（Fermentas公司）合成cDNA用于后續(xù)試驗。

根據(jù)Phytozome（https：//phytozome.jgi.doe.gov）木薯數(shù)據(jù)庫和木薯基因組數(shù)據(jù)庫[9]，獲得2條與脯氨酸合成相關(guān)的序列，二者均包含了一個完整的開放閱讀框（ORF）。根據(jù)ORF采用生物學(xué)軟件Primer Premier 5.0設(shè)計引物（表1），進(jìn)行PCR擴(kuò)增。反應(yīng)程序：94 ℃ 2 min；94 ℃ 30 s，55 ℃ 30 s，72 ℃ 2 min，30個循環(huán)；72 ℃ 10 min。將擴(kuò)增得到的片段連接到pMD19-T載體后，轉(zhuǎn)化DH5α感受態(tài)細(xì)胞，挑選陽性克隆。經(jīng)菌落PCR和質(zhì)粒酶切驗證后送上海生工生物工程有限公司測序。

1.3 生物信息學(xué)分析

將測序得到的核苷酸序列和推導(dǎo)的氨基酸序列采用BLAST軟件（http：//blast.ncbi.nlm.nih.gov/Blast.cgi）進(jìn)行序列搜索；采用CDD數(shù)據(jù)庫（http：//www.ncbi.nlm.nih.gov/cdd/）進(jìn)行保守結(jié)構(gòu)域預(yù)測；采用Clustal X軟件進(jìn)行序列比對；利用MEGA 6.0軟件內(nèi)置的Neighbor-Joining法（bootstrap=1 000）構(gòu)建系統(tǒng)進(jìn)化樹。

1.4 基因表達(dá)分析

為了研究干旱條件下不同組織中基因表達(dá)量的變化，用20%PEG處理木薯0、3、24 h后，分別采集了不同葉片（包括未展開葉、第一片完全展開葉和老葉）和根的樣本，分別用RNA Plant Plus Reagent試劑提取總RNA，用Revert Aid First Strand cDNA Synthesis Kit試劑盒合成cDNA，進(jìn)行qRT-PCR分析，引物見表1。

qRT-PCR采用SYBR GreenⅠ試劑盒（TaKaRa公司），按照操作說明和推薦的反應(yīng)體系在 Mx 3005P型熒光定量PCR儀（Stratagene公司）上進(jìn)行，反應(yīng)程序：95 ℃預(yù)變性30 s；95℃ 10 s，55℃ 10 s， 72 ℃ 20 s，40個循環(huán)。每個樣品3次重復(fù)，表達(dá)量按照2-ΔΔCt進(jìn)行計算[11]。

2 結(jié)果與分析

2.1 木薯葉片脯氨酸含量在不同干旱脅迫下的動態(tài)變化

為了研究木薯干旱脅迫下脯氨酸含量的動態(tài)變化，在5個不同的PEG濃度（0%、20%、30%、40%、50%）、6個不同的時間點（0、3、6、12、24和48 h）處理下分別測定了木薯葉片中的脯氨酸含量。如圖1所示，在不同PEG濃度處理下，脯氨酸含量在處理后0～24 h均呈線性增長，24～48 h在20%、30%、40%PEG處理下，脯氨酸含量隨著PEG濃度升高而明顯升高；然而當(dāng)PEG濃度升至50%時，脯氨酸含量下降，表明木薯在高濃度（如50%）PEG下處理較長時間（如>24 h）可能會破壞其細(xì)胞滲透調(diào)節(jié)的能力。干旱脅迫能誘導(dǎo)木薯葉片中脯氨酸含量的增加；在以后的木薯干旱試驗中，應(yīng)該選擇適當(dāng)?shù)腜EG濃度（如20%、30%）來模擬較長時間的干旱脅迫。

2.2 MeP5CS和MeP5CR基因克隆與序列分析

以木薯SC124 cDNA為模板，用基因特異性引物進(jìn)行PCR擴(kuò)增，克隆目的基因的全長cDNA序列（圖2）。將擴(kuò)增產(chǎn)物連接pMD19-T載體，經(jīng)測序、序列比對后，分別命名為MeP5CS和MeP5CR。

MeP5CS基因全長編碼區(qū)序列為2 217 bp，編碼738個氨基酸（圖3A）。將其與文庫中獲得的參考序列比對，存在9個堿基差異，其中1個為非同義突變，這可能是由于不同木薯品種之間的差異造成的。保守結(jié)構(gòu)域分析顯示，MeP5CS含有P5CS保守結(jié)構(gòu)域（PLN02418，圖3B），證明克隆的就是MeP5CS基因。BLASTN分析顯示，MeP5CS與麻風(fēng)樹（XM_012233139）和蓖麻（XM_002524184）P5CS基因的同源性較高，核苷酸序列相似性分別達(dá)到89%和88%。用蛋白質(zhì)序列構(gòu)建MeP5CS基因的進(jìn)化樹，分析結(jié)果也表明，木薯MeP5CS與麻風(fēng)樹（Jatropha curcas，XP_012088529）和蓖麻（Ricinus communis，XP_002524230）P5CS的親緣關(guān)系較近（圖3C），而與水稻（Oryza sativa，BAA19916）和小麥（Triticum aestivum，AAX35536）的親緣關(guān)系較遠(yuǎn)。

MeP5CR基因含有一個完整的開放讀碼框，長度為828 bp，編碼275個氨基酸（圖4A）。將其與數(shù)據(jù)庫中獲得的木薯參考序列比對發(fā)現(xiàn)，兩者之間存在9個堿基差異，但無氨基酸差異。BLASTN分析表明，MeP5CR與蓖麻P5CR基因（XM_002531450）的同源性最高，核苷酸相似性為84%。保守結(jié)構(gòu)域分析顯示，MeP5CR含有P5CR保守結(jié)構(gòu)域（PLN02688，圖4B），證明克隆的就是MeP5CR基因。用蛋白質(zhì)序列構(gòu)建MeP5CR基因的進(jìn)化樹，分析結(jié)果表明，MeP5CR與蓖麻（Ricinus communis，XP_002531496）、麻風(fēng)樹（Jatropha curcas，NP_001303280）P5CR基因的親緣關(guān)系較近，而與玉米（Zea mays，NP_001105787）、黑麥草（Lolium perenne，AGS78291）和小麥（Triticum aestivum，AAW82908）的親緣關(guān)系較遠(yuǎn)（圖4C）。

2.3 MeP5CS和MeP5CR基因表達(dá)分析

為了研究MeP5CS和MeP5CR對干旱脅迫的響應(yīng)，用20%PEG處理木薯0、3和24 h后，分別考察了MeP5CS和MeP5CR在木薯不同葉片（包括未展開葉、第一片完全展開葉和老葉）和根中表達(dá)量的變化情況（圖5）。

PEG處理3 h后，MeP5CS在第一片完全展開葉中的表達(dá)量稍有下降，在老葉中有輕微上升，然而在PEG處理24 h后，其表達(dá)量均明顯上升。相比之下，MeP5CS在未展開葉和根中的表達(dá)量變化不大，表明MeP5CS主要在第一片完全展開葉和老葉中起作用。

在未展開葉中，MeP5CR的表達(dá)量在3、24 h都明顯下降；在第一片完全展開葉中，MeP5CR的表達(dá)量在3 h明顯下降，然后在24 h時升高；在老葉中，MeP5CR的表達(dá)量呈現(xiàn)相反的表達(dá)趨勢，在3 h先上升然后在24 h下降；在根中，MeP5CR的表達(dá)量在3、24 h均輕微下降。上述結(jié)果表明，MeP5CS和MeP5CR在轉(zhuǎn)錄水平受到滲透脅迫的調(diào)控。

3 小結(jié)與討論

在干旱等逆境因素脅迫下，許多植物通常會大量積累脯氨酸、可溶性糖、可溶性蛋白等滲透物質(zhì)。脯氨酸不僅是植物體內(nèi)的氮源、碳源，作為滲透調(diào)節(jié)物質(zhì)，它還在植物抵御滲透脅迫中起著重要作用[12]。有研究表明，脯氨酸累積與植物對干旱和鹽脅迫適應(yīng)性之間呈正相關(guān)[8，13]。本研究中，在PEG處理前期（0～24 h）木薯葉片中脯氨酸含量隨著處理時間的延長呈線性增長，而在處理后期（24～48 h）其含量隨著PEG濃度（20%～40%）升高而明顯升高，表現(xiàn)出較好的相關(guān)性，這與前人的研究結(jié)果[8，13]一致。

植物體內(nèi)，脯氨酸的合成是通過谷氨酸途徑和鳥氨酸途徑進(jìn)行的。有研究認(rèn)為滲透脅迫下，兩種途徑在脯氨酸積累中都發(fā)揮重要作用[14，15]；也有研究認(rèn)為在滲透脅迫與氮缺乏時谷氨酸途徑是植物體內(nèi)脯氨酸合成的主要途徑，而在非脅迫與高氮水平條件下鳥氨酸途徑占優(yōu)勢[16]。P5CS和P5CR是脯氨酸合成谷氨酸途徑中的兩個關(guān)鍵基因，這兩個基因的克隆、分離和遺傳轉(zhuǎn)化的研究已取得很大進(jìn)展[17-20]。本研究從木薯葉片中分別克隆了MeP5CS和MeP5CR基因。表達(dá)分析表明，MeP5CS在第一片完全展開葉和老葉中的表達(dá)量受到滲透脅迫誘導(dǎo)，其表達(dá)量變化趨勢與脯氨酸含量上升趨勢較為一致，表明MeP5CS與脯氨酸的合成有關(guān)。相比之下，MeP5CR在各葉片中的表達(dá)量不但沒有上升，反而下降了，與脯氨酸含量的變化趨勢相反，表明MeP5CR的表達(dá)量受到脯氨酸含量的反饋抑制。盡管這一現(xiàn)象在其他物種中也有報道[21]，但其相關(guān)機(jī)制還有待進(jìn)一步探討。

參考文獻(xiàn)：

[1] 王曉慶.耐旱木薯種質(zhì)資源的鑒定與評價[D].海口：海南大學(xué)，2012.

[2] 李合生.植物生理生化實驗原理和技術(shù)[M].北京：高等教育出版社，2000.

[3] OKOGBENIN E，SETTER T L，F(xiàn)ERGUSON M，et al. Phenotypic approaches to drought in cassava：Review[J]. Fron Physiol，2013，4：1-15.

[4] KAVIKISHOR P B，HONG Z，MIAO G H，et al. Overexpression of △1-pyrroline-5-carboxylate synthestase increases proline production and confers osmotolerance in transgenic plant[J]. Plant Physiol，1995，108：1387-1394.

[5] BATES L S，WALDREN R P，TEARE I D. Rapid determination of free proline for water-stress studies[J].Plant and Soil，1973， 39：205-207.

[6] DELAUNEY A J，HU C A A，KISHOR P B K，et al. Cloning of ornithine-aminotransferase cDNA from Vigna aconitifolia by trans-complementation in Escherichia coli and regulation of proline biosynthesis[J]. J Biol Chem，1993，268：18673-18678.

[7] 王麗媛，丁國華，黎 莉.脯氨酸代謝的研究進(jìn)展[J].哈爾濱師范大學(xué)自然科學(xué)學(xué)報，2010，26（2）：84-89.

[8] MOLINARI H B C，MARUR C J，DAROS E，et al. Evaluation of the stress-inducible production of proline in transgenic sugarcan（Saccharum spp.）：Osmotic adjustment，chlorophyll fluorescence and oxidative stress[J]. Physiol Plantarum，2007，130（2）：218-229.

[9] WANG W Q，F(xiàn)ENG B X，XIAO J F，et al. Cassava genome from a wild ancestor to cultivated varieties[J]. Nat Commun， 2014，5：1-9.

[10] XU J，DUAN X G，YANG J，et al. Enhanced reactive oxygen species scavenging by overproduction of superoxide dismutase and catalase delays postharvest physiological deterioration of cassava storage roots[J]. Plant Physiol，2013，161（3）：1517-1528.

[11] 胡 偉，顏 彥，韋運謝，等.木薯MeASR基因克隆及表達(dá)分析[J].生物技術(shù)通報，2015，31（10）：125-130.

[12] MOLINARI H B C，MARUR C J，F(xiàn)ILHO J C B，et al. Osmotic adjustment in transgenic citrus rootstock Carrizo citrange （Citrus sinensis Osb.×Poncirus trifoliate L. Raf.） overproducting proline[J]. Plant Sci，2004，167（6）：1375-1381.

[13] 李 玲，余光輝，曾富華.水分脅迫下植物脯氨酸累積的分子機(jī)理[J].華南師范大學(xué)學(xué)報（自然科學(xué)版），2003（1）：126-134.

[14] HONG Z，LAKKINENI K，ZHANG Z，et al. Removal of feedback inhibition of △1-pyrroline-5-carboxylate synthetase results in increased proline accumulation and protection of plants from osmotic stress[J]. Plant Physiol，2000，122：1129-1136.

[15] 楊 善，葉昌輝，莫俊杰，等.不同磷肥水平對甘蔗脯氨酸合成積累及關(guān)鍵基因表達(dá)的效應(yīng)[J].華南農(nóng)業(yè)大學(xué)學(xué)報，2015， 36（6）：68-75.

[16] HUA X J，VAN D C B，VAN M M，et al. Develpmental regulation of pyrroline-5-carboxylate reductase gene expression in Arabidopsis[J]. Plant Physiol，1997，114：1215-1224.

[17] 徐 晟，江宜龍，蔣明敏，等.石蒜LrP5CS基因的克隆與表達(dá)特性分析[J].西北植物學(xué)報，2015，35（10）：1925-1933.

[18] 周精華，邢虎成，揭雨成，等.苧麻Δ1-吡咯啉-5-羧酸合成酶（P5CS）基因的克隆和表達(dá)分析[J].作物學(xué)報，2012，38（3）：549-555.

[19] 曹 麗，孫振元，義鳴放，等.多年生黑麥草P5CS基因的cDNA克隆、表達(dá)及亞細(xì)胞定位[J].園藝學(xué)報，2010，37（9）：1477-1484.

[20] 馮昌增，謝月輝，曹 毅，等.果蠅吡咯啉-5-羧酸還原酶的表達(dá)、純化及理化特點[J].中國生物工程雜志，2012，32（6）：27-34.

[21] P？魪REZ-ARELLANO I，CARMONA-？魣LVAREZ F，MART？魱NEZ A，et al.Pyrroline-5-carboxylate synthase and proline biosynthesis：From osmotolerance to rare metabolic disease[J].Protein Sci，2010，19（3）：372-382.