位明明++李維國

摘 要 miRNAs（microRNA，微RNA）作為植物體內(nèi)一類重要的內(nèi)源性非編碼小分子RNA，能夠顯著調(diào)控細(xì)胞發(fā)育的多個生物學(xué)過程，對植物生長起重要調(diào)控作用。諸多研究表明：在植物生長發(fā)育過程中miRNA的表達(dá)具有組織特異性和時空特異性。因此，檢測和分析不同組織或樣品中miRNA及其靶基因的種類與表達(dá)特征，可為研究miRNA的生物學(xué)功能提供重要信息和依據(jù)。本文綜述了近年來植物中miRNA及其靶基因的檢測與鑒定方法及其存在的優(yōu)缺點(diǎn)，為進(jìn)一步開展miRNA及其靶基因的功能研究提供借鑒和參考。

關(guān)鍵詞 miRNA ；靶基因 ；檢測方法

分類號 Q522

Detection and Identification of miRNAs and Their Target Genes in Plants

WEI Mingming LI Weiguo

（State Centre for Rubber Breeding / Ministry of Agriculture Key Laboratory of Biology and Genetic Resources of Rubber Tree / Rubber Research Institute， CATAS，Danzhou，Hainan 571737）

Abstract miRNAs， an important class of endogenous non-coding small RNA molecules in plants， can significantly regulate multiple biological processes of cell development and play an important role in plant growth regulation. Many studies show that the expression of miRNA in the process of plant growth has spatial and temporal specificity. Therefore，detection and analysis the characterization of miRNAs and their target genes in different tissues or samples can provide important information and basis for studying the biological function of miRNA. This paper discusses the detection and identification methods of miRNAs and their target genes in plants， as well as the advantages and disadvantages of different methods in recent years，so to provide the reference and instruction for further studying the function of miRNAs and their target genes.

Keywords miRNA ； target gene ； detection methods

對于小分子RNA的關(guān)注，始于2000年6月人類基因組計劃的完成，該計劃的完成向人類揭示了高等哺乳動物基因組中蛋白質(zhì)編碼序列不足2%，而98%以上的序列都是非蛋白質(zhì)編碼DNA[1]。此后，隨著多種模式植物基因組測序的完成，科學(xué)家發(fā)現(xiàn)，非編碼序列占絕大多數(shù)的情況同樣存在于植物基因組中。起初許多學(xué)者將此類非蛋白編碼基因認(rèn)為是“垃圾DNA”，但隨后研究發(fā)現(xiàn)，生物細(xì)胞所表達(dá)的RNA中，除了前人熟悉的tRNA、rRNA、mRNA外，還蘊(yùn)藏著大量不編碼蛋白質(zhì)的非編碼RNA分子（noncoding RNA，ncRNA），這些在生物體基因組上廣泛分布的非蛋白質(zhì)編碼區(qū)并非“垃圾DNA”，而是包含了大量轉(zhuǎn)錄調(diào)控元件的非編碼RNA基因[2]，它們能夠調(diào)控真核細(xì)胞發(fā)育的多個方面，如影響功能基因的表達(dá)、調(diào)控細(xì)胞周期和個體發(fā)育等多種行為[3]。ncRNA的發(fā)現(xiàn)，打破了長期以來對RNA僅限于DNA的“轉(zhuǎn)錄信使”這一角色的認(rèn)識，使科學(xué)家重新審視對細(xì)胞功能及其進(jìn)化的認(rèn)識。非編碼RNA分子在多個物種中不斷被發(fā)現(xiàn)和鑒定，揭示出了生物體內(nèi)存在著一個巨大的且不為人知的“非編碼RNA世界”。同時，基于不同物種體內(nèi)非編碼RNA分子的研究也拓生出了一門新的科學(xué)——RNA組學(xué)（RNomics）。目前，RNA組學(xué)研究已成為后基因組學(xué)時代的一個重要科學(xué)前沿[4]，尤其是微RNA（microRNA，miRNA）的鑒定和功能研究，已成為繼小干擾RNA（short interfering RNA，siRNA）之后新的研究熱點(diǎn)之一。本研究對miRNA的驗(yàn)證方法、miRNA靶基因的鑒定方法進(jìn)行了概述，為進(jìn)一步開展植物miRNA的功能研究提供借鑒和參考。

1 miRNA研究進(jìn)展

1.1 miRNA的發(fā)現(xiàn)

miRNA于1993年首次在秀麗新小桿線蟲（caenorkabditi solegans）中被發(fā)現(xiàn)，當(dāng)時研究人員在對線蟲突變體進(jìn)行遺傳分析時，發(fā)現(xiàn)了一類長度為22個核苷酸的小分子RNA，它們能夠調(diào)控細(xì)胞發(fā)育的時序，并將其命名為lin-14 RNA[5]，進(jìn)一步研究發(fā)現(xiàn)lin-14 RNA雖不編碼蛋白質(zhì)，卻能與靶基因 lin-14的3′非翻譯區(qū)（3′untranslated region，3′UTR）以堿基互補(bǔ)配對的方式結(jié)合，并抑制Lin-14基因蛋白的翻譯，從而調(diào)控線蟲的細(xì)胞發(fā)育由第一期向第二期轉(zhuǎn)化。然而，當(dāng)Lin-14 RNA首次被發(fā)現(xiàn)時，并未引起科研人員太多的關(guān)注，直到2000年另一個類似的具有細(xì)胞時序調(diào)節(jié)功能的小分子RNA——let-7 RNA被發(fā)現(xiàn)后，這些非編碼RNA才逐步得到人們的關(guān)注[6]。由于lin-14 RNA和let-7 RNA具有顯著的細(xì)胞時序調(diào)控功能，起初被命名為小時序RNA（small temporal RNA，stRNA）。隨著研究的深入，在人和小鼠等哺乳動物中同樣發(fā)現(xiàn)了大量具有轉(zhuǎn)錄后調(diào)節(jié)功能的小分子RNA。2001年世界不同國家的科研人員將這一類具有轉(zhuǎn)錄后調(diào)節(jié)功能的小分子RN再統(tǒng)一命名為miRNA[7-9]。

植物中miRNA于2002年首次在擬南芥小分子RNA文庫中被發(fā)現(xiàn)。盡管植物中小分子RNA發(fā)現(xiàn)的時間比動物晚了將近10 a，但植物中小分子RNA的鑒定方法要比動物中完善。目前，不僅在擬南芥和水稻等模式植物中發(fā)現(xiàn)和驗(yàn)證了大量的miRNA，而且在其它植物中也發(fā)現(xiàn)和驗(yàn)證了數(shù)量可觀的miRNA。隨著生物體內(nèi)發(fā)現(xiàn)的miRNA家族成員的不斷增多，為了更好地對不同物種中發(fā)現(xiàn)的miRNA進(jìn)行登記和命名，國際上建立了miRNA文庫—miRBase[10-11]。截至2015年5月，在miRNA數(shù)據(jù)庫 （http：microrna.sanger.ac.uk/sequenees/）中注冊提交的miRNA成熟序列已達(dá)28 654條，極大地豐富了人們對生物體內(nèi)非編碼小分子RNA的認(rèn)識。

