沈婕 董世滿 李淑霞 彭明

摘? 要：低溫脅迫嚴(yán)重影響木薯的生長、產(chǎn)量以及塊根品質(zhì)。長鏈非編碼RNAs（（lncRNAs）lncRNAs） 作為關(guān)鍵的調(diào)節(jié)因子，在植物響應(yīng)非生物脅迫過程中發(fā)揮重要的功能。然而，lncRNAs在木薯中的功能和調(diào)控機(jī)制尚無報(bào)道。因此，為了研究木薯lncRNAs在低溫脅迫應(yīng)答中的功能，本研究以木薯主栽品種‘60444為材料，利用低溫處理前后的葉片轉(zhuǎn)錄組數(shù)據(jù)，鑒定到一個(gè)受低溫脅迫調(diào)控的lncRNA（cold-responsive lncRNA5，CRR5）。通過對(duì)其序列比對(duì)進(jìn)行同源性分析， RNA二級(jí)結(jié)構(gòu)預(yù)測(cè)和共表達(dá)分析鑒定來分析lncRNA-CRR5響應(yīng)低溫脅迫的分子功能。并根據(jù)基因的共表達(dá)分析推測(cè)CRR5的順式作用和反式作用靶基因。這些結(jié)果為更進(jìn)一步研究CRR5的功能提供了理論依據(jù)。

關(guān)鍵詞：木薯;低溫脅迫;lncRNA;功能分析;同源性分析;RNA二級(jí)結(jié)構(gòu)預(yù)測(cè);共表達(dá)基因分析;靶基因

中圖分類號(hào)：S533? ? ? 文獻(xiàn)標(biāo)識(shí)碼：A

Abstract： Low temperature stress seriously affects the growth and yield of cassava and the quality of roots. Long non-coding RNAs （lncRNAs） function as major regulators， have been shown to play essential roles in plants in response to abiotic stresses. However， the function and regulation mechanism in cassava has not been reported yet. Therefore， in order to investigate the function of lncRNAs in low-temperature response， we identified a cold-responsive lncRNA 5 （CRR5） from the transcriptome data using cassava cultivars ‘60444as the material. We analyzed the molecular function of lncRNA-CRR5 in response to low temperature stress by homology analysis， RNA secondary structure prediction and co-expression analysis. The cis-acting target gene and trans-acting target gene of CRR5 were predicted based on the analysis of gene co-expression. The results would provide a new perspective theoretical basis for further exploring the response of CRR5 to low temperature.

Keywords： cassava; cold stress; lncRNA; functional analysis; homology analysis; RNA secondary structure; co-expression gene analysis; target genes

DOI： 10.3969/j.issn.1000-2561.2021.01.006

木薯（Manihot esculenta Crantz）是熱帶、亞熱帶地區(qū)重要的經(jīng)濟(jì)作物，在中國廣東、廣西、海南等省份被廣泛種植。木薯塊根用途廣泛，可食用、飼用及作為工業(yè)原料生產(chǎn)乙醇和制取淀粉，即便是生產(chǎn)后的木薯渣殘，亦可用作飼料、肥料被開發(fā)與利用。然而，中國熱帶地區(qū)地處熱帶北緣，經(jīng)常遭受寒流侵襲。低溫是一種常見的影響植物地理分布、生長發(fā)育和品質(zhì)產(chǎn)量的非生物脅迫因素，嚴(yán)重時(shí)會(huì)導(dǎo)致植物死亡[1-2]。木薯對(duì)低溫非常敏感，在受到低溫脅迫時(shí)，頂端生長受抑制、植株萎蔫、莖桿壞死，在薯塊形成后受低溫脅迫則導(dǎo)致塊根膨大和淀粉累積受限、產(chǎn)量下降[3-4]。近年來，木薯雖為中國熱區(qū)種植面積最大的薯類作物，但受寒流侵襲和區(qū)域性氣候的限制，中國熱區(qū)生產(chǎn)的木薯自給嚴(yán)重不足，遠(yuǎn)遠(yuǎn)無法滿足國內(nèi)的產(chǎn)業(yè)發(fā)展需要，需大量從國外進(jìn)口。因此，開展耐低溫育種基礎(chǔ)性研究，培育耐低溫、抗旱木薯品種，使木薯種植北移，滿足我國和世界能源工業(yè)對(duì)木薯原材料的需求，是我國木薯產(chǎn)業(yè)發(fā)展的迫切需求。

