桑樹miRNA的鑒定及其在蠶?；プ髦械墓δ苎芯?pdf

1、研究利用寡食性昆蟲-家蠶和桑葉這對組合，探究家蠶攝食桑葉后，桑樹的miRNA是否進入了家蠶的體內。利用Illumina高通量測序方法對野生桑（川桑）、栽培桑（湖桑32號和倫教109）葉片以及家蠶血淋巴小RNA進行大規(guī)模測序，分析鑒定桑葉的miRNA，利用生物信息學方法分析家蠶血淋巴中是否存在桑樹源miRNA。采用莖環(huán)PCR、T-A克隆、Sanger測序驗證該結果，并探究了桑樹源miRNA是否廣泛存在家蠶的組織中。利用微滴數字PCR技術檢

2、測人工合成的miR166b是否也能進入家蠶的體內。另外，采用轉錄組測序、RT-PCR以及表型調查等方法探究桑樹源miR166b對家蠶的生理活動是否有影響。獲得的主要結果如下:
　　1.利用高通量測序技術鑒定川桑的保守miRNA和新miRNA
　　研究采用solexa技術對川桑3個組織（葉、皮和雄花）的小RNA文庫進行高通量測序，利用生物信息學方法，在川桑的3個組織中共鑒定到85個保守miRNA，分屬于31個家族，262個新m

3、iRNA，分別位于川?；蚪M的371個位點上。采用miRNA靶基因預測的方法，20個保守miRNA家族和113個新miRNA分別預測得到89和332個候選靶基因。
　　川桑保守miRNA的進化保守性分析表明川桑大部分miRNA同其它雙子葉和單子葉植物都較為保守。表達水平分析表明川桑保守miRNA家族的表達水平范圍較廣，具有組織傾向性，葉、皮和雄花組織傾向性表達的miRNA家族分別是5個、5個和8個。RT-PCR結果顯示部分保守mi

4、RNA在葉、皮和雄花的表達譜模式同高通量測序結果一致。
　　川桑新miRNA的特點分析表明，新miRNA前體長度范圍與其它已報道的植物相似;最小折疊自由能平均值低于tRNA和rRNA;鑒定到90個新miRNA對應miRNA*序列，進一步表明了新預測miRNA的準確性。表達水平分析表明新miRNA的表達水平也較為寬廣，總體表達水平遠低于保守miRNA，也具有組織傾向性，其中58個、84個和72個分別具有葉、皮和雄花組織傾向性表達模式

5、。選取部分新miRNA進行RT-PCR驗證，結果顯示其在3個組織的表達模式與高通量測序的結果一致。
　　本部分研究通過高通量測序技術鑒定了川桑保守miRNA、新miRNA以及其候選靶基因，對桑樹miRNA的功能分析和家蠶體內是否有桑樹源miRNA的研究提供了基礎數據。
　　2.栽培桑miRNA的鑒定及其與野生桑miRNA的比較分析
　　之前我們對川桑（野生桑）葉片進行了miRNA的鑒定，但是家蠶攝食的是栽培桑桑葉，于是

6、本部分研究采用solexa測序對兩個栽培桑品系（湖桑32號和倫教109）葉片miRNA進行鑒定并比較分析了野生桑和栽培桑之間的差異。
　　結合miRBase的第21版釋放數據，共鑒定到兩個栽培桑品系以及川桑的保守miRNA分別是192個、180個和151個，湖桑32號和倫教109的新miRNA分別是177個和167個。三者共有的保守miRNA是137個，表明野生桑和栽培桑的miRNA種類差異不大。表達水平差異分析表明野生桑（川桑）

7、和栽培桑（湖桑32號和倫教109）的差異保守miRNA共有100個，其中75個miRNA在栽培桑中上調表達，25個下調表達;差異新miRNA共有266個，其中89個新miRNA在栽培桑中上調表達，177個下調表達，雖然差異表達的miRNA數目較多，但是表達量高的差異表達miRNA較少，RPM值在100以上的只有6個。說明了野生桑（川桑）和兩個栽培桑miRNA整體的表達水平大同小異。同時我們也分析了這些差異miRNA可能與川桑和栽培桑之間

8、的關系。差異表達miRNA對應靶基因的KEGG分析表明，野生桑和栽培桑參與碳水化合物和氨基酸代謝的差異miRNA對應靶基因數目最多，暗示了二者在糖類和蛋白質含量可能有所差異。
　　本部分研究表明栽培桑和野生桑在miRNA的種類和表達水平上大同小異，說明用川桑數據進行下一步的分析是可行的。
　　3.家蠶血淋巴小RNA種類的鑒定和分析
　　為探究家蠶體內是否有桑樹源miRNA，本部分研究采用高通量測序技術對家蠶血淋巴小RN

9、A文庫進行測序，通過分析共鑒定的小RNA種類包括:桑樹源miRNA，家蠶的miRNA、tsRNA、rRNA、snRNA以及snoRNA。
　　家蠶血淋巴小RNA和桑樹miRNA的比對分析結果表明，家蠶血淋巴中有8個桑樹源miRNA，包括mno-miR166b，mno-miR166c，mno-miR167e，mno-miR396b, mno-miR156c,mno-miR398, mno-miR162, mno-miR159a，其總