1.2 植物miRNA的生成過程

研究發(fā)現(xiàn)大部分植物miRNA在基因組上具有一個獨(dú)立的轉(zhuǎn)錄單位，只有少數(shù)miRNA在基因組上呈串聯(lián)重復(fù)存在，在轉(zhuǎn)錄過程中由一個轉(zhuǎn)錄單位可加工形成多個miRNA。miRNA在植物體內(nèi)的生成過程是一個復(fù)雜的、涉及多種酶和功能蛋白參與的過程[12]，具體步驟如圖1所示：首先植物基因組上具有特定位點(diǎn)的MIR基因，在RNA聚合酶的作用下轉(zhuǎn)錄成一個或多個較長的初級轉(zhuǎn)錄產(chǎn)物（pri-miRNA），然后初級轉(zhuǎn)錄物在Drosha酶的作用下[13]，進(jìn)一步加工形成一個具有特定二級莖環(huán)（stem-loop）結(jié)構(gòu)的次級轉(zhuǎn)錄物（pre-miRNA），這一系列過程是在細(xì)胞核內(nèi)完成的；接著pre-miRNA要在nuclear export的作用下從細(xì)胞核運(yùn)輸?shù)郊?xì)胞質(zhì)中，由Dicer酶通過兩步法進(jìn)行酶切加工，隨后酶切產(chǎn)物在HYL1（HYPONASTICLEAVES1）、HEN1（HUAENHANCER1）、HST（HASTY）等蛋白的作用下，形成成熟的miRNA[14-15]。此外，成熟的單鏈miRNA必須與Argonaute（AGO）蛋白相互作用形成一個RNA誘導(dǎo)沉默復(fù)合物（RNA-induced silencing complex，RISC），才能與目標(biāo)mRNA結(jié)合而發(fā)揮作用[16-17]。

1.3 miRNA的特點(diǎn)

miRNA是生物體內(nèi)一類5′端帶磷酸基團(tuán)，3′端帶羥基基團(tuán)的內(nèi)源性小分子RNA。在生物體內(nèi)miRNA是由非蛋白編碼基因轉(zhuǎn)錄物形成的莖環(huán)結(jié)構(gòu)加工而來，并以堿基互補(bǔ)配對的方式靶向mRNAs，進(jìn)而影響目標(biāo)mRNAs的翻譯和穩(wěn)定性。從目前的研究結(jié)果來看，miRNA主要有以下特征：（1）miRNA廣泛存在于生物體內(nèi)。不僅在模式植物擬南芥[18]和水稻[19]中發(fā)現(xiàn)了大量的miRNA，而且在一些重要農(nóng)作物如玉米[20]、小麥[21]、棉花[22]、及大豆[23]中發(fā)現(xiàn)了大量的miRNA，甚至在苔蘚、蕨類等低等植物中也發(fā)現(xiàn)了大量的miRNA存在[24]。（2）miRNA沒有開放閱讀框（ORF）及蛋白質(zhì)編碼基因。在動植物中絕大多數(shù)miRNA產(chǎn)生于基因間隔和非蛋白編碼區(qū)，只有極少數(shù)miRNA產(chǎn)生于基因的啟動子區(qū)，它們是由不同于mRNA的轉(zhuǎn)錄單元表達(dá)產(chǎn)生的[25]。（3）miRNA具有典型的長度范圍，長度約為20～24 nt。大多數(shù)植物miRNA的長度都集中在21～23 nt，但miRNA的前體長度變化較大，有時長度甚至超過1 kb[26-27]。（4）不同物種中miRNA基因位點(diǎn)通常成簇存在。如研究發(fā)現(xiàn)，一些miRNA在生成過程中通常使用一個共同的前體結(jié)構(gòu)，即來源于物種基因組上的同一個miRNA基因簇，并由該前體結(jié)構(gòu)加工形成多個不同的成熟miRNA[28]。（5）miRNA具有典型的序列保守性。與蛋白編碼基因相比，miRNA在進(jìn)化上更具保守性。如在擬南芥中發(fā)現(xiàn)的42個miRNA家族也廣泛存在于水稻、玉米和高粱等其它作物中[29]。同時，miRNA的序列保守性，也為從其它物種中鑒定并克隆miRNA提供了方便。（6）miRNA的表達(dá)具有時序性和組織特異性的特征[30-31]。miRNA在不同組織或器官中的表達(dá)量差異很大，在同一器官的不同發(fā)育時期表達(dá)量也有很大差異，表明在生物體的不同發(fā)育階段，需要開啟和關(guān)閉不同的miRNA來調(diào)控其發(fā)育進(jìn)程。（7）大多數(shù)miRNA的成熟體都是從前體5′端或3′端的一條臂上加工而來，只有少數(shù)miRNA的成熟體是由其前體的兩條壁上同時加工產(chǎn)生[32]。（8）成熟miRNA的 5′端有一個磷酸基團(tuán)，3′端為羥基，且miRNA的起始堿基多為尿嘧啶核苷酸[33-34]。（9）與動物miRNA的加工和成熟過程需要Drosha、Dicer和Exportin-5等蛋白所不同的是，植物miRNA的生成過程必須經(jīng)過Dicer-likel（DCL1）、Hasty（HST）、Hyponastie leaves l（HYLI），以及Hua enhancerl（HENI） 等蛋白的加工過程[35-37]，這些蛋白在植物miRNA的加工和成熟過程中起著至關(guān)重要的作用，它們可直接影響植物體內(nèi)miRNA的生成、種類、以及分布。（10）大多數(shù)植物miRNA與其靶基因幾乎完全互補(bǔ)，而動物miRNA與其靶基因則不完全配對[38]。此外，植物miRNA與靶基因的結(jié)合位點(diǎn)不僅限于靶基因的3′非翻譯區(qū)（untranslated region，UTR），還可以位于轉(zhuǎn)錄區(qū)域[39-40]。植物中miRNA及靶基因的高度互補(bǔ)性也為相關(guān)miRNA的靶基因鑒定和生物信息學(xué)預(yù)測提供了很大便利。

2 miRNA的鑒定方法

前人研究表明，在哺乳動物中miRNA的異常表達(dá)與癌癥[41]、神經(jīng)紊亂[42]、以及心臟疾病的發(fā)生密切相關(guān)[43-44]，如在大腸癌和乳腺癌細(xì)胞中，miRNA的表達(dá)量較正常組織細(xì)胞有較大改變，人們推測miRNA可能在癌變過程中起到原癌基因和抑癌基因的作用[45]。在植物體內(nèi)miRNA能夠顯著調(diào)控其多個生物學(xué)進(jìn)程，對植物的生長發(fā)育至關(guān)重要。如與miRNA代謝及其作用機(jī)制相關(guān)的幾個重要基因的突變體，均體現(xiàn)出發(fā)育過程中多向性的缺陷[46-47]。因此，對生物體不同組織或樣本中miRNA的表達(dá)特征進(jìn)行檢測和分析，將有助于人們進(jìn)一步了解相關(guān)miRNA與生物體發(fā)育的關(guān)系，并為相關(guān)miRNA的功能研究提供信息和依據(jù)。

