Patterns of Gene Expression in Western Corn Rootworm (Diabrotica virgifera virgifera) Neonates, Challenged with Cry34Ab1, Cry35Ab1 and Cry34/35Ab1, Based on Next-Generation Sequencing

A draft Diabrotica virgifera virgifera genome: insights into control and host plant adaption by a major maize pest insect

BACKGROUND

Adaptations by arthropod pests to host plant defenses of crops determine their impacts on agricultural production. The larval host range of western corn rootworm, Diabrotica virgifera virgifera (Coleoptera: Chrysomelidae), is restricted to maize and a few grasses. Resistance of D. v. virgifera to crop rotation practices and multiple insecticides contributes to its status as the most damaging pest of cultivated maize in North America and Europe. The extent to which adaptations by this pest contributes to host plant specialization remains unknown.

RESULTS

A 2.42 Gb draft D. v. virgifera genome, Dvir_v2.0, was assembled from short shotgun reads and scaffolded using long-insert mate-pair, transcriptome and linked read data. K-mer analysis predicted a repeat content of ≥ 61.5%. Ortholog assignments for Dvir_2.0 RefSeq models predict a greater number of species-specific gene duplications, including expansions in ATP binding cassette transporter and chemosensory gene families, than in other Coleoptera. A majority of annotated D. v. virgifera cytochrome P450s belong to CYP4, 6, and 9 clades. A total of 5,404 transcripts were differentially-expressed between D. v. virgifera larvae fed maize roots compared to alternative host (Miscanthus), a marginal host (Panicum virgatum), a poor host (Sorghum bicolor) and starvation treatments; Among differentially-expressed transcripts, 1,908 were shared across treatments and the least number were between Miscanthus compared to maize. Differentially-expressed transcripts were enriched for putative spliceosome, proteosome, and intracellular transport functions. General stress pathway functions were unique and enriched among up-regulated transcripts in marginal host, poor host, and starvation responses compared to responses on primary (maize) and alternate hosts.

CONCLUSIONS

Manual annotation of D. v. virgifera Dvir_2.0 RefSeq models predicted expansion of paralogs with gene families putatively involved in insecticide resistance and chemosensory perception. Our study also suggests that adaptations of D. v. virgifera larvae to feeding on an alternate host plant invoke fewer transcriptional changes compared to marginal or poor hosts. The shared up-regulation of stress response pathways between marginal host and poor host, and starvation treatments may reflect nutrient deprivation. This study provides insight into transcriptomic responses of larval feeding on different host plants and resources for genomic research on this economically significant pest of maize.

Towards Understanding the Function of Aegerolysins

Aegerolysins are remarkable proteins. They are distributed over the tree of life, being relatively widespread in bacteria and fungi, but also present in some insects, plants, protozoa, and viruses. Despite their abundance in cells of certain developmental stages and their presence in secretomes, only a few aegerolysins have been studied in detail. Their function, in particular, is intriguing. Here, we summarize previously published findings on the distribution, molecular interactions, and function of these versatile aegerolysins. They have very diverse protein sequences but a common fold. The machine learning approach of the AlphaFold2 algorithm, which incorporates physical and biological knowledge of protein structures and multisequence alignments, provides us new insights into the aegerolysins and their pore-forming partners, complemented by additional genomic support. We hypothesize that aegerolysins are involved in the mechanisms of competitive exclusion in the niche.

Co-Transcriptomic Analysis of the Maize-Western Corn Rootworm Interaction

