白明興,莊澤龍,姬祥卓,彭云玲

(甘肅農(nóng)業(yè)大學(xué)農(nóng)學(xué)院,甘肅省干旱生境作物學(xué)重點實驗室/甘肅省作物遺傳改良與種質(zhì)創(chuàng)新重點實驗室, 甘肅 蘭州 730070)

干旱作為全球最重要的非生物脅迫之一，極大地影響了作物的生長發(fā)育和產(chǎn)量。玉米作為我國主要的糧飼兼用作物，是典型的干旱敏感作物，各生育期需水量較大[1-3]。研究表明，干旱缺水不僅嚴(yán)重影響玉米的生長發(fā)育，造成玉米大幅減產(chǎn)，而且會阻礙玉米葉片對CO2的吸收，進而影響其光合性能[4-5]。此外，干旱脅迫會導(dǎo)致植物體內(nèi)活性氧(ROS)不斷積累，ROS在植物生長發(fā)育中有著十分重要的作用，正常條件下主要參與調(diào)節(jié)植物的生長發(fā)育，維持動態(tài)平衡狀態(tài)[6]；當(dāng)遭受非生物脅迫時，植物體內(nèi)的ROS快速積累，多余的ROS不能及時清除，動態(tài)平衡遭到破壞，從而影響植物的發(fā)育[7]。活性氧是氧分子沒有完全還原的產(chǎn)物，其毒性很強，在非生物脅迫下，它不僅可以破壞細(xì)胞的結(jié)構(gòu)和功能，而且是信號轉(zhuǎn)導(dǎo)的重要調(diào)節(jié)因子[8]。由于ROS具有極高的氧化活性，植物體內(nèi)蛋白質(zhì)、磷脂和核酸很容易被氧化，造成細(xì)胞中蛋白質(zhì)降解和核酸裂解，以及磷脂過氧化產(chǎn)生丙二醛，破壞細(xì)胞膜結(jié)構(gòu)[9]。且大量研究也表明，ROS 作為信號分子在植物細(xì)胞防御響應(yīng)和細(xì)胞程序性死亡等方面扮演著極其重要的角色[10]。當(dāng)植物遭受輕度干旱脅迫時，會導(dǎo)致植物體內(nèi)活性氧清除啟動一系列復(fù)雜的清除反應(yīng)[11]，促使超氧化物歧化酶(SOD)、過氧化氫酶(CAT)和過氧化物酶(POD)等酶活力提高，保護植物免受環(huán)境脅迫的危害。根系作為植物吸收和運輸土壤中水分及養(yǎng)料的主要器官，在遭遇非生物脅迫時，根系一般最先感受到逆境脅迫，從而改變根系的形態(tài)以及分布以適應(yīng)不利環(huán)境。因此，深入剖析玉米苗期根系對干旱脅迫的分子基礎(chǔ)及響應(yīng)機制，對改善干旱和半干旱地區(qū)作物的早期生長狀況及促進全球糧食生產(chǎn)安全有著極其重要的作用[12]。

