Comparative Transcriptome Analysis of Bombyx mori (Lepidoptera) Larval Midgut Response to BmNPV in Susceptible and Near-Isogenic Resistant Strains

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.

Proteomics Analysis to Explore the Resistance Genes of Silkworm to Bombyx mori Nuclear Polyhedrosis Virus

The resistance of silkworms to Bombyx mori nuclear polyhedrosis virus (BmNPV) is controlled by a major dominant gene and multiple modifying genes. Given the presence of modified genes, it is difficult to determine the main gene by positional cloning. In this study, the main anti-BmNPV gene of BmNPV-resistant silkworm variety N was introduced into the susceptible variety Su to breed the near-isogenic line SuN with BmNPV resistance. The infection process of BmNPV in the hemolymph of Su and SuN was analyzed using the cell analysis system TissueFAXS PLUS. According to the law of infection and proliferation, hemolymph was extracted every 6 h for two-dimensional electrophoresis (2-DE) analysis and quantitative real-time polymerase chain reaction (qRT-PCR). Seven DEPs were found in comparisons between Su and SuN by 2-DE analysis. Among them, acid phosphatase, storage protein, and phenoloxidase can prevent pathogen invasion, which may play a role against BmNPV. Polyamine oxidase plays an important role in energy metabolism, which may be indirectly involved in the process of resisting BmNPV. Most of the transcriptional expression profiles of the seven DEPs were consistent with the 2-DE results. This study can provide a reference for the identification of anti-BmNPV genes and the breeding of BmNPV-resistant silkworm varieties.

Transcriptomics and interactomics during the primary infection of an SfNPV baculovirus on Spodoptera frugiperda larvae

The fall armyworm (FAW), Spodoptera frugiperda, has been the most devastating pest of corn as well as of other crops in America, and more recently in Africa and Asia. The development of resistance to chemical insecticides led the search for environmentally friendly biological alternatives such as baculoviruses. This study focuses on the primary infection of the baculovirus SfNPV-Ar in the FAW's midgut epithelium, by analyzing the differential expression of transcripts in excised midguts at 6, 12, and 24 h post-infection (hpi), and predicted their interactions. Interaction of viral factors with the infected midgut tissue could alters various cellular processes, such as the apoptotic system due to the up-regulation observed of FABP at 6 hpi and of HSP90 at 24 hpi, along with the down-regulated PRX at 6 hpi and FABP transcripts between 12 and 24 hpi. Changes in transcript regulation could affect the cellular architecture of infected cells due to up-regulation of ARP 2/3 at 6 and 12 hpi, followed by down-regulation at 24 hpi. In relation to protein folding proteins, HSP90 was up-regulated at 24 hpi and PDI was down-regulated between 6 and 12 hpi. With respect to metabolism and cellular transport, AcilBP and ATPS0 were up regulated at 6 hpi and 12 hpi, respectively. In reference to transcription and translation up-regulation of RPL11 at 6 hpi and of FPN32 and RPL19 at 24 hpi was detected, as well as the down-regulation of RPL19 at 6 hpi, of PDI and RPL7 at 12 hpi, and of FABP at 24 hpi. In conclusion, gene regulation induced by viral infection could be related to the cytoskeleton and cellular metabolism as well as to oxidative stress, apoptosis, protein folding, translation, and ribosomal structure. The results presented in this work are an approach to understanding how the virus takes control of the general metabolism of the insect host during the primary infection period.

Bombyx mori Ecdysone Receptor B1 May Inhibit BmNPV Infection by Triggering Apoptosis

Bombyx mori nucleopolyhedrovirus (BmNPV) is a serious threat to sericulture. Nevertheless, no effective control strategy is currently available. The innate immunity of silkworm is critical in the antiviral process. Exploring its molecular mechanism provides theoretical support for the prevention and treatment of BmNPV. Insect hormone receptors play an essential role in regulating host immunity. We found a correlation between Bombyx mori ecdysone receptor B1 (BmEcR-B1) and BmNPV infection, whereas the underlying mechanism remains unclear. In this study, the expression patterns and sequence characteristics of BmEcR-B1 and its isoform, BmEcR-A, were initially analyzed. BmEcR-B1 was found to be more critical than BmEcR-A in silkworm development and responses to BmNPV. Moreover, RNAi and an overexpression in BmN cells showed BmEcR-B1 had antiviral effects in the presence of 20-hydroxyecdysone (20E); Otherwise, it had no antiviral activity. Furthermore, BmEcR-B1 was required for 20E-induced apoptosis, which significantly suppressed virus infection. Finally, feeding 20E had no significant negative impacts on larval growth and the cocoon shell, suggesting the regulation of this pathway has practical value in controlling BmNPV in sericulture. The findings of this study provide important theoretical support for understanding the mechanism of the silkworm innate immune system in response to BmNPV infection.

Identification of Host Molecules Involved in the Proliferation of Nucleopolyhedrovirus in Bombyx mori

The Bombyx mori nucleopolyhedrovirus (BmNPV), a foodborne infectious virus, is the pathogen causing nuclear polyhedrosis and high lethality in the silkworm. In this study, we characterized the molecules involved in BmNPV-silkworm interaction by RNA sequencing of the fat body isolated from the virus-susceptible strain P50. Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation showed that the upregulated differentially expressed genes (DEGs) were mainly involved in translation, signal transduction, folding, sorting, and degradation, as well as transport and catabolism, while the downregulated DEGs were predominantly enriched in the metabolism of carbohydrates, amino acids, and lipids at 72 h post BmNPV infection. Knockout of the upregulated somatomedin-B and thrombospondin type-1 domain-containing protein, probable allantoicase, trifunctional purine biosynthetic protein adenosine-3, and Psl and pyoverdine operon regulator inhibited the proliferation of BmNPV, while knockout of the downregulated clip domain serine protease 3 and carboxylesterase clade H, member 1 promoted it. The molecules herein identified provide a foundation for developing strategies and designing drugs against BmNPV.

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