長鏈非編碼RNA（long non-coding RNA，lncRNA）是一類長度在200 nt以上的非編碼RNA。根據(jù)其在基因組上的位置可將其分為三3種類型：基因間長鏈非編碼RNA（long intergenetic noncoding RNA，linc RNA）、天然反向轉(zhuǎn)錄本（natural antisense transcript，NAT）以及內(nèi)含子長鏈非編碼 RNA（intronic long noncoding RNA）[5-6]。由于lncRNA不具備蛋白質(zhì)編碼功能，因此，很長一段時(shí)間，人們認(rèn)為它是轉(zhuǎn)錄過程中的“噪音”。然而，越來越多的研究結(jié)果表明，lncRNAs作為功能調(diào)控元件在真核生物的基因表達(dá)調(diào)控中起著非常重要的作用[6]。目前關(guān)于長鏈非編碼RNA的工作主要集中在動(dòng)物中，哺乳動(dòng)物lncRNAs參與調(diào)控多種生物學(xué)過程，包括X染色體沉默、基因組印跡、染色體修飾、轉(zhuǎn)錄激活、轉(zhuǎn)錄干擾、核內(nèi)運(yùn)輸?shù)萚7-8]。在植物長鏈非編碼RNA的生物學(xué)功能和作用機(jī)制的研究主要集中在模式作物，如擬南芥、水稻等。相對(duì)于研究較多的非編碼小RNA，目前大部分lncRNAs的功能尚不清楚。近年來，有研究發(fā)現(xiàn)植物lncRNA參與了春化作用、根瘤形成、光形態(tài)建成及逆境響應(yīng)等生物學(xué)過程[9-12]。例如：植物春化開花途徑中的關(guān)鍵調(diào)控基因FLC（Flowering locus C），其mRNA的轉(zhuǎn)錄受到兩個(gè)lncRNA（COOLAIR和COLDAIR）的調(diào)控[13-14]。在豆科植物中，人們鑒定到一個(gè)非常保守的lncRNA，命名為ENOD40（EARLY NODULIN 40）。研究發(fā)現(xiàn)ENOD40參與到共生根瘤的形成過程，在結(jié)瘤過程中，ENOD40作為RNA分子可以影響Mt RBP1（RNA binding protein 1）蛋白的核質(zhì)轉(zhuǎn)運(yùn)[15]。

得益于高通量測(cè)序技術(shù)的發(fā)展，使大量植物基因組信息的解讀成為可能，同時(shí)通過多組學(xué)數(shù)據(jù)的整合，研究者可以清晰地認(rèn)識(shí)基因的功能，研究基因的轉(zhuǎn)錄和調(diào)控機(jī)制。目前，利用組學(xué)數(shù)據(jù)，人們?cè)诒姸嘀参镂锓N中鑒定到大量的lncRNA，這些lncRNA或參與了植物的發(fā)育調(diào)控、或介導(dǎo)了逆境脅迫響應(yīng)過程[16-18]。例如：擬南芥中已經(jīng)發(fā)現(xiàn)1832個(gè)lncRNAs在干旱、低溫、高鹽或脫落酸處理后表達(dá)量發(fā)生顯著改變[19-20]。小麥中鑒定到125個(gè)脅迫（白粉菌感染和熱脅迫）響應(yīng)的lncRNAs[21]。Lv等[22]在玉米中鑒定到7245個(gè)lncRNA，其中637個(gè)與氮營養(yǎng)脅迫相關(guān)。本研究中，前期利用組學(xué)數(shù)據(jù)篩選到一個(gè)低溫響應(yīng)的lncRNA，命名為CRR5（Cold-responsive lncRNA 5）。運(yùn)用序列比對(duì)進(jìn)行同源性分析以及RNA二級(jí)結(jié)構(gòu)預(yù)測(cè)，對(duì)CRR5的潛在功能和靶基因進(jìn)行了預(yù)測(cè)和分析。結(jié)果表明，木薯中CRR5對(duì)低溫脅迫作出了應(yīng)答，通過對(duì)lncRNAs作用機(jī)制的研究，從而篩選特異的功能基因，將為木薯的轉(zhuǎn)基因遺傳育種提供新的思路。

