筛选柑橘木虱应对球孢白僵菌侵染的免疫应答及其网络调控基因,以进一步探讨柑橘木虱对球孢白僵菌的免疫防御机制。采用Illumina高通量测序平台对感染球孢白僵菌24h、48h、72h与健康柑橘木虱转录组进行了测序,利用生物信息学软件对筛选到的差异表达基因进行了功能注释、分类以及参与的信号通路分析。组装得到138 313条不可延长的非冗余Unigene,其N50和N90分别为2 532bp和413bp,平均长度为1 191.26bp。在CK vs.S24h、CK vs.S48h和CK vs.S72h三个转录组里分别获得了405、614和542个显著差异表达基因,GO富集分析表明56个、83个和60个GO terms分别有显著性富集现象,KEGG pathway分析结果显示分别有98、333和247个差异表达基因显著富集到10、27和25条代谢通路。差异表达基因主要富集在能量代谢、离子运输、转录和翻译调控、生殖和发育调控以及免疫防御反应等相关通路,大部分编码潜在的与免疫识别及调控相关的基因,并从中筛选到5个显著上调表达的免疫相关基因,为开展柑橘木虱免疫应答虫生真菌的研究奠定理论基础。