段珊 ，竇玲玲 ，程帥帥 ，郭亞寧 ，龐朝友 ，馬啟峰 ，魏恒玲 ，范術(shù)麗 ，王寒濤，喻樹(shù)迅*

（1.西北農(nóng)林科技大學(xué)農(nóng)學(xué)院，陜西楊陵712100；2.中國(guó)農(nóng)業(yè)科學(xué)院棉花研究所/棉花生物學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，河南安陽(yáng)455000）

植物microRNA （miRNA）是1類長(zhǎng)度約為21 nt的內(nèi)源性無(wú)編碼功能的單鏈小分子RNA[1]，通過(guò)降解靶基因或抑制其翻譯在轉(zhuǎn)錄后的水平上調(diào)節(jié)基因表達(dá)[2]，從而廣泛參與植物生長(zhǎng)發(fā)育的一系列重要進(jìn)程，包括調(diào)控根的形成和發(fā)育、葉形態(tài)發(fā)生和極性及植物幼年期到生殖期轉(zhuǎn)變。同時(shí)，miRNA在植物對(duì)外界環(huán)境應(yīng)答和非生物脅迫方面發(fā)揮重要的調(diào)控作用[3]。

棉花是我國(guó)最重要的經(jīng)濟(jì)作物之一,葉片是棉花最重要的源器官，葉片早衰嚴(yán)重影響棉產(chǎn)量和纖維品質(zhì)，進(jìn)而造成棉花減產(chǎn)[4]。葉片衰老是葉片發(fā)育的最后階段，多種激素參與了這一過(guò)程的調(diào)控，乙烯[5]、脫落酸[6]、茉莉酸[7]和水楊酸[8]可加速葉片衰老，細(xì)胞分裂素[9]、生長(zhǎng)素[10]和赤霉素[11]則延緩葉片衰老。

miR160在植物的發(fā)育過(guò)程中起到了重要的調(diào)控作用。在擬南芥中，miR160通過(guò)靶向調(diào)控ARF10、ARF16和ARF17，參與調(diào)控植物生長(zhǎng)發(fā)育及生物脅迫等[12]。在大豆中，過(guò)表達(dá)gm-MIR160A可增加葉綠素含量，同時(shí)其靶基因GmARF和衰老標(biāo)記基因下調(diào)，表明gma-miR160通過(guò)抑制靶基因表達(dá)對(duì)葉片的衰老進(jìn)程進(jìn)行負(fù)調(diào)控[13]。將普通野生稻MIR160f前體轉(zhuǎn)入擬南芥，可使蓮座葉數(shù)目減少，首次抽薹時(shí)間提前，表明MIR160f可能參與了植物花期的調(diào)控[14]。而在龍眼中，miR160及其靶基因不僅參與激素信號(hào)傳導(dǎo)，還參與體細(xì)胞胚胎發(fā)生過(guò)程[15]。番茄的miR160a通過(guò)抑制其靶基因SLARF10可影響葉片的長(zhǎng)出和早期果實(shí)的發(fā)育[16]。

棉屬A組、D組和AD組基因組測(cè)序已經(jīng)完成[17-19]，可通過(guò) psRNATarget（http://plantgrn.noble.org/psRNATarget/）預(yù)測(cè)miR160的靶基因。經(jīng)分析其為核酸外切體組分RRP41-like（Ribosomal RNA-processing protein 41），它的功能需要進(jìn)一步驗(yàn)證。

本研究利用擬南芥內(nèi)源miR319a前體骨架[20]，構(gòu)建了棉花miR160的過(guò)表達(dá)載體，通過(guò)農(nóng)桿菌介導(dǎo)法轉(zhuǎn)入擬南芥等方法，初步明確了ghr-miR160的作用，以期為棉花葉片衰老相關(guān)基因的研究提供新思路。

1　材料與方法

1.1　材料

棉花miR160及其靶基因GhRRP41L在葉片中的相對(duì)表達(dá)量分析所用材料為早熟早衰品種中棉所10號(hào) （CCRI 10）和早熟不早衰品種遼4086（Liao 4086），種植于中國(guó)農(nóng)業(yè)科學(xué)院棉花研究所（河南省安陽(yáng)市）溫室。擬南芥為哥倫比亞型，于本實(shí)驗(yàn)室溫室種植，生長(zhǎng)條件為22℃，16 h光照、8 h黑暗。pMD19-T載體購(gòu)于TAKARA公司，pBI121載體為本實(shí)驗(yàn)室保存，大腸桿菌DH5α購(gòu)于北京全式金有限公司，農(nóng)桿菌菌株LBA4404A為本實(shí)驗(yàn)室保存。

