孫廣鑫,欒雨時,崔娟娟

大連理工大學生命科學與技術(shù)學院,大連 116024

MicroRNA(miRNA)普遍存在于生物體內(nèi),是由內(nèi)源基因編碼的長度為22nt左右的單鏈非編碼小分子 RNA[1],可與 RNA誘導沉默復合體(RNA-induced silencing complex,RISC) 結(jié)合識別其靶mRNA。植物中miRNA序列高度保守并呈現(xiàn)出獨特的表達時序性和組織特異性。在植物細胞核中,內(nèi)源miRNA基因先形成前體轉(zhuǎn)錄本pri-miRNA,然后被加工形成較為穩(wěn)定的初級莖環(huán)結(jié)構(gòu),再進一步形成具有發(fā)夾結(jié)構(gòu)的 pre-miRNA,最后經(jīng) DCL1酶剪切形成 miRNA/miRNA*雙鏈復合體,從細胞核進入細胞質(zhì)形成成熟的 miRNA[2],與 RISC結(jié)合后作用于靶基因,發(fā)揮其翻譯阻遏或者降解作用[3]。miRNA參與植物生長發(fā)育的多種生理生化過程以及各種生物與非生物脅迫反應[4]等。

番茄作為重要的園藝作物,不僅具有很高的經(jīng)濟價值,而且在生命科學研究中也占有舉足輕重的地位。由于番茄具有易雜交、繁殖系數(shù)高、遺傳資源豐富、突變體庫較多、遺傳圖譜廣泛、遺傳轉(zhuǎn)化體系高效穩(wěn)定等特性,使其成為作物研究中的模式植物之一。在植物與病原物互作的研究方面,由于番茄基因組復雜且病原種類豐富,陸續(xù)得到了眾多的理論信息[5]。2009年,Feng等[6]率先檢測了與病毒相關(guān)的番茄 miRNA,其后的研究表明 miRNA在黃瓜花葉病毒(Cucumber mosaic virus,CMV)、番茄不孕病毒(Tomato aspermy virus,TAV)、番茄卷葉新德里病毒(Tomato leaf curl new delhi virus,ToLCNDV)、番茄灰霉病等病害的防御和發(fā)生中均發(fā)揮著重要的調(diào)控作用。因此,對miRNA的抗逆功能研究可為番茄種質(zhì)資源改良提供一定的理論依據(jù)。

目前對 miRNA的挖掘主要有直接克隆(Direct cloning)、正向遺傳學(Forward genetics)篩選和生物信息學(Bioinformatics)分析等方法。直接克隆法比較簡單,第一條植物miRNA就是用該方法發(fā)現(xiàn)的,很多實驗室已經(jīng)利用這種方法從不同生物中發(fā)現(xiàn)了眾多miRNA[7～9],該方法更傾向于發(fā)現(xiàn)中、高量表達的miRNA,但很難發(fā)現(xiàn)一些表達量低、組織特異性強或受誘導表達的miRNA。近年來,高通量測序技術(shù)的發(fā)展極大地完善了該方法,并使得大規(guī)模鑒定miRNA成為可能[10,11],很多植物的miRNA都是通過高通量測序獲得的[12～14]。正向遺傳學篩選是在獲得突變個體后,將突變基因分離出來,從而得到一個新的產(chǎn)物。最初的lin-4和let-7就是利用正向遺傳學篩選發(fā)現(xiàn)的[15,16],該方法可以清楚地了解miRNA的功能,但效率很低,只能鑒定到少量的miRNA,因此,正向遺傳學篩選方法無法成為發(fā)現(xiàn)miRNA的主要手段,只是在 miRNA研究早期使用較多。生物信息學分析方法可篩選出經(jīng)各種測序方法獲得的大量數(shù)據(jù)中的miRNA,它能夠彌補直接克隆法的不足,發(fā)現(xiàn)一些豐度低、組織特異性強或受誘導表達的miRNA。隨著一些物種基因組序列、大量EST(Expressed sequence tag)和GSS(Genome survey sequences)等序列的公布,也可以不通過測序,而是根據(jù)miRNA的保守性和獨有的特征,利用生物信息學方法在基因組、EST和 GSS中預測 miRNA前體,進而找到保守的 miRNA,但通過該方法獲得的miRNA還需經(jīng)過實驗驗證以確定其真實性[17,18]。生物信息學分析方法快速、簡單、高效,但是很難發(fā)現(xiàn)與現(xiàn)有miRNA無同源性的新的miRNA。