Bombyx mori nucleopolyhedrovirus (BmNPV) is one of the most serious pathogens restricting the sustainable development of the sericulture industry. Currently, there is no effective treatment for BmNPV infection in silkworms, and the mechanism underlying BmNPV resistance in silkworms is also not clear. In this study, comparative transcriptome analyses were carried out in midguts of two silkworm varieties, namely BaiyuN, which is a resistance variety, and Baiyu, which is a susceptible variety, at five different time points (i.e., 0, 1, 3, 6, and 9 h) post-BmNPV infection to detect the early-stage transcriptional changes in these silkworms. In total, 1911 and 1577 differentially expressed genes (DEGs) were identified in the Baiyu and BaiyuN varieties, respectively, involving a total of 48 metabolic pathways. Of these pathways, eight were shared by the Baiyu and BaiyuN varieties in response to BmNPV infection. Notably, four genes (i.e., BGIBMGA08815, BGIBMGA003935, BGIBMGA003571, BGIBMGA010059) were upregulated in the Baiyu variety while downregulated in the BaiyuN variety. The inhibited expression of these four genes in the resistant variety highlighted their potential roles in the resistance of early-stage viral replication. Thus, our study provided a new avenue for the further study of the mechanism underlying BmNPV infection in silkworms and the potential treatment of BmNPV infection.

Bombyx mori C-Type Lectin (BmIML-2) Inhibits the Proliferation of B. mori Nucleopolyhedrovirus (BmNPV) through Involvement in Apoptosis

C-type lectins (CTLs) are widely distributed in mammals, insects, and plants, which act as pattern recognition receptors (PRRs) to recognize pathogens and initiate immune responses. In this study, we identified a C-type lectin gene called BmIML-2 from the silkworm Bombyx mori. Its open reading frame (ORF) encodes 314 amino acids, which contain dual tandem C-type lectin-like domain (CTLD). BmIML-2 is highly expressed in the fat body and is significantly induced at 24 h after BmNPV infection. Moreover, overexpression of BmIML-2 dramatically inhibited the proliferation of BmNPV, and knockdown assay via siRNA further validated the inhibition of BmIML-2 on viral proliferation. In addition, transcript level detection of apoptosis-related genes and observation of apoptosis bodies implied that overexpression of BmIML-2 promoted BmNPV-induced apoptosis. Immunofluorescence analysis indicated that BmIML-2 distributed throughout the cytoplasm and was slightly concentrated in the cell membrane. Taken together, our results suggest that BmIML-2 could inhibit in the proliferation of BmNPV by facilitating cell apoptosis.

CRISPR/Cas9-Mediated Disruption of the lef8 and lef9 to Inhibit Nucleopolyhedrovirus Replication in Silkworms

Bombyx mori nucleopolyhedrovirus (BmNPV) is a pathogen that causes severe disease in silkworms. In a previous study, we demonstrated that by using the CRISPR/Cas9 system to disrupt the BmNPV ie-1 and me53 genes, transgenic silkworms showed resistance to BmNPV infection. Here, we used the same strategy to simultaneously target lef8 and lef9, which are essential for BmNPV replication. A PCR assay confirmed that double-stranded breaks were induced in viral DNA at targeted sequences in BmNPV-infected transgenic silkworms that expressed small guide RNAs (sgRNAs) and Cas9. Bioassays and qPCR showed that replication of BmNPV and mortality were significantly reduced in the transgenic silkworms in comparison with the control groups. Microscopy showed degradation of midgut cells in the BmNPV-infected wild type silkworms, but not in the transgenic silkworms. These results demonstrated that transgenic silkworms using the CRISPR/Cas9 system to disrupt BmNPV lef8 and lef9 genes could successfully prevent BmNPV infection. Our research not only provides more alternative targets for the CRISPR antiviral system, but also aims to provide new ideas for the application of virus infection research and the control of insect pests.

Unraveling the innate immune responses of Bombyx mori hemolymph, fat body, and midgut to Bombyx mori nucleopolyhedrovirus oral infection by metabolomic analysis

Bombyx mori nucleopolyhedrovirus (BmNPV) infection causes a series of physiological and pathological changes in Bombyx mori (B. mori). Here, a metabolomic study of the innate immunity organs including hemolymph, fat body, and midgut of the silkworm strain Dazao following BmNPV challenge was conducted to reveal the metabolic variations in B. mori. Compared to the control, 4964 and 4942 features with 4077 and 4327 high-quality features were generated under positive and negative modes, respectively, from BmNPV-infected larvae. The principal component analysis and supervised learning method using partial least squares discrimination analysis demonstrated good analytical stability and experimental reproducibility of the metabolic profiles. Based on database annotations, a total of 296, 108, and 215 differential expressed metabolites (DEMs) were identified from BmNPV-infected group of hemolymph, fat body, and midgut, respectively, which were all mainly grouped into carboxylic acids and derivatives, fatty acyls, and glycerophospholipids. Kyoto Encyclopedia of Genes and Genomes Database enrichment analysis of the DEMs showed that amino acid metabolism was increased at 24 h after BmNPV infection. BmNPV induction was adopted to significantly alter a series of immune-related pathways including phospholipase D signaling pathway, FoxO signaling pathway, metabolism of xenobiotics by cytochrome P450, melanogenesis, membrane transport, carbohydrate metabolism, and lipid metabolism. The different levels of expression of several DEMs including l-glutamate, naphthalene, 3-succinoylpyridine 1-acyl-sn-glycerol 3-phosphate, and l-tyrosine which were involved in those pathways exhibited the immune responses of B. mori to BmNPV infection. Our findings are valuable for a better understanding of the antiviral mechanism of B. mori underlying the interaction between the silkworm and BmNPV.

