【目的】长链非编码RNA(lncRNA)在真核生物的基因表达、表观遗传和细胞周期调控等方面发挥重要功能。本研究旨在探究意大利蜜蜂(Apis mellifera ligustica,简称意蜂)工蜂中肠发育过程中lncRNA的表达谱及其作用。【方法】利用RNA-seq技术和链特异性建库方法对意蜂7和10日龄工蜂中肠(Am7、Am10)进行深度测序,下机的原始数据经过Perl脚本过滤得到高质量有效读段。利用bowtie工具将有效读段比对核糖体数据库,进一步利用Top Hat2软件将未比对到核糖体数据库上的数据比对到参考基因组。利用CPC和CNCI软件对转录本的编码能力进行预测。通过RT-PCR对部分lncRNA进行鉴定。利用edgeR软件进行差异表达lncRNA(DElncRNA)分析,进而预测lncRNA的上下游基因,并对上下游基因进行GO及KEGG代谢通路富集分析。联用RNAhybrid、Miranda和Target Scan软件预测DElncRNA靶向结合的mi RNA及mi RNA靶向结合的靶基因,并通过Cytoscape软件构建DElncRNAs-mi RNAs-m RNAs的调控网络。最后,通过RT-qPCR验证测序数据的可靠性。【结果】Am7和Am10的深度测序分别获得134 802 058和147 051 470条原始读段,经严格过滤分别得到134 166 157和146 293 288条有效读段;共得到3 890个DElncRNA,包括2 005个上调lncRNA与1 885个下调lncRNA。RT-PCR验证结果显示共有8个lncRNA能扩增出符合预期的目的片段,表明预测出的lncRNA真实存在。DElncRNA的上下游基因数为1 793个,它们涉及42个GO条目,包括代谢进程、发育进程、细胞进程、应激反应和免疫系统进程等;这些上下游基因还涉及251条代谢通路,包括碳代谢、嘌呤代谢和脂肪酸的生物合成等物质代谢通路,硫代谢、甲烷代谢和氧化磷酸化等能量代谢通路,Hippo信号通路、Wnt信号通路和Notch信号通路等信号通路,溶酶体、内吞作用和泛素介导的蛋白水解等细胞免疫通路,以及MAPK信号通路、Jak-STAT信号通路和NF-kappa B信号通路等体液免疫通路,上述结果表明DElncRNA在意蜂中肠发育过程中参与物质和能量代谢、细胞生命活动和免疫调控。进一步分析发现TCONS_00020918可通过调控西方蜜蜂王浆主蛋白1编码基因在意蜂工蜂中肠的营养吸收、级型分化中发挥功能。DElncRNA的调控网络分析结果显示DElncRNA与mi RNA、m RNA间存在复杂的调控关系,部分DElncRNA处于调控网络的中心位置且能靶向结合较多的mi RNA,也有部分mi RNA可被多个DElncRNA共同靶向,表明这些DElncRNA可能在中肠发育中发挥重要作用。随机挑取5个DElncRNA进行RT-qPCR验证,结果显示它们的表达量变化趋势与测序结果一致,证实了本研究测序数据的可靠性。【结论】差异表达长链非编码RNA(DElncRNA)广泛参与意蜂工蜂中肠的新陈代谢、细胞活动和免疫调控并作为竞争性内源RNA(ce RNA)发挥作用,研究结果为关键lncRNA的筛选和功能研究提供了必要的数据支持。

【目的】环状RNA(circular RNA,circRNA)在可变剪接、转录调控和来源基因的表达调控等方面具有重要功能。本研究旨在探究意大利蜜蜂(Apis mellifera ligustica,简称意蜂)工蜂中肠发育过程中circRNA的表达谱及其发育过程中的差异表达circRNA(differentially expressed circRNA,DEcircRNA),进而在转录组水平探究DEcircRNA在中肠发育中的作用。【方法】基于前期获得的意蜂7和10日龄工蜂中肠样品(Am7和Am10)的全转录组数据,利用find_circ软件从质控后的数据中预测circRNA。采用RPM算法归一化处理得到circ RNA的表达量。利用DEGseq软件对circ RNA进行差异分析,按照差异倍数(fold change)≥2、P

