杜 宇，熊翠玲，史秀丽，郑燕珍，付中民，徐细建，陈大福，郭 睿，*

(1.福建农林大学 蜂学学院，福建 福州 350002； 2.新疆维吾尔自治区蜂业发展中心，新疆 乌鲁木齐 830000)

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

杜 宇1，熊翠玲1，史秀丽2，郑燕珍1，付中民1，徐细建1，陈大福1，郭 睿1，*

(1.福建农林大学 蜂学学院，福建 福州 350002； 2.新疆维吾尔自治区蜂业发展中心，新疆 乌鲁木齐 830000)

为检测蜜蜂球囊菌(Ascosphaeraapis)在胁迫意大利蜜蜂(Apismelliferaligustica，简称意蜂)幼虫后期的基因表达谱，利用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个基因上调表达，说明该通路在球囊菌胁迫后期被激活。研究结果不仅为揭示球囊菌在胁迫意蜂幼虫后期的病原-宿主互作提供了重要信息，也为阐明球囊菌致病的分子机制奠定了基础。

球囊菌；意大利蜜蜂；幼虫肠道；差异表达基因；RNA-seq

蜜蜂是世界范围内最重要的授粉昆虫，其在经济创收和生态保护等方面均发挥着重要作用[1]。……