雖然木薯在塊根作物中研究基礎(chǔ)較好，但仍遠(yuǎn)遠(yuǎn)落后于水稻等糧食作物，挖掘木薯特異功能基因，進(jìn)行轉(zhuǎn)基因木薯品質(zhì)改良，是木薯品種改良的有效方法之一。

1? 材料與方法

1.1? 材料

實(shí)驗(yàn)材料為木薯主栽品種，來自中國熱帶農(nóng)業(yè)科學(xué)院熱帶生物技術(shù)研究所研究基地提供。

將組培瓶中的木薯莖去葉片留腋芽和頂芽切成1 cm左右，插入MS培養(yǎng)基中，置于溫室培養(yǎng)，培養(yǎng)溫度為（（26±2）26±2）℃，光照16 h和黑暗8 h。培養(yǎng)1個(gè)月后，分別收集葉片、莖尖、莖段和根等組織部位用于RNA抽提。對(duì)于低溫處理組，把木薯組培苗放于4 ℃分別處理1、3、6、12、24、48、72 h，并且和不做低溫處理的MS固體培養(yǎng)為空白對(duì)照，3次重復(fù)，分別取處理組和對(duì)照組的葉片，液氮速凍后儲(chǔ)存在?80 ℃冰箱用于提取組織總RNA。

1.2? 方法

1.2.1? 葉片總RNA的提取? 取木薯葉片，加入液氮后，迅速充分的地將葉片研磨成粉末狀，再將磨完后的粉末裝入2 mL離心管中，加入500 μL裂解液SL，立即渦旋劇烈震蕩混勻。12000 r/min離心2 min。將上清液轉(zhuǎn)移至過濾柱CS上，12000 r/min離心2 min，小心吸取收集管中的上清液至新的RNase-Free的離心管中，吸頭盡量避免接觸收集管中的細(xì)胞碎片沉淀。緩慢加入0.4倍上清液體積的無水乙醇，混勻，將得到的溶液和沉淀一起轉(zhuǎn)入CR3中，12000 r/min離心15 s，倒掉收集管中的廢液，將吸附柱CR3放回收集管中。向吸附柱CR3中加入350 μL去蛋白液RW1，12000 r/min離心15 s，倒掉收集管中的廢液，將吸附柱CR3放回收集管中。向吸附柱CR3中央加入80 μL的DNase 1工作液，室溫放置15 min（DNase 1工作液的配置：取10 μL DNase 1儲(chǔ)存液放入新的RNA-Free離心管中，加入70 μL RDD緩沖液，輕柔混勻）。向吸附柱CR3中加入350 μL去蛋白液RW1，12000 r/min離心15 s，倒掉收集管中的廢液，將吸附柱CR3放回收集管中。向吸附柱CR3中加入500 μL漂洗液RW，12000 r/min離心15 s，倒掉收集管中的廢液，將吸附柱CR3放回收集管中，并重復(fù)1次。12000 r/min離心2 min，將吸附柱CR3放入一個(gè)新的RNase-Free的離心管中，向吸附柱中央部位懸空滴加50 μL RNase-Free ddH2O，室溫放置2 min，12000 r/min離心1 min，得到RNA溶液。

1.4 1. 2.2? RNA反轉(zhuǎn)錄成cDNA

cDNA反轉(zhuǎn)錄體系（50 ng～-2 μg）總RNA可建立20 μL體系。先進(jìn)行基因組DNA去除體系配制混合液，徹底混勻簡短離心，并置于42 ℃孵育3 min，然后放到冰上。（1）gDNA去除反應(yīng)體系：每個(gè)樣品使用5 gDNA Buffer 2 μL，將Total RNA補(bǔ)足50 ng，隨后使用RNase-Free ddH2O補(bǔ)足到10 μL。（2）反轉(zhuǎn)錄反應(yīng)體系：

每個(gè)樣品使用10 Fast RT Buffer 2 μL，RT Enzyme Mix 1 μL，RNase-Free ddH2O 7 μL。