【目的】前期已对意大利蜜蜂(Apis mellifera ligustica,简称意蜂)7日龄工蜂中肠(Am7)、10日龄工蜂中肠(Am10)进行全转录组测序,本研究基于高质量的组学数据探究中肠发育过程中的差异基因表达谱及其调控网络,以期解析意蜂工蜂中肠发育的分子机理。【方法】根据FPKM(fragments per kilobase of transcript per million mapped reads)算法计算基因表达量,并以|log2 fold change|≥1且P≤0.05作为标准筛选得到差异表达基因(differentially expressed gene,DEG)。利用TargetFinder软件预测ame-miR-6001-3p的靶mRNA。利用相关生物信息学软件,对全部DEG进行GO和KEGG数据库注释。筛选出与AMPK、P13K-Akt、Wnt、cAMP、FoxO、Hippo、mTOR、Jak-STAT、Toll-like受体、TGF-beta、Notch、MAPK和NF-κB 13条信号通路存在富集关系的DEG,以及与ame-miR-6001-3p存在靶向结合关系的DEG,通过Cytoscape软件构建调控网络将其富集关系与调控关系可视化。利用茎环反转录PCR(Stem loop RT-PCR)和实时荧光定量PCR(RT-qPCR)验证ame-miR-6001-3p以及DEG在Am7和Am10中的差异表达情况。【结果】Am7 vs Am10比较组中共有1 038个DEG,包括515个上调基因和523个下调基因。这些DEG涉及细胞进程、代谢进程和催化活性等功能条目,并显著富集在氧化磷酸化、氨基糖与核苷酸糖代谢、脂肪酸代谢和嘌呤代谢等能量和物质代谢通路,表明工蜂中肠内存在旺盛细胞生命与新陈代谢活动。表达量聚类分析发现分别有20、18、15和14个DEG富集在AMPK信号通路、P13K-Akt信号通路、内吞作用和Hippo信号通路。57个DEG与P13K-Akt、Wnt、Jak-STAT等上述13条与生长发育和免疫防御相关的信号通路存在富集关系,且1个DEG可与多条信号通路存在富集关系。调控网络分析结果显示,分别有54个上调基因和44个下调基因可被ame-miR-6001-3p靶向结合;上调基因富集在磷酸肌醇代谢、胰岛素信号通路、Hippo信号通路和谷胱甘肽代谢等43条代谢通路,而下调基因富集在Hippo信号通路、新陈代谢途径、谷胱甘肽代谢和花生四烯酸代谢等20条代谢通路。RT-qPCR结果显示随机挑选的6个DEG的表达量变化趋势与测序数据一致,证实了本研究中基因差异表达真实可靠。此外ame-miR-6001-3p在Am7和Am10内均真实表达,并在Am10中的相对表达量显著较低。【结论】对意蜂工蜂中肠发育过程的DEG表达谱和DEG与ame-miR-6001-3p之间的调控关系,以及DEG的潜在作用进行深入分析和探讨,发现DEG可参与TGF-beta、Wnt、Hippo、Notch、PI3K-Akt、mTOR、AMPK和NF-κB等各类信号通路进而影响中肠的生长发育和免疫防御,DEG可通过与显著下调表达的ame-miR-6001-3p形成复杂的调控网络参与意蜂工蜂中肠发育过程中胰岛素信号通路等多条代谢途径。

