Comparative analysis of Bombyx mori nucleopolyhedrovirus responsive genes in fat body and haemocyte of B. mori resistant and susceptible strains

Comprehensive transcriptome sequencing of silkworm Midguts: Uncovering extensive isoform diversity and alternative splicing in BmNPV-Sensitive and BmNPV-resistant strains

The silkworm, Bombyx mori, stands out as one of the few economically valuable insects within the realm of model organisms. However, Bombyx mori nucleopolyhedrovirus (BmNPV) poses a significant threat, decreasing the quality and quantity of silkworm cocoons. Over the past few decades, a multitude of researchers has delved into the mechanisms that underlie silkworm resistance to BmNPV, employing diverse methodologies and approaching the problem from various angles. Despite this extensive research, the role of alternative splicing (AS) in the silkworm's response to BmNPV infection has been largely unexplored. This study leveraged both third-generation (Oxford Nanopore Technologies) and second-generation (Illumina) high-throughput sequencing technologies to meticulously identify and analyze AS patterns in the context of BmNPV response, utilizing two distinct silkworm strains-the susceptible strain 306 and the resistant strain NB. Consequently, we identified five crucial genes (Dsclp, LOC692903, LOC101743583, LOC101742498, LOC101743809) that are linked to the response to BmNPV infection through AS and differential expression. Additionally, a thorough comparative analysis was conducted on their diverse transcriptomic expression profiles, including alternative polyadenylation, simple sequence repeats, and transcription factors.

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.

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.

A Protein Asteroid with PIN Domain in Silkworm Bombyx mori Is Involved in Anti-BmNPV Infection

Nuclease is a type of protein that degrades nucleic acids, which plays an important role in biological processes, including RNA interference efficiency and antiviral immunity. However, no evidence of a link between nuclease and Bombyx mori nucleopolyhedrovirus (BmNPV) infection in silkworm B. mori has been found. In this study, a protein asteroid (BmAst) containing the PIN domain and XPG domain was identified in silkworm B. mori. BmAst gene was highest expressed in hemocytes and fat body of the 5th instar larvae, and high expression in the pupa stage. The transcriptional levels of the BmAst gene in 5th instar larvae were significantly induced by BmNPV or dsRNA. After knocking down BmAst gene expression by specific dsRNA, the proliferation of BmNPV in B. mori was increased significantly, whereas the survival rate of larvae was significantly lower when compared with the control. Our findings indicate that BmAst is involved in silkworm resistance to BmNPV infection.

The Mechanisms of Silkworm Resistance to the Baculovirus and Antiviral Breeding

Silkworm (Bombyx mori) is not only an economic insect but also a model organism for life science research. Bombyx mori nucleopolyhedrovirus (BmNPV) disease is a major infectious disease in the world's sericulture industry. The cocoon loss caused by this disease accounts for more than 60% of the total loss caused by all silkworm diseases. To date, there has been no effective solution for preventing and treating this disease. The most effective measure is to breed disease-resistant varieties. The quickest way to breed disease-resistant varieties is to apply genetic modification. However, this requires that we obtain disease resistance genes and know the mechanism of disease resistance. Since the discovery of disease-resistant resources in 1989, scholars in the sericulture industry around the world have been inspired to search for resistance genes. In the past two decades, with the help of multi-omics technologies, screening of resistance genes, gene localization, protein modification, virus-host interactions, etc., researchers have found some candidate genes that have been proposed to function at the cellular or individual level. Several disease-resistant varieties have been obtained and used in production through hybrid breeding, RNA interference, and genetic modification. This article summarizes and reviews the discovery of and research advances related to silkworm resistance to BmNPV. It is anticipated that the review will inspire scientific researchers to continue searching for disease resistance genes, clarify the molecular mechanism of silkworm disease resistance, and promote disease-resistant silkworm breeding.

The antiviral role of NF-κB-mediated immune responses and their antagonism by viruses in insects

The antiviral role of innate immune responses mediated by the NF-κB family of transcription factors is well established in vertebrates but was for a long time less clear in insects. Insects encode two canonical NF-κB pathways, the Toll and Imd ('immunodeficiency') pathways, which are best characterised for their role in antibacterial and antifungal defence. An increasing body of evidence has also implicated NF-κB-mediated innate immunity in antiviral responses against some, but not all, viruses. Specific pattern recognition receptors (PRRs) and molecular events leading to NF-κB activation by viral pathogen-associated molecular patterns (PAMPs) have been elucidated for a number of viruses and insect species. Particularly interesting are recent findings indicating that the cyclic GMP-AMP synthase (cGAS)-stimulator of interferon genes (STING) pathway detects viral RNA to activate NF-κB-regulated gene expression. We summarise the literature on virus-NF-κB pathway interactions across the class Insecta, with a focus on the dipterans Drosophila melanogaster and Aedes aegypti. We discuss potential reasons for differences observed between different virus-host combinations, and highlight similarities and differences between cGAS-STING signalling in insects versus vertebrates. Finally, we summarise the increasing number of known molecular mechanisms by which viruses antagonise NF-κB responses, which suggest that NF-κB-mediated immunity exerts strong evolutionary pressures on viruses. These developments in our understanding of insect antiviral immunity have relevance to the large number of insect species that impact on humans through their transmission of human, livestock and plant diseases, exploitation as biotechnology platforms, and role as parasites, pollinators, livestock and 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.

