Identification of immunity-related genes in Plutella xylostella in response to fungal peptide destruxin A: RNA-Seq and DGE analysis

A PLA2 deletion mutant using CRISPR/Cas9 coupled to RNASeq reveals insect immune genes associated with eicosanoid signaling

Eicosanoids mediate insect immune responses and synthesized by the catalytic activity of phospholipase A2 (PLA2). A uniquely encoded secretory PLA2 (sPLA2) is associated with immune responses of a lepidopteran insect, Spodoptera exigua. Its deletion mutant was generated using a CRISPR/Cas9 genome editing technology. Both wild and mutant lines were then immune-challenged, and the resulting transcripts were compared with their naïve transcripts by RNASeq using the Illumina-HiSeq platform. In total, 12,878 unigenes were further analyzed by differentially expressed gene tools. Over 69% of the expressed genes in S. exigua larvae are modulated in their expression levels by eicosanoids, recorded from CRISPR/Cas9 mutagenesis against an eicosanoid-synthetic gene, Se-sPLA2. Further, about 36% of the immune-associated genes are controlled by the eicosanoids in S. exigua. Indeed, the deletion mutant suffered significant immunosuppression in both cellular and humoral responses in response to bacterial challenge as well as severely reduced developmental and reproductive potentials.

Comparative transcriptome analysis reveals disruption of Plutella xylostella immune system by fungal peptide cyclosporin C

The diamondback moth (Plutella xylostella), a major threat to crucifers across the globe, has developed resistance against the majority of insecticides enhancing the need for alternate control measures against this pest. Recently cyclosporin C, a secondary metabolite produced by the insect pathogenic fungus Purpeocillium lilacinum, has been reported to induce lethal and sub-lethal effects against P. xylostella. To date, little is known about the molecular mechanisms of interaction between cyclosporin C and P. xylostella immune systems. This study reports the transcriptome-based immune response of P. xylostella to cyclosprin C treatment. Our results showed differential expression of 322, 97, and 504 differentially expressed genes (DEGS) in P. xylostella treated with cyclosporin C compared to control 24, 48, and 72 h post-treatment, respectively. Thirteen DEGs were commonly expressed at different time intervals in P. xylostella larvae treated with cyclosporin C compared to control. Cyclosporin C treatment induced the down-regulated expression of majority of immune-related genes related to pattern recognition responses, signal modulation, Toll and IMD pathways, antimicrobial peptides and antioxidant responses confirming the ability to suppress immune response of P. xylostella. These results will further improve our knowledge of the infection mechanism and complex biochemical processes involved in interaction between cyclosporin C and insect immune systems.

Immune-Related Gene Profiles and Differential Expression in the Grey Garden Slug Deroceras reticulatum Infected with the Parasitic Nematode Phasmarhabditis hermaphrodita

The grey garden slug (Deroceras reticulatum), a common terrestrial slug native to Europe with a global distribution including North America, is commonly considered the most severe slug pest in agriculture. The nematode Phasmarhabditis hermaphrodita, which has been used in the U.K. and Europe as a commercial biocontrol agent since 1994, has also recently been collected in Oregon and California and has long been considered a candidate biocontrol agent for slug management in the U.S. In this study, we report differential gene expressions in nematode-infected slugs using RNA-seq to identify slug immune-related genes against nematodes. Comparison of gene expression levels between the whole bodies of a nematode-infected slug (N-S) and an uninfected control slug (C-S) revealed that there were a total of 39,380 regulated unigenes, of which 3084 (3%) were upregulated and 6761 (6%) were downregulated at greater than 2-fold change (FC > 2) in the nematode-infected slug. To further investigate the biological functions of differentially expressed genes (DEGs), gene ontology (GO) and functional enrichment analysis were performed to map the DEGs to terms in the GO, eukaryotic ortholog groups of proteins (KOG) and Kyoto Encyclopedia of Genes and Genome Pathway (KEGG) databases. Among these DEGs, approximately 228 genes associated with immunity or immune-related pathways were upregulated 2-fold or more in the N-S compared to C-S. These genes include toll, Imd, JNK, scavenger receptors (SCRs), C-type lectins (CTLs), immunoglobulin-like domains, and JAK/STAT63 signaling pathways. From the RNA-seq results, we selected 18 genes and confirmed their expression levels by qRT-PCR. Our findings provide insights into the immune response of slugs during nematode infection. These studies provide fundamental information that will be valuable for the development of new methods of pest slug control using pathogenic nematodes in the field.

