西鄂爾多斯珍稀植物綿刺的避旱相關(guān)基因研究

陶彥彤, 馬 娟, 張穎娟,3

(1.內(nèi)蒙古師范大學(xué) 生命科學(xué)與技術(shù)學(xué)院,內(nèi)蒙古 呼和浩特 010022;2.哈爾濱市香坊中學(xué),黑龍江 哈爾濱 150036; 3.內(nèi)蒙古師范大學(xué) 學(xué)報(bào)編輯部,內(nèi)蒙古 呼和浩特 010022)

西鄂爾多斯是荒漠化草原到草原化荒漠的過(guò)渡地區(qū),年降水量稀少,蒸發(fā)量遠(yuǎn)大于降水量。該地區(qū)是亞洲中部荒漠區(qū)古老特有植物分布最集中、最豐富的區(qū)域,分布著多種國(guó)家級(jí)重點(diǎn)保護(hù)珍稀瀕危植物,有第三紀(jì)孑遺植物“避難所”之稱[1-3]。綿刺(Potaniniamongolica)為荒漠強(qiáng)旱生小灌木,是薔薇科(Rosaceae)綿刺屬(Potaninia)的古地中海孑遺植物,分布于西鄂爾多斯和阿拉善荒漠[4]。綿刺對(duì)干旱氣候有特殊適應(yīng)性,若生長(zhǎng)季降水稀少,則葉、花、果實(shí)全部脫落,并以假死方式進(jìn)入休眠; 若秋季有降水補(bǔ)給,綿刺則迅速吸水,進(jìn)行二次生長(zhǎng)繁殖[4]。目前對(duì)于綿刺的研究多集中于分類、形態(tài)、生理生態(tài)[5-8]等方面,對(duì)其分子生物學(xué)方面的研究較少。研究人員利用AFLP分子標(biāo)記法對(duì)8個(gè)綿刺種群遺傳結(jié)構(gòu)進(jìn)行分析,發(fā)現(xiàn)隨機(jī)遺傳漂變不是影響綿刺種群遺傳多樣性的主要過(guò)程[9-10]; 采用兼并PCR(polymerase chain reaction)和RACE(rapid amplification of cDNA ends)技術(shù)分離了綿刺肌動(dòng)蛋白(actin)基因cDNA序列,驗(yàn)證了PmActin基因可作為分子內(nèi)標(biāo)的可靠性[11]。

隨著高通量測(cè)序技術(shù)的快速發(fā)展,轉(zhuǎn)錄組測(cè)序因能夠快速有效的從整體水平探究基因功能與基因結(jié)構(gòu),在生物學(xué)各個(gè)領(lǐng)域得到廣泛應(yīng)用,已成為研究基因差異表達(dá)的重要手段[12]。2017年Khan等[13]利用轉(zhuǎn)錄組測(cè)序技術(shù)對(duì)自然條件下及干旱脅迫下的耐旱型蠶豆Hosa wi-2進(jìn)行了轉(zhuǎn)錄組分析,獲得了Hosa wi-2相關(guān)抗旱基因。段娜[14]通過(guò)對(duì)干旱脅迫下的唐古特白刺葉片進(jìn)行轉(zhuǎn)錄組學(xué)分析,發(fā)現(xiàn)干旱脅迫下白刺差異表達(dá)基因的主要相關(guān)酶。

本研究利用轉(zhuǎn)錄組Illumina HiSeq X-ten測(cè)序平臺(tái),對(duì)不同時(shí)期綿刺進(jìn)行轉(zhuǎn)錄組測(cè)序,篩選出生長(zhǎng)期與休眠期耐旱相關(guān)的差異表達(dá)的功能基因和調(diào)控基因,并對(duì)基因的相關(guān)信息進(jìn)行了分析,旨在為西鄂爾多斯地區(qū)珍稀瀕危植物的耐旱研究提供依據(jù)。

1 材料與方法

1.1 材料來(lái)源