將反轉(zhuǎn)錄反應(yīng)中的Mix加到gDNA去除步驟的反應(yīng)液中，充分混勻，42 ℃孵育15 min，95 ℃孵育3 min之后放于冰上，得到的cDNA可用于后續(xù)的實(shí)驗(yàn)。

1.51.2.3? 實(shí)時(shí)熒光定量PCR的檢測(cè)

（1）20 μL反應(yīng)體系：

2 SYBRⅡ 10 μL，Primer（正反向） 1.6 μL，ddH2O 7 μL，cDNA 1 μL，SYBR ROX DYe 50 0.4 μL。

（2）反應(yīng)步驟：第一步，95 ℃預(yù)變性30 s，第二步（45個(gè)循環(huán)），95 ℃：5 s，60 ℃：30 s，第三步，95 ℃：15 s，60 ℃：1 min，95 ℃：15 s。

1.61.2.4? CRR5二級(jí)結(jié)構(gòu)預(yù)測(cè)

目前有多種算法可以用來預(yù)測(cè)RNA的二級(jí)結(jié)構(gòu)，其中最小自由能算法（Minimum Free Energy，MFE）作為典型代表使用最為廣泛。本文利用RNA fold軟件，利用最小自由能算法對(duì)CRR5進(jìn)行二級(jí)結(jié)構(gòu)特征進(jìn)行分析。

2? 結(jié)果與分析

2.1? CRR5的篩選與分析

利用木薯莖尖及幼嫩葉片的轉(zhuǎn)錄組數(shù)據(jù)，鑒定到682個(gè)lncRNA，其中318個(gè)受低溫調(diào)控表達(dá)[23]，并重點(diǎn)篩選出lncRNA-CRR5展開后續(xù)功能研究。CRR5全長976 bp，含有3個(gè)外顯子，位一號(hào)染色負(fù)鏈表達(dá)，上下游鄰近蛋白質(zhì)編碼基因分別是Manes.01G155600和Manes.01G155500。通過設(shè)計(jì)特異性PCR引物，以木薯cDNA為模板，經(jīng)過特異性的PCR擴(kuò)增，最終克隆到CRR5全長序列。為了驗(yàn)證CRR5的組織表達(dá)特征，分別采集了在溫室培養(yǎng)30 d，溫室溫度為（26±2）℃，光照16 h和黑暗8 h的木薯葉、莖段、根等不同的組織部位，抽提RNA后檢測(cè)CRR5的表達(dá)水平。結(jié)果發(fā)現(xiàn)CRR5基因具有較為明顯的組織表達(dá)特異性，在根中表達(dá)最低、莖次之、莖尖和葉片中最高，葉柄這種表達(dá)量與根中接近（圖1）。

2.2? CRR5同源序列比對(duì)以及系統(tǒng)進(jìn)化樹的構(gòu)建

為研究木薯CRR5的序列保守程度，從而在進(jìn)化層面上對(duì)其進(jìn)行功能注釋。以CRR5的全長序列為檢索對(duì)象搜索NCBI數(shù)據(jù)庫，選取了12個(gè)與CRR5同源性較高的木薯及其他植物的序列，用MEGA6通過Neighbor joining（NJ）法構(gòu)建了系統(tǒng)發(fā)育樹（圖2）。結(jié)果顯示，CRR5與木薯中的組氨酸?；D(zhuǎn)移酶的同源程度較高，且已有研究結(jié)果表明，組蛋白乙?；揎椩趨⑴c并調(diào)控植物的生長發(fā)育過程中起著重要的調(diào)控作用。此外， CRR5和木薯中一個(gè)功能位置的蛋白質(zhì)編碼基因（Manes.01G194100）同源性最為接近。與此同時(shí)，其在橡膠樹、赤豆中也存在同源基因。

2.3? CRR5二級(jí)結(jié)構(gòu)預(yù)測(cè)