【目的】微小RNA(microRNA,miRNA)是一类在转录后水平对mRNA进行负调控的关键调控因子。本研究旨在通过分析意大利蜜蜂(Apis mellifera ligustica,简称意蜂)幼虫肠道发育过程中差异表达miRNA(differentially expressed miRNA,DEmiRNA)及其调控网络,提供miRNA的表达谱和差异表达信息,揭示DEmiRNA在幼虫肠道发育中的作用。【方法】利用small RNA-seq(sRNA-seq)技术对意蜂4、5和6日龄幼虫肠道样品(Am4、Am5和Am6)进行测序,将质控后的数据与西方蜜蜂(Apis mellifera)参考基因组进行比对,然后将比对上的序列标签(tags)注释到miRBase数据库,并利用TPM(tags per million)算法归一化处理所得miRNA的表达量,再通过相关生物信息学软件对miRNA进行表达量聚类、前体二级结构预测及差异表达分析。利用TargetFinder软件预测DEmiRNA的靶基因,使用Blast软件对靶基因进行GO和KEGG数据库注释,进而通过Cytoscape软件构建mi RNA-m RNA的调控网络。采用茎环实时荧光定量PCR(Stem-loop RT-qPCR)验证测序数据的可靠性。【结果】意蜂幼虫肠道样品的测序分别得到10 841 644、12 037 678和9 230 496条有效序列标签;Am4 vs Am5比较组包含16个上调和10个下调mi RNA,Am5 vs Am6比较组包含5个上调和7个下调mi RNA。其中,novel-m0031-3p为两个比较组所共有,并结合5个与蜕皮激素诱导蛋白相关的靶基因;二者的特有DEmi RNA数分别为25和11个。Am4 vs Am5的26个DEmi RNA结合5 742个靶基因,其中2 725个靶基因可注释到GO数据库中的46个GO term,并主要富集在结合、细胞进程、代谢进程和单组织进程等;Am5 vs Am6的12个DEmi RNA结合3 733个靶基因,且其中2 725个靶基因富集在结合、细胞进程、单组织进程和代谢进程等41个GO term;此外,两个比较组中的DEmi RNA分别有1 046和676个靶基因可注释到116和92条KEGG代谢通路,且Am4 vs Am5比较组的DEmi RNA的靶基因富集在Wnt信号通路、Hippo信号通路、嘌呤代谢和内吞作用等通路上的数量均高于Am5 vs Am6比较组。进一步分析结果显示Am4 vs Am5中的上调和下调mi RNA可分别结合611和85个靶基因,其中ame-mi R-6052结合的靶基因数最多,可通过结合5个靶基因,参与对细胞色素P450的调控;mi R-281-x结合49个靶基因,并间接调控组氨酸代谢、TGF-β信号通路以及Hippo信号通路;Am5 vs Am6中的上调和下调mi RNA可分别结合43和431个靶基因,其中mi R-iab-4-x结合靶基因数量最多,并广泛参与调控背腹轴的形成、Hippo信号通路、Wnt信号通路、Fox O信号通路、Notch信号通路以及m TOR信号通路等与生长发育相关的通路。调控网络分析结果表明,DEmi RNA与靶基因间形成较为复杂的调控网络,DEmi RNA居于调控网络的中心位置,而m RNA位于调控网络的外周。最后,通过茎环实时荧光定量PCR对随机选取的3个DEmiRNA进行验证,结果证实了测序数据的可靠性。【结论】在全基因组水平对意蜂幼虫肠道的DEmi RNA及其靶基因进行预测和分析,并对DEmi RNA的调控网络进行构建及分析,发现意蜂幼虫通过调节包括ame-mi R-6052、miR-iab-4-x、mi R-281-x和novel-m0031-3p在内的多个mi RNA的表达水平对肠道生长和发育进行调控,研究结果不仅提供了意蜂肠道发育过程的mi RNA表达谱及差异表达信息,也为阐明意蜂幼虫肠道的发育机理打下基础。