Proteomic analysis of Laodelphax striatellus in response to Rice stripe virus infection reveal a potential role of ZFP36L1 in restriction of viral proliferation

Persistent plant viruses multiply and circulate inside insect vectors following the route of midgut-hemolymph-salivary gland. Currently, how viruses interact with insect vectors after they are released into hemolymph is not entirely clear. In this study, we found that the hemolymph and fat body (HF) contained the highest Rice stripe virus (RSV) levels. Proteomic analysis on RSV-free and RSV-infected HF identified 156 differentially expressed proteins (DEPs), with the majority of them participating in metabolism, transportation, and detoxification. The RNA binding protein esf2 was the most downregulated protein. Knocking down the expression of esf2 did not influence the RSV burden, but caused the lethal effect to L. striatellus. In contrast, the mRNA decay protein ZFP36L1 was 69% more abundant upon RSV infection, and suppression of ZFP36L1 significantly increased the RSV burden. Our results reveal the potential role of ZFP36L1 in restricting the viral proliferation, and provide valuable clues for unravelling the interaction between RSV and L. striatellus in HF. SIGNIFICANCE: More than 76% of plant viruses are transmitted by insect vectors. For persistent propagative transmission, plant viruses multiply and circulate inside insects following the route of midgut-hemolymph-salivary gland. However, how viruses interact with vector insects after they are released into hemolymph is not entirely clear. Our study investigated the influence of rice stripe virus (RSV) on insect hemolymph and fat body by iTRAQ labeling method. Among the 156 differentially expressed proteins (DEPs) identified, two proteins associated with mRNA metabolism were selected for function analysis. We found that the mRNA decay activator protein ZFP36L1 influenced the RSV proliferation, and RNA binding protein esf2 caused the lethal effect to L. striatellus. Our results provide valuable clues for unveiling the interaction between RSV and L. striatellus, and might be useful in pest management.

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.

Integrative analysis of circRNA/miRNA/mRNA regulatory network reveals the potential immune function of circRNAs in the Bombyx mori fat body

Bombyx mori nucleopolyhedrosis virus (BmNPV) is one of the greatest threats to sustainable development of the sericulture industry. Circular RNA (circRNA), a type of non-coding RNA, has been shown to play important roles in gene expression regulation, immune response, and diseases. The fat body is a tissue with both metabolic and immune functions. To explore the potential immune function of circRNAs, we analyzed differentially expressed (DE)circRNAs, microRNAs(miRNAs), and mRNAs in the B. mori fat body in response to BmNPV infection using high-throughput RNA sequencing. A total of 77 DEcircRNAs, 32 DEmiRNAs, and 730 DEmRNAs that are associated with BmNPV infection were identified. We constructed a DEcircRNA/DEmiRNA/DEmRNA and DEcircRNA/DEmiRNA/BmNPV gene regulatory network and validated the differential expression of circ_0001432 and its corresponding miRNA (miR-2774c and miR-3406-5p) and mRNA (778467 and 101745232) in the network. Tissue-specific expression of circ_0001432 and its expression at different time points were also examined. KEGG pathway analysis of DEmRNAs, target genes of DEmiRNAs, and host genes of DEcircRNAs in the network showed that these genes were enriched in several metabolic pathways and signaling pathways, which could play important roles in insect immune responses. Our results suggest that circRNA could be involved in immune responses of the B. mori fat body and help in understanding the molecular mechanisms underlying silkworm-pathogen interactions.

Antimicrobial Peptides as Potential Antiviral Factors in Insect Antiviral Immune Response

Antimicrobial peptides (AMPs) with antiviral activity (antiviral peptides: AVPs) have become a research hotspot and already show immense potential to become pharmaceutically available antiviral drugs. AVPs have exhibited huge potential in inhibiting viruses by targeting various stages of their life cycle. Insects are the most speciose group of animals that inhabit almost all ecosystems and habitats on the land and are a rich source of natural AMPs. However, insect AVP mining, functional research, and drug development are still in their infancy. This review aims to summarize the currently validated insect AVPs, explore potential new insect AVPs and to discuss their possible mechanism of synthesis and action, with a view to providing clues to unravel the mechanisms of insect antiviral immunity and to develop insect AVP-derived antiviral drugs.

