崔 翠，王利鵑，周清元，譚尊飛，曲存民，張正圣

(西南大學(xué), 農(nóng)學(xué)與生物科技學(xué)院，重慶 400715)

低溫脅迫下烤煙幼苗葉片光合作用和抗氧化能力基因差異表達(dá)譜

崔 翠*，王利鵑，周清元，譚尊飛，曲存民，張正圣

(西南大學(xué), 農(nóng)學(xué)與生物科技學(xué)院，重慶 400715)

對(duì)低溫(5—7 ℃)脅迫下烤煙“K326”幼苗葉片光合指標(biāo)、膜氧化水平及其抗氧化指標(biāo)進(jìn)行測(cè)定，并利用數(shù)字化基因表達(dá)譜技術(shù)進(jìn)行基因差異表達(dá)分析。低溫脅迫后烤煙幼苗葉綠素含量、光合能力顯著下降，脯氨酸含量、丙二醛含量上升，超氧化物歧化酶活性、過(guò)氧化氫酶活性、抗壞血酸含量和谷胱甘肽含量均顯著上升。低溫脅迫后有2357個(gè)基因發(fā)生了顯著差異表達(dá)，其中1673個(gè)基因表達(dá)上調(diào)、684個(gè)基因表達(dá)下調(diào)，其分子功能、細(xì)胞位置和主要代謝過(guò)程均涉及光系統(tǒng)、膜氧化系統(tǒng)和抗氧化系統(tǒng)。對(duì)涉及到的代謝過(guò)程進(jìn)行分析，結(jié)果表明:光合天線蛋白調(diào)控基因表達(dá)量均顯著下降、光合作用的主要調(diào)控基因表達(dá)量多數(shù)表現(xiàn)為顯著下調(diào)、而與氧化能力相關(guān)的谷胱甘肽代謝差異表達(dá)基因大多數(shù)顯著上調(diào)?；虿町惐磉_(dá)譜分析結(jié)果和低溫脅迫后葉片光合能力、抗氧化能力生理生態(tài)指標(biāo)測(cè)定結(jié)果基本一致，為進(jìn)一步研究低溫脅迫對(duì)作物的生態(tài)影響和研究基因克隆與功能提供基礎(chǔ)。

烤煙；低溫脅迫；數(shù)字化基因表達(dá)譜；抗氧化能力

低溫冷害是限制農(nóng)業(yè)生產(chǎn)的主要因子之一。低溫會(huì)使葉綠素分解，細(xì)胞受到傷害，質(zhì)膜相對(duì)透性增大[1]，葉綠素含量降低[1- 4]。葉綠體是低溫脅迫的直接目標(biāo)[5]，低溫削弱植物通過(guò)光合作用利用光能的能力[6- 8]。植物在一定低溫范圍內(nèi)通過(guò)提高保護(hù)酶含量來(lái)維持體內(nèi)自由基產(chǎn)生與清除間的動(dòng)態(tài)平衡。若低溫脅迫嚴(yán)重，就會(huì)導(dǎo)致自由基含量增高，造成膜脂過(guò)氧化，丙二醛含量上升，細(xì)胞膜系統(tǒng)受到傷害[9]。溫度逆境也誘導(dǎo)植物本身具有的抗氧化物質(zhì)，如過(guò)氧化物酶、過(guò)氧化氫酶、超氧化物歧化酶、谷胱甘肽及抗壞血酸等含量上升來(lái)抵御逆境傷害[1,10]。

生物對(duì)各種生物、理化及致病因子的應(yīng)答，本質(zhì)上都是基因差異表達(dá)的結(jié)果。比較同一類細(xì)胞或不同類細(xì)胞在不同生理狀態(tài)或生長(zhǎng)發(fā)育階段下的基因表達(dá)差異，對(duì)研究生物生長(zhǎng)發(fā)育調(diào)控及代謝機(jī)制具有重要意義[11]。近年來(lái)，利用mRNA 差異顯示技術(shù)和基因芯片技術(shù)等研究逆境對(duì)作物生長(zhǎng)發(fā)育的影響取得了重要進(jìn)展，對(duì)作物逆境生理和遺傳育種研究發(fā)揮了積極作用[12- 16]。高通量測(cè)序的數(shù)字化基因表達(dá)譜(Digital gene expression profile, DGE)技術(shù)是新發(fā)展起來(lái)的分析生物轉(zhuǎn)錄組基因表達(dá)的方法, 能更加真實(shí)、全面地評(píng)價(jià)細(xì)胞中全部基因的表達(dá)情況，也能夠檢測(cè)未知基因的表達(dá),是發(fā)現(xiàn)新基因的一種有效方法[17- 18]。在轉(zhuǎn)錄水平上對(duì)外界非生物脅迫下植物應(yīng)答機(jī)制的研究主要集中在擬南芥、水稻等植物上，脅迫類型以干旱、冷害、鹽脅迫等為主[19- 23]。關(guān)于煙草在低溫脅迫下基因的特異表達(dá)研究較少，本研究模擬南方煙區(qū)“倒春寒”，研究低溫脅迫對(duì)烤煙幼苗光合能力、膜氧化及抗氧化能力的影響，再通過(guò)高通量的數(shù)字化基因表達(dá)譜技術(shù)，分析低溫脅迫下烤煙幼苗抗寒相關(guān)基因的差異表達(dá)譜，尋找低溫脅迫后與煙草光合作用、抗氧化能力有關(guān)的功能基因，為深入分析煙草低溫脅迫傷害機(jī)理、挖掘抗寒相關(guān)基因，提高煙草適應(yīng)性栽培等方面提供理論指導(dǎo)。

1 材料與方法

1.1 材料培養(yǎng)與處理

試驗(yàn)材料為生產(chǎn)上大面積栽培的Msk326(K326)，由云南煙草科學(xué)研究所選育。選取大小一致且飽滿的烤煙種子浸種催芽出苗后，分別移栽于裝有營(yíng)養(yǎng)土的塑料缽內(nèi)，在3葉期選取大小基本一致煙苗，每缽定苗3株，置于光照培養(yǎng)箱中培養(yǎng)生長(zhǎng)。培養(yǎng)條件：溫度為23—25 ℃；光照強(qiáng)度為6000 lx，每日光照12 h。煙苗長(zhǎng)至5—6片真葉時(shí)，5—7 ℃低溫脅迫處理3 d(標(biāo)記為B)，常溫23—25 ℃培養(yǎng)為對(duì)照(標(biāo)記為A)。取其相同部位葉片進(jìn)行RNA提取和相關(guān)生理指標(biāo)測(cè)定。

總RNA提取試劑盒為WASON(RNAex Reagent System Ⅳ)試劑盒，上海華舜生物工程有限公司生產(chǎn)，其它均采用國(guó)產(chǎn)分析純?cè)噭?/p>

1.2 試驗(yàn)方法

1.2.1 生理生化指標(biāo)測(cè)定

取處理(B)和對(duì)照(A)葉片，參考崔翠方法[24]，分別測(cè)定總?cè)~綠素Chl 含量，Chl a 及Chl b含量，Pn，Gs，Tr，Ci，電解質(zhì)滲透率，Pro及MDA含量，POD、CAT和SOD活性，GSH及Vc含量等生理指標(biāo)。每處理4次重復(fù)，測(cè)定數(shù)據(jù)利用Micorsoft excel 軟件進(jìn)行處理，利用DPS統(tǒng)計(jì)軟件進(jìn)行差異顯著性檢驗(yàn)。

1.2.2 低溫脅迫后烤煙幼苗葉片基因差異表達(dá)譜的分析