綿刺樣品取自西鄂爾多斯自然保護(hù)區(qū)內(nèi)(106°45′24″E,39°56′15″N),本區(qū)地處亞非荒漠東部邊緣,生態(tài)環(huán)境脆弱。年均溫8 ℃,極端最高溫度40.2 ℃,地面極端最高氣溫可達(dá)68.5 ℃,≥10 ℃的活動(dòng)積溫3 200 ℃左右; 年降水量為190～250 mm,年蒸發(fā)量約3 000 mm,濕潤(rùn)系數(shù)小于0.13。氣候干旱少雨,風(fēng)大沙多,春遲秋旱,為中溫帶大陸性氣候。

對(duì)自然生境下具有休眠復(fù)蘇特性的珍稀植物綿刺,隨機(jī)選取長(zhǎng)勢(shì)一致的成年植株90株,設(shè)計(jì)3個(gè)階段植物的處理,每10株作為一個(gè)處理,每個(gè)處理重復(fù)3次。90株分為三個(gè)階段并標(biāo)記, (1) 生長(zhǎng)期(Z): 綿刺正常生長(zhǎng)的時(shí)期; (2) 休眠期(X): 綿刺植株下部葉片發(fā)黃達(dá)到50%; (3) 復(fù)蘇期(F): 綿刺第二次生長(zhǎng)。將每個(gè)階段處理的植株選取頂端嫩莖上的葉片,用鋁箔包好,迅速放入布袋中用液氮冷凍保存帶回實(shí)驗(yàn)室,隨后送到北京百邁克公司進(jìn)行RNA提取以及轉(zhuǎn)錄組測(cè)序。

1.2 方法

1.2.1 綿刺總RNA提取、文庫(kù)構(gòu)建及轉(zhuǎn)錄組測(cè)序 提取樣品總RNA并檢測(cè)合格后用帶有Oligo(dT)的磁珠富集真核生物mRNA,加入Fragmentation Buffer將mRNA進(jìn)行隨機(jī)打斷; 以被打斷的mRNA為模板,用六堿基隨機(jī)引物(random hexamers)合成第一條cDNA鏈,然后加入緩沖液、dNTPs、RNase H和DNA polymerase I合成第二條cDNA鏈,利用AMPure XP beads純化cDNA; 純化得到的雙鏈cDNA再進(jìn)行末端修復(fù)、加A尾并連接測(cè)序接頭,使用AMPure XP beads進(jìn)行片段大小選擇,最后通過(guò)PCR富集得到cDNA文庫(kù)。北京百邁克公司使用Illumina HiSeq X-ten高通量測(cè)序平臺(tái)完成測(cè)序。

1.2.2 Unigene功能注釋與差異表達(dá)基因分析 使用BLAST軟件將Unigene序列與NR、GO、Swiss-Prot、COG、KOG、eggNOG4.5、KEGG的數(shù)據(jù)庫(kù)進(jìn)行比對(duì),使用KOBAS2.0技術(shù)鑒定Unigene在KEGG中的KEGG Orthology結(jié)果,使用HMMER軟件與Pfam數(shù)據(jù)庫(kù)比對(duì),獲得Unigene的注釋信息。采用Benjamini-Hochberg法篩選差異表達(dá)基因,并進(jìn)行差異基因的GO和KEGG富集分析。

2 結(jié)果

2.1轉(zhuǎn)錄組測(cè)序與組裝結(jié)果分析

綿刺生長(zhǎng)期(Z)、休眠期(X)和復(fù)蘇期(F)轉(zhuǎn)錄組測(cè)序共得到75.63 Gb Clean Data,每個(gè)樣品Clean Data質(zhì)量值≥30的堿基百分比均不小于92.84%。利用Trinity軟件經(jīng)對(duì)過(guò)de novo組裝后,Clean reads被組裝得到435 276條Transcript(轉(zhuǎn)錄本)和141 608條Unigene,Transcript和Unigene平均長(zhǎng)度分別是1 974.70 bp和702.73 bp; Transcript與Unigene的N50分別為3 269和1 103,由表1可以看出Unigene分布變化隨著長(zhǎng)度區(qū)間的增加而減少,組裝完整性較高。

表1 組裝結(jié)果統(tǒng)計(jì)Tab.1 Assemble result table

2.2 Unigene功能注釋