Functional analysis of Dicer-2 gene in Bombyx mori resistance to BmNPV virus

Bombyx mori nucleopolyhedrovirus (BmNPV) is a primary pathogen in silkworm, and the molecular mechanism of B. mori defense to BmNPV infection is still unclear. RNA interference (RNAi) is well-known as an intracellular conserved mechanism that is critical in gene regulation and cell defense. The antiviral RNAi pathway processes viral double-stranded RNA (dsRNA) into viral small interfering RNAs that guide the recognition and cleavage of complementary viral target RNAs. In this study, a Dicer-2 (Dcr2) gene was identified in B. mori and its antiviral function was explored. Dcr2 messenger RNA (mRNA) expression was the highest in hemocytes and expressed in all stages of silkworm growth. After infection with BmNPV, the expression of Dcr2 mRNA was significantly increased after infection in midgut and hemocytes. The expression of Dcr2 was significantly upregulated by injecting dsRNA (dsBmSPH-1) into silkworm after 48 hr. Knocking down the expression level of Dcr2 using specific dsRNA in silkworm, which modestly enhanced the production of viral genomic DNA. Our results suggested that the Dcr2 gene in B. mori plays an important role in against BmNPV invasion.

Increased expression of Suppressor of cytokine signaling 2 (BmSOCS2) is correlated with suppression of Bombyx mori nucleopolyhedrovirus replication in silkworm larval tissues and cells

The resistance of silkworm to infection by Bombyx mori nuclear polyhedrosis virus (BmNPV) is a main focus of sericultural research. Previously, a BmNPV-resistant strain, NB, was identified among a collection of Chinese silkworm strains in our lab. To better understand the molecular mechanism of NB strain resistance, the patterns of host immune response gene transcription in resistant (NB) and susceptible (306) strains were examined. Quantative real-time PCR (qRT-PCR) revealed that multiple insect innate immune signaling pathways (Toll, Imd and JAK/STAT) were strongly activated upon infection with BmNPV. Notably, Suppressor of cytokine signaling 2 (BmSOCS2) mRNA expression was significantly up-regulated in midgut tissues of the resistant NB strain, suggesting that the BmSOCS2 gene product may be involved in host immune defense against BmNPV infection. A significant inhibition of BmNPV replication was also observed in BmN cells transfected with a vector encoding BmSOCS2. The results suggest that BmSOCS2 is a key gene involved in the resistance of the NB silkworm strain to BmNPV infection.

Adenosine Receptor Modulates Permissiveness of Baculovirus (Budded Virus) Infection via Regulation of Energy Metabolism in Bombyx mori

Although the modulation of host physiology has been interpreted as an essential process supporting baculovirus propagation, the requirement of energy supply for host antivirus reactions could not be ruled out. Our present study showed that metabolic induction upon AcMNPV (budded virus) infection of Bombyx mori stimulated virus clearance and production of the antivirus protein, gloverin. In addition, we demonstrated that adenosine receptor signaling (AdoR) played an important role in regulating such metabolic reprogramming upon baculovirus infection. By using a second lepidopteran model, Spodoptera frugiperda Sf-21 cells, we demonstrated that the glycolytic induction regulated by adenosine signaling was a conservative mechanism modulating the permissiveness of baculovirus infection. Another interesting finding in our present study is that both BmNPV and AcMNPV infection cause metabolic activation, but it appears that BmNPV infection moderates the level of ATP production, which is in contrast to a dramatic increase upon AcMNPV infection. We identified potential AdoR miRNAs induced by BmNPV infection and concluded that BmNPV may attempt to minimize metabolic activation by suppressing adenosine signaling and further decreasing the host's anti-baculovirus response. Our present study shows that activation of energy synthesis by adenosine signaling upon baculovirus infection is a host physiological response that is essential for supporting the innate immune response against infection.

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.

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.

Resistant silkworm strain block viral infection independent of melanization

Melanization mediated by the prophenoloxidase-activating system (proPO) is an important immune response in invertebrates. However, the role of melanization on viral infection has not been wildly revealed in Bombyx mori (B. mori), silkworm. Here, we investigated the extent of melanization of susceptible (871) and resistant (near-isogenic line 871C) B. mori strains. The result showed that the extent of melanization was significantly higher in the susceptible strain than in the resistant strain. A majority of Serpins were up-regulated in the resistant strain through iTRAQ-based quantitative proteomics, comparing with susceptible strain. Our data further identified that Serpin-5, Serpin-9 and Serpin-19 reduced PO activity, indicating that the menlanization pathway was inhibited in the resistant strain. Moreover, our results indicated that the hemolymph of 871 reduced viral infection in the presence of PTU, indicating that melanization of 871 strain hemolymph blocked viral infection. However, viral infection was significantly suppressed in the hemolymph of 871C strain regardless of the presence of PTU or not, which implied that the resistant strain inhibited viral infection independent of the melanization pathway. Taken together, our findings indicated that the melanization pathway was inhibited in resistant strain. These results expend the analysis of melanization pathway in insects and provide insights into understanding the antiviral mechanism.

