Identification of Key Genes Involved in Resistance to Early Stage of BmNPV Infection in Silkworms

Study on anti-BmNPV mechanism of branched-chain amino acid aminotransferases in silkworm

Bombyx mori nucleopolyhedrovirus (BmNPV) is the most important virus that threatens sericulture industry. At present, there is no effective treatment for BmNPV infection in silkworms, and lncRNA plays an important role in biological immune response and host-virus interaction, but there are relatively few studies in silkworms. In this study, the four midgut tissue samples of the resistance strain NB (NB) and susceptible strain 306 (306) and the NB and 306 continuously infected with BmNPV for 96 h are used for whole transcriptome sequencing to analyze the differences in the genetic background of NB and 306 and the differences after inoculation of BmNPV, and the significantly different mRNA, miRNA and lnRNA between NB and 306 after BmNPV inoculation were screened. By comparing NB and 306, 2651 significantly different mRNAs, 57 significantly different miRNAs and 198 significantly different lncRNAs were screened. By comparing NB and 306 after BmNPV inoculation, 2684 significantly different mRNAs, 39 significantly different miRNAs and 125 significantly different lncRNAs were screened. According to the significantly different mRNA, miRNA and lncRNA screened from NB and 306 and NB and 306 after virus inoculation, the mRNA-miRNA-lncRNA regulatory network was constructed before and after virus inoculation, and the BmBCAT-Bomo_chr7_8305-MSTRG.3236.2 regulatory axis was screened from them, and it was found that BmBCAT was not Bomo_chr7_8305 regulated in the genetic background, after viral infection, MSTRG.3236.2 competes for binding Bomo_chr7_8305 regulates BmBCAT. The whole transcriptome sequencing results were verified by qPCR and the time-series expression analysis was performed to prove the reliability of the regulatory network. The BmBCAT-Bomo_chr7_8305-MSTRG.3236.2 regulatory axis may play a potential role in the interaction between silkworms and BmNPV. These results provide new insights into the interaction mechanism between silkworms and BmNPV.

The Bombyx mori G protein β subunit 1 (BmGNβ1) gene inhibits BmNPV infection

G protein β subunit 1 (GNβ1) has several functions, including cell growth regulation, the control of second messenger levels, and ion channel switching. Previous transcriptome analyses in our laboratory have shown that BmGNβ1 transcription is reduced following infection with Bombyx mori nucleopolyhedrovirus (BmNPV), but it is unknown what role this gene may have in the host response to BmNPV infection. In this study, the BmGNβ1 gene was cloned using the RACE method. After BmNPV infection, BmGNβ1 was downregulated in Baiyu strains in tissues such as the hemolymph and midgut. Indirect immunofluorescence showed that BmGNβ1 was localized to the cytoplasm. We further constructed a BmGNβ1-pIZ/V5-His-mCherry overexpression plasmid and designed siRNA to evaluate the role of BmGNβ1 in host response to infection. The results showed that BmGNβ1 overexpression inhibited BmNPV proliferation, while knockdown of BmGNβ1 was correlated with increased BmNPV proliferation. The siRNA-mediated reduction of BmGNβ1 was correlated with an increase in BmNPV infection of BmN cells, increased BmNPV vp39 transcription, and reduced survival time of BmNPV-infected B. mori. Overexpression of BmGNβ1 in BmN cells was also correlated with apoptosis and a modification in transcript levels of genes involved in host response to BmNPV infection (PI3K, AKT, Bmp53, BmFOXO, Caspase-1, Bmp21, BmPKN and BmCREB), suggesting that BmGNβ1 may influence the apoptotic host response of infected B. mori through the PI3K-AKT pathway. This study provides potential targets and theoretical support for breeding BmNPV-resistant silkworm varieties.

Transcriptome analysis indicates the mechanisms of BmNPV resistance in Bombyx mori midgut

