意大利蜜蜂6日龄幼虫肠道内球囊菌的差异表达基因分析

摘要

The purified Ascosphaera apis spores and the 6-day-old larval gut of Apis mellifera ligustica under the stress of A.apis were sequenced using RNA-seq technology to detect gene expression patterns of the fungal pathogen.A total of 21 361 differentially expressed genes (DEGs) were obtained, including 15 306 up-regulated genes and 6 055 down-regulated genes.GO enrichment analysis showed that the up-regulated genes were enriched in 39 terms, among them the mostly enriched ones were cellular process (2 416 unigenes), cell (2 408 unigenes) and cell part (2 408 unigenes);the down-regulated genes were enriched in 38 terms, and the mostly enriched ones were metabolic process (1 227 unigenes), cellular process (1 185 unigenes) and cell (1 131 unigenes).KEGG enrichment analysis displayed that the up-regulated genes were enriched in 119 pathways, and the largest groups were ribosome (221 unigenes), protein processing in endoplasmic reticulum (190 unigenes) and endocytosis (177 unigenes);the down-regulated genes were enriched in 113 pathways, and the mostly enriched ones were ribosome (144 unigenes), RNA transport (113 unigenes) as well as biosynthesis of amino acid (108 unigenes).Furthermore, 64 up-regulated genes were found to be enriched in MAPK signaling pathway, suggesting that it was induced to a large extent.The data not only provided key information for uncovering the pathogen-host interaction during the late stage of the stress of A.apis on A.m.ligustica larvae, but also laid some foundation for clarifying molecular mechanisms regulating the pathogenesis of A.apis.%为检测蜜蜂球囊菌(Ascosphaera apis)在胁迫意大利蜜蜂(Apis mellifera ligustica,简称意蜂)幼虫后期的基因表达谱,利用RNA-seq技术对纯培养的球囊菌孢子及球囊菌胁迫的意蜂6日龄幼虫肠道进行深度测序.通过比较处理组和对照组得到21 361个差异表达基因(DEGs),包括15 306个上调基因和6 055个下调基因.GO富集分析结果显示,上调基因富集于39个GO条目(term),基因富集数最多的为细胞进程(2 416 unigenes)、细胞(2 408 unigenes)和细胞组件(2 408 unigenes);下调基因富集于38个GO term,基因富集数最多的为代谢进程(1 227 unigenes)、细胞进程(1 185 unigenes)和细胞(1 131 unigenes).KEGG代谢通路富集分析结果显示,上调基因富集在119个通路,其中基因富集数最多的是核糖体(221 unigenes),其次为内质网中蛋白加工(190 unigenes)和内吞作用(177 unigenes);下调基因富集在113个通路上,基因富集数最多的是核糖体(144 unigenes),其次为RNA转运(113 unigenes)和氨基酸生物合成(108 unigenes).进一步分析发现,富集在MAPK信号通路上的64个基因上调表达,说明该通路在球囊菌胁迫后期被激活.研究结果不仅为揭示球囊菌在胁迫意蜂幼虫后期的病原-宿主互作提供了重要信息,也为阐明球囊菌致病的分子机制奠定了基础.