西方蜜蜂几丁质合成酶基因的生物信息学分析及功能研究

【目的】本研究旨在明确西方蜜蜂Apis mellifera几丁质合成酶(chitin synthase, CHS)基因的分子特性,并揭示CHS在西方蜜蜂工蜂幼虫响应蜜蜂球囊菌Ascosphaera apis胁迫的免疫应答中的功能。【方法】利用相关生物信息学软件预测和分析西方蜜蜂CHS蛋白的分子特性、保守基序和结构域。使用MEGA X软件对西方蜜蜂和其他昆虫CHSs的氨基酸序列进行系统进化分析。采用体外转录法合成CHS和egfp的dsRNA,饲喂蜜蜂球囊菌胁迫的西方蜜蜂工蜂3日龄幼虫以进行RNAi,并利用RT-qPCR检测西方蜜蜂工蜂5日龄幼虫肠道中CHS及宿主响应蜜蜂球囊菌侵染的免疫相关基因abaecin,apidaecin,birc5,defensin-1和PGRP-S2的表达量。【结果】西方蜜蜂CHS蛋白含有1 572个氨基酸,属于20种氨基酸,其中正电荷氨基酸与负电荷氨基酸分别为177和169个,分子量约为178.77 kD,等电点为6.65;含纤维素合酶CESA3催化结构域。在西方蜜蜂和东方蜜蜂A.cerana等共13种昆虫的CHSs中均鉴定到3个Motif (Motif 1,...     

意大利蜜蜂幼虫肠道在球囊菌胁迫后期的差异表达微小RNA及其靶基因分析

【目的】蜜蜂球囊菌Ascosphaera apis(简称球囊菌)是一种特异性侵染蜜蜂幼虫肠道的致死性真菌病原。微小RNA(microRNA,mi RNA)是一种在转录后水平对m RNA进行负调控的关键调控因子。本研究旨在对意大利蜜蜂Apis mellifera ligustica(简称意蜂)幼虫肠道在球囊菌胁迫后期的差异表达mi RNAs(DEmiRNAs)及其靶基因进行深入分析,为揭示DEmiRNAs在胁迫应答中的作用提供重要信息。【方法】利用small RNA-seq(s RNA-seq)技术对正常及球囊菌侵染的意蜂6日龄幼虫肠道(分别表示为Am CK和Am T)进行测序,通过相关生物信息学软件对DEmiRNAs进行预测和分析。利用TargetFinder软件预测DEmiRNAs的靶基因,然后利用Blast软件对靶基因进行GO和KEGG数据库的功能注释。通过Cytoscape软件构建DEmiRNAs与其靶m RNAs的调控网络。通过茎环实时荧光定量PCR(stem-loop RT-qPCR)验证测序数据的可靠性。【结果】Am CK和Am T的测序分别得到9 230 496和10 823 667条有效序列标签; Am CK vs Am T比较组中包含15个上调和6个下调mi RNAs,分别结合3 503和3 252个靶基因,它们可分别注释到40和39个GO terms以及104和99条代谢通路;进一步分析发现调控网络中的17个DEmiRNAs靶向结合116个与丝氨酸蛋白酶相关的m RNAs,14个DEmiRNAs靶向结合54个与泛素介导的蛋白水解相关的m RNAs。Stem-loop RT-qPCR验证结果显示随机选择的4个DEmiRNAs (mi R-251-x,mi R-9277-y,mi R-1672-x和mi R-4968-y)的表达量变化趋势与测序结果一致,表明测序数据真实可信。【结论】本研究率先对意蜂幼虫肠道在球囊菌胁迫后期的DEmiRNAs及其靶基因进行预测与分析,提供了mi RNAs的表达谱和差异表达信息。ame-miR-927b,mi R-429-y和mi R-8440-y等DEmiRNAs可能参与对丝氨酸蛋白酶的调控和对泛素介导的蛋白水解的调控,DEmiRNAs与靶基因之间存在复杂的调控网络关系。DEmiRNAs可参与到意蜂幼虫肠道与球囊菌间的互作。     

环状RNA ame_circ_000115调节蜜蜂球囊菌胁迫下意大利蜜蜂幼虫肠道基因表达

环状RNA (circular RNA, circRNA)是一类新型非编码RNA,已被证实可通过竞争性内源RNA(competing endogenous RNA,ceRNA)调控网络等多种方式参与调节昆虫基因表达和免疫应答。目前,circRNA在蜜蜂免疫应答中的作用尚不清楚。本研究对意大利蜜蜂(Apis mellifera ligustica,简称意蜂)ame＿circ＿000115(简称ac115)的反向剪切(back splicing,BS)位点进行验证。采用RT-qPCR检测蜜蜂球囊菌(Ascosphaera apis)胁迫下意蜂幼虫肠道内ac115的表达谱,利用双荧光素酶报告基因试验验证ac115与ame-miR-13b的结合关系。通过饲喂特异性siRNA对蜜蜂球囊菌胁迫的意蜂幼虫肠道内ac115进行干扰,进而检测干扰ac115对宿主免疫应答相关的6个基因表达量的影响。结果表明,ac115的BS位点真实存在;相较于未接种组,蜜蜂球囊菌接种组4日龄幼虫肠道内ac115的表达量极显著上调(P<0.000 1),5和6日龄幼虫肠道内ac115显著上调表达(P<0.05)...     

意大利蜜蜂工蜂中肠piRNA的鉴定与分析

【目的】Piwi蛋白互作RNA(PIWI-interacting RNA, piRNA)在昆虫的发育和免疫等重要生物学过程发挥重要的调控作用。本研究旨在丰富西方蜜蜂Apis mellifera的piRNA信息,并为进一步探究差异表达piRNA(differentially expressed piRNA, DEpiRNA)调控意大利蜜蜂Apis mellifera ligustica工蜂中肠发育的分子机理提供基础。【方法】基于已获得的意大利蜜蜂7日龄(Am7)和10日龄(Am10)工蜂中肠的small RNA(sRNA)组学数据,对piRNA进行预测和分析。将质控后的sRNA数据比对西方蜜蜂参考基因组,再将比对上的序列标签(tags)进一步比对数据库以滤除rRNA和tRNA等小非编码RNA(non-coding RNA,ncRNA),进而根据piRNA的长度特征鉴定piRNA。采用TPM(tags per million)算法对piRNA表达量进行计算和归一化处理。根据|log₂ fold change|≥1且P≤0.05的标准筛选Am7 vs Am10 比较组的DEpiRNA。通过相关软件预测DEpiRNA的靶mRNA并进行GO和KEGG数据库注释。利用Cytoscape软件对DEpiRNA-mRNA调控网络进行可视化。通过Stem-loop RT-PCR对随机选取的Am7与Am10两组共有的6个piRNA的表达进行验证。利用RT-qPCR对随机挑选的6个DEpiRNA的表达趋势进行验证。【结果】从意大利蜜蜂工蜂中肠中共鉴定到596个piRNA,长度介于24~33 nt,且Am7和Am10组中不同长度分布范围的piRNA的首位碱基偏向性具有明显差异。在Am7 vs Am10比较组共筛选出41个DEpiRNA,其中piR-ame-11093, piR-ame-1111451,piR-ame-190949和piR-ame-932156可分别靶向1 195, 1 018, 4 040和1 063条mRNA。上述靶mRNA可分别注释到45个功能条目和45条通路。Stem-loop RT-PCR验证结果显示6个piRNA(piR-ame-1084826, piR-ame-11093, piR-ame-14476, piR-ame-24995, piR-ame-39500和piR-ame-774987)均真实表达。RT-qPCR结果显示共有6个DEpiRNA (piR-ame-1084826,piR-ame-11093, piR-ame-14476, piR-ame-24995, piR-ame-39500和piR-ame-774987)的表达趋势与测序数据中的表达趋势一致,证实了本研究中sRNA-seq数据的可靠性和piRNA差异表达趋势的真实性。【结论】本研究从意大利蜜蜂工蜂中肠中鉴定到596个piRN,长度介于24~33 nt。不同长度分布范围的piRNA的首位碱基偏向性具有明显差异;piR-ame-11093, piR-ame-1111451,piR-ame-190949和piR-ame-932156具有靶向调控基因表达参与工蜂中肠发育过程的潜力。     

意大利蜜蜂幼虫肠道的miRNAs的生物信息学预测及分析

【目的】蜜蜂肠道是食物消化和营养吸收的主要部位,同时也是抵御病原侵染的重要场所。本研究旨在对意大利蜜蜂Apismelliferaligustica(简称意蜂)幼虫肠道的microRNAs(miRNAs)及其靶基因进行深入分析,进而解析miRNAs在幼虫肠道发育和生长过程中的作用。【方法】利用small RNA-seq(sRNA-seq)技术对意蜂4、5和6日龄幼虫肠道进行测序,通过相关生物信息学软件对意蜂幼虫肠道的miRNAs进行预测及分析。利用TargetFinder软件预测miRNAs的靶基因,然后将靶基因通过BLAST软件进行GO和KEGG数据库的功能注释。利用Cytoscape软件构建miRNAs与其靶向结合的mRNAs的调控网络。【结果】意蜂幼虫肠道样品的测序共获得96 329 456条有效标签序列(Clean tags),预测出560个miRNAs,包括45个novelmiRNAs。上述miRNAs的长度主要分布在17-27nt之间,且不同长度的miRNAs的首位碱基偏向性具有明显差异。表达量聚类分析结果显示,331个miRNAs为4、5和6日龄幼虫肠道所共有且表达量较为稳定,随着发育时间延长,特有miRNAs的表达水平呈总体增加的趋势。利用TargetFinder软件共预测出16479个意蜂幼虫肠道的靶基因,其中有8132和3361个可分别注释到GO和KEGG数据库,进一步分析发现有224个靶基因注释在Wnt信号通路,分别有2个和27个靶基因与保幼激素和蜕皮激素的调控密切相关。【结论】本研究在全基因组水平对意蜂幼虫肠道的miRNAs进行预测和分析,研究结果不仅丰富了意蜂的miRNAs信息,也为深入研究意蜂幼虫肠道的生长发育机理打下了初步基础。     

中华蜜蜂6日龄幼虫响应蜜蜂球囊菌胁迫的环状RNA应答

【目的】蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)专性侵染蜜蜂幼虫而导致白垩病,危害蜜蜂健康和养蜂生产。本研究旨在探究中华蜜蜂(Apis cerana cerana,简称中蜂) 6日龄幼虫响应球囊菌胁迫的环状RNA(circular RNA,circRNA)差异表达谱及差异表达circRNA (differentially expressed circRNA,DEcircRNA)在宿主胁迫应答中的潜在功能。【方法】利用去除线性RNA的circRNA-seq技术对正常和球囊菌侵染的中蜂6日龄幼虫肠道(AcCK和AcT)进行测序。利用find_circ软件鉴定circRNA,统计circRNA的长度和环化类型。根据|log_2(Fold change)|≥1和P≤0.05的标准筛选DEcircRNA。将DEcircRNA的来源基因比对Gene ontology (GO)数据库和Kyoto Encyclopedia of Genes and Genomes (KEGG)数据库,从而获得功能及通路(pathway)注释。随机挑选3个DEcircRNA进行RT-qPCR验证。【结果】AcCK和AcT的circRNA-seq分别得到76342570和68269362条原始读段(raw reads),经严格质控得到74524108和66974392条有效读段(clean reads),Q30分别为92.75%和94%,GC含量分别为54.31%和54.90%。比对上东方蜜蜂(Apis cerana)参考基因组的短序列读段(anchor reads)共计23648400条。AcCK和AcT中分别鉴定到805和702个circRNA,长度均介于201–1000 nt,数量最多的环化类型均为已注释外显子circRNA,但分布在不同长度、不同环化类型的circRNA数量存在差异。AcCK vs AcT比较组共有494个DEcircRNA,包括257个上调circRNA和237个下调circRNA;上调和下调幅度最大的circRNA分别为novel_circ_000123和novel_circ_000726。上述DEcircRNA的来源基因可注释到11条生物学进程相关条目,9条分子功能相关条目,9条细胞组分相关条目,以及73条通路。进一步分析发现,部分DEcircRNA的来源基因注释到7条细胞免疫通路和3条体液免疫通路。【结论】中蜂6日龄幼虫响应球囊菌胁迫的过程中可能通过改变分布在不同长度和环化类型的circRNA数量,以及特异性表达一些circRNA和调节部分circRNA的表达量对病原产生应答;novel_circ_000027、novel_circ_000127、novel_circ_000312等DEcircRNA在宿主的胁迫应答过程中可能通过调控氧化磷酸化、细胞和体液免疫等通路发挥特殊作用。研究结果为深入理解中蜂幼虫对球囊菌的胁迫应答机制及二者的相互作用机制提供了新见解。     