由于成熟miRNA具有片段小、細(xì)胞中表達(dá)水平普遍偏低、且同一miRNA家族成員間僅有1-2個堿基的差別等特點(diǎn)，這些不利因素給miRNA的檢測帶來了較大困難和挑戰(zhàn)。盡管如此，前人還是根據(jù)miRNA的表達(dá)特征發(fā)明了以下5種miRNA檢測方法：Northern blotting法[48]、微陣列法[49]、克隆和測序法[50]、以及實(shí)時熒光定量PCR法[51]。根據(jù)miRNA的檢測過程中是否需要擴(kuò)增樣本，又可將以上5種方法大體歸為兩類：一類是不需要樣本擴(kuò)增的探針直接雜交法，如Northern blotting法；另一類是基于樣本miRNA擴(kuò)增的方法，如實(shí)時熒光定量PCR法?？傊瑹o論何種檢測方法，一種理想的miRNA定量檢測方法都應(yīng)具有靈敏度高、特異性強(qiáng)、通量高和成本低等優(yōu)點(diǎn)，才能適應(yīng)生物學(xué)的發(fā)展需求。

2.1 基于探針雜交技術(shù)的Northern blotting法

基于探針雜交技術(shù)的Northern blotting法是目前最為經(jīng)典的一種miRNA檢測方法，該方法可以直接對樣品中的miRNA進(jìn)行檢測，而不需要對樣本中miRNA進(jìn)行預(yù)擴(kuò)增。它的基本原理是先將標(biāo)記過的miRNA探針固定在尼龍膜上，然后再將樣本RNA與經(jīng)過標(biāo)記的探針進(jìn)行雜交，最后洗滌多余的雜交探針后進(jìn)行信號檢測。該方法在miRNA的檢測中最為普遍，其優(yōu)點(diǎn)是可同時對多個不同組織中的miRNA進(jìn)行檢測，且檢測成本較低。由于經(jīng)典的Northern blotting方法是基于固相探針雜交技術(shù)，所以該檢測方法最大的缺點(diǎn)是雜交速度慢、靈敏度低；同時，該方法在檢測過程需要耗費(fèi)大量的RNA樣本，對于那些RNA產(chǎn)率很低和容易降解的樣品不太實(shí)用。

為了克服Northern blotting探針雜交法特異性差和靈敏度低的缺點(diǎn)，人們發(fā)明了用鎖核酸（LNA）修飾的雜交探針來取代傳統(tǒng)探針的方法，即改進(jìn)的Northern blotting方法[52]。由于LNA是一類寡核苷酸衍生物，它對核酸分子具有較高的親和性，可滲入到核酸序列的任何位置，因而可以提高LNA探針與目標(biāo)miRNA分子結(jié)合后的雙鏈熱穩(wěn)定性，使得檢測的靈敏度和特異性顯著提高[53]。

2.2 基于微陣列芯片的高通量檢測方法

雖然傳統(tǒng)的Northern blotting法在miRNA檢測中較為常見，但該方法操作步驟繁瑣，不能同時對樣品中多個不同的miRNA進(jìn)行檢測分析，其miRNA的檢測通量明顯不夠。為了克服這一缺點(diǎn)，人們又發(fā)明了miRNA微陣列芯片技術(shù)（Microarray），它可以在同一塊芯片上固定多個與miRNA互補(bǔ)的探針，通過對樣本中相關(guān)miRNA進(jìn)行雜交后，最后進(jìn)行相關(guān)miRNA的信號檢測[54]。盡管該方法實(shí)現(xiàn)了miRNA的高通量檢測分析，但該方法只能檢測小RNA數(shù)據(jù)庫中已經(jīng)提交的miRNA序列，不能對樣品中新的未知miRNA進(jìn)行檢測，導(dǎo)致樣品中一些miRNA表達(dá)信息的缺失，這是該方法的顯著缺陷。

2.3 克隆測序法

基于探針技術(shù)的Northern blotting法和微陣列芯片法都是對已知序列的miRNA進(jìn)行定量檢測。而克隆測序法則可以對樣品中各種不同類型的miRNA進(jìn)行檢測分析?？寺y序法的基本原理是：首先富集樣品中的小分子RNA，然后通過反轉(zhuǎn)錄構(gòu)建樣品相關(guān)的小分子RNA的cDNA文庫，接著對cDNA文庫進(jìn)行PCR擴(kuò)增，最后將PCR擴(kuò)增產(chǎn)物連接到表達(dá)載體上進(jìn)行克隆測序，并根據(jù)文庫中序列出現(xiàn)的豐度對相關(guān)miRNA進(jìn)行定量分析。由于克隆法在測序前必須對樣本進(jìn)行細(xì)菌克隆，然后再對每一個克隆進(jìn)行測序，因而需要消耗較大的初始樣本量（幾百微克RNA），并投入大量的人力和財力[55]。此外，基于克隆測序的方法通常是Sanger測序法，不僅測序效率低，而且費(fèi)時費(fèi)工，這些都給深入研究miRNA的功能帶來了諸多不便。

2.4 實(shí)時定量PCR法

前面幾種方法在檢測過程中都不需要對樣本進(jìn)行擴(kuò)增，而實(shí)時定量PCR法則需要擴(kuò)增樣本。目前采用實(shí)時定量PCR對miRNA進(jìn)行定量檢測主要有以下3種方法：第一種方法稱之為引物延伸定量PCR法（PE-qPCR），它的基本原理是在常規(guī)引物延伸技術(shù)（PE）的基礎(chǔ)上，用一個加尾的特異性引物將miRNA反轉(zhuǎn)錄成cDNA，然后利用一個鎖核酸（LNA）修飾的miRNA特異性反向引物和一個加尾的通用引物序列，對相關(guān)miRNA進(jìn)行PCR擴(kuò)增[56]；第二種方法稱之加尾法，它的基本原理是先用poly（A）聚合酶對樣品總RNA進(jìn)行處理，使得樣品中miRNA的3′端加上poly（A）尾巴，然后對其進(jìn)行反轉(zhuǎn)錄，使得cDNA序列加上poly（T）接頭，最后再用通用引物和miRNA序列特異引物序列，進(jìn)行miRNA的PCR擴(kuò)增[57]。由于該方法在PCR擴(kuò)增前，需要對miRNA進(jìn)行加尾反應(yīng)，因此該方法被稱為加尾法；第三種方法稱之為stem-loop RT-qPCR檢測法，其原理是首先利用一條特殊的莖-環(huán)狀引物將miRNA反轉(zhuǎn)錄成cDNA雙鏈，再利用與miRNA互補(bǔ)的特異性正反向引物和探針，對目標(biāo)miRNA進(jìn)行定量檢測。

