基于RNA-seq的谷子萌芽期抗旱相關(guān)基因挖掘與分析

代小冬 朱燦燦 宋迎輝 王春義 代書桃 秦 娜 李君霞

(河南省農(nóng)業(yè)科學(xué)院糧食作物研究所，河南鄭州 450002)

干旱是限制作物生長發(fā)育的主要非生物逆境災(zāi)害[1]。隨著氣候變化和世界人口的不斷增長，水資源短缺的問題更加突顯，干旱發(fā)生的頻率和嚴(yán)重程度將會進(jìn)一步加劇[2]。因此，解析作物抗旱機(jī)制、培育抗旱作物品種對減輕干旱災(zāi)害及其造成的損失具有重要的意義。

谷子[Setaria italica(L.) Beauv]起源于中國，迄今已有8 700 多年的栽培歷史[3]。谷子是二倍體自花授粉植物，且擁有較小的基因組(約430 Mb)。谷子具有生育期短、單株籽粒多、適應(yīng)性廣、抗旱、耐瘠薄等特點(diǎn)，使其漸漸成為抗旱研究的模式作物。在抗旱性鑒定方法方面，前人相繼確立了谷子萌芽期[4－5]、孕穗期[6]和全生育期[7]的抗旱性鑒定指標(biāo)，并進(jìn)行了抗旱性鑒定。在抗旱相關(guān)QTL 定位方面，Qie 等[8]利用豫谷1 號與青狗尾草構(gòu)建的重組自交系(recombinant inbred lines，RIL)群體定位到18 個(gè)與谷子抗旱相關(guān)的QTL。在參與谷子抗旱的功能基因與蛋白挖掘方面，Tang 等[9]對抗旱品種豫谷1 號和干旱敏感品種安04 在干旱條件下的轉(zhuǎn)錄組進(jìn)行分析，并結(jié)合抗旱相關(guān)QTL 定位結(jié)果，確定了20 個(gè)抗旱候選基因。趙晉鋒等[10]研究認(rèn)為SiCIPK19 基因參與了谷子拔節(jié)期、抽穗期和灌漿期對干旱脅迫的響應(yīng)，推測該基因參與谷子對非生物逆境的應(yīng)答，尤其在抽穗期和灌漿期干旱脅迫應(yīng)答中發(fā)揮重要作用。Wang 等[11]通過對谷子中149 個(gè)bHLH 轉(zhuǎn)錄因子的研究，發(fā)現(xiàn)其可能通過調(diào)控下游抗旱基因的表達(dá)參與谷子的抗旱反應(yīng)。Feng 等[12]通過分析谷子中39 個(gè)核轉(zhuǎn)錄因子的功能，發(fā)現(xiàn)SiNF-YA1 和SiNF-YB8 可以激活相關(guān)抗旱基因，提升谷子生理狀態(tài)，提高谷子抗旱能力。Pan 等[13]利用定量蛋白質(zhì)組學(xué)分析的方法發(fā)現(xiàn)了321個(gè)與谷子抗旱性相關(guān)的蛋白。

目前雖然在谷子抗旱性研究方面取得了一定的進(jìn)展，但是由于谷子抗旱機(jī)制復(fù)雜和研究起步較晚，對谷子抗旱分子調(diào)節(jié)機(jī)制的認(rèn)識還比較有限，因此在谷子抗旱關(guān)鍵基因發(fā)掘和干旱分子響應(yīng)機(jī)制方面的研究有待進(jìn)一步加強(qiáng)。本研究以谷子抗旱品種山西2010 和干旱敏感品種K359*M4－1 為材料，應(yīng)用轉(zhuǎn)錄組測序RNA-seq 技術(shù)對兩個(gè)品種干旱脅迫前后萌發(fā)期種子進(jìn)行轉(zhuǎn)錄組測定，獲得抗旱相關(guān)差異表達(dá)基因(differentially expressed genes，DEG)，并通過GO(gene ontology)富集和KEGG(kyoto encyclopedia of genes and genomes)分析對DEG 進(jìn)行功能分類，以期挖掘谷子抗旱相關(guān)基因，為解析谷子抗旱遺傳機(jī)制和指導(dǎo)谷子抗旱育種提供理論依據(jù)。

1 材料與方法

1.1 試驗(yàn)材料

供試材料K359*M4－1 和山西2010 分別來源于河北省農(nóng)林科學(xué)院谷子研究所和山西省農(nóng)業(yè)科學(xué)院經(jīng)濟(jì)作物研究所。