Bombyx mori β-1,3-Glucan Recognition Protein 4 (BmβGRP4) Could Inhibit the Proliferation of B. mori Nucleopolyhedrovirus through Promoting Apoptosis

β-1,3-glucan recognition proteins (βGRPs) as pattern recognition receptors (PRRs) play an important role in recognizing various pathogens and trigger complicated signaling pathways in insects. In this study, we identified a Bombyx mori β-1,3-glucan recognition protein gene named BmβGRP4, which showed differential expression, from a previous transcriptome database. The full-length cDNA sequence was 1244 bp, containing an open reading frame (ORF) of 1128 bp encoding 375 amino acids. BmβGRP4 was strongly expressed in the larval stages and highly expressed in the midgut of B. mori larvae in particular. After BmNPV infection, the expression of BmβGRP4 was reduced significantly in the midgut. Furthermore, a significant increase in the copy number of BmNPV was observed after the knockdown of BmβGRP4 in 5th instar larvae, while the overexpression of BmβGRP4 suppressed the proliferation of BmNPV in BmN cells. Subsequently, the expression analysis of several apoptosis-related genes and observation of the apoptosis morphology demonstrated that overexpression of BmβGRP4 facilitated apoptosis induced by BmNPV in BmN cells. Moreover, BmβGRP4 positively regulated the phosphatase and tensin homolog gene (BmPTEN), while expression of the inhibitor of apoptosis gene (BmIAP) was negatively regulated by BmβGRP4. Hence, we hypothesize that BmNPV infection might suppress BmPTEN and facilitate BmIAP to inhibit cell apoptosis by downregulating the expression of BmβGRP4 to escape host antiviral defense. Taken together, these results show that BmβGRP4 may play a role in B. mori response to BmNPV infection and lay a foundation for studying its functions.

The Antiviral Molecule 5-Pyridoxolactone Identified Post BmNPV Infection of the Silkworm, Bombyx mori

Bombyx mori nucleopolyhedrovirus (BmNPV) is a pathogen that causes great economic losses in sericulture. Many genes play a role in viral infection of silkworms, but silkworm metabolism in response to BmNPV infection is unknown. We studied BmE cells infected with BmNPV. We performed liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS)-based non-targeted metabolomics analysis of the cytosolic extract and identified 36, 76, 138, 101, 189, and 166 different molecules at 3, 6, 12, 24, 48, and 72 h post BmNPV infection (hpi) compared with 0 hpi. Compounds representing different areas of metabolism were increased in cells post BmNPV infection. These areas included purine metabolism, aminoacyl−tRNA biosynthesis, and ABC transporters. Glycerophosphocholine (GPC), 2-hydroxyadenine (2-OH-Ade), gamma-glutamylcysteine (γ-Glu-Cys), hydroxytolbutamide, and 5-pyridoxolactone glycerophosphocholine were continuously upregulated in BmE cells post BmNPV infection by heat map analysis. Only 5-pyridoxolactone was found to strongly inhibit the proliferation of BmNPV when it was used to treat BmE cells. Fewer infected cells were detected and the level of BmNPV DNA decreased with increasing 5-pyridoxolactone in a dose-dependent manner. The expression of BmNPV genes ie1, helicase, GP64, and VP39 in BmE cells treated with 5-pyridoxolactone were strongly inhibited in the BmNPV infection stage. This suggested that 5-pyridoxolactone may suppress the entry of BmNPV. The data in this study characterize the metabolism changes in BmNPV-infected cells. Further analysis of 5-pyridoxolactone, which is a robust antiviral molecule, may increase our understanding of antiviral immunity.

Bmcas-1 plays an important role in response against BmNPV infection in vitro

Apoptosis, as one kind of innate immune system, is involved in host response against pathogens innovation. Caspases play a vital role in the execution stage of host cell apoptosis. It has been reported that Bmcaspase-1 (Bmcas-1) has a close relationship with Bombyx mori nucleopolyhedrovirus (BmNPV) infection for its differentially expressed patterns after viral infection. However, its underlying response mechanism is still unclear. The significant differential expression of Bmcas-1 in different tissues of differentially resistant strains revealed its vital role in BmNPV infection. To further validate its role in BmNPV infection, budded virus (BV)-eGFP was analyzed after knockdown and overexpression of Bmcas-1 by small interfering RNA and the pIZT-mCherry vector, respectively. The reproduction of BV-eGFP obviously increased at 72 h after knockdown of Bmcas-1, and decreased after overexpression in BmN cells. Moreover, the conserved functional domain of Cas-1 among different species and the closed evolutionary relationship of Cas-1 in Lepidoptera hinted that Bmcas-1 might be associated with apoptosis, and this was also validated by the apoptosis inducer, Silvestrol, and the inhibitor, Z-DEVD-FMK. Therefore, Bmcas-1 plays an essential antiviral role by activating apoptosis, and this result lays a fundament for clarifying the molecular mechanism of silkworm in response against BmNPV infection and breeding of resistant strains.

The validation of the role of several genes related to Bombyx mori nucleopolyhedrovirus infection in vivo

Bombyx mori nucleopolyhedrovirus (BmNPV) is one of primary silkworm pathogens and causes a serious damage of cocoon losses every year. Recent years, many works have been done to clarify the silkworm anti-BmNPV mechanism, and a significant progress has been made in screening and studying of genes and proteins related to BmNPV infection, but several of them lacked the proofs in vivo. In this study, to further validate the function of seven newly reported genes in vivo, including BmAtlatin-n, Bmferritin-heavy chain (BmFerHCH), Bmthymosin (BmTHY), Bmseroin1, Bmseroin2, Bmnuclear hormone receptors 96 (BmNHR96), and BmE3 ubiquitin-protein ligase SINA-like 10 (BmSINAL10), the response of them in the midgut, fat body, and hemolymph of differentially resistant strains (resistant strain YeA and susceptible strain YeB) at 48 h following BmNPV infection were analyzed. The results showed that the relative stable or upregulated expression level of BmAtlatin-n, BmTHY, Bmseroin1, and Bmseroin2 in YeA resistant strain following BmNPV infection further indicated their antiviral role in vivo, compared with susceptible YeB strain. Moreover, the significant downregulation of BmFerHCH, BmNHR96, and BmSINAL10 in both strains following BmNPV infection revealed their role in benefiting virus infection, as well as the upregulation of BmFerHCH in YeB midgut and BmSINAL10 in YeB hemolymph. These data could be used to complementary the proofs of the function of these genes in response to BmNPV infection.