選擇BLAST參數(shù)E-value≤10-5和HMMER參數(shù)E-value≤10-10,最終獲得87 707個(gè)有注釋信息的Unigene(表2)。比對(duì)到eggNOG和NR數(shù)據(jù)庫(kù)里的Unigene數(shù)較多,分別占All Unigene的90.92%和89.54%,其他依次為Pfam、KOG、Swissprot、GO、KEGG、COG數(shù)據(jù)庫(kù),分別占All Unigene的59.95%、59.21%、55.63%、41.16%、37.97%、29.36%,其中每個(gè)Unigene可能被比對(duì)到多個(gè)數(shù)據(jù)庫(kù)中。

表2 Unigene注釋結(jié)果Tab.2 Results of Unigene annotations

2.3 差異表達(dá)基因的篩選結(jié)果

綿刺三個(gè)時(shí)期兩兩之間的差異表達(dá)基因篩選結(jié)果見(jiàn)表3。生長(zhǎng)期與休眠期差異表達(dá)基因數(shù)共999條,其中上調(diào)基因數(shù)635條,下調(diào)基因數(shù)364條; 生長(zhǎng)期與復(fù)蘇期差異表達(dá)基因數(shù)共2 596條,其中上調(diào)基因數(shù)1 631條,下調(diào)基因數(shù)965條; 休眠期與復(fù)蘇期的差異表達(dá)基因數(shù)共476條,上調(diào)基因數(shù)204條,下調(diào)基因數(shù)272條。生長(zhǎng)期與復(fù)蘇期差異基因表達(dá)的數(shù)最多,說(shuō)明采用的試驗(yàn)樣品時(shí)間的間段增加,參與基因調(diào)控的數(shù)目增多。

表3 差異表達(dá)基因數(shù)目統(tǒng)計(jì)結(jié)果Tab.3 The number of DEGs statistical result

2.4 差異表達(dá)基因功能篩選

基于基因在不同樣品中的表達(dá)量,對(duì)識(shí)別到的差異表達(dá)基因進(jìn)行功能注釋,各差異表達(dá)基因注釋的統(tǒng)計(jì)結(jié)果見(jiàn)表4。注釋到各功能數(shù)據(jù)庫(kù)數(shù)量最多的是生長(zhǎng)期與復(fù)蘇期,為2 475條; 其次是生長(zhǎng)期與休眠期的938條; 最少的是休眠期與復(fù)蘇期的451條。不同差異表達(dá)基因的各功能數(shù)據(jù)庫(kù)注釋到的基因數(shù)目,與總數(shù)變化呈一致性。

表4 注釋的差異表達(dá)基因數(shù)量Tab.4 The number of DEGs annotation

2.5 差異表達(dá)基因GO功能富集

綿刺GO功能富集在生長(zhǎng)期與休眠期有340條樣本基因注釋到“生物學(xué)過(guò)程”,主要聚集在代謝過(guò)程(metabolic process)、細(xì)胞過(guò)程(cellular process)和單生物過(guò)程(single-organism process)等(圖略)。注釋到“分子功能”中有356條差異表達(dá)基因被富集,主要富集條數(shù)最多的為催化活性(catalytic activity)和物質(zhì)結(jié)合(binding)等。注釋到“細(xì)胞組分”中共246條,注釋到的差異表達(dá)基因數(shù)最多的條目為細(xì)胞(cell)、細(xì)胞成分(cell part)和細(xì)胞器(organelle)等,這些富集差異表達(dá)基因數(shù)最多的條目主要與細(xì)胞、細(xì)胞成分有關(guān),表現(xiàn)在干旱脅迫下,細(xì)胞及細(xì)胞成分的性狀發(fā)生變化。在綿刺生長(zhǎng)期與復(fù)蘇期得到注釋的共有1 224條。其中綿刺基因注釋到“生物學(xué)過(guò)程”中共有989條,數(shù)量多于生長(zhǎng)期和休眠期間的注釋,主要富集差異基因較多的仍是代謝過(guò)程、細(xì)胞過(guò)程和單生物過(guò)程等。在“分子功能”中同樣富集條數(shù)最多的為催化活性和物質(zhì)結(jié)合等,這與植物的耐旱功能緊密相關(guān)?！凹?xì)胞組分”中共677條差異基因被富集,細(xì)胞、細(xì)胞成分和細(xì)胞器等得到大量聚集。休眠期與復(fù)蘇期有200條注釋到GO數(shù)據(jù)庫(kù)。其中綿刺基因富集到“生物學(xué)過(guò)程”中共有154條,富集較多的過(guò)程同生長(zhǎng)期與休眠期相同,但數(shù)量較小?！胺肿庸δ堋敝杏?67條差異表達(dá)基因被富集,富集條數(shù)最多的同生長(zhǎng)期與休眠期一致?！凹?xì)胞組分”中共105條差異表達(dá)基因被富集,注釋到基因數(shù)較多的同生長(zhǎng)期與休眠期相符。這些功能表達(dá)說(shuō)明在綿刺耐旱性中起重要作用。