1.2 谷子萌芽期干旱脅迫處理

選用0、10%、15%、20%、25%、30% 6 個(gè)聚乙二醇6000(polyethylene glycol 6000，PEG－6000)濃度進(jìn)行預(yù)試驗(yàn)，篩選出谷子萌芽期抗旱鑒定的PEG 6000 適宜濃度為20%。利用濃度為20%的PEG 6000 對山西2010 和K359*M4－1 進(jìn)行萌芽期抗旱性鑒定。每份材料選擇飽滿的種子，首先用5%次氯酸鈉消毒15 min，然后用ddH2O 沖洗干凈。選擇50 粒滅過菌的種子放入直徑為9 cm 的培養(yǎng)皿中，再將培養(yǎng)皿放入培養(yǎng)箱，28℃暗培養(yǎng)。干旱脅迫組培養(yǎng)皿中加入5 mL 20%PEG 6000，對照組培養(yǎng)皿中加入5 mL ddH2O。每個(gè)試驗(yàn)設(shè)3 次重復(fù)。

分別調(diào)查第2、第4、第6、第8 天種子的發(fā)芽數(shù)，并將第2、第4、第6、第8 天種子的萌發(fā)率記為nd2、nd4、nd6、nd8。參考Bouslama 等[14]的方法計(jì)算萌發(fā)指數(shù)(promptness index，PI)和萌發(fā)抗旱指數(shù)(sprout index of drought resisting，SIDR)，PI ＝1.00nd2 ＋0.75nd4 ＋0.50nd6 ＋0.25nd8，SIDR＝處理萌發(fā)指數(shù)/對照萌發(fā)指數(shù)。

1.3 RNA 文庫構(gòu)建和測序

文庫構(gòu)建和轉(zhuǎn)錄組測序由北京諾禾致源生物信息科技有限公司完成。

1.4 數(shù)據(jù)的處理與篩選

序列數(shù)據(jù)由Illumina HiSeq 2500 測序得到的原始圖像數(shù)據(jù)經(jīng)base calling 轉(zhuǎn)化后形成FASTQ 格式的文件，該序列數(shù)據(jù)包含reads 的序列及堿基的測序質(zhì)量，在經(jīng)過去除和過濾含adapter、N 過多或低質(zhì)量堿基的dirty raw reads 后獲得高質(zhì)量的clean reads，用于后續(xù)的信息分析。

利用轉(zhuǎn)錄組數(shù)據(jù)比對軟件TopHat2 將所有的clean reads 與谷子(豫谷1 號)參考基因組序列進(jìn)行比對。利用比對的reads 計(jì)算基因的表達(dá)量以確保結(jié)果的高準(zhǔn)確性。用FPKM (fragments per kilobase of transcript per million mapped reads)值衡量計(jì)算表達(dá)量。分別利用GOseq R 和KOBAS 軟件包進(jìn)行差異表達(dá)基因的GO 和KEGG 富集分析(偽發(fā)現(xiàn)率false discovery rate，F(xiàn)DR<0.05)。

本研究報(bào)道的原始數(shù)據(jù)已存入中國科學(xué)院北京基因組研究所大數(shù)據(jù)中心基因組序列檔案庫，accession ID:CRA002537。公眾可訪問網(wǎng)址:https:/ /bigd.big.ac.cn/gsa 獲取數(shù)據(jù)。

1.5 實(shí)時(shí)熒光定量PCR 分析

為了驗(yàn)證谷子萌芽期干旱脅迫轉(zhuǎn)錄組結(jié)果的準(zhǔn)確性，從DEG 數(shù)據(jù)庫中隨機(jī)挑選28 個(gè)基因，利用Primer 6.0 設(shè)計(jì)引物(表1)，分析實(shí)時(shí)熒光定量PCR(realtime quantitative PCR，RT-qPCR)結(jié)果是否與轉(zhuǎn)錄組測序結(jié)果一致。

表1 RT-qPCR 引物序列Table 1 RT-qPCR primer sequence

2 結(jié)果與分析

2.1 干旱處理對谷子發(fā)芽的影響

以濃度為20%的PEG 6000 對91 份谷子材料的萌芽期抗旱性進(jìn)行評價(jià)。由圖1可知，山西2010 的萌發(fā)抗旱指數(shù)(sprout index of drought resisting，SIDR)達(dá)到92%，抗旱性較強(qiáng)，而K359*M4－1 的萌發(fā)抗旱指數(shù)僅為64%，抗旱性較弱。