基于前期已获得的意大利蜜蜂(Apis mellifera ligustica)工蜂中肠转录组数据,在全转录组水平上对意大利蜜蜂的可变剪接基因(alternatively spliced gene, ASG)进行鉴定和分析,旨在揭示ASG在宿主响应东方蜜蜂微孢子虫(Nosema ceranae)胁迫中的作用.在正常意大利蜜蜂7和10日龄工蜂中肠(Am7CK、Am10CK)和东方蜜蜂微孢子虫胁迫的意大利蜜蜂7和10日龄工蜂中肠(Am7T、Am10T)样品中检测到可变剪接事件数分别为73 383、85 618、79 813和74 693次,对应的基因数分别为9 586、10 063、9 829和9 622个,平均每个ASG上发生的可变剪接事件数分别为7.66、8.50、8.12和7.76次.各样品的可变剪接事件中外显子跨越的数量最多,分别为4 411、5 247、4 993和4 629个.Venn分析结果显示,Am7CK、Am10CK、Am7T和Am10T中的共有ASG数为9 056个,特有ASG数分别为100、153、275和95个.GO分类结果显示,这些共有ASG分布在50个功能条目,特有ASG分别分布于26、30、33和24个功能条目.KEGG代谢通路富集分析结果显示,上述共有ASG富集在325个代谢通路,包括内质网蛋白加工、核糖体和RNA转运等;特有ASG分别富集在31、104、126和22个代谢通路.研究结果提供了意大利蜜蜂响应东方蜜蜂微孢子虫胁迫过程中的ASG数量和种类信息,揭示宿主ASG在胁迫响应中参与了新陈代谢和细胞免疫,为关键ASG的筛选和功能研究建立了数据基础.

【目的】对原种意大利蜜蜂(意蜂)和原种意大利蜜蜂的王浆高产品系(浆蜂)工蜂蛹期头部的蛋白质组进行比较,以探明浆蜂和意蜂3个日龄(13、15、17日龄)蛋白质表达调控方面的差异,为阐明蜜蜂发育生物学提供一定理论基础。【方法】采用双向电泳建立浆蜂和意蜂3个日龄蛹期头部蛋白质组表达谱,获得图谱中蛋白质的种类、表达量、等电点和分子量等信息,并通过质谱分析、数据库检索,鉴定部分差异蛋白。【结果】在13日龄蛹头部,浆蜂表达279个蛋白,意蜂表达268个蛋白,浆蜂和意蜂共有蛋白为212个,而浆蜂和意蜂的特异蛋白分别为67个和56个;到15日龄时,浆蜂表达296个蛋白,意蜂为277个,浆蜂和意蜂共有蛋白点为2...

根据中华蜜蜂(Apis cerana cerana)基因组序列,鉴定抗氧化酶系统的主要组分,并用RNA-Seq数据分析了工蜂行为发育中抗氧化基因在脑组织中的转录水平.研究表明,共有47个抗氧化基因,包括34个抗氧化酶基因和13个硫氧还蛋白基因.比较东方蜜蜂(Apis cerana)、西方蜜蜂(Apis mellifera)和黑腹果蝇(Drosophila melanogaster)的抗氧化基因,抗氧化酶系统的主要组分在进化上保守,但抗氧化酶基因家族的成员存在物种间差异.如东方蜜蜂和黑腹果蝇的SOD基因家族含3个成员,而西方蜜蜂有2个;果蝇有43个GST基因,而东方蜜蜂和西方蜜蜂仅有14和13个...

【目的】比较王浆高产蜜蜂(Apis mellifera L.浆蜂)与原种意大利蜜蜂(Apis mellifera L.原意)雄蜂卵期3d发育阶段蛋白质组的特点。【方法】根据双向电泳所得表达谱,对两个蜂种雄蜂卵期发育的蛋白质的种类、等电点、分子量和表达量进行分析。【结果】浆蜂雄蜂卵期所表达的蛋白数(332、377和339)显著的高于相应日龄原意雄蜂卵期3个日龄所表达的蛋白数(283、305和293)(P<0.05),其中在两蜂种3个日龄共同表达的蛋白数为254个。通过雄蜂卵期3d发育过程蛋白质表达谱的比较,两个蜂种分别检测到274、298和285个共有蛋白。在这些共有蛋白中,浆蜂雄蜂分别有54、...