Circular RNA alterations in the Bombyx mori midgut following B. mori nucleopolyhedrovirus infection

Thus far, no systematic studies have examined circRNA expression profiles in the silkworm following B.mori nucleopolyhedrovirus (BmNPV) infection. To explore the expression patterns of circRNAs in the silkworm midgut following BmNPV infection, circRNAs in normal midguts and BmNPV-infected midguts were analyzed by high-throughput sequencing. A total of 353 circRNAs were significantly differentially expressed, of which 241 were upregulated and 112 were downregulated following infection. GO annotation and KEGG pathways analyses of these circRNAs showed that many key immunity pathways and metabolism pathways were enriched in the BmNPV-infected midguts. The potential roles of the predicted targets of the miRNAs that interacted with the circRNAs showed that ubiquitin, apoptosis, and endocytosis signaling pathways were enriched significantly by BmNPV infection.

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.

De novo characterization of Phenacoccus solenopsis transcriptome and analysis of gene expression profiling during development and hormone biosynthesis

The cotton mealybug Phenacoccus solenopsis is a devastating pest of cotton causing tremendous loss in the yield of crops each year. Widespread physiological and biological studies on P. solenopsis have been carried out, but the lack of genetic information has constrained our understanding of the molecular mechanisms behind its growth and development. To understand and characterize the different developmental stages, RNA-Seq platform was used to execute de-novo transcriptome assembly and differential gene expression profiling for the eggs, first, second, third instar and adult female stages. About 182.67 million reads were assembled into 93,781 unigenes with an average length of 871.4 bp and an N50 length of 1899 bp. These unigenes sequences were annotated and classified by performing NCBI non-redundant (Nr) database, Kyoto Encyclopedia of Genes and Genomes (KEGG) and Clusters of Orthologous Groups (COG), Gene ontology (GO), the Swiss-Prot protein database (Swiss-Prot), and nearest related organism Acyrthosiphon pisum (pea aphid) database. To get more information regarding the process of metamorphosis, we performed a pairwise comparison of four developmental stages and obtained 29,415 differentially expressed genes. Some of the differentially expressed genes were associated with functional protein synthesis, anti-microbial protection, development and hormone biosynthesis. Functional pathway enrichment analysis of differentially expressed genes showed the positive correlation with specific physiological activities of each stage, and these results were confirmed by qRT-PCR experiments. This study gives a valuable genomics resource of P. solenopsis covering all its developmental stages and will promote future studies on biological processes at the molecular level.

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.

Discovery of anti-viral molecules and their vital functions in Bombyx mori

The silkworm Bombyx mori (B. mori), a lepidopteran model organism, has become an important model for molecular biology researches with its genome completely sequenced. Silkworms confront different types of virus diseases, mainly including those caused by Bombyx mori nucleopolyhedrovirus (BmNPV), Bombyx mori densovirus type 1 (BmDNV-1), Bombyx mori bidesovirus (BmBDV) which was termed as Bombyx mori densovirus type 2 (BmDNV-2) or Bombyx mori parvo-like virus (BmPLV) before in sericulture. B. mori offers excellent models to study the molecular mechanisms of insect innate immune responses to viruses. A variety of molecules and pathways have been identified to be involved in the immune responses in the silkworm to viruses, such as the antimicrobial peptides, prophenoloxidase-activating system, apoptosis, ROS, small RNA and related molecules. Here in this review, we summarize the current research advances in molecules involved in silkworm anti-virus pathways. Moreover, taking BmNPV as an example, we proposed a schematic model of molecules and pathways involved in silkworm immune responses against virus infection. We hope this review can facilitate further study of antiviral mechanisms in silkworm, and provide a reference for virus diseases in other organisms.

Describing the Diapause-Preparatory Proteome of the Beetle Colaphellus bowringi and Identifying Candidates Affecting Lipid Accumulation Using Isobaric Tags for Mass Spectrometry-Based Proteome Quantification (iTRAQ)

