Epistasis confers resistance to Bt toxin Cry1Ac in the cotton bollworm

Resistance of Lepidopteran Pests to Bacillus thuringiensis Toxins: Evidence of Field and Laboratory Evolved Resistance and Cross-Resistance, Mode of Resistance Inheritance, Fitness Costs, Mechanisms Involved and Management Options

Bacillus thuringiensis (Bt) toxins are potential alternatives to synthetic insecticides for the control of lepidopteran pests. However, the evolution of resistance in some insect pest populations is a threat and can reduce the effectiveness of Bt toxins. In this review, we summarize the results of 161 studies from 20 countries reporting field and laboratory-evolved resistance, cross-resistance, and inheritance, mechanisms, and fitness costs of resistance to different Bt toxins. The studies refer mainly to insects from the United States of America (70), followed by China (31), Brazil (19), India (12), Malaysia (9), Spain (3), and Australia (3). The majority of the studies revealed that most of the pest populations showed susceptibility and a lack of cross-resistance to Bt toxins. Factors that delay resistance include recessive inheritance of resistance, the low initial frequency of resistant alleles, increased fitness costs, abundant refuges of non-Bt, and pyramided Bt crops. The results of field and laboratory resistance, cross-resistance, and inheritance, mechanisms, and fitness cost of resistance are advantageous for predicting the threat of future resistance and making effective strategies to sustain the effectiveness of Bt crops.

Genomic signals of local adaptation across climatically heterogenous habitats in an invasive tropical fruit fly (Bactrocera tryoni)

Local adaptation plays a key role in the successful establishment of pest populations in new environments by enabling them to tolerate novel biotic and abiotic conditions experienced outside their native range. However, the genomic underpinnings of such adaptive responses remain unclear, especially for agriculturally important pests. We investigated population genomic signatures in the tropical/subtropical Queensland fruit fly, Bactrocera tryoni, which has an expanded range encompassing temperate and arid zones in Australia, and tropical zones in the Pacific Islands. Using reduced representation sequencing data from 28 populations, we detected allele frequency shifts associated with the native/invasive status of populations and identified environmental factors that have likely driven population differentiation. We also determined that precipitation, temperature, and geographic variables explain allelic shifts across the distribution range of B. tryoni. We found spatial heterogeneity in signatures of local adaptation across various climatic conditions in invaded areas. Specifically, disjunct invasive populations in the tropical Pacific Islands and arid zones of Australia were characterised by multiple significantly differentiated single nucleotide polymorphisms (SNPs), some of which were associated with genes with well-understood function in environmental stress (e.g., heat and desiccation) response. However, invasive populations in southeast Australian temperate zones showed higher gene flow with the native range and lacked a strong local adaptive signal. These results suggest that population connectivity with the native range has differentially affected local adaptive patterns in different invasive populations. Overall, our findings provide insights into the evolutionary underpinnings of invasion success of an important horticultural pest in climatically distinct environments.

Down-regulation of HaABCC3, potentially mediated by a cis-regulatory mechanism, is involved in resistance to Cry1Ac in the cotton bollworm, Helicoverpa armigera

Evolution of resistance to Cry proteins in multiple pest insects has been threatening the sustainable use of Bacillus thuringiensis (Bt)-transgenic crops. Better understanding about the mechanism of resistance to Cry proteins in insects is needed. Our preliminary study reported that the transcription of HaABCC3 was significantly decreased in a near-isogenic line (LFC2) of a Cry1Ac-resistant strain (LF60) of the global pest Helicoverpa armigera. However, the causality between HaABCC3 downregulation and resistance to Cry1Ac remains to be verified, and the regulatory mechanism underlying the HaABCC3 downregulation is still unclear. In this study, our data showed that both HaABCC3 and HaABCC3 downregulation were genetically linked to resistance to Cry1Ac in LF60. However, no InDels were observed in the coding sequence of HaABCC3 from LF60. Furthermore, F₁ offspring from the cross of LF60 and a HaABCC2/3-knockout mutant exhibited moderate resistance to Cry1Ac toxin; this indicated that the high resistance to Cry1Ac toxin in LF60 may have resulted from multiple genetic factors, including HaABCC2 mis-splicing and HaABCC3 downregulation. Results from luciferase reporter assays showed that promoter activity of HaABCC3 in LF60 was significantly lower than that in the susceptible strain, which indicated that HaABCC3 downregulation was likely mediated by promoter variation. Consistently, multiple variations of the GATA- or FoxA-binding sites in the promoter region of HaABCC3 were identified. Collectively, all results in this study suggested that the downregulation of HaABCC3 observed in the H. armigera LF60 strain, which is resistant to Cry1Ac, may be mediated by a cis-regulatory mechanism.

