Comparative Subcellular Proteomics Analysis of Susceptible and Near-isogenic Resistant Bombyx mori (Lepidoptera) Larval Midgut Response to BmNPV infection

Effects of baculovirus infection on intestinal microflora of BmNPV resistant and susceptible strain silkworm

Bombyx mori L. (Lepidoptera: Bombycidae) nucleopolyhedrovirus (BmNPV) is a serious pathogen causing huge economic losses to sericulture. There is growing evidence that the gut microbiota of silkworms plays a critical role in shaping host responses and interactions with viral infection. However, little is known about the differences in the composition and diversity of intestinal microflora, especially with respect to silkworm strain differences and BmNPV infection-induced changes. Here, we aim to explore the differences between BmNPV-resistant strain A35 and susceptible strain P50 silkworm and the impact of BmNPV infection on intestinal microflora in different strains. The 16S rDNA sequencing analysis revealed that the fecal microbial populations were distinct between A35 and P50 and were significantly changed post BmNPV infection in both strains. Further analysis showed that the BmNPV-resistant strain silkworm possessed higher bacterial diversity than the susceptible strain, and BmNPV infection reduced the diversity of intestinal flora assessed by feces in both silkworm strains. In response to BmNPV infection, the abundance of Muribaculaceae increased in P50 and decreased in A35, while the abundance of Enterobacteriaceae decreased in P50 and increased in A35. These results indicated that BmNPV infection had various effects on the abundance of fecal microflora in different silkworm strains. Our findings not only broadened the understanding of host-pathogen interactions but also provided theoretical help for the breeding of resistant strains and healthy rearing of silkworms based on symbiotic bacteria.

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

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

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

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

Multi-omics study and ncRNA regulation of anti-BmNPV in silkworms, Bombyx mori: an update

Bombyx mori silkworm is an important economic insect which has a significant contribution to the improvement of the economy. Bombyx mori nucleopolyhedrovirus (BmNPV) is a vitally significant purulent virus that impedes the sustainable and stable development of the silkworm industry, resulting in substantial economic losses. In recent years, with the development of biotechnology, transcriptomics, proteomics, metabolomics, and the related techniques have been used to select BmNPV-resistant genes, proteins, and metabolites. The regulatory networks between viruses and hosts have been gradually clarified with the discovery of ncRNAs, such as miRNA, lncRNA, and circRNA in cells. Thus, this paper aims to highlight the results of current multi-omics and ncRNA studies on BmNPV resistance in the silkworm, providing some references for resistant strategies in the silkworm to BmNPV.

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.

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.

Genomic Variation in the Tea Leafhopper Reveals the Basis of Adaptive Evolution

Tea green leafhopper (TGL), Empoasca onukii, is of biological and economic interest. Despite numerous studies, the mechanisms underlying its adaptation and evolution remain enigmatic. Here, we use previously untapped genome and population genetics approaches to examine how the pest adapted to different environmental variables and thus has expanded geographically. We complete a chromosome-level assembly and annotation of the E. onukii genome, showing notable expansions of gene families associated with adaptation to chemoreception and detoxification. Genomic signals indicating balancing selection highlight metabolic pathways involved in adaptation to a wide range of tea varieties grown across ecologically diverse regions. Patterns of genetic variations among 54 E. onukii samples unveil the population structure and evolutionary history across different tea-growing regions in China. Our results demonstrate that the genomic changes in key pathways, including those linked to metabolism, circadian rhythms, and immune system functions, may underlie the successful spread and adaptation of E. onukii. This work highlights the genetic and molecular bases underlying the evolutionary success of a species with broad economic impacts, and provides insights into insect adaptation to host plants, which will ultimately facilitate more sustainable pest management.

Metabolic characterization of hemolymph in Bombyx mori varieties after Bombyx mori nucleopolyhedrovirus infection by GC-MS-based metabolite profiling

