| 摘要: |
| 为探究中华蜜蜂（Apis cerana cerana，以下简称中蜂）与意大利蜜蜂（Apis mellifera ligustica，以下简称意蜂）幼虫的球囊菌（Ascospahaera apis）抗性差异产生的原因，利用Venn分析、GO分类分析、KEGG代谢通路分析以及Ka/Ks分析对二者的转录组进行比较研究。Venn分析结果显示，中蜂与意蜂的同源基因有17 656个，归属于8 111个基因家族，二者的特有基因分别有1 158和241个，分别归属于468和86个基因家族。GO分类结果显示，中蜂和意蜂的特有基因分别富集在38和28个GO条目。KEGG代谢通路富集分析结果显示，中蜂和意蜂的特有基因分别富集于96和21个代谢通路；进一步分析发现中蜂的特有基因富集在内吞作用、溶酶体、泛素介导的蛋白水解和黑化作用等4个细胞免疫通路，以及MAPK信号通路和Toll-like受体信号通路等2个体液免疫通路，而意蜂仅有1个特有基因富集在MAPK信号通路。Ka/Ks分析结果显示受到强烈正向选择、弱正向选择、中性选择和负向选择的单拷贝同源基因分别有37、281、2 630和3 269个。对受到强烈正向选择的基因进行GO分类和KEGG代谢通路富集分析，GO结果显示中蜂和意蜂的单拷贝同源基因富集的GO条目类型相同；KEGG结果显示中蜂的单拷贝基因仅富集在氧化磷酸化。上述结果表明免疫相关基因数量的差异是二者的球囊菌抗性差异产生的重要原因之一，中蜂在与球囊菌的协同进化过程中可能通过提高能量利用从而限制球囊菌的增殖。本研究结果可为阐明中蜂及意蜂幼虫的球囊菌抗性差异产生的分子机制提供参考。 |
| 关键词: 中华蜜蜂 意大利蜜蜂 幼虫 比较转录组学 球囊菌 |
| DOI：10.11841/j.issn.1007-4333.2019.05.13 |
| 投稿时间：2018-09-20 |
| 基金项目:国家自然科学基金项目（31702190）；现代农业产业技术体系建设专项资金（CARS-44-KXJ7）；福建省教育厅中青年教师教育科研项目（JAT170158）；福建农林大学科技创新专项基金（CXZX2017343）；福建省大学生创新创业训练计划项目（201610389053） |
|
| Unraveling the mechanism regulating the Ascosphaera apis-resistance difference between Apis cerana cerana and Apis mellifera ligustica larvae based on comparative transcriptome analysis |
|
XIONG Cuiling, DU Yu, WANG Hongquan, ZHENG Yanzhen, FU Zhongmin, WANG Haipeng, ZHANG Lu, CHEN Dafu, GUO Rui
|
| (College of Bee Science, Fujian Agriculture and Forestry University, Fuzhou 350002, China) |
| Abstract: |
| To explore the factors regulating the Ascosphaera apis-resistance difference between Apis cerana cerana and Apis mellifera ligustica larvae,Venn analysis,GO classification analysis and KEGG pathway enrichment analysis and Ka/Ks analysis were adopted in this study.The Venn analysis showed there were 17 656 homologous genes belonging to 8 111 genes families between A.c.cerana and A.m.ligustica,and 1 158 and 241 unique genes belonging to 468 and 86 gene families,respectively.The GO classification suggested these unique genes were involved in 38 and 28 GO terms,respectively.The KEGG pathway enrichment analysis demonstrated the unique genes were enriched in 96 and 21 pathways.Further analysis showed the unique genes of A.c.cerana were enriched in 4 cellular immune pathways including endocytosis,and 2 humoral immune pathways including MAPK signaling pathway,while only one unique gene of A.m.ligustica was enriched in MAPK signaling pathway.The Ka/Ks analysis indicated there were respectively 37,281,2 630 and 3 269 homologous genes with a single copy under strong positive selection,weak positive selection,neutral selection and negative selection.GO classification and KEGG pathway enrichment analysis were further performed for homologous gene with a single copy under strong positive selection,and the result suggested that the enriched GO terms were the same between A.c.cerana and A.m.ligustica,and only one single-copy homologous gene of the latter was enriched in oxidative phosphorylation.The above results demonstrated that the difference in the quantity of immunity-related genes might be a significant factor leading to the A.apis-resistance difference between A.c.cerana and A.m.ligustica larvae.A.c.cerana might enhance the utilization of energy to limit the proliferation of A.apis.This study would provide basice foundation to clarify the molecular mechanism underlying the A.apis-resistance difference between A.c.cerana and A.m.ligustica larvae. |
| Key words: Apis cerana cerana Apis mellifera ligustica larvae comparative transcriptome analysis Ascospaera apis |