1.2　取樣方法

棉花子葉展平后，開(kāi)始進(jìn)行第1次取樣，每隔1周取1次葉片，共取4次樣，3個(gè)重復(fù)，于液氮中速凍后存于－80℃冰箱備用。

1.3　DNA及RNA提取

采用改良的CTAB法提取擬南芥基因組DNA[21]，并將其保存于－20℃冰箱備用。使用Small RNA Purification Kit（Sigma）進(jìn)行 microRNA的提取，按照TransScript miRNA First-Strand cDNA Synthesis SuperMix（北京全式金公司）說(shuō)明書進(jìn)行反轉(zhuǎn)錄，microRNA及cDNA存于－20℃冰箱備用。RNA的提取參照多糖多酚類RNAprep pure plant kit（天根生化科技有限公司）說(shuō)明書，應(yīng)用微量分光光度計(jì)和1%（質(zhì)量分?jǐn)?shù)）瓊脂糖凝膠進(jìn)行濃度和質(zhì)量檢測(cè)，合格的樣品使用 PrimeScriptRT reagent kit（TaKaRa）合成 cDNA模板后，于－20℃保存或直接用于反轉(zhuǎn)錄實(shí)時(shí)熒光定量聚合酶鏈?zhǔn)椒磻?yīng) （Quantitative reverse transcription polymerase chainreaction，qRT-PCR）。

1.4　MicroRNA160過(guò)表達(dá)載體的構(gòu)建

從 miRBase （http://www.mirbase.org/，release 21）獲取擬南芥pre-miR319a序列，將ath-miR319a和ath-miR319a*替換為ghr-miR160和ghr-miR 160*，再?gòu)臄M南芥pre-miR319a上下游各取80 bp堿基，新的序列送于蘇州金唯智生物科技有限公司合成。以合成序列為模板，用AtmiR319a引物及 Mighty Amp○RDNA Polymerase Ver.2（TaKaRa）酶獲得所需序列后，回收目的片段并連接pMD19-T載體，轉(zhuǎn)入大腸桿菌感受態(tài)細(xì)胞DH5α，用菌液PCR方法鑒定獲得陽(yáng)性克隆，送至金維智公司進(jìn)行測(cè)序。

對(duì)序列正確的菌液提取質(zhì)粒檢測(cè)后與經(jīng)過(guò)雙酶切 （酶切位點(diǎn)分別是BamHⅠ和SacⅠ）的pBI121載體進(jìn)行重組，轉(zhuǎn)化大腸桿菌DH5α，用測(cè)序驗(yàn)證正確的大腸桿菌重組質(zhì)粒轉(zhuǎn)化農(nóng)桿菌。

1.5　MicroRNA160過(guò)表達(dá)載體轉(zhuǎn)化擬南芥

利用Bechtold等建立的蘸花法轉(zhuǎn)化擬南芥[22]。待種子成熟后收獲，為T0種子，消毒滅菌后種于具卡那霉素抗性的1/2 MS培養(yǎng)基，于4℃冰箱春化2 d后，移至16 h光照、8 h黑暗，21℃的培養(yǎng)箱中，生長(zhǎng)14 d后，將其移栽至營(yíng)養(yǎng)土中。待擬南芥生長(zhǎng)約40 d，取蓮座葉進(jìn)行PCR檢測(cè)，分單株收獲陽(yáng)性株，為T1種子，以此類推，直至T3純合株系。選取陽(yáng)性純合株系，進(jìn)行g(shù)hr-miR160的轉(zhuǎn)錄水平分析及表型觀察。

表1　引物序列Table 1　The primer sequences in this study

1.6　qRT-PCR分析

利用軟件Oligo7設(shè)計(jì)熒光定量引物（表1），以棉花葉片和擬南芥葉片的cDNA為模板，使用康為世紀(jì)的試劑盒UltraSYBR mixture with ROX在ABI7500 Real-time PCR system熒光定量?jī)x上進(jìn)行qRT-PCR分析，反應(yīng)條件為3步法，擴(kuò)增程序?yàn)椋侯A(yù)變性 95℃ 10 min；95℃ 10 s；60℃ 30 s；72℃32 s；40個(gè)循環(huán)。每個(gè)樣品技術(shù)重復(fù)3次，采用 2-△△Ct法，內(nèi)參基因?yàn)镚hU6 和GhActin。