【目的】通过对中华蜜蜂(中蜂)和意大利蜜蜂(意蜂)雄蜂胚胎发育期3个日龄差异蛋白质组和磷酸化蛋白质组进行分析,以探明中蜂、意蜂雄蜂胚胎发育特征并了解雄蜂胚胎发育生物学,为蜜蜂遗传改良提供理论依据。【方法】采用双向电泳方法建立中蜂、意蜂雄蜂胚胎发育3个日龄的蛋白质组及磷酸化蛋白质组表达谱,获得图谱中蛋白质的分子量、等电点和表达量等信息,并对部分差异蛋白质进行质谱分析、功能鉴定和生物信息学分析。【结果】在雄蜂胚胎发育的第1、2、3日龄,中蜂表达332、362、340个蛋白质,其中111、117、112个蛋白质发生了磷酸化修饰;意蜂表达302、331、311个蛋白质,其中有107、110、106个发...

【目的】对王浆高产蜜蜂(浆蜂)(A.m.ligustica)和原种意大利蜜蜂(原意)(A.m.ligustica)1、3、6日龄工蜂咽下腺进行蛋白质组研究,揭示咽下腺在此阶段的发育特征。【方法】采用双向电泳方法对浆蜂和原意咽下腺进行蛋白质组研究,通过与已鉴定蛋白功能的咽下腺和蜂王浆蛋白质组图谱比较,推断部分蛋白的功能。【结果】浆蜂咽下腺在3个日龄表达的蛋白数(210、192、230)分别显著的高于原意(169、188、212),两蜂种均在6日龄表达蛋白最多(P<0.05)。浆蜂咽下腺3个日龄表达的共有蛋白数为119个,其中21个蛋白表达量随咽下腺的发育显著上调,14个蛋白显著下调;原意咽下腺3...

温度是影响蜜蜂发育的关键因子之一.为研究低温胁迫对西方蜜蜂幼虫基因表达的影响,利用高通量转录组测序技术分别对最适发育温度和低温处理的西方蜜蜂工蜂6日龄幼虫的基因表达情况进行分析,获得差异表达基因(DEGs),并对DEGs进行GO功能分类和KEGG富集分析.结果显示,低温处理后共检测到98个DEGs.GO功能分类显示,DEGs富集在21条二级GO功能注释上,富集最多DEGs的GO功能注释为结合、细胞进程、代谢进程、催化活性;KEGG富集分析显示,上调和下调DEGs分别富集在25和3条通路上,主要包括甘氨酸-丝氨酸-苏氨酸代谢、核黄素代谢、昆虫激素生物合成等代谢通路,以及与信号转导相关的Hippo、Notch等信号通路.qPCR结果显示,转录组测序获得的DEGs IP3K、Ex、CYP18A1、GOD、GCDH表达差异显著.本研究结果揭示西方蜜蜂工蜂6日龄幼虫主要通过提高抗氧化能力和减慢发育速度应对低温.

【目的】20世纪90年代中国成功培育出的王浆高产蜜蜂(Apis mellifera L.)是当今世界最优良的蜂种,使中国的王浆年产量高达2600多吨,占世界总产量的90%以上,稳居世界第一,通过对该蜂种不同发育日龄工蜂幼虫的蛋白质组进行研究,以探明其发育机理。【方法】采用双向电泳法对王浆高产蜜蜂(Apis mellifera L.)不同发育期的工蜂幼虫进行蛋白质组研究。【结果】在幼虫期6d的发育过程中,2日龄、4日龄、6日龄幼虫中分别检测到了262、418、194个蛋白点。这些蛋白的分子量在12.2～88.2kD的较大范围之间,等电点在pH4.00～9.24之间。有84个蛋白在幼虫期整个发育过...