Prior to entering diapause, insects must prepare themselves physiologically to withstand the stresses of arresting their development for a lengthy period. While studies describing the biochemical and cellular milieu of the maintenance phase of diapause are accumulating, few studies have taken an “omics” approach to describing molecular events during the diapause preparatory phase. We used isobaric tags and mass spectrometry (iTRAQ) to quantitatively compare the expression profiles of proteins identified during the onset of diapause preparation phase in the heads of adult female cabbage beetles, Colaphellus bowringi. A total of 3,175 proteins were identified, 297 of which were differentially expressed between diapause-destined and non-diapause-destined female adults and could therefore be involved in diapause preparation in this species. Comparison of identified proteins with protein function databases shows that many of these differentially expressed proteins enhanced in diapause destined beetles are involved in energy production and conversion, carbohydrate metabolism and transport, and lipid metabolism. Further hand annotation of differentially abundant peptides nominates several associated with stress hardiness, including HSPs and antioxidants, as well as neural development. In contrast, non-diapause destined beetles show substantial increases in cuticle proteins, suggesting additional post-emergence growth. Using RNA interference to silence a fatty acid-binding protein (FABP) that was highly abundant in the head of diapause-destined females prevented the accumulation of lipids in the fat body, a common product of diapause preparation in this species and others. Surprisingly, RNAi against the FABP also affected the transcript abundance of several heat shock proteins. These results suggest that the identified differentially expressed proteins that play vital roles in lipid metabolism may also contribute somehow to enhanced hardiness to environmental stress that is characteristic of diapause.

In vitro production of baculoviruses: identifying host and virus genes associated with high productivity

Baculoviruses are recognized as viral workhorses of biotechnology, being used for production of vaccines, complex recombinant proteins, gene delivery vectors' and safe biological pesticides. Improving production yields and understanding the interactions of the virus and its host cell are important aspects of ensuring baculovirus-based processes are commercially competitive. This study aims at potential optimization of host cells used in in vitro virus production by systemically investigating changes in host gene expression in response to virus replication and transcription inside host cells. The study focuses on in vitro interactions of the Helicoverpa armigera virus with Helicoverpa zea insect cells. We used 22 genome-wide microarrays to simultaneously measure both virus and host genes in infected cells in multiple batch suspension cultures, representing high- and low-producing infection conditions. Among 661 differentially expressed genes, we identified a core set of 59 host genes consistently overexpressed post infection, with strong overrepresentation of genes involved in retrotransposition, protein processing in the endoplasmic reticulum, and ubiquitin-mediated proteolysis. Applying a whole genome correlation network analysis to link gene expression to productivity, we revealed 18 key genes significantly associated to virus yield. In addition, this study is among the first to perform a genome-wide expression study for a major baculovirus group II strain, the H. armigera virus, extending current understanding of baculovirus-insect interactions, which mainly focuses on group I viruses.

Digital gene expression analysis of Helicoverpa armigera in the early stage of infection with Helicoverpa armigera nucleopolyhedrovirus

Helicoverpa armigera single nucleocapsid nucleopolyhedrovirus (HearNPV) is an obligatory and lethal parasite of the cotton bollworm and has been extensively used in China for the control of this notorious pest. Digital gene expression (DGE) analysis was adopted for an overall comparison of transcriptome profiling between HearNPV-infected and control healthy Helicoverpa armigera larvae during an early stage post-inoculation. A total of 908 differentially expressed genes (DEGs) were identified, of which 136 were up-regulated and 597 were down-regulated. GO category and KEGG pathway analysis demonstrated that the identified DEGs involved in ribosome biogenesis, aminoacyl-tRNA biosynthesis, protein processing in endoplasmic reticulum, biosynthesis of valine, leucine, isoleucine and the spliceosome were significantly down-regulated, whereas genes involved in pancreatic secretion, protein digestion and absorption and salivary secretion showed obviously up-regulated transcription. The DEGs were verified by quantitative real-time PCR, and genes that participated in defensive response, nutritional digestion and developmental regulation exhibited specific expression patterns in a continuous time-course assessment. These results provide basic data for future research on the molecular mechanism of HearNPV infection and the interactions between lepidopteran hosts and their specific NPV parasites.

Proteomics analysis of digestive juice from silkworm during Bombyx mori nucleopolyhedrovirus infection

Previous studies have analyzed the midgut transcriptome and proteome after challenge with Bombyx mori nucleopolyhedrovirus (BmNPV), however little information is available on the digestive juice proteome after BmNPV challenge. This study investigated BmNPV infection-induced protein changes in the digestive juice of silkworms using shotgun proteomics and MS sequencing. From the digestive juice of normal third-day, fifth-instar silkworm larvae, 75 proteins were identified, 44 of which were unknown; from larvae 6 h after inoculation with BmNPV, 106 proteins were identified, of which 39 were unknown. After BmNPV challenge, more secreted proteins appeared that had antiviral and digestive features. GO annotation analysis clustered most proteins in the lumen into catalytic, binding, and metabolic processes. Numerous proteins were reported to have BmNPV interactions. Hsp70 protein cognate, lipase-1, and chlorophyllide A-binding protein precursor were upregulated significantly after BmNPV challenge. Levels of trypsin-like serine protease, beta-1,3-glucanase, catalase, and serine protease transcripts decreased or were not significantly change after BmNPV challenge. Taken together, these findings provided insights into the interaction between host and BmNPV and revealed potential functions of digestive juice after per os BmNPV infection.