The Western corn rootworm (WCR; Diabrotica virgifera virgifera) is an economically important belowground pest of maize. Belowground feeding by WCR is damaging because it weakens the roots system, diminishes nutrient uptake, and creates entry points for fungal and bacterial pathogens and increases lodging, all of which can significantly suppress maize yields. Previously, it was demonstrated that belowground herbivory can trigger plant defense responses in the roots and the shoots, thereby impacting intraplant communication. Although several aspects of maize-WCR interactions have been reported, co-transcriptomic remodeling in the plant and insect are yet to be explored. We used a maize genotype, Mp708, that is resistant to a large guild of herbivore pests to study the underlying plant defense signaling network between below and aboveground tissues. We also evaluated WCR compensatory transcriptome responses. Using RNA-seq, we profiled the transcriptome of roots and leaves that interacted with WCR infestation up to 5 days post infestation (dpi). Our results suggest that Mp708 shoots and roots had elevated constitutive and WCR-feeding induced expression of genes related to jasmonic acid and ethylene pathways, respectively, before and after WCR feeding for 1 and 5 days. Similarly, extended feeding by WCR for 5 days in Mp708 roots suppressed many genes involved in the benzoxazinoid pathway, which is a major group of indole-derived secondary metabolites that provides resistance to several insect pests in maize. Furthermore, extended feeding by WCR on Mp708 roots revealed several genes that were downregulated in WCR, which include genes related to proteolysis, neuropeptide signaling pathway, defense response, drug catabolic process, and hormone metabolic process. These findings indicate a dynamic transcriptomic dialog between WCR and WCR-infested maize plants.

Toxicometabolomic profiling of resistant and susceptible western corn rootworm larvae feeding on Bt maize seedlings

The western corn rootworm (WCR), Diabrotica virgifera virgifera LeConte, is the most serious pest of maize (Zea mays L.) in the U.S. Corn Belt and parts of Europe. Transgenic maize hybrids expressing at least one of the four currently available insecticidal toxins from Bacillus thuringiensis (Bt) Berliner, currently the most widely adopted control method in continuous maize, have faltered due to the emergence of resistance. The resistance mechanisms of WCR to Bt toxins are not fully understood. We identified metabolic profiles of susceptible and resistant WCR larvae fed on maize hybrids expressing each of three available Cry3 proteins (eCry3Ab1, mCry3A, and Cry3Bb1) targeting corn rootworms and a control non-Bt maize via an untargeted metabolomics approach. Over 580 unique metabolites found in WCR larvae were classified into different pathways (amino acids, carbohydrates, cofactors and vitamins, energy, lipid, nucleotide, peptide, and xenobiotics). By exploring shifts in WCR larval metabolome exclusively by Bt toxins, several candidate metabolites and metabolic pathways were identified in susceptible and resistant larvae that may be involved in defense against or recovery from Bt ingestion by these larvae. These findings would provide mechanistic insights into altered metabolic pathways associated with the resistance mechanisms of WCR to Bt toxins.

Up-regulation of apoptotic- and cell survival-related gene pathways following exposures of western corn rootworm to B. thuringiensis crystalline pesticidal proteins in transgenic maize roots

BACKGROUND

Resistance of pest insect species to insecticides, including B. thuringiensis (Bt) pesticidal proteins expressed by transgenic plants, is a threat to global food security. Despite the western corn rootworm, Diabrotica virgifera virgifera, being a major pest of maize and having populations showing increasing levels of resistance to hybrids expressing Bt pesticidal proteins, the cell mechanisms leading to mortality are not fully understood.

RESULTS

Twenty unique RNA-seq libraries from the Bt susceptible D. v. virgifera inbred line Ped12, representing all growth stages and a range of different adult and larval exposures, were assembled into a reference transcriptome. Ten-day exposures of Ped12 larvae to transgenic Bt Cry3Bb1 and Gpp34/Tpp35Ab1 maize roots showed significant differential expression of 1055 and 1374 transcripts, respectively, compared to cohorts on non-Bt maize. Among these, 696 were differentially expressed in both Cry3Bb1 and Gpp34/Tpp35Ab1 maize exposures. Differentially-expressed transcripts encoded protein domains putatively involved in detoxification, metabolism, binding, and transport, were, in part, shared among transcripts that changed significantly following exposures to the entomopathogens Heterorhabditis bacteriophora and Metarhizium anisopliae. Differentially expressed transcripts in common between Bt and entomopathogen treatments encode proteins in general stress response pathways, including putative Bt binding receptors from the ATP binding cassette transporter superfamily. Putative caspases, pro- and anti-apoptotic factors, as well as endoplasmic reticulum (ER) stress-response factors were identified among transcripts uniquely up-regulated following exposure to either Bt protein.