2　結(jié)果與分析

2.1　MicroRNA160家族分析

在miRBase中共獲得pre-miR160數(shù)據(jù)141條，miR160數(shù)據(jù)150條，分布在19個(gè)物種共41種植物中。利用DNAMAN對(duì)150條miR160數(shù)據(jù)進(jìn)行多重序列比對(duì) （圖1-A），以棉花miR160為例，第11位、14位和16位的C、G和A堿基在所有miR160中高度保守。利用DNAMAN和i-TOL（http://itol.embl.de/index.shtml） 對(duì) 141 條pre-miR160數(shù)據(jù)進(jìn)行系統(tǒng)進(jìn)化分析 （圖1-B），結(jié)果表明，棉花 miR160 和甜橙（Citrus sinensis）、木薯（Manihot esculenta）、葡萄（Vitis vinifera）、毛果楊（Populus trichocarpa）及番木瓜（Carica papaya）的miR160具有較高的同源性。

2.2　棉花miR160的表達(dá)模式分析

葉片衰老最顯著的變化是顏色由綠變黃，且葉片顏色越黃，衰老程度越嚴(yán)重[23]。觀察CCRI 10和 Liao 4086 在 4 個(gè)時(shí)期（7 d、21 d、35 d、49 d）的子葉表型，CCRI 10子葉生長(zhǎng)35 d時(shí)顏色明顯變黃，而同時(shí)期Liao 4086子葉并無(wú)明顯的黃化表型（圖2-A）。因此，此時(shí)期CCRI 10的衰老程度比Liao 4086更為嚴(yán)重。qRT-PCR測(cè)定ghr-miR160（圖 2-B）及其靶基因GhRRP41L（圖2-C）的相對(duì)表達(dá)量表明，在35 d的Liao 4086子葉中g(shù)hr-miR160表達(dá)量上調(diào)，GhRRP41L表達(dá)量下調(diào)，說(shuō)明ghr-miR160可能通過(guò)作用于靶基因延緩子葉的衰老。

圖1　miR160多重序列比對(duì)及進(jìn)化關(guān)系分析Fig.1　Multiple sequence alignment and phylogenetic analysis of miR160 family

圖2　棉花miR160及其靶基因GhRRP41L在CCRI 10和Liao 4086不同時(shí)期子葉中的相對(duì)表達(dá)量Fig.2　qRT-PCR analysis of ghr-miR160 andGhRRP41Lin CCRI 10 and Liao 4086 at different development stages

2.3　棉花miR160的功能驗(yàn)證分析

提取轉(zhuǎn)ghr-miR160基因擬南芥株系LR-1、LR-2、LR-3和野生型擬南芥葉片的 DNA和RNA，進(jìn)行PCR驗(yàn)證和表達(dá)量分析。在DNA水平，轉(zhuǎn)ghr-miR160基因擬南芥的條帶和陽(yáng)性質(zhì)粒條帶位置相同，而陰性對(duì)照和野生型擬南芥無(wú)條帶（圖3-A），說(shuō)明ghr-miR160基因已經(jīng)插入擬南芥基因組中；在RNA水平，ghr-miR160在轉(zhuǎn)該基因擬南芥株系中表達(dá)量顯著上調(diào)45～390倍（圖3-B），說(shuō)明ghr-miR160已在擬南芥中轉(zhuǎn)錄表達(dá)。表型觀察結(jié)果表明，轉(zhuǎn)ghr-miR160基因擬南芥比野生型擬南芥要晚抽薹，晚開(kāi)花，葉片衰老延緩（圖 4-A）。 分別取生長(zhǎng) 20 d、30 d、40 d和 50 d的轉(zhuǎn)ghr-miR160基因和野生型擬南芥蓮座葉，測(cè)定擬南芥衰老相關(guān)基因SAG12的表達(dá)量，同時(shí)測(cè)定葉綠素。結(jié)果表明，在擬南芥衰老進(jìn)程中，SAG12的表達(dá)量總體呈逐漸上升趨勢(shì)，在擬南芥生長(zhǎng)的40 d和50 d，轉(zhuǎn)該基因擬南芥SAG12的表達(dá)量低于野生型擬南芥（圖4-B）。然而，擬南芥葉綠素含量隨著葉片衰老總體呈下降趨勢(shì)，在生長(zhǎng)40 d和50 d時(shí)，轉(zhuǎn)ghr-miR160基因擬南芥的葉綠素含量高于野生型擬南芥（圖4-C）。