利用以上3種方法對miRNA進(jìn)行PCR擴(kuò)增時，主要采用2種信號檢測模式，即熒光染料法和探針法。利用熒光染料法對miRNA進(jìn)行PCR擴(kuò)增時，由于不能進(jìn)行多重PCR反應(yīng)，同時引物二聚體、單鏈二級結(jié)構(gòu)以及錯誤擴(kuò)增產(chǎn)物的出現(xiàn)會增加熒光值，進(jìn)而影響檢測結(jié)果的準(zhǔn)確性，限制了其在miRNA定量檢測中的應(yīng)用；而特異性熒光定量PCR中TaqMan熒光探針法則沒有上述缺點(diǎn)，在實(shí)際操作中被廣泛采納[58-59]。目前，從不同物種中miRNA的檢測效果來看，由于熒光探針法在miRNA的PCR擴(kuò)增中不僅表現(xiàn)出特異性強(qiáng)，靈敏度高，樣品消耗少，定量線性范圍寬等特點(diǎn)[60-62]，能準(zhǔn)確區(qū)分某一miRNA家族中高度同源的不同miRNA序列，甚至還能檢測出只有一個堿基差異的不同miRNA表達(dá)水平，因而在miRNA的定量檢測中得以廣泛應(yīng)用。然而，盡管探針法在miRNA的定量檢測中表現(xiàn)出較大優(yōu)勢，但在實(shí)際應(yīng)用過程中每一條miRNA序列都需要設(shè)計相應(yīng)的探針，費(fèi)用相對較高，限制了該技術(shù)的大規(guī)模應(yīng)用。

2.5 高通量測序技術(shù)

植物中蘊(yùn)藏著巨大數(shù)量的小分子RNA世界，而傳統(tǒng)的Sanger測序法在測序中程序復(fù)雜，效率較低，且測序深度有限，在很大程度上制約著植物中miRNA的發(fā)現(xiàn)速度。為了克服傳統(tǒng)Sanger測序法的缺點(diǎn)，更好的挖掘物種中未知的miRNA，自2005年科研人員發(fā)明了以大規(guī)模平行測序（Massively Parallel Signature Sequencing，MPSS）為代表的第二代高通量測序技術(shù)，來挖掘植物中蘊(yùn)藏的miRNA種類[63-65]。盡管第二代高通量測序技術(shù)獲得的序列較短，但它們在測序過程中可以對上萬個樣本同時進(jìn)行測序，極大地提高了測序效率，具有速度快、成本低、覆蓋度深、產(chǎn)出巨大等優(yōu)點(diǎn)[66-67]。因此，非常適合小分子RNA測序，尤其是對那些拷貝數(shù)低的小分子RNA。

尤其值得一提的是，二代測序技術(shù)中Solexa測序技術(shù)較其它二代測序技術(shù)表現(xiàn)出了前所未有的優(yōu)勢。它能一次檢測上億個核苷酸片段，且測定成本僅為常規(guī)毛細(xì)管電泳測序技術(shù)成本的1%，是二代測序技術(shù)中高質(zhì)量、高通量、低成本的測序法，具有很高的應(yīng)用前景。該測序技術(shù)的基本流程為：先提取高質(zhì)量的樣品總RNA，然后用聚丙烯酰胺凝膠電泳分離大小為18～30 nt之的小RNA片段，接著用T4 RNA連接酶將該范圍內(nèi)的小分子RNA加上5′端接頭和3′端接頭，并對兩端加上接頭的小RNA分子進(jìn)行反轉(zhuǎn)錄和PCR擴(kuò)增。最后將PCR擴(kuò)增后得到的cDNA放置在flow cell上，進(jìn)行上機(jī)測序，從而直接獲得所有miRNA的序列信息。

由于植物中miRNA的表達(dá)具有時空特異性，并隨著諸多生物學(xué)過程的變化而發(fā)生變化，而且在植物的不同發(fā)育時期、不同生長部位其作用機(jī)理表現(xiàn)出較大差異，這就要求miRNA的檢測過程必須具備規(guī)?；?、高通量、假陽性率低等特征，而新一代高通量測序技術(shù)能很好的解決以上問題，較以往的miRNA檢測方法表現(xiàn)出較大優(yōu)勢（表1）。目前，小RNA數(shù)據(jù)庫中提交的絕大多數(shù)miRNA序列都是通過新一代高通量測序技術(shù)發(fā)現(xiàn)和鑒定而來。同時，高通量測序技術(shù)在植物學(xué)研究領(lǐng)域的廣泛應(yīng)用，不僅提升了植物中miRNA的發(fā)掘力度，促進(jìn)了植物中miRNA的功能研究，提高了人們對不同物種中小分子RNA進(jìn)化的認(rèn)識，也為全面、深入研究植物小分子RNA的生物學(xué)功能提供最基本的素材。

3 miRNA靶基因鑒定方法

為了更好地了解miRNA的功能及其下游靶基因的相互作用關(guān)系，對相關(guān)miRNA進(jìn)行靶基因鑒定就顯得十分必要。目前，主要采用以下3種方法來鑒定相關(guān)miRNA的靶基因。

3.1 利用生物信息學(xué)預(yù)測miRNA的靶基因

植物體內(nèi)miRNA通常與靶基因進(jìn)行完全或接近完全的配對來引起mRNA的剪切，從而達(dá)到調(diào)控基因表達(dá)的目的。與動物miRNA相比，植物miRNA與靶基因的互補(bǔ)程度較高，這就為植物miRNA的靶基因預(yù)測帶來很大便利?；趍iRNA與靶基因的互補(bǔ)配對關(guān)系，前人發(fā)明了利用生物信息學(xué)方法來預(yù)測相關(guān)miRNA的靶基因。根據(jù)前人的研究和總結(jié)，miRNA靶基因的生物信息學(xué)預(yù)測過程中必須遵循一定的流程和規(guī)則[68-70]，具體規(guī)則如下：（1）miRNA與靶基因不得超過4個堿基錯配（G-U配對認(rèn)為0.5個錯配）；（2）miRNA與靶基因的匹配過程中不能有超過2處的相鄰位點(diǎn)堿基錯配；（3）miRNA與靶基因復(fù)合體中，從miRNA的5′端起第2～12個位點(diǎn)不得有相鄰位點(diǎn)發(fā)生錯配；（4） miRNA與靶基因的匹配過程中不得發(fā)生第10或11位點(diǎn)堿基錯配。大多數(shù)miRNA都在第10或第11位點(diǎn)對mRNA進(jìn)行切割，因此這兩個位點(diǎn)不允許有堿基錯配；（5）miRNA與靶基因匹配中，從miRNA的5′端起，第1～12個位點(diǎn)的種子序列之間不得有超過2.5個堿基錯配；（6）miRNA與靶基因匹配后的最小結(jié)合自由能（MFE），不能小于該miRNA與靶基因完全互補(bǔ)結(jié)合時最小接合自由能的75%。