Bmapaf-1 is Involved in the Response against BmNPV Infection by the Mitochondrial Apoptosis Pathway

Simple Summary

Apaf-1 is involved in the apoptosis pathway and Bmapaf-1 showed a significant response to BmNPV infection in our previous transcriptome data. In this study, the underlying mechanism of Bmapaf-1 in response to BmNPV infection was studied. To preliminarily determine the relationship of Bmapaf-1 with BmNPV, the expression pattern of Bmapaf-1 was analyzed in different tissues of differentially resistant silkworm strains following virus infection. To further define the role of Bmapaf-1 in BmNPV infection, the alteration of BmNPV infection in BmN cells and the expression patterns of Bmcas-Nc and Bmcas-1 were analyzed following knockdown and overexpression of Bmapaf-1 using siRNA and the pIZT/V5-His-mCherry insect vector, respectively. Furthermore, to analyze whether Bmapaf-1 is involved in BmNPV infection by apoptosis, the inducer NSC348884 and inhibitor Z-DEVD-FMK were used.

Abstract

Discovery of the anti-BmNPV (Bombyx mori nuclearpolyhedrovirus) silkworm strain suggests that some kind of antiviral molecular mechanism does exist but is still unclear. Apoptosis, as an innate part of the immune system, plays an important role in the response against pathogen infections and may be involved in the anti-BmNPV infection. Several candidate genes involved in the mitochondrial apoptosis pathway were identified from our previous study. Bombyx mori apoptosis protease-activating factor-1 (Bmapaf-1) was one of them, but the antiviral mechanism is still unclear. In this study, sequences of BmApaf-1 were characterized. It was found to contain a unique transposase_1 functional domain and share high CARD and NB-ARC domains with other species. Relatively high expression levels of Bmapaf-1 were found at key moments of embryonic development, metamorphosis, and reproductive development. Further, the significant difference in expression of Bmapaf-1 in different tissues following virus infection indicated its close relationship with BmNPV, which was further validated by RNAi and overexpression in BmN cells. Briefly, infection of budded virus with enhanced green fluorescent protein (BV-EGFP) was significantly inhibited at 72 h after overexpression of Bmapaf-1, which was confirmed after knockdown of Bmapaf-1 with siRNA. Moreover, the downstream genes of Bmapaf-1, including Bmnedd2-like caspase (BmNc) and Bmcaspase-1 (Bmcas-1), were upregulated after overexpression of Bmapaf-1 in BmN cells, which was consistent with the RNAi results. Furthermore, the phenomenon of Bmapaf-1 in response to BmNPV infection was determined to be related to apoptosis using the apoptosis inducer NSC348884 and inhibitor Z-DEVD-FMK. Therefore, Bmapaf-1 is involved in the response against BmNPV infection by the mitochondrial apoptosis pathway. This result provides valuable data for clarifying the anti-BmNPV mechanism of silkworms and breeding of resistant silkworm strains.

Identification and Functional Analysis of BmNPV-Interacting Proteins From Bombyx mori (Lepidoptera) Larval Midgut Based on Subcellular Protein Levels

Bombyx mori nucleopolyhedrovirus (BmNPV) is a major pathogen causing severe economic loss. However, the molecular mechanism of silkworm resistance to BmNPV and the interactions of this virus with the host during infection remain largely unclear. To explore the virus-binding proteins of silkworms, the midgut subcellular component proteins that may interact with BmNPV were analyzed in vitro based on one- and two-dimensional electrophoresis and far-western blotting combined with mass spectrometry (MS). A total of 24 proteins were determined to be specifically bound to budded viruses (BVs) in two subcellular fractions (mitochondria and microsomes). These proteins were involved in viral transportation, energy metabolism, apoptosis and viral propagation, and they responded to BmNPV infection with different expression profiles in different resistant strains. In particular, almost all the identified proteins were downregulated in the A35 strain following BmNPV infection. Interestingly, there were no virus-binding proteins identified in the cytosolic fraction of the silkworm midgut. Two candidate proteins, RACK1 and VDAC2, interacted with BVs, as determined with far-western blotting and reverse far-western blotting. We speculated that the proteins interacting with the virus could either enhance or inhibit the infection of the virus. The data provide comprehensive useful information for further research on the interaction of the host with BmNPV.

Identification of FerLCH, isolation of ferritin and functional analysis related to interaction with pathogens in Eri-silkworm, Samia cynthia ricini