Generation of Human Domain Antibody Fragments as Potential Insecticidal Agents against Helicoverpa armigera by Cadherin-Based Screening

New insecticidal genes and approaches for pest control are a hot research area. In the present study, we explored a novel strategy for the generation of insecticidal proteins. The midgut cadherin of Helicoverpa armigera (H. armigera) was used as a target to screen materials that have insecticidal activity. After three rounds of panning, the phage-displayed human domain antibody B1F6, which not only binds to the H. armigera cadherin CR9-CR11 but also significantly inhibits Cry1Ac toxins from binding to CR9-CR11, was obtained from a phage-displayed human domain antibody (DAb) library. To better analyze the relevant activity of B1F6, soluble B1F6 protein was expressed by Escherichia coli BL21 (DE3). The cytotoxicity assays demonstrated that soluble B1F6 induced Sf9 cell death when expressing H. armigera cadherin on the cell membrane. The insect bioassay results showed that soluble B1F6 protein (90 μg/cm²) caused 49.5 ± 3.3% H. armigera larvae mortality. The midgut histological results showed that soluble B1F6 caused damage to the midgut epithelium of H. armigera larvae. The present study explored a new strategy and provided a basic material for the generation of new insecticidal materials.

Pyramiding of cry toxins and methanol producing genes to increase insect resistance in cotton

The idea of enhanced methanol production from cell wall by pectin methyl esterase enzymes (PME) combined with expression of cry genes from Bacillus thuringiensis as a strategy to improve insect pest control in cotton is presented. We constructed a cassette containing two cry genes (cry1Fa and Cry32Aa) and two pme genes, one from Arabidopsis thaliana (AtPME), and other from Aspergillus. niger (AnPME) in pCAMBIA1301 plant expression vector using CAMV-35S promoter. This construction was transformed in Eagle-2 cotton variety by using shoot apex-cut Agrobacterium-mediated transformation. Expression of cry genes and pme genes was confirmed by qPCR. Methanol production was measured in control and in the cry and pme transformed plants showing methanol production only in transformed plants, in contrast to the non-transgenic cotton plants. Finally, insect bioassays performed with transgenic plants expressing cry and pme genes showed 100% mortality for Helicoverpa armigera (cotton bollworm) larvae, 70% mortality for Pectinophora gossypiella (pink bollworm) larvae and 95% mortality of Earias fabia, (spotted bollworm) larvae, that was higher than the transgenic plants expressing only cry genes that showed 84%, 49% and 79% mortality, respectively. These results demonstrate that Bt. cry-genes coupled with pme genes are an effective strategy to improve the control of different insect pests.

Evaluating Cross-Resistance to Cry and Vip Toxins in Four Strains of Helicoverpa armigera With Different Genetic Mechanisms of Resistance to Bt Toxin Cry1Ac

Evolution of resistance by pests has diminished the efficacy of transgenic crops producing insecticidal proteins from Bacillus thuringiensis (Bt). In China, where transgenic cotton producing Bt toxin Cry1Ac has been planted since 1997, field control failures have not been reported but the frequency of resistance to Cry1Ac has increased in the cotton bollworm, Helicoverpa armigera. This provides incentive to switch to multi-toxin Bt cotton, which is grown in many other countries. Previous work created four laboratory strains of H. armigera with >100-fold resistance to Cry1Ac, with the genetic basis of resistance known in all but the LF256 strain. Here, we analyzed the genetic basis of resistance in Cry1Ac in LF256 and evaluated cross-resistance of all four strains to three toxins produced by widely planted multi-toxin Bt cotton: Cry1Fa, Cry2Ab, and Vip3Aa. DNA sequencing revealed that LF256 lacked the mutations in three genes (HaTSPAN1, HaABCC2, and HaABCC3) that confer resistance to Cry1Ac in two other strains of H. armigera we analyzed. Together with previous results, the data reported here show that each of the four strains examined has a different genetic basis of resistance to Cry1Ac. Significant positive cross-resistance occurred to Cry1Fa in three of the four strains tested but not to Cry2Ab or Vip3Aa in any strain. Thus, Cry2Ab and Vip3Aa are likely to be especially valuable for increasing the efficacy and durability of Bt cotton against H. armigera populations that have some resistance to Cry1Ac.

Synergistic selection of a Helicoverpa armigera cadherin fragment with Cry1Ac in different cells and insects

The midgut cadherin fragments were extensively studied as Bt synergists in insects, while their synergistic selection modes with Bt toxins in different mechanisms of resistance or insects have never been determined. Here, a soluble Helicoverpa armigera cadherin fragment which corresponds to the Cry1Ac binding region (HaCad-TBR) was expressed in Escherichia coli and its synergism with Cry1Ac toxin in H. armigera and Plutella xylostella larvae as well as Sf9 cells expressing different cadherins was tested. HaCad-TBR exhibited higher synergism factor in P. xylostella larvae (4.84-fold) than in H. armigera larvae (2.45-fold). Among the cells expressing HaCad alleles, HaCad-TBR enhanced the Cry1Ac toxicity only in the cells expressing the mutant lacking the extracellular domain. Moreover, HaCad-TBR had a weak enhancement of Cry1Ac toxicity in Sf9 cells expressing the P. xylostella cadherin. Further researches revealed that the enhancement of toxicity in Sf9 cells was correlated with increased toxin binding. These results suggested that cadherin fragments which have high binding level with Cry1Ac are more likely to enhance toxin toxicity well against the cells or larvae where the cadherin has lower binding level with Cry1Ac, especially in the cases lacking the toxin binding domain.

Shared and Independent Genetic Basis of Resistance to Bt Toxin Cry2Ab in Two Strains of Pink Bollworm

Evolution of pest resistance threatens the benefits of crops genetically engineered to produce insecticidal proteins from Bacillus thuringiensis (Bt). Field populations of the pink bollworm (Pectinophora gossypiella), a global pest of cotton, have evolved practical resistance to transgenic cotton producing Bt toxin Cry2Ab in India, but not in the United States. Previous results show that recessive mutations disrupting an autosomal ATP-binding cassette gene (PgABCA2) are associated with pink bollworm resistance to Cry2Ab in field-selected populations from India and in one lab-selected strain from the United States (Bt4-R2). Here we discovered that an independently derived, lab-selected Cry2Ab-resistant pink bollworm strain from the United States (BX-R) also harbors mutations that disrupt PgABCA2. Premature stop codons introduced by mis-splicing of PgABCA2 pre-mRNA were prevalent in field-selected larvae from India and in both lab-selected strains. The most common mutation in field-selected larvae from India was also detected in both lab-selected strains. Results from interstrain crosses indicate BX-R has at least one additional mechanism of resistance to Cry2Ab that does not involve PgABCA2 and is not completely recessive or autosomal. We conclude that recessive mutations disrupting PgABCA2 are the primary, but not the only, mechanism of resistance to Cry2Ab in pink bollworm.

Identification of a Cry1Fa binding site of cadherin in Plutella xylostella through fragment exchanging and molecular docking methods

Binding to the cadherin in target pests is the primary step in the action mechanism of Cry toxins, but little is known regarding the interaction of Cry1Fa with cadherin. Our previous study suggested that a Plutella xylostella cadherin fragment (PxCad-TBR) can bind to Cry1Fa, while its homologous fragment (HaCad-TBR) in Helicoverpa armigera cannot. In this study, we expressed two cadherin fragments that combine parts of PxCad-TBR and HaCad-TBR in Escherichia coli and tested their binding to the Cry1Fa. The results showed that the fragment containing amino acids T1202-A1341 of P. xylostella cadherin showed binding ability to Cry1Fa. Furthermore, two regions (V1219-E1233 and D1326-F1337) were predicted as hot spot regions that are involved in the interaction of Cry1Fa and PxCad-TBR with computer-aided molecular docking. We then constructed two PxCad-TBR mutations by fragment exchanging based on the molecular docking results and verified the mutations' binding abilities to the Cry1Fa. The results showed that the region that contains amino acids D1326-F1337 was one important binding site to Cry1Fa in P. xylostella cadherin. These results suggested that a combination of computer-aided molecular docking and fragment exchanging is an effective way to locate the key binding sites of Bt toxins in receptors.

A Combinational Approach of Enhanced Methanol Production and Double Bt Genes for Broad Spectrum Insect Resistance in Transgenic Cotton

The prevalence of insect resistance against Bt toxins has led to the idea of enhancing demethylation from cell wall pectin by pectin methylesterase enzyme for overproduction of methanol which is toxic to insects pests. The AtPME and AnPME fragments ligated into pCAMBIA1301 vector were confirmed through restriction digestion with EcoR1 and BamH1. Excision of 3363 bp fragment from 11,850 bp vector confirmed the ligation of both fragments into pCAMBIA1301 vector. Transformation of pectin methylesterase-producing genes, i.e., AtPME and AnPME from Arabidopsis thaliana and Aspergillus niger cloned in plant expression vector pCAMBIA1301 under 35S promoter into cotton variety CEMB-33 harboring two Bt genes Cry1Ac and Cry2A, respectively, was done by using shoot apex-cut Agrobacterium-mediated transformation method. The plantlets were screened on MS medium supplemented with hygromycin on initial basis. Amplification of 412 and 543 bp, respectively, through gene-specific primer has been obtained which confirmed the successful introduction of pCAMBIA AtPME and AnPME genes into cotton variety CEMB 33. Relative expression of AtPME and AnPME genes through real-time PCR determined the expression level of both gene ranges between 3- and 3.5-fold in different transgenic cotton lines along with quantity of methanol ranging from 0.8 to 0.9% of maximum while 0.5% to 0.6% of minimum but no expression was obtained in negative non-transgenic control cotton plant with least quantity of methanol, i.e., 0.1%. Almost 100% mortality was observed in insect bioassay for Helicoverpa armigera on detached leaves bioassay and 63% for Pink Bollworm (Pectinophora gossypiella) on growing transgenic cotton bolls as compared to positive control transgenic cotton with double Bt genes where mortality was found to be 82% for H. armigera and 50% for P. gossypiella while 0% in negative control non-transgenic plants.

Roles of Midgut Cadherin from Two Moths in Different Bacillus thuringiensis Action Mechanisms: Correlation among Toxin Binding, Cellular Toxicity, and Synergism

The midgut cadherin has been described as one of the main functional receptors for Bacillus thuringiensis (Bt) toxins. Plutella xylostella (P. xylostella) and Helicoverpa armigera (H. armigera) are two major target pests of Bt toxins in China, and the roles of their cadherins in the action of Bt toxins have been only partially studied. Here, we expressed the two cadherins in Sf9 cells and their partial extracellular domains in Escherichia coli and tested them for Bt toxin binding, cellular toxicity, and synergism with toxins. Our results suggested that PxCad might function as a Cry1Ac receptor, although it showed lower binding levels to Cry1Ac and reduced cytotoxicity compared with HaCad. PxCad and HaCad are not receptors for Cry2A, Cry1B, Cry1C, and Cry1F toxins, although some of them can bind to the cadherins. The PxCad-TBR exhibits higher enhancement of Cry1Ac and weak enhancement of Cry1F toxicity in P. xylostella larvae, although it is not the receptor of Cry1F.

Newly identified APN splice isoforms suggest novel splicing mechanisms may underlie circRNA circularization in moth

Circular RNA (circRNA) have long been considered by-products of splicing errors, but the coordination of RNA transcription and exon circularization events remains poorly understood. Here, we investigated this question using genes encoding aminopeptidases N (APNs), which are receptors of Bacillus thuringiensis toxins, in the cotton bollworm, Helicoverpa armigera. We cloned and sequenced the cDNA of ten APN genes (HaAPN1-10) located in the same APN gene cluster, and detected 20 and 14 novel splicing isoforms with exon skipping in HaAPN1 and HaAPN3, respectively, whereas no or very few variants were found in the remaining genes. Further study identified 14 and 6 circular RNA (circRNA) in HaAPN1 and HaAPN3, respectively. Neither novel splicing isoforms nor circRNA were detected in HaAPN2 and HaAPN5. Distinct from the conventional GT/AG splicing signal, short co-directional repeats were involved in the splicing of the linear and circular isoforms of HaAPN1 and HaAPN3. Identification of the splice sites revealed that the linear isoforms may be related in some way to the circularization. Moreover, phylogenetic analysis and detection of circRNA of the APN gene of the diamondback moth, Plutella xylostella (PxAPN3), showed that circRNA formation is relatively conserved during the lepidopteran evolutionary process. These results contribute to an improved understanding of lepidopteran APNs and this novel class of insect circRNA.