近年來，隨著全基因組測序的快速發(fā)展，不同學(xué)者在研究植物抗逆性方面也逐漸采用轉(zhuǎn)錄組和蛋白質(zhì)組學(xué)技術(shù)來揭示逆境脅迫中植物的生理、生化和分析響應(yīng)機制，從而挖掘大量與抗逆相關(guān)的基因。ZENDA等[13]為了揭示不同玉米自交系對干旱脅迫的耐受性，采用耐旱自交系YE8112和干旱敏感自交系MO17進行比較轉(zhuǎn)錄組分析，通過分析獲得了一些耐旱基因和調(diào)控干旱脅迫的主要途徑，闡述了不同玉米自交系對干旱脅迫的響應(yīng)機制，并依據(jù)目前研究構(gòu)建了玉米幼苗抗旱的分子模型。LIU等[14]通過對兩個雜交品種(耐旱ND476和旱敏ZX978)在干旱脅迫下進行轉(zhuǎn)錄組分析發(fā)現(xiàn)，玉米對干旱脅迫的適應(yīng)是一種品種特異性的反應(yīng)，也是一種階段性的反應(yīng)過程，并且該研究也挖掘了一些調(diào)控玉米干旱脅迫反應(yīng)的品種特異性基因和生育期特異性基因。CHEN等[15]通過對兩份耐鹽性不同的玉米自交系根系進行蛋白質(zhì)組學(xué)分析發(fā)現(xiàn)，當(dāng)遭遇鹽脅迫時，不同耐鹽性自交系差異表達蛋白(DEPs)間雖存在明顯差異，但是不同自交系間對鹽脅迫的響應(yīng)也存在一些相同的趨勢。JIANG等[16]利用iTRAQ(Isobaric tags for relative and absolute quantification)技術(shù)發(fā)現(xiàn)玉米幼苗在干旱脅迫下存在多重轉(zhuǎn)錄后調(diào)控和翻譯后修飾，并提出在干旱脅迫下玉米幼苗會通過提高其抗旱性來應(yīng)對干旱的策略。ZHENG等[17]采用耐旱材料28M和旱敏感材料753F從基因調(diào)控水平揭示了287M抗旱性的機制，表明287M和753F的抗旱性差異可能與清除ROS能力、信號相互作用網(wǎng)絡(luò)和某些轉(zhuǎn)錄因子有關(guān)。HAO等[18]利用RNA測序技術(shù)對耐旱性不同的兩自交系根系干旱相關(guān)基因進行全基因組鑒定和比較分析，并表明玉米的抗旱性與ABA信號通路、海藻糖生物合成、活性氧清除和轉(zhuǎn)錄因子相關(guān)基因的表達模式有關(guān)。以上研究試圖從某一組學(xué)水平去解析一些生物學(xué)現(xiàn)象，不能深入剖析其響應(yīng)機理。例如，轉(zhuǎn)錄組只能關(guān)注一些關(guān)鍵編碼轉(zhuǎn)錄本的表達量變化，而不能了解其轉(zhuǎn)錄調(diào)控以及對蛋白翻譯的影響。因此，利用多組學(xué)聯(lián)合分析剖析一些生物學(xué)問題十分必要。轉(zhuǎn)錄組和蛋白質(zhì)組學(xué)可以從兩個不同層面反映基因的表達情況，所以要了解轉(zhuǎn)錄組與蛋白質(zhì)組的相互調(diào)控作用，需要對mRNA 與蛋白質(zhì)的表達進行同步監(jiān)測。關(guān)聯(lián)分析的目的之一是實現(xiàn)數(shù)據(jù)互補，得到生物體更加完整的表達信息。

因此，為了進一步剖析玉米幼苗對干旱脅迫的響應(yīng)機制，本研究以耐旱自交系昌7-2和干旱敏感自交系TS141的苗期根系為試驗材料，結(jié)合轉(zhuǎn)錄組測序和iTRAQ技術(shù)探究不同玉米自交系對干旱的響應(yīng)機制并篩選一些耐旱候選基因，為進一步了解玉米耐旱的分子機制提供新的思路。

1 材料與方法

1.1 材料處理

用耐旱自交系昌7-2(Chang 7-2)和干旱敏感自交系TS141為試材，對兩自交系分別進行蒸餾水(CK)和20%PEG(D)處理。挑選整齊、一致的種子用0.5%NaClO表面滅菌10 min并用無菌蒸餾水洗滌5次。將種子浸泡在處理溶液中24 h后點播在營養(yǎng)缽(15 cm×13 cm)中，消毒蛭石(溶液∶蛭石=1∶5 v·m-1)用與各處理相同的處理溶液攪拌均勻，每缽10粒，試驗設(shè)置3次重復(fù)。將材料置于人工氣候箱，晝夜溫度分別為25℃、20℃，光照和黑暗時間分別為16、8 h，光照強度為600 μmol·s-1·m-2。對相應(yīng)材料每隔1 d加入50 mL蒸餾水，每3 d施加50 mL對應(yīng)處理液，處理12 d后收集根系樣品進行轉(zhuǎn)錄組測序和蛋白質(zhì)組學(xué)測序。

1.2 轉(zhuǎn)錄組測序和蛋白質(zhì)組測序分析