Ferritin is a ubiquitous and conserved iron storage protein that plays a significant role in host detoxification, iron storage, and immune response. Although ferritin has been studied in many species, little is known about its role in the Eri-silkworm (Samia cynthia ricini). In this study, the ferritin light-chain subunit gene, named ScFerLCH, was identified from S. c. ricini. The full-length gene, ScFerLCH, was 1,155 bp and encoded a protein consisting of 231 amino acids with a deduced molecular weight of 26.38 kDa. Higher ScFerLCH expression levels were found in the midgut, silk gland, and fat body by quantitative reverse-transcription polymerase chain reaction and western blot analysis. Injection of Staphylococcus aureus and Pseudomonas aeruginosa could induce upregulation of ScFerLCH in the hemolymph, fat body, and midgut, indicating that ScFerLCH may contribute to the host defense against invading pathogens. In addition, the native ferritin protein was isolated from S. c. ricini by native polyacrylamide gel electrophoresis and its two subunits, ferritin heavy-chain subunit (ScFerHCH) and ferritin light-chain subunit (ScFerLCH), were identified by mass spectrometry. Specifically, we found that recombinant ferritin subunits could self-assemble into a protein complex in vitro; moreover, both recombinant subunits and the protein complex were found to bind different bacteria, including Escherichia coli, P. aeruginosa, S. aureus, and Bacillus subtilis. However, bactericidal tests showed that the protein complex could not inhibit the growth of bacteria directly. Taken together, our results suggest that ScFerritin might play an important role in mediating molecular interaction with pathogens.

Transcriptome Analyses of Diaphorina citri Midgut Responses to Candidatus Liberibacter Asiaticus Infection

The Asian citrus psyllid (ACP), Diaphorina citri Kuwayama (Hemiptera: Liviidae), is an important transmission vector of the citrus greening disease Candidatus Liberibacter asiaticus (CLas). The D. citri midgut exhibits an important tissue barrier against CLas infection. However, the molecular mechanism of the midgut response to CLas infection has not been comprehensively elucidated. In this study, we identified 778 differentially expressed genes (DEGs) in the midgut upon CLas infection, by comparative transcriptome analyses, including 499 upregulated DEGs and 279 downregulated DEGs. Functional annotation analysis showed that these DEGs were associated with ubiquitination, the immune response, the ribosome, endocytosis, the cytoskeleton and insecticide resistance. KEGG enrichment analysis revealed that most of the DEGs were primarily involved in endocytosis and the ribosome. A total of fourteen DEG functions were further validated by reverse transcription quantitative PCR (RT-qPCR). This study will contribute to our understanding of the molecular interaction between CLas and D. citri.

A Novel Digestive Proteinase Lipase Member H-A in Bombyx mori Contributes to Digestive Juice Antiviral Activity Against B. mori Nucleopolyhedrovirus

Previous studies have revealed that some proteins in Bombyx mori larvae digestive juice show antiviral activity. Here, based on the label-free proteomics data, BmLipase member H-A (BmLHA) was identified as being involved in the response to BmNPV infection in B. mori larvae digestive juice. In the present study, a gene encoding the BmLHA protein in B. mori was characterized. The protein has an open reading fragment of 999 bp, encoding a predicted 332 amino acid residue-protein with a molecular weight of approximately 35.9 kDa. The phylogenetic analysis revealed that BmLHA shares a close genetic distance with Papilio xuthus Lipase member H-A. BmLHA was highly expressed in the middle part of the B. mori gut, and the expression level increased with instar rising in larvae. There was higher expression of BmLHA in A35 than in P50 strains, and it was upregulated in both A35 and P50 strains, following BmNPV infection. The expression level of VP39 decreased significantly in appropriate recombinant-BmLHA-treated groups compared with the PBS-treated group in B. mori larvae and BmN cells. Meanwhile, overexpression of BmLHA significantly reduced the infectivity of BmNPV in BmN cells. These results indicated that BmLHA did not have digestive function but had anti-BmNPV activity. Taken together, our work provides valuable data for the clarification of the molecular characterization BmLHA and supplements research on proteins of anti-BmNPV activity in B. mori.

Bmserpin2 Is Involved in BmNPV Infection by Suppressing Melanization in Bombyx mori

Melanization, an important defense response, plays a vital role in arthropod immunity. It is mediated by serine proteases (SPs) that convert the inactive prophenoloxidase (PPO) to active phenoloxidase (PO) and is tightly regulated by serine protease inhibitors (serpins) which belong to a well distributed superfamily in invertebrates, participating in immune mechanisms and other important physiological processes. Here, we investigated the Bmserpin2 gene which was identified from a transcriptome database in response to Bombyx mori nucleopolyhedrovirus (BmNPV) infection. Quantitative real-time polymerase chain reaction (qRT-PCR) results showed that Bmserpin2 was expressed in all tissues, with maximum expression in fat body. Upon BmNPV infection, the expression of Bmserpin2 was up-regulated in P50 (susceptible strain) and BC9 (resistant strain) in haemocytes, fat body and the midgut. However, up-regulation was delayed in BC9 (48 or 72 h), in contrast to P50 (24 h), after BmNPV infection. Meanwhile, Bmserpin2 could delay or inhibit melanization in silkworm haemolymph. Significant increased PO activity can be observed in Bmserpin2-depleted haemolymph under NPV infection. Furthermore, the viral genomic DNA copy number was decreased in Bmserpin2-depleted haemolymph. We conclude that Bmserpin2 is an inducible gene which might be involved in the regulation of PPO activation and suppressed melanization, and have a potential role in the innate immune system of B. mori.

Quantitative label-free proteomic analysis reveals differentially expressed proteins in the digestive juice of resistant versus susceptible silkworm strains and their predicted impacts on BmNPV infection