早期miRNA的功能研究中，許多miRNA的靶基因都是通過生物信息學(xué)預(yù)測鑒定而來，當(dāng)時為相關(guān)miRNA的功能研究帶來了一定便利。由于傳統(tǒng)的利用生物信息學(xué)預(yù)測miRNA靶基因的方法主要依靠計算機(jī)模擬進(jìn)行驗(yàn)證，該方法存在一定的局限性，如相關(guān)miRNA的靶基因鑒定結(jié)果會出現(xiàn)數(shù)據(jù)海量，假陽性高，準(zhǔn)確性差等弊端。因此，該方法只適合對物種間保守的miRNA進(jìn)行靶基因鑒定，而對于那些未知的新miRNA進(jìn)行靶基因鑒定時效果很差，這些不利因素限制了miRNA的后續(xù)研究進(jìn)程。

3.2 利用5′ RLM-RACE技術(shù)鑒定miRNA的靶基因

為了從實(shí)驗(yàn)水平直接驗(yàn)證miRNA的靶基因，人們隨后又發(fā)明了基于RNA連接反應(yīng)的5′ Race技術(shù)來鑒定相關(guān)miRNA的靶基因[71]。該技術(shù)的原理是將提取的RNA連接上5′ Race接頭，在反轉(zhuǎn)錄酶和GeneRacer Oligo dT 引物的作用下合成cDNA單鏈，然后cDNA單鏈在GeneRacer 5′ Primer（5′-CGACTGGAGCACGAGGACACTGA-3′）和GeneRacer 3′ Primer（5′-GCTGTCAACGATACGCTACGTAACG-3′）的作用下產(chǎn)生非特異性基因的5′ Race擴(kuò)增產(chǎn)物文庫。接著，使用GeneRacer 5′ Nested Primer和Gene-specific primers進(jìn)行相關(guān)特異性基因的5′ Race產(chǎn)物擴(kuò)增。最后，對特異性基因的5′ Race擴(kuò)增產(chǎn)物進(jìn)行膠回收、克隆、測序分析。

該方法雖能準(zhǔn)確鑒定某一特定組織中相關(guān)miRNA的靶基因，但在實(shí)際應(yīng)用過程中，仍存在一些不利因素[72-73]。如該技術(shù)在操作過程中需要消耗大量的樣品RNA，對于那些采集比較困難的樣品中miRNA的靶基因鑒定不太現(xiàn)實(shí)。其次，5′Race試劑盒的價格偏高，致使相關(guān)miRNA靶基因的鑒定成本過高。此外，在實(shí)際操作過程中，每次僅能從樣品中鑒定有限數(shù)目的靶基因，而不能對該組織中的全部miRNA進(jìn)行靶基因鑒定，使得組織樣品中很多重要的靶基因切割位點(diǎn)信息遺失。因此，該技術(shù)只適合小范圍對樣品中的miRNA進(jìn)行靶基因鑒定，而不能對樣品中絕大多數(shù)miRNA的靶基因進(jìn)行分析，這些不利因素限制了該技術(shù)在miRNA功能研究中的大規(guī)模應(yīng)用。

3.3 利用降解組測序技術(shù)鑒定miRNA的靶基因

考慮到生物信息學(xué)預(yù)測miRNA靶基因時存在假陽性偏高、準(zhǔn)確性差等弊端，而5′ Race技術(shù)又存在成本偏高，不能對組織中全部miRNA的靶基因進(jìn)行高通量分析的缺點(diǎn)，這些不利因素困擾著miRNA的研究進(jìn)程。近年來隨著新一代高通量測序（High-throughput Sequencing） 技術(shù)的出現(xiàn)和發(fā)展，出現(xiàn)了一種新的稱為降解組測序的方法（Degradome Sequencing） ，該技術(shù)結(jié)合了高通量測序技術(shù)與生物信息學(xué)分析的優(yōu)勢，能夠很好的解決miRNA靶基因鑒定過程中存在的弊端和缺陷[74-76]，使之在miRNA的研究中得以大規(guī)模應(yīng)用。

降解組測序的原理是基于植物體內(nèi)絕大多數(shù)miRNA是通過剪切作用來實(shí)現(xiàn)對靶基因的調(diào)控，且miRNA對靶基因的剪切位點(diǎn)通常發(fā)生在與mRNA互補(bǔ)區(qū)域的第10或第11位堿基位點(diǎn)上。如果mRNA受到相關(guān)miRNA的剪切作用，mRNA被剪切后就會產(chǎn)生兩個斷片——即含有帽子結(jié)構(gòu)的5′端剪切片段和帶有poly A尾巴的3′端剪切片段。由于mRNA被剪切后產(chǎn)生的3′端剪切片段含有5′單磷酸自由基團(tuán)和3′poly A尾巴，該斷片可被RNA連接酶連接，其連接產(chǎn)物可用于下一步的高通量測序；而含有帽子結(jié)構(gòu)的5′端剪切片段因缺少5′單磷酸基團(tuán)，無法被RNA連接酶連接，因而不能進(jìn)入下一步的測序?qū)嶒?yàn)。緊接著使用HiSeq2000測序儀對連接產(chǎn)物進(jìn)行測序，得到50 nt長度的原始數(shù)據(jù)。原始數(shù)據(jù)再經(jīng)去接頭、去污染、去低質(zhì)量等步驟后獲得“干凈”的序列——即稱為降解序列。降解序列經(jīng)過一系列的注釋之后，就可得到mRNA的降解片段。最后，獲得的mRNA降解片段通過與miRNA的比對分析后推測出其互補(bǔ)序列，并借助于生物信息學(xué)工具target finder做一個反向預(yù)測來尋找mRNA-miRNA的配對關(guān)系。綜合這兩方面的結(jié)果，可以直觀地發(fā)現(xiàn)在mRNA序列的某個位點(diǎn)會出現(xiàn)一個波峰，而該峰值對應(yīng)的mRNA位點(diǎn)正是候選miRNA的剪切位點(diǎn)。

利用降解組測序鑒定miRNA的靶基因具有以下特點(diǎn)：（1）高通量性：一次測序可以檢測到1 000萬條以上的序列信息；（2）高準(zhǔn)確性：從幾個到幾十萬個拷貝的精確計數(shù)；（3）高分辨率：可以精確檢測到單堿基的差異；（4）重復(fù)性好：深度測序保證了檢測的隨機(jī)性，不需要技術(shù)重復(fù)等特點(diǎn)。目前，該方法已被成功應(yīng)用于水稻、棉花、玉米等重要農(nóng)作物miRNA靶基因的篩選和鑒定[77-80]。同時，該方法的問世也為miRNA靶基因的鑒定和功能研究帶來了革命性的變革，有利地推動了動植物miRNA的研究進(jìn)程。