轉(zhuǎn)錄組測序參照CHEN等[19]的方法，對兩不同玉米自交系處理后12 d的幼苗根系進行總RNA提取，并進行轉(zhuǎn)錄組測序。測序所獲得的原始讀段(Raw reads)進行過濾，得到干凈讀段(Clean reads)，后續(xù)分析都基于Clean reads，并采用HISAT和RSEM軟件以第三版B73玉米參考基因組對Clean reads進行比對(Version ZmB73_5a.59)(http://ftp.maizesequence.org/)和基因表達水平檢測。通過RPKM計算12個樣品的基因表達(每百萬Kbp讀數(shù))方法，其測序的數(shù)據(jù)已提交到NCBI的SRA(登錄號：PRJNA698088)。為了提高DEG數(shù)據(jù)的準(zhǔn)確性，差異基因篩選標(biāo)準(zhǔn)采用|log2FC(Fold change)|≥1且FDR≤0.001。而蛋白質(zhì)組測序的具體方法及相關(guān)數(shù)據(jù)來自ZENG等[20]。采用WEGO軟件[21]對篩選出的DEGs在GO數(shù)據(jù)庫中的功能注釋進行分類。通路顯著性富集分析基于KEGG(Kyoto encyclopedia of genes and genomes)公共數(shù)據(jù)庫(http:/iwww.genome.jp/kegg/pathway.html)進行，對顯著富集的通路進行鑒定。其中把Q value≤0.05的Pathway定義為在DEGs中顯著富集的Pathway。

1.3 關(guān)聯(lián)分析

干旱脅迫下，利用LAN等[22]的方法計算了昌7-2和TS141蛋白質(zhì)水平與轉(zhuǎn)錄水平的Pearson相關(guān)系數(shù)，確定蛋白質(zhì)與mRNA的相關(guān)性，并對每個品種關(guān)聯(lián)了兩個生物學(xué)重復(fù)，其中Chang 7-2 D1vsCK1比較組表示耐旱自交系昌7-2轉(zhuǎn)錄組和蛋白組第一個生物學(xué)重復(fù)所關(guān)聯(lián)的結(jié)果，Chang 7-2 D2vsCK2比較組表示耐旱自交系昌7-2轉(zhuǎn)錄組和蛋白組第二個生物學(xué)重復(fù)所關(guān)聯(lián)的結(jié)果；而TS141 D1vsCK1和TS141 D2vsCK2比較組分別表示旱敏感自交系TS141中兩個重復(fù)的關(guān)聯(lián)結(jié)果。

1.4 qRT-PCR(Quantitative real time-PCR)驗證

采用RNA提取試劑盒獲取兩自交系不同處理下幼苗根系的總RNA，隨后利用FastKing RT Kit(天根生物技術(shù)，北京，中國)對提取的總RNA反轉(zhuǎn)錄為cDNA，備用。利用NCBI(https://www.ncbi.nlm.nih.gov/)官網(wǎng)進行引物設(shè)計，本試驗根據(jù)不同基因的表達量共設(shè)計了9個基因的特異性PCR引物對轉(zhuǎn)錄組數(shù)據(jù)和蛋白數(shù)據(jù)進行驗證，實時熒光定量(qRT-PCR)分析的具體方法參考CHEN等[19]，所選基因的相對表達水平通過2-ΔΔCT方法計算，并以Actin基因的表達水平進行歸一化，每個處理使用3個生物重復(fù)。

2 結(jié)果與分析

2.1 轉(zhuǎn)錄組和蛋白質(zhì)組數(shù)據(jù)分析

2.1.1 轉(zhuǎn)錄組測序分析 采用BGISEQ-500測序平臺對兩份玉米自交系根系樣品進行RNA-Seq測序。結(jié)果表明測序樣品平均產(chǎn)生116025890.2條Raw reads，經(jīng)過剔除低質(zhì)量的Reads后, 每個樣品平均得到107905100.5條高質(zhì)量的Clean reads。采用HISAT軟件將獲得的Clean reads與玉米參考基因組B73(http://ftp.maizesequence.org/)進行比對，結(jié)果表明不同樣品在參考基因組中的平均比對率約68.77%。此外，以|log2Ratio(FC, Fold Change)|≥1,Pvalue≤0.001作為兩樣品間差異表達基因(DEGs)的篩選條件對不同樣品間進行差異表達基因篩選，結(jié)果表明，在昌7-2中鑒定出7483個DEGs，其中3162個DEGs上調(diào)表達，4321個DEGs下調(diào)表達；在TS141中鑒定出6829個DEGs，其中3290個DEGs上調(diào)表達，3539個DEGs下調(diào)表達(圖1)。

2.1.2 蛋白質(zhì)組數(shù)據(jù)分析 采用iTRAQ技術(shù)對抗旱性不同的兩自交系幼苗根系進行蛋白質(zhì)組學(xué)分析，結(jié)果共產(chǎn)生了97374張光譜、39371個肽以及30148個獨特的肽，最終鑒定到7723個蛋白質(zhì)。當(dāng)?shù)鞍踪|(zhì)豐度比達到0.83倍以下或1.2倍以上，且P<0.05時，認(rèn)為該蛋白是兩個不同樣品間的差異表達蛋白(DEPs)。結(jié)果表明，在Chang 7-2-DvsCK中鑒定出1243個DEPs，其中572個DEPs上調(diào)表達，671個DEPs下調(diào)表達；TS141-DvsCK中鑒定出419個DEPs，其中172個DEPs上調(diào)表達，247個DEPs下調(diào)表達(圖1)。

圖1 干旱處理下兩玉米自交系的DEGs和DEPs

2.2 轉(zhuǎn)錄組和蛋白質(zhì)組學(xué)關(guān)聯(lián)性分析

2.2.1 轉(zhuǎn)錄組和蛋白質(zhì)組差異表達基因與差異表達蛋白的關(guān)聯(lián)分析 對轉(zhuǎn)錄組和蛋白質(zhì)組數(shù)據(jù)進行關(guān)聯(lián)分析，當(dāng)某一個蛋白質(zhì)在轉(zhuǎn)錄組水平有表達量時，被認(rèn)為關(guān)聯(lián)到。關(guān)聯(lián)結(jié)果表明(表1)，在鑒定、定量和顯著差異表達3個方面，Chang 7-2 D1vsCK1中7612個基因在mRNA和蛋白質(zhì)水平被同時鑒定，并定量到7609個基因，有157個DEGs 同DEPs關(guān)聯(lián)。Chang 7-2 D2vsCK2中7621個基因在mRNA和蛋白質(zhì)水平被同時鑒定，并定量到7619個基因，有98個DEGs 同DEPs關(guān)聯(lián)。TS141 D1vsCK1中7609個基因在mRNA和蛋白質(zhì)水平被同時鑒定，并定量到7606個基因，有76個DEGs 同DEPs關(guān)聯(lián)。TS141 D2vsCK2中7619個基因在mRNA和蛋白質(zhì)水平被同時鑒定，并定量到7616個基因，有116個DEGs 同DEPs關(guān)聯(lián)。

表1 鑒定、定量和顯著差異中被關(guān)聯(lián)的基因和蛋白質(zhì)數(shù)量

2.2.2 轉(zhuǎn)錄組和蛋白質(zhì)組相關(guān)性分析 基于mRNA水平和蛋白質(zhì)水平的表達結(jié)果，使關(guān)聯(lián)到的基因在mRNA表達水平和蛋白質(zhì)水平上存在多種類型，即mRNA和蛋白表達趨勢相同、mRNA和蛋白表達趨勢相反、mRNA表達有差異而蛋白表達無差異、mRNA表達無差異而蛋白表達有差異以及mRNA和蛋白表達均無差異5種關(guān)聯(lián)類型。本研究發(fā)現(xiàn)，Chang 7-2 D1vsCK1、Chang 7-2 D2vsCK2、TS141 D1vsCK1和TS141 D2vsCK2等4個比較組中所關(guān)聯(lián)到DEGs和DEPs的相關(guān)性系數(shù)分別為0.3063、0.2208、0.2295和-0.0203(圖2)。其中，mRNA和蛋白表達趨勢相同的基因分別為98(重復(fù)基因2個)、64、55(重復(fù)基因2個)、62個(重復(fù)基因1個)(表2)，去除組內(nèi)重復(fù)基因后分別為96、64、53、61；而mRNA和蛋白表達趨勢相反的基因分別為59(重復(fù)基因1個)、34(重復(fù)基因2個)、21(重復(fù)基因1個)、54個(表2), 去除組內(nèi)重復(fù)基因后分別為58、32、20、54。

圖2 顯著差異蛋白質(zhì)和顯著差異基因表達關(guān)聯(lián)圖

為了進一步剖析兩自交系在干旱脅迫下基因表達的調(diào)控情況，對每個自交系D1vsCK1和D2vsCK2的關(guān)聯(lián)結(jié)果進行了整合，以獲得在干旱脅迫下穩(wěn)定表達的基因信息。結(jié)果表明，Chang 7-2的D1vsCK1和D2vsCK2樣品中mRNA和蛋白表達趨勢相同的基因有19個(表3；圖3A)，mRNA和蛋白表達趨勢相反的基因有9個(表3；圖3B)；TS141 D1vsCK1和D2vsCK2中mRNA和蛋白表達趨勢相同的基因有11個(表3)，mRNA和蛋白表達趨勢相反的基因在兩次重復(fù)中沒檢索到較為穩(wěn)定表達的基因(圖3B)。

圖3 兩玉米自交系在不同處理下DEGs和DEPs的Venn圖

2.2.3 GO富集關(guān)聯(lián)分析 對兩個組學(xué)中差異表達蛋白和基因進行GO功能注釋分析，結(jié)果表明Chang 7-2 D1vsCK1、Chang 7-2 D2vsCK2、TS141 D1vsCK1和TS141 D2vsCK2等4個比較組的GO功能富集總體分為細(xì)胞成分(Cellular component, C)、分子功能(Molecular function, F)和生物過程(Biological process , P)，而干旱處理下，比較組Chang 7-2 D1vsCK1和Chang 7-2 D2vsCK2所關(guān)聯(lián)DEGs/DEPs的主要細(xì)胞組分包括細(xì)胞內(nèi)非膜-b(Intracellular non-membrane-b)、非膜結(jié)合細(xì)胞器(Non-membrane-bounded organel)、細(xì)胞外圍(Cell periphery)；生物學(xué)過程主要有氧化-還原過程(Oxidation-reduction process)、脅迫反應(yīng)(Response to stress)、碳水化合物代謝過程(Carbohydrate metabolic proce)；分子功能主要包括金屬離子結(jié)合(Metal ion binding)、陽離子結(jié)合(Cation binding)、氧化還原酶活性(Oxidoreductase activity)。比較組TS141 D1vsCK1和TS141 D2vsCK2所關(guān)聯(lián)到的DEGs/DEPs主要細(xì)胞組分包括Intracellular non-membrane-b、Non-membrane-bounded organel和核糖核蛋白復(fù)合物(Ribonucleoprotein complex)；分子功能主要包括Oxidoreductase activity和Metal ion binding；其生物學(xué)過程主要有脅迫反應(yīng)(Response to stress)和刺激反應(yīng)(Response to stimulus)。此外，我們還發(fā)現(xiàn)兩玉米自交系幼苗根系所關(guān)聯(lián)到的DEGs/DEPs均富集到脅迫響應(yīng)(Response to stress)過程，即遭遇干旱脅迫后，其幼苗根系會立即對干旱脅迫做出響應(yīng)，并通過促進根系內(nèi)過氧化物酶活性、氧化還原酶活性等多種酶活性來抵御干旱。

2.2.4 Pathway富集關(guān)聯(lián)分析 利用KEGG pathway數(shù)據(jù)庫對Chang 7-2 D1vsCK1、Chang 7-2 D2vsCK2、TS141 D1vsCK1和TS141 D2vsCK2等4個比較組的轉(zhuǎn)錄組和蛋白質(zhì)組數(shù)據(jù)以及所關(guān)聯(lián)到的DEGs/DEPs進行通路富集分析，結(jié)果表明，差異基因和差異蛋白在各代謝通路中富集顯著性分為3大類：轉(zhuǎn)錄組與蛋白質(zhì)組都顯著富集、轉(zhuǎn)錄組與蛋白質(zhì)組任意一個組學(xué)數(shù)據(jù)顯著富集和轉(zhuǎn)錄組與蛋白質(zhì)組均無顯著富集。而我們僅對轉(zhuǎn)錄組與蛋白組兩個組學(xué)都顯著富集的通路進行了分析，發(fā)現(xiàn)4個比較組所關(guān)聯(lián)到的DEGs/DEPs主要富集的通路包括苯丙素的生物合成(Phenylpropanoid biosynthesis)、次生代謝產(chǎn)物的生物合成(Biosynthesis of secondary metabolites)和核糖體(Ribosome)(表4)。由此發(fā)現(xiàn)，玉米幼苗在應(yīng)對干旱脅迫時會發(fā)生一系列復(fù)雜的代謝活動，而苯丙素的生物合成、次生代謝產(chǎn)物的生物合成和核糖體等相關(guān)代謝通路可能是根系響應(yīng)干旱脅迫的主要通路。

2.2.5 qRT-PCR驗證 為了驗證干旱脅迫下兩玉米自交系轉(zhuǎn)錄組和蛋白質(zhì)組數(shù)據(jù)對干旱脅迫響應(yīng)的可靠性，我們通過qRT-PCR驗證了在干旱處理下昌7-2和TS141中表達水平不同的9個基因(表5、圖4)。我們發(fā)現(xiàn)這9個基因在兩個自交系中的基因表達水平與RNA-Seq和iTRAQ結(jié)果呈現(xiàn)相同趨勢，這表明兩組學(xué)測序數(shù)據(jù)是可靠的。

圖4 qRT-PCR驗證所選基因的相對表達量

表5 用于qRT-PCR分析的基因

3 討 論

為了進一步剖析不同玉米自交系對干旱脅迫的響應(yīng)機制，采用耐旱自交系昌7-2和干旱敏感自交系TS141作為試驗材料，在干旱處理12 d后對其幼苗根系進行轉(zhuǎn)錄組學(xué)分析和蛋白質(zhì)組學(xué)分析，并對其根系差異基因與差異蛋白組表達趨勢相同的部分基因進行了進一步分析。

3.1 耐旱玉米自交系昌7-2根系對干旱的響應(yīng)機制

Jasmonate O-methyltransferase作為催化茉莉酸合成MeJA最重要的限速步驟[23]，是MeJA合成途徑遺傳改造的最重要位置。研究表明，茉莉酸類化合物作為信號分子參與調(diào)控植物的生長發(fā)育與代謝途徑，響應(yīng)低溫、干旱和鹽等多種非生物脅迫[24]。此外，茉莉酸類物質(zhì)在植物對逆境防御及脅迫傳遞信號等方面有著極其重要的作用。MA等[25]和WU等[26]研究結(jié)果也表明，茉莉酸甲酯緩解物質(zhì)通過調(diào)節(jié)小麥(Triticumaestivum)和甘藍(Brassicaoleracea)氣孔運動及相關(guān)基因表達來調(diào)節(jié)干旱缺水對作物所帶來的損傷，提高植物的抗旱性。本研究在昌7-2中也發(fā)現(xiàn)一個編碼Jasmonate O-methyltransferase的基因(103629478)在轉(zhuǎn)錄組和蛋白質(zhì)組中均上調(diào)表達，推測其在干旱脅迫下通過調(diào)節(jié)幼苗根系的氣孔開度抗氧化酶活性來抵御干旱。

幾丁質(zhì)酶(Endochitinase)是一種能夠?qū)锥≠|(zhì)分解成幾丁質(zhì)單糖—N-乙?；咸烟堑拿?，幾乎在植物各器官中都有分布。有研究表明，當(dāng)植物遭受干旱、低溫、凍害及重金屬離子等多重脅迫時，其組織內(nèi)幾丁質(zhì)酶的活性會快速升高并水解病原菌細(xì)胞壁的幾丁質(zhì)，抑制病原菌在植物體內(nèi)的生長繁殖[27]。近年來，關(guān)于幾丁質(zhì)酶基因的研究在不同植物中也不斷被報道，陳鵬等[28]通過對小麥進行干旱脅迫后發(fā)現(xiàn)，隨著干旱脅迫程度的增加，其體內(nèi)幾丁質(zhì)酶的活性也在不斷上升，復(fù)水處理后又逐漸恢復(fù)至正常水平，表明在干旱脅迫下幾丁質(zhì)酶可以提高植物抗旱性。MUOKI等[29]通過對茶樹進行環(huán)境脅迫后發(fā)現(xiàn)，幾丁質(zhì)酶基因Chi與Thaumatin蛋白可共同響應(yīng)干旱脅迫，進而促進茶樹的抗旱性，減弱干旱對茶樹的損害。在本研究中，我們也發(fā)現(xiàn)當(dāng)昌7-2在遭遇干旱脅迫后其根系內(nèi)幾丁質(zhì)酶相關(guān)基因出現(xiàn)差異上調(diào)表達。

此外，當(dāng)昌7-2遭受干旱脅迫時，其體內(nèi)與熱激蛋白(100283821,Hsp20)、滲透蛋白(100276695,Osmotin-like protein OSM34)及過氧化物酶(100285894、100273479、541674、107548102)相關(guān)基因及蛋白也發(fā)生上調(diào)或下調(diào)表達。研究表明，熱激蛋白作為植物響應(yīng)非生物脅迫的重要蛋白之一，在植物耐寒耐旱、保護植物相關(guān)蛋白完成正確的折疊、翻譯和聚集及結(jié)構(gòu)功能方面發(fā)揮著重要作用[30]，目前已發(fā)現(xiàn)多種類型熱激蛋白，例HSP60、HSP70和HSP90等家族[31-32]。另外，干旱脅迫下，滲透調(diào)節(jié)蛋白基因和過氧化物酶(POD)對植物應(yīng)對逆境脅迫有非常重要的作用。Osm基因(即ProBA基因)是一種可增強植物細(xì)胞抗?jié)B透脅迫能力的蛋白質(zhì)，其主要功能為在逆境脅迫下促進脯氨酸在細(xì)胞內(nèi)的合成和累積[33]，而在逆境脅迫下不同類型的POD相互作用可有效清除植物體內(nèi)活性氧。因此，我們推測昌7-2受到干旱脅迫后，其幼苗根系會通過調(diào)節(jié)體內(nèi)熱激蛋白和滲透蛋白相關(guān)基因上調(diào)，同時通過調(diào)控不同類型的過氧化物酶上或下調(diào)來減輕干旱所帶來的損傷。