2.6 差異表達(dá)基因KEGG功能富集

綿刺生長(zhǎng)期與休眠期有316條和KEGG富集通路相關(guān)的差異表達(dá)基因,被注釋到通路中的有88個(gè),其中富集顯著性可靠且參考價(jià)值較大包括真核生物核糖體起源(ribosome biogenesis ineukaryotes)、葡萄糖醛酸轉(zhuǎn)換(pentose andglucuronate interconversions)、光合作用(photosynthesis)等。生長(zhǎng)期與復(fù)蘇期有1 162條和此通路相關(guān)的差異表達(dá)基因,被注釋到通路中的有114個(gè),其中富集顯著性可靠且參考價(jià)值較大為光合作用、光合作用-天線蛋(photosynthesis-antenna proteins)等。休眠期與復(fù)蘇期有225條和此通路相關(guān)的差異表達(dá)基因有78個(gè),其中光合作用-天線蛋白、精氨酸與脯氨酸代謝(arginine and proline metabolism)等富集顯著性可靠且參考價(jià)值較大,其他通路的富集顯著可靠性較低。

2.7 綿刺耐旱相關(guān)的功能基因

隨著干旱脅迫的逐漸增加,綿刺在生長(zhǎng)期與休眠期的差異表達(dá)基因中,耐旱相關(guān)的基因有36條Unigenes預(yù)測(cè)為功能基因(表5),主要有滲透調(diào)節(jié)物質(zhì)、內(nèi)源激素類物質(zhì)、活性氧清除類物質(zhì)、果膠酶和保護(hù)生物大分子物質(zhì)等相關(guān)的基因。綿刺在水分脅迫下,植物體內(nèi)相關(guān)的基因均會(huì)發(fā)揮重要的調(diào)節(jié)作用,篩選出6條差異表達(dá)基因與植物激素相關(guān),其中生長(zhǎng)素(auxin)5條占比最多,4條上調(diào)表達(dá)基因,1條為下調(diào)表達(dá)。

表5 干旱誘導(dǎo)與功能相關(guān)的差異表達(dá)基因的表達(dá)情況Tab.5 Expression of differentially expressed genes related to drought and function

續(xù)表5 干旱誘導(dǎo)與功能相關(guān)的差異表達(dá)基因的表達(dá)情況Continude Tab.5 Expression of differentially expressed genes related to drought and function

果膠酶分為果膠酯酶(pectinesterase,PME)和多聚半乳糖醛酸酶(polygalacturonase,PG)兩種,注釋到的差異表達(dá)基因中篩選到10個(gè)基因與果膠酶相關(guān),干旱脅迫下植物體內(nèi)的果膠酯酶有6條均上調(diào)表達(dá)。滲透調(diào)節(jié)類有兩類相關(guān)的差異基因,脯氨酸(proline)和海藻糖(trehalose)類均有1個(gè)上調(diào)表達(dá)的差異基因,有11個(gè)差異表達(dá)基因與活性氧清除類物質(zhì)相關(guān),其中與過(guò)氧化物酶(peroxidase,POD)相關(guān)的基因有5條,均呈上調(diào)表達(dá)。與谷胱甘肽(glutathione)相關(guān)的基因有3條,2條為上調(diào)表達(dá)1條下調(diào)表達(dá)。與抗壞血酸氧化酶(ascorbic acid oxidase,AO)相關(guān)的3條基因均呈上調(diào)表達(dá)。經(jīng)篩選共有6條保護(hù)生物大分子物質(zhì),其中與LEA蛋白相關(guān)的基因有3條,2條上調(diào)表達(dá),1條下調(diào)表達(dá)。與水通道蛋白(aquaporin)相關(guān)的基因有3條,均呈上調(diào)表達(dá)。