The "Huakang 2" silkworm variety, bred by the Sericulture Research Institute of the Chinese Academy of Agricultural Sciences, is highly resistant to Bombyx mori nucleopolyhedrovirus (BmNPV) and effectively solves the issue of frequent Bombyx mori nuclear polyhedrosis in sericultural production. The molecular mechanism of its resistance to BmNPV, however, is still unknown. The purpose of the present study was therefore to identify these anti-BmNPV mechanisms by using metabolomics in combination with transcriptomics after subcutaneous injection of budded virus (BV) with high concentrations of BmNPV from specimens of the Baiyu N variety (which is highly resistant to BmNPV) and the Baiyu variety (which is sensitive to BmNPV). A total of 375 differential metabolites were identified, which mainly included sugars, acids, amines, alcohols, glycosides, and other small molecules. KEGG enrichment analysis and functional clustering of differential metabolites identified possible metabolic pathways, including tyrosine metabolism, oxidative phosphorylation, and alanine, aspartate, and glutamate metabolism. The differentially expressed genes (DEGs) identified by transcriptome analysis were annotated in KEGG. Association analysis showed that the metabolic pathways of different silkworm varieties are affected differently by BmNPV infection, triggering a series of complex physiological and biochemical changes in the organism. In particular, oxidative phosphorylation might be an essential pathway involved in regulation of disease resistance.

Integration of Transcriptomic and Proteomic Analyses Reveals New Insights into the Regulation of Immune Pathways in Midgut of Samia ricini upon SariNPV Infection

Simple Summary

SariNPV is one of the main pathogens of Samia ricini and its infection of Samia ricini sericulture has caused significant economic losses to society. In this study, we aim to reveal the molecular mechanism of pathogen–host interactions in SariNPV-infected S. ricini through transcriptomic and proteomic analyses. Using RNA-sequencing and iTRAQ, we mapped the differentially expressed genes (DEGs) and proteins (DEPs) that are involved in the immune responses of S. ricini upon virus invasion. Based on our analyses, we identified specific DEGs and DEPs that are involved in various essential biological signaling pathways and immune-related pathways upon SariNPV invasion. These DEGs and DEPs play an important role in triggering host immune responses to pathogens. Our study provides molecular insights into host immune responses regarding pathogen invasion, in particular, the immune response mechanism and network of S. ricini in response to SariNPV infection.

Abstract

Samia ricini nucleopolyhedrovirus (SariNPV) is one of the main pathogens of S. ricini sericulture and its infection causes severe impacts on economic sericulture development. To explore and reveal the molecular mechanisms of S. ricini in response to SariNPV infection, we employed RNA sequencing (RNA-seq), adopting isobaric tags for relative and absolute quantitation (iTRAQ), and carried out combination analysis of the obtained differentially expressed genes (DEGs) and proteins (DEPs). Through transcriptome sequencing, a total of 2535 DEGs were detected, and with iTRAQ, 434 DEPs with significant expression difference were identified. Through correlation analysis, we found that the expression trends of 116 differentially expressed proteins were the same as those of differentially expressed genes (including 106 up-regulated and 10 down-regulated). Twenty-five key differentially expressed genes (proteins) involved in several signaling and immune related pathways (mainly involving Toll, Imd, Jak-STAT and Wnt signaling pathways, as well as other immune related pathways) were screened through real-time quantitative PCR. Our results, not only provide insights into the pathogenic mechanism of SariNPV infection in ricin silkworm and the immune response mechanism within the host, but also provide a significant contribution for identifying and preventing diseases caused by SariNPV.

A Review on Transcriptional Responses of Interactions between Insect Vectors and Plant Viruses

This review provides a synopsis of transcriptional responses pertaining to interactions between plant viruses and the insect vectors that transmit them in diverse modes. In the process, it attempts to catalog differential gene expression pertinent to virus–vector interactions in vectors such as virus reception, virus cell entry, virus tissue tropism, virus multiplication, and vector immune responses. Whiteflies, leafhoppers, planthoppers, and thrips are the main insect groups reviewed, along with aphids and leaf beetles. Much of the focus on gene expression pertinent to vector–virus interactions has centered around whole-body RNA extraction, whereas data on virus-induced tissue-specific gene expression in vectors is limited. This review compares transcriptional responses in different insect groups following the acquisition of non-persistent, semi-persistent, and persistent (non-propagative and propagative) plant viruses and identifies parallels and divergences in gene expression patterns. Understanding virus-induced changes in vectors at a transcriptional level can aid in the identification of candidate genes for targeting with RNAi and/or CRISPR editing in insect vectors for management approaches.

The digestive proteinase trypsin, alkaline A contributes to anti-BmNPV activity in silkworm (Bombyx mori)

Bombyx mori nucleopolyhedrovirus (BmNPV) is a serious pathogenic microorganism that causes tremendous loss to sericulture. Previous studies have found that some proteins of serine protease family in the digestive juice of B. mori larvae have anti-BmNPV activity. In our previous publication about proteome analysis of the digestive juice of B. mori larvae, the digestive enzyme trypsin, alkaline A (BmTA) was filtered as a differentially expressed protein possibly involved in BmNPV resistance. Here, the biological characteristics and anti-BmNPV functions of BmTA were comprehensively analysed. The cDNA sequence of BmTA had an ORF of 768 nucleotides encoding 255 amino acid residues. Domain architecture analysis showed that BmTA contained a signal peptide and a typical Tryp_SPc domain. Quantitative real-time PCR analysis showed that BmTA was highly expressed in the larval stages and specifically expressed in the midgut of B. mori larvae. The expression level of BmTA in BmNPV resistant strain A35 was higher than that in susceptible strain P50. After BmNPV infection, the expression of BmTA increased in both strains from 24 to 72 h. Virus amplification analysis showed that the relative levels of VP39 in B. mori larvae and BmN cells infected with the appropriate concentration of recombinant-BmTA-treated BmNPV were significantly lower than in the control groups. Moreover, overexpression of BmTA in BmN cells significantly inhibited the amplification of BmNPV. Taken together, the results of this study indicated that BmTA possessed anti-BmNPV activity in B. mori, which broadens the horizon for virus-resistant breeding of silkworms.

Proteomic Identification of Bombyx mori Organelles Using the Engineered Ascorbate Peroxidase APEX and Development of Silkworm Organelle Proteome Database (SilkOrganPDB)

Silkworm Bombyx mori is an economically important insect and a lepidopteran model. Organelle proteome is vital to understanding gene functions; however, it remains to be identified in silkworm. Here, using the engineered ascorbate peroxidase APEX, we constructed transgenic B. mori embryo cells (BmE) expressing APEX-NLS, COX4-APEX, APEX-Rev, and APEX-KDEL in nucleus, mitochondrial matrix (MM), cytosol, and endoplasmic reticulum (ER), and isolated the biotin-labeled proteins using streptavidin-affinity purification, respectively. The isolated proteins were determined using LC-MS/MS and annotated by searching B. mori genomes downloaded from GenBank, SilkBase, SilkDB 2.0, and SilkDB 3.0, resulting in 842, 495, 311, and 445 organelle proteins identified, respectively. We mapped the 296 MM proteins annotated in the GenBank data to mitochondrial protein databases of the fly, human, and mouse, and found that 140 (47%) proteins are homologous to 80 fly proteins, and 65 (22%) proteins match to 31 and 29 human and mouse proteins, respectively. Protein orthology was predicted in multiple insects using OrthoMCL, producing 460 families containing 839 proteins we identified. Out of 460 families, 363 were highly conserved and found in all insects, leaving only three proteins without orthology in other insects, indicating that the identified proteins are highly conserved and probably play important roles in insects. A gene ontology enrichment analysis by clusterProfiler revealed that the nucleus proteins significantly enriched in cellular component terms of nucleus and nucleolus, the MM proteins markedly enriched in molecular function terms of nucleotide binding, and the cytosol proteins mainly enriched in biological process terms of small molecule metabolism. To facilitate the usage and analysis of our data, we developed an open-access database, Silkworm Organelle Proteome Database (SilkOrganPDB), which provides multiple modules for searching, browsing, downloading, and analyzing these proteins, including BLAST, HMMER, Organelle Proteins, Protein Locations, Sequences, Gene Ontology, Homologs, and Phylogeny. In summary, our work revealed the protein composition of silkworm BmE organelles and provided a database resource helpful for understanding the functions and evolution of these proteins.

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.

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.

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.

Plasma Membrane Proteomic Profile Discovers Macrophage-capping Protein Related to Latent HIV-1

BACKGROUND

Due to the persistence of latent HIV-infected cellular reservoirs, HIV virus can not be eradicated completely.

OBJECTIVE

To identify proteins related to HIV latency, we performed a subcellular proteomic study in HIV latent cell lines.

METHODS

An established HIV-1 latent cell model (J-Lat Tat-GFP Clone A7 cells, A7 cells) and its parental cell line (Jurkat cells) were used. The plasma membrane (PM) fraction from cultured cells was enriched through aqueous two-phase partition. PM proteins were extracted and then separated using two-dimensional electrophoresis (2DE). Differentially expressed proteins were identified by mass spectrometry, and verified by western blotting.

RESULTS

Thirteen non-redundant proteins were identified to be differentially expressed in the A7 PM fraction compared to those in the Jurkat PM. Eight had a PM location through Gene Ontology (GO) analysis. A differential protein network of CAPG-ACTR3-CD3D was detected to have interactions with HIV Vpr, Tat, gp160, etc. through STRING software analysis. One of the differential proteins (Macrophage-capping protein (CAPG)) was verified by western blotting to be down- regulated in two cell lines and HIV resting CD4+ T cells negatively selected from patients.

CONCLUSION

We identified 13 proteins in A7 compared to Jurkat cells. CAPG may be a potential biomarker related to HIV latency.

Comparative Proteomic Analysis Reveals Immune Competence in Hemolymph of Bombyx mori Pupa Parasitized by Silkworm Maggot Exorista sorbillans

The silkworm maggot, Exoristasorbillans, is a well-known larval endoparasitoid of the silkworm Bombyxmori that causes considerable damage to the silkworm cocoon crop. To gain insights into the response mechanism of the silkworm at the protein level, we applied a comparative proteomic approach to investigate proteomic differences in the hemolymph of the female silkworm pupae parasitized by E. sorbillans. In total, 50 differentially expressed proteins (DEPs) were successfully identified, of which 36 proteins were upregulated and 14 proteins were downregulated in response to parasitoid infection. These proteins are mainly involved in disease, energy metabolism, signaling pathways, and amino acid metabolism. Eight innate immune proteins were distinctly upregulated to resist maggot parasitism. Apoptosis-related proteins of cathepsin B and 14-3-3 zeta were significantly downregulated in E. sorbillans-parasitized silkworm pupae; their downregulation induces apoptosis. Quantitative PCR was used to further verify gene transcription of five DEPs, and the results are consistent at the transcriptional and proteomic levels. This was the first report on identification of possible proteins from the E. bombycis-parasitized silkworms at the late stage of parasitism, which contributes to furthering our understanding of the response mechanism of silkworms to parasitism and dipteran parasitoid biology.

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.

Differences in gene expression in whitefly associated with CYSDV-infected and virus-free melon, and comparison with expression in whiteflies fed on ToCV- and TYLCV-infected tomato

Background

Cucurbit yellow stunting disorder virus (CYSDV; genus Crinivirus, Closteroviridae) is transmitted in a semipersistent manner by the whitefly, Bemisia tabaci, and is efficiently transmitted by the widely prevalent B. tabaci cryptic species, MEAM1. In this study, we compared transcriptome profiles of B. tabaci MEAM1, after 24 h, 72 h and 7 days of acquisition feeding on melon plants infected with CYSDV (CYSDV-whiteflies) with those fed on virus-free melon, using RNA-Seq technology. We also compared transcriptome profiles with whiteflies fed on tomato plants separately infected with Tomato chlorosis virus (ToCV), a crinivirus closely related to CYSDV, and Tomato yellow leaf curl virus (TYLCV), a member of the genus Begomovirus, which has a distinctly different mode of transmission and their respective virus-free controls, to find common gene expression changes among viruliferous whiteflies feeding on different host plants infected with distinct (TYLCV) and related (CYSDV and ToCV) viruses.

Results

A total of 275 differentially expressed genes (DEGs) were identified in CYSDV-whiteflies, with 3 DEGs at 24 h, 221 DEGs at 72 h, and 51 DEGs at 7 days of virus acquisition. Changes in genes encoding orphan genes (54 genes), phosphatidylethanolamine-binding proteins (PEBP) (20 genes), and AAA-ATPase domain containing proteins (10 genes) were associated with the 72 h time point. Several more orphan genes (20 genes) were differentially expressed at 7 days. A total of 59 common DEGs were found between CYSDV-whiteflies and ToCV-whiteflies, which included 20 orphan genes and 6 lysosomal genes. A comparison of DEGs across the three different virus-host systems revealed 14 common DEGs, among which, eight showed similar and significant up-regulation in CYSDV-whiteflies at 72 h and TYLCV-whiteflies at 24 h, while down-regulation of the same genes was observed in ToCV-whiteflies at 72 h.

Conclusions

Dynamic gene expression changes occurred in CYSDV-whiteflies after 72 h feeding, with decreased gene expression changes associated with 7 days of CYSDV acquisition. Similarities in gene expression changes among CYSDV-whiteflies, ToCV-whiteflies and TYLCV-whiteflies suggest the possible involvement of common genes or pathways for virus acquisition and transmission by whiteflies, even for viruses with distinctly different modes of transmission.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5999-0) contains supplementary material, which is available to authorized users.

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.

A 1H NMR based study of hemolymph metabonomics in different resistant silkworms, Bombyx mori (Lepidotera), after BmNPV inoculation

Bombyx mori nucleopolyhedrovirus (BmNPV) is a primary silkworm pathogen, and the molecular mechanism of silkworm defense to BmNPV infection is still unclear. Herein, comparative metabolomics was adopted to analyze the variations in the hemolymph metabolites of different resistant silkworm strains following BmNPV inoculation using a ¹H NMR method. Trehalose, as an instant source of energy, plays a crucial role in the response to pathogen infections in insects. The level of trehalose was persistently upregulated in the hemolymph of the resistant silkworm strain YeA following infection with BmNPV, compared to that of the susceptible strain YeB, indicating that trehalose metabolism plays a vital role in the response to BmNPV infection. The significant upregulation of TCA cycle relevant metabolites, including malate, fumarate, citrate, succinate, and α-ketoglutarate, was identified at 0 h, 12 h, 48 h, and 96 h post-infection in YeA hemolymph, whereas a significant upregulation in YeB hemolymph was only detected at an early stage of infection (0 h-24 h). The expression level of selected key metabolic enzymes, determined using RT-qPCR, validated the differences in trehalose and TCA cycle relevant metabolite levels. The variations in branched-chain amino acid (BCAA) pathway relevant metabolites in resistant silkworm strains following BmNPV infection showed a regular undulation at different times after infection. A significant accumulation of phenylalanine and tyrosine was observed in YeA following BmNPV infection compared to YeB. The glycolysis and gluconeogenesis pathways showed a relatively low activity in YeA following BmNPV infection. Moreover, the levels of other metabolites related to fat metabolism, transamination, energy metabolism, and glycometabolism, such as glycine, threonine, glutamine, and glutamate, were unstable in the two silkworm strains following BmNPV infection. Thus, our study provides an overview of the metabolic response of the silkworm in response to BmNPV infection, which lays the foundation for clarifying the mechanism of silkworm resistance to BmNPV infection.

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.

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.

Molecular and biochemical responses in the midgut of the silkworm, Bombyx mori, infected with Nosema bombycis

BACKGROUND

Microsporidia are a group of eukaryotic intracellular parasites that infect almost all vertebrates and invertebrates. However, there is little information available of how microsporidia obtain nutrients and energy from host cells. The purpose of this study was to investigate the energy and material requirements of Nosema bombycis for the invasion procedure through analyzing the global variation of the gene expression, protein abundance, fatty acids level and ATP flux induced by the microsporidia N. bombycis infection in the midgut of the silkworm Bombyx mori.

METHODS

A suppression subtractive hybridization (SSH) and quantitative real-time PCR (qPCR) analysis were performed to identify the genes upregulated in the midgut of B. mori 48 h following N. bombycis infection. Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were used to annotate and summarize the differentially expressed genes, according to the categories 'molecular function', 'cellular component' and 'biological process'. To evaluate the nutrition material and energy costs in B.mori infected by N. bombycis, biochemical analysis was performed to determine the variation of protein abundance, fatty acid levels and ATP flux with or without the microsporidia N. bombycis infection in the midgut of the silkworm B. mori.

RESULTS

A total of 744 clones were obtained, 288 clones were randomly selected for sequencing, and 110 unigenes were generated. Amongst these, 49.21%, 30.16% and 14.29% genes were involved in 19 molecular functions, 19 biological processes and nine cellular components, respectively. A total of 11 oxidative phosphorylation- and eight proton-coupled ATP synthesis-related genes were upregulated. Seven protein degradation-, three fat degradation-related genes were upregulated, and no genes related to the de novo synthesis of amino acids and fatty acids were significantly upregulated. The data from the biochemical analysis showed the contents of total protein and ATP of B. mori midgut tissues decreased significantly, whereas the fatty acid content did not significantly change after four days of N. bombycis infection. Microsporidia N. bombycis infection upregulated the expression level of genes involved in host ATP synthesis, protein and fat degradation, which eventually causes the obvious decline of protein content and ATP synthesis in the host midgut, whereas the fatty acids content did not change significantly.