降解組測序技術(shù)的出現(xiàn)突破了傳統(tǒng)miRNA靶基因鑒定技術(shù)的局限性，為科研人員深入開展miRNA的功能研究提供了更加靈活和有效的分析方法，該技術(shù)檢測到的靶基因數(shù)量和準(zhǔn)確率相對于生物信息學(xué)預(yù)測有較大提高，而且能大大縮小后續(xù)實(shí)驗(yàn)驗(yàn)證的范圍，其實(shí)際效果是傳統(tǒng)的研究方法所無法比擬的（表2）。在實(shí)際操作中，借助于生物學(xué)實(shí)驗(yàn)手段對樣本中miRNA的靶基因進(jìn)行高通量測序，可以快速確認(rèn)上游miRNA的剪切位點(diǎn)及其相互作用的靶基因種類，進(jìn)而得知樣本中miRNA與下游靶基因的調(diào)控通徑。同時，結(jié)合不同手段的生物信息學(xué)分析，如基因本體論分析（Gene Ontology，GO）、代謝通路（Pathway）分析與表型性狀分析，可以從不同角度對植物miRNA的調(diào)控機(jī)制進(jìn)行整體分析。目前，降解組測序技術(shù)已被廣泛應(yīng)用于植物生長發(fā)育調(diào)控以及逆境脅迫相關(guān)的生物學(xué)響應(yīng)機(jī)制研究中，成為研究miRNA及其靶基因相互作用關(guān)系，進(jìn)一步探究植物生長發(fā)育調(diào)控和脅迫相關(guān)的調(diào)節(jié)通路、以及從基因表達(dá)的轉(zhuǎn)錄水平和轉(zhuǎn)錄后水平構(gòu)建調(diào)控網(wǎng)絡(luò)的一種全新而有效的方法。

4 展望

miRNAs的發(fā)現(xiàn)是生命科學(xué)領(lǐng)域的一個重大突破，它拓展和豐富了人們對生物體內(nèi)基因表達(dá)和小分子RNA調(diào)控功能的認(rèn)識，促使生物學(xué)家重新反思對細(xì)胞進(jìn)化的理解和認(rèn)知。自1993年動物中發(fā)現(xiàn)第一個miRNA lin4以來，miRNA領(lǐng)域的相關(guān)研究得以飛速發(fā)展，并逐步成為生物學(xué)研究領(lǐng)域追蹤的熱點(diǎn)之一。現(xiàn)已證實(shí)生物體細(xì)胞中存在著數(shù)量龐大的miRNAs，它們是基因表達(dá)調(diào)控網(wǎng)絡(luò)中豐富而重要的組成部分。在動物研究領(lǐng)域，諸多研究表明，miRNA與人類一些重大疾病的產(chǎn)生密切相關(guān)，理解miRNA介導(dǎo)的基因沉默機(jī)制，可對諸如癌癥和心腦血管疾病產(chǎn)生的根本原因提供新的見解與治療方案；在植物研究領(lǐng)域，miRNA與一些重要農(nóng)作物產(chǎn)量形成的分子機(jī)制密切相關(guān)，通過對一些物種中重要農(nóng)藝性狀相關(guān)miRNA及其靶基因進(jìn)行鑒定和功能分析，對于詮釋一些重要器官中特定miRNA的生物學(xué)功能與代謝機(jī)制，破解重要農(nóng)作物產(chǎn)量形成中的分子調(diào)控機(jī)制具有重要意義。隨著miRNA及其靶基因鑒定技術(shù)的進(jìn)一步完善和深入，以及大規(guī)模高通量測序技術(shù)與平臺的廣泛應(yīng)用，必將促使人們開始對各個物種中特有的miRNA、不同發(fā)育狀態(tài)下的miRNA、以及各種脅迫條件的miRNA進(jìn)行挖掘和鑒定，以此來全面探究該物種中miRNA在不同條件及不同發(fā)育時期的功能。同時，隨著miRNA及靶基因研究手段的進(jìn)一步完善，借助于轉(zhuǎn)基因技術(shù)過量表達(dá)特定的miRNA或獲得靶基因缺失突變體對miRNA的功能進(jìn)行驗(yàn)證技術(shù)的進(jìn)一步完善，這些研究必將推動整個生物學(xué)研究領(lǐng)域的飛速發(fā)展。

縱觀近年來miRNA及其靶基因檢測與鑒定技術(shù)的進(jìn)展，雖然鑒定并提交至小RAN數(shù)據(jù)庫中的miRNA數(shù)量呈指數(shù)增加，但許多miRNA的功能仍然未知，一些提交至小RNA數(shù)據(jù)庫中miRNA的真實(shí)性仍需進(jìn)一步確認(rèn)。在miRNA的鑒定與功能研究領(lǐng)域，仍存在許多未解的困惑和未知的領(lǐng)域?？偟膩碚f，問題主要包括以下6個方面：其一、關(guān)于物種中miRNA種類上限的問題。迄今人們不清楚任何一種物種中miRNA的種類究竟有多少，某一特定物種中發(fā)現(xiàn)的miRNA種類距離該物種中全部miRNA種類的上限還有多少？ 其二、關(guān)于miRNA的靶基因問題。研究發(fā)現(xiàn)很多已經(jīng)鑒定的miRNA根本找不到其作用的靶基因，或者它的靶基因功能未知。其三、關(guān)于miRNA與靶基因的接合效率和作用方式問題。在miRNA與靶基因的匹配過程中，究竟什么因素影響了miRNA與UTR位點(diǎn)的接近性和效率？ 在什么條件下，miRNA對靶基因進(jìn)行裂解或者抑制翻譯？ 其四、關(guān)于miRNA與靶基因的接合位點(diǎn)問題。研究發(fā)現(xiàn)很多情況下，多個不同的miRNAs會作用于同一靶基因，這些miRNA結(jié)合在同一靶基因位點(diǎn)上，它們?nèi)绾纹鹱饔茫?其五、關(guān)于miRNA靶基因功能冗余和疊加的問題。研究發(fā)現(xiàn)，許多miRNA具有相似的功能，它們的靶基因功能也相似。在進(jìn)化中這些miRNA靶基因功能冗余和疊加的本質(zhì)和意義是什么？ 其六，關(guān)于miRNA與其它幾種內(nèi)源性小分子RNA作用機(jī)制的區(qū)別與相互聯(lián)系問題。如生物體內(nèi)miRNAs與siRNAs都可以導(dǎo)致基因表達(dá)沉默的產(chǎn)生，但究竟在多大程度上，miRNAs生成通路與siRNAs生產(chǎn)通路使用相似的同源進(jìn)化基因，何種因素介導(dǎo)了miRNA復(fù)合體由細(xì)胞核向細(xì)胞質(zhì)轉(zhuǎn)運(yùn)，這些問題都有待進(jìn)一步研究和探討。

參考文獻(xiàn)

[1] International Human Genome Sequencing Consortium. Finishing the euchromatic sequence of the human genome[J]. Nature， 2004， 31（7011）： 931-945.

