摘要：【目的】通过对意大利蜜蜂(Apis mellifera ligustica,简称意蜂)工蜂中肠响应东方蜜蜂微孢子虫(Nosema ceranae)胁迫的差异表达基因(differentially expressed gene,DEG)及免疫通路进行深入细致的分析,揭示宿主的细胞和体液免疫应答,为关键免疫应答基因的筛选及功能研究打下基础。【方法】基于前期获得的正常及东方蜜蜂微孢子虫胁迫的意蜂7日龄和10日龄工蜂中肠(Am7CK、Am7T、Am10CK、Am10T)转录组数据,按照FDR≤1,P≤0.05和|log2 fold change|≥1的标准筛选出各组的DEG,利用相关生物信息学软件对DEG进行Pearson相关性分析、Venn分析、GO分类和KEGG代谢通路富集分析,进而对免疫通路富集基因进行统计和分析,利用实时荧光定量PCR(real-time quantitative PCR,RT-qPCR)验证转录组数据的可靠性。【结果】差异分析结果显示,Am7CK vs Am7T比较组包含472个上调基因和385个下调基因;Am10CK vs Am10T比较组包含611个上调基因和360个下调基因。Venn分析结果显示,两个比较组特有的DEG分别为739和853个,共有的DEG为118个。GO分类结果显示,Am7CK vs Am7T比较组中上调和下调基因分别涉及23和29个功能条目,其中富集上调基因数最多的前5位分别为结合、催化活性、代谢进程、细胞进程和单组织进程;富集下调基因数最多的前5位分别为代谢进程、单组织进程、催化活性、细胞进程和结合。Am10CK vs Am10T比较组中上调和下调基因分别涉及36和26个功能条目,其中富集上调基因数最多的前5位分别为单组织进程、结合、细胞进程、催化活性和代谢活性;富集下调基因数最多的前5位分别为结合、细胞进程、催化活性、代谢进程和单组织进程。KEGG代谢通路富集分析结果显示,Am7CK vs Am7T比较组中上调和下调基因分别富集在38和33条代谢通路,富集上调基因数最多的前5条分别为胆汁分泌、内质网蛋白加工、泛素介导的蛋白水解、PI3K-Akt