A hypothetical model of crossing Bombyx mori nucleopolyhedrovirus through its host midgut physical barrier

Bombyx mori nucleopolyhedrovirus (BmNPV) is a primary pathogen of silkworm (B. mori) that causes severe economic losses each year. However, the molecular mechanisms of silkworm-BmNPV interactions, especially the silkworm proteins that can interact with the virus, are still largely unknown. In this study, the total and membrane proteins of silkworm midguts were displayed using one- and two-dimensional electrophoresis. A virus overlay assay was used to detect B. mori proteins that specifically bind to BmNPV particles. Twelve proteins were located and identified using mass spectrometry, and the different expression of the corresponding genes in BmNPV susceptible and resistant silkworm strains also indicated their involvement in BmNPV infection. The 12 proteins are grouped based on their potential roles in viral infection, for example, endocytosis, intracellular transportation, and host responses. Based on these results, we hypothesize the following: I) vacuolar ATP synthase catalytic subunit A and subunit B may be implicated in the process of the membrane fusion of virus and the release of the nucleocapsid into cytoplasm; II) actin, enolase and phosphoglycerate kinase are cytoskeleton associated proteins and may play an important role in BmNPV intracellular transportation; III) mitochondrial prohibitin complex protein 2, ganglioside-induced differentiation-associated protein, calreticulin, regucalcin-like isoform X1 and 60 kDa heat shock protein are involved in cell apoptosis regulation during BmNPV infection in larvae midguts; IV) ribosomal P0 may be associated with BmNPV infection by regulating gene expression of BmNPV; V) arginine kinase has a role in the antiviral activities against BmNPV. Our work should prove informative by providing multiple protein targets and a novel direction to investigate the molecular mechanisms of the interactions between silkworms and BmNPV.

Genomic insights into the serine protease gene family and expression profile analysis in the planthopper, Nilaparvata lugens

Background

The brown planthopper (Nilaparvata lugens) is one of the most destructive rice plant pests in Asia. N. lugens causes extensive damage to rice by sucking rice phloem sap, which results in hopper burn (complete death of the rice plants). Despite its importance, little is known about the digestion, development and defense mechanisms of this hemimetabolous insect pest. In this study, we aim to identify the serine protease (SP) and serine protease homolog (SPH) genes, which form a large family in eukaryotes, due to the potential for multiple physiological roles. Having a fully sequenced genome for N. lugens allows us to perform in-depth analysis of the gene structures, reveal the evolutionary relationships and predict the physiological functions of SP genes.

Results

The genome- and transcriptome-wide analysis identified 90 putative SP (65) and SPH (25) genes in N. lugens. Detailed gene information regarding the exon-intron organization, size, distribution and transcription orientation in the genome revealed that many SP/SPH loci are closely situated on the same scaffold, indicating the frequent occurrence of gene duplications in this large gene family. The gene expression profiles revealed new findings with regard to how SPs/SPHs respond to bacterial infections as well as their tissue-, development- and sex-specific expressions.

Conclusions

Our findings provide comprehensive gene sequence resources and expression profiles of the N. lugens SP and SPH genes, which give insights into clarifying the potentially functional roles of these genes in the biological processes including development, digestion, reproduction and immunity.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-507) contains supplementary material, which is available to authorized users.

bmnpv-miR-3 facilitates BmNPV infection by modulating the expression of viral P6.9 and other late genes in Bombyx mori

During the last decade, microRNAs (miRNAs) have emerged as fine tuners of gene expression in various biological processes including host-pathogen interactions. Apart from the role of host encoded miRNAs in host-virus interactions, recent studies have also indicated the key role of virus-encoded miRNAs in the regulation of host defense responses. In the present study, we show that bmnpv-miR-3, a Bombyx mori nucleopolyhedrovirus (BmNPV) encoded miRNA, regulates the expression of DNA binding protein (P6.9) and other late genes, vital for the late stage of viral infection in the host, Bombyx mori. We have performed both cell culture and in vivo experiments to establish the role of bmnpv-miR-3 in the infection cycle of BmNPV. Our findings showed that bmnpv-miR-3 expresses during early stage of infection, and negatively regulates the expression of P6.9. There was an upregulation in P6.9 expression upon blocking of bmnpv-miR-3 by Locked Nucleic Acid (LNA), whereas overexpression of bmnpv-miR-3 resulted in a decreased expression of P6.9. Besides, a remarkable enhancement and reduction in the viral loads were observed upon blocking and overexpression of bmnpv-miR-3, respectively. Furthermore, we have also assessed the host immune response using one of the Lepidoptera-specific antimicrobial proteins, Gloverin-1 upon blocking and overexpression of bmnpv-miR-3, which correlated viral load with the host immune response. All these results together; clearly imply that bmnpv-miR-3-mediated controlled regulation of BmNPV late genes in the early stage of infection helps BmNPV to escape the early immune response from the host.