[2] Reinhart B J， Bartel D P. Small RNAs correspond to centromere heterochromatic repeats[J]. Science， 2002， 297（5 588）： 1 831-1 831.

[3] Malone C D， Hannon G J. Small RNAs as Guardians of the Genome[J]. Cell， 2009， 136： 656-668.

[4] 屈良鵠. RNA 組學(xué)： 后基因組時代的科學(xué)前沿[J].中國科學(xué)C 輯，2009，39：1-2.

[5] Lee R C， Feinbaum R L， Ambros V， et al. Elegtms heterochronic gene lin-4 encodes small RNAs with antisense complementarity to Lin-l4[J]. Cell， 1993， 75（5）： 843-854.

[6] Reinhart B J， Bartel D P. Small RNAs correspond to centromere heterochromatic repeats[J]. Science， 2002， 297（5 588）： 1 831-1 831.

[7] Lagos-Quintana M， Rauhut R， Yalcin A， et al. Identification of tissue-specific microRNAs from mouse[J]. Curr Biol， 2002， 12（9）： 735-739.

[8] Pasquinelli A E， Reinhart B J， Slack F， et al. Conservation of the sequence and temporal expression of let-7 heterochronic regulatory RNA[J]. Nature， 2000， 408（6 808）： 86-89.

[9] Lee R C， Ambros V. An extensive class of small RNAs in Caenorhabditis elegans[J]. Science， 2001， 294（5 543）： 862-864.

[10] Griffiths-Jones S， Grocock RJ， Van Dongen S， et al. miRBase： microRNA sequences， targets and gene nomenclature[J]. Nucleic acids research， 2006， 34（1）： 140-144.

[11] Griffiths-Jones S， Saini H K， van Dongen S， et al. miRBase： tools for microRNA genomics[J].Nucleic acids research， 2008， 36（1）： 154-158.

[12] Wierzbicki AT， Haag JR， Pikaard CS. Noncoding transcription by RNA polymerase Pol IVb/Pol V mediates transcriptional silencing of overlapping and adjacent genes[J]. Cell， 2008， 135： 635–648.

[13] Rehwinkel J， Natalin P， Stark A， et al. Genome-wide analysis of mRNAs regulated by Drosha and Argonaute proteins in Drosophila melanogaster[J]. Mol Cell Biol， 2006， 26： 2 965-2 975.

[14] Bernstein E， Caudy A A， Hammond S M， Hannon G J. Role of a bidentate ribonuclease in the initiation step of RNA interference[J]. Nature， 2001， 409（6 818）： 363-366.

[15] Jones-Rhoades M W， Bartel D P， Bartel B. MicroRNAs and their regulatory roles in plants[J]. Annu Rev Plant Biol， 2006， 57： 19-53.

[16] Hammond S C， Bernstein E， Beach D， et al. An RNA-directed nuclease mediates Post transcriptional gene silencing in Drosophila cells[J]. Nature， 2000， 404（6 775）： 293-296

[17] Hutvagner G， Zamore P D. RNAi： nature abhors a double strand[J]. Curr Opin Genet Dev， 2002， 12（2）： 225-232.

[18] Sunkar R， Zhu J K. Novel and stress-regulated micro-RNAs and other small RNAs from Arabidopsis[J]. Plant Cell， 2004， 16（8）： 2 001-2 019.

[19] Sunkar R， Girke T， Jain PK， et al. Cloning and characterization of microRNAs from rice[J]. Plant Cell， 2005， 17（5）： 1 397-1 411.

[20] Mica E， Gianfranceschi L， Pe M E. Characterization of five microRNA families in maize[J]. J Exp Bot， 2006， 57（11）： 2 601-2 612.

[21] Yao Y Y， Guo G G， Ni Z Y， et al. Cloning and characterization of microRNAs from wheat （Triticum aestivum L.） [J]. Genome Biol， 2007， 8（6）： 1-13.

[22] Kwak P B， Wang Q Q， Chen X S， et al. Enrichment of a set of microRNAs during the cotton fiber development[J] .BMC Genomics， 2009， 10（1）： 1-11.

[23] Zhang B， Pan X， Stellwag E J. Identification of soybean microRNAs and their targets[J]. Planta， 2008， 229（1）： 161-182.

[24] Arazi T， Talmor NM， Stav R， et al. Cloning and characterization of micro-RNAs from moss[J].Plant J， 2005， 43（6）： 837-848.

[25] Tang G， Reinhart B J， Bartel D P， et al. A biochemical frame work for RNA silencing in plants[J]. Genes & development， 2003， 17（1）： 49-63.

[26] Allen E， Xie Z， Gustafson AM， et al. microRNA-directed phasing during trans-acting siRNA biogenesis in plants[J]. Cell， 2005， 121（2）： 207-221.

[27] Jung J H， Seo P J， Park C M. MicroRNA biogenesis and function in higher plants[J]. Plant Bio technol Rep， 2009， 3（2）： 111-126.

[28] Sunkar R， Zhu J K. Novel and stress-regulated micro-RNAs and other small RNAs from Arabidopsis[J]. Plant Cell， 2004， 16（8）： 2 001-2 019.

[29] Axtell M J， Bartel D P. Antiquity of microRNAs and their targets in land plants[J]. Plant Cell， 2005， 17（6）： 1 658-1 673.

[30] Bartel D P. MicroRNAs： genomics， biogenesis， mechanism， and function[J].Cell， 2004， 116（2）： 281-297.

[31] Wienholds E， Kloosterman W P， Miska E， et al. MicroRNA expression in zebrafish embryonic development[J]. Science， 2005， 309（5 732）： 310-311.

[32] Ambros V. The functions of animal microRNAs[J]. Nature， 2004， 431（7 006）： 350-355.

[33] Lau N C， Lim L P， Weinstein E G， et al. An abundant class of tiny RNAs with probable regulatory roles in Caenorhabditis elegans[J]. Science， 2001， 294（5 543）： 858-862.

[34] Ambros V， Bartel B， Bartel DP， et al. A uniform system for microRNA annotation[J].RNA， 2003， 9（3）： 277-279.

[35] Papp I， Mette MF， Aufsatz W， et al. Evidence for nuclear processing of plant micro RNA and short interfering RNA precursors[J]. Plant physi， 2003， 132（3）： 1 382-1 390.

[36] Bartel B， Bartel DP. MicroRNAs： at the root of plant development[J]. Plant Physiology， 2003， 132（2）： 709-717.

[37] Lund E， Güttinger S， Calado A， et al. Nuclear export of microRNA precursors[J].Science， 2004， 303（5 654）： 95-98.

[38] Liave C， Xie Z， Kasschau KD， et al. Cleavage of Scarecrow-like mRNA targets directed by a class of Arabidopsis miRNA[J].Science， 2002， 297： 2 053-2 056.

[39] Floyd S K， Bowman J L. Gene regulation： ancient microRNA target sequences in plants[J]. Nature， 2004， 428（6 982）： 485-486.