圖3　轉(zhuǎn)基因擬南芥陽(yáng)性檢測(cè)和ghr-miR160表達(dá)量變化Fig.3　The checking of transgenicArabidopsisby PCR and qRT-PCR of ghr-miR160

3　討論

葉片衰老是植物葉片發(fā)育的最后階段，在衰老過(guò)程中，葉綠素喪失，大分子物質(zhì)如蛋白質(zhì)、脂類、核酸等被降解，導(dǎo)致光合作用減弱，同時(shí)衰老組織的營(yíng)養(yǎng)物質(zhì)被運(yùn)輸?shù)接啄劢M織和生殖器官中[24]。棉花葉片衰老是1個(gè)普遍性的問(wèn)題，發(fā)生于各國(guó)的產(chǎn)棉區(qū)，導(dǎo)致棉花減產(chǎn)降質(zhì)[4]。目前對(duì)于棉花葉片衰老的研究主要集中在生理生化指標(biāo)上；因此，研究棉花早衰過(guò)程中的關(guān)鍵基因具有重要意義。

根據(jù)Zhang等[25]對(duì)71個(gè)物種37個(gè)miRNA家族研究表明，至少存在于10個(gè)不同物種的miRNA家族是高度保守的。本研究發(fā)現(xiàn)植物miR160基因在進(jìn)化上高度保守，這與Zhang等[25]的研究一致。保守的miRNA基因往往形成復(fù)雜的調(diào)控網(wǎng)絡(luò)，功能具有多樣性[26]。在所有miR160分布中，玉米（Zea mays）、亞麻（Linum usitatissimum）、 小立碗蘚 （Physcomitrella patens）、 木薯（Manihot esculenta）中分別發(fā)現(xiàn) 13 個(gè)、10 個(gè)、9 個(gè)和 8 個(gè)；而在菊芋（Helianthus tuberosus）、向日葵（Helianthus annuus）和小麥（Triticum aestivum）中分別發(fā)現(xiàn)1個(gè)。由此可見(jiàn)，miR160在低等植物和高等植物中都存在，是1個(gè)相對(duì)古老的家族；但是各個(gè)植物之間的存在數(shù)量有明顯差別，說(shuō)明miR160基因在進(jìn)化中，有些舊的基因消失，也有些新的基因出現(xiàn)，導(dǎo)致其在基因組呈分散分布，也進(jìn)一步說(shuō)明miR160家族調(diào)控的復(fù)雜性和研究?jī)r(jià)值。

為了探究棉花miR160的表達(dá)模式，本研究利用qRT-PCR分析了CCRI 10和Liao 40864各時(shí)期子葉的ghr-miR160及靶基因GhRRP41L的表達(dá)量。ghr-miR160在Liao 4086在子葉生長(zhǎng)35 d中表達(dá)量顯著高于CCRI 10。Xu等[27]也發(fā)現(xiàn)在水稻灌漿后期，miR160在不早衰品種的葉片中優(yōu)勢(shì)表達(dá)。這表明ghr-miR160可能調(diào)控棉花葉片的衰老。

目前，對(duì)于miR160的研究主要集中于其在種子萌發(fā)、根冠發(fā)育及響應(yīng)逆境等方面的作用[12]。為了進(jìn)一步驗(yàn)證miR160在葉片衰老中的功能，用ghr-miR160過(guò)表達(dá)載體轉(zhuǎn)化擬南芥。在葉片衰老過(guò)程中，一些基因被激活，表達(dá)增強(qiáng)，此類基因?yàn)樗ダ舷嚓P(guān)基因[28]。qRT-PCR檢測(cè)擬南芥衰老相關(guān)基因SAG12表達(dá)量結(jié)果顯示，生長(zhǎng)40 d的野生型擬南芥的SAG12表達(dá)量約為轉(zhuǎn)ghr-miR160基因擬南芥的2倍，而至50 d時(shí)高達(dá)5倍。并且擬南芥生長(zhǎng)40 d、50 d時(shí)，轉(zhuǎn)ghr-miR160基因擬南芥葉綠素含量約為野生型擬南芥的18倍。以上結(jié)果表明，轉(zhuǎn)ghr-miR160基因擬南芥的衰老速度遠(yuǎn)遠(yuǎn)低于野生型擬南芥。張停停[29]發(fā)現(xiàn)miR160的靶基因ARF16突變體可使擬南芥提前衰老。以上研究表明，miR160負(fù)調(diào)控葉片衰老。