球囊菌侵染对意大利蜜蜂幼虫肠道免疫与解毒相关基因表达及蛋白质、脂质和糖类含量的影响

[目的]球囊菌Ascosphaera apis是一种真菌性病原体,可以侵染意大利蜜蜂Apis mellifera ligustica幼虫,导致蜂群患白垩病.本研究通过球囊菌侵染意大利蜜蜂幼虫,探究对意大利蜜蜂幼虫肠道免疫解毒相关基因表达及蛋白质、脂质和糖类含量的影响,分析球囊菌侵染胁迫下意大利蜜蜂免疫应答反应与营养物质...     

中华蜜蜂6日龄幼虫响应蜜蜂球囊菌侵染的长链非编码RNA应答研究

[目的]本研究旨在探究长链非编码RNA (long non-coding RNA,lncRNA)在中华蜜蜂(Apis cerana cerana,简称中蜂)6日龄幼虫应答蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)侵染过程中的差异表达谱及调控作用.[方法]利用链特异性cDNA建库的RNA-seq技术对未被...     

中华蜜蜂和意大利蜜蜂感染囊状幼虫病毒后免疫及发育相关基因表达比较

【目的】本文旨在探究囊状幼虫病毒(sacbrood virus, SBV)对中华蜜蜂(Apis cerana cerana,简称中蜂)和意大利蜜蜂(Apis mellifera ligustica,简称意蜂)工蜂幼虫发育和免疫基因、营养代谢基因、抗病毒基因、细胞发育及代谢相关基因表达的影响。【方法】从蜂群中移取2日龄的中蜂和意蜂幼虫,在培养箱(34℃, RH 85%)进行人工饲养,3日龄时接种SBV病毒,每天观察记录死亡情况,并通过实时定量PCR (real-time quantitative polymerase chain reaction, RT-qPCR)检测4日龄和7日龄幼虫体内SBV基因相对表达量,及免疫基因(Apidaecin、Abaecin、Hymenoptaecin、Denfensin、Lys-1、Pgrp-lc、Kenny、Domeless)、营养代谢基因(Ilp1、Hex110、Vg)、抗病毒基因(Dis3、Dicer、Ago1)、细胞组成及发育调控基因(Vhdl、Co-1-iv)以及细胞代谢和调控基因(Mta1)的表达水平。【结果】通过分析比较发现,感染相同剂量...     

