Transcriptome analysis of the immune response of silkworm at the early stage of Bombyx mori bidensovirus infection

Suppression of Adipokinetic hormones enhances Ascovirus HvAV-3h killing speed in Helicoverpa armigera

The emergence of pesticide resistance in Helicoverpa armigera and other pests represents a challenge, necessitating continued development of innovative insecticides. Ascoviruses are a potential insecticide for H. armigera. These double-stranded DNA viruses cause cell death and lethal disease in infected larvae. However, this disease is chronic and prolongs the lifespan of the infected larvae by several weeks. The slow killing speed limits the success of these pathogens in the insecticidal market. Here, we demonstrate that the dsRNA silencing of H. armigera neuropeptide Adipokinetic Hormones (HaAKH) accelerates the killing speed of Heliothis virescens ascovirus 3h (HvAV-3h) against the third-instar larvae of H. armigera. The LT50 was reduced to 3.96 days in larvae with silenced HaAKH1 and 4.7 days in those with silenced HaAKH3. Moreover, the histopathological examinations revealed the destruction of the host's fat body and epidermal tissue shrinkage after HaAKH silencing during HvAV-3h infection. Examining detoxification and antioxidant enzyme activity in HvAV-3h infected larvae showed reduced detoxification mechanisms after HaAKH gene silencing. Furthermore, the silencing of HaAKH resulted in an overall reduction in the fold changes of proline dehydrogenase. In conclusion, this study demonstrates that the ascovirus killing speed can be accelerated by interfering with the host neuropeptide-related gene expression. Moreover, the silencing of H. armigera HaAKH1 and HaAKH3 decreased the antiviral immunity against HvAV-3h.

Identification of long noncoding RNAs of silkworm at the early stage of Bombyx mori bidensovirus infection

Bombyx mori bidensovirus (BmBDV) is one of the most important pathogens of silkworm. It mainly infects midgut cells of silkworm and causes losses to the sericulture industry. Long noncoding RNAs (lncRNAs) have been reported to play an important role in the regulation of antiviral immune response in silkworm. To explore whether lncRNAs are involved in BmBDV infection and immune response of silkworm, we performed a comparative transcriptome analysis to identify the lncRNAs and mRNAs between the BmBDV infected and noninfected silkworm larvae at the early stage. A total of 16,069 genes and 974 candidate lncRNAs were identified, among which 142 messenger RNA (mRNAs) and four lncRNAs were differentially expressed (DE). Target gene prediction revealed that 142 DEmRNAs were coexpressed with four DElncRNAs, suggesting that the expression of mRNA is mainly affected through trans-regulation activities. A regulatory network of DElncRNAs and DEmRNAs was constructed, showing that many genes targeted by different DElncRNAs are involved in metabolism and immunity, which implies that these genes and lncRNAs play an important role in the replication of BmBDV. Our results will help us to improve our understanding of lncRNA-mediated regulatory roles in BmBDV infection, providing a new perspective for further exploring the interaction between host and BmBDV.

Transcriptome analysis reveals that knocking out BmNPV iap2 induces apoptosis by inhibiting the oxidative phosphorylation pathway

Apoptosis is essential for the normal growth, development, and immunity defense of living organisms, and its function and mechanisms have been intensively studied. When viral infection occurs, apoptosis is triggered, causing programmed death of the infected cells. Meanwhile, viruses have also evolved countermeasures to inhibit apoptosis in host cells. We previously constructed a transgenic silkworm line with significantly improved resistance to Bombyx mori nucleopolyhedrovirus (BmNPV) by knocking out the BmNPV inhibitor of apoptosis 2 (iap2) gene. However, the mechanism of how IAP2 induces apoptosis still needs to be further investigated. Here, the transcriptomes of Cas9(-)/sgiap2 (-) and Cas9(+)/sgiap2(+) strains were analyzed at 48 h after BmNPV infection, and a total of 709 differential genes were obtained. A KEGG analysis revealed that the differentially expressed genes were enriched in the oxidative phosphorylation, proteasome, and ribosome pathways. In the oxidative phosphorylation pathway, 41 differentially expressed genes were downregulated, and 12 of these genes were verified by qRT-PCR. More importantly, the knockout of BmNPV iap2 led to the inhibition of the oxidative phosphorylation pathway, followed by activated oxidative stress triggered apoptosis, thereby inhibiting the replication of BmNPV in vitro and vivo. The results provide a basis for the analysis of the initiation of apoptosis that can inhibit virus proliferation, and the study presents new ideas for the subsequent creation of resistant material.

Comparative transcriptome analysis reveals regional specialization of gene expression in larval silkworm (Bombyx mori) midgut

Insect midgut plays a central role in food digestion and nutrition absorption. Larval silkworm midgut could be divided into 3 distinct regions based on their morphological colors. However, it remains rudimentary of regional gene expression and physiological function in larval silkworm midgut. Through transcriptome sequencing of 3 midgut compartments, a comprehensive analysis of gene expression atlas along the anterior-posterior axis was conducted. Posterior midgut was found transcriptionally divergent from anterior and middle midgut. Differentially expressed gene analysis revealed the regional specialization of digestive enzyme production, transmembrane transport, chitin metabolism, and hormone regulation in different midgut regions. In addition, gene subsets of pan-midgut and region-specific transcription factors (TFs) along the length of midgut were also identified. The results suggested that homeobox TFs might play an essential role in transcriptional variations across the midgut. Altogether, our study provided the first fundamental resource to investigate physiological function and regulation mechanism in larval midgut compartmentalization.

Assembly-free rapid differential gene expression analysis in non-model organisms using DNA-protein alignment

BACKGROUND

RNA-seq is being increasingly adopted for gene expression studies in a panoply of non-model organisms, with applications spanning the fields of agriculture, aquaculture, ecology, and environment. For organisms that lack a well-annotated reference genome or transcriptome, a conventional RNA-seq data analysis workflow requires constructing a de-novo transcriptome assembly and annotating it against a high-confidence protein database. The assembly serves as a reference for read mapping, and the annotation is necessary for functional analysis of genes found to be differentially expressed. However, assembly is computationally expensive. It is also prone to errors that impact expression analysis, especially since sequencing depth is typically much lower for expression studies than for transcript discovery.

RESULTS

We propose a shortcut, in which we obtain counts for differential expression analysis by directly aligning RNA-seq reads to the high-confidence proteome that would have been otherwise used for annotation. By avoiding assembly, we drastically cut down computational costs - the running time on a typical dataset improves from the order of tens of hours to under half an hour, and the memory requirement is reduced from the order of tens of Gbytes to tens of Mbytes. We show through experiments on simulated and real data that our pipeline not only reduces computational costs, but has higher sensitivity and precision than a typical assembly-based pipeline. A Snakemake implementation of our workflow is available at: https://bitbucket.org/project_samar/samar .

CONCLUSIONS

The flip side of RNA-seq becoming accessible to even modestly resourced labs has been that the time, labor, and infrastructure cost of bioinformatics analysis has become a bottleneck. Assembly is one such resource-hungry process, and we show here that it can be avoided for quick and easy, yet more sensitive and precise, differential gene expression analysis in non-model organisms.

Genome-wide identification of lipases in silkworm (Bombyx mori) and their spatio-temporal expression in larval midgut

Lipases play crucial roles in food digestion by degrading dietary lipids into free fatty acids and glycerols. The domesticated silkworm (Bombyx mori) has been widely used as an important Lepidopteran model for decades. However, little is known about the lipase gene family in the silkworm, especially their hydrolytic activities as digestive enzymes. In this study, a total of 38 lipase genes were identified in the silkworm genome. Phylogenetic analysis indicated that they were divided into three major groups. Twelve lipases were confirmed to be expressed in the midgut at both transcriptional and translational levels. They were grouped into the same gene cluster, suggesting that they could have similar physiological functions. Quantitative real-time PCR (qRT-PCR) analyses indicated that lipases were mainly expressed in anterior and middle midgut regions, and their expression levels varied greatly along the length of midgut. A majority of lipases were down-regulated in the midgut when larvae stopped feeding. However, a unique lipase gene (Bmlip10583) showed low expression level during feeding stage, but it was significantly up-regulated during the larvae-pupae transition. These results demonstrated that expression of silkworm lipases was spatially and temporally regulated in the midgut during larval development. Taken together, our results provide a fundamental research of the lipase gene family in the silkworm.

High-Quality de novo Chromosome-Level Genome Assembly of a Single Bombyx mori With BmNPV Resistance by a Combination of PacBio Long-Read Sequencing, Illumina Short-Read Sequencing, and Hi-C Sequencing

The reference genomes of Bombyx mori (B. mori), Silkworm Knowledge-based database (SilkDB) and SilkBase, have served as the gold standard for nearly two decades. Their use has fundamentally shaped model organisms and accelerated relevant studies on lepidoptera. However, the current reference genomes of B. mori do not accurately represent the full set of genes for any single strain. As new genome-wide sequencing technologies have emerged and the cost of high-throughput sequencing technology has fallen, it is now possible for standard laboratories to perform full-genome assembly for specific strains. Here we present a high-quality de novo chromosome-level genome assembly of a single B. mori with nuclear polyhedrosis virus (BmNPV) resistance through the integration of PacBio long-read sequencing, Illumina short-read sequencing, and Hi-C sequencing. In addition, regular bioinformatics analyses, such as gene family, phylogenetic, and divergence analyses, were performed. The sample was from our unique B. mori species (NB), which has strong inborn resistance to BmNPV. Our genome assembly showed good collinearity with SilkDB and SilkBase and particular regions. To the best of our knowledge, this is the first genome assembly with BmNPV resistance, which should be a more accurate insect model for resistance studies.

Profiling of MicroRNAs in Midguts of Plutella xylostella Provides Novel Insights Into the Bacillus thuringiensis Resistance

The diamondback moth (DBM), Plutella xylostella, one of the most destructive lepidopteran pests worldwide, has developed field resistance to Bacillus thuringiensis (Bt) Cry toxins. Although miRNAs have been reported to be involved in insect resistance to multiple insecticides, our understanding of their roles in mediating Bt resistance is limited. In this study, we constructed small RNA libraries from midguts of the Cry1Ac-resistant (Cry1S1000) strain and the Cry1Ac-susceptible strain (G88) using a high-throughput sequencing analysis. A total of 437 (76 known and 361 novel miRNAs) were identified, among which 178 miRNAs were classified into 91 miRNA families. Transcripts per million analysis revealed 12 differentially expressed miRNAs between the Cry1S1000 and G88 strains. Specifically, nine miRNAs were down-regulated and three up-regulated in the Cry1S1000 strain compared to the G88 strain. Next, we predicted the potential target genes of these differentially expressed miRNAs and carried out GO and KEGG pathway analyses. We found that the cellular process, metabolism process, membrane and the catalytic activity were the most enriched GO terms and the Hippo, MAPK signaling pathway might be involved in Bt resistance of DBM. In addition, the expression patterns of these miRNAs and their target genes were determined by RT-qPCR, showing that partial miRNAs negatively while others positively correlate with their corresponding target genes. Subsequently, novel-miR-240, one of the differentially expressed miRNAs with inverse correlation with its target genes, was confirmed to interact with Px017590 and Px007885 using dual luciferase reporter assays. Our study highlights the characteristics of differentially expressed miRNAs in midguts of the Cry1S1000 and G88 strains, paving the way for further investigation of miRNA roles in mediating Bt resistance.

Host Response against Virus Infection in an Insect: Bidensovirus Infection Effect on Silkworm (Bombyx mori)

Silk cocoons obtained from silkworms are the primary source of commercial silk, making the silkworm an economically important insect. However, the silk industry suffers significant losses due to various virus infections. Bombyx mori bidensovirus (BmBDV) is one of the pathogens that cause flacherie disease in silkworms. Most silkworm strains die after BmBDV infection. However, certain silkworm strains show resistance to the virus, which is determined by a single recessive gene, nsd-2. The +^nsd-2 gene (allele of nsd-2; the susceptibility gene) encodes a putative amino acid transporter expressed only in the insect’s midgut, where BmBDV can infect, suggesting that this membrane protein may function as a receptor for BmBDV. Interestingly, the expression analysis revealed no changes in the +^nsd-2 gene expression levels in virus-uninfected silkworms, whereas the gene expression drastically decreased in the virus-infected silkworm. This condition indicates that the host factor’s expression, the putative virus receptor, is affected by BmBDV infection. It has recently been reported that the expression levels of some host genes encoding cuticle, antioxidant, and immune response-related proteins were significantly regulated by BmBDV infection. In this review, we discuss the host response against virus infection based on our knowledge and long-term research experience in this field.

Insights Into the Antiviral Pathways of the Silkworm Bombyx mori

The lepidopteran model silkworm, Bombyx mori, is an important economic insect. Viruses cause serious economic losses in sericulture; thus, the economic importance of these viruses heightens the need to understand the antiviral pathways of silkworm to develop antiviral strategies. Insect innate immunity pathways play a critical role in the outcome of infection. The RNA interference (RNAi), NF-kB-mediated, immune deficiency (Imd), and stimulator of interferon gene (STING) pathways, and Janus kinase/signal transducer and activator of transcription (JAK/STAT) pathway are the major antiviral defense mechanisms, and these have been shown to play important roles in the antiviral immunity of silkworms. In contrast, viruses can modulate the prophenol oxidase (PPO), phosphatidylinositol 3-kinase (PI3K)/protein kinase B (Akt), and the extracellular signal-regulated kinase (ERK) signaling pathways of the host to elevate their proliferation in silkworms. In this review, we present an overview of the current understanding of the main immune pathways in response to viruses and the signaling pathways modulated by viruses in silkworms. Elucidation of these pathways involved in the antiviral mechanism of silkworms furnishes a theoretical basis for the enhancement of virus resistance in economic insects, such as upregulating antiviral immune pathways through transgenic overexpression, RNAi of virus genes, and targeting these virus-modulated pathways by gene editing or inhibitors.

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.