RNA提取方法參考試劑盒操作手冊(cè)進(jìn)行，分別取對(duì)照(A)和處理(B)葉片混合成兩個(gè)樣品池，從兩個(gè)樣品池中各取混合煙葉葉片組織0.5 g左右，提取6 μg總RNA，-78 ℃保存。測(cè)序由深圳華大基因科技有限公司承擔(dān)。采用邊合成邊測(cè)序法(Sequencing by synthesis，SBS)測(cè)序[25]。測(cè)序所得原始圖像數(shù)據(jù)經(jīng)base calling 轉(zhuǎn)化為序列數(shù)據(jù)。對(duì)測(cè)序得到的原始標(biāo)簽去除3′adaptor 序列、低質(zhì)量序列(含未知堿基)以及長(zhǎng)度過(guò)小或過(guò)大，拷貝數(shù)為1的tag，確定Clean tag。對(duì)初步篩選獲得的Clean tags進(jìn)行標(biāo)準(zhǔn)化處理[26]，統(tǒng)計(jì)不同tag的分布特征，分析測(cè)序飽和度。將Clean tags 與基因序列進(jìn)行比對(duì)，確定Clean tag 代表的基因。通過(guò)tags 統(tǒng)計(jì)相應(yīng)基因的表達(dá)量。

從NCBI網(wǎng)站檢索mRNA Unigene序列上所有CATG位點(diǎn)，生成CATG+17堿基的參考標(biāo)簽數(shù)據(jù)庫(kù)。將全部clean tags與參考標(biāo)簽數(shù)據(jù)庫(kù)比對(duì)，允許最多1個(gè)堿基錯(cuò)配，對(duì)其中唯一比對(duì)到1個(gè)基因的標(biāo)簽(Unambiguous tags)進(jìn)行基因注釋，統(tǒng)計(jì)每個(gè)基因?qū)?yīng)的原始的clean tag數(shù)目，然后對(duì)原始clean tag數(shù)做標(biāo)準(zhǔn)化處理，獲得標(biāo)準(zhǔn)化的基因表達(dá)量，衡量基因的表達(dá)水平。標(biāo)準(zhǔn)化的方法為：每個(gè)基因包含的原始clean tag 數(shù)/該樣本中總clean tag 數(shù)×1000000[26]。參照Audic S.等發(fā)表在Genome Research上的數(shù)字化基因表達(dá)譜差異基因檢測(cè)方法[17]。

參考李余良等方法[15]把所有差異表達(dá)基因向GO(Gene ontology)數(shù)據(jù)庫(kù)(http：//www.genetntology.org/)的各個(gè)項(xiàng)目映射，找出與整個(gè)Tag庫(kù)背景相比，在差異表達(dá)基因中顯著富集的GO條目和顯著性富集的Pathway，分析其中與光合及抗氧化能力相關(guān)的差異表達(dá)基因。

2 結(jié)果與分析

2.1 低溫脅迫處理對(duì)烤煙幼苗葉片部分光合能力和抗氧化指標(biāo)的影響

葉綠素與光合作用密切相關(guān)，同時(shí)也是衡量植物抗寒力高低的一項(xiàng)重要指標(biāo)。結(jié)果如表1所示，低溫處理3 d后的葉綠素含量，無(wú)論是葉綠素a (Chl a)還是葉綠素b (Chl b)含量均顯著下降，說(shuō)明低溫脅迫嚴(yán)重阻礙了葉綠素合成；低溫處理后烤煙葉片光合能力降低，低溫導(dǎo)致凈光合速率、氣孔導(dǎo)度、蒸騰速率顯著下降。氧自由基的積累對(duì)細(xì)胞存在著巨大傷害，處理3 d后烤煙幼苗葉片電解質(zhì)滲透率、脯氨酸含量及丙二醛含量均顯著上升，說(shuō)明低溫脅迫已對(duì)烤煙幼苗葉片細(xì)胞膜造成一定傷害，但其氧自由基生成速率(形成速率)顯著下降；另外，膜保護(hù)酶中CAT及SOD活性顯著上升，而POD活性則顯著下降；抗氧化物質(zhì)Vc和GSH含量均顯著上升。從結(jié)果來(lái)看，低溫脅迫3 d后，K326幼苗葉片光合能力下降，細(xì)胞膜氧化，但其自身抗氧化能力顯著提升。低溫處理3d后烤煙K326幼苗在光合能力、膜氧化及抗氧化能力方面均發(fā)生顯著變化。

表1 低溫脅迫后烤煙幼苗葉片部分生理指標(biāo)

*表示差異顯著(P<0.05)

2.2 基因差異表達(dá)譜分析

2.2.1 RNA質(zhì)量檢測(cè)和測(cè)序質(zhì)量評(píng)估

分別提取低溫脅迫(B)及常溫培養(yǎng)(A)的烤煙幼苗葉片的總RNA，樣品經(jīng)NANODrop 2000C 核酸檢測(cè)儀檢測(cè)，OD260/OD280均在1.8—2.2范圍內(nèi)，純度較高。按照檢測(cè)的濃度，將相同處理的總RNA等量混合。用DNase I降解基因組DNA，電泳顯示RNA均保持完整性?？俁NA用電泳檢測(cè)表明，28S rRNA和18S rRNA條帶清晰可見(jiàn)，無(wú)拖帶現(xiàn)象，說(shuō)明RNA完整性較好，可滿足實(shí)驗(yàn)要求。

測(cè)序得到對(duì)照樣品A 和處理樣品B 原始序列數(shù)據(jù)，經(jīng)去除雜質(zhì)后得到高質(zhì)量tag(Clean tag)分別有3490430條和3506495條，占總tag數(shù)的95.91%和94.69%；僅有5%左右的tag拷貝數(shù)小于2或只有接頭序列。通過(guò)測(cè)序分析，最終得到的clean tag的數(shù)量均達(dá)到總測(cè)序量90%以上，說(shuō)明樣品制備和測(cè)序質(zhì)量良好。對(duì)2個(gè)樣品Clean tag拷貝數(shù)分布進(jìn)行分析，顯示拷貝數(shù)大于100的高表達(dá)Tags 在對(duì)照和處理中分別占63.79%和55.50%，在數(shù)量上占絕對(duì)優(yōu)勢(shì)； 拷貝數(shù)小于5的低表達(dá)Tags雖然在總拷貝數(shù)中分別只占5.94%和7.24%，但是在種類上非常豐富, 均達(dá)到57%以上，從整體上評(píng)估測(cè)序數(shù)據(jù)正常。

從NCBI網(wǎng)站上(http：//www.ncbi.nlm.nih.gov/)下載和煙草相關(guān)的參考基因有24432條，參考tag有74261條。通過(guò)軟件檢索Unigene上所有CATG位點(diǎn)，生成CATG+17堿基的參考標(biāo)簽數(shù)據(jù)庫(kù)，其中參考基因中有CATG位點(diǎn)的基因有22263個(gè)，唯一的參考標(biāo)簽有72259個(gè)。將clean tags分別與參考標(biāo)簽數(shù)據(jù)庫(kù)比對(duì)，并統(tǒng)計(jì)每個(gè)基因?qū)?yīng)clean tag的數(shù)目。對(duì)照(A)一個(gè)tag完全比對(duì)上一個(gè)基因的標(biāo)簽(允許1個(gè)堿基錯(cuò)配)數(shù)目為1066364，占總數(shù)的30.55%，種類為16983種，占總數(shù)的13.47%；而處理(B)一個(gè)tag完全比對(duì)上一個(gè)基因的標(biāo)簽(允許1個(gè)堿基錯(cuò)配)數(shù)目為889559，占總數(shù)的25.37%，種類為18619，占總數(shù)的12.03%。

2.2.2 低溫脅迫后烤煙幼苗葉片基因表達(dá)譜分析

經(jīng)過(guò)Unigene標(biāo)簽庫(kù)進(jìn)行基因注釋后，通過(guò)基因比對(duì)，設(shè)定FDR≤0.001[25]，且倍數(shù)差異在2倍以上(∣log2ratio∣≥1)作為閾值來(lái)分析判斷基因在2個(gè)溫度狀態(tài)下是否差異表達(dá)，篩選出差異表達(dá)基因。低溫脅迫后(處理B)與對(duì)照A相比較(圖1)，有2357個(gè)基因發(fā)生了顯著的差異表達(dá)，其中有1673個(gè)基因表達(dá)上調(diào)，log2ratio值大于10的有62個(gè)，5.00 — 9.99之間的有185個(gè)，1.00 — 4.99之間有1426個(gè)；684個(gè)基因表達(dá)下調(diào)，log2ratio值小于 -10的有11個(gè)， -9.99— -5.00之間的有51個(gè)， -4.99— -1.00之間的有622個(gè)。

圖1 低溫脅迫下烤煙幼苗差異表達(dá)基因Fig.1 Differentially expressed genes(DEGs)of leaves of tobacco seedling under chilling stress

通過(guò)GO功能富集分析(Gene ontology，簡(jiǎn)稱GO)，對(duì)差異表達(dá)顯著的基因進(jìn)行功能分類注釋(圖2)。圖2結(jié)果表明，從差異表達(dá)基因所處的細(xì)胞位置來(lái)看，低溫脅迫后烤煙幼苗與對(duì)照煙苗發(fā)生顯著差異表達(dá)的基因主要集中在光合系統(tǒng)(Photosystem)、光合膜(Photosynthetic Membrane)、類囊體(Thylakoid Part)、葉綠體(Chloroplast)、質(zhì)體(Plastid Part)、細(xì)胞質(zhì)(Cell Plastid)、細(xì)胞膜(Membrane)等位置，其中與細(xì)胞質(zhì)有關(guān)基因發(fā)生差異表達(dá)最為顯著，有524個(gè)基因發(fā)生差異表達(dá)，占總差異表達(dá)基因的45.1%；與生物膜有關(guān)的差異表達(dá)基因達(dá)到420個(gè)，占總差異表達(dá)基因的36.1%；細(xì)胞器部分，差異表達(dá)基因達(dá)到358個(gè)，占總差異表達(dá)基因的30.8%；此外，質(zhì)體及細(xì)胞內(nèi)細(xì)胞器部分等位置差異表達(dá)基因亦均達(dá)到27%以上。從分子功能分析結(jié)果看，低溫脅迫導(dǎo)致烤煙幼苗基因的氧化還原能力、跨膜運(yùn)輸能力、結(jié)構(gòu)分子活力發(fā)生顯著改變。其中與氧化還原能力相關(guān)基因差異最為顯著，差異表達(dá)基因達(dá)到196個(gè)，占總差異表達(dá)基因的16.3%；與跨膜運(yùn)輸能力相關(guān)的差異表達(dá)基因?yàn)?04個(gè)，占總差異表達(dá)基因的8.6%；在結(jié)構(gòu)分子活力方面，差異表達(dá)基因達(dá)到74個(gè)，占總差異表達(dá)基因的6.1%。發(fā)生顯著差異表達(dá)的生物過(guò)程主要集中在脅迫響應(yīng)、代謝物質(zhì)及能源的初始反應(yīng)、光合及光反應(yīng)，對(duì)激素刺激的細(xì)胞響應(yīng)、激素調(diào)節(jié)的信號(hào)途徑等方面。其中脅迫刺激響應(yīng)的生物過(guò)程差異最為顯著，差異表達(dá)基因達(dá)到426個(gè)，占總差異表達(dá)基因的36.3%；與溫度脅迫響應(yīng)相關(guān)差異表達(dá)基因?yàn)?23個(gè)，占總差異表達(dá)基因的19.0%；與先導(dǎo)代謝物及能源相關(guān)差異表達(dá)基因74個(gè)，占總差異表達(dá)基因的6.3%。因此，從差異基因所處的位置、分子功能和影響的生物過(guò)程來(lái)看，調(diào)控烤煙幼苗光合作用和抗氧化能力的基因在低溫脅迫后發(fā)生了較大的差異表達(dá)，為進(jìn)一步分析提供了基礎(chǔ)。