3.2 玉米干旱敏感自交系TS141根系對干旱的響應(yīng)機制

當(dāng)對TS141進行干旱處理后，我們發(fā)現(xiàn)其mRNA和蛋白表達趨勢相同的基因主要編碼6-磷酸葡萄糖酸內(nèi)酯酶(100285843)、核苷二磷酸激酶I(100217209)、脯氨酸-4-羥化酶(100272381)及丙酮酸脫氫酶(541919、542376)等一系列基因上、下調(diào)表達來減輕其對干旱脅迫的響應(yīng)，而這些基因的功能多數(shù)與糖酵解生成丙酮酸的途徑相關(guān)聯(lián)。丙酮酸脫羧酶(Pyruvate decarboxylase, PDC)是一種以焦磷酸硫胺素為輔酶的羧基裂解酶，是發(fā)酵途徑中的關(guān)鍵酶。PDC通過與PDH競爭代謝丙酮酸參與植物體內(nèi)能量代謝和物質(zhì)代謝，響應(yīng)植物多種非生物脅迫，如缺氧、高鹽及低溫等[34-35]。研究表明，在逆境脅迫下丙酮酸對植物的呼吸作用有交替作用，而植物保持較高的交替呼吸本身便是一種抗逆機制，對植物抗旱等方面具有重要意義[36-37]。因此，我們推測當(dāng)TS141感受到干旱脅迫后，主要通過誘導(dǎo)其糖酵解途徑來促進根系進行交替呼吸，進而抵御干旱所帶來的損傷。

4 結(jié) 論

本研究以兩份耐旱性不同的玉米自交系昌7-2和TS141為試材，對其幼苗進行12 d干旱脅迫后進行轉(zhuǎn)錄組分析和蛋白質(zhì)組學(xué)分析，以剖析不同自交系對干旱脅迫的響應(yīng)機制及其抵御反應(yīng)。結(jié)果表明，當(dāng)玉米幼苗根系遭遇干旱脅迫后，為抵御干旱脅迫其體內(nèi)會發(fā)生一系列生理生化反應(yīng)，且不同自交系的根系在轉(zhuǎn)錄和翻譯變化方面也存在明顯的差異。此外，在干旱脅迫下，我們在昌7-2和TS141中分別鑒定了19、11個mRNA和蛋白表達趨勢相同的基因、9和0個趨勢相反的基因，并且這些基因主要涉及苯丙素的生物合成、次生代謝產(chǎn)物的生物合成和核糖體等代謝通路，可為今后了解玉米耐旱性的分子機制提供一定的參考。