根管充填即刻質(zhì)量評價(jià)標(biāo)準(zhǔn)[1]（1）適充：根充材料距解剖根尖≤2mm,根管三維充填封閉嚴(yán)密；（2）超充：根管充填材料超出根尖孔；（3）欠充：根充材料距解剖根尖＞2mm或者根管三維充填不嚴(yán)密，X線顯示根充物不致密。

圖1 干旱脅迫對山西2010 和K359*M4－1 發(fā)芽的影響Fig.1 Effects of drought stress on germination of Shanxi 2010 and K359*M4－1

2.2 干旱處理萌芽期谷子種子的轉(zhuǎn)錄組分析

利用RNA-seq 技術(shù)對山西2010 和K359*M4－1萌芽期種子進(jìn)行轉(zhuǎn)錄組測序。將干旱處理前后山西2010 樣品分別命名為SC 和ST，將干旱處理前后K359*M4－1 樣品分別命名為KC 和KT。將校正后P<0.05 作為標(biāo)準(zhǔn)篩選每個(gè)比較組的DEG。結(jié)果顯示，山西2010 干旱處理前后(SC vs ST)共鑒定到2 300 個(gè)DEG，其中上調(diào)基因有1 002 個(gè)，下調(diào)基因有1 298 個(gè)。K359*M4－1 干旱處理前后(KC vs KT)共鑒定到3 652個(gè)DEG，其中上調(diào)基因1 861 個(gè)，下調(diào)基因1 791個(gè)。綜上可知，在干旱脅迫下，與抗旱品種山西2010相比，干旱敏感品種K359*M4－1 的代謝發(fā)生了更大的變化。

韋恩圖分析顯示，干旱脅迫處理前后山西2010 和K359*M4－1 之間存在1 268 個(gè)共有DEG，在山西2010 和K359*M4－1 2 個(gè)品種中分別存在1 032 個(gè)和2 384 個(gè)特有的DEG(圖2)。這些基因參與了谷子對干旱脅迫的響應(yīng)，并可以用來解釋2 個(gè)谷子品種的抗旱性差異。

圖2 DEG 表達(dá)韋恩圖Fig.2 Venn diagrams of DEGs among different comparisons

2.3 DEG 功能分析

為了進(jìn)一步探索谷子對干旱脅迫的響應(yīng)機(jī)制，將矯正后P<0.05、閾值(｜Log 2 fold change ｜≥2)作為標(biāo)準(zhǔn)篩選每個(gè)比較組的DEG。分別在山西2010 和K359*M4－1 中鑒定到82 個(gè)(40 個(gè)上調(diào)基因，42 個(gè)下基因調(diào))和112 個(gè)(51 個(gè)上調(diào)基因，61 個(gè)下調(diào)基因)DEG(表2)。DEG 主要為編碼類鋅誘導(dǎo)的促進(jìn)因子蛋白、類萌發(fā)素蛋白、蛋白質(zhì)磷酸酶、轉(zhuǎn)運(yùn)蛋白、胚胎發(fā)育晚期豐富蛋白、轉(zhuǎn)錄因子、過氧化物酶等的基因。

表2 部分谷子萌芽期響應(yīng)干旱脅迫DEGTable 2 Part of DEGs response to drought stress at germination stage of foxtail millet under drought stress

2.4 DEG 功能富集分析

通過對不同比較組鑒定到的DEG 進(jìn)行GO 和KEGG 代謝途徑富集分析。根據(jù)GO 功能注釋，SC vs ST 和KC vs KT 中的DEG 分別富集在52 個(gè)和21 個(gè)生物學(xué)過程。主要集中在催化活性(GO:0003824)、單個(gè)有機(jī)體代謝過程(GO:0044710)、氧化還原過程(GO:0055114)、血紅素結(jié)合(GO:0020037)、四吡咯結(jié)合(GO:0046906)等生物學(xué)過程(表3)。

表3 干旱脅迫下谷子DEG 的GO 富集分析Table 3 GO enrichment analysis of s DEGs in foxtail millet under drought stress

KEGG 富集分析表明，不同比較組鑒定到的DEG主要富集在淀粉和蔗糖代謝，植物激素信號轉(zhuǎn)導(dǎo)，光合作用－天線蛋白，苯丙素生物合成，角質(zhì)、亞氨酸和蠟生物合成，次生代謝物的生物合成等途徑(圖3)。