2.2.3 低溫脅迫對(duì)烤煙幼苗葉片光合能力的影響

在生物體內(nèi)，不同基因相互協(xié)調(diào)行使其生物學(xué)功能，基于Pathway的分析有助于進(jìn)一步了解基因的生物學(xué)功能。KEGG是有關(guān)Pathway的主要公共數(shù)據(jù)庫(kù)[27]，以KEGG Pathway為單位，應(yīng)用超幾何檢驗(yàn)，進(jìn)行Pathway顯著性富集分析，在Qvalue≤0.05時(shí)，低溫對(duì)光合天線蛋白(捕光色素蛋白復(fù)合體，light-harvesting chlorophyll protein complex，簡(jiǎn)寫(xiě)LHC)和光合能力影響差異顯著。參照Kanehisa laboratories的光合作用天線蛋白示意圖3 (a)，標(biāo)記顯著差異基因于圖3(b)，并將對(duì)應(yīng)的差異表達(dá)基因列于表2。從圖3(b)和表2來(lái)看，光系統(tǒng)Ⅰ(photosystem Ⅰ，簡(jiǎn)寫(xiě)PSⅠ)和光系統(tǒng)Ⅱ(photosystem Ⅱ，簡(jiǎn)寫(xiě)PSⅡ)中部分捕光葉綠素a/b

圖2 低溫脅迫下烤煙幼苗葉片表達(dá)基因的功能分類Fig.2 Classification on functional of differentially expressed genes of leaves of tobacco seedling under chilling stress

圖3 KEGG數(shù)據(jù)庫(kù)中光合作用-天線蛋白示意圖(a); 低溫脅迫下烤煙幼苗與光合天線蛋白相關(guān)差異表達(dá)基因(b)Fig.3 The skechmap of antenna proteins related photosynthesis in KEGG database(a); Differential expression genes related photosynthesis-antenna proteins of tobacco seedling under chilling stress(b)綠色方框代表表達(dá)基因下調(diào)

色素天線蛋白(Lhca1、Lhca2、Lhca3、Lhca4、Lhca5、Lhcb1、Lhcb2、Lhcb3、Lhcb4、Lhcb5)等均在低溫脅迫下發(fā)生顯著下調(diào)。說(shuō)明低溫脅迫導(dǎo)致烤煙幼苗對(duì)光電子的捕獲能力下降。

表2 與光合作用-天線蛋白相關(guān)的差異表達(dá)基因相關(guān)信息

