基于前期已获得的意大利蜜蜂(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的筛选和功能研究建立了数据基础.

【目的】通过对意大利蜜蜂(Apis mellifera ligustica,简称意蜂)工蜂中肠响应东方蜜蜂微孢子虫(Nosema ceranae)胁迫的差异表达基因(differentially expressed gene,DEG)及免疫通路进行深入细致的分析,揭示宿主的细胞和体液免疫应答,为关键免疫应答基因的筛选及功能研究打下基础。【方法】基于前期获得的正常及东方蜜蜂微孢子虫胁迫的意蜂7日龄和10日龄工蜂中肠(Am7CK、Am7T、Am10CK、Am10T)转录组数据,按照FDR≤1,P≤0.05和|log2 fold change|≥1的标准筛选出各组的DEG,利用相关生物信息学软件对DEG进行Pearson相关性分析、Venn分析、GO分类和KEGG代谢通路富集分析,进而对免疫通路富集基因进行统计和分析,利用实时荧光定量PCR(real-time quantitative PCR,RT-q PCR)验证转录组数据的可靠性。【结果】差异分析结果显示,Am7CK vs Am7T比较组包含472个上调基因和385个下调基因;Am10CK vs Am10T比较组包含611个上调基因和360个下调基因。Venn分析结果显示,两个比较组特有的DEG分别为739和853个,共有的DEG为118个。GO分类结果显示,Am7CK vs Am7T比较组中上调和下调基因分别涉及23和29个功能条目,其中富集上调基因数最多的前5位分别为结合、催化活性、代谢进程、细胞进程和单组织进程;富集下调基因数最多的前5位分别为代谢进程、单组织进程、催化活性、细胞进程和结合。Am10CK vs Am10T比较组中上调和下调基因分别涉及36和26个功能条目,其中富集上调基因数最多的前5位分别为单组织进程、结合、细胞进程、催化活性和代谢活性;富集下调基因数最多的前5位分别为结合、细胞进程、催化活性、代谢进程和单组织进程。KEGG代谢通路富集分析结果显示,Am7CK vs Am7T比较组中上调和下调基因分别富集在38和33条代谢通路,富集上调基因数最多的前5条分别为胆汁分泌、内质网蛋白加工、泛素介导的蛋白水解、PI3K-Akt信号通路和神经营养因子信号通路;富集下调基因数最多的前5条分别为胞质DNA传感通路、嘌呤代谢、嘧啶代谢、RNA聚合酶和核糖体;涉及泛素介导的蛋白水解等3条细胞免疫通路,以及PI3K-Akt信号通路等7条体液免疫通路。Am10CK vs Am10T比较组中上调和下调基因分别富集在54和43条代谢通路,富集上调基因数最多的前5条分别为Hippo信号通路、药物代谢-细胞色素P450、P450对外源物质的代谢、泛素介导的蛋白水解和鞘脂类代谢;富集下调基因数最多的前5条分别为m RNA监测、鞘脂类信号通路、果糖和甘露糖代谢、半乳糖代谢和鞘脂类代谢;涉及泛素介导的蛋白水解等7条细胞免疫通路,以及NF-κB信号通路等2条体液免疫通路。RT-q PCR验证结果显示6个随机挑选的DEG的表达量变化趋势与测序结果一致,证实了本研究中测序数据的可靠性。进一步分析发现意蜂7日龄和10日龄工蜂中肠的NF-κB信号通路均被东方蜜蜂微孢子虫激活,随即启动了3种抗菌肽apidaecin、defensin-1和hymenoptaecin的合成,体现了它们在宿主抵御东方蜜蜂微孢子虫入侵中的重要性。【结论】在转录组水平解析了意蜂工蜂中肠对东方蜜蜂微孢子虫的免疫应答,揭示宿主在胁迫前期同时作出细胞和体液免疫应答,前者可能在抵御病原入侵方面发挥主要作用;宿主的细胞免疫在胁迫后期持续增强,但体液免疫较大程度地减弱;泛素介导的蛋白水解通路及富集DEG,NF-κB信号通路及富集DEG,以及抗菌肽编码基因apidaecin、defensin-1和hymenoptaecin可能在意蜂工蜂对东方蜜蜂微孢子虫的免疫应答中起到关键作用。

基于前期低温胁迫3日龄封盖子不同时间的转录组数据,对可变剪接基因(ASG)作数量、种类统计和表达谱等分析,进一步解释低温对蜜蜂化蛹影响的分子机制;同时利用rMATS软件及Omicshare在线工具对T18、T36处理组(在20℃低温下分别处理18和36 h)及CK(没有低温处理)的差异基因进行可变剪接(AS)事件分析.结果表明,在各组别(CK与T18、CK与T36、T18与T36)中鉴定出4 971个ASG,产生58 017个AS事件,以外显子跳跃、外显子选择性跳跃和可变3′端剪接3种剪接类型为主.各组别间的共有ASG为4 027个,其基因本体(GO)分类结果显示,这些ASG富集在41条GO条目上;KEGG代谢通路富集分析显示,这些ASG共富集在132条代谢通路上,极显著富集的有Wnt、Hippo、TGF-beta和胞吞作用等信号通路.上述结果表明,低温胁迫可通过ASG介导代谢、生长发育和免疫等多重功能调控蜜蜂的化蛹过程.

温度是影响蜜蜂发育的关键因子之一.为研究低温胁迫对西方蜜蜂幼虫基因表达的影响,利用高通量转录组测序技术分别对最适发育温度和低温处理的西方蜜蜂工蜂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日龄幼虫主要通过提高抗氧化能力和减慢发育速度应对低温.

【目的】通过对东方蜜蜂微孢子虫和其他物种的Rad18-like重组与DNA修复蛋白(Rad18-like recombination and DNA repair protein, RRDRP)进行系统进化分析，为进一步研究其功能提供理论依据。【方法】利用Protparam、ProtScale和SWISS-model等软件分析RRDRP的理化性质、亲水性和三级结构。利用SignalP 4.1 Server、NetPhos 3.1 Server、TMHMM和SOPMA等软件预测RRDRP的信号肽、磷酸化位点、跨膜结构域及二级结构。使用PSORTⅡ软件预测RRDRP蛋白的亚细胞定位。使用MEME软件预测东方蜜蜂微孢子虫和其他物种RRDRP的保守基序。通过Mega 11.0软件对东方蜜蜂微孢子虫和其他物种的RRDRP进行氨基酸序列多重比对，并采用邻接法构建系统进化树。【结果】东方蜜蜂微孢子虫RRDRP基因含有2 928个核苷酸，可编码975个氨基酸，分子质量约11.50 ku,脂溶系数93.58,等电点6.63,分子式C_(5135)H_(8253)N_(1377)O_(1545)S_(34),平均亲水系数为-0.627,可同时定位于细胞核、细胞质、过氧化物酶体和细胞骨架；RRDRP包含55个磷酸化位点、678个α-螺旋、99条延长链、22个β-转角及176个无规则卷曲；RRDRP不含信号肽与跨膜结构域。东方蜜蜂微孢子虫、兔脑炎微孢子虫、颗粒病微粒子虫等11个物种的RRDRP中均含有5个相同的保守基序(Motif 1、Motif 2、Motif 3、Motif 4和Motif 5)。东方蜜蜂微孢子虫与颗粒病微粒子虫在进化树上聚为一支，同源性较高。【结论】东方蜜蜂微孢子虫RRDRP是潜在的亲水性蛋白和胞内蛋白，在东方蜜蜂微孢子虫以及其他微孢子虫和真菌中高度保守。

蜜蜂蜂子发育具有狭温性,偏离最适温度影响蜂子正常发育.为了揭示温度胁迫对蜂子能量消耗的影响,本研究测定了29、31、35、37、38℃条件下,意大利蜜蜂工蜂封盖子不同发育阶段的能量消耗.结果表明:在29℃条件下,封盖子能量消耗受到抑制,但完成封盖期发育会增加糖、蛋白和脂肪的消耗量;31℃条件下,封盖子发育过程中的总糖和总蛋白的消耗也受到抑制,但完成封盖期发育消耗的总能量与35℃相同;37℃条件对封盖子能量消耗无影响;38℃高温胁迫抑制封盖子能量的消耗.

利用相关生物信息学软件对东方蜜蜂微孢子虫(Nosema ceranae)的DNA拓扑异构酶3(DTⅢ)进行了理化性质、分子特性、保守基序、结构域和系统进化的分析。结果表明：DTⅢ包含805个氨基酸，分子式为C_(4267)H_(6719)N_(1117)O_(1240)S_(34),分子质量约为94.60 ku,脂溶系数为81.84,等电点为9.08,平均亲水系数为-0.595,亲水氨基酸多于疏水氨基酸，不含典型信号肽；含有41个丝氨酸磷酸化位点，23个酪氨酸磷酸化位点及14个苏氨酸磷酸化位点；包含310个α-螺旋，117个β-折叠，37个β-转角及341个无规则卷曲；可同时定位于细胞核、细胞质和线粒体；东方蜜蜂微孢子虫与其他10种微孢子虫的DTⅢ均含有5个相同的保守基序(Motif 1、Motif 2、Motif 3、Motif 4和Motif 5)和2个相同的结构域(Toprim和Topoisom＿bac);通过Mega 11.0软件构建东方蜜蜂微孢子虫和其他物种的DTⅢ的系统进化树，结果显示东方蜜蜂微孢子虫与绒鳌蟹肝孢虫(Hepatospora eriocheir)的DTⅢ聚为一支，进化距离最近。研究结果明确了东方蜜蜂微孢子虫DTⅢ的理化性质和分子特性，并揭示东方蜜蜂微孢子虫与绒鳌蟹肠孢虫的DTⅢ亲缘关系最近，DTⅢ在东方蜜蜂微孢子虫和其他微孢子虫中具有较高的保守性。

【目的】长链非编码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的筛选和功能研究提供了必要的数据支持。

【目的】研究意大利蜜蜂(Apis mellifera)工蜂幼虫饲料中适宜色氨酸水平,为探明意大利蜜蜂工蜂幼虫发育阶段的色氨酸营养需要提供理论依据。【方法】选用1日龄工蜂幼虫1 008只,随机分为7组,每组3个重复,每个重复48只幼虫,分别饲喂色氨酸水平为7.84、8.84、9.84、10.84、11.84、12.84和13.84 mg·g-1的7种日粮。取5日龄工蜂幼虫虫体测定色氨酸羟化酶(tryptophan hydroxylase,TPH)基因与5羟色胺受体(5-hydroxytryptamine receptors,5-HTR)1、2α、2β、7基因表达量;取6日龄工蜂幼虫虫体测定总蛋白...

对东方蜜蜂微孢子虫RCDP 42基因进行分子克隆，并通过生物信息学解析RCDP 42蛋白的分子特性，进而对东方蜜蜂微孢子虫和其他物种的RCDP 42蛋白进行保守基序和结构域预测及分析。结果表明，成功克隆到东方蜜蜂微孢子虫RCDP 42基因的CDS区，RCDP 42基因在孢子中真实表达，含有540个核苷酸，可编码179个氨基酸；RCDP 42蛋白的分子质量约为20.3 ku,分子式为C_(915)H_(1416)N_(230)O_(274)S_(10),脂溶系数为86.03,等电点为6.57;包含的亲水氨基酸数量多于疏水氨基酸，但不含信号肽与跨膜结构域，说明RCDP 42可能为亲水性蛋白、胞内蛋白和非跨膜蛋白，可能包含16个磷酸化位点、63个α-螺旋、43条延长链、11个β-转角及62个无规则卷曲。在东方蜜蜂微孢子虫、蚱蜢脑炎微孢子虫和肠脑炎微孢子虫的RCDP 42蛋白中均鉴定到3个保守基序(motif 1、motif 2和motif 3)与1个结构域(Rho)。

【目的】前期已对意大利蜜蜂(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形成复杂的调控网络参与意蜂工蜂中肠发育过程中胰岛素信号通路等多条代谢途径。

利用前期获得的中华蜜蜂(Apis cerana cerana)幼虫肠道转录组数据对东方蜜蜂(Apis cerana)基因组的已注释基因进行结构优化、对未注释的新基因进行预测、分析和鉴定:利用Cufflink软件对reads进行转录本重构,通过与参考转录本序列进行比较,共对3 366个东方蜜蜂的基因结构进行了优化;利用Cuffcompare软件将重构转录本与参考基因组进行比对,共预测出527个东方蜜蜂的新基因,其中234个新基因在Nr数据库中存在注释信息;随机选取12个新基因进行PCR鉴定,有11个能够扩增出符合预期的目的片段,表明多数预测出的新基因真实存在;GO分类和KEGG代谢通路分析结果显示上述新基因可能在东方蜜蜂的生长发育、物质代谢、遗传信息传递等方面具有重要功能。本研究结果丰富和完善了东方蜜蜂基因组的基因结构与功能注释信息,也为新基因的功能研究打下基础。

【目的】明确东方蜜蜂微孢子虫(Nosema ceranae)对密林熊蜂(Bombus patagiatus)的侵染性及致病机理。【方法】采用传统生物学和电镜超微结构观察的方法,结合qPCR定量分析对东方蜜蜂微孢子虫在密林熊蜂上的致病机理进行探讨。【结果】感染初期工蜂除取食减少和行动迟缓外无明显外观染病特征,感染后期工蜂萎靡、衰弱、飞翔无力。解剖后镜检发现中肠仅存少量孢子,但充满大量细菌;熊蜂肠道组织切片发现东方蜜蜂微孢子虫侵染中肠上皮细胞后导致核膨大并变形、线粒体体积变小甚至解体,内质网紊乱,但孢子只侵染寄主细胞质而不侵入细胞核,最终导致线粒体解体,细胞破裂;qPCR定量分析得出接种后3—4 ...

【目的】白垩病是困扰养蜂生产的顽疾。目前,尚无利用二代测序技术研究中华蜜蜂(Apis cerana cerana,简称中蜂)幼虫白垩病的报道。本研究利用RNA-seq技术对健康(Ac CK)及球囊菌(Ascosphaera apis)胁迫的中蜂4日龄幼虫肠道(Ac T)进行深度测序,在转录组水平研究中蜂幼虫在球囊菌胁迫早期的胁迫应答。【方法】通过Illumina Hi Seq 2500平台对Ac CK和Ac T进行双端(PE125)测序,首先对测序数据进行质控和评估,利用edge R软件进行差异表达基因(

【目的】利用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参与调节发育和免疫相关途径。