圖3 干旱脅迫下山西2010(A)和K359*M4－1(B)中DEG 的KEGG 富集圖Fig.3 KEGG enrichments map of the DEGs in Shanxi 2010(A) and K359*M4－1 (B) under drought stress

2.5 實(shí)時(shí)熒光定量PCR 驗(yàn)證

圖4 實(shí)時(shí)熒光定量PCR 結(jié)果Fig.4 Confirmation of the transcriptomic profiles of selected genes by RT-qPCR

3 討論

RNA-seq 可以全面地獲得某一物種特定細(xì)胞或組織在某一狀態(tài)下的幾乎所有的轉(zhuǎn)錄本及基因序列，是一種快速、有效鑒定差異表達(dá)基因的方法，是揭示植物抗逆機(jī)理以及篩選抗逆基因的主要研究技術(shù)[15－20]。本研究利用RNA-seq 在全基因組表達(dá)水平上鑒定了谷子萌芽期的抗旱基因。通過對干旱處理前后抗旱品種和干旱敏感品種的轉(zhuǎn)錄組比較，在山西2010 和K359*M4－1 中分別發(fā)現(xiàn)了2 300 個(gè)和3 652 個(gè)可能參與干旱脅迫響應(yīng)的DEG。這些基因的表達(dá)和調(diào)控能使谷子抵御干旱脅迫的危害。進(jìn)一步比較發(fā)現(xiàn)，K359*M4－1 中的DEG 明顯多于山西2010 中的DEG，表明在干旱脅迫下，與抗旱品種相比，干旱敏感品種的代謝發(fā)生了更大的變化。上述結(jié)果與前人研究結(jié)果一致[21－22]。

本研究中，DEG 功能富集分析結(jié)果表明，DEG 顯著富集于氧化還原相關(guān)的生物學(xué)過程，表明植物可以通過調(diào)控活性氧代謝來保護(hù)細(xì)胞膜的穩(wěn)定性，從而提高植物的抗旱性[23－24]。苯丙素生物合成、次生代謝物生物合成、淀粉和蔗糖代謝、苯丙氨酸代謝等途徑顯著富集在兩個(gè)比較組中。植物次生代謝產(chǎn)物是在長期進(jìn)化過程中植物與生物和非生物因素相互作用的結(jié)果，在植物對環(huán)境脅迫的適應(yīng)過程中發(fā)揮重要作用[25]。當(dāng)植物處于干旱環(huán)境時(shí)，植物次生代謝產(chǎn)物會通過減緩淀粉代謝來保存能量以延長細(xì)胞的存活時(shí)間[15]。說明富集在這些途徑上的基因在谷子響應(yīng)干旱脅迫方面發(fā)揮了重要作用。

本研究獲得了很多與干旱脅迫相關(guān)的DEG，如編碼類鋅誘導(dǎo)促進(jìn)因子蛋白的基因Seita.3G089200 在兩個(gè)比較組中的表達(dá)都表現(xiàn)明顯上調(diào)。前人研究表明，類鋅誘導(dǎo)促進(jìn)因子可以通過調(diào)節(jié)氣孔的關(guān)閉提高植物的抗旱性[26]，因此，推測Seita.3G089200 參與了谷子對干旱的響應(yīng)。

在植物中，鉀不僅是生長發(fā)育所必需的宏量營養(yǎng)素，而且是介導(dǎo)植物響應(yīng)脅迫的主要信號[27]。鉀吸收透性酶/高親和力鉀轉(zhuǎn)運(yùn)體/鉀轉(zhuǎn)運(yùn)體(K＋uptake permease/high-affinity K＋transporter/K＋transporter，KUP/HAK/KT)轉(zhuǎn)運(yùn)蛋白在調(diào)控甘蔗響應(yīng)低鉀、干旱和鹽脅迫方面發(fā)揮了重要作用[28]。本研究發(fā)現(xiàn)2 個(gè)編碼KUP/HAK/KT 轉(zhuǎn)運(yùn)蛋白的基因(Seita.9G182400和Seita.9G182500)在干旱脅迫下上調(diào)表達(dá)。表明這些基因可能是參與谷子抗旱的關(guān)鍵基因。ABC 轉(zhuǎn)運(yùn)蛋白參與跨膜運(yùn)輸、信號傳導(dǎo)和細(xì)胞解毒等多種重要的生理過程，使植物能夠適應(yīng)環(huán)境的變化，并調(diào)節(jié)植物對生物和非生物脅迫的反應(yīng)[29]。在擬南芥(Arabidopsis thaliana)中超表達(dá)AtABCG36 可以顯著增強(qiáng)相株的抗旱性[30]。

