朱 昊 ，張 博 ，張薈薈 ，王玉祥 ，柯 梅 ，李學森 *

（1.新疆農業(yè)大學草業(yè)與環(huán)境科學學院，烏魯木齊 830052；2.新疆畜牧科學院草業(yè)研究所，烏魯木齊 830011）

我國的牧草種質資源極為豐富，目前的研究主要集中在以全球變化與農業(yè)多樣性、農業(yè)系統(tǒng)固碳減排、草地農業(yè)決策系統(tǒng)及林草生態(tài)等宏觀方面[1]；同時以分子生物學為基礎的微觀方面也有許多研究，主要包括DNA水平上的遺傳多樣性和分子標記輔助選擇育種等方面[2-3]。利用新一代高通量測序技術進行牧草種質資源創(chuàng)新及抗逆基因挖掘方面的研究相對較少[4]。

1 RNA-seq技術原理

RNA-seq即轉錄組測序技術?？死锟耍‵rancis Harry Compton Crick,1916-2004）于1958年提出“中心法則”，并于1970年在Nature上的一篇文章中進行了重申。中心法則指出遺傳信息從DNA傳遞給RNA，再從RNA傳遞給蛋白質，即完成遺傳信息的轉錄和翻譯的過程。依據(jù)中心法則，隨著二代測序技術的興起和成熟，轉錄組學研究得到了蓬勃發(fā)展。簡而言之，轉錄組學（transcriptomics）就是從RNA水平研究基因表達的情況，是一門在整體水平上研究細胞中基因轉錄的情況及轉錄調控規(guī)律的學科。

轉錄組（transcriptome）廣義上指某一生理條件下，細胞內所有轉錄產物的集合，包括mRNA、rRNA、tRNA及Non-coding RNA；狹義上指所有mRNA的集合。 RNA-seq，就是把mRNA、small RNA及Noncoding RNA等或者其中一些用高通量測序技術把它們的序列測出來。反映出它們的表達水平。不管物種有無基因序列信息，該技術都能夠全面檢測其所有基因的表達，具有通量高、成本低、定量準、靈敏度高、重復性好以及檢測范圍廣等特點。

2 RNA-seq一般流程

轉錄組測序分析技術的流程如下（圖1）。

圖1 轉錄組測序及后續(xù)生物信息學分析流程圖

轉錄組測序需要注意的事項：（1）取樣，每個處理至少得有3個生物學重復，這個對于無參考基因組序列更為重要；（2）RNA的抽提檢測，理想的RNA樣品是試驗取得成功的一切的前提；（3）測序得到Uigene序列后的生物信息學分析是關鍵，結合試驗設計的目的對數(shù)據(jù)進行深度挖掘。

3 植物應答環(huán)境脅迫方面的RNA-seq研究

植物對逆境脅迫的應答是通過細胞分子水平修飾、自身的生理生化改變來順應不利的生存環(huán)境?；虮磉_具有時空性，抗逆基因的表達尤其明顯。對植物的不同組織器官、不同生長發(fā)育階段、不同環(huán)境脅迫因子響應時的差異表達的基因進行分析篩選，挖掘關鍵基因和抗性之間的關聯(lián)關系，將有助于從轉錄水平上了解脅迫因子的傷害機理及植物適應逆境脅迫的分子機制。

3.1 溫度脅迫

溫度是植物生長發(fā)育過程中十分重要的環(huán)境影響因子。溫度脅迫引起植物生理上的變化主要體現(xiàn)在光合作用、細胞膜結構、滲透調節(jié)物質、抗氧化物質等方面。吳月燕等［5］以杜鵑為研究材料，人工模擬低溫脅迫，進行轉錄組測序以探索杜鵑抗寒的分子機制。從頭測序組裝得到Unigenes序列122 654個，有大量的Unigenes注釋到同源數(shù)據(jù)庫中，部分Unigenes未得到數(shù)據(jù)庫的注釋信息，這可能是杜鵑的特有的基因序列。其轉錄組數(shù)據(jù)分析結果表明，差異表達的基因共有10 829個，這些差異基因為揭示杜鵑抗寒機理，提供了大量可靠的轉錄組水平的基本研究數(shù)據(jù)，還有待后續(xù)研究的進一步挖掘分析。顏朗等［6］通過比較轉錄組學分析，初步揭示了低溫處理下皇竹草的抗寒分子機制。篩選出冬夏兩季皇竹草差異表達基因11 443個，挖掘到在冬季低溫條件下高度表達的耐寒基因112個，且17個耐寒基因在冬季低溫材料中特意表達。韓超等［7］研究了梭梭對高溫脅迫的轉錄組響應情況，分析了差異表達基因的功能和蛋白代謝通路，結果表明，梭梭的同化枝在應答高溫單一環(huán)境因子脅迫時，多個基因、多個生物過程協(xié)同調控。因此，溫度這個單一的環(huán)境因子變化引起梭梭啟動多基因表達系統(tǒng)來應答。