當番茄受病原物侵染時,體內(nèi)的多種miRNA及相關(guān)靶基因的表達量都會發(fā)生變化,說明miRNA在番茄響應生物脅迫中發(fā)揮著關(guān)鍵的調(diào)控作用[19,20]。為了篩選出番茄體內(nèi)與致病密切相關(guān)的miRNA,本研究采用生物信息學手段,從Sanger miRNA數(shù)據(jù)庫——miRBase(Release 20:June 2013,http://www.mirbase.org/)和已報道的文獻中找到34個miRNA的成熟序列,預測了它們的靶基因,利用Cytoscope軟件構(gòu)建這些 miRNA之間的調(diào)控網(wǎng)絡,從中篩選出miR169、miR482、miR5300、miR6024、miR6026和 miR6027,進行了進一步的分析及實時定量 PCR驗證,為深入研究miRNAs的作用機制奠定基礎(chǔ)。

1 材料和方法

1.1 材料

野生醋栗番茄(Solanum pimpinellifolium)材料L3708為本實驗室保藏; 致病疫霉菌種P12103為西北農(nóng)林科技大學單衛(wèi)星教授惠贈。通過NCBI、Web of Knowledge及文獻檢索等途徑,統(tǒng)計出番茄中與致病相關(guān)的 miRNA共 34個(表1),其中 miR156、miR159、miR160、miR162、miR166、miR167、miR169、miR171、miR172、miR319、miR395、miR397、miR399、miR482、miR5300、miR6022、miR6023、miR6024、miR6026和 miR6027是在miRBase中已經(jīng)登錄的。利用Primer5軟件分析設(shè)計其中 5個 miRNA的上游引物序列及內(nèi)參基因Actin(GenBank:U60478.1)的上、下游引物序列(表 2)。

1.2 方法

1.2.1 致病相關(guān)miRNA序列的獲得及靶基因預測

從 Sanger miRNA數(shù)據(jù)庫——miRBase(http://www.mirbase.org/search.shtml)(release 20)和已報道的文獻及其附件中收集到的34個miRNA全部家族成員的成熟序列。

將上述 miRNA利用軟件 psRNATarget(http://plantgrn.noble.org/psRNATarget) 進行靶基因預測。在Solanum Lycopersicon(tomato)、transcript、cDNA library version2.3數(shù)據(jù)庫中預測番茄miRNA的靶基因。預測參數(shù)除Maximum expectation設(shè)為5.0外,其余均使用程序默認參數(shù)。

1.2.2 致病密切相關(guān)miRNA的篩選

利用Cytoscope軟件構(gòu)建病原誘導的miRNA及其靶基因的調(diào)控網(wǎng)絡。根據(jù)各靶基因的個數(shù)確定對應圓點的大小,從中篩選出與致病密切相關(guān)的miRNA作為主要的研究和分析對象。

表1 番茄中與致病相關(guān)的34個miRNA

表2 設(shè)計的5個miRNA的上游引物及Actin的上、下游引物序列

1.2.3 致病密切相關(guān)miRNA的靶基因功能分析

根據(jù)上述靶基因預測結(jié)果,針對篩選出的若干miRNA的靶基因名稱和 GO注釋號在 Uni-protKB(http://www.uniprot.org/help/uniprotkb)蛋白數(shù)據(jù)庫中分析各靶基因的功能。統(tǒng)計出與病害相關(guān)的靶基因。

1.2.4 致病密切相關(guān)miRNA的啟動子分析

針對篩選出的若干 miRNA 的基因 5′端上游1 500 bp的啟動子片段,利用 PlantCARE(http://bioinformatics.psb.ugent.be/webtools/plantcare/html/)在線程序進行上游啟動子預測及分析。找出與病害相關(guān)的啟動子元件。

1.2.5 致病密切相關(guān)miRNA的實時定量PCR驗證

選取長勢相同的番茄五葉期幼苗為試材,以106個孢子/mL的致病疫霉孢子懸浮液噴施葉片,分別于處理后 0 h、3h、6 h、12 h、24 h、48 h、72 h、96 h時取樣,用RNAiso plus(TaKaRa,Japan)試劑盒提取總RNA,再用miRNA反轉(zhuǎn)錄試劑盒One Step PrimeScript miRNA cDNA Synthesis Kit(TaKaRa,Japan)對樣品中的 miRNA進行加尾,加上一個長為60 bp左右的接頭,通過實時定量PCR進行驗證。本研究使用的各miRNA的上游引物見表2,下游引物為試劑盒中提供的通用引物Uni-miR qPCR Primer。反應體系為 25 μL,包括：2 × SYBR?Premix ExTaqTM12.5 μL,cDNA 模板 2 μL,正反向引物各 1 μL,ddH2O 8.5 μL。反應條件：95℃預變性10 s; 95℃變性5 s; 60℃復性20 s,40個循環(huán)。用番茄Actin作為內(nèi)參基因,熒光數(shù)據(jù)在每個循環(huán)的復性末期采集。對得到的結(jié)果用 2-DDCt法分析繪出柱狀圖,并對實時定量PCR產(chǎn)物進行4%瓊脂糖凝膠電泳檢測。

2 結(jié)果與分析

2.1 致病密切相關(guān)miRNA的靶基因預測及篩選

整理分析psRNATarget在線預測出的miRNA及其靶基因。先構(gòu)建miRNA與靶基因間的Cytoscope調(diào)控網(wǎng)絡。發(fā)現(xiàn)在整個網(wǎng)絡中有多個miRNA攜帶它的靶基因單獨存在,還有一部分miRNA是通過它們的靶基因與少數(shù)幾個miRNA相關(guān)聯(lián),而只有一部分miRNA(miR156、miR157、miR159、miR167、miR169、miR171、miR319、miR403、miR482、miR5300、miR6024、miR6026和 miR6027)是通過它們的靶基因與多個 miRNA相關(guān)聯(lián)形成了致密網(wǎng)絡。再以這13個分布較集中的 miRNA為研究對象構(gòu)建得到它們相互間的Cytoscope調(diào)控網(wǎng)絡(圖1)。從圖1可知miR169、miR482、miR5300、miR6024、miR6026和miR6027的靶基因個數(shù)較多且與其他miRNA的關(guān)聯(lián)性較大(至少有4個相關(guān)聯(lián)的miRNA)。

2.2 致病密切相關(guān)miRNA的靶基因功能分析

圖1 致病密切相關(guān)各 miRNA之間的調(diào)控網(wǎng)絡(圓的大小與靶基因數(shù)目呈正相關(guān))

用psRNATarget在線軟件預測miR169、miR482、miR5300、miR6024、miR6026和 miR6027的靶基因,經(jīng)篩選后各獲得多條靶基因。根據(jù)各靶基因的名稱和GO注釋號在Uni-protKB蛋白數(shù)據(jù)庫中查詢每條靶基因的功能,可將它們分為 5類：(1)編碼與脅迫響應相關(guān)的蛋白質(zhì);(2)編碼與代謝相關(guān)的酶蛋白;(3)編碼與跨膜運輸相關(guān)的蛋白質(zhì);(4)編碼轉(zhuǎn)錄因子;(5)未知蛋白。統(tǒng)計發(fā)現(xiàn) miR169、miR482、miR5300、miR6024、miR6026和 miR6027編碼與脅迫響應相關(guān)蛋白質(zhì)的靶基因分別占各自總靶基因的22.45%、80.30%、71.43%、95.45%、59.09%和 100.00%?？梢?預測到的靶基因大部分與番茄的脅迫響應相關(guān),其中有很多參與致病過程(表3)。

表3 與致病密切相關(guān)的各miRNA的靶基因

2.3 致病密切相關(guān)miRNA的啟動子分析

應用PlantCARE在線程序?qū)iR169、miR482、miR5300、miR6024、miR6026和miR6027前體基因5′端上游1 500 bp的序列進行啟動子預測分析。結(jié)果表明,除 miR6027外,其他各成員的啟動子區(qū)域不僅存在典型的具有轉(zhuǎn)錄起始功能的 TATA-box和CAAT-box,還存在多種脅迫相關(guān)的順式作用元件(包括激素響應、光響應、生物脅迫響應、厭氧響應、干旱響應和溫度響應元件等)以及與生物脅迫相關(guān)的啟動子元件 TC-rich repeats; 同時 miR169、miR482、miR5300和miR6026中還存在與生物脅迫相關(guān)的啟動子元件Box-W1。由于不能確定miR6027在染色體上的具體位置,所以未對其進行啟動子分析。

2.4 致病密切相關(guān)miRNA的實時定量PCR驗證

5個miRNA的實時定量PCR檢測結(jié)果見圖2。致病疫霉處理后,番茄葉片中的miR169、miR482、miR6024、miR6026和 miR6027的表達量均在 3 h內(nèi)迅速降低,在24 h時普遍有所回升,在48 h時又降到最低點。以上結(jié)果說明這5個miRNA都能對致病疫霉處理做出響應。且通過瓊脂糖凝膠電泳檢測(圖 3)表明,miR169、miR482、miR6024、miR6026和miR6027的擴增產(chǎn)物都在80 bp左右,說明符合預期長度。

圖2 致病疫霉處理后5個miRNA的表達特性

圖3 實時定量PCR產(chǎn)物的瓊脂糖凝膠電泳檢測

3 討 論

近年來,miRNA作為基因調(diào)控因子備受矚目,已迅速成為生命科學領(lǐng)域的研究熱點。眾所周知miRNA是植物對脅迫響應過程中不可或缺的調(diào)控因子,它通過互補作用靶向mRNA抑制或降解相應基因的表達,參與調(diào)控生物體內(nèi)的各項生命活動。本研究采用生物信息學分析方法構(gòu)建 Cytoscape網(wǎng)絡圖,欲篩選與致病密切相關(guān)的 miRNA,對選出的miR169、miR482、miR5300、miR6024、miR6026和miR6027進行靶基因功能分析和啟動子分析。結(jié)果發(fā)現(xiàn),這6個miRNA的靶基因和啟動子基本都與生物脅迫相關(guān),表明它們極可能參與了番茄對脅迫響應的調(diào)控過程。為了探究各miRNA在致病疫霉處理下的表達情況,本文又對其中5個miRNA進行了表達特性分析。發(fā)現(xiàn)這5個miRNA的表達量都有顯著變化,且趨勢基本一致,均在48 h時降到最低水平,說明這5個miRNA都能對致病疫霉的接種做出響應,且它們之間可能存在某些相似的通路。由于miR5300的引物特異性不理想,所以本文未對其進行實驗分析。究竟在脅迫下各miRNA與其靶基因存在怎樣的相互關(guān)系,后者又是如何表達的,這些問題還有待于進一步的實驗驗證。

番茄作為重要的栽培植物,在受到真菌、細菌、病毒以及線蟲和昆蟲的侵害時常造成大量減產(chǎn)[37],甚至絕收,給農(nóng)業(yè)生產(chǎn)帶來嚴重影響。研究發(fā)現(xiàn),植物抗病(R)基因能夠識別病原物并引起過敏性反應,它作為防衛(wèi)系統(tǒng)的一員,在沒有病原物侵染時表達量很低,但當病原物侵染后,能夠識別激發(fā)子,誘導抗病基因的過量表達,并激活抗病信號傳導途徑,產(chǎn)生防衛(wèi)反應[38],從而提高作物對病害的抗性。所以,克隆抗性基因?qū)ρ芯恐参锖筒≡锘プ?、提高植物抗病性具有重要意義。近年來已從不同植物中克隆出了各種R基因。多數(shù)R基因具有NBS-LRR結(jié)構(gòu)域,其中NBS(Nucleotide binding site)是抗病基因核苷酸結(jié)合位點。而NBS-LRR編碼基因是抗病基因中的最大家族,大約 3/4的抗病基因都來自于這一家族。根據(jù)N端信號肽的不同,又將NBS-LRR基因分為 TIR-NBS-LRR和CC-NBS-LRR兩類[39,40]。這些抗病基因大量存在于植物基因組中,它的作用對象囊括了所有病原物[41]。經(jīng)靶基因功能分析發(fā)現(xiàn)篩選出的6個miRNA的靶基因中均含有NBS-LRR抗性蛋白,推測植物在受到病原物侵染時,這 6個miRNA的表達量也會發(fā)生相應的變化,從而說明了我們篩選的miRNA的確是與致病密切相關(guān)的。并且已有實驗證明,當番茄植株受到CMV、TAV、Tolcv和灰霉病菌侵染時,用基因芯片技術(shù)檢測到miR169表達量逐漸上調(diào),而其靶基因的表達量下調(diào),說明miR169參與了病原物的脅迫響應。miR482、miR6024、miR6026和 miR6027也被預測到參與植物的免疫應答過程,靶向植物的免疫受體[24]。

由于啟動子控制著基因表達的起始時間和表達程度,而某些miRNA的表達又具有時空特異性和組織特異性。所以,對啟動子的研究將成為miRNA表達特性分析的重要工具[42]。研究發(fā)現(xiàn),大部分miRNA都含有獨立的啟動子,是 RNA聚合酶Ⅱ的產(chǎn)物[43],經(jīng)過轉(zhuǎn)錄而形成。所以,它的表達很有可能受到精確的調(diào)控。為進一步剖析miRNA在脅迫下的表達模式,本研究利用 PlantCARE[44]對番茄miR169、miR482、miR5300、miR6024 和 miR6026基因上游啟動子區(qū)域進行分析,發(fā)現(xiàn)在這些區(qū)域中存在共同的植物病原物侵染和損傷響應順式作用元件TC-rich repeats和Box-W1。表明這些生物脅迫相關(guān)miRNA可能參與了相應的逆境脅迫過程,也暗示植物在響應生物脅迫時可能存在某些共同的通路。研究發(fā)現(xiàn)Box-W1存在于如芪合成酶(STS)和病程相關(guān)蛋白PR-1、Ypr-10中[45,46],位于多病原體侵染和非生物脅迫響應基因 VvALDH2B8的啟動子區(qū)域,能特異的結(jié)合WRKY轉(zhuǎn)錄因子[47],參與脅迫響應基因的表達調(diào)控。所以,這些脅迫響應順式元件的存在為 miRNA參與植物抗逆反應提供了證據(jù)。此外,番茄作為一種模式植物,挖掘和分析其與病害相關(guān)的 miRNA還有助于揭示植物與病原物互作的分子機制。

[1]Carrington JC,Ambros V.Role of microRNAs in plant and animal development.Science,2003,301(5631):336–338.

[2]Jones-Rhoades MW,Bartel DP,Bartel B.MicroRNAs and their regulatory roles in plants.Annu Rev Plant Biol,2006,57:19–53.

[3]金龍國,王川,劉進元.植物 MicroRNA.中國生物化學與分子生物學報,2006,22(8):609–614.

[4]Kulcheski FR,De Oliveira LFV,Molina LG,Almer?o MP,Rodrigues FA,Marcolino J,Barbosa JF,Stolf-Moreira R,Nepomuceno AL,Marcelino-Guimar?es FC,Abdelnoor RV,Nascimento LC,Carazzolle MF,Pereira GAG,Margis R.Identification of novel soybean microRNAs involved in abiotic and biotic stresses.BMC Genomics,2011,12(1):307.

[5]Arie T,Teraoka T,Takahashi H,Kodama M.Tomato as a model plant for plant-pathogen interactions.Plant Biotechnology,2007,24(1):135–147.

[6]Feng JL,Wang K,Liu X,Chen SN,Chen JS.The quantification of tomato microRNAs response to viral infection by stem-loop real-time RT-PCR.Gene,2009,437(1–2):14–21.

[7]Sunkar R,Girke T,Jain PK,Zhu JK.Cloning and characterization of microRNAs from rice.Plant Cell,2005,17(5):1397–1411.

[8]Tang GL,Reinhart BJ,Bartel DP,Zamore PD.A biochemical framework for RNA silencing in plants.Genes Dev,2003,17(1):49–63.

[9]Llave C,Xie ZX,Kasschau KD,Carrington JC.Cleavage ofScarecrow-likemRNA targets directed by a class ofArabidopsismiRNA.Science,2002,297(5589):2053–2056.

[10]Lu C,Tej SS,Luo SJ,Haudenschild CD,Meyers BC,Green PJ.Elucidation of the small RNA component of the transcriptome.Science,2005,309(5740):1567–1569.

[11]Nakano M,Nobuta K,Vemaraju K,Tej SS,Skogen JW,Meyers BC.Plant MPSS databases:signature-based transcriptional resources for analyses of mRNA and small RNA.Nucleic Acids Res,2006,34(S1):D731–D735.

[12]Zhao CZ,Xia H,Frazier TP,Yao YY,Bi YP,Li AQ,Li MJ,Li CS,Zhang BH,Wang XJ.Deep sequencing identifies novel and conserved microRNAs in peanuts(Arachis hypogaeaL.).BMC Plant Biol,2010,10(1):3.

[13]Li BS,Qin YR,Duan H,Yin WL,Xia XL.Genome-wide characterization of new and drought stress responsive microRNAs inPopulus euphratica.J Exp Bot,2011,62(11):3765–3779.

[14]Zhou ZS,Zeng HQ,Liu ZP,Yang ZM.Genome-wide identification ofMedicago truncatulamicroRNAs and their targets reveals their differential regulation by heavy metal.Plant Cell Environ,2012,35(1):86–89.

[15]Lee RC,Feinbaum RL,Ambros V.The C.elegans heterochronic genelin-4encodes small RNAs with antisense complementarity tolin-14.Cell,1993,75(5):843– 854.

[16]Reinhart BJ,Slack FJ,Basson M,Pasquinelli AE,Bettinger JC,Rougvie AE,Horvitz HR,Ruvkun G.The 21 nucleotidelet-7RNA regulates developmental timing inCaenorhabditis elegans.Nature,2000,403(6772):901–906.

[17]張志明,宋銳,彭華,羅茂,沈亞歐,劉麗,趙茂俊,潘光堂.用生物信息學挖掘玉米中的 microRNAs及其靶基因.作物學報,2010,36(8):1324–1335.

[18]Zhang JG,Zeng R,Chen JS,Liu X,Liao QS.Identification of conserved microRNAs and their targets fromSolanum lycopersicumMill.Gene,2008,423(1):1–7.

[19]Feng JL,Lai LY,Lin RH,Jin CZ,Chen JS.Differential effects ofcucumber mosaic virussatellite rnas in the perturbation of microrna-regulated gene expression in tomato.Mol Biol Rep,2012,39(1):775–784.

[20]Diermann N,Matousek J,Junge M,Riesner D,Steger G.Characterization of plant mirnas and small rnas derived from potato spindle tuber viroid(pstvd) in infected tomato.Biol Chem,2010,391(12):1379–1390.

[21]Feng JL,Lin RH,Chen JS.Alteration of tomato micro-RNAs expression during fruit development uponCucumber mosaic virusandTomato aspermy virusinfection.Mol Biol Rep,2013,40(5):3713-3722.

[22]Feng JL,Lai LY,Lin RH,Jin CZ,Chen JS.Differential effects ofCucumber mosaic virussatellite RNAs in the perturbation of microRNA-regulated gene expression in tomato.Mol Biol Rep,2012,39(1):775–784.

[23]Feng JL,Liu X,Lai LY,Chen JS.Spatio-temporal expression of miRNAs in tomato tissues uponCucumber mosaic virusandTomato aspermy virusinfections.Acta Biochim Biophys Sin,2011,43(4):258–266.

[24]Naqvi AR,Haq QMR,Mukherjee SK.MicroRNA profiling ofTomato leaf curl new delhi virus(tolcndv) infected tomato leaves indicates that deregulation of mir159/319 and mir172 might be linked with leaf curl disease.Virol J,2010,7(1):281.

[25]Vu TV,Choudhury NR,Mukherjee SK.Transgenic tomato plants expressing artificial microRNAs for silencing the pre-coat and coat proteins of a begomovirus,Tomato leaf curl New Delhi virus,show tolerance to virus infection.Virus Res,2013,172(1-2):35–45.

[26]Várallyay E,Válóczi A,ágyi A,Burgyán J,Havelda Z.Plant virus-mediated induction of miR168 is associated with repression of ARGONAUTE1 accumulation.EMBO J,2010,29(20):3507–3519.

[27]Stav R,Hendelman A,Buxdorf K,Arazi T.Transgenic expression ofTomato bushy stunt virussilencing suppressor P19 via the pOp/LhG4 transactivation system induces viral-like symptoms in tomato.Virus Genes,2010,40(1):119–129.

[28]Martín R,Arenas C,Daròs JA,Covarrubias A,Reyes JL,Chua NH.Characterization of small RNAs derived fromCitrus exocortis viroid(CEVd) in infected tomato plants.Virology,2007,367(1):135–146.

[29]Pilcher RLR,Moxon S,Pakseresht N,Moulton V,Manning K,Seymour G,Dalmay T.Identification of novel small RNAs in tomato(Solanum lycopersicum).Planta,2007,226(3):709–717.

[30]Feng JL,Wang K,Liu X,Chen SN,Chen JS.The quantification of tomato microRNAs response to viral infection by stem-loop real-time RT-PCR.Gene,2009,437(1-2):14–21.

[31]Feng JL,Wang YW,Lin RH,Chen JS.Altered expression of microRNAs and target mRNAs in tomato root and stem tissues upon different viral infection.J Phytopathol,2013,161(2):107–119.

[32]Cillo F,Mascia T,Pasciuto MM,Gallitelli D.Differential effects of mild and severeCucumber mosaic virusstrains in the perturbation of microRNA-regulated gene expression in tomato map to the 3′ sequence of RNA 2.Mol Plant Microbe Interact,2009,22(10):1239–1249.

[33]Naqvi AR,Choudhury NR,Mukherjee SK,Haq QMR.In silicoanalysis reveals that several tomato micro-RNA/microRNA* sequences exhibit propensity to bind toTomato leaf curl virus(ToLCV) associated genomes and most of their encoded open reading frames(ORFs).Plant Physiol Biochem,2011,49(1):13–17.

[34]Shivaprasad PV,Chen HM,Patel K,Bond DM,Santos BACM,Baulcombe DC.A microRNA superfamily regulates nucleotide binding site–leucine-rich repeats and other mRNAs.Plant Cell,2012,24(3):859–874.

[35]Goswami S,Sahana N,Pandey V,Doblas P,Jain RK,Palukaitis P,Canto T,Praveen S.Interference in plant defense and development by non-structural protein NSs ofGroundnut bud necrosis virus.Virus Res,2012,163(1):368–373.

[36]Li F,Pignatta D,Bendix C,Brunkard JO,Cohn MM,Tung J,Sun H,Kumar P,Baker B.MicroRNA regulation of plant innate immune receptors.Proc Natl Acad Sci USA,2012,109(5):1790–1795.

[37]van Kan JAL.Licensed to kill:the lifestyle of a necrotrophic plant pathogen.Trends Plant Sci,2006,11(5):247–253.

[38]Thurau T,Kifle S,Jung C,Cai DG.The promoter of the nematode resistance gene Hs1pro-1activates a nematode-responsive and feeding site-specific gene expression in sugar beet(Beta vulgarisL.) andArabidopsis thaliana.Plant Mol Biol,2003,52(3):643–660.

[39]Dangl JL,Jones JDG.Plant pathogens and integrated defence responses to infection.Nature,2001,411(6839):826–833.

[40]van Der Biezen EA,Jones JDG.The NB-ARC domain:a novel signalling motif shared by plant resistance gene products and regulators of cell death in animals.Curr Biol,1998,8(7):R226–R227.

[41]He LM,Du CG,Covaleda L,Xu ZY,Robinson AF,Yu JZ,Kohel RJ,Zhang HB.Cloning,characterization,and evolution of the NBS-LRR-encoding resistance gene analogue family in polyploid cotton(Gossypium hirsutumL.).Mol Plant Microbe Interact,2004,17(11):1234–1241.

[42]Sunkar R,Kapoor A,Zhu JK.Posttranscriptional induction of two Cu/Zn superoxide dismutase genes inArabidopsisis mediated by downregulation of miR398 and important for oxidative stress tolerance.Plant Cell,2006,18(8):2051–2065.

[43]Johnston RJ Jr,Hobert O.A microRNA controlling left/right neuronal asymmetry inCaenorhabditis elegans.Nature,2003,426(6968):845–849.

[44]Lescot M,Dehais P,Thijs G,Marchal K,Moreau Y,van de Peer Y,Rouz P,Rombauts S.PlantCARE,a database of plantcis-acting regulatory elements and a portal to tools forin silicoanalysis of promoter sequences.Nucleic Acids Res,2002,30(1):325–327.

[45]Pühringer H,Moll D,Hoffmann-Sommergruber K,Watillon B,Katinger H,Machado MLD.The promoter of an appleYpr10gene,encoding the major allergen Mal d 1,is stress-and pathogen-inducible.Plant Sci,2000,152(1):35–50.

[46]Pape S,Thurow C,Gatz C.The ArabidopsisPR-1promoter contains multiple integration sites for the coactivator NPR1 and the repressor SNI1.Plant Physiol,2010,154(4):1805–1818.

[47]Rushton PJ,Torres JT,Parniske M,Wernert P,Hahlbrock K,Somssich IE.Interaction of elicitor-induced DNA-binding proteins with elicitor response elements in the promoters of parsleyPR1genes.EMBO J,1996,15(20):5690–5700.