Bombyx mori nucleopolyhedrovirus (BmNPV) is a major pathogen causing severe economic loss. Previous studies have revealed that some proteins in silkworm digestive juice show antiviral activity. In this study, antiviral activity examination of different resistant strains showed that the digestive juice of the resistant strain (A35) had higher inhibition to virus than the susceptible strain (P50). Subsequently, the label-free quantitative proteomics was used to study the midgut digestive juice response to BmNPV infection in P50 and A35 strains. A total of 98 proteins were identified, of which 80 were differentially expressed proteins (DEPs) with 54 enzymes and 26 nonenzymatic proteins by comparing the proteomes of infected and non-infected P50 and A35 silkworms. These DEPs are mainly involved in metabolism, proteolysis, neuroactive ligand receptor interaction, starch and sucrose metabolism and glutathione metabolism. After removing the genetic background and individual immune stress response proteins, 9 DEPs were identified potentially involved in resistance to BmNPV. Further studies showed that a serine protease, an alkaline phosphatase and serine protease inhibitor 2 isoform X1 were differentially expressed in A35 compared to P50 or post BmNPV infection. Taken together, these results provide insights into the potential mechanisms for silkworm digestive juice to provide resistance to BmNPV infection. Signifcance: Bombyx mori nucleopolyhedrovirus (BmNPV) is highly pathogenic, which has a great impact on the sericulture. BmNPV entered the midgut lumen and exposed to digestive juices after oral infection. Previous studies have revealed that some proteins in silkworm digestive juice show antiviral activity, however, current information on the digestive juice proteome of high resistant silkworm strain after BmNPV challenge compared to susceptible strain is incomprehensive. Here, we combined label-free quantification method, bioinformatics, RT-qPCR and western blot analysis and found that BmNPV infection causes some protein changes in the silkworm midgut digestive juice. The DEPs were identified in the digestive juices of different resistant strains following BmNPV infection, and screened out some proteins potentially related to resistance to BmNPV. Three important differentially expression proteins were validated by independent approaches. These findings uncover the potential role of silkworm digestive juice in providing resistance to BmNPV and supplemented the profile of the proteome of the digestive juices in B. mori.

Study on the Role of Cytc in Response to BmNPV Infection in Silkworm, Bombyx mori (Lepidoptera)

Bombyx mori nucleopolyhedrovirus (BmNPV) is one of the primary pathogens of the silkworm. Cytochrome c (cytc) showed a significant response to BmNPV infection in our previous transcriptome study. However, little is known about the role of Bombyx mori cytc (Bmcytc) in resistance to BmNPV infection. In this study, the expression levels analysis of Bmcytc showed stable expression levels in selected tissues of the resistant strain AN following BmNPV infection, while there was downregulation in the susceptible strain p50, except in the malpighian tubule. To further study the role of Bmcytc in viral infection, Bmcytc was knocked down with siRNA in vitro, resulting in significant downregulation of selected downstream genes of the mitochondrial pathway, including Bmapaf, Bmcaspase-Nc, and Bmcaspase-1; this was also confirmed by overexpression of Bmcytc using the pIZT/V5-His-mCherry insect vector, except Bmcaspase-1. Moreover, knockdown of Bmcytc significantly promoted the infection process of BmNPV in vitro, while the infection was inhibited by overexpression of Bmcytc at the early stage and subsequently increased rapidly. Based on these results, we concluded that Bmcytc plays a vital role in BmNPV infection by regulating the mitochondrial apoptosis pathway. Our work provides valuable data for the clarification of the mechanism of silkworm resistance to BmNPV infection.

Specific genes related to nucleopolyhedrovirus in Bombyx mori susceptible and near-isogenic resistant strains through comparative transcriptome analysis

Bombyx mori nucleopolyhedrovirus (BmNPV) is one of the primary pathogens that causes severe economic losses to sericulture. Comparative transcriptomics analysis has been widely applied to explore the antiviral mechanism in resistant strains. Here, to identify genes involved in BmNPV infection, we identified differentially expressed genes (DEGs) and performed weighted gene co-expression network analysis (WGCNA) between two Bombyx mori strains: strain 871 (susceptible to BmNPV infection) and the near-isogenic strain 871C (resistant to BmNPV). Our results showed that 400 genes were associated with resistance in strain 871C, and 76 genes were related to susceptibility in strain 871. In addition, the correlation analysis of DEGs and WGCNA showed that 40 genes related to resistance were highly expressed in the resistant strain. Among them, gene BGIBMGA004291 was the most noticeable. We further identified the effect of gene BGIBMGA004291, which encoded a multiprotein bridge factor 2 (MBF2) family member (MBF2-10), on viral infection in cells. Our data suggested that MBF2-10 inhibited viral infection. Taken together, this study showed specific module trait correlations related to viral infection in strains 871 and 871C, and we identified a resistance-related gene. These findings suggested promising candidate genes with antiviral activity, aiding in the analysis of the antiviral molecular mechanisms in resistant strains.

The regulation of crecropin-A and gloverin 2 by the silkworm Toll-like gene 18 wheeler in immune response

The innate immune system is conserved among different insect species in its response to microorganism infection. The transmembrane receptors of the Toll superfamily play an important role in activating immune response, however, the function of silkworm Toll family member 18 Wheeler (18 W) remained unclear. Here, the 18w gene in silkworm was characterized. A relatively high transcription level of Bm18w mRNA was found in Malpighian tubules, and in eggs, larvae pre-molt to fourth instar, pupae and adults. When silkworm larvae were infected with E. coli or S. aureus, Bm18w showed a significant response, especially to E. coli, but did not have antibacterial activity. To further identify the downstream antimicrobial peptide genes of Bm18w, expression of Bm18w was knocked down with siRNA in vitro, resulting in significant decreases of cecropin-A, gloverin 2, and moricin B3. The overexpression of Bm18w was carried out using pIZT/V5-His-mCherry insect vector in BmN cells and significant upregulation of cecropin-A and gloverin 2 was detected, as well as upregulation of attacin and defensin. Based on the results, we concluded that Bm18w is involved in response to bacterial infection by selectively inducing the expression of antimicrobial peptide genes, especially cecropin-A and gloverin 2. This study provides valuable data to supplement understanding of the immune pathway of the silkworm.

Knockdown of BmTCP-1β Delays BmNPV Infection in vitro