2.8 綿刺耐旱相關(guān)調(diào)控基因

干旱脅迫的環(huán)境條件會(huì)在植物體內(nèi)誘導(dǎo)一系列轉(zhuǎn)錄因子基因,干旱脅迫信號(hào)轉(zhuǎn)導(dǎo)過(guò)程中,轉(zhuǎn)錄因子通過(guò)啟動(dòng)多條途徑從多個(gè)層面調(diào)節(jié)、降低脅迫對(duì)植物的傷害[15],對(duì)植物在逆境下的生長(zhǎng)發(fā)育起到重要作用。綿刺生長(zhǎng)期與休眠期得到注釋的差異表達(dá)基因中,耐旱相關(guān)的轉(zhuǎn)錄因子共8條,主要來(lái)自AP2/EREBP、NAC、bHLH、MYB/MYC四個(gè)轉(zhuǎn)錄因子家族的轉(zhuǎn)錄因子。篩選的每個(gè)轉(zhuǎn)錄因子注釋的信息和調(diào)控表達(dá)見(jiàn)表6。屬于NAC的基因共2條,基因c66600.graph_c2和c58089.graph_c0均為下調(diào)表達(dá)。屬于AP2/EREBP的基因共有4條,c57840.graph_c0、c40137.graph_c0和c69583.graph_c0基因上調(diào)表達(dá),AP2/EREBP家族中ERF類轉(zhuǎn)錄因子c65442.graph_c0表現(xiàn)出下調(diào)表達(dá),同是下調(diào)表達(dá)基因的有bHLH的基因c68718.graph_c0。只有一條屬于MYB/MYC的基因c41021.graph_c0為上調(diào)表達(dá)。

表6 干旱誘導(dǎo)與調(diào)控相關(guān)的差異表達(dá)基因的表達(dá)情況Tab.6 Expression of differentially expressed genes related to drought induction and regulation

3 討論

轉(zhuǎn)錄組測(cè)序適用于缺乏基因信息的研究對(duì)象,該技術(shù)可以對(duì)在單核苷酸水平上檢測(cè)指定材料的整體轉(zhuǎn)錄活動(dòng),且不需要參考基因組來(lái)獲得有用的轉(zhuǎn)錄信息[16]。本研究通過(guò)高通量測(cè)序?qū)ψ匀簧诚碌木d刺在轉(zhuǎn)錄組水平分析其耐旱響應(yīng)。在對(duì)綿刺不同時(shí)期的差異表達(dá)基因(DEG)的研究中,三組的差異表達(dá)基因的數(shù)目共有4 071條,而成功被注釋到的差異表達(dá)基因數(shù)量為3 864條這一結(jié)果與其他植物轉(zhuǎn)錄組測(cè)序結(jié)果一致[17]。由于物種基因信息不足,不是所有Unigenne都能得到注釋。但基于被檢測(cè)出來(lái)的DEG,可篩選出多條富集通路和過(guò)程,從而為研究綿刺對(duì)干旱脅迫的響應(yīng)機(jī)制提供理論基礎(chǔ)。在DEG的KEGG富集分析中,核糖體生物合成“ribosome biogenesis in eukaryotes”在生長(zhǎng)期與休眠期的比較中富集顯著性最高,生長(zhǎng)期與復(fù)蘇期中“光合作用”的富集顯著性最高,而休眠期與復(fù)蘇期的“光合作用-天線蛋白”富集顯著最高。說(shuō)明隨著環(huán)境的變化,激發(fā)了各通路的相關(guān)基因的表達(dá),而富集顯著性較高的基因,就可能在綿刺適應(yīng)脅迫過(guò)程中起重要的作用。