【目的】微小RNA(microRNA,miRNA)是一类重要的基因表达调控因子,可影响宿主与病原间的互作过程。蜜蜂球囊菌(Ascosphaera apis)是一种特异性侵染蜜蜂幼虫的致死性真菌病原。本研究旨在对意大利蜜蜂(Apis mellifera ligustica,简称意蜂)幼虫肠道在球囊菌胁迫前期的差异表达miRNA(differentially expressed miRNA,DEmiRNA)及其靶基因进行深入分析,在miRNA组学水平探究意蜂幼虫在球囊菌侵染前期的胁迫应答,并通过构建显著DEmiRNA的调控网络筛选出与宿主应答相关的关键miRNA。【方法】利用small RNA-seq(sRNA-seq)技术对正常及球囊菌侵染的意蜂4日龄幼虫肠道(AmCK和AmT)进行高通量测序,首先对原始数据进行质控和评估,随后将过滤后的数据与西方蜜蜂(Apis mellifera)参考基因组进行比对;将比对上的序列标签(tags)注释到miRBase数据库,得出已知miRNA的表达量;通过TPM(tags per million)算法对各样本中miRNA的表达量进行归一化处理,以|log_2 fold change|≥1且P≤0.05作为标准筛选得到显著DEmiRNA;利用TargetFinder软件预测显著DEmiRNA的靶基因,并对其进行GO和KEGG代谢通路(pathway)富集分析。利用Cytoscape软件对miRNA-mRNA调控网络进行可视化。最后,利用茎环反转录PCR(Stem-loop RT-PCR)和荧光定量PCR(qPCR)验证测序数据的可靠性。【结果】AmCK和AmT样品的测序分别得到13 553 302和10 777 534条原始读段(raw reads),经严格过滤后得到的有效读段(clean reads)数分别为13 186 921和10 480 913条。各样品的生物学重复间的Pearson相关性系数分别在0.9822和0.9508以上。共有10个显著DEmiRNA,包括4个上调miRNA和6个下调miRNA。显著DEmiRNA在AmT的整体表达水平低于AmCK。10个显著DEmiRNA可靶向结合3 788个靶基因,其中上调miRNA的1 240个靶基因可注释到GO数据库中的39个GO条目,主要富集在结合、细胞进程、代谢进程和应激反应等;下调miRNA的749个靶基因可注释到34个GO条目,主要富集在细胞进程、结合、代谢进程和应激反应等。KEGG数据库注释结果显示,上调miRNA和下调miRNA的靶基因分别注释到95和66条代谢通路,富集基因数最多的分别是Wnt信号通路、Hippo信号通路、光传导以及内吞作用、磷脂酰肌醇信号系统、嘌呤代谢。对于上调和下调miRNA,分别有31和52个靶基因注释到内吞作用,15和7个靶基因注释到泛素介导的蛋白水解,11和5个靶基因注释到Jak-STAT信号通路,1和3个靶基因注释到MAPK信号通路。显著DEmiRNA与靶mRNA之间形成复杂的调控网络,7个显著DEmiRNA靶向结合96个与Wnt信号通路相关的mRNA,8个显著DEmiRNA靶向结合55个与内吞作用相关的mRNA。Stem-loop RT-PCR和qPCR结果验证了测序数据的可靠性。【结论】对意蜂幼虫肠道在球囊菌侵染前期的DEmiRNA及其靶基因进行预测和分析,并构建和分析了DEmiRNA-mRNA调控网络,研究结果提供了宿主miRNA的表达谱和差异表达信息,揭示了DEmiRNA通过调控细胞生命活动、新陈代谢以及部分细胞和体液免疫等生物学过程参与宿主的胁迫应答。miR-4331-y、miR-4968-y、miR-8440-y、novel-m0023-5p和novel-m0025-3p共同参与了宿主的Wnt信号通路和内吞作用的调控,可作为白垩病治疗的潜在分子靶标。

【目的】对原种意大利蜜蜂(Apis mellifera ligustica Spinola,1806)(意蜂)雄蜂的3d卵期发育进行蛋白质组研究,以探明该蜂种雄蜂卵期蛋白质表达调控方面的一些特点,从而揭示其发育的分子机理。【方法】采用双向电泳法对意蜂雄蜂卵期蛋白质组进行研究。【结果】3d的卵期发育中,按顺序分别检测到200、242、233个蛋白,这些蛋白的分子量在12.42～169.60kD,等电点为4.50～9.00。164个蛋白在3d卵期发育过程中均有表达,其中79个蛋白在3个日龄的表达量无显著差异(P>0.05),7个共有蛋白表达量随卵期发育呈显著上调(P<0.05)趋势,而4个共有蛋白...

【目的】利用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 Prism7软件统计各组肠道样品中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参与调节发育和免疫相关途径。