LncRNA13164通过ace-miR-4968-y调节中华蜜蜂幼虫对蜜蜂球囊菌侵染的免疫应答

【目的】长链非编码RNA (long non-coding RNA,lncRNA)可通过竞争性内源RNA(competing endogenous RNA,ce RNA)等多种方式发挥重要的调控作用,但蜜蜂lncRNA的功能研究至今依然缺失,lncRNA在蜜蜂免疫应答中的作用尚不清楚。本研究旨在揭示中华蜜蜂(Apis cerana cerana,简称中蜂)幼虫响应蜜蜂球囊菌(Ascosphaera apis)侵染的免疫应答中lncRNA的调控功能及作用机制。笔者团队前期通过深度测序和生物信息学分析发现,lncRNA13164与ace-miR-4968存在靶向关系且潜在参与在中蜂幼虫对蜜蜂球囊菌侵染的应答。【方法】利用实时荧光定量PCR (real-time fluorescence quantitative polymerase chain reaction,RT-qPCR)检测蜜蜂球囊菌接种后中蜂幼虫肠道内lncRNA13164的表达谱。采用ILoc-Lnc RNA软件预测lncRNA13164的亚细胞定位。联用RNAhybrid、Miranda和Target Scan软件预测lnc...     

球囊菌胁迫中华蜜蜂幼虫肠道过程中病原的转录组学研究