本研究發(fā)現(xiàn)4 個(gè)編碼ABC 轉(zhuǎn)運(yùn)蛋白的基因(Seita.9G001100、Seita.9G447200、Seita.3G042600 和Seita.2G290000)參與了谷子對干旱脅迫的響應(yīng)。類萌發(fā)素蛋白基因受多種環(huán)境信號誘導(dǎo)，在植物生長發(fā)育及對非生物脅迫、生物脅迫的響應(yīng)中起著重要的作用[31－34]。超表達(dá)大豆(Glycine max)的GmGLP7 可以提高擬南芥的抗旱性和抗氧化能力[35]。本研究有10個(gè)編碼類萌發(fā)素蛋白的基因在干旱脅迫下差異表達(dá)，推測這些基因與谷子抗旱性相關(guān)。蛋白質(zhì)磷酸酶是許多信號轉(zhuǎn)導(dǎo)途徑的調(diào)控因子，在植物抗氧化、耐冷、抗旱等過程中起關(guān)鍵作用[36]。如煙草(Nicotiana tabacum)中的NtPP2C1 的表達(dá)受干旱誘導(dǎo)，在干旱脅迫條件下的表達(dá)顯著提高[37]。本研究中，在干旱條件下，谷子中3 個(gè)編碼蛋白質(zhì)磷酸酶的基因(Seita.7G021400、Seita.6G005300 和Seita.3G218800)表達(dá)上調(diào)，可能參與了谷子對干旱脅迫的應(yīng)答。轉(zhuǎn)錄因子可通過調(diào)控多個(gè)下游靶基因的表達(dá)，誘導(dǎo)多種非生物脅迫的保護(hù)機(jī)制[38－40]。如NAC、MYB、bHLH、AP2/ERF 等轉(zhuǎn)錄因子直接或間接參與抗旱響應(yīng)[41－43]。本研究發(fā)現(xiàn)編碼NAC 轉(zhuǎn)錄因子基因Seita.2G028200、編碼MYB 相關(guān)蛋白基因Seita.1G374700、bH035 的同源基因Seita.7G045200 的表達(dá)在干旱脅迫前后差異顯著，推測這些基因可能通過調(diào)節(jié)相關(guān)基因表達(dá)參與谷子對干旱的響應(yīng)。LEA 蛋白參與植物對各種非生物和生物脅迫響應(yīng)，包括鹽害[44]、干旱[45]、熱害[46]、冷害[47]等。提高LEA 蛋白基因的表達(dá)可以提高植物的耐旱性[48－50]。本研究中，編碼LEA 蛋白的Seita.2G437900、Seita.9G420500、Seita.5G287800 和Seita.3G162400 在干旱脅迫下表達(dá)顯著上調(diào)，表明這些基因可能與谷子抗旱相關(guān)。當(dāng)植物遭受干旱脅迫時(shí)，其通過抗氧化酶系統(tǒng)保護(hù)組織免受由于活性氧積累導(dǎo)致的氧化損傷[51]。過氧化物酶是一種抗氧化酶，在植物抗旱脅迫過程中發(fā)揮重要的作用[52]。本研究發(fā)現(xiàn)，Seita.7G247400 和Seita.3G052900 在谷子受到干旱脅迫時(shí)，表達(dá)顯著上升，推測它們可能參與了干旱脅迫下谷子中活性氧的清除過程。

4 結(jié)論

本研究通過對干旱脅迫前后2 個(gè)谷子品種的萌芽期種子進(jìn)行轉(zhuǎn)錄組分析，獲得了一批包括編碼類萌發(fā)素蛋白、蛋白磷酸酶、LEA 蛋白、轉(zhuǎn)錄因子、ABC 轉(zhuǎn)運(yùn)蛋白等與干旱脅迫相關(guān)的基因。DEG 顯著富集在淀粉和蔗糖代謝，植物激素信號轉(zhuǎn)導(dǎo)，光合作用－天線蛋白，苯丙素生物合成，角質(zhì)、亞氨酸和蠟生物合成，次生代謝物生物合成等途徑，表明這些途徑參與了谷子的抗旱反應(yīng)。本研究結(jié)果為解析谷子抗旱機(jī)制、篩選谷子抗旱品種奠定了很好的理論基礎(chǔ)。