參照Kanehisa laboratories的光合作用示意圖4(a)，將差異表達(dá)基因分別標(biāo)記在圖4(b)，并將對(duì)應(yīng)基因列于表3。表3結(jié)果表明，低溫脅迫后烤煙幼苗與光合作用有關(guān)的差異表達(dá)基因有35個(gè)，其中11個(gè)基因表達(dá)上調(diào)，24個(gè)基因表達(dá)下調(diào)。從圖4(b)和表3來(lái)看，在PSⅡ上，上調(diào)表達(dá)的基因主要有影響PSⅡ中核心編碼蛋白CP47及分子量為44kDa、22kDa等蛋白質(zhì)相關(guān)的基因，這些基因主要控制光合作用中葉綠體類囊體上PsbC、PsbB、PsbS和PsbZ等蛋白；在下調(diào)表達(dá)的基因中，部分調(diào)節(jié)PSⅡ中水裂解酶、水氧化蛋白及分子量為10kDa的多肽等相關(guān)基因表達(dá)下調(diào)，導(dǎo)致PSⅡ中6個(gè)核心編碼蛋白(PsbO、PsbP、PsbQ、PsbR、PsbW和PsbY)表達(dá)下調(diào)；下調(diào)表達(dá)基因部分降低PSⅠ亞基(PSⅠ-D1、PSⅠ-H、Ⅲ)的前體及PSⅠ反應(yīng)中心亞基(Ⅱ、IX A、X psaK、XI、XIB)的生成，因此，導(dǎo)致PSⅠ中有9個(gè)核心編碼蛋白(PsaD、PsaE、PsbF、PsaG、PsaH、PsaK、PsbL、PsaO和PsaN)相關(guān)基因均在轉(zhuǎn)錄表達(dá)水平上顯著下調(diào)。細(xì)胞色素b/f復(fù)合物 (Cytochrome b6-f complex) 是類囊體膜上可分離出來(lái)的多亞基膜蛋白。在本研究中，控制光合電子傳遞的細(xì)胞色素b6/f復(fù)合物的部分基因Cytb6發(fā)生了上調(diào)表達(dá)(PetB蛋白表達(dá)上調(diào))、與光合電子傳遞鏈有關(guān)的蛋白PetF、PetH上調(diào)，而PetE蛋白下調(diào)。此外，ATP合成酶復(fù)合體有關(guān)的基因c蛋白基因上調(diào)，b蛋白基因下調(diào)及ATP合成酶(如ATP合酶δ鏈，ATP合酶c鏈)發(fā)生上調(diào)表達(dá)?？傊?，在光合作用中下調(diào)表達(dá)的24個(gè)基因中，涉及到光合作用的原初反應(yīng)、光合電子鏈的傳遞、光合磷酸化、光合碳同化等主要的代謝通路；部分基因的上調(diào)表達(dá)，也提高了烤煙幼苗對(duì)低溫脅迫的適應(yīng)。

圖4 KEGG數(shù)據(jù)庫(kù)中光合作用示意圖(a)；低溫脅迫下烤煙幼苗與光合作用相關(guān)的差異表達(dá)基因(b)Fig.4 The skechmap of photosynthesis in KEGG database(a); Differential expression genes related photosynthesis of tobacco seedling under chilling stress(b)紅色方框表示基因上調(diào)

基因登IDgeneIDlog2ratio假定功能Putativefunction圖中位置PositionO04397.1|FENR2_TOBAC1.109ferredoxin-NADPreductase,root-typeisozyme,chloroplastic1.18.1.2XP_002271630.1-1.1predicted:hypotheticalprotein3.6.3.14/bBAH11228.12.725ATPsynthaseCF0Cchain3.6.3.14/cXP_002523293.12.194ATPsynthasedeltachain,putative3.6.3.14/DeltAP32980.1|ATPD_TOBAC-1.9ATPsynthasedeltachain,chloroplastic3.6.3.14/DeltAYP_001936513.12.333photosystemII44kDaproteinPsbCADD30625.111.42photosystemIICP47chlorophyllapoproteinPscBAAP03871.1-3.54oxygenevolvingcomplex33kDaphotosystemIIproteinPsbOCAA45699.1-4.4223kDapolypeptideofwater-oxidizingcomplexofpho-tosystemIIPsbPCAA44292.1-4.3123-kDaploypeptideofphotosystemIIoxygen-evolvingcomplexPsbPP18212.2|PSBP2_TOBAC-2.43oxygen-evolvingenhancerprotein2-2,chloroplastic;PsbPXP_002279556.1-1.38predicted:hypotheticalproteinPsbPBAD97359.1-3.33PsbQPsbQAAU03361.1-2.65photosystemIIoxygen-evolvingcomplexprotein3PsbQQ40519.1|PSBR_TOBAC-1.78photosystemII10kDapolypeptide,chloroplastic;SPsbRQ9SMB4.1|PSBS_TOBAC1.039recName:Full=PhotosystemII22kDaprotein,chloroplasticPsbSACU13646.1-3.49unknownPsbWXP_002285325.1-2.02predicted:hypotheticalproteinPsbYYP_001109497.13.512hypotheticalproteinPoptr_cp018PsbZBAA02871.1-3.74PSⅠ-D1precursorPsaDP29302.1|PSAD_NICSY-2.64photosystemⅠreactioncentersubunitⅡ,chloroplas-tic;PsaDQ41228.1|PSAEA_NICSY-3.43photosystemⅠreactioncentersubunitchloroplastic;PsaEQ41229.1|PSAEB_NICSY-2.71photosystemⅠreactioncentersubunitⅣB,chloro-plasticPsaEAAP03872.1-4.11putativephotosystemⅠsubunitⅢprecursorPsaFCBI39855.3-2.2unnamedproteinproductPsaGBAA04634.1-3.62PSⅠ-HprecursorPsaHAAP03873.1-3.15photosystemIreactioncentersubunitXpsaKPsaKAAO85557.1-2.87photosystemIsubunitⅪPsaLACZ72945.1-3photosystemIreactioncentersubunitPsaNCAB75430.1-5.3putative16kDamembraneproteinPsaOP06247.2|CYB6_TOBAC4.362cytochromeb6PetBACV32157.1-2.64chloroplastplastocyaninprecursorpetEACU15698.15.601unknownPetFXP_002275771.11.1predicted:hypotheticalproteinisoform2PetFO04397.1|FENR2_TOBAC1.1ferredoxin-NADPreductase,root-typeisozyme,chloroplasticPetH

2.2.4 低溫脅迫對(duì)烤煙幼苗葉片抗氧化能力的影響

5—7 ℃低溫脅迫后K326烤煙幼苗代謝通路中，谷胱甘肽代謝途徑是基因顯著差異表達(dá)顯著富集的最主要生化代謝途徑之一。結(jié)合KEGG數(shù)據(jù)庫(kù)中谷胱甘肽代謝圖將差異基因標(biāo)注在圖5，并將相關(guān)基因信息列于表4。表4結(jié)果表明，低溫誘導(dǎo)烤煙幼苗谷胱甘肽代謝過(guò)程差異表達(dá)基因共有26個(gè)，其中上調(diào)表達(dá)的基因有23個(gè)，3個(gè)基因下調(diào)表達(dá)。上調(diào)表達(dá)的基因主要是與谷胱甘肽-s-轉(zhuǎn)移酶(GSTs)、抗壞血酸過(guò)氧化物酶、磷脂過(guò)氧化氫酶、谷胱甘肽過(guò)氧化物酶、6-磷酸葡萄糖酸鹽脫氫酶等酶蛋白有關(guān)的基因。其中，與谷胱甘肽S-轉(zhuǎn)移酶(GST)相關(guān)的基因有13個(gè)發(fā)生了上調(diào)表達(dá)，占總上調(diào)表達(dá)基因的50.00%；谷胱甘肽二硫化物和谷胱甘肽之間的代謝過(guò)程中差異基因8個(gè)，占總差異基因的29.63%，其中，APIC基因(基因登錄號(hào)為P46440.1|GSTF2_TOBAC)上調(diào)表達(dá)最顯著，較對(duì)照上升近1000倍，C- 7(基因登錄號(hào)為CAA45741.1)上調(diào)表達(dá)達(dá)到500倍，均在清除細(xì)胞內(nèi)活性氧及保護(hù)細(xì)胞抵御活性親電物質(zhì)時(shí)起重要作用。