【目的】本研究利用RNA-seq技术对球囊菌胁迫的中华蜜蜂(中蜂)幼虫肠道进行深度测序,经趋势分析得到差异表达基因(DEGs)的显著表达模式,进而对胁迫过程中的球囊菌进行转录组学分析。【方法】利用Illumina HiSeq 2500平台对球囊菌胁迫的中蜂幼虫肠道进行深度测序,并利用相关软件进行了深入分析。最后,通过RT-qPCR对RNA-seq数据进行了验证。【结果】本研究共得到球囊菌的41133932条高质量clean reads。22865个DEGs共聚类为8个基因表达模式,其中,16769个DEGs聚类为2个显著上调趋势与2个显著下调趋势。GO富集分析结果显示,显著上调与显著下调趋势中的DEGs分别富集于40与37个GO term,基因富集数最多的为细胞进程(2486 unigenes)。KEGG代谢通路(pathway)富集分析结果显示,显著上调与显著下调趋势中的DEGs分别富集于119和112个pathway,基因富集数最多的分别是氨基酸生物合成(127 unigenes)与核糖体(98 unigenes)。进一步分析表明球囊菌在胁迫中蜂幼虫肠道的过程中通过提高物质合成促进其增殖,而宿主通过抑制球囊菌的蛋白合成抵御病原入侵。富集在MAPK信号通路的11个DEGs的表达水平随着胁迫时间的延长而逐渐下降,推测中蜂幼虫通过抑制该通路而阻遏球囊菌增殖。【结论】本研究不仅为揭示白垩病过程中的球囊菌-中蜂幼虫互作提供了重要信息,也为阐明不同抗性蜂种的球囊菌抗性差异奠定了基础。     