對(duì)CRR5的一級(jí)序列的分析，促使我們獲得其部分信息，而對(duì)二級(jí)結(jié)構(gòu)的研究可以進(jìn)一步探索其具有的生物學(xué)功能。研究表明，單鏈的RNA自身折疊配對(duì)形成螺旋區(qū)和環(huán)，對(duì)RNA二級(jí)結(jié)構(gòu)的預(yù)測(cè)是生物信息學(xué)領(lǐng)域的研究焦點(diǎn)問題。CRR5的二級(jí)結(jié)構(gòu)預(yù)測(cè)結(jié)果如下圖3所示，從圖中可以看出，CRR5二級(jí)結(jié)構(gòu)是多莖環(huán)發(fā)夾結(jié)構(gòu)，其中，3個(gè)主要的分支環(huán)包含12個(gè)莖環(huán)、2個(gè)凸環(huán)結(jié)構(gòu)。藍(lán)色區(qū)域是5端，紅色的則是3端。藍(lán)色區(qū)域表示堿基對(duì)的匹配程度較高，自由能值低，穩(wěn)定性好（圖3）。以往研究表明，RNA的莖環(huán)結(jié)構(gòu)在與RNA結(jié)合蛋白相互作用過程中至關(guān)重要[5]。因此，對(duì)CRR5二級(jí)結(jié)構(gòu)的預(yù)測(cè)，可以為研究其功能和作用機(jī)制等提供重要數(shù)據(jù)。

2.4? 對(duì)不同低溫時(shí)間處理下木薯CRR5表達(dá)的實(shí)時(shí)熒光定量分析

為了研究CRR5響應(yīng)低溫脅迫的動(dòng)態(tài)表達(dá)情況，本實(shí)驗(yàn)用熒光定量PCR的方法對(duì)CRR5進(jìn)行驗(yàn)證。結(jié)果發(fā)現(xiàn)，低溫處理后CRR5的表達(dá)量呈現(xiàn)出整體的上升趨勢(shì)。在低溫脅迫1 h后，處理組中CRR5的相對(duì)表達(dá)量約為對(duì)照組的7倍;3、6、12 h處理后的相對(duì)表達(dá)量約為對(duì)照組的10倍;在低溫處理24 h后，差異表達(dá)最顯著，相對(duì)表達(dá)量約為處理前的20倍。處理48 h和72 h后CRR5的相對(duì)表達(dá)量略有下降（圖4）。

2.5? CRR5共表達(dá)基因鑒定

基于lncRNA與鄰近基因可能存在協(xié)同表達(dá)或調(diào)控關(guān)系，因此，著重分析了CRR5上下游鄰近蛋白質(zhì)編碼基因Manes.01G155600、Manes. 01G155500的表達(dá)。結(jié)果表明，二者在低溫條件下表達(dá)變化不顯著，因此推測(cè)CRR5和臨近的上下游基因不存在共表達(dá)或調(diào)控關(guān)系。而基于基因共表達(dá)的方法，通過分析轉(zhuǎn)錄組中與CRR5的表達(dá)關(guān)聯(lián)系數(shù)大于0.9的蛋白編碼基因的功能?？梢詫?duì)CRR5進(jìn)行功能注釋。利用該方法，在木薯基因組中鑒定了29個(gè)蛋白質(zhì)編碼基因，它們的表達(dá)量與CRR5高度相關(guān)，將這些編碼基因進(jìn)行功能注釋。其中29個(gè)基因參與代謝過程。6個(gè)基因與分子運(yùn)輸相關(guān)，15個(gè)基因參與轉(zhuǎn)錄調(diào)節(jié)。8個(gè)基因參與激素應(yīng)答與信號(hào)轉(zhuǎn)導(dǎo)（（表1）表1）。從而推測(cè)位于編碼基因上下游的這些lncRNA可能參與轉(zhuǎn)錄調(diào)節(jié)、代謝、激素應(yīng)答和信號(hào)轉(zhuǎn)導(dǎo)等生物過程來響應(yīng)低溫脅迫。