CONCLUSIONS

Our study suggests that the up-regulation of genes involved in ER stress management and apoptotic progression may be important in determining cell fate following exposure of susceptible D. v. virgifera larvae to Bt maize roots. This study provides novel insights into insect response to Bt intoxication, and a possible framework for future investigations of resistance mechanisms.

Detection of alternative splicing in western corn rootworm (Diabrotica virgifera virgifera LeConte) in association with eCry3.1Ab resistance using RNA-seq and PacBio Iso-Seq

Alternative splicing is a common feature in eukaryotes that not only increases the transcript diversity, but also has functional consequences. In insects, alternative splicing has been found associated with resistance to pesticides and Bt toxins. Up to date, the alternative splicing in western corn rootworm (Diabrotica virgifera virgifera LeConte) has not been studied. To investigate its alternative splicing pattern and relation to Bt resistance, we carried out single-molecule real-time (SMRT) transcript sequencing and Iso-seq analysis on resistant, eCry3.1Ab-selected and susceptible, unselected, western corn rootworm neonate midguts which fed on seedling maize with and without eCry3.1Ab for 12 and 24 h. We present transcriptome-wide alternative splicing patterns of western corn rootworm midgut in response to feeding on eCry3.1Ab-expressing corn using a comprehensive approach that combines both RNA-seq and SMRT transcript sequencing techniques. The results showed genes in western corn rootworm are highly alternatively spliced, which happens on 67.73% of multi-exon genes. One of the alternative splicing events we identified was a novel peritrophic matrix protein with two alternative splicing isoforms. Analysis of differential exon usage between resistant and susceptible colonies showed that in eCry3.1Ab-resistant western corn rootworm, expression of one isoform was significantly higher than in the susceptible colony, while no significant differences between colonies were observed with the other isoform. Our results provide the first survey of alternative splicing in western corn rootworm and suggest that the observed alternatively spliced isoforms of peritrophic matrix protein may be associated with eCry3.1Ab resistance in western corn rootworm.

Dissecting Out the Molecular Mechanism of Insecticidal Activity of Ostreolysin A6/Pleurotolysin B Complexes on Western Corn Rootworm

Ostreolysin A6 (OlyA6) is a protein produced by the oyster mushroom (Pleurotus ostreatus). It binds to membrane sphingomyelin/cholesterol domains, and together with its protein partner, pleurotolysin B (PlyB), it forms 13-meric transmembrane pore complexes. Further, OlyA6 binds 1000 times more strongly to the insect-specific membrane sphingolipid, ceramide phosphoethanolamine (CPE). In concert with PlyB, OlyA6 has potent and selective insecticidal activity against the western corn rootworm. We analysed the histological alterations of the midgut wall columnar epithelium of western corn rootworm larvae fed with OlyA6/PlyB, which showed vacuolisation of the cell cytoplasm, swelling of the apical cell surface into the gut lumen, and delamination of the basal lamina underlying the epithelium. Additionally, cryo-electron microscopy was used to explore the membrane interactions of the OlyA6/PlyB complex using lipid vesicles composed of artificial lipids containing CPE, and western corn rootworm brush border membrane vesicles. Multimeric transmembrane pores were formed in both vesicle preparations, similar to those described for sphingomyelin/cholesterol membranes. These results strongly suggest that the molecular mechanism of insecticidal action of OlyA6/PlyB arises from specific interactions of OlyA6 with CPE, and the consequent formation of transmembrane pores in the insect midgut.

