Comparative Transcriptome Analysis of Bombyx mori (Lepidoptera) Larval Hemolymph in Response to Autographa californica Nucleopolyhedrovirus in Differentially Resistant Strains

Abstract

Bombyx mori nucleopolyhedrovirus (BmNPV) is a kind of pathogen that causes huge economic losses to silkworm production. Although Autographa californica nucleopolyhedrovirus (AcMNPV) and BmNPV are both baculoviruses, the host domains of these two viruses have almost no intersection in nature. Recently, it has been found that some silkworms could be infected by recombinant AcMNPV through a puncture, which provided valuable material for studying the infection mechanism of baculovirus to silkworm. In this study, comparative transcriptomics was used to analyse the hemolymph of two differentially resistant strains following AcMNPV inoculation. There were 678 DEGs in p50 and 515 DEGs in C108 following viral infection. Among them, the upregulation and downregulation of DEGs were similar in p50; however, the upregulated DEGs were nearly twice as numerous as the downregulated DEGs in C108. The DEGs in different resistant strains differed by GO enrichment. Based on KEGG enrichment, DEGs were mainly enriched in metabolic pathways in p50 and the apoptosis pathway in C108. Moreover, 13 genes involved in metabolic pathways and 11 genes involved in the apoptosis pathway were analysed. Among the DEGs involved in apoptosis, the function of BmTex261 in viral infection was analysed. The BmTex261 showed the highest expression in hemolymph and a significant response to viral infection in the hemolymph of C108, indicating that it is involved in anti-AcMNPV infection. This was further validated by the significantly decreased expression of viral gene lef3 after overexpression of BmTex261 in BmN cells. The results provide a theoretical reference for the molecular mechanism of resistance to BmNPV in silkworms.