在低溫脅迫下，植物細(xì)胞至少有3種溫度感受及信號(hào)傳導(dǎo)途徑：膜流動(dòng)性的改變，葉綠體能量轉(zhuǎn)換失衡和蛋白磷酸化的植物冷響應(yīng)[1， ， 24]。鈣離子作為細(xì)胞內(nèi)的第二信使會(huì)引起一系列代謝及基因表達(dá)的變化，激活相應(yīng)的轉(zhuǎn)錄因子，如MYB、WRKY等來誘導(dǎo)基因的表達(dá)[25-26]。此外，蛋白磷酸化酶和蛋白激酶等也是信號(hào)傳導(dǎo)分子。在轉(zhuǎn)錄水平上的調(diào)控途徑，闡述最清楚的是CBF基因家族的調(diào)控，經(jīng)過轉(zhuǎn)錄組分析，ICE1可結(jié)合在CBF3的啟動(dòng)子MYC元件區(qū)域激活基因轉(zhuǎn)錄，轉(zhuǎn)錄的產(chǎn)物再結(jié)合在COR基因的啟動(dòng)子元件上，激活了冷響應(yīng)有關(guān)基因的轉(zhuǎn)錄從而增強(qiáng)植物的耐寒性[27]。在低溫脅迫下，木薯內(nèi)CRR5與一些重要的蛋白激酶基因存在協(xié)同表達(dá)，比如谷胱甘肽S-轉(zhuǎn)移酶，AUX/IAA蛋白等，還有與RNA轉(zhuǎn)錄相關(guān)的蛋白（（表1）表1）。

3? 討論

低溫和寒害對(duì)熱帶作物木薯的影響巨大。因此，解析低溫響應(yīng)關(guān)鍵基因的功能和作用機(jī)理，將會(huì)為木薯抗逆性研究提供分子基礎(chǔ)，有助于木薯更好地適應(yīng)寒冷環(huán)境。在過去的10年里，諸多證據(jù)表明，lncRNA在植物開花調(diào)控、根瘤形成、逆境脅迫響應(yīng)等方面發(fā)揮重要功能[5]。近年來，隨著新一代測(cè)序技術(shù)的進(jìn)步，許多新的非編碼RNA在不同的植物物種中被發(fā)現(xiàn)，如水稻、玉米等[18， 22， 28]。然而，由于植物對(duì)低溫脅迫表現(xiàn)出復(fù)雜的生理、生化和分子反應(yīng)，以及實(shí)驗(yàn)技術(shù)的局限性，其功能研究一直未得到深入驗(yàn)證，尤其是在木薯中。近年來，科研人員對(duì)lncRNA的研究的大部分還停留在對(duì)不同植物的識(shí)別或調(diào)控現(xiàn)象等方面。根據(jù)lncRNA的一級(jí)結(jié)構(gòu)推測(cè)其二級(jí)結(jié)構(gòu)，是功能蛋白質(zhì)與RNA互作的重要依據(jù)，或成為lncRNA的功能研究的途徑之一。

采用鏈特異性RNA-seq的方法，尋找在低溫脅迫下表達(dá)顯著變化的lncRNA，根據(jù)嚴(yán)格的篩選標(biāo)準(zhǔn)，篩選到了候選lncRNA-CRR5。本研究從木薯葉片中克隆得到CRR5的cDNA全長，基于CRR5的一級(jí)序列，發(fā)現(xiàn)其在植物物種間同源性較低，但系統(tǒng)發(fā)育樹結(jié)果顯示，在木薯中存在12個(gè)與之同源性較高的蛋白質(zhì)編碼基因，這一結(jié)果暗示CRR5可能與這些基因的功能高度相關(guān)。另外CRR5的二級(jí)結(jié)構(gòu)預(yù)測(cè)其存在12個(gè)穩(wěn)定的莖環(huán)結(jié)構(gòu)，通過對(duì)CRR5二級(jí)結(jié)構(gòu)的預(yù)測(cè)，將為尋找與其存在相互作用的蛋白、蛋白質(zhì)復(fù)合體提供參考和理論依據(jù)。