The molecular mechanism of silkworm resistance to Bombyx mori nucleopolyhedrovirus (BmNPV) infection remains unclear. The chaperonin containing t-complex polypeptide 1 (TCP-1) is essential for the folding of tubulin and actin to produce stable and functional competent protein conformation. However, little is known about this protein in silkworm. In the present study, a gene encoding the TCP-1β protein in silkworm was characterized, which has an open reading fragment of 1,611 bp encoding a predicted 536 amino acid residue-protein with a molecular weight of approximately 57.6 kDa containing a Cpn60_TCP1 functional domain. The sequence conservation is 81.52%. The highest level of BmTCP-1β mRNA expression was found in the midgut, while the lowest was in the hemolymph. To further study the function of BmTCP-1β, expression was knocked down with siRNA in vitro, resulting in significant downregulation of the selected cytoskeletal-related genes, actin and tubulin, which was also confirmed by overexpression of BmTCP-1β in BmN cells using the pIZT/V5-His-mCherry insect vector. Moreover, knockdown of BmTCP-1β significantly prolonged the infection process of BmNPV in BmN cells, which was also verified by overexpression of BmTCP-1β in BmN cells. Based on the results of the present study, we concluded that BmTCP-1β plays a vital role in BmNPV infection by regulating the expression of tubulin and actin. Taken together, our work provides valuable data for the clarification of the molecular mechanism of silkworm resistance to BmNPV infection.

Molecular Characterization of Two Mitogen-Activated Protein Kinases: p38 MAP Kinase and Ribosomal S6 Kinase From Bombyx mori (Lepidoptera: Bombycidae), and Insight Into Their Roles in Response to BmNPV Infection

Proteins p38 map kinase and ribosomal S6 kinase (S6K) as members of mitogen-activated protein kinases (MAPKs) play important roles against pathogens. In this study, Bmp38 and BmS6K were identified as differentially expressed proteins from iTRAQ database. Bmp38 and BmS6K were expressed, and recombinant proteins were purified. The bioinformatics analysis showed that both proteins have serine/threonine-protein kinases, catalytic domain (S_TKc) with 360 and 753 amino acids, respectively. The real-time quantitative polymerase chain reaction (RT-qPCR) results suggest that Bmp38 and BmS6K had high expression in the midgut and hemolymph. The comparative expression level of Bmp38 and BmS6K in BC9 was upregulated than in P50 in the midgut after Bombyx mori nucleopolyhedrovirus (BmNPV) infection. Western bolt results showed a positive correlation between RT-qPCR and iTRAQ data for Bmp38, but BmS6K data showed partial correlation with iTRAQ. Injection of anti-Bmp38 and anti-BmS6K serum suggested that Bmp38 may be involved against BmNPV infection, whereas BmS6K may require phosphorylation modification to inhibit BmNPV infection. Taken together, our results suggest that Bmp38 and BmS6k might play an important role in innate immunity of silkworm against BmNPV.

C-lysozyme contributes to antiviral immunity in Bombyx mori against nucleopolyhedrovirus infection

Lysozymes is a ubiquitous immune effector that is widely distributed in both vertebrates and invertebrates. Previous reports have shown that lysozymes significantly inhibit viral infections in vertebrates. However, the antiviral effects of lysozymes in invertebrates remain unclear. Here, we investigated the role of lysozymes in Bombyx mori (B. mori) response to viral infection by overexpressing B. mori C-lysozyme (BmC-LZM) in larvae and cells. We found that BmC-LZM was up-regulated in cells in response to viral infection. Indeed, the overexpressing of BmC-LZM significantly inhibited viral replication in cells during late-stage infection. However, this effect was reversed by BmC-LZM mRNA. BmC-LZM was successfully overexpressed in B. mori strain 871 using Baculovirus Expression Vector System (BEVS). This overexpression markedly reduced viral proliferation and increased larval survival percentage. Thus, BmC-LZM inhibited viral replication both in vivo and in vitro, indicating that BmC-LZM is involved in the insect immune response to viral infection. Our results provide a basis for further applications of lysozymes.

Comparative transcriptome analysis of silkworm, Bombyx mori colleterial gland suggests their functional role in mucous secretion

Colleterial glands (CG) present in the body of adult female of Bombyx mori, which can help adhere eggs on the surface of the host plants. Although this organ has been known for centuries, only morphology and its secretions have been studied. Their gene expression profiles and physiological roles remain largely unknown. Aided by high-throughput next generation sequencing (NGS), we reported the comparative transcriptome analysis of CG isolated from the H9 and the P50 strains of Bombyx mori. A total of 19,896,957 and 20,446,366 clean reads were obtained from CG of H9 and the P50 strains, respectively; then differential expression analysis was performed, and 1,509 differentially expressed genes (DEGs) were identified. Among them, 1,001 genes are up-regulated and 508 genes are down-regulated in P50 individuals compared with H9 individuals. The enrichment of GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) of DEGs confirmed that many DEGs were associated with "Amino acid transport and metabolism", "Nucleotide transport and metabolism", and "Inorganic ion transport and metabolism", 25 of the DEGs related to the "ECM-receptor interaction passway", "sphingolipid metabolism passway", and "amino sugar and nucleotide sugar metabolism passway" were potentially involved in the process of CG development and mucus secretion. According to these data, we hypothesized that CG play an important role in providing favorable physiological environment for the glue secretion formation. In addition, GO enrichment and differential expression analysis of the DEGs in the CG indicate that this gland may be involved in the transporting of small solutes such as sugars, ions, amino acids and nucleotide sugar to the CG. Our findings lay the foundation for further research on CG function.

Response analysis of host Spodoptera exigua larvae to infection by Heliothis virescens ascovirus 3h (HvAV-3h) via transcriptome