4　結(jié)論

通過(guò)對(duì)ghr-miR160進(jìn)行生物信息學(xué)分析、棉花子葉表達(dá)模式分析及過(guò)表達(dá)載體轉(zhuǎn)化擬南芥驗(yàn)證分析，結(jié)果表明：其在進(jìn)化過(guò)程中高度保守，且在35 d的早熟不早衰品種遼4086子葉中優(yōu)勢(shì)表達(dá)；轉(zhuǎn)該基因擬南芥出現(xiàn)晚抽薹、晚開(kāi)花和葉片衰老延遲現(xiàn)象。測(cè)定擬南芥葉片SAG12相對(duì)表達(dá)量和葉綠素含量結(jié)果進(jìn)一步表明，轉(zhuǎn)該基因擬南芥出現(xiàn)衰老延緩現(xiàn)象。以上結(jié)果初步證明棉花miR160基因可能對(duì)葉片衰老過(guò)程進(jìn)行負(fù)調(diào)控，為miRNA在葉片衰老進(jìn)程中發(fā)揮調(diào)控作用機(jī)制提供了新的思路，也為棉花分子育種提供了潛在的候選基因。

圖4　轉(zhuǎn)ghr-miR160基因擬南芥衰老延緩表型及相關(guān)指標(biāo)測(cè)定Fig.4　PhenotypefordelayedleafsenescenceandresultsofrelateddetectioninArabidopsisoverexpressedghr-miR160

參考文獻(xiàn)：

[1]Bartel D P,Bartel B.MicroRNAs and their regulatory roles in plant[J].Annual Review of Plant Biology,2006,57:19-53.

[2]Bartel D P.MicroRNAs:Genomics,biogenesis,mechanism,and function[J].Cell,2004,116(2):281-297.

[3]Mallory A C,Vaucheret H.Functions of microRNAs and related small RNAs in plants[J].Nature Genetics,2006,38(Suppl):31-36.

[4]王丹.棉花葉片衰老相關(guān)基因的數(shù)字化表達(dá)譜分析[D].楊凌:西北農(nóng)林科技大學(xué),2011.Wang Dan.Identification and analysis of senescence-associated genes in cotton leaves[D].Yangling:Northwest A&F University,2011.

[5]Abeles F B,Dunn L J,Morgens P,et al.Induction of 33-kD and 60-kD peroxidases during ethylene-induced senescence of cucumber cotyledons[J].Plant Physiology,1988,87(3):609-615.

[6]Zeevaart J A,Creelman R A D.Metabolism and physiology of abscisic acid[J].Annual Review of Plant Biology,1988,39:439-473.

[7]He Y,Fukushige H,Hildebrand D F,et al.Evidence supporting a role of jasmonic acid inArabidopsisleaf senescence[J].Plant Physiology,2002,128(3):876-884.

[8]Van der Graaff E,Schwacke R,Schneider A,et al.Transcription analysis ofArabidopsismembrane transporters and hormone pathways during developmental and induced leaf senescence[J].Plant Physiology,2006,141(2):776-792.

[9]Balibrea L M E,Gonzalez G M C,Fatima T,et al.Extracellular invertaseis an essential component of cytokinin mediated delay of senescence[J].Plant Cell,2004,16(5):1276-1287.

[10]Lim P O,Kim H J,Nam H G.Leaf senescence[J].Annual Review of Plant Biology,2007,58:115-136.

[11]Li J R,Yu K,Wei J R,et al.Gibberellin retards chlorophyll degradation during senescence ofParis polyphylla[J].Biologia Plantarum,2010,54(2):395-399.

[12]Rhoades M W,Reinhart B J,Lim L P,et al.Prediction of plant microRNA targets[J].Cell,2002,110(4):513-520.

[13]李小平,曾慶發(fā),張根生,等.大豆 microRNA基因Gm-MIR160A負(fù)調(diào)控植物葉片衰老進(jìn)程[J].廣西植物,2015,35(1):84-91.Li Xiaoping,Zeng Qingfa,Zhang Gensheng,et al.Gm-MIR160A,a class of soybean microRNA gene,negetively regulates progress of leaf senescence[J].Guihaia,2015,35(1):84-91.

