[目的]本研究旨在通过差异表达基因(differentially expressed gene,DEG)分析以及毒力因子和其他侵染相关因子分析,在转录组水平揭示东方蜜蜂微孢子虫Nosema ceranae侵染意大利蜜蜂Apis mellifera ligustica的分子机制.[方法]基于前期已获得高质量的东方蜜蜂微孢子虫纯化孢子(NcCK)及侵染意大利蜜蜂工蜂7和10 d的东方蜜蜂微孢子虫(分别为NcT1和NcT2)转录组数据,根据P≤0.05且|log2(Fold change)|≥1的标准,通过比较分析筛选出NcCK vs NcT1,NcCK vsNcT2和NcT1 vs NcT2比较组的DEG.通过相关生物信息学软件对上述DEG进行Venn分析、GO分类和KEGG代谢通路富集分析.根据Nr和KEGG数据库注释信息和相关文献进行对东方蜜蜂微孢子虫的毒力因子和侵染相关因子的统计和分析.通过RT-qPCR验证转录组数据及DEG表达趋势.[结果]从NcCK vs NcT1,NcCK vs NcT2和NcT1 vs NcT2比较组分别鉴定出1 397,1 497和52个DEG.Venn分析结果显示各比较组共有的上调和下调基因分别为10和1个.GO分类结果显示,NcCK vs NcT1和NcCK vs NcT2中DEG富集数最多的功能条目为代谢进程、细胞进程、单组织进程、细胞、细胞组件、细胞器、催化活性和结合,而NcT1 vs NcT2中DEG富集数最多的是代谢进程、细胞进程、单组织进程、催化活性和结合.KEGG代谢通路

作者：耿四海；周丁丁；范小雪；蒋海宾；祝智威；王杰；范元婵；王心蕊；熊翠玲；郑燕珍；付中民；陈大福；郭睿

来源：昆虫学报 2020 年 63卷 3期