Functional validation of DvABCB1 as a receptor of Cry3 toxins in western corn rootworm, Diabrotica virgifera virgifera

Western corn rootworm (WCR), Diabrotica virgifera virgifera (Coleoptera: Chrysomelidae), is a serious insect pest in the major corn growing areas of North America and in parts of Europe. WCR populations with resistance to Bacillus thuringiensis (Bt) toxins utilized in commercial transgenic traits have been reported, raising concerns over their continued efficacy in WCR management. Understanding the modes of action of Bt toxins is important for WCR control and resistance management. Although different classes of proteins have been identified as Bt receptors for lepidopteran insects, identification of receptors in WCR has been limited with no reports of functional validation. Our results demonstrate that heterologous expression of DvABCB1 in Sf9 and HEK293 cells conferred sensitivity to the cytotoxic effects of Cry3A toxins. The result was further validated using knockdown of DvABCB1 by RNAi which rendered WCR larvae insensitive to a Cry3A toxin. However, silencing of DvABCB2 which is highly homologous to DvABCB1 at the amino acid level, did not reduce the sensitivity of WCR larvae to a Cry3A toxin. Furthermore, our functional studies corroborate different mode-of-actions for other insecticidal proteins including Cry34Ab1/35Ab1, Cry6Aa1, and IPD072Aa against WCR. Finally, reduced expression and alternatively spliced transcripts of DvABCB1 were identified in a mCry3A-resistant strain of WCR. Our results provide the first clear demonstration of a functional receptor in the molecular mechanism of Cry3A toxicity in WCR and confirmed its role in the mechanism of resistance in a mCry3A resistant strain of WCR.

Aegerolysins from the fungal genus Pleurotus - Bioinsecticidal proteins with multiple potential applications

The aegerolysin proteins ostreolysin A6, pleurotolysin A2 and erylysin A are produced by mushrooms of the genus Pleurotus. These aegerolysins can interact specifically with sphingolipid-enriched membranes. In particular, they strongly bind insect cells and to artificial lipid membranes that contain physiologically relevant concentrations of the main invertebrate-specific sphingolipid, ceramide phosphoethanolamine. Moreover, the aegerolysins permeabilise these membranes when combined with their protein partner pleurotolysin B, which contains a membrane-attack-complex/perforin domain. These aegerolysin/ pleurotolysin B complexes show strong and selective toxicity towards western corn rootworm larvae and adults and Colorado potato beetle larvae. Their insecticidal activities arise through aegerolysin binding to ceramide phosphoethanolamine in the insect midgut. This mode of membrane binding is different from those described for similar aegerolysin-based complexes of bacterial origin (e.g., Cry34Ab1/Cry35Ab1), or other Bacillus thuringiensis proteinaceous crystal toxins, which associate with protein receptors. The ability of Pleurotus aegerolysins to specifically interact with sphingolipid-enriched domains in mammalian cells can be further exploited to visualize lipid rafts in living cells, and to treat certain types of tumours and metabolic disorders. Finally, these proteins can strongly enhance fruiting initiation of P. ostreatus even when applied externally. In this review, we summarise the current knowledge of the potential biotechnological and biomedical applications of the Pleurotus aegerolysins, either alone or when complexed with pleurotolysin B, with special emphasis on their bioinsecticidal effects.

Midgut de novo transcriptome analysis and gene expression profiling of Spodoptera exigua larvae exposed with sublethal concentrations of Cry1Ca protein