研究表明，lncRNA的表達(dá)具有組織特異性，本研究發(fā)現(xiàn)，CRR5在莖尖、葉片中表達(dá)豐度相對(duì)較高，表明其在莖尖分生組織的維持、葉片分化及形態(tài)建成過程中可能發(fā)揮功能。木薯的莖頂端和葉片是感受溫度變化的重要部位，對(duì)低溫最為敏感。低溫處理后，莖頂端和葉片分別呈現(xiàn)出壞死和萎蔫的癥狀。CRR5在這兩2個(gè)部位的高表達(dá)，表明其在低溫脅迫響應(yīng)中擔(dān)任重要角色。利用基因組共表達(dá)方法，注釋了29個(gè)與CRR5共表達(dá)的蛋白質(zhì)編碼基因，并根據(jù)這些基因的功能，推測(cè)CRR5主要參與了轉(zhuǎn)錄調(diào)節(jié)、代謝、激素應(yīng)答和信號(hào)轉(zhuǎn)導(dǎo)等生物過程。本研究中的發(fā)現(xiàn)為未來開展CRR5的功能與分子機(jī)制的研究及提高木薯抗逆性奠定基礎(chǔ)。

CRR5與鄰近蛋白質(zhì)編碼基因的表達(dá)并不存在協(xié)同或相關(guān)性，表明其并不通過順式作用方式調(diào)控靶基因。而共表達(dá)基因網(wǎng)絡(luò)分析顯示CRR5與真菌素金屬肽酶（M36），Di-葡萄糖結(jié)合內(nèi)質(zhì)網(wǎng)，不分生頂端蛋白，短鏈脫氫酶，木瓜蛋白酶家族半胱氨酸蛋白酶，WRKY轉(zhuǎn)錄因子的共表達(dá)系數(shù)較高，我們推測(cè)CRR5與這些基因的表達(dá)有關(guān)。這些功能基因與激素的調(diào)節(jié)，蛋白質(zhì)的轉(zhuǎn)運(yùn)等功能息息相關(guān)。研究這些基因一方面可以了解低溫脅迫下CRR5如何發(fā)揮功能，也可以通過驗(yàn)證這些功能基因的功能推測(cè)木薯抵御低溫脅迫的分子遺傳通路，為揭示植物響應(yīng)低溫脅迫提供新思路。目前CRR5介導(dǎo)木薯抗低溫的分子機(jī)制尚不清楚，隨后我們還將遺傳轉(zhuǎn)化技術(shù)，構(gòu)建CRR5高/低表達(dá)轉(zhuǎn)基因植株，并通過分析轉(zhuǎn)基因株系在低溫脅迫下的抗寒表型檢驗(yàn)，以及低溫脅迫前后各項(xiàng)生理指標(biāo)的變化確定CRR5的生物學(xué)功能。并通過RNA熒光原位雜交技術(shù)確定其行使功能的具體部位、結(jié)合轉(zhuǎn)錄組學(xué)分析該lncRNA的分子機(jī)制。

參考文獻(xiàn)

Xiong L， Schumaker K S， Zhu J K. Cell signaling during cold， drought， and salt stress[J]. Plant Cell， 2002， 14（Suppl）： S165-183.

Mahajan S， Tuteja N. Cold， salinity and drought stresses： An overview[J]. Archives of Biochemistry & Biophysics， 2005， 444（2）： 139-158.

El-Sharkawy M A. Cassava biology and physiology[J]. Plant Molecular Biology， 2004， 56（4）： 481-501.

An D， Yang J， Zhang P. Transcriptome profiling of low temperature-treated cassava apical shoots showed dynamic responses of tropical plant to cold stress[J]. BMC Genomics， 2012， 13（1）： 1-25.

Zhang Y C， Chen Y Q. Long noncoding RNAs： New regulators in plant development[J]. Biochemical and Biophysical Research Communications， 2013， 436（2）： 111-114.

Chekanova J A. Long non-coding RNAs and their functions in plants[J]. Current Opinion in Plant Biology， 2015（27）： 207-216.

Loda A， Heard E. Xist RNA in action： Past， present， and future[J]. PLoS Genetics， 2019， 15（9）： e1008333.

Zhang X， Hong R， Chen W， et al. The role of long noncoding RNA in major human disease[J]. Bioorganic Chemistry， 2019， 92： 103214.

Heo J B， Sung S. Vernalization-mediated epigenetic silencing by a long intronic noncoding RNA[J]. Science， 2011. 331（6013）： 76-79.

Ding J， Lu Q， Ouyang Y， et al. A long noncoding RNA regulates photoperiod-sensitive male sterility， an essential component of hybrid rice[J]. Proceedings of the National Academy of Sciences， 2012， 109（7）： 2654-2659.