由于2017年綿刺復(fù)蘇期返青較晚(10月上旬),平均氣溫接近10 ℃,返青前有降雪,影響綿刺葉片的生長(zhǎng),差異表達(dá)基因篩選只比較了生長(zhǎng)期與休眠期,共有37條Unigenes預(yù)測(cè)為抗旱相關(guān)功能基因?；钚匝跚宄愇镔|(zhì)在植物中廣泛存在,綿刺在逆境脅迫加劇的環(huán)境下,呈差異表達(dá)基因只發(fā)現(xiàn)POD、GSH和AO,并且三種酶的基因表達(dá)量變化大部分呈現(xiàn)上調(diào)表達(dá)。與過(guò)氧化物酶(POD)相關(guān)的基因含有5條,均是上調(diào)表達(dá)。這與段麗[14]對(duì)白刺干旱脅迫下POD相關(guān)差異基因表達(dá)呈上調(diào)的研究結(jié)果相一致; 研究發(fā)現(xiàn)與GSH相關(guān)的基因有3條,其中1條下調(diào)表達(dá)2條上調(diào)表達(dá),而與AO基因相關(guān)的有3條,均是上調(diào)表達(dá),說(shuō)明在休眠期干旱脅迫加劇AO相關(guān)基因表達(dá)量的增加使綿刺耐旱性增強(qiáng)。

滲透調(diào)節(jié)物質(zhì)中與脯氨酸和海藻糖相關(guān)的差異基因各有一條,均呈上調(diào)表達(dá)。研究表明,脯氨酸的過(guò)表達(dá)可增強(qiáng)植物的耐旱能力。杜俊瑞[18]對(duì)白沙蒿相關(guān)抗旱基因研究發(fā)現(xiàn),控制白沙蒿脯氨酸合成的基因表現(xiàn)為上調(diào)表達(dá),與本文研究結(jié)果一致,說(shuō)明脯氨酸相關(guān)基因的上調(diào)表達(dá)有助于綿刺的耐旱性。在綿刺中參與海藻糖相關(guān)的差異表達(dá)基因僅有一條上調(diào)表達(dá)與多數(shù)研究相一致。黃曉鈺[19]對(duì)干旱脅迫下的檸條錦雞兒進(jìn)行轉(zhuǎn)錄組分析結(jié)果表明,海藻糖合成途徑相關(guān)酶基因的表達(dá)量出現(xiàn)不同程度地上調(diào),表明海藻糖合成途徑在檸條錦雞兒應(yīng)對(duì)干旱脅迫時(shí)發(fā)揮重要作用。

姜揚(yáng)[20]對(duì)長(zhǎng)春在逆脅迫下的變化研究發(fā)現(xiàn),干旱脅迫使PME基因表達(dá)上調(diào),有助于保持細(xì)胞水分,增強(qiáng)抗逆能力; 在逆脅迫處理下番茄的PG基因表達(dá)顯著增加[21]。本研究結(jié)果也表明,隨干旱脅迫加劇,綿刺體內(nèi)的果膠酯酶有6條均上調(diào)表達(dá),PG相關(guān)的差異表達(dá)基因有4條均是上調(diào)表達(dá)。綿刺休眠期葉片脫落,PG可能起到了一定的作用,PME和PG可能增強(qiáng)綿刺耐旱能力。

LEA蛋白是一類種子胚胎發(fā)育后期富集的脫水保護(hù)蛋白。Veeranagamallaiah等[22]通過(guò)SDS-PAGE研究發(fā)現(xiàn),隨LEA蛋白的表達(dá)增加來(lái)減弱水分脅迫導(dǎo)致的蛋白質(zhì)凝聚現(xiàn)象,表明LEA蛋白可間接增強(qiáng)植物抗逆性。Zhang等[23]將沙漠牧草隱子草(Cleistgenessongorica)中的脫水蛋白CsLEA轉(zhuǎn)化到紫花苜蓿(MedicagosativaL.),發(fā)現(xiàn)CsLEA在紫花苜蓿中的過(guò)表達(dá)顯著提高了轉(zhuǎn)基因植株的抗旱性和耐鹽性。本研究也發(fā)現(xiàn),綿刺在干旱脅迫加劇的休眠期,2個(gè)LEA蛋白基因呈上調(diào)表達(dá),僅有一個(gè)是下調(diào)表達(dá),上調(diào)的LEA蛋白基因,可能對(duì)綿刺適應(yīng)干旱脅迫起較重要的調(diào)控作用。