Heliothis virescens ascovirus 3 h (HvAV-3h), a dsDNA insect virus, belonging to the family Ascoviridae, can infect caterpillars of several Noctuidae species by ovipositing parasitoid wasps. In order to provide a comprehensive overview of the interactive responses of host larvae after infection by the ascovirus, a transcriptome analysis of Spodoptera exigua to HvAV-3h was conducted from 6 to 168 hours post infection (hpi). Approximately 101.64 Gb of RNA sequencing (RNA-seq) data obtained from infected and uninfected S. exigua larvae were used to perform a de novo transcriptome assembly, which generated approximately 62,258 S. exigua unigenes. Using differential gene expression analysis, it was determined that the majority of host transcripts were down-regulated beginning at 6 hpi and continuing throughout the infection period, although there was an increase in up-regulated unigene number during the 12 to 72 hpi stage. It is noteworthy that the most abundantly enriched pathways in KEGG annotation were Metabolism terms, indicating that the host larval metabolic mechanisms were highly influenced post HvAV-3h infection. In addition, the host cuticle protein encoding unigenes were highly down-regulated in most of the situations, suggesting that the host larval cuticle synthesis were inhibited by the viral infection.

Comparative transcriptome analysis reveals significant metabolic alterations in eri-silkworm (Samia cynthia ricini) haemolymph in response to 1-deoxynojirimycin

Samia cynthia ricini (Lepidoptera: Saturniidae) is an important commercial silk-producing insect; however, in contrast to the silkworm, mulberry leaves are toxic to this insect because the leaves contain the component 1-deoxynojirimycin (DNJ). A transcriptomic analysis of eri-silkworm haemolymph was conducted to examine the genes related to different metabolic pathways and to elucidate the molecular mechanism underlying eri-silkworm haemolymph responses to DNJ. Eight hundred sixty-five differentially expressed genes (DEGs) were identified, among which 577 DEGs were up-regulated and 288 DEGs were down-regulated in the 2% DNJ group compared to control (ddH₂O) after 12h. Based on the results of the functional analysis, these DEGs were associated with ribosomes, glycolysis, N-glycan biosynthesis, and oxidative phosphorylation. In particular, according to the KEGG analysis, 138 DEGs were involved in energy metabolism, glycometabolism and lipid metabolism, and the changes in the expression of nine DEGs were confirmed by reverse transcription quantitative PCR (RT-qPCR). Thus, DNJ induced significant metabolic alterations in eri-silkworm haemolymph. These results will lay the foundation for research into the toxic effects of DNJ on eri-silkworm as a model and provide a reference for the exploitation of new drugs in humans.

SWATH-based quantitative proteomics reveals the mechanism of enhanced Bombyx mori nucleopolyhedrovirus-resistance in silkworm reared on UV-B treated mulberry leaves

Bombyx mori nucleopolyhedrovirus (BmNPV) is one of the most acute infectious diseases in silkworm, which has led to great economic loss in sericulture. Previous study showed that the content of secondary metabolites in mulberry leaves, particularly for moracin N, was increased after UV-B irradiation. In this study, the BmNPV resistance of silkworms reared on UV-B treated and moracin N spread mulberry leaves was improved. To uncover the mechanism of enhanced BmNPV resistance, silkworm midguts from UV-B treated mulberry leaves (BUM) and moracin N (BNM) groups were analyzed by SWATH-based proteomic technique. Of note, the abundance of ribosomal proteins in BUM and BNM groups was significantly changed to maintain the synthesis of total protein levels and cell survival. While, cytochrome c oxidase subunit II, calcium ATPase and programmed cell death 4 involved in apoptotic process were up-regulated in BNM group. Expressions of lipase-1, serine protease precursor, Rab1 protein, and histone genes were increased significantly in BNM group. These results suggest that moracin N might be the main active component in UV-B treated mulberry leaves which could improve the BmNPV-resistance of silkworm through promoting apoptotic cell death, enhancing the organism immunity, and regulating the intercellular environment of cells in silkworm. It also presents an innovative process to reduce the mortality rate of silkworms infected with BmNPV.

iTRAQ-based quantitative proteomics analysis of molecular mechanisms associated with Bombyx mori (Lepidoptera) larval midgut response to BmNPV in susceptible and near-isogenic strains

Bombyx mori nucleopolyhedrovirus (BmNPV) has been identified as a major pathogen responsible for severe economic loss. Most silkworm strains are susceptible to BmNPV, with only a few highly resistant strains thus far identified. Here we investigated the molecular basis of silkworm resistance to BmNPV using susceptible (the recurrent parent P50) and resistant (near-isogenic line BC9) strains and a combination of iTRAQ-based quantitative proteomics, reverse-transcription quantitative PCR and Western blotting. By comparing the proteomes of infected and non-infected P50 and BC9 silkworms, we identified 793 differentially expressed proteins (DEPs). By gene ontology and KEGG enrichment analyses, we found that these DEPs are preferentially involved in metabolism, catalytic activity, amino sugar and nucleotide sugar metabolism and carbon metabolism. 114 (14.38%) DEPs were associated with the cytoskeleton, immune response, apoptosis, ubiquitination, translation, ion transport, endocytosis and endopeptidase activity. After removing the genetic background and individual immune stress response proteins, we identified 84 DEPs were found that are potentially involved in resistance to BmNPV. Further studies showed that a serine protease was down-regulated in P50 and up-regulated in BC9 after BmNPV infection. Taken together, these results provide insights into the molecular mechanism of silkworm response to BmNPV.

BIOLOGICAL SIGNIFICANCE

Bombyx mori nucleopolyhedrovirus (BmNPV) is highly pathogenic, causing serious losses in sericulture every year. However, the molecular mechanisms of BmNPV infection and host defence remain unclear. Here we combined quantitative proteomic, bioinformatics, RT-qPCR and Western blotting analyses and found that BmNPV invasion causes complex protein alterations in the larval midgut, and that these changes are related to cytoskeleton, immune response, apoptosis, ubiquitination, translation, ion transport, endocytosis and endopeptidase activity. Five important differentially expression proteins were validation by independent approaches. These finding will help address the molecular mechanisms of silkworm resistance to BmNPV and provide a molecular target for resisting BmNPV.