3.2 干旱脅迫

水是生命之源。在干旱脅迫下，植物體通過改變自身組織器官的生理生化狀態(tài)來適應不利的生存環(huán)境[8]。植物對干旱脅迫響應的分子機制通常是通過信號轉導、以不同代謝通路上的基因轉錄和翻譯來調控應答。周炎等[9]分析了烤煙品系LY1360在干旱脅迫下的轉錄組，通過GO和KEGG分析比較干旱脅迫下供試材料的基因和代謝通路差異，結果表明激素信號轉導途徑的關鍵基因在干旱脅迫下表達量均顯著上調。張春榮等[10]對干旱脅迫下的甘草根進行轉錄組測序，發(fā)現(xiàn)適度干旱促進抗逆相關基因的表達和甘草根有效成分積累。GO富集分析顯示天冬酰胺酰內肽酶、β-木糖苷酶、GDP-L-巖藻糖合酶基因下調并顯著富集，這表明干旱脅迫可能抑制干草根初生壁的降解與程序性細胞死亡。孫愛清等［11］對花生干旱處理前后轉錄表達譜變化進行了分析比較，篩選出935個差異表達基因，其中323個基因涉及次生代謝、大分子代謝、能量代謝等代謝途徑。在干旱脅迫條件下參與類黃酮代謝途徑的9個基因的表達量均顯著上調，推斷花生在干旱脅迫下類黃酮對抵御干旱起到十分重要的作用。

3.3 鹽脅迫

鹽脅迫對植物生長的威脅是干旱之外的第二大環(huán)境因子脅迫。有關農作物及鹽生植物對鹽脅迫的轉錄組研究相對比較深入，現(xiàn)已鑒定出部分相關耐鹽基因，涉及到多個代謝途徑。吉福桑等［12］研究了巴西蕉葉片在中性單鹽NaCl脅迫下的轉錄組情況，結果表明脅迫12 h差異基因的上調和下調表達數(shù)均高于脅迫24 h。KEGG分析表明差異基因主要涉及類酮類和丙酮素生物合成途徑以及光合作用等代謝通徑。qRT-PCR驗證了轉錄組測序結果可靠可信，為香蕉耐鹽機理分析提供了科學參考依據(jù)。彭振等［13］以耐鹽和敏鹽的陸地棉為研究材料，分別研究2個材料在NaCl鹽脅迫下幼苗葉片的轉錄因子的表達情況。通過差異基因篩選，鑒定出26個與耐鹽相關的轉錄因子家族，高表達的11個轉錄因子基因與陸地棉的耐鹽密切相關。張麗麗［14］以荒漠鹽堿地常見的鹽生植物鹽穗木為研究材料，分別對在低鹽和高鹽下脅迫3 h的同化枝進行了轉錄組測序。篩選出高鹽和低鹽脅迫下共差異基因1 254個，有關滲透調節(jié)和活性氧清除的基因大多表現(xiàn)為上調基因，說明鹽穗木在短期鹽脅迫下以增強滲透調節(jié)和活性氧清除為應答策略。

4 RNA-seq在牧草抗逆基因挖掘中的研究現(xiàn)狀

張純（2017）[15]梳理了許多植物在逆境脅迫下的轉錄組應答的研究報告，發(fā)現(xiàn)主要集中在水稻、棉花、玉米、大豆、馬鈴薯、丹參、黃瓜、牡丹等這些作物和藥用植物及花卉植物方面；有關牧草植物的研究相對少，僅見星星草、高羊茅及紫羊茅等物種。郭云（2016）[16]基于RNA-seq技術研究了苜蓿根蘗性狀發(fā)生的下調基因情況。齊曉（2017）[17]對低溫脅迫下苜蓿轉錄因子的鑒定和表達進行了分析。Zhipeng Liu（2013）[18]在國外報到了同源四倍體紫花苜蓿“金皇后”從種子萌發(fā)到開花生長發(fā)育的階段的轉錄組表達情況，基于此轉錄組數(shù)據(jù)開發(fā)了紫花苜蓿的EST-SSR分子標記引物。此外，冷暖（2017）[19]研究報到了草地早熟禾在干旱脅迫下的轉錄組應答。通過查閱文獻可知，目前國內有關牧草的轉錄組學研究相對滯后，即使作為牧草之王的紫花苜蓿在轉錄組方面的研究也很少。

