1、素有“活化石”之稱的銀杏（Ginkgo biloba L.），被認(rèn)為是地球上繼蕨類物種之后出現(xiàn)最早的種子植物之一。銀杏起源歷史悠久，其生殖過(guò)程表現(xiàn)出獨(dú)特的性狀，尤其是雌性生殖器官胚珠的生長(zhǎng)與發(fā)育，更是具有許多獨(dú)特的原始特征。目前大多數(shù)研究集中在有關(guān)銀杏胚珠的微觀形態(tài)結(jié)構(gòu)、胚珠組織形成和生理生化變化等方面，而有關(guān)調(diào)控銀杏胚珠生長(zhǎng)發(fā)育的分子機(jī)制方面鮮有報(bào)道。MicroRNAs（miRNAs），是一類內(nèi)源性的、單鏈的、長(zhǎng)度約為22個(gè)核苷酸(n

2、t)的非編碼小RNA，在真核生物中進(jìn)化具有高度保守性，通過(guò)直接對(duì)靶基因切割或抑制靶基因翻譯的手段，參與到植物生長(zhǎng)發(fā)育、信號(hào)轉(zhuǎn)導(dǎo)過(guò)程和環(huán)境脅迫響應(yīng)等生命進(jìn)程的各個(gè)方面。因此，對(duì)銀杏胚珠開(kāi)展miRNA的鑒定與功能研究具有重要的意義，可為進(jìn)一步研究銀杏胚珠生長(zhǎng)發(fā)育的分子機(jī)理提供有力幫助。
　　本研究中，以銀杏胚珠為實(shí)驗(yàn)材料，采用高通量測(cè)序技術(shù)與生物信息學(xué)分析技術(shù)相結(jié)合的方法，鑒定和分析了銀杏胚珠的miRNA及其靶基因。取得主要結(jié)果如下:

3、
　　(1)通過(guò)測(cè)序，從銀杏胚珠中總共得到12，121，385個(gè)高質(zhì)量的sRNA序列，去除3'接頭缺失序列、5'接頭污染序列、無(wú)插入片段、小于18nt的短序列及Poly A序列等低質(zhì)量序列后，最終得到11，919，364（1，991，395 unique）條干凈序列。
　　(2)對(duì)銀杏胚珠小RNA文庫(kù)中的clean reads進(jìn)行分析與分類注釋。文庫(kù)中的小RNA長(zhǎng)度集中在19-24 nt之間，符合典型miRNA長(zhǎng)度分布。通過(guò)

4、比對(duì)基因組，其中包含核糖體RNA(rRNA)有519，513條序列(31，039 unique)，占到sRNA總數(shù)的4.36％(1.56％ unique);小核RNA（small nuclear RNA，snRNA）5，769條序列(1，014unique)，占到sRNA的0.05％(0.05％ unique);核仁小RNA(small Nucleolus RNA，snoRNA)1，858條序列(376 unique)，占總數(shù)的0.02

5、％(0.02％ unique);轉(zhuǎn)運(yùn)RNA(transfer ribonucleicacid，tRNA)203，471條序列(6，636unique)，占總數(shù)的1.71％;未注釋(unannatation)的sRNA8，627，585條序列(1，932，584unique)，占到總數(shù)的72.38％(97.05％ unique)。
　　(3)將胚珠的sRNA序列比對(duì)到miRBase數(shù)據(jù)庫(kù)（miRBase20.0）中所有植物的miRN

6、A成熟序列以及前體序列，最終比對(duì)得到19，746條唯一同源序列，經(jīng)過(guò)嚴(yán)格篩選之后，最終鑒定得到保守序列144條，隸屬于35個(gè)miRNA家族。
　　(4)將小RNA注釋分析去除各種小RNA序列以后得到的未注釋序列，用來(lái)進(jìn)行新的miRNA預(yù)測(cè)。得到新的miRNA種數(shù)為99種，總表達(dá)量為132，843。在對(duì)新的miRNA進(jìn)行首位堿基偏向性分析中發(fā)現(xiàn)，銀杏胚珠中19-24nt的首位堿基對(duì)堿基U具有偏向性，符合成熟序列首位堿基對(duì)堿基U的偏向

7、性要求。
　　(5)應(yīng)用psRNATarget軟件，對(duì)銀杏胚珠的miRNA進(jìn)行靶基因預(yù)測(cè)，銀杏胚珠中已知的445個(gè)miRNA，預(yù)測(cè)得到2，062個(gè)靶基因，靶基因位點(diǎn)數(shù)為2，487個(gè);123個(gè)新的miRNA，預(yù)測(cè)得到靶基因898個(gè)，共有靶基因位點(diǎn)數(shù)1，331個(gè)。GO富集分析發(fā)現(xiàn)這些靶基因的功能參與銀杏胚珠的代謝過(guò)程、細(xì)胞進(jìn)程等生長(zhǎng)發(fā)育的各個(gè)生命活動(dòng)進(jìn)程。
　　(6)對(duì)銀杏胚珠的靶基因功能注釋分析發(fā)現(xiàn)，大部分靶基因?yàn)橹参镛D(zhuǎn)錄因子

8、，我們從中挑選出具有前體序列，并且與光合作用、代謝合成、抗病以及信號(hào)轉(zhuǎn)導(dǎo)等功能相關(guān)的14個(gè)已知的miRNA和表達(dá)量比較高的11個(gè)新的miRNA，應(yīng)用亞克隆技術(shù)，進(jìn)行前體驗(yàn)證，最終得到10個(gè)miRNA前體與Illumina測(cè)序結(jié)果完全匹配，7個(gè)低于3個(gè)堿基錯(cuò)配。
　　(7)對(duì)銀杏胚珠與銀杏雌株葉片中的miRNA的表達(dá)量進(jìn)行差異分析，發(fā)現(xiàn)已知的miRNA中胚珠相對(duì)于雌株葉片有300個(gè)上調(diào)，171個(gè)下調(diào)，73個(gè)無(wú)顯著差異;新的miRNA

9、中胚珠相對(duì)于雌株葉片有46個(gè)上調(diào)，47個(gè)下調(diào)，21個(gè)無(wú)顯著差異。
　　(8)根據(jù)上面靶基因功能以及前體驗(yàn)證結(jié)果，我們從中一共挑選出19個(gè)miRNA前體序列進(jìn)行熒光定量PCR實(shí)驗(yàn)，以兩個(gè)不同時(shí)期的銀杏胚珠與銀杏雌株葉片為實(shí)驗(yàn)材料，分析miRNA在同一時(shí)期不同組織中以及同一組織不同時(shí)期的表達(dá)差異性。
　　綜上，本研究通過(guò)對(duì)銀杏胚珠中miRNA的鑒定與分析，揭示了銀杏胚珠中存在較多數(shù)量的小RNA，并且這些小RNA參與到了銀杏的生長(zhǎng)