Bombyx mori nucleopolyhedrovirus (BmNPV) caused serious economic losses in sericulture. Analyzing the molecular mechanism of silkworms (B. mori) resistance to BmNPV is of great significance for the prevention and control of silkworm virus diseases and the biological control of agricultural lepidopteran pests. In order to clarify the defense mechanisms of silkworms against BmNPV, we constructed a near isogenic line BC8 with high resistance to BmNPV through the highly BmNPV-resistant strain NB and the highly BmNPV-susceptible strain 306. In this study, RNA-Seq technique was used to analyze the transcriptome level differences in the midgut of BC8 and 306 following BmNPV infection. A total of 1350 DEGs were identified. Clustering analysis showed that these genes could be divided into 8 clusters with different expression patterns. Functional annotations based on GO and KEGG analysis indicated that they were involved in various metabolism pathways. Finally, 32 BmNPV defense responsive genes were screened. They were involved in metabolism, reactive oxygen species (ROS), signal transduction and immune response, and insect hormones. The further verification shows that HSP70 should participate in resistance responses of anti-BmNPV. These findings have paved the way in further functional characterization of candidate genes and subsequently can be used in breeding of BmNPV resistance dominant silkworms.

Bmo-miR-3351 modulates glutathione content and inhibits BmNPV proliferation by targeting BmGSTe6 in Bombyx mori

MicroRNAs (miRNAs) are small non-coding RNAs that play pivotal roles in the host response to invading pathogens. Among these pathogens, Bombyx mori nucleopolyhedrovirus (BmNPV) is one of the main causes of substantial economic losses in sericulture, and there are relatively few studies on the specific functions of miRNAs in the B. mori-BmNPV interaction. Therefore, we conducted transcriptome sequencing to identify differentially expressed (DE) messenger RNAs (mRNAs) and miRNAs in the midgut of 2 B. mori strains (BmNPV-susceptible strain P50 and BmNPV-resistant strain A35) after BmNPV infection. Through correlation analysis of the miRNA and mRNA data, we identified a comprehensive set of 21 miRNAs and 37 predicted target mRNAs. Notably, miR-3351, which has high expression in A35, exhibited remarkable efficacy in suppressing BmNPV proliferation. Additionally, we confirmed that miR-3351 binds to the 3' untranslated region (3' UTR) of B. mori glutathione S-transferase epsilon 6 (BmGSTe6), resulting in its downregulation. Conversely, BmGSTe6 displayed an opposite expression pattern to miR-3351, effectively promoting BmNPV proliferation. Notably, BmGSTe6 levels were positively correlated with glutathione S-transferase activity, consequently influencing intracellular glutathione content in the infected samples. Furthermore, our investigation revealed the protective role of glutathione against BmNPV infection in BmN cells. In summary, miR-3351 modulates glutathione content by downregulating BmGSTe6 to inhibit BmNPV proliferation in B. mori. Our findings enriched the research on the role of B. mori miRNAs in the defense against BmNPV infection, and suggests that the antiviral molecule, glutathione, offers a novel perspective on preventing viral infection in sericulture.

Single-nucleus sequencing of silkworm larval midgut reveals the immune escape strategy of BmNPV in the midgut during the late stage of infection

The midgut is an important barrier against microorganism invasion and proliferation, yet is the first tissue encountered when a baculovirus naturally invades the host. However, only limited knowledge is available how different midgut cell types contribute to the immune response and the clearance or promotion of viral infection. Here, single-nucleus RNA sequencing (snRNA seq) was employed to analyze the responses of various cell subpopulations in the silkworm larval midgut to B. mori nucleopolyhedrovirus (BmNPV) infection. We identified 22 distinct clusters representing enteroendocrine cells (EEs), enterocytes (ECs), intestinal stem cells (ISCs), Goblet cell-like and muscle cell types in the BmNPV-infected and uninfected silkworm larvae midgut at 72 h post infection. Further, our results revealed that the strategies for immune escape of BmNPV in the midgut at the late stage of infection include (1) inhibiting the response of antiviral pathways; (2) inhibiting the expression of antiviral host factors; (3) stimulating expression levels of genes promoting BmNPV replication. These findings suggest that the midgut, as the first line of defense against the invasion of the baculovirus, has dual characteristics of "resistance" and "tolerance". Our single-cell dataset reveals the diversity of silkworm larval midgut cells, and the transcriptome analysis provides insights into the interaction between host and virus infection at the single-cell level.