摘要  使用对家蚕Bombyx mori具有高、低毒力的球孢白僵菌Beauveria bassiana菌株侵染家蚕，分析家蚕应对球孢白僵菌浸染的免疫应答基因。采用新一代Illumina高通量测序技术对接种高、低毒力球孢白僵菌家蚕样品与未接种的家蚕样品分别进行转录组测序，筛选出差异表达基因，利用生物信息学软件分析差异表达基因的功能、分类及涉及的信号通路等。结果表明，接种高毒力菌株后诱导924个家蚕基因发生差异表达（FDR≤0.001，log₂ Ratio≥1），包括252个上调基因和672个下调基因；KEGG分析表明，557个差异表达基因被富集到31个功能通路。低毒力菌株侵染家蚕后诱导2 697个基因显著差异表达（FDR≤0.001，log₂ Ratio≥1），包括778个上调基因和1 919个下调基因；KEGG分析表明，1 545个差异表达基因被富集到48个功能通路。分析发现，差异表达基因编码蛋白包括RNA聚合酶、丝氨酸蛋白酶、丝氨酸蛋白酶抑制剂以及ECM受体作用蛋白免疫相关蛋白等。另外，膜结合碱性磷酸酶样蛋白、E4酯酶、肽聚糖识别蛋白、黑芥子酶等的编码基因在球孢白僵菌侵染过程中也高表达，推测这些基因可能参与家蚕应对球孢白僵菌侵染的免疫反应。研究结果为进一步研究家蚕和其他昆虫相关免疫基因在抵抗病原物侵染中的功能奠定了基础。