表4 低溫脅迫烤煙幼苗谷胱甘肽代謝差異表達(dá)基因的相關(guān)信息

圖5 低溫脅迫下烤煙幼苗的谷胱甘肽代謝過(guò)程中的差異表達(dá)基因Fig.5 Differential expression genes related to glutathione metabolism of tobacco seedling under chilling stress

3 討論

煙草生長(zhǎng)最適溫度為25— 28 ℃。目前，烤煙大多采用工廠化集約漂浮育苗，雖然育苗小環(huán)境有所改善，但受自然氣溫影響仍然很大，當(dāng)溫度降至12 ℃以下會(huì)造成出苗慢、降低出苗率及苗期生長(zhǎng)緩慢。南方很多煙區(qū)在早春育苗階段常有“倒春寒”等低溫危害發(fā)生，一般3 d左右，有時(shí)甚至達(dá)1周，影響煙苗質(zhì)量及移栽后還苗期延長(zhǎng)，最終影響煙葉生產(chǎn)中、上等煙的比例及優(yōu)質(zhì)特色煙葉生產(chǎn)。晉艷等研究了低溫脅迫0.5—1.5 h后2個(gè)煙草品種(云煙85和K326)幼苗的膜透性、脯氨酸含量、膜保護(hù)酶和葉綠素含量變化特點(diǎn)[5]，模擬南方煙區(qū)“倒春寒”特點(diǎn)，低溫脅迫3 d后測(cè)定烤煙幼苗光合能力、膜氧化和抗氧化能力，研究結(jié)果對(duì)烤煙壯苗培育具有重要參考意義。本研究結(jié)果表明，低溫脅迫后其葉綠素含量、光合能力顯著下降，細(xì)胞膜受損，結(jié)果與短時(shí)間處理[5]結(jié)果一致，但膜保護(hù)酶活性和抗氧化物質(zhì)(谷胱甘肽和抗壞血酸)含量顯著上升，抗氧化能力加強(qiáng)，植物本身耐低溫能力得到提高。植物在逆境脅迫響應(yīng)下受到多個(gè)基因的調(diào)控，在逆境條件下大量基因表達(dá)發(fā)生變化，這些基因?qū)ο鄳?yīng)逆境耐性的獲得起到重要作用[15]?；跍y(cè)序技術(shù)的數(shù)字化基因表達(dá)譜具有高通量和高靈敏性的優(yōu)勢(shì)，不需要獲得目標(biāo)基因序列，可以從全基因組突破同時(shí)研究成千上萬(wàn)個(gè)基因的表達(dá)，能很好反映生物在生長(zhǎng)發(fā)育過(guò)程中整個(gè)基因組的基因表達(dá)變化情況[27- 28]。因此，我們采用該基因表達(dá)譜分析技術(shù)分析低溫處理3 d后烤煙幼苗葉片的基因差異表達(dá)情況，從中進(jìn)一步驗(yàn)證基因表達(dá)差異和光合能力、抗氧化能力之間的一致性，發(fā)現(xiàn)大量的上調(diào)和下調(diào)差異表達(dá)基因，其中部分基因表達(dá)量較高，為進(jìn)一步利用這些基因提供了理論依據(jù)。

在高等植物中大部分的葉綠素分子都結(jié)合在PSⅠ和PSⅡ的天然色素蛋白復(fù)合物上，類囊體膜是光合作用能量轉(zhuǎn)換所需的脂質(zhì)雙層膜，含有捕光及將光能轉(zhuǎn)化為化學(xué)能所必需的色素蛋白及酶系統(tǒng)。其中，PSⅡ捕光色素蛋白復(fù)合物在類囊體膜上的含量最為豐富,它所結(jié)合的葉綠素約占類囊體膜上色素量的50%[29]。通常PSⅠ捕光色素蛋白 (LHC)含有Lhca1、Lhca2、Lhca3 和Lhca4等種蛋白質(zhì)，PSⅡ捕光色素蛋白 (LHCII)含有Lhcb1、Lhcb2、Lhcb3、Lhcb4、Lhcb5等蛋白質(zhì)，捕光色素蛋白分子量在20 — 29 kD 之間[30]。內(nèi)囊體膜中有兩種捕光天線葉綠素a/b色素蛋白復(fù)合體，PSⅠ捕光天線葉綠素a/b 色素蛋白復(fù)合體 (LHCI) 和PSⅡ捕光天線葉綠素a/b色素蛋白復(fù)合體 (包括LHCII、CP29 、CP26 、CP24 和CP22)。由于LHCII 結(jié)合有葉片葉綠素總量的50%，在光合作用光能的吸收中具有重要作用。CP29、CP26、CP24 和CP22是PSⅡ捕光天線葉綠素a/b色素蛋白復(fù)合體中4種主要的亞復(fù)合物，含有β-胡蘿卜素和葉黃素，參與葉黃素循環(huán)，分散過(guò)多的激發(fā)能，對(duì)抗光抑制，起到保護(hù)反應(yīng)中心的作用，其中CP22和CP29更靠近PSⅡ核心復(fù)合物。本研究結(jié)果表明，低溫脅迫后葉綠素a/b結(jié)合蛋白編碼基因，無(wú)論是CP22還是CP29，表達(dá)都下調(diào)，說(shuō)明低溫脅迫后葉片葉綠素捕光能力下降。光合作用的原初反應(yīng)，即捕光色素分子吸收光能，將能量快速高效地傳遞給光合作用反應(yīng)中心，進(jìn)行電荷分離，固定激發(fā)能，經(jīng)過(guò)一系列電子傳遞步驟，將能量轉(zhuǎn)化為電化學(xué)自由能。PSⅡ反應(yīng)中心負(fù)責(zé)光能的轉(zhuǎn)換、電子和質(zhì)子的產(chǎn)生以及分子氧的釋放過(guò)程的啟動(dòng)。葉綠體PsbA和PsbD 基因編碼的D1和D2多肽是PSⅡ RC 的主要功能蛋白，它們結(jié)合著P680，脫鎂葉綠素a和 β-胡羅卜素，并且含有PSⅡ電子傳遞鏈中最重要的組分Z 和D，還為與水裂解相關(guān)的錳簇提供結(jié)合位點(diǎn)[31]。CP47 和CP43 與反應(yīng)中心的D1 和D2 蛋白緊密相連，除了將外周天線葉綠素a/b結(jié)合蛋白(LHCII、CP29、CP27、CP24、CP22)捕獲的激發(fā)能匯集給反應(yīng)中心外，還具有核心天線的作用。本研究發(fā)現(xiàn)，低溫脅迫后煙草CP47(PsbB)和CP43 (PsbC) 相關(guān)基因表達(dá)明顯上調(diào)，可能是為了適應(yīng)天線蛋白基因下調(diào)的結(jié)果，起到一種補(bǔ)償作用。光合作用的開(kāi)始是從PSⅠ開(kāi)始。PSⅠ是電子傳遞鏈過(guò)程，從初級(jí)接受者(Primary acceptor)開(kāi)始，經(jīng)過(guò)鐵氧化還原蛋白(Fd)、細(xì)胞色素復(fù)合體(Cytochrome Complex)、質(zhì)體藍(lán)蛋白(含銅蛋白質(zhì))(Pc)后再回到PSⅠ，本文中PSⅠ相關(guān)調(diào)控差異基因均表現(xiàn)為下調(diào)表達(dá)，說(shuō)明低溫影響電子鏈的傳遞。光合電子傳遞是由反應(yīng)中心啟動(dòng)，由電荷分離所導(dǎo)致的，由PSⅡ和PSⅠ串聯(lián)起來(lái)進(jìn)行的，其中細(xì)胞色素b6/f復(fù)合物位于PSⅡ和PSⅠ兩者之間。本文中低溫導(dǎo)致烤煙細(xì)胞色素b6/f復(fù)合物上調(diào)，光電子傳遞部分基因上調(diào)，部分基因下調(diào)，這也可能是烤煙幼苗低溫適應(yīng)的結(jié)果。葉綠素含量高低變化是光合作用重要的影響因素，也會(huì)影響葉片對(duì)光能的吸收。低溫脅迫3d后烤煙K326多數(shù)葉綠素相關(guān)基因處于下調(diào)，同時(shí)，葉綠素含量較對(duì)照顯著下降，研究結(jié)果表明葉綠素總含量變化和作為兩種捕光天線葉綠素a/b色素蛋白復(fù)合體基因下調(diào)是相對(duì)一致的。同時(shí)，K326在低溫處理后，其凈光合速率、蒸騰速率、葉片氣孔導(dǎo)度等均顯著下降，說(shuō)明其光合能力大大下降，與基因表達(dá)譜分析結(jié)果一致。

