昆虫学报 ›› 2021, Vol. 64 ›› Issue (3): 363-373.doi: 10.16380/j.kcxb.2021.03.008

• 研究论文 • 上一篇    下一篇

蜜蜂球囊菌菌丝和孢子中全长转录本的差异表达分析

杜宇1,#, 蒋海宾1,#, 王杰1, 范小雪1, 王秀娜3,4, 冯睿蓉1, 张文德1, 隆琦1, 熊翠玲1, 郑燕珍1,2, 陈大福1,2, 郭睿1,2,*   

  1. (1. 福建农林大学动物科学学院(蜂学学院), 福州 350002; 2. 福建农林大学, 蜂产品加工与应用教育部工程研究中心, 福州 350002; 3. 福建农林大学生命科学学院, 福州 350002; 4. 福建农林大学, 福建省病原真菌与真菌毒素重点实验室, 福州 350002)
  • 出版日期:2021-03-20 发布日期:2021-04-20

Differential expression analysis of full-length transcripts in mycelia and spores of Ascosphaera apis 

DU Yu1,#, JIANG Hai-Bin1,#, WANG Jie1, FAN Xiao-Xue1, WANG Xiu-Na3,4, FENG Rui-Rong1, ZHANG Wen-De1, LONG Qi1, XIONG Cui-Ling1, ZHENG Yan-Zhen1, CHEN Da-Fu1,2, GUO Rui1,2,*   

  1.  (1. College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, Fuzhou 350002, China; 2. Engineering Research Center of Processing and Application of Bee Products of Ministry of Education, Fujian Agriculture and Forestry University, Fuzhou 350002, China; 3. College of Life Sciences, Fujian Agriculture and Forestry University, Fuzhou 350002, China; 4. Key Laboratory of Pathogenic Fungi and Mycotoxins of Fujian Province, Fujian Agriculture and Forestry University, Fuzhou 350002, China)
  • Online:2021-03-20 Published:2021-04-20

摘要:

【目的】本研究旨在探究蜜蜂球囊菌Ascosphaera apis(简称“球囊菌”)菌丝(Aam)和孢子(Aas)中的全长转录本与菌丝生长、孢子萌发以及有性生殖的相关性。【方法】将前期获得的Aam 和Aas 的Nanopore长读段测序数据中的有效读段与球囊菌参考基因组注释的已知转录本进行序列比对,获取全长转录本与已知转录本的对应信息。利用生物信息学方法对Aam和Aas的全长转录本进行差异表达、功能注释以及结构特征分析。【结果】Aam和Aas共有的全长转录本数量为19 966条,特有的全长转录本数量分别为1 273和2 856条。Aas vs Aam比较组含有3 230条差异表达转录本(differentially expressedtranscripts, DETs),包含3 072条上调和158条下调。GO功能注释结果显示,这些DETs涉及代谢进程、细胞、催化活性等GO条目。KEGG注释结果显示,这些DETs注释到内吞作用、MAPK信号通路、糖酵解/糖异生、碳代谢以及氨基酸的生物合成等相关通路。进一步分析发现,部分全长转录本的剪接异构体在Aam和Aas中具有不同的表达量和结构。【结论】本研究发现球囊菌在菌丝和孢子两个不同阶段伴随着转录本表达量和结构的变化,研究结果为深入探究不同剪接异构体在球囊菌的菌丝生长、孢子萌发和有性生殖中的分子功能提供了理论依据和数据基础。

关键词: 蜜蜂球囊菌, 菌丝, 孢子, 差异表达转录本, 全长转录本, 三代测序

Abstract:

【Aim】 This study aims to investigate the correlation of the full-length transcripts in Ascosphaera apis mycelium (Aam) and spore (Aas) with the mycelium growth, spore germination and sexual reproduction. 【Methods】 The congruent information about the full-length transcripts and known transcripts was obtained by aligning the clean reads generated from the previously gained Nanopore long-read sequencing data of Aam and Aas to the known transcripts annotated in the reference genome of A. apis. Bioinformatics methods were used to analyze the differential expression pattern, functional annotation and structural characteristics of the full-length transcripts in Aam and Aas. 【Results】There were 19 966 full-length transcripts in both Aam and Aas. The specific full-length transcripts of Aam and Aas were 1 273 and 2 856, respectively. In the Aas vs Aam comparison group, 3 230 differentially expressed transcripts (DETs) including 3 072 up-regulated and 158 down-regulated transcripts were identified. These DETs involve such GO items as metabolic process, cell, catalytic activity, etc. KEGG pathway annotation analysis demonstrated that these DETs are involved in pathways including endocytosis, MAPK signaling pathway, glycolysis/gluconeogenesis, carbon metabolism, and amino acid biosynthesis. Further investigation indicated that the splicing isoforms of some full-length transcripts had different expression levels and structure in Aam and Aas. 【Conclusion】 In this study it was discovered that the expression levels and structure of transcripts in A. apis vary in the two stages of mycelium and spore. These results provide a theoretical and data foundation for deeply investigating the molecular function of various splicing isoforms involved in the mycelium growth, spore germination and sexual reproduction of A. apis.

Key words: Ascosphaera apis, mycelium, spore, differentially expressed transcript, full-length transcript, third-generation sequencing