Destruxin A inhibits the hemocytin-mediated hemolymph immunity of host insects to facilitate Metarhizium infection

Insects have an effective innate immune system to protect themselves against fungal invasion. Metarhizium employs a toxin-based strategy using a nonribosomal peptide called destruxin A (DA) to counteract the host immune response. However, the mechanism by which DA inhibits insect immunity is still unclear. Here, we identified 48 DA-binding proteins in silkworm hemolymph, with the binding affinity (KD) ranging from 2 to 420 μM. Among these proteins, hemocytin, an important immune factor, was determined to be the strongest DA-binding protein. DA binds to hemocytin and regulates its conformation in a multisite manner. Furthermore, DA exerts a significant inhibitory effect on hemocytin-mediated hemocyte aggregation. By disrupting the interaction between hemocytin, actin A3, and gelsolin, DA prevents the transformation of granules into vesicles in hemocytes. These vesicles are responsible for storing, maturing, and exocytosing hemocytin. Therefore, hemocytin secretion is reduced, and the formation of structures that promote aggregation in outer hemocytes is inhibited.

Immune responses of the Asian onion moth, Acrolepiopsis sapporensis, and their genetic factors from RNA-Seq analysis

A nonmodel insect, Acrolepiopsis sapporensis, has been analyzed in immune responses. The total hemocytes in the fifth instar larvae were 2.33 × 106  cells/mL. These hemocytes comprised at least five different types and different relative ratios: 47% granulocytes, 26% plasmatocytes, 11% oenocytoid, 8% prohemocytes, and 5% spherulocytes. Upon bacterial challenge, some of the hemocytes exhibited typical hemocyte-spreading behaviors, such as focal adhesion, and filopodial and lamellipodial cytoplasmic extensions. The hemocyte behaviors induced cellular immune responses demonstrated by nodule formation. In addition, the plasma collected from the immune-challenged larvae exhibited humoral immune responses by bacterial growth inhibition along with enhanced phenoloxidase enzyme activity. These cellular and humoral immune responses were further analyzed by determining the immune-associated genes from a transcriptome generated by RNA-Seq. A total of about 12 Gb sequences led to about 218,116 contigs, which were predicted to encode about 46,808 genes. Comparative expression analysis showed 8392 uniquely expressed genes in the immune-challenged larvae. Differentially expressed gene (DEG) analysis among the commonly expressed genes indicated that 782 genes were upregulated and 548 genes were downregulated in the expressions after bacterial challenge. These immune-associated genes included pattern recognition receptors, immune mediation/signaling genes, and various immune effectors. Specifically, the genetic components of the Toll, IMD, and JAK/STAT immune signaling pathways were included in the DEG database. These results demonstrate the immune responses of A. sapporensis larvae and suggest the genes associated with the immune responses in this nonmodel insect.

Binding proteins of destruxin A from Metarhizium against insect cell

Destruxin A (DA) is a cyclo-hexadepsipeptidic insecticidal mycotoxin isolated from the entomopathogenic fungi, Metarhizium spp. However, its mode of action is unknown. In this study, we isolated 149 candidate DA-binding proteins by drug affinity response target stability, and determined the interactions of 80 canditates with DA in vitro by surface plasmon resonance. The affinity coefficients (K_D) ranged from 24 to 469 μM. Binding proteins were functionally diverse and included cytoskeletal components and cell motility, protein transcription and translation pathways, ubiquitin dependent protein metabolic processes, nucleus pore entry and exit, and endoplasmic reticulum vesicle transport and etc. Electron microscopy revealed that DA damaged the cytoskeleton and multiple organelles, disrupted cell adhesion and motility, and led to cell death. DA appeared to have a multi-targeted approach to cellular structures and multiple life processes, leading to cell death. The results of this study could provide molecular evidence for the analysis of the insecticidal toxicology of DA and further improve the study of the pathogenic insect mechanism of Metarhizium.

Antibiotics accelerate growth at the expense of immunity

Antibiotics have long been used in the raising of animals for agricultural, industrial or laboratory use. The use of subtherapeutic doses in diets of terrestrial and aquatic animals to promote growth is common and highly debated. Despite their vast application in animal husbandry, knowledge about the mechanisms behind growth promotion is minimal, particularly at the molecular level. Evidence from evolutionary research shows that immunocompetence is resource-limited, and hence expected to trade off with other resource-demanding processes, such as growth. Here, we ask if accelerated growth caused by antibiotics can be explained by genome-wide trade-offs between growth and costly immunocompetence. We explored this idea by injecting broad-spectrum antibiotics into wood tiger moth (Arctia plantaginis) larvae during development. We follow several life-history traits and analyse gene expression (RNA-seq) and bacterial (r16S) profiles. Moths treated with antibiotics show a substantial depletion of bacterial taxa, faster growth rate, a significant downregulation of genes involved in immunity and significant upregulation of growth-related genes. These results suggest that the presence of antibiotics may aid in up-keeping the immune system. Hence, by reducing the resource load of this costly process, bodily resources may be reallocated to other key processes such as growth.

Insecticidal action of the botanical insecticide wilforine on Mythimna separata (Walker) related with the changes of ryanodine receptor expression

The detailed molecular mechanism of wilforine, a novel botanical insecticidal component, remains unclear, except for the knowledge that it affects the calcium signaling pathway. The aim of the current study was to examine the underlying molecular mechanism of wilforine in Mythimna separata (Walker) by transcriptome and RNA interference (RNAi), with chlorantraniliprole as control. RNA sequencing showed that the relative expression of genes related to the calcium signaling pathway and muscle contraction in M. separata treated with wilforine significantly changed and was further validated by qRT-PCR. Interestingly, the expression level of the ryanodine receptor (MsRyR) gene was downregulated by wilforine at relatively high concentrations and long treatment time, contrary to that observed using chlorantraniliprole. Furthermore, a putative MsRyR was cloned using a 16,258-bp contiguous sequence containing a 308-bp 5'-untranslated region and 578-bp 3'-untranslated region by RT-PCR and RACE. The results of the RNAi experiment showed that injection of dsMsRyR significantly reduced MsRyR mRNA levels, and growth and development were inhibited. Importantly, silencing of the MsRyR gene resulted in decreased susceptibility to both wilforine and chlorantraniliprole. Together with the results of our previous studies on toxic symptoms and muscle tissue lesions between wilforine and chlorantraniliprole, we propose that RyR Ca²⁺ release channel dysfunction is closely related with significant lethal mechanisms of wilforine.

sPLA2 behaves like a prophylactic agent and mediates cellular and humoral immune responses in Plutella xylostella

Most immune effectors are inducible to microbial pathogen infection while some are already present to act as prophylactic immunity against as yet unseen infection. This study identified secretory phospholipase A₂ (sPLA₂ ) as a prophylactic factor in diamondback moth, Plutella xylostella. Western blotting using a polyclonal antibody raised against other lepidopteran sPLA₂ reacted specifically with ∼25 kDa protein, which was present at approximately 0.4 mM in the plasma of naïve larvae. Interrogation of P. xylostella transcriptomes revealed an open-reading frame for sPLA₂ (Px-sPLA₂ ), exhibiting high homology with other Group III sPLA₂ s. Px-sPLA₂ was expressed in all developmental stages. In the larval stage, bacterial challenge induced its expression in hemocytes and fat body but not in gut or epidermis. RNA interference (RNAi) suppressed Px-sPLA₂ messenger RNA level and sPLA₂ activity in plasma. An inhibition zone assay showed that Px-sPLA₂ exhibited antibacterial activities against different species, because specific RNAi knockdown impaired the activity. The RNAi treatment also suppressed the cellular immune response assessed by hemocyte nodule formation and humoral immune response assessed by antimicrobial peptide gene expression. Finally, benzylideneacetone (BZA, a specific sPLA₂ inhibitor) treatment inhibited plasma sPLA₂ activity of naive larvae in a dose-dependent manner. An addition of BZA significantly increased the bacterial virulence of an entomopathogen, Bacillus thuringiensis. These results suggest that Px-sPLA₂ is an immune-associated factor of P. xylostella and its relatively high level of concentration in the plasma of naive larvae strongly suggests its role as a prophylactic factor in defending against pathogens at early infection stages.

Role of serine protease inhibitors in insect-host-pathogen interactions

Serine protease inhibitors (serpins), evolutionary old, structurally conserved molecules, are a superfamily of proteins found in almost all living organisms. Serpins are relatively large, typically 350-500 amino acids in length, with three β-sheets and seven to nine α-helices folding into a conserved tertiary structure with a reactive center loop. Serpins perform various physiological functions in insects, including development, digestion, host-pathogen interactions, and innate immune response. In insects, the innate immune system is characterized as the first and major defense system against the invasion of microorganisms. Serine protease cascades play a critical role in the initiation of innate immune responses, such as melanization and the production of antimicrobial peptides, and are strictly and precisely regulated by serpins. Herein, we provide a microreview on the role of serpins in the insect-host-pathogen interactions, emphasizing their role in immune responses, particularly in diamondback moth (Plutella xylostella), highlighting the important discoveries and also the gaps that remain to be explored in future studies.

Effects of Destruxin A on Silkworm's Immunophilins

Destruxin A (DA), a major secondary metabolite of Metarhizium anisopliae, has anti-immunity to insects. However, the detailed mechanism and its interactions with target proteins are elusive. Previously, three immunophilins, peptidyl–prolyl cis–trans isomerase (BmPPI), FK506 binding-protein 45 (BmFKBP45) and BmFKBP59 homologue, were isolated from the silkworm, Bombyx mori Bm12 cell line following treatment with DA, which suggested that these proteins were possible DA-binding proteins. To validate the interaction between DA and the three immunophilins, we performed bio-layer interferometry (BLI) assay, and the results showed that DA has interaction with BmPPI, whose affinity constant value is 1.98 × 10⁻³ M and which has no affinity with FKBP45 and FKBP59 homologue in vitro. Furthermore, we investigated the affinity between DA and human PPI protein (HsPPIA) and the affinity constant (K_D) value is 2.22 × 10⁻³ M. Additionally, we compared the effects of silkworm and human PPI proteins produced by DA and immunosuppressants, cyclosporine A (CsA), and tacrolimus (FK506), by employing I2H (insect two-hybrid) in the SF-9 cell line. The results indicated that in silkworm, the effects created by DA and CsA were stronger than FK506. Furthermore, the effects created by DA in silkworm were stronger than those in humans. This study will offer new thinking to elucidate the molecular mechanism of DA in the immunity system of insects.

BmTudor-sn Is a Binding Protein of Destruxin A in Silkworm Bm12 Cells

Destruxin A (DA), a hexa-cyclodepsipeptidic mycotoxin secreted by the entomopathogenic fungus Metarhizium anisopliae, was reported to have an insecticidal effect and anti-immunity activity. However, its molecular mechanism of action remains unclear. Previously, we isolated several potential DA-affinity (binding) proteins in the Bombyx mori Bm12 cell line. By docking score using MOE2015, we selected three proteins—BmTudor-sn, BmPiwi, and BmAGO2—for further validation. First, using Bio-Layer Interferometry in vitro, we found that BmTudor-sn had an affinity interaction with DA at 125, 250, and 500 µM, while BmPiwi and BmAGO2 had no interaction signal with DA. Second, we employed standard immunoblotting to verify that BmTudor-sn is susceptible to DA, but BmPiwi and BmAGO2 are not. Third, to verify these findings in vivo, we used a target engagement strategy based on shifts in protein thermal stability following ligand binding termed the cellular thermal shift assay and found no thermal stability shift in BmPiwi and BmAGO2, whereas a shift was found for BmTudor-sn. In addition, in BmTudor-sn knockdown Bm12 cells, we observed that cell viability increased under DA treatment. Furthermore, insect two-hybrid system results indicated that the key site involved in DA binding to BmTudor-sn was Leu704. In conclusion, in vivo and in vitro experimental evidence indicated that BmTudor-sn is a binding protein of DA in silkworm Bm12 cells at the 100 µM level, and the key site of this interaction is Leu704. Our results provide new perspectives to aid in elucidating the molecular mechanism of action of DA in insects and developing new biopesticide.

