浙江农业学报 ›› 2017, Vol. 29 ›› Issue (7): 1119-1128.DOI: 10.3969/j.issn.1004-1524.2017.07.09

• 动物科学 • 上一篇    下一篇

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

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

  1. 1.福建农林大学 蜂学学院,福建 福州 350002;
    2.新疆维吾尔自治区蜂业发展中心,新疆 乌鲁木齐 830000
  • 收稿日期:2016-03-08 出版日期:2017-07-20 发布日期:2017-07-24
  • 通讯作者: 郭睿,E-mail:ruiguo@fafu.edu.cn
  • 作者简介:杜宇(1994—),男,河北唐山人,本科生,专业为蜂学,E-mail:599195720@qq.com
  • 基金资助:
    现代农业产业技术体系建设专项资金(CARS-45-KXJ7); 福建农林大学科技发展资金(KF2015123); 国家级大学生创新创业项目(201610389016)

Transcriptome analysis of differentially expressed genes in Ascosphaera apis stressing the 6-day-old larval gut of Apis mellifera ligustica

DU Yu1, XIONG Cuiling1, SHI Xiuli2, ZHENG Yanzhen1, FU Zhongmin1, XU Xijian1, CHEN Dafu1, GUO Rui1, *   

  1. 1. College of Bee Science, Fujian Agriculture and Forestry University, Fuzhou 350002, China;
    2. Apicultural Development Centre in Xinjiang Uygur Autonomous Region, Urumqi 830000, China
  • Received:2016-03-08 Online:2017-07-20 Published:2017-07-24

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

关键词: 球囊菌, 意大利蜜蜂, 幼虫肠道, 差异表达基因, RNA-seq

Abstract: 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.

Key words: Ascosphaera apis, Apis mellifera ligustica, larval gut, DEGs, RNA-seq

中图分类号: