Variability in the cadherin gene in an Ostrinia nubilalis strain selected for Cry1Ab resistance

Recent progress on the interaction between insects and Bacillus thuringiensis crops

Extensive use of chemical pesticides poses a great threat to the environment and food safety. The discovery of Bacillus thuringiensis (Bt) toxins with effective insecticidal activity against pests and the development of transgenic technology of plants opened a new era of pest control. Transgenic Bt crops, including maize, cotton and soya bean, have now been produced and commercialized to protect against about 30 major coleopteran and lepidopteran pests, greatly benefiting the environment and the economy. However, with the long-term cultivation of Bt crops, some target pests have gradually developed resistance. Numerous studies have indicated that mutations in genes for toxins activation, toxin-binding and insect immunization are important sources in Bt resistance. An in-depth exploration of the corresponding Bt-resistance mechanisms will aid in the design of new strategies to prevent and control pests. Future research will focus on Bt crops expressing new genes and multiple genes to control a broader range of pests as part of an integrated pest management programme.

This article is part of the theme issue ‘Biotic signalling sheds light on smart pest management’.

An Overview of Mechanisms of Cry Toxin Resistance in Lepidopteran Insects

Arthropods have the capacity to evolve resistance to insecticides and insecticidal traits in genetically modified crops. Resistance development among Lepidoptera is a common phenomenon, and a repertoire of resistance mechanisms to various Cry toxins have been identified from laboratory, greenhouse, and field studies in this insect order. Elucidation of such resistance mechanisms is crucial for developing IRM (insect resistance management) strategies to ensure sustainable use of genetically modified crops. This mini review provides a comprehensive overview of mechanisms of resistance that have been reported for lepidopteran pests. This study demonstrated that resistance mechanisms are highly complex, and the most common mechanism of resistance is altered binding sites. It is yet to be established whether all these altered binding sites are regulated by an MAPK signaling pathway, which might suggest a universal mechanism of resistance in lepidopterans.

Field-Evolved Mode 1 Resistance of the Fall Armyworm to Transgenic Cry1Fa-Expressing Corn Associated with Reduced Cry1Fa Toxin Binding and Midgut Alkaline Phosphatase Expression

Insecticidal protein genes from the bacterium Bacillus thuringiensis (Bt) are expressed by transgenic Bt crops (Bt crops) for effective and environmentally safe pest control. The development of resistance to these insecticidal proteins is considered the most serious threat to the sustainability of Bt crops. Resistance in fall armyworm (Spodoptera frugiperda) populations from Puerto Rico to transgenic corn producing the Cry1Fa insecticidal protein resulted, for the first time in the United States, in practical resistance, and Bt corn was withdrawn from the local market. In this study, we used a field-collected Cry1Fa corn-resistant strain (456) of S. frugiperda to identify the mechanism responsible for field-evolved resistance. Binding assays detected reduced Cry1Fa, Cry1Ab, and Cry1Ac but not Cry1Ca toxin binding to midgut brush border membrane vesicles (BBMV) from the larvae of strain 456 compared to that from the larvae of a susceptible (Ben) strain. This binding phenotype is descriptive of the mode 1 type of resistance to Bt toxins. A comparison of the transcript levels for putative Cry1 toxin receptor genes identified a significant downregulation (>90%) of a membrane-bound alkaline phosphatase (ALP), which translated to reduced ALP protein levels and a 75% reduction in ALP activity in BBMV from 456 compared to that of Ben larvae. We cloned and heterologously expressed this ALP from susceptible S. frugiperda larvae and demonstrated that it specifically binds with Cry1Fa toxin. This study provides a thorough mechanistic description of field-evolved resistance to a transgenic Bt crop and supports an association between resistance and reduced Cry1Fa toxin binding and levels of a putative Cry1Fa toxin receptor, ALP, in the midguts of S. frugiperda larvae.

A single type of cadherin is involved in Bacillus thuringiensis toxicity in Plutella xylostella

Cadherins have been described as one the main functional receptors for the toxins of the entomopathogenic bacterium, Bacillus thuringiensis (Bt). With the availability of the whole genome of Plutella xylostella, different types of cadherins have been annotated. In this study we focused on determining those members of the cadherin-related proteins that potentially play a role in the mode of action of Bt toxins. For this, we mined the genome of P. xylostella to identify these putative cadherins. The genome screening revealed 52 genes that were annotated as cadherin or cadherin-like genes. Further analysis revealed that six of these putative cadherins had three motifs common to all Bt-related cadherins: a signal peptide, cadherin repeats and a transmembrane domain. From the six selected cadherins, only P. xylostella cadherin 1 (PxCad1) was expressed in the larval midgut and only the silencing of this gene by RNA interference (double-stranded RNA feeding) reduce toxicity and binding to the midgut of the Cry1Ac type toxin from Bt. These results indicate that from the whole set of cadherin-related genes identified in P. xylostella, only PxCad1 is associated with the Cry1Ac mode of action.

Generation of a Transcriptome in a Model Lepidopteran Pest, Heliothis virescens, Using Multiple Sequencing Strategies for Profiling Midgut Gene Expression

Heliothine pests such as the tobacco budworm, Heliothis virescens (F.), pose a significant threat to production of a variety of crops and ornamental plants and are models for developmental and physiological studies. The efforts to develop new control measures for H. virescens, as well as its use as a relevant biological model, are hampered by a lack of molecular resources. The present work demonstrates the utility of next-generation sequencing technologies for rapid molecular resource generation from this species for which lacks a sequenced genome. In order to amass a de novo transcriptome for this moth, transcript sequences generated from Illumina, Roche 454, and Sanger sequencing platforms were merged into a single de novo transcriptome assembly. This pooling strategy allowed a thorough sampling of transcripts produced under diverse environmental conditions, developmental stages, tissues, and infections with entomopathogens used for biological control, to provide the most complete transcriptome to date for this species. Over 138 million reads from the three platforms were assembled into the final set of 63,648 contigs. Of these, 29,978 had significant BLAST scores indicating orthologous relationships to transcripts of other insect species, with the top-hit species being the monarch butterfly (Danaus plexippus) and silkworm (Bombyx mori). Among identified H. virescens orthologs were immune effectors, signal transduction pathways, olfactory receptors, hormone biosynthetic pathways, peptide hormones and their receptors, digestive enzymes, and insecticide resistance enzymes. As an example, we demonstrate the utility of this transcriptomic resource to study gene expression profiling of larval midguts and detect transcripts of putative Bacillus thuringiensis (Bt) Cry toxin receptors. The substantial molecular resources described in this study will facilitate development of H. virescens as a relevant biological model for functional genomics and for new biological experimentation needed to develop efficient control efforts for this and related Noctuid pest moths.

MAPK signaling pathway alters expression of midgut ALP and ABCC genes and causes resistance to Bacillus thuringiensis Cry1Ac toxin in diamondback moth

Insecticidal crystal toxins derived from the soil bacterium Bacillus thuringiensis (Bt) are widely used as biopesticide sprays or expressed in transgenic crops to control insect pests. However, large-scale use of Bt has led to field-evolved resistance in several lepidopteran pests. Resistance to Bt Cry1Ac toxin in the diamondback moth, Plutella xylostella (L.), was previously mapped to a multigenic resistance locus (BtR-1). Here, we assembled the 3.15 Mb BtR-1 locus and found high-level resistance to Cry1Ac and Bt biopesticide in four independent P. xylostella strains were all associated with differential expression of a midgut membrane-bound alkaline phosphatase (ALP) outside this locus and a suite of ATP-binding cassette transporter subfamily C (ABCC) genes inside this locus. The interplay between these resistance genes is controlled by a previously uncharacterized trans-regulatory mechanism via the mitogen-activated protein kinase (MAPK) signaling pathway. Molecular, biochemical, and functional analyses have established ALP as a functional Cry1Ac receptor. Phenotypic association experiments revealed that the recessive Cry1Ac resistance was tightly linked to down-regulation of ALP, ABCC2 and ABCC3, whereas it was not linked to up-regulation of ABCC1. Silencing of ABCC2 and ABCC3 in susceptible larvae reduced their susceptibility to Cry1Ac but did not affect the expression of ALP, whereas suppression of MAP4K4, a constitutively transcriptionally-activated MAPK upstream gene within the BtR-1 locus, led to a transient recovery of gene expression thereby restoring the susceptibility in resistant larvae. These results highlight a crucial role for ALP and ABCC genes in field-evolved resistance to Cry1Ac and reveal a novel trans-regulatory signaling mechanism responsible for modulating the expression of these pivotal genes in P. xylostella.

Biopesticide and transgenic crops based on Bacillus thuringiensis (Bt) Cry toxins are widely used worldwide, yet the development of field resistance seriously threatens their sustainability. Unraveling these resistance mechanisms are of great importance for delaying insect field resistance evolution. The diamondback moth was the first insect to evolve field resistance to Bt biopesticides and it is an excellent model for the study of Bt resistance mechanisms. In this work, we present strong empirical evidence supporting that (1) field-evolved resistance to Bt in P. xylostella is tightly associated with differential expression of a membrane-bound alkaline phosphatase (ALP) and a suite of ATP-binding cassette transporter subfamily C (ABCC) genes, and (2) a constitutively transcriptionally-activated upstream gene (MAP4K4) in the MAPK signaling pathway is responsible for this trans-regulatory signaling mechanism. These findings identify key resistance genes and provide the first comprehensive mechanistic description responsible for the field-evolved Bt resistance in P. xylostella. Given that expression alterations of multiple receptor genes result in Bt resistance in many other insects, it can now be tested to determine whether the previously unidentified trans-regulatory mechanism characterized in this study is also involved in these cases.

The midgut cadherin-like gene is not associated with resistance to Bacillus thuringiensis toxin Cry1Ac in Plutella xylostella (L.)

The Gram-positive bacterium Bacillus thuringiensis (Bt) produces Cry toxins that have been used to control important agricultural pests. Evolution of resistance in target pests threatens the effectiveness of these toxins when used either in sprayed biopesticides or in Bt transgenic crops. Although alterations of the midgut cadherin-like receptor can lead to Bt Cry toxin resistance in many insects, whether the cadherin gene is involved in Cry1Ac resistance of Plutella xylostella (L.) remains unclear. Here, we present experimental evidence that resistance to Cry1Ac or Bt var. kurstaki (Btk) in P. xylostella is not due to alterations of the cadherin gene. The bona fide P. xylostella cadherin cDNA sequence was cloned and analyzed, and comparisons of the cadherin cDNA sequence among susceptible and resistant P. xylostella strains confirmed that Cry1Ac resistance was independent of mutations in this gene. In addition, real-time quantitative PCR (qPCR) indicated that cadherin transcript levels did not significantly differ among susceptible and resistant P. xylostella strains. RNA interference (RNAi)-mediated suppression of cadherin gene expression did not affect larval susceptibility to Cry1Ac toxin. Furthermore, genetic linkage assays using four cadherin gDNA allelic biomarkers confirmed that the cadherin gene is not linked to resistance against Cry1Ac in P. xylostella. Taken together, our findings demonstrate that Cry1Ac resistance of P. xylostella is independent of the cadherin gene.

Downregulation and mutation of a Cadherin gene associated with Cry1Ac resistance in the Asian Corn Borer, Ostrinia furnacalis (Guenée)

Development of resistance in target pests is a major threat to long-term use of transgenic crops expressing Bacillus thuringiensis (Bt) Cry toxins. To manage and/or delay the evolution of resistance in target insects through the implementation of effective strategies, it is essential to understand the basis of resistance. One of the most important mechanisms of insect resistance to Bt crops is the alteration of the interactions between Cry toxins and their receptors in the midgut. A Cry1Ac-selected strain of Asian corn borer (ACB), Ostrinia furnacalis, a key pest of maize in China, evolved three mutant alleles of a cadherin-like protein (OfCAD) (MPR-r1, MPR-r2 and MPR-r3), which mapped within the toxin-binding region (TBR). Each of the three mutant alleles possessed two or three amino acid substitutions in this region, especially Thr¹⁴⁵⁷→Ser. In highly resistant larvae (ACB-Ac200), MPR-r2 had a 26-amino acid residue deletion in the TBR, which resulted in reduced binding of Cry1Ac compared to the MPR from the susceptible strain, suggesting that the number of amino acid deletions influences the level of resistance. Furthermore, downregulation of OfCAD gene (ofcad) transcription was observed in the Cry1Ac resistant strain, ACB-Ac24, suggesting that Cry1Ac resistance in ACB is associated with the downregulation of the transcript levels of the cadherin-like protein gene. The OfCAD identified from ACB exhibited a high degree of similarity to other members of the cadherin super-family in lepidopteran species.

Alternative splicing and highly variable cadherin transcripts associated with field-evolved resistance of pink bollworm to bt cotton in India

Evolution of resistance by insect pests can reduce the benefits of insecticidal proteins from Bacillus thuringiensis (Bt) that are used extensively in sprays and transgenic crops. Despite considerable knowledge of the genes conferring insect resistance to Bt toxins in laboratory-selected strains and in field populations exposed to Bt sprays, understanding of the genetic basis of field-evolved resistance to Bt crops remains limited. In particular, previous work has not identified the genes conferring resistance in any cases where field-evolved resistance has reduced the efficacy of a Bt crop. Here we report that mutations in a gene encoding a cadherin protein that binds Bt toxin Cry1Ac are associated with field-evolved resistance of pink bollworm (Pectinophora gossypiella) in India to Cry1Ac produced by transgenic cotton. We conducted laboratory bioassays that confirmed previously reported resistance to Cry1Ac in pink bollworm from the state of Gujarat, where Bt cotton producing Cry1Ac has been grown extensively. Analysis of DNA from 436 pink bollworm from seven populations in India detected none of the four cadherin resistance alleles previously reported to be linked with resistance to Cry1Ac in laboratory-selected strains of pink bollworm from Arizona. However, DNA sequencing of pink bollworm derived from resistant and susceptible field populations in India revealed eight novel, severely disrupted cadherin alleles associated with resistance to Cry1Ac. For these eight alleles, analysis of complementary DNA (cDNA) revealed a total of 19 transcript isoforms, each containing a premature stop codon, a deletion of at least 99 base pairs, or both. Seven of the eight disrupted alleles each produced two or more different transcript isoforms, which implicates alternative splicing of messenger RNA (mRNA). This represents the first example of alternative splicing associated with field-evolved resistance that reduced the efficacy of a Bt crop.

Comparative binding of Cry1Ab and Cry1F Bacillus thuringiensis toxins to brush border membrane proteins from Ostrinia nubilalis, Ostrinia furnacalis and Diatraea saccharalis (Lepidoptera: Crambidae) midgut tissue

The European (Ostrinia nubilalis Hübner) and Asian corn borers (Ostrinia furnacalis Guenée) are closely related and display similar sensitivity to Cry1 toxins. In this study, we compared the binding patterns of Cry1Ab and Cry1F toxins between both Ostrinia spp., as well as the expression of putative cadherin- and aminopeptidase-N (APN)-like protein receptors. Additionally, cDNA sequences of these putative toxin receptors from both Ostrinia species were compared. Ligand blots for both species indicated a similar binding pattern for Cry1Ab with the strongest immunoreactive band at 260 kDa in both species. In addition, similar expression of the putative cadherin- and APN-like protein receptors were observed at 260 and 135 kDa, respectively. A high degree of similarity (98% amino acid sequence identity) of cDNA sequences for both putative receptor sequences was observed. The Cry1F ligand blot revealed that O. furnacalis and O. nubilalis BBMV exhibited slightly different binding patterns, with strong binding to putative proteins at 150 and 140 kDa, respectively. Both proteins appeared to also bind Cry1Ab, although the signal intensity was much reduced with Cry1Ab. O. furnacalis showed an additional but weaker band at 210 kDa relative to the 150 kDa band. Diatraea saccharalis (Fabricius), which was used as an outgroup species, exhibited different binding patterns than either Ostrinia species, with both Cry1Ab and Cry1F toxins binding to a 210 kDa protein. These results support the previous experiments indicating that O. nubilalis and O. furnacalis share similar patterns of susceptibility to Cry toxins.

A cadherin-like protein influences Bacillus thuringiensis Cry1Ab toxicity in the oriental armyworm, Mythimna separata

Cadherins comprise a family of calcium-dependent cell adhesion proteins that act in cell-cell interactions. Cadherin-like proteins (CADs) in midguts of some insects act as receptors that bind some of the toxins produced by the Bacillus thuringiensis (Bt). We cloned a CAD gene associated with larval midguts prepared from Mythimna separata. The full-length cDNA (MsCAD1, GenBank Accession No. JF951432) is 5642 bp, with an open reading frame encoding a 1757 amino acid and characteristics typical of insect CADs. Expression of MsCAD1 is predominantly in midgut tissue, with highest expression in the 3rd- to 6th-instars and lowest in newly hatched larvae. Knocking-down MsCAD1 decreased Cry1Ab susceptibility, indicated by reduced developmental time, increased larval weight and reduced larval mortality. We expressed MsCAD1 in E. coli and recovered the recombinant protein, rMsCAD1, which binds Cry1Ab toxin. Truncation analysis and binding experiments revealed that a contiguous 209-aa, located in CR11 and CR12, is the minimal Cry1Ab binding region. These results demonstrate that MsCAD1 is associated with Cry1Ab toxicity and is one of the Cry1Ab receptors in this insect. The significance of this work lies in identifying MsCAD1 as a Cry1Ab receptor, which helps understand the mechanism of Cry1Ab toxicity and of potential resistance to Bt in M. separata.

Quantitative genetic analysis of Cry1Ab tolerance in Ostrinia nubilalis Spanish populations

Tolerance to Bacillus thuringiensis Cry1Ab toxin in Spanish Ostrinia nubilalis populations was analyzed by quantitative genetic techniques, using isolines established from field-derived insects. F1 offspring was tested for susceptibility to trypsin activated Cry1Ab using a concentration that caused a mean larval mortality of 87% (±17% SD). The progeny of the most tolerant isolines (that had shown mortalities lower than 60%) was crossed to obtain the F2 generation that was exposed to the same Cry1Ab concentration. A clear reduction in mortality (62±17% SD) was observed. The upper limit for heritability was estimated to range between 0.82 and 0.90, suggesting that a high part of phenotypic variation in tolerance to Cry1Ab was attributable to genetic differences. An estimate of the minimum number of segregating factors indicated that the loci involved in tolerance to Cry1Ab were at least two. The role of the cadherin gene, which is a B. thuringiensis resistance gene in Lepidoptera, was assessed in the most tolerant isolines by using an EPIC-PCR marker specifically developed for this study. Association between cadherin and tolerance was obtained in one tolerant isoline; however it could be not confirmed by segregation analysis in the F2 progeny because F2 offspring was not viable. Our results indicate that the tolerance trait is common in Spanish field populations. Quantitative genetic techniques may be helpful for estimating the influence of genetic factors to Cry1Ab tolerance in O. nubilalis.

Down regulation of a gene for cadherin, but not alkaline phosphatase, associated with Cry1Ab resistance in the sugarcane borer Diatraea saccharalis

The sugarcane borer, Diatraea saccharalis, is a major target pest of transgenic corn expressing Bacillus thuringiensis (Bt) proteins (i.e., Cry1Ab) in South America and the mid-southern region of the United States. Evolution of insecticide resistance in such target pests is a major threat to the durability of transgenic Bt crops. Understanding the pests' resistance mechanisms will facilitate development of effective strategies for delaying or countering resistance. Alterations in expression of cadherin- and alkaline phosphatase (ALP) have been associated with Bt resistance in several species of pest insects. In this study, neither the activity nor gene regulation of ALP was associated with Cry1Ab resistance in D. saccharalis. Total ALP enzymatic activity was similar between Cry1Ab-susceptible (Cry1Ab-SS) and -resistant (Cry1Ab-RR) strains of D. saccharalis. In addition, expression levels of three ALP genes were also similar between Cry1Ab-SS and -RR, and cDNA sequences did not differ between susceptible and resistant larvae. In contrast, altered expression of a midgut cadherin (DsCAD1) was associated with the Cry1Ab resistance. Whereas cDNA sequences of DsCAD1 were identical between the two strains, the transcript abundance of DsCAD1 was significantly lower in Cry1Ab-RR. To verify the involvement of DsCAD1 in susceptibility to Cry1Ab, RNA interference (RNAi) was employed to knock-down DsCAD1 expression in the susceptible larvae. Down-regulation of DsCAD1 expression by RNAi was functionally correlated with a decrease in Cry1Ab susceptibility. These results suggest that down-regulation of DsCAD1 is associated with resistance to Cry1Ab in D. saccharalis.

Quantitative real-time PCR with SYBR Green detection to assess gene duplication in insects: study of gene dosage in Drosophila melanogaster (Diptera) and in Ostrinia nubilalis (Lepidoptera)

BACKGROUND

The accurate determination of the number of copies of a gene in the genome (gene dosage) is essential for a number of genetic analyses. Quantitative real time PCR (qPCR) with TaqMan detection has shown advantages over traditional Southern-blot and FISH techniques, however the high costs of the required labeled probes is an important limitation of this method. qPCR with SYBR Green I detection is a simple and inexpensive alternative, but it has never been applied to the determination of the copy number of low copy number genes in organisms with high allelic variability (as some insects), where a very small margin of error is essential.

FINDINGS

We have tested the suitability of the qPCR with SYBR Green I detection methodology for the detection of low copy number genes in two insects: the genetically well characterized Drosophila melanogaster (Diptera) and the poor genetically characterized Ostrinia nubilalis (Lepidoptera). The system was applied to determine the copy number of: (1) the O. nubilalis cadherin gene, involved in the mode of action of Bacillus thuringiensis toxins, which showed indirect evidence of duplication, and (2) the D. melanogaster BarH1 and BarH2 genes, located within the Bar region of the X chromosome, to clearly determine whether they both are covered by the tandem duplication in the classical Bar (B1) mutant. Our results showed that the O. nubilalis cadherin gene is an autosomal single copy gene and that BarH1, but not BarH2, is duplicated in the Drosophila B1 mutant.

CONCLUSIONS

This work shows that qPCR with SYBR Green I detection can be specific and accurate enough to distinguish between one and two gene copies per haploid genome of genes with high allelic variability. The technique is sensitive enough to give reliable results with a minimum amount of sample (DNA from individual thoraxes) and to detect gene duplications in tandem.

Candidate genes revealed by a genome scan for mosquito resistance to a bacterial insecticide: sequence and gene expression variations

Background

Genome scans are becoming an increasingly popular approach to study the genetic basis of adaptation and speciation, but on their own, they are often helpless at identifying the specific gene(s) or mutation(s) targeted by selection. This shortcoming is hopefully bound to disappear in the near future, thanks to the wealth of new genomic resources that are currently being developed for many species. In this article, we provide a foretaste of this exciting new era by conducting a genome scan in the mosquito Aedes aegypti with the aim to look for candidate genes involved in resistance to Bacillus thuringiensis subsp. israelensis (Bti) insecticidal toxins.

Results

The genome of a Bti-resistant and a Bti-susceptible strains was surveyed using about 500 MITE-based molecular markers, and the loci showing the highest inter-strain genetic differentiation were sequenced and mapped on the Aedes aegypti genome sequence. Several good candidate genes for Bti-resistance were identified in the vicinity of these highly differentiated markers. Two of them, coding for a cadherin and a leucine aminopeptidase, were further examined at the sequence and gene expression levels. In the resistant strain, the cadherin gene displayed patterns of nucleotide polymorphisms consistent with the action of positive selection (e.g. an excess of high compared to intermediate frequency mutations), as well as a significant under-expression compared to the susceptible strain.

Conclusion

Both sequence and gene expression analyses agree to suggest a role for positive selection in the evolution of this cadherin gene in the resistant strain. However, it is unlikely that resistance to Bti is conferred by this gene alone, and further investigation will be needed to characterize other genes significantly associated with Bti resistance in Ae. aegypti. Beyond these results, this article illustrates how genome scans can build on the body of new genomic information (here, full genome sequence and MITE characterization) to finally hold their promises and help pinpoint candidate genes for adaptation and speciation.