侵染中华蜜蜂6日龄幼虫的蜜蜂球囊菌的微小RNA差异表达谱及调控网络

【目的】蜜蜂球囊菌(Ascosphaeraapis,简称球囊菌)是专性侵染蜜蜂幼虫的致死性真菌病原。MicroRNA(miRNA)作为一类重要的基因表达调控因子,能够广泛参与真菌及其宿主的相互作用过程。本研究通过比较分析球囊菌孢子(AaCK)和侵染中华蜜蜂(Apis cerana cerana,简称中蜂) 6日龄幼虫肠道内的球囊菌(AaT)的smallRNA(sRNA)组学数据对球囊菌的差异表达miRNA(differentiallyexpressed miRNA,DEmiRNA)、靶mRNA及二者间的调控网络进行全面解析,旨在揭示miRNA介导的球囊菌对中蜂幼虫的侵染机制。【方法】对于球囊菌侵染的中蜂6日龄幼虫肠道的small RNA-seq (sRNA-seq)数据,利用BLAST工具连续比对东方蜜蜂(Apiscerana)和球囊菌的参考基因组筛滤得到AaT的sRNA组学数据。分别将AaCK和AaT的sRNA组学数据比对miRBase数据库,对球囊菌侵染宿主前后miRNA的数量和结构特征进行分析。联用RNAhybrid+svm_light、Miranda和Tar...     

东方蜜蜂微孢子虫侵染意大利蜜蜂工蜂过程的nce-miR-23928及其靶基因表达谱

【背景】东方蜜蜂微孢子虫(Nosema ceranae)专性侵染成年蜜蜂中肠上皮细胞而导致的微孢子虫病给养蜂业造成严重损失。【目的】检测东方蜜蜂微孢子虫nce-miR-23928及其靶基因在侵染意大利蜜蜂(Apis mellifera ligustica)工蜂过程的表达谱,为深入探究nce-miR-23928在东方蜜蜂微孢子虫侵染中的功能及调控机制提供依据。【方法】通过RNAhybrid、miRanda和TargetScan软件预测nce-miR-23928的靶基因。使用BLAST工具将上述靶基因比对到基因本体论(geneontology,GO)、京都基因和基因组百科全书(Kyoto encyclopedia of genes and genomes, KEGG)、Nr和Swiss-Prot数据库以获得相应注释。采用实时荧光定量PCR(realtimequantitativePCR,RT-qPCR)技术检测nce-miR-23928及其靶基因在东方蜜蜂微孢子虫侵染意蜂工蜂过程中的相对表达量。【结果】相较于接种后1 d (1 day post infection, 1 dpi),nce-...     

基于small RNA组学分析揭示意大利蜜蜂响应东方蜜蜂微孢子虫胁迫的免疫应答机制