Expression analysis of several antiviral related genes to BmNPV in different resistant strains of silkworm, Bombyx mori

Bombyx mori L. (Lepidoptera: Bombycidae) nucleopolyhedrovirus (BmNPV) is a highly pathogenic virus in the sericultural industry, often causing severe damage leading to large economic losses. The immune mechanisms of B. mori against this virus remain obscure. Previous studies had demonstrated Bmlipase-1, BmNox and Bmserine protease-2 showing antiviral activity in vitro, but data on the transcription levels of these proteins in different resistant strains were not reported. In order to determine the resistance level of the four different strains (P50, A35, A40, A53) and gain a better understanding of the mechanism of resistance to BmNPV in B. mori, the relative expression level of the genes coding the three antiviral proteins in larval haemolymph and midgut of different B. mori strains resistant to BmNPV was determined. The results showed that these genes expressed significantly higher in the resistant strains compared to the susceptible strain, and the differential expression levels were consistent with the LC50 values in different strains. The transcription level of the target genes almost all up-regulated in the larvae midgut and down-regulated in the haemolymph. The results indicate the correlation of these genes to BmNPV resistance in B. mori.

The progress and future of enhancing antiviral capacity by transgenic technology in the silkworm Bombyx mori

Bombyx mori is a common lepidopteran model and an important economic insect for silk production. B. mori nucleopolyhedrovirus (BmNPV) is a typical pathogenic baculovirus that causes serious economic losses in sericulture. B. mori and BmNPV are a model of insect host and pathogen interaction including invasion of the host by the pathogen, host response, and enhancement of host resistance. The antiviral capacity of silkworms can be improved by transgenic technology such as overexpression of an endogenous or exogenous antiviral gene, RNA interference of the BmNPV gene, or regulation of the immune pathway to inhibit BmNPV at different stages of infection. Antiviral capacity could be further increased by combining different methods. We discuss the future of an antiviral strategy in silkworm, including possible improvement of anti-BmNPV, the feasibility of constructing transgenic silkworms with resistance to multiple viruses, and the safety of transgenic silkworms. The silkworm model could provide a reference for disease control in other organisms.

Genome scale transcriptomics of baculovirus-insect interactions

Baculovirus-insect cell technologies are applied in the production of complex proteins, veterinary and human vaccines, gene delivery vectors' and biopesticides. Better understanding of how baculoviruses and insect cells interact would facilitate baculovirus-based production. While complete genomic sequences are available for over 58 baculovirus species, little insect genomic information is known. The release of the Bombyx mori and Plutella xylostella genomes, the accumulation of EST sequences for several Lepidopteran species, and especially the availability of two genome-scale analysis tools, namely oligonucleotide microarrays and next generation sequencing (NGS), have facilitated expression studies to generate a rich picture of insect gene responses to baculovirus infections. This review presents current knowledge on the interaction dynamics of the baculovirus-insect system' which is relatively well studied in relation to nucleocapsid transportation, apoptosis, and heat shock responses, but is still poorly understood regarding responses involved in pro-survival pathways, DNA damage pathways, protein degradation, translation, signaling pathways, RNAi pathways, and importantly metabolic pathways for energy, nucleotide and amino acid production. We discuss how the two genome-scale transcriptomic tools can be applied for studying such pathways and suggest that proteomics and metabolomics can produce complementary findings to transcriptomic studies.

Gene cloning and expression patterns of two prophenoloxidases from Catantops pinguis (Orthoptera: Catantopidae)

In insect, fat body plays major roles in insect innate immunity. Phenoloxidase (PO) is an important component in insect innate immunity and is necessary for acclimatization. In our study, two prophenoloxidase (PPO) subunits were obtained from fat body of Catantops pinguis (Stål). The full-length cDNA sequence of one PPO (CpPPO1) consisted of 2347 bp with an open reading frame (ORF) of 2187 bp encoding 728 amino acids, while the other subunit (CpPPO2) had a full length of 2445 bp, encoding 691 amino acids. Both the PPO gene products are predicted to possess all the structural features of other PPO members, including two putative tyrosinase copper-binding motifs with six highly conserved histidine residues and a thiolester-like motif. Tissue distribution analysis showed that both PPO mRNAs were abundantly expressed in the fat body among 11 tissues examined, and they were transiently up-regulated after Escherichia coli infection, consistent with them being immune-responsive genes. Total levels of CpPPO1 and CpPPO2 mRNA transcripts were much higher in first instar larvae and adults. A much higher transcript level of CpPPO1 was detected in several months, while there were extremely high mRNA expression levels of CpPPO2 in January, July, October, and December. The above results suggested that PPO from fat body might also bring significant function during the processes of development and acclimatization for C. pinguis.