5 展 望

牧草種質資源的挖掘大多集中于優(yōu)良抗性種質篩選方面，相對作物來講，有關基因資源分子水平的研究起步晚、研究相對比較淺，有關抗性基因的分析很少。隨著二代測序成本的進一步下降及三代測序的興起，通過RNA-seq技術在牧草優(yōu)良種質資源基因挖掘方面的應用前景十分廣闊。

此外，轉錄組學可以結合代謝組學和蛋白組學一起來研究，多組學的關聯(lián)可為牧草種質資源挖掘和優(yōu)良基因克隆及抗逆機制等方面提供分子水平研究的大數(shù)據(jù)平臺。

參考文獻

[1]張自和．西部大開發(fā)與草業(yè)發(fā)展［J］．草業(yè)科學，2010，27（8）：1-4．

[2]黃文達，趙學勇，趙昕，等．分子標記在種群遺傳學研究中的應用［J］．草業(yè)科學，2010，27（11）：115-120．

[3]張昆，李明娜，曹世豪，等.牧草與草坪草耐鹽基因工程研究進展［J］．中國草地學報，2017，39（2）：96-104．

[4]湯海港，黃艷華，路海博，等.轉錄組測序技術及其在能源草基因挖掘和品種選育中的應用前景分析[J].草地學報，2016，（4）：731-737.

[5]吳月燕，明萌，何靜雯，等．低溫脅迫下杜鵑雜交后代優(yōu)良株系的轉錄組分析［J］．華北農學報，2017，32（5）：52-60．

[6]顏朗，吳燕，陳鼎，等．皇竹草轉錄組分析揭示其抗凍特性［J］．應用與環(huán)境生物學報，2017（3）：0459-0466．

[7]韓超，謝文華，李建貴，等．夏季高溫干旱逆境下梭梭枝條中轉錄組表達譜分析［J］．基因組學與應用生物學，2015 （3 ）：607-614．

[8]HU H，XIONG L．Genetic engineering and breeding of drought-resistant crops［J］．Annual Review of Plant Biology，2014，（65）:715-741．

[9]周炎，丁松爽，趙莉，等 烤煙品系LY1306在PEG干旱脅迫下的生理響應及轉錄組學分析［J］．中國煙草學報，2017，（6）：36-44．

[10]張春榮，桑雪雨，渠萌，等．基于轉錄組測序揭示適度干旱脅迫對甘草根基因表達的調控［J］．中國中藥雜志，2015，40 （24）:4817-4823．

[11]孫愛清，張杰道，萬勇善，等．花生干旱脅迫響應基因的數(shù)字表達譜分析［J］．作物學報，2013，39（6）:1045-1053．

[12]吉福桑，李元元，唐露，等．香蕉葉片相應鹽脅迫轉錄組分析［J］．分子植物育種，2017，15（3）:875-882．

[13]彭振，何守樸，龔文芳，等．陸地棉幼苗NaCl脅迫下轉錄因子的轉錄組學分析［J］．作物學報，2017，43（3）:354-370．

[14]張麗麗，張富春．短期鹽脅迫下鹽穗木的轉錄組分析［J］．植物研究，2018，38（1） :91-99．

[15]張純，唐承晨，王吉永，等.轉錄組學在植物應答逆境脅迫中的研究進展［J］．生物學雜志，2017，34（2）：86-90．

[16]郭云，王鐵梅.基于 RNA-seq 技術的苜蓿根蘗性狀發(fā)生相關下調基因［J］．草業(yè)科學，2016，33（1）：75-85．

[17]齊曉，閔學陽，張正社，等.紫花苜蓿響應低溫脅迫轉錄因子的鑒定及表達分析［J］．草業(yè)科學，2017，34（9）：1824-1829．

[18]Zhipeng Liu,Tianlong Chen,Lichao Ma，et al．Global transcriptome sequencing using the Illumina platform and the development of EST-SSR markers in autotetraploid alfalfa[J]．PLOS ONE,2013，8（12）：1-13．

[19]冷暖，劉曉巍，張娜，等．草地早熟禾干旱脅迫轉錄組差異性分析［J］．草業(yè)學報，2017，26（12）：128-137．