中华蜜蜂幼虫肠道中微小RNA的鉴定及分析

【目的】利用small RNA-seq（sRNA-seq）技术和生物信息学方法对中华蜜蜂（Apis cerana cerana,简称中蜂）幼虫肠道的微小RNA（microRNA,miRNA）进行全转录组鉴定和分析,旨在丰富中蜂的miRNA信息,并为深入研究miRNA调控中蜂幼虫肠道发育的分子机理提供依据。【方法】利用sRNA-seq技术对中蜂4、5和6日龄幼虫肠道样品（Ac1、Ac2和Ac3）进行测序,通过数据质控获得有效标签序列（clean tags）。采用Blast工具将clean tags连续比对东方蜜蜂（Apis cerana）基因组和miRBase数据库,以鉴定保守miRNA和新miRNA。采用TPM法对miRNA的表达量进行归一化处理。通过GraphPad Prism 7软件统计各组肠道样品中sRNA占比、miRNA长度分布及首位碱基偏向性。利用相关软件预测上述miRNA靶向的mRNA并进行GO和KEGG数据库注释。进一步根据靶向结合关系构建和分析注释到发育和免疫相关通路的基因及其靶向miRNA的调控网络,并利用Cytoscape软件进行可视化。利用茎环反转录PCR（Stem-loop RT-PCR）、分子克隆和Sanger测序验证miRNA的表达和序列的真实性。【结果】共鉴定到中蜂的371个保守miRNA和64个新miRNA;这些miRNA的长度介于18—25 nt且首位碱基主要偏向于U;上述miRNA共靶向14 750条mRNA,涉及离子结合、金属离子结合、细胞膜、细胞膜组件和单一有机体进程等2 270个GO条目,以及内吞作用、细胞凋亡、mTOR信号通路、RNA转运和昆虫激素的生物合成等332条KEGG通路。进一步分析结果显示156个miRNA与注释到Wnt、Hippo、Notch和mTOR等生长发育相关通路的67个靶基因存在调控关系,145个miRNA与注释到Toll、Imd/JNK、Jak-STAT和抗菌效应因子等免疫相关途径的21个靶基因存在调控关系。Stem-loop RT-PCR 结果显示miR-8-y、miR-9-z、miR-14-y、miR-281-y、miR-283-x和miR-306-x均能扩增出预期的特异性片段;Sanger测序结果显示上述6个miRNA的序列与深度测序结果一致。【结论】提供了中蜂miRNA的数量、结构特征和表达谱;揭示中蜂幼虫肠道miRNA潜在调控诸多生命进程与细胞活动;中蜂幼虫肠道的部分miRNA可通过靶向结合相应的mRNA参与调节发育和免疫相关途径。