Bacillus thuringiensis Suppresses the Humoral Immune System to Overcome Defense Mechanism of Plutella xylostella

Background: Plutella xylostella has become a notorious pest of cruciferous crops all over the world. Delta-endotoxins of Bacillus thuringiensis are widely used insecticidal proteins for controlling P. xylostella. However, the interaction mechanism of B. thuringiensis with the immune system of P. xylostella, at the genomic level, is still unclear. This study explored the immune response of P. xylostella to B. thuringiensis, at different time intervals, 6 h, 12 h, 18 h, 24 h, and 36 h, by using RNA-Sequencing (RNA-Seq) and RT-qPCR. Results: In total, 167 immunity-related genes were identified and placed into different families, including pattern recognition receptors (PRRs), signal modulators, immune pathways (Toll, IMD, and JAK/STAT), and immune effectors. It is worth mentioning that the analyses of the differentially expressed immunity-related genes revealed that most of the differentially expressed genes (DEGs) (87, 56, 76, 67, and 73 genes) were downregulated in P. xylostella following B. thuringiensis oral infection at 6 h, 12 h, 18 h, 24 h, and 36 h. Interestingly, our RNA-Seq analysis also revealed reduced expression of antimicrobial peptides, that play a vital role in the humoral immune system of P. xylostella. Conclusion: This study demonstrates that B. thuringiensis plays a novel role in controlling P. xylostella, by suppressing the immune system.

Genome-Wide Identification of Destruxin A-Responsive Immunity-Related MicroRNAs in Diamondback Moth, Plutella xylostella

Plutella xylostella, a global key pest, is one of the major lepidopteran pests of cruciferous vegetables owing to its strong ability of resistance development to a wide range of insecticides. Destruxin A, a mycotoxin of the entomopathogenic fungus, Metarhizium anisopliae, has broad-spectrum insecticidal effects and has been used as an alternative control strategy to reduce harmful effects of insecticides. However, microRNA (miRNA)-regulated reactions against destruxin A have not been elucidated yet. Therefore, here, to identify immunity-related miRNAs, we constructed four small RNA libraries from destruxin A-injected larvae of P. xylostella at three different time courses (2, 4, and 6 h) with a control, and sequenced by Illumina. Our results showed that totally 187 known and 44 novel miRNAs were identified in four libraries by bioinformatic analysis. Interestingly, among differentially expressed known miRNAs, some conserved miRNAs, such as miR-263, miR-279, miR-306, miR-2a, and miR-308, predicted to be involved in regulating immunity-related genes, were also identified. Worthy to mention, miR-306 and miR-279 were also listed as common abundantly expressed miRNA in all treatments. The Kyoto Encyclopedia of Genes and Genomes pathway analysis also indicated that differentially expressed miRNAs were involved in several immunity-related signaling pathways, including toll signaling pathway, IMD signaling pathway, JAK-STAT signaling pathway, and cell adhesion molecules signaling pathway. To the best of our knowledge, this is the first comprehensive report of destruxin A-responsive immunity-related miRNAs in P. xylostella. Our findings will improve in understanding the role of destruxin A-responsive miRNAs in the host immune system and would be useful to develop biological control strategies for controlling P. xylostella.

Genome-Wide Profiling of Plutella xylostella Immunity-Related miRNAs after Isaria fumosorosea Infection

The development of resistance by Plutella xylostella to almost all insecticides is of significant concern all over the world. Entomopathogenic fungi such as Isaria fumosorosea have been used as an alternative to insecticides. However, the knowledge of miRNA-regulated reactions against entomopathogenic fungi is still in its infant stage. In the present study, P. xylostella was challenged with I. fumosorosea at four different time points (12, 18, 24, and 36 h) including a control, to build miRNA libraries by Illumina sequencing. The results of differential expression analysis exhibited that 23 miRNAs were differentially expressed, compared to control, in all treatments. It is worth mentioning, of these, some conserved miRNAs such as miR-2, miR-9a, miR-745, miR-7b, and miR-2767, known to play critical roles in host-pathogen interaction, were also identified. Furthermore, differentially expressed miRNAs were validated by RT-qPCR. Our results provide an essential information for further functional studies of the interaction between I. fumosorosea and P. xylostella at the post-transcriptional level.