Spodoptera exigua (Hübner) is a polyphagous pest on agricultural crops, whose control is based mainly on the application of chemical insecticides. Bacillus thuringiensis (Bt) is one of the most important biological agents that have been successfully applied as a biological control, and Cry1Ca protein is considered to be active against S. exigua. Therefore, to understand the response of S. exigua to Cry1Ca protein, high-throughput sequencing was used to analyse the S. exigua larval midgut after treatment with sublethal concentrations of Cry1Ca protein. Transcriptome data showed that a total of 98,571 unigenes with an N50 value of 1135 bp and a mean length of 653 bp were obtained. Furthermore, 2962 differentially expressed genes (DEGs) were identified after Cry1Ca challenge, including 1508 up-regulated and 1454 down-regulated unigenes. Among these DEGs, detoxification (CYP, CarE, and GST) and Bt resistance (ALP, APN, and ABC transporter)-related genes were differentially expressed in the midgut of S. exigua after Cry1Ca treatment. However, most DEGs of protective enzymes were down-regulated, while most DEGs related with serine protease and REPAT were up-regulated. Furthermore, almost all DEGs related to the immune signaling pathway, antimicrobial protein, and lysozyme were up-regulated by Cry1Ca treatment. These results indicated that the detoxification enzyme, protective enzymes, Bt resistance-related genes, serine protease, REPAT, and the immune response might have been involved in the response of S. exigua to Cry1Ca protein. In summary, analysis of the transcriptomal expression of genes involved in Cry1Ca protein against S. exigua provided potential clues for elucidating the host response processes and defensive mechanisms underlying Cry1Ca toxicity.

Characterization of the Spectrum of Insecticidal Activity for IPD072Aa: A Protein Derived from Pseudomonas chlororaphis with Activity Against Diabrotica virgifera virgifera (Coleoptera: Chrysomelidae)

Western corn rootworm (Diabrotica virgifera virgifera LeConte) presents significant pest management challenges for farmers in both North America and Europe. IPD072Aa, a protein derived from Pseudomonas chlororaphis, has previously been shown to have activity against western corn rootworm. In the current study, the spectrum of activity of IPD072Aa was evaluated in controlled laboratory diet bioassays. IPD072Aa was fed at high concentrations in subchronic or chronic bioassays to 11 different insect species, representing 4 families within Coleoptera, and an additional 4 species representing four families of Lepidoptera. No adverse effects were noted in the Lepidoptera species. Within the order Coleoptera, western corn rootworm was the most sensitive species tested. A range of responses was observed within each of the four families of Coleoptera evaluated that included either no-observed effects or reduced growth, developmental delays, and/or reduced survival. These data will help inform the environmental risk assessment of genetically modified plants that express the IPD072Aa protein for western corn rootworm control.

Differential gene expression in response to eCry3.1Ab ingestion in an unselected and eCry3.1Ab-selected western corn rootworm (Diabrotica virgifera virgifera LeConte) population

Diabrotica virgifera virgifera LeConte, the western corn rootworm (WCR) is one of the most destructive pests in the U.S. Corn Belt. Transgenic maize lines expressing various Cry toxins from Bacillus thuringiensis have been adopted as a management strategy. However, resistance to many Bt toxins has occurred. To investigate the mechanisms of Bt resistance we carried out RNA-seq using Illumina sequencing technology on resistant, eCry3.1Ab-selected and susceptible, unselected, whole WCR neonates which fed on seedling maize with and without eCry3.1Ab for 12 and 24 hours. In a parallel experiment RNA-seq experiments were conducted when only the midgut of neonate WCR was evaluated from the same treatments. After de novo transcriptome assembly we identified differentially expressed genes (DEGs). Results from the assemblies and annotation indicate that WCR neonates from the eCry3.1Ab-selected resistant colony expressed a small number of up and down-regulated genes following Bt intoxication. In contrast, unselected susceptible WCR neonates expressed a large number of up and down-regulated transcripts in response to intoxication. Annotation and pathway analysis of DEGs between susceptible and resistant whole WCR and their midgut tissue revealed genes associated with cell membrane, immune response, detoxification, and potential Bt receptors which are likely related to eCry3.1Ab resistance. This research provides a framework to study the toxicology of Bt toxins and mechanism of resistance in WCR, an economically important coleopteran pest species.