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

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

作　　者：杜宇 蒋海宾 王杰 范小雪 王秀娜[3,4] 冯睿蓉 张文德 隆琦 熊翠玲 郑燕珍 陈大福 郭睿 DU Yu;JIANG Hai-Bin;WANG Jie;FAN Xiao-Xue;WANG Xiu-Na;FENG Rui-Rong;ZHANG Wen-De;LONG Qi;XIONG Cui-Ling;ZHENG Yan-Zhen;CHEN Da-Fu;GUO Rui(College of Animal Sciences(College of Bee Science),Fujian Agriculture and Forestry University,Fuzhou 350002,China;Engineering Research Center of Processing and Application of Bee Products of Ministry of Education,Fujian Agriculture and Forestry University,Fuzhou 350002,China;College of Life Sciences,Fujian Agriculture and Forestry University,Fuzhou 350002,China;Key Laboratory of Pathogenic Fungi and Mycotoxins of Fujian Province,Fujian Agriculture and Forestry University,Fuzhou 350002,China)

机构地区：[1]福建农林大学动物科学学院(蜂学学院),福州350002 [2]福建农林大学,蜂产品加工与应用教育部工程研究中心,福州350002 [3]福建农林大学生命科学学院,福州350002 [4]福建农林大学,福建省病原真菌与真菌毒素重点实验室,福州350002

出　　处：《昆虫学报》2021年第3期363-373,共11页Acta Entomologica Sinica

基　　金：国家现代农业产业技术体系建设专项资金项目(CARS-44-KXJ7);福建省自然科学基金项目(2018J05042);福建省教育厅中青年教师教育科研项目(JAT170158);福建农林大学杰出青年科研人才计划项目(xjq201814);福建省病原真菌与真菌毒素重点实验室(福建农林大学)开放课题;福建农林大学优秀硕士学位论文资助基金(杜宇)。

摘　　要：【目的】本研究旨在探究蜜蜂球囊菌Ascosphaera apis(简称“球囊菌”)菌丝(Aam)和孢子(Aas)中的全长转录本与菌丝生长、孢子萌发以及有性生殖的相关性。【方法】将前期获得的Aam和Aas的Nanopore长读段测序数据中的有效读段与球囊菌参考基因组注释的已知转录本进行序列比对,获取全长转录本与已知转录本的对应信息。利用生物信息学方法对Aam和Aas的全长转录本进行差异表达、功能注释以及结构特征分析。【结果】Aam和Aas共有的全长转录本数量为19966条,特有的全长转录本数量分别为1273和2856条。Aas vs Aam比较组含有3230条差异表达转录本(differentially expressed transcripts,DETs),包含3072条上调和158条下调。GO功能注释结果显示,这些DETs涉及代谢进程、细胞、催化活性等GO条目。KEGG注释结果显示,这些DETs注释到内吞作用、MAPK信号通路、糖酵解/糖异生、碳代谢以及氨基酸的生物合成等相关通路。进一步分析发现,部分全长转录本的剪接异构体在Aam和Aas中具有不同的表达量和结构。【结论】本研究发现球囊菌在菌丝和孢子两个不同阶段伴随着转录本表达量和结构的变化,研究结果为深入探究不同剪接异构体在球囊菌的菌丝生长、孢子萌发和有性生殖中的分子功能提供了理论依据和数据基础。【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 19966 full-length transcripts in both Aam and Aas.The specific full-length transcripts of Aam and Aas were 1273 and 2856,respectively.In the Aas vs Aam comparison group,3230 differentially expressed transcripts(DETs)including 3072 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.

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

分 类 号：Q895.1[生物学]