Genome scale analysis of differential mRNA expression of Helicoverpa zea insect cells infected with a H. armigera baculovirus

Knowledge of baculovirus-insect host interactions at a genome-scale level is important for developing a number of baculovirus-based applications, but the gathering of such knowledge is hindered by the lack of genomic sequences in most insect hosts. In this study, expression kinetics of 24,206 Helicoverpa zea insect transcripts and 134 Helicoverpa armigera nucleopolyhedrovirus (HearNPV) genes at 0, 12, 24 and 48 h post-infection (hpi) were simultaneously analyzed using microarrays, which were developed from sequences obtained by deep transcriptome sequencing. Host genes in pathways important for infection such as those for energy generation, anti-viral peptides, apoptosis, detoxification, DNA polymerase activities, RNA polymerase activities, translation initiation, protein processing and cell cycle arrest were identified. Differential expression was linked to changes in the number of intracellular and extracellular viral genomes and occlusion bodies. The first comprehensive elucidation of HearNPV-H. zea expression kinetics was obtained.

Microarray analysis of gene expression profile in resistant and susceptible Bombyx mori strains reveals resistance-related genes to nucleopolyhedrovirus

To investigate the molecular mechanism of silkworm resistance to BmNPV infection, we constructed a near-isogenic line (BC8) with BmNPV resistance using highly resistant (NB) and highly susceptible parental strains (306). We investigated variations in the gene expression in the midguts of BmNPV-infected BC8 and 306 at 12 h pi using the microarray. 92 differentially expressed genes were identified. Real-time qPCR analysis confirmed that 10 genes were significantly up-regulated or down-regulated in the midguts of BC8 and NB compared to 306. To our knowledge, we first defined the role of the amino acid transporter and 26S proteasome in insect antiviral. However, serine protease was not completely consistent with data of reported previously in insect antiviral. The role of the 5 genes (Bm123, Bm122, COP β', aquaporin, glycoside hydrolases) was also demonstrated in insect antiviral. Our results provided new insights into the molecular mechanism of the Bombyx mori immune response against BmNPV infection.

Dynamic interactions between Bombyx mori nucleopolyhedrovirus and its host cells revealed by transcriptome analysis

Although microarray and expressed sequence tag (EST)-based approaches have been used to profile gene expression during baculovirus infection, the response of host genes to baculovirus infection and the interaction between baculovirus and its host remain largely unknown. To determine the host response to Bombyx mori nucleopolyhedrovirus infection and the dynamic interaction between the virus and its host, eight digital gene expression libraries were examined in a Bm5 cell line before infection and at 1.5, 3, 6, 12, 24, 48, and 96 h postinfection. Gene set enrichment analysis of differentially expressed genes at each time point following infection showed that gene sets including cytoskeleton, transcription, translation, energy metabolism, iron ion metabolism, and the ubiquitin-proteasome pathway were altered after viral infection. In addition, a time course depicting protein-protein interaction networks between the baculovirus and the host were constructed and revealed that viral proteins interact with a multitude of cellular machineries, such as the proteasome, cytoskeleton, and spliceosome. Several viral proteins, including IE2, CG30, PE38, and PK-1/2, were predicted to play key roles in mediating virus-host interactions. Based on these results, we tested the role of the ubiquitin-proteasome pathway and iron ion metabolism in the viral infection cycle. Treatment with a proteasome inhibitor and deferoxamine mesylate in vitro and in vivo confirmed that these pathways regulate viral infection. Taken together, these findings provide new insights into the interaction between the baculovirus and its host and identify molecular mechanisms that can be used to block viral infection and improve baculovirus expression systems.

Low cost antiviral activity of Plodia interpunctella haemolymph in vivo demonstrated by dose dependent infection

Given the ubiquity of infectious disease it is important to understand the way in which hosts defend themselves and any costs that they may pay for this defence. Despite this, we know relatively little about insect immune responses to viruses when compared to their well-characterized responses to other pathogens. In particular it is unclear whether there is significant haemocoelic response to viral infection. Here we directly examine this question by examining whether there is a dose-dependency in infection risk when a DNA virus is injected directly into the haemocoel. Infection from direct injection into the haemocoel showed a clear dose dependency that is indicative of an active intrahaemocoelic immune response to DNA viruses in insects. In contrast to the natural oral infection route, we found no measurable sublethal effects in the survivors from direct injection. This suggests that the immune responses in the haemocoel are less costly than those that occur earlier.