AQP是一種跨膜輸水蛋白,能增加生物膜對(duì)水分的透性,調(diào)節(jié)細(xì)胞滲透平衡,脅迫條件下調(diào)節(jié)水分在植物內(nèi)長(zhǎng)距離或短距離運(yùn)輸[24]。其在逆脅迫下的響應(yīng)表達(dá)不同,且每個(gè)AQP也有單獨(dú)的意義。韋興炎[25]對(duì)箭舌豌豆的水通道蛋白相關(guān)基因進(jìn)行研究發(fā)現(xiàn),在根、莖、葉組織中有8個(gè)VsAQP基因在干旱處理過(guò)程中呈現(xiàn)出不同的表達(dá)模式,反映了它們?cè)诩Q豌豆響應(yīng)干旱脅迫中可能具備不同的功能。本研究結(jié)果發(fā)現(xiàn),綿刺在休眠期綿刺有3條與AQP有關(guān)的基因呈上調(diào)表達(dá),具體如何參加綿刺耐旱調(diào)控需要進(jìn)一步研究。

干旱脅迫信號(hào)轉(zhuǎn)導(dǎo)過(guò)程中,轉(zhuǎn)錄因子在水分脅迫中起到不同作用,表達(dá)量的變化也不同。張翠梅[26]在紫花苜蓿響應(yīng)干旱脅迫的研究中發(fā)現(xiàn),NAC轉(zhuǎn)錄因子苜蓿抵抗干旱脅迫過(guò)程中起負(fù)調(diào)控因子。王一航對(duì)舟山新木姜子(Neolitseasericea)幼苗葉片的耐旱性進(jìn)行研究,發(fā)現(xiàn)bHLH家族的轉(zhuǎn)錄因子AtbHLH112可能參與調(diào)控干旱脅迫的ABA信號(hào)轉(zhuǎn)導(dǎo)[27]。MYB44的過(guò)表達(dá)通過(guò)抑制蛋白PHOSPHA Tase 2C(PP2C)基因賦予擬南芥植物抗旱性[28]。TaPIMP1表達(dá)上調(diào),轉(zhuǎn)基因煙草在干旱脅迫期間表現(xiàn)出增強(qiáng)的耐性[29]。與本文篩選綿刺的AP2/EREBP、NAC、bHLH、MYB/MYC四類轉(zhuǎn)錄因子的DEG變化不完全一致,轉(zhuǎn)錄因子在干旱脅迫加劇的休眠期,呈現(xiàn)不同的調(diào)控表達(dá),表明轉(zhuǎn)錄因子對(duì)綿刺的作用相對(duì)復(fù)雜多變,具體轉(zhuǎn)錄因子如何提高綿刺耐旱能力,需要更進(jìn)一步的研究驗(yàn)證。

4 結(jié)論

本研究對(duì)自然生境下不同時(shí)期的綿刺進(jìn)行轉(zhuǎn)錄組測(cè)序,共得到75.63 Gb Clean Data。共有87 707條Unigenes被注釋到NR、GO、COG、KEGG、Swissprot、eggNOG和Pfam數(shù)據(jù)庫(kù)。綿刺不同時(shí)期GO功能富集發(fā)現(xiàn)注釋到GO數(shù)據(jù)庫(kù)中的DEG有440條,在三大類別中富集DEG較多,分別為代謝過(guò)程、催化活性細(xì)胞、細(xì)胞成分等。不同時(shí)期KEGG功能富集分析,主要集中于核糖體生物合成,光合作用,光合作用-天線蛋白等代謝通路等。發(fā)現(xiàn)有37條Unigene與綿刺抗旱相關(guān)的功能基因主要有滲透調(diào)節(jié)物質(zhì)、內(nèi)源激素、果膠酶等相關(guān)基因。在生長(zhǎng)期與休眠期的注釋到的差異表達(dá)基因中,統(tǒng)計(jì)與耐旱相關(guān)的轉(zhuǎn)錄因子共8條,主要來(lái)自AP2/EREBP、NAC、bHLH。