[40] Wu S， Huang S， Ding J， et al. Multiple microRNAs modulate p21Cip1/Waf1 expression by directly targeting its 3′ untranslated region[J]. Oncogene， 2010， 29（15）： 2 302-2 308.

[41] Cho W C S. OncomiRs： the discovery and progress of microRNAs in cancers[J]. Mol Cancer， 2007， 6（1）： 60.

[42] Schetter A J， Leung S Y， Sohn J J， et al. MicroRNA expression profiles associated with prognosis and therapeutic outcome in colon adenocarcinoma[J]. Jama， 2008， 299（4）： 425-436.

[43] Zimmerman AL， Wu S. MicroRNAs， cancer and cancer stem cells[J]. Cancer letters， 2011， 300（1）： 10-19.

[44] Kong Y W， Ferland-McCollough D， Jackson TJ， et al. microRNAs in cancer management[J].The lancet oncology， 2012， 13（6）： 249-258.

[45] Di Leva G， Garofalo M， Croce CM. MicroRNAs in cancer[J]. Annual review of pathology， 2014， 60： 167-179.

[46] Lu C， Fedoroff N. A mutation in the Arabidopsis HYL1 gene encoding a dsRNA binding protein affects responses to abscisic acid， auxin， and cytokinin[J]. Plant Cell， 2000， 12（12）： 2 351-2 365.

[47] Schauer S E， Jacobsen S E， Meinke D W， et al. DICER-LIKE1： blind men and elephants in Arabidopsis development[J].Trends in plant science， 2002， 7（11）： 487-491.

[48] Pall G S， Codony-Servat C， Byrne J， et al. Carbodiimide-mediated cross-linking of RNA to nylon membranes improves the detection of siRNA， miRNA and piRNA by northern blot[J]. Nucleic Acids Res， 2007， 35（8）： 1-9.

[49] Liu C G， Calin G A， Meloon B， et al. An oligonucleotide microchip for genome-wide microRNA profiling in human and mouse tissues[J]. Proc Natl Acad Sci USA， 2004， 101（26）： 9 740-9 744.

[50] Chen C， Ridzon D A， Broomer AJ， et al. Real-time quantification of microRNAs by stem–loop RT-PCR[J]. Nucleic acids research， 2005， 33（20）： 179-187.

[51] Shi R， Chiang V L. Facile means for quantifying microRNA expression by real-time PCR[J]. Biotechniques， 2005， 39（4）： 519-525

[52] Varallyay E， Burgyan J， Havelda Z. MicroRNA detection by northern blotting using locked nucleic acid probes[J]. Nat Protoc， 2008， 3（2）： 190-196.

[53] Kubota K， Ohashi A， Imachi H， et al. Improved in situ hybridization efficiency with locked-nucleic-acid-incorporated DNA probes[J]. Appl Environ Microbiol， 2006， 72（8）： 5 311-5 317

[54] Castoldi M， Schmidt S， Benes V， et al. miChip： an array-based method for microRNA expression profiling using locked nucleic acid capture probes[J]. Nat Protoc， 2008， 3（2）： 321-329.

[55] Lagos-Quintana M， Rauhut R， Meyer J， et al. New microRNAs from mouse and human[J]. Rna， 2003， 9（2）： 175-179.

[56] Allawi H T， Dahlberg J E， Olson S， et al. Quantitation of microRNAs using a modified Invader assay[J]. Rna， 2004， 10（7）： 1 153-1 161.

[57] Raymond C K， Roberts B S， Garrett-Engele P， et al. Simple， quantitative primer-extension PCR assay for direct monitoring of microRNAs and short-interfering RNAs[J].RNA， 2005， 11（11）： 1 737-1 744.

[58] Chen C， Ridzon D A， Broomer A J， et al. Real-time quantification of microRNAs by stem-loop RT-PCR[J]. Nucleic acids research， 2005， 33（20）： 179-179.

[59] Varkonyi-Gasic E， Wu R， Wood M， et al. Protocol： a highly sensitive RT-PCR method for detection and quantification of microRNAs[J]. Plant methods， 2007， 3（1）： 1-12.

[60] Kramer M F. Stem-Loop RT-qPCR for miRNAs[J]. Current protocols in molecular biology， 2011， 15（10）： 1-15.

[61] Hurley J， Roberts D， Bond A， et al. Next-generation microRNA expression profiling technology[J]. Humana Press， 2012： 33-52.

[62] Salone V， Rederstorff M. Stem-Loop RT-PCR Based Quantification of Small Non-Coding RNAs[J]. Small Non-Coding RNAs： Methods and Protocols， 2015， 8（10）： 103-108.

[63] Brennecke J， Aravin A， Stark A， et al. Discrete Small RNA-generating loci as master regulators of transposon activity in Drosophila[J]. Cell， 2007， 128（6）： 1 089-1 103.

[64] Margulies M， Egholm M， Altman WE， et al. Genome sequencing in microfabricated high-density picolitre reactors[J]. Nature， 2005， 437（7 057）： 376-380.

[65] Cokus S J， Feng S， Zhang X， et al. Shotgun bisulphite sequencing of the Arabidopsis genome reveals DNA methylation patterning[J]. Nature， 2008， 452（7 184）： 215-219.

[66] Sánchez-León N， Arteaga-Vázquez M， Alvarez-Mejía C， et al. Transcriptional analysis of the Arabidopsis ovule by massively parallel signature sequencing[J]. J Exp Bot， 2012， 63（10）： 3 829-3 842.

[67] Nuttle X， Itsara A， Shendure J， et al. Resolving genomic disorder-associated breakpoints within segmental DNA duplications using massively parallel sequencing[J].Nature protocols， 2014， 9（6）： 1 496-1 513.

[68] Qiu CX， Xie FL， Zhu YY， et al. Computational identification of microRNAs and their targets in Gossypium hirsutum expressed sequence tags[J]. Gene， 2007， 395： 49-61.

[69] Rhoades MW， Reinhart BJ， Lim LP， et al. Prediction of plant microRNA targets[J]. Cell， 2002， 110： 513-520.

[70] Meyers BC， Axtell MJ， Bartel B， et al. Criteria for annotation of plant microRNAs[J]. Plant Cell， 2008， 20（12）： 3 186-3 190.

[71] Liave C， Xie Z， Kasschau KD， et al. Cleavage of Scarecrow-like mRNA targets directed by a class of Arabidopsis miRNA.Science， 2002， 297： 2 053-2 056.

[72] Wang C， Han J， Korir NK， et al. Characterization of target mRNAs for grapevine microRNAs with an integrated strategy of modified RLM-RACE， newly developed PPM-RACE and qPCRs[J]. Journal of plant physiology， 2013， 170（10）： 943-957.

[73] Shangguan L， Song C， Han J， et al. Characterization of regulatory mechanism of Poncirus trifoliata microRNAs on their target genes with an integrated strategy of newly developed PPM-RACE and RLM-RACE[J]. Gene， 2014， 535（1）： 42-52.

（下轉(zhuǎn)第36頁）