CONCLUSIONS

This study suggested to some extent that N. bombycis invasion can activate the host protein degradation and accelerate the production of host ATP. Microsporidia of N. bombycis show preference for proteins rather than fatty acids from the host to ensure the material preparation required by their parasitic life-cycle. Requirements of N. bombycis for energy were also mainly dependent on the host ATP production. This study provides a new data that may help our understanding of the molecular mechanisms of obtaining energy and nutrients from the host by the microsporidium N. bombycis.

Immune function of a Rab-related protein by modulating the JAK-STAT signaling pathway in the silkworm, Bombyx mori

The Rab-family GTPases mainly regulate intracellular vesicle transport, and play important roles in the innate immune response in invertebrates. However, the function and signal transduction of Rab proteins in immune reactions remain unclear in silkworms. In this study, we analyzed a Rab-related protein of silkworm Bombyx mori (BmRABRP) by raising antibodies against its bacterially expressed recombinant form. Tissue distribution analysis showed that BmRABRP mRNA and protein were high expressed in the Malpighian tubule and fat body, respectively. However, among the different stages, only the fourth instar larvae and pupae showed significant BmRABRP levels. After challenge with four pathogenic microorganisms (Escherichia coli, BmNPV, Beauveria bassiana, Micrococcus luteus), the expression of BmRABRP mRNA in the fat body was significantly upregulated. In contrast, the BmRABRP protein was significantly upregulated after infection with BmNPV, while it was downregulated by E. coli, B. bassiana, and M. luteus. A specific dsRNA was used to explore the immune function and relationship between BmRABRP and the JAK-STAT signaling pathway. After BmRABRP gene interference, significant reduction in the number of nodules and increased mortality suggested that BmRABRP plays an important role in silkworm's response to bacterial challenge. In addition, four key genes (BmHOP, BmSTAT, BmSOCS2, and BmSOCS6) of the JAK-STAT signaling pathway showed significantly altered expressions after BmRABRP silencing. BmHOP and BmSOCS6 expressions were significantly decreased, while BmSTAT and BmSOCS2 were significantly upregulated. Our results suggested that BmRABRP is involved in the innate immune response against pathogenic microorganisms through the JAK-STAT signaling pathway in silkworm.

Systematic identification of mitochondrial lysine succinylome in silkworm (Bombyx mori) midgut during the larval gluttonous stage

Lysine succinylation is a newly identified protein post-translational modification (PTM) of lysine residues. Increasing evidences demonstrate that this modification is prevalent in mitochondria and regulates many vital cellular processes, especially metabolism. Here, we determined the succinylome of the silkworm (Bombyx mori) midgut mitochondria during the larval gluttonous stage (the fifth instar) using succinylated peptides enrichment coupled with nano HPLC/MS/MS. A total of 1884 lysine succinylation sites on 373 mitochondrial proteins were identified. The bioinformatic analysis reveal that succinylated proteins are significantly enriched in central metabolic processes and mitochondrial protein synthesis. Several apoptosis and detoxification related enzymes or proteins are succinylated. The findings suggest the crucial role of lysine succinylation in silkworm midgut metabolism and resistance. Our data provide a rich resource for further analysis of lysine succinylation in silkworm.

SIGNIFICANCE

Insect midgut is the vital tissue for nutrient metabolism and also for xenobiotic metabolism. There is a growing body of knowledge on regulation of midgut function at the gene or protein levels in silkworm, however, the regulation at post-translation modification level remains largely unknown. We provide a first global analysis of the mitochondrial lysine succinylome in silkworm midgut. A total of 1884 lysine succinylation sites on 373 mitochondrial proteins were identified. Bioinformatics results suggest an important role of this modification in regulating metabolism and mitochondrial protein synthesis. Our data greatly expand the catalog of lysine succinylation substrates and sites in insects, and represents an important resource for understanding the physiological function of lysine succinylation in insect midgut.