[14]楊松楠,王姣,陳宗祥,等.普通野生稻miR160f的克隆和功能分析[J].生物技術(shù)通報(bào),2014(11):114-118.Yang Songnan,Wang Jiao,Chen Zongxiang,et al.Clone and function analysis miR160f in common wild rice(Oryza rufipogonGriff.)[J].Biotechnology Bulletin,2014(11):114-118.

[15]Lin Yuling,Lai Zhongxiong,Tian Qilin,et al.Endogenous target mimics down-regulate miR160 mediation ofARF10,-16,and-17 cleavage during somatic embryogenesis inDimocarpus longanLour.[J/OL].Frontiers in Plant Science,2015,6:956(2015-11-25)[2016-12-01].https://doi.org/10.3389/fpls.2015.00956.

[16]Hendelmen A,Buxdorf K,Stav R,et al.Inhibition of lamina outgrowth followingSolanum lycopersicumAUXIN RESPONSE FACTOR10(SIARF10)derepression[J].Plant Molecular Biology,2012,78(6):561-576.

[17]Wang Kunbo,Wang Zhiwen,Li Fuguang,et al.The draft genome of a diploid cotton[J].Nature Genetics,2012,44(10):1098-1103.

[18]Li Fuguang,Fan Guangyi,Wang Kunbo,et al.Genome sequence of the cultivated cottonGossypium arboreum[J].Nature Genetics,2014,46(6):567-572.

[19]Li Fuguang,Fan Guangyi,Lu Cairui,et al.Genome sequence of cultivated upland cotton(Gossypium hirsutumTM-1)provides insights into genome evolution[J].Nature Biotechnology,2015,33(5):524-530.

[20]李文超,趙淑清.人工microRNAs對(duì)擬南芥At1g13770和At2g23470基因的特異沉默[J].遺傳,2012,34(3):348-355.Li Wenchao,ZhaoShuqing.Specific gene silencing ofAt1g13770 andAt2g23470 by artificial microRNAs inArabidopsis[J].Hereditas(Beijing),2012,34(3):348-355.

[21]宋國(guó)立,崔榮霞.改良CTAB法快速提取棉花DNA[J].棉花學(xué)報(bào),1998,10(5):273-275.Song Guoli,Cui Rongxia.A rapid improved CTAB method for extraction of cotton genomic DNA[J].Cotton Science,1998,10(5):273-275.

[22]Bechtold N.In planta agrobacteriummediated gene transfer by infiltration of adultArabidopsis thalianaplants[J].Comptes Rendus De Lacademie Des Sciences Serie Sciences De La Vie,1993,316(10):1194-1199.

[23]張寶來(lái).水稻葉片衰老的研究進(jìn)展[J].天津農(nóng)業(yè)科學(xué),2013,19(4):19-24.Zhang Baolai.Advances of research on leaf senescence in rice[J].Tianjin Agricultural Sciences,2013,19(4):19-24.

[24]Lim P O,Kim Y,Breeze E,et al.Overexpression of a chromatin architecture-controlling AT-hook protein extends leaf longevity and increases the post-harvest storage life of plants[J].The Plant Journal,2007,52(6):1140-1153.

[25]Zhang Baohong,Pan Xiaoping,Charles H C,et al.Conservation and divergence of plant microRNA genes[J].The Plant Journal,2006,46(2):243-259.

[26]魏強(qiáng),梁永宏,李廣林.植物miRNA的進(jìn)化[J].遺傳,2013,35(3):315-323.Wei Qiang,Liang Yonghong,Li Guanglin.Evolution of miRNA in plant[J].Hereditas(Beijing),2013,35(3):315-323.

[27]Xu Xiangbin,Bai Haiqi,Liu Chaoping,et al.Genome-wide analysis of microRNAs and their target genes related to leaf senescence of rice[J/OL].PLoS ONE,2014,9(12):e114313(2014-12-05)[2016-12-01].https://doi.org/10.1371/journal.pone.0114313.

[28]Gan S,Amasino R M.Making sense of senescence(molecular genetic regulation and manipulation of leaf senescence)[J].Plant Physiology,1997,113(2):313-319.

[29]張停停.擬南芥生長(zhǎng)素響應(yīng)轉(zhuǎn)錄因子ARF16生物學(xué)功能的初步探究[D].淮北:淮北師范大學(xué),2013.Zhang Tingting.The function of auxin response transcription factorARF16 inArabidopsis thaliana[D].Huaibei:Huaibei Normal University,2013.