【目的】通过深度测序和组学分析在small RNA组学层面揭示意大利蜜蜂(Apis mellifera ligustica,简称意蜂)响应东方蜜蜂微孢子虫(Nosema ceranae)的免疫应答机制。【方法】利用small RNA-seq技术对正常及N. ceranae胁迫7 d和10 d的意蜂工蜂中肠(Am7CK、Am7T、Am10CK和Am10T)进行深度测序。利用相关生物信息学软件对测序数据进行质控、已知micro RNA (mi RNA)鉴定、新mi RNA预测及mi RNA的结构特征分析。通过Stem-loop RT-PCR验证新mi RNA的表达。按照|log2(Fold change)|≥1和P≤0.05的标准筛选出Am7CKvs.Am7T和Am10CKvs.Am10T比较组的差异表达mi RNA(differentially expressed mi RNA,DEmi RNA)。利用软件预测DEmi RNA靶向结合的m RNA并进行GO和KEGG数据库注释,根据注释信息对细胞和体液免疫相关通路及富集靶m RNA进行统计和分析。根据靶向结合关系构建DEmi RNA和免疫通路相关差异表达m RNA (differentially expressed m RNA,DEm RNA)的调控网络。采用RT-q PCR对数据的可靠性和DEmi RNA的差异表达进行验证。【结果】共获得165895574条原始读段和132028990条有效序列标签,各组的组内Pearson相关性平均在87.92%及以上。共鉴定到928个已知mi RNA和56个新mi RNA。这些mi RNA的长度介于18–28nt,多数的长度为18nt和22nt且首位碱基主要偏向U。验证了12个新mi RNA的真实表达。Am7CKvs.Am7T比较组包含48个上调mi RNA和36个下调mi RNA;Am10CK vs. Am10T比较组包含56个上调mi RNA和51个下调mi RNA。两个比较组的DEmi RNA可分别靶向结合9827个和10720个m RNA。这些靶m RNA可分别注释到50和47条功能条目,以及138和135条KEGG通路。DEmi RNA与免疫通路相关靶m RNA的调控网络分析结果显示,Am7CK vs. Am7T中有26个DEmi RNA靶向与内吞作用等免疫通路相关的10个DEm RNA;Am10CK vs. Am10T中有15个DEmi RNA靶向与MAPK信号通路等免疫通路相关的10个DEm RNA。验证了测序数据和4个DEmi RNA差异表达的可靠性。【结论】研究结果揭示了宿主DEmi RNA可能通过调控物质和能量代谢、细胞和体液免疫对N. ceranae产生应答,但DEmi RNA不参与抗菌肽基因的表达调控;mi R-1-z可能参与宿主的细胞增殖、细胞凋亡和免疫进程;氧化磷酸化通路可能在宿主免疫应答及宿主-病原互作中发挥特殊作用。     

中华蜜蜂幼虫肠道响应球囊菌胁迫的microRNA应答分析

【目的】蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)是一种能够侵染中华蜜蜂(Apis cerana cerana,简称中蜂)幼虫的致死性真菌病原。微小RNA(microRNA,miRNA)可通过在转录后水平靶向抑制或降解mRNA而参与宿主与病原互作过程。本研究旨在对球囊菌胁迫的中蜂6日龄幼虫肠道的差异表达miRNA(DEmiRNA)及其靶基因进行深入分析,进而揭示DEmiRNA在中蜂响应球囊菌胁迫应答过程中的作用。【方法】利用Illumina MiSeq平台对正常及球囊菌胁迫的中蜂6日龄幼虫肠道(AcCK和AcT)进行测序,通过相关生物信息学软件预测DEmiRNA及其靶基因。通过Blast将靶基因注释到GO和KEGG数据库。利用Cytoscape软件构建DEmiRNA与其靶mRNA的调控网络。通过Stem-loop RT-PCR和qPCR验证测序数据的可靠性。【结果】本研究共预测出537个miRNA,其长度分布介于16–35 nt之间,且不同长度的miRNA首位碱基偏向性差异明显。通过Stem-loop RT-PCR证实了10个novel miRNA的表达。AcCK vs AcT比较组共有54个DEmiRNA,包含31个上调和23个下调miRNA,可分别靶向结合6170和8199个靶基因。GO分类结果显示上调和下调miRNA的靶基因分别涉及47和47个条目,富集基因数最多的皆为结合细胞进程和催化活性。KEGG代谢通路(pathway)富集分析结果表明上调和下调miRNA的靶基因分别富集在134和126条pathway,富集基因数最多的均为内吞作用和内质网中的蛋白质加工。调控网络分析结果表明,DEmiRNA及其靶mRNA形成十分复杂的调控关系;31个DEmiRNA可靶向结合51个与泛素介导的蛋白水解相关的mRNA,18个DEmiRNA可靶向结合14个与Jak-STAT信号通路相关的mRNA;miR-1277-x、miR-26-x、miR-27-y、miR-30-x、miR-6052-x等16个miRNA共同参与了上述两条免疫通路的调控。最后,随机挑选3个DEmiRNA进行qPCR验证,结果证明了测序数据的可靠性。【结论】本研究提供了中蜂幼虫肠道在球囊菌胁迫后期的miRNA的表达谱和差异表达信息,揭示了球囊菌与宿主之间在miRNA组学水平存在复杂的互作。miR-6052-x和miR-1277-x作为调控网络的核心可能通过影响细胞凋亡参与宿主的免疫防御,miR-26-x和miR-30-x可能通过调控Jak-STAT信号通路参与宿主的胁迫应答。本研究筛选出的关键DEmiRNA有望作为治疗白垩病的分子靶标。