谷胱甘肽是一種在生物體內(nèi)廣泛存在的活性三肽，具有作為藥物、機(jī)體異物的解毒劑、機(jī)體抗氧化劑、氨基酸運(yùn)輸物質(zhì)、輔酶等各種生化作用[32]。尤其是可自發(fā)地或在GPX(谷胱甘肽過(guò)氧化物酶)催化下來(lái)清除體內(nèi)的活性氧，也可通過(guò)谷胱甘肽硫轉(zhuǎn)移酶(GST)催化形成GSX來(lái)清除活性物質(zhì)以保護(hù)細(xì)胞，降低氧化損傷[33]。低溫處理K326后進(jìn)行基因差異表達(dá)分析結(jié)果顯示，葡萄糖- 6-磷酸脫氫酶、精胺合酶、葉綠體類囊體29kDa蛋白等蛋白質(zhì)有關(guān)的基因均處于下調(diào)表達(dá)，說(shuō)明有機(jī)體有一套強(qiáng)大的自我保護(hù)及調(diào)控機(jī)制來(lái)應(yīng)對(duì)環(huán)境脅迫。谷胱甘肽相關(guān)基因的上調(diào)及下調(diào)表達(dá)是其更充分發(fā)揮抵御氧化壓力及解毒功能的作用，這種調(diào)控與自我保護(hù)機(jī)制密切相關(guān)。本文通過(guò)測(cè)定K326煙草幼苗葉片中自由基生成速率發(fā)現(xiàn)，低溫處理3 d后自由生產(chǎn)速率下降，從差異基因表達(dá)情況來(lái)看，谷胱甘肽相關(guān)基因多處于上調(diào)(表3)，這和自由基生產(chǎn)速率下降一致。本文谷胱甘肽上調(diào)基因中，主要包含谷胱甘肽S轉(zhuǎn)移酶、多胺(亞精胺合酶，腐胺氨丙基轉(zhuǎn)移酶)及抗壞血酸過(guò)氧化物酶等類型，均與抵御逆境脅迫有關(guān)。谷胱甘肽和Vc含量在低溫處理后顯著上升，說(shuō)明K326烤煙幼苗通過(guò)相關(guān)基因上調(diào)，提高了其清除自由基能力。林元震[34]通過(guò)將克隆的葡萄糖- 6-磷酸脫氫酶(G6PDH)基因轉(zhuǎn)移到煙草中進(jìn)行低溫脅迫后發(fā)現(xiàn)，隨著脅迫時(shí)間延長(zhǎng)，轉(zhuǎn)化煙草的SOD、POD活性升高而MDA含量下降，說(shuō)明該基因在煙草中的表達(dá)有助于膜保護(hù)性酶SOD和POD活性提高及細(xì)胞膜穩(wěn)定，進(jìn)而提高了轉(zhuǎn)基因煙草的耐低溫能力[34]。通過(guò)低溫脅迫處理K326后，其CAT酶活性上升，SOD酶活性顯著上升，POD酶活性則表現(xiàn)為顯著下降，這可能與差異基因表達(dá)譜中G- 6-PD基因中部分上調(diào)、部分處于下調(diào)有關(guān)。

植物對(duì)低溫脅迫的響應(yīng)是植物體內(nèi)一個(gè)復(fù)雜的過(guò)程，主要包括感應(yīng)低溫信號(hào)，信號(hào)轉(zhuǎn)導(dǎo)，及引起轉(zhuǎn)錄因子與相關(guān)功能性基因的差異表達(dá)。植物主要通過(guò)氣孔關(guān)閉、滲透性物質(zhì)積累、活性氧清除以及對(duì)膜和蛋白結(jié)構(gòu)的保護(hù)等一系列反應(yīng)作為對(duì)低溫脅迫的耐受響應(yīng)。本研究差異表達(dá)基因涉及到烤煙幼苗對(duì)低溫脅迫響應(yīng)機(jī)制的各個(gè)方面。此外，還有大量功能未知基因，這些基因參與多種代謝途徑，功能表達(dá)多表現(xiàn)為上調(diào)，少數(shù)表達(dá)下調(diào)，部分呈現(xiàn)超過(guò)10倍的表達(dá)量，這些基因與烤煙幼苗在低溫脅迫下的生理生態(tài)反應(yīng)將是下一步研究的主要內(nèi)容之一。本研究也獲得了一批與耐寒相關(guān)的基因，尤其是抗氧化相關(guān)的超量表達(dá)基因，為耐寒基因克隆及分子機(jī)理研究奠定了基礎(chǔ)。

Reference:

[1] Liang H M, Xia Y, Du F, Zhang P J. Effect of low temperature stress on physiological process of kentucky bluegrass. Acta Agrestia Sinica, 2001, 9(4): 283- 286.

[2] Wang Y, Yang H F, Li S D. Studies on chilling injury and cold hardiness of horticultural crops. Acta Horticulturae Sinica, 1994, 21(3): 239- 244.

[3] Yan Q, Ma Y S, Shi J J, Wang Y L, Yang S H. Effect of low temperature stress on physiological process of three kinds of herbage grasses in the seedling stage. Journal of Qinghai University: Natural Science, 2007, 25(1): 54- 57.

[4] Strand M, Oquist G. Effects of frost hardening, dehardening and freezing stress on in vivo fluorescence of seedlings of Scots pine (PinsylvestrisL.). Plant Cell and Environment, 1988, 11(4): 231- 238.

[5] Jin Y, Yang Y H, Hua S J, Duan Y Q, Cheng X Y. Effects of low temperature stress on the protective enzymes and contents of nitrogen and carbon compounds of tobacco seedlings. Journal of Southwest China Normal University: Natural Science Edition, 2007, 32(3): 74- 79.

[6] Hart J J, Stemler A. High light-induced reduction and low light-enhanced recovery of photon yield in triazine-resistantBrassicanapusL. Plant Physiology, 1990, 94(3): 1301- 1307.

[7] Zhou J, Yang L F, Hao F G, You Y. Photosynthesis and chlorophyll-fluorescence ofmagnoliagrandifloraseedlings under low temperature stress. Acta Botanica Boreali-Occidentalia Sinica, 2009, 29(1): 136- 142.