10、體的表達量非常低，最高的reads數只有11。該結果表明家蠶攝食桑葉后，部分桑樹miRNA能夠穿過家蠶的消化道進入家蠶的體內，但是由于檢測到的桑樹源miRNA含量較低，也不能排除這一結果是由于高通量測序帶來的誤差或是樣品收集過程中桑葉帶來的污染。
　　家蠶血淋巴小RNA和已報道的家蠶miRNA的比對結果表明，家蠶血淋巴中共有133個已知miRNA和6個新miRNA，其整體表達水平遠遠低于其在家蠶的其它組織（例如脂肪體、絲腺等），推

11、測其原因可能是血淋巴中的miRNA大部分來源于血細胞中的miRNA和組織細胞分泌的miRNA。
　　此外，我們還發(fā)現，血淋巴小RNA長度分布結果表明血淋巴中含量最高的小RNA是tRNA衍生的小RNA(tsRNA)，分別占據兩個小RNA文庫總量的86.99％和92.04％，長度分布主要在25-28nt區(qū)間。tsRNA分類結果表明，家蠶血淋巴中含量前2位分別是tsRNA-Asp(GAC)和tsRNA-His(CAC)。以上兩種tsRN

12、A在其對應的tRNA的位置分布結果表明幾乎99％都位于tRNA的5'端。此外，tsRNA-Asp(GAC)和tsRNA-His(CAC)在家蠶體液和組織中的表達水平分析表明，這兩個tsRNA僅在家蠶血淋巴中高量表達，在全蠶、脂肪體和絲腺中極微量表達，以上結果暗示家蠶的tsRNA-Asp(GAC)和tsRNA-His(CAC)可能在血淋巴中起著特殊的作用。
　　該研究利用二代測序法和生物信息學方法對家蠶血淋巴的小RNA種類進行了分析

13、，對家蠶血淋巴中miRNA、tsRNA的功能研究和桑樹源miRNA跨界轉運到家蠶的研究提供了數據支撐。
　　4.桑樹miR166b對家蠶的跨界調節(jié)研究
　　家蠶血淋巴小RNA和桑樹miRNA的比對分析結果顯示，家蠶血淋巴中有8個桑樹源miRNA。本研究對桑樹源的miRNA是否真正進入了家蠶體內進行了驗證，并研究了進入家蠶體內的桑樹源miRNA對家蠶的生理功能是否有作用。
　　莖環(huán)PCR、T-A克隆和Sanger測序結果

14、顯示:7個桑樹源的miRNA(mno-miR166b,mo-miR166c,mno-miR167e,mno-miR396b,mno-miR156c,mno-miR398和mno-miR159a)進入了家蠶的血淋巴中，通過與對照miR169a的測序頻率比較，證明了這7個桑樹源miRNA不是來自大規(guī)模測序的污染，也非樣品收集過程中桑葉帶來的污染，是真實存在家蠶的血淋巴中。此外，克隆測序結果表明桑樹源的4個miRNA(mno-miR166b，

15、mno-miR166c，mno-miR167e和mno-miR396b)能進入家蠶的9個組織（脂肪體、絲腺、腦、卵巢、中腸、涎腺、精巢、前胸腺和馬氏管）。微滴數字化PCR(ddPCR)結果表明人工合成的miR166b在被家蠶攝食后，能進入家蠶的血淋巴和脂肪體中，然而該miRNA僅能在家蠶體內持續(xù)保持約3小時便被降解。
　　RT-PCR結果表明miR166b的11個候選靶基因在喂食人工合成miR166b前后的表達水平沒有顯著變化，暗

16、示了miR166b在家蠶中沒有功能或是這11個基因不是miR166b的靶基因。轉錄組測序的進一步分析結果表明，5齡2天的雄蠶被喂食人工合成的miR166b之后與其被喂食前相比較，全蠶的整個轉錄組共有30個差異表達基因，其中17個上調基因，13個下調基因?；蚍诸惐砻鞴灿?2個基因涉及免疫和響應壓力，占據30個差異表達基因的40％，我們推測通過攝食進入家蠶體內的人工合成miR166b可能激發(fā)了家蠶抵抗非自我分子。6個差異表達基因(DEGs

17、)的RT-PCR結果顯示其中有3個DEGs在喂食了人工合成miR166c和/或miR167e后上調和/或下調表達，該結果暗示了轉錄組測序得到的12個與免疫和壓力相關的DEGs表達水平變化的發(fā)生不是由于喂食的miR166b的序列引起，而是miR166b這種核酸物質造成的。表型調查表明喂食miR166b和沒有喂食miR166b的家蠶在幼蟲重量的增加倍數、進入熟蠶期的比率、致死率、蛹和繭殼重量等方面沒有顯著差異。此外，我們還預測了miR166

18、b在30個差異表達基因中是否有作用的靶位點，結果表明沒有重要的靶位點。以上結果都表明了桑樹miR166b對家蠶的生理活動基本沒有影響。
　　綜合上文高通量測序、克隆測序實驗、ddPCR檢測以及喂食人工合成miR166b的轉錄組分析、RT-PCR檢測和表型調查的結果，我們認為桑葉中的miRNAs能通過家蠶攝食桑葉的方式進入家蠶的血淋巴和9個組織中，人工合成miR166b也能進入家蠶的血淋巴和脂肪體，然而喂食的人工合成miR166b對