Flynn R A， Chang H Y. Long noncoding RNAs in cell-fate programming and reprogramming[J]. Cell Stem Cell， 2014， 14（6）： 752-761.

Wang Y， Fan X， Lin F， et al. Arabidopsis noncoding RNA mediates control of photomorphogenesis by red light[J]. Proceedings of the National Academy of Sciences of the United States of America， 2014， 111（28）： 10359-10364.

Csorba T， Questa J I， Sun Q W， et al. Antisense COOLAIR mediates the coordinated switching of chromatin states at FLC during vernalization[J]. Proceedings of the National Academy of Sciences of the United States of America， 2014， 111（45）： 16160-16165.

Marquardt S， Raitskin O， Wu Z， et al. Functional consequences of splicing of the antisense transcript COOLAIR on FLC transcription[J]. Molecular Cell， 2014， 54（1）： 156-165.

Campalans A， Kondorosi A， Crespi1 M. Enod40， a short open reading frame–containing mRNA， induces cytoplasmic localization of a nuclear RNA binding protein in Medicago truncatula[J]. Plant Cell， 2017， 29（4）： 1047-1059.

Li L， Eichten S R， Shimizu R， et al. Genome-wide discovery and characterization of maize long non-coding RNAs[J]. Gary Journal Muehlbauer， 2014， 15（2）： R40.

Shuai P， Liang D， Tang S， et al. Genome-wide identification and functional prediction of novel and drought-responsive lincRNAs in Populus trichocarpa[J]. Journal of Experimental Botany， 2014， 65（17）： 4975-4983.

Wang H， Chung P J， Liu J， et al. Genome-wide identification of long noncoding natural antisense transcripts and their responses to light in Arabidopsis[J]. Genome Research， 2014， 24（3）： 444-453.

Liu J， Jung C， Xu J， et al. Genome-wide analysis uncovers regulation of long intergenic noncoding RNAs in Arabidopsis[J]. The Plant Cell， 2012， 24（11）： 4333-4345.

Amor B B， Wirth S， Merchan F， et al. Novel long non-protein coding RNAs involved in Arabidopsis differentiation and stress responses[J]. Genome Research， 2009， 19（1）： 57-69.

Xin M M， Wang Y， Yao Y Y， et al. Identification and characterization of wheat long non-protein coding RNAs responsive to powdery mildew infection and heat stress by using microarray analysis and SBS sequencing[J]. BMC Plant Biology， 2011， 11（1）： 61.

Lv Y， Liang Z K， Ge M， et al. Genome-wide identification and functional prediction of nitrogen-responsive intergenic and intronic long non-coding RNAs in maize （Zea mays L.）[J]. BMC genomics， 2016， 17（1）： 350.

Li S X， Yu X， Lei N， et al. Genome-wide identification and functional prediction of cold and/or drought-responsive lncRNAs in cassava[J]. Scientific Reports， 2017， 7（1）： 45981.

Yamaguchi-Shinozaki K， Shinozaki K. Transcriptional regulatory networks in cellular responses and tolerance to dehydration and cold stresses[J]. Annual Review of Plant Biology， 2006， 57（1）： 781-803.

Saijo Y， kinoshita N， Hata S， et al. A Ca2+-dependent protein kinase that endows rice plants with cold- and salt-stress tolerance functions in vascular bundles[J]. Plant Cell Physiol， 2001， 42（11）： 1228-1233.

Huang G T， Ma S L， Bai L P， et al. Signal transduction during cold， salt， and drought stresses in plants[J]. Molecular Biology Reports， 2012， 39（2）： 969-987.

Agarwal P K， Agarwal P， Reddy M K， et al. Role of DREB transcription factors in abiotic and biotic stress tolerance in plants[J]. Plant Cell Reports， 2006， 25（12）： 1263-1274.

Zhang Y C， Liao J Y， Li Z Y， et al. Genome-wide screening and functional analysis identify a large number of long noncoding RNAs involved in the sexual reproduction of rice[J]. Genome Biology， 2014， 15（12）： 1-16.

責(zé)任編輯：白? 凈