[8] Wu X X, Chen J L, Zha D S. Effects of low temperature stress on chlorophyll fluorescence characteristics and excitation energy dissipation in eggplant seedling leaves. Plant Nutrition and Fertilizer Science, 2009, 15(1): 164- 169.

[9] Chen H H, Li D H. Plant Cold Hardness and Freezing Stress. Vol. 2. New York: Academy Press, 1982: 5- 15.

[10] Prasad T K. Role of catalase in inducing chilling tolerance in pre-emergent maize seedling. Plant Physiology, 1997, 114(4): 1369- 1376.

[11] Liang P, Pardee A B. Differential display of eukaryotic messenger RNA by means of the polymerase chain reaction. Science, 1992, 257(5072): 967- 971.

[12] Nikiforova V, Freitag J, Kempa S, Adamik M, Hesse H, Hoefgen R. Transcriptome analysis of sulfur depletion inArabidopsisthaliana: interlacing of biosynthetic pathways provides response specificity. The Plant Journal, 2003, 33(4): 633- 650.

[13] Wang R, Guegler K, LaBrie S T, Crawford N M. Genomic analysis of a nutrient response in Arabidopsis reveals diverse expression patterns and novel metabolic and potential regulatory genes induced by nitrate. The Plant Cell, 2000, 12(8): 1491- 1509.

[14] Zou Z W, Fang W P, Zhang D, Duan Y S, Li X H. Analysis of differential expression genes in cold-induced tea plant. Journal of Tea Science, 2008, 28(4): 249- 254.

[15] Li Y L, Liu J H, Zheng J R, Hu J G. Gene expression profile of sweet corn ears under heat stress. Acta Agronomica Sinica, 2013, 39(2): 269- 279.

[16] Zhu S Q, Ji B H, Chen M W, Chen Y H, Zhou R, Liang J S. Expression profiling of rice transcription factor genes under salt or low temperature stress. Bulletin of Science and Technology, 2010, 26(6): 844- 852.

[17] Audic S, Claverie J M. The significance of digital gene expression profiles. Genome Research, 1997, 7(10): 986- 995.

[18] Ansorge W J. Next-generation DNA sequencing techniques. New Biotechnology, 2009, 25(4): 195- 203.

[19] Fowler S, Thomashow M F. Arabidopsis transcriptome profiling indicates that multiple regulatory pathways are activated during cold acclimation in addition to the CBF cold response pathway. Plant Cell, 2002, 14(8): 1675- 1690.

[20] Kaplan F, Kopka J, Sung D Y, Zhao W, Popp M, Porat R, Guy C L. Transcript and metabolite profiling during cold acclimation of Arabidopsis reveals an intricate relationship of cold-regulated gene expression with modifications in metabolite content. The Plant Journal, 2007, 50(6): 967- 981.

[21] Li L L, Wang F R, Han J M, Dong J G. Screening of drought-sensitive mutant inArabidopsisthalianaand responses of drought-sensitive mutant to drought stress. Journal of Agricultural University of Hebei, 2011, 34(4): 35- 40.

[22] Seki M, Narusaka M, Abe H, Kasuga M, Yamaguchi-Shinozaki K, Carninci P, Hayashizaki Y, Shinozaki K. Monitoring the expression pattern of 1300 Arabidopsis genes under drought and cold stresses by using a full-length cDNA microarray. Plant Cell, 2001, 13(1): 61- 72.

[23] Yamaguchi T, Nakayama K, Hayashi T, Yazaki J, Kishi-moto N, Kikuchi S, Koike S. cDNA microarray analysis of rice anther genes under chilling stress at the microsporo-genesis stage revealed two genes with DNA transponson castaway in 5′-flanking region. Bioscience Biotechnology, Biochemmistry, 2004, 68(6): 1315- 1323.

[24] Cui C. Studies on physiological ecological response and gene difference expression of flue-cured tobacco seedlings in chilling stress [D]. Chongqing: Southwest University, 2012.

[25] Benjamini Y, Yekutieli D. The control of the false discovery rate in multiple testing under dependency. The Annals of Statistics, 2001, 29(4): 1165- 1188.

[26] Morrissy A S, Morin R D, Delaney A, Zeng T, McDonald H, Jones S, Zhao Y, Hirst M, Marra M A. Next-generation tag sequencing for cancer gene expression profiling. Genome Research, 2009, 19(10): 1825- 1835.

[27] Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, Yamanishi Y. KEGG for linking genomes to life and the environment. Nucleic Acids Research, 2008, 36(Database issue): D480- 484.

[28] Hao Q N, Zhou X A, Sha A H, Wang C, Zhou R, Chen S L. Identification of genes associated with nitrogen-use efficiency by genome-wide transcriptional analysis of two soybean genotypes. BMC Genomics, 2011, 12: 525- 525.

[29] Xiang T H, Wang L L, Pang J L. Cloning and characterization of a full-length cab gene encoding the lightharvesting chlorophyll a/b-binding proteins in rice (OryzaSativaL.). Acta Agronomica Sinica, 2005, 31(9): 1227- 1232.

[30] Green B R, Pichersky E, Kloppstech K. Chlorophyll a/b binding proteins: an extended family. Trends Biochemical Sciences, 1991, 16(5): 181- 186.

[31] Nanba O, Satoh K. Isolation of a photosystem II reaction center consisting of D- 1 and D- 2 polypeptides and cytochrome b- 559. Proceedings of the National Academy of Sciences of the United States of America, 1987, 84(1): 109- 112.

[32] Mao Z, Qiu J P. Main physiological function and metabolic modulation of glutathione in yeast. Journal of Microbiology, 2005, 25(1): 94- 96.

[33] Jelinsky S A, Samson L D. Global response of saccharomyces cerevisiaeto an alkylating agent. Proceedings of the National Academy of Sciences of the United States of America, 1999, 96(4): 1486- 1491.

[34] Lin Y Z, Zhang Z Y, Guo H, Liu C X, Chen X Y. Isolation and analysis of glucose-6-phosphate dehydrogenase (G6PDH) promoter from poplar. Genomics and Applied Biology, 2009, 28(3): 445- 449.

[1] 梁慧敏，夏陽(yáng)，杜峰，張普金.低溫脅迫對(duì)草地早熟禾抗性生理生化指標(biāo)的影響.2001,9(4):283- 286.

[3] 嚴(yán)青，馬玉壽，施建軍，王彥龍，楊時(shí)海.低溫脅迫對(duì)3種牧草幼苗抗性生理指標(biāo)的影響.青海大學(xué)學(xué)報(bào)(自然科學(xué)版),2007,25(1):54- 57.

[5] 晉艷,楊宇虹,華水金,段玉琪,程新宇.低溫脅迫對(duì)煙草保護(hù)性酶類及氮和碳化合物的影響.西南師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2007,32(3): 74- 79.

[7] 周建，楊立峰，郝峰鴿，尤揚(yáng).低溫脅迫對(duì)廣玉蘭幼苗光合及葉綠素?zé)晒馓匦缘挠绊?西北植物學(xué)報(bào)，2009,29(1)：136- 42.

[8] 吳雪霞, 陳建林, 查丁石.低溫脅迫對(duì)茄子幼苗葉片葉綠素?zé)晒馓匦院湍芰亢纳⒌挠绊?植物營(yíng)養(yǎng)與肥料學(xué)報(bào)，2009,15(1): 164- 169.

[14] 鄒中偉, 房婉萍, 張定, 段云裳, 黎星輝. 低溫脅迫下茶樹(shù)基因表達(dá)的差異分析. 茶葉科學(xué), 2008, 28(4): 249- 254.

[15] 李余良, 劉建華, 鄭錦榮, 胡建廣. 高溫脅迫下甜玉米雌穗發(fā)育基因差異表達(dá)譜分析. 作物學(xué)報(bào), 2013, 39(2): 269- 279.

[16] 朱素琴, 季本華, 陳名蔚, 陳艷紅, 周蓉, 梁建生. 高鹽低溫脅迫下水稻轉(zhuǎn)錄因子基因表達(dá)譜分析. 科技通報(bào), 2010, 26(6): 844- 852.

[21] 李琳琳, 王鳳茹, 韓建民, 董金皋. 擬南芥干旱敏感突變體的篩選及其對(duì)干旱脅迫的響應(yīng). 河北農(nóng)業(yè)大學(xué)學(xué)報(bào), 2011, 34(4): 35- 40.

[24] 崔翠. 低溫脅迫下烤煙幼苗的生理生態(tài)響應(yīng)及基因差異表達(dá)研究 [D]. 重慶: 西南大學(xué), 2012.

[29] 向太和, 王利琳, 龐基良. 水稻(OryzasativaL.)捕光葉綠素a/b結(jié)合蛋白基因全長(zhǎng)cDNA的克隆和特性分析. 作物學(xué)報(bào), 2005, 31(9): 1227- 1232.

[32] 毛珍, 裘娟萍. 酵母菌中谷胱甘肽的主要生理功能及其代謝調(diào)控. 微生物學(xué)雜志, 2005, 25(1): 94- 96.

[34] 林元震, 張志毅, 郭海, 劉純鑫, 陳曉陽(yáng). 楊樹(shù)葡萄糖- 6- 磷酸脫氫酶(G6PDH)基因啟動(dòng)子的克隆與分析. 基因組學(xué)與應(yīng)用生物學(xué), 2009, 28(3): 445- 449.

Expression profiling of genes related to photosynthesis and antioxidant capacity in flue-cured tobacco seedlings subjected to chilling stress

CUI Cui*, WANG Lijuan, ZHOU Qingyuan, TAN Zunfei, QU Cunmin, ZHANG Zhengsheng

CollegeofAgronomyandBiotechnology,SouthwestUniversity,Chongqing400715,China

Tobacco is an important economical leaf crop and complex model organism that is cultivated worldwide. Low temperature is the one of major factors causing abiotic damage to flue-cured seedlings; this damage can affect the growth and development of flue-cured tobacco seedlings and can decrease the yield and quality of flue-cured tobacco leaves. Currently, technologies for cultivating seedlings by using water and oxygen have been widely adopted in southern tobacco-cultivating regions. However, temperature is dependent on climate, which is difficult to control artificially. Low temperatures during the early spring season limit the culture of strong seedlings in southern tobacco-cultivation areas. The objective of this study was to analyze physiological and ecological adaptations of flue-cured tobacco seedlings by measuring parameters related to levels of photosynthetic, oxidant, and antioxidant factors in the cell membranes of leaves. We also used digital gene expression profiling technology to analyze the differential expression of genes in flue-cured tobacco seedlings after chilling stress. Seedlings of flue-cured tobacco strain K326 were used as experiment materials. Tobacco seedlings with 5—6 true leaves were divided in two groups. One group of flue-cured seedlings was placed in a light incubator at a low temperature range of 5 ℃ (night) to 7 ℃ (day) for 3 days as chilling stress treatment. Seedlings in the other (control) group remained at the initial temperature and illumination conditions and grew at a normal temperature range of 23 ℃ (night) to 25 ℃ (day), with other similar conditions. After chilling stress for 3 days, total chlorophyll (Chl), chlorophyll a (Chl a), and chlorophyll b (Chl b) contents decreased significantly by 2.88%—6.82%, and the net photosynthetic rate (Pn), stomatal conductance (Gs), and transpiration (Tr) decreased significantly by 14.14%—68.50% compared with that of seedlings grown under the favorable temperature. Intercellular CO2concentration (Ci) showed no obvious changes. In terms of the membrane oxidant levels of leaves from flue-cured tobacco seedlings, proline and malondialdehyde (MDA) contents and electrolyte permeability increased significantly by 3.88%—144.22% compared with that of control seedlings, while the oxyradical generation rate decreased significantly by 10.66%. Additionally, the membrane antioxidant capacity of flue-cured tobacco seedlings and the activities of superoxide dismutase (SOD) and catalase (CAT) increased significantly by 6.07% and 45.64%, respectively, compared with that grown under the favorable temperature, while peroxidase (POD) activity was significantly reduced. The vitamin C (Vc) and glutathione (GSH) contents increased significantly by 197.36% and 14.15%, respectively, compared with that of the control. Under chilling stress conditions, 2357 genes from K326 seedlings were differentially expressed significantly compared with that of seedlings grown under the favorable temperature. Of these, 1673 genes were upregulated, and 684 genes downregulated. Gene ontology analysis revealed that these differentially expressed genes were mainly involved in transcription factors, transmembrane transporter proteins, antioxidant systems, signaling pathways, and others. Gene ontology analysis revealed relationships with the photosystem, photosynthetic membrane, chloroplast, plastid, cytoplasm, membrane, ribosomes, and more. In particular, 524 differentially expressed genes were associated with plastid expression. The molecular functions of these differentially expressed genes were mainly associated with oxidoreductase activity and the antioxidant system of the cell membrane of flue-cured seedlings. Differential expression was also observed for genes involved in biological processes such as response to stress, photosynthesis, light reaction of photosynthesis, hormone-mediated signaling pathways, and more. After analysis of metabolism pathways related to differential gene expression, our results showed that genes related to photosynthesis-antenna proteins and photosynthesis were downregulated, while genes involved in oxidant capacity, such as Vc, GSH, and proline metabolism, were upregulated. Leaves of flue-cured tobacco seedlings suffered oxidant damage under chilling stress condition, but the promotion of antioxidant ability by chilling stress had an active protective effect on flue-cured tobacco seedlings. Therefore, flue-cured seedlings had some capacity to adapt to chilling stress by adjustment of gene expression, which altered pathways involved in photosynthetic, oxidant, and antioxidant metabolism to reduce the damage resulting from the stress conditions. Importantly, the results of gene expression analysis and physiological ecological adaptation were consistent. Hence, it is possible to analyze ecological adaptation and differential gene expression of crops under stress, and such studies will facilitate further analysis of gene ontology functions and gene expression relationships.

flue-cured tobacco; chilling stress; digital gene expression profile; differentially expressed genes

高等學(xué)校學(xué)科創(chuàng)新引智計(jì)劃(B12006); 重慶市自然科學(xué)基金項(xiàng)目(cst2011jjA80023); 西南大學(xué)博士啟動(dòng)基金(SWU112074); 中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)(2362014xk09, XDJK2014C150)

2013- 10- 24;

2014- 07- 03

10.5846/stxb201310242565

*通訊作者Corresponding author.E-mail: cuigreeny@163.com

崔翠，王利鵑，周清元，譚尊飛，曲存民，張正圣.低溫脅迫下烤煙幼苗葉片光合作用和抗氧化能力基因差異表達(dá)譜.生態(tài)學(xué)報(bào),2014,34(21):6076- 6089.

Cui C, Wang L J, Zhou Q Y, Tan Z F, Qu C M, Zhang Z S.Expression profiling of genes related to photosynthesis and antioxidant capacity in flue-cured tobacco seedlings subjected to chilling stress.Acta Ecologica Sinica,2014,34(21):6076- 6089.