Comparison of the performance of multiple whole-genome sequence-based tools for the identification of Bacillus cereus sensu stricto biovar Thuringiensis

The Bacillus cereus sensu stricto (s.s.) species comprises strains of biovar Thuringiensis (Bt) known for their bioinsecticidal activity, as well as strains with foodborne pathogenic potential. Bt strains are identified (i) based on the production of insecticidal crystal proteins, also known as Bt toxins, or (ii) based on the presence of cry, cyt, and vip genes, which encode Bt toxins. Multiple bioinformatics tools have been developed for the detection of crystal protein-encoding genes based on whole-genome sequencing (WGS) data. However, the performance of these tools is yet to be evaluated using phenotypic data. Thus, the goal of this study was to assess the performance of four bioinformatics tools for the detection of crystal protein-encoding genes. The accuracy of sequence-based identification of Bt was determined in reference to phenotypic microscope-based screening for the production of crystal proteins. A total of 58 diverse B. cereus sensu lato strains isolated from clinical, food, environmental, and commercial biopesticide products underwent WGS. Isolates were examined for crystal protein production using phase contrast microscopy. Crystal protein-encoding genes were detected using BtToxin_Digger, BTyper3, IDOPS (identification of pesticidal sequences), and Cry_processor. Out of 58 isolates, the phenotypic production of crystal proteins was confirmed for 18 isolates. Specificity and sensitivity of Bt identification based on sequences were 0.85 and 0.94 for BtToxin_Digger, 0.97 and 0.89 for BTyper3, 0.95 and 0.94 for IDOPS, and 0.88 and 1.00 for Cry_processor, respectively. Cry_processor predicted crystal protein production with the highest specificity, and BtToxin_Digger and IDOPS predicted crystal protein production with the highest sensitivity. Three out of four tested bioinformatics tools performed well overall, with IDOPS achieving high sensitivity and specificity (>0.90).IMPORTANCEStrains of Bacillus cereus sensu stricto (s.s.) biovar Thuringiensis (Bt) are used as organic biopesticides. Bt is differentiated from the foodborne pathogen Bacillus cereus s.s. by the production of insecticidal crystal proteins. Thus, reliable genomic identification of biovar Thuringiensis is necessary to ensure food safety and facilitate risk assessment. This study assessed the accuracy of whole-genome sequencing (WGS)-based identification of Bt compared to phenotypic microscopy-based screening for crystal protein production. Multiple bioinformatics tools were compared to assess their performance in predicting crystal protein production. Among them, identification of pesticidal sequences performed best overall at WGS-based Bt identification.

A chromosome-level genome assembly of the soybean pod borer: insights into larval transcriptional response to transgenic soybean expressing the pesticidal Cry1Ac protein

BACKGROUND

Genetically modified (GM) crop plants with transgenic expression of Bacillus thuringiensis (Bt) pesticidal proteins are used to manage feeding damage by pest insects. The durability of this technology is threatened by the selection for resistance in pest populations. The molecular mechanism(s) involved in insect physiological response or evolution of resistance to Bt is not fully understood.

RESULTS

To investigate the response of a susceptible target insect to Bt, the soybean pod borer, Leguminivora glycinivorella (Lepidoptera: Tortricidae), was exposed to soybean, Glycine max, expressing Cry1Ac pesticidal protein or the non-transgenic parental cultivar. Assessment of larval changes in gene expression was facilitated by a third-generation sequenced and scaffolded chromosome-level assembly of the L. glycinivorella genome (657.4 Mb; 27 autosomes + Z chromosome), and subsequent structural annotation of 18,197 RefSeq gene models encoding 23,735 putative mRNA transcripts. Exposure of L. glycinivorella larvae to transgenic Cry1Ac G. max resulted in prediction of significant differential gene expression for 204 gene models (64 up- and 140 down-regulated) and differential splicing among isoforms for 10 genes compared to unexposed cohorts. Differentially expressed genes (DEGs) included putative peritrophic membrane constituents, orthologs of Bt receptor-encoding genes previously linked or associated with Bt resistance, and those involved in stress responses. Putative functional Gene Ontology (GO) annotations assigned to DEGs were significantly enriched for 36 categories at GO level 2, respectively. Most significantly enriched cellular component (CC), biological process (BP), and molecular function (MF) categories corresponded to vacuolar and microbody, transport and metabolic processes, and binding and reductase activities. The DEGs in enriched GO categories were biased for those that were down-regulated (≥ 0.783), with only MF categories GTPase and iron binding activities were bias for up-regulation genes.

CONCLUSIONS

This study provides insights into pathways and processes involved larval response to Bt intoxication, which may inform future unbiased investigations into mechanisms of resistance that show no evidence of alteration in midgut receptors.

Mutational analysis of the transmembrane α4-helix of Bacillus thuringiensis mosquito-larvicidal Cry4Aa toxin

Cry4Aa, produced by Bacillus thuringiensis subsp. israelensis, exhibits specific toxicity to larvae of medically important mosquito genera. Cry4Aa functions as a pore-forming toxin, and a helical hairpin (α4-loop-α5) of domain I is believed to be the transmembrane domain that forms toxin pores. Pore formation is considered to be a central mode of Cry4Aa action, but the relationship between pore formation and toxicity is poorly understood. In the present study, we constructed Cry4Aa mutants in which each polar amino acid residues within the transmembrane α4 helix was replaced with glutamic acid. Bioassays using Culex pipiens mosquito larvae and subsequent ion permeability measurements using symmetric KCl solution revealed an apparent correlation between toxicity and toxin pore conductance for most of the Cry4Aa mutants. In contrast, the Cry4Aa mutant H178E was a clear exception, almost losing its toxicity but still exhibiting a moderately high conductivity of about 60% of the wild-type. Furthermore, the conductance of the pore formed by the N190E mutant (about 50% of the wild-type) was close to that of H178E, but the toxicity was significantly higher than that of H178E. Ion selectivity measurements using asymmetric KCl solution revealed a significant decrease in cation selectivity of toxin pores formed by H178E compared to N190E. Our data suggest that the toxicity of Cry4Aa is primarily pore related. The formation of toxin pores that are highly ion-permeable and also highly cation-selective may enhance the influx of cations and water into the target cell, thereby facilitating the eventual death of mosquito larvae.

Bacillus thuringiensis Cry9Aa Insecticidal Protein Domain I Helices α3 and α4 Are Two Core Regions Involved in Oligomerization and Toxicity

Bacillus thuringiensis Cry9 proteins show high insecticidal activity against different lepidopteran pests. Cry9 could be a valuable alternative to Cry1 proteins because it showed a synergistic effect with no cross-resistance. However, the pore-formation region of the Cry9 proteins is still unclear. In this study, nine mutations of certain Cry9Aa helices α3 and α4 residues resulted in a complete loss of insecticidal activity against the rice pest Chilo suppressalis; however, the protein stability and receptor binding ability of these mutants were not affected. Among these mutants, Cry9Aa-D121R, Cry9Aa-D125R, Cry9Aa-D163R, Cry9Aa-E165R, and Cry9Aa-D167R are unable to form oligomers in vitro, while the oligomers formed by Cry9Aa-R156D, Cry9Aa-R158D, and Cry9Aa-R160D are unstable and failed to insert into the membrane. These data confirmed that helices α3 and α4 of Cry9Aa are involved in oligomerization, membrane insertion, and toxicity. The knowledge of Cry9 pore-forming action may promote its application as an alternative to Cry1 insecticidal proteins.

Coffee Cell Suspensions as a Platform for Transient Gene Expression Analysis

Coffea arabica L. is a crucial crop globally, but its genetic homogeneity leads to its susceptibility to diseases and pests like the coffee berry borer (CBB). Chemical and cultural control methods are difficult due to the majority of the CBB life cycle taking place inside coffee beans. One potential solution is the use of the gene cyt1Aa from Bacillus thuringiensis as a biological insecticide. To validate candidate genes against CBB, a simple, rapid, and efficient transient expression system is necessary. This study uses cell suspensions as a platform for expressing the cyt1Aa gene in the coffee genome (C. arabica L. var. Catuaí) to control CBB. The Agrobacterium tumefaciens strain GV3101::pMP90 containing the bar and cyt1Aa genes are used to genetically transform embryogenic cell suspensions. PCR amplification of the cyt1Aa gene is observed 2, 5, and 7 weeks after infection. This chapter describes a protocol that can be used for the development of resistant varieties against biotic and abiotic stresses and CRISPR/Cas9-mediated genome editing.

A soybean trypsin inhibitor reduces the resistance to transgenic maize in a population of Spodoptera frugiperda (Lepidoptera: Noctuidae)

Lepidopteran pests have been successfully managed by the adoption of insect resistant transgenic plants expressing Cry and/or Vip insecticidal proteins derived from Bacillus thuringiensis (Bt plants). Among such pests, Spodoptera frugiperda (Smith, 1797) (Lepidoptera: Noctuidae) is highlighted for its destructive potential in maize crops and for cases of field-evolved resistance to Bt plants. Cry insecticidal proteins expressed in Bt plants are known for their interaction with insect midgut receptors and subsequent midgut cell disruption that leads to target pest death. In the midgut of lepidopteran larval pests such as S. frugiperda, serine proteases are important in dietary protein digestion and activation or degradation of insecticidal proteins. This work was conducted to evaluate if the use of a soybean trypsin inhibitor (SBTI) could disrupt the development of a Bt-susceptible and a Bt-resistant population of S. frugiperda ingesting Bt (expressing Cry1F, Cry1A.105, and Cry2Ab2 Cry proteins) and non-Bt maize plants. The SBTI was produced and purified using recombinant expression in E. coli followed by purification in Ni-Sepharose. Bioassays using non-Bt maize leaves indicated that the development of susceptible and resistant populations of S. frugiperda was not influenced by the ingestion of SBTI. However, when the resistant population consumed Bt maize plants amended with SBTI, high mortality along with a reduction in larval weight and reduced activity of digestive trypsins were observed. Although the mode of action was not elucidated, it is possible that the consumption of SBTI increased susceptibility to Bt maize in the resistant population of S. frugiperda.

Channel Formation in Cry Toxins: An Alphafold-2 Perspective

Bacillus thuringiensis (Bt) strains produce pore-forming toxins (PFTs) that attack insect pests. Information for pre-pore and pore structures of some of these Bt toxins is available. However, for the three-domain (I-III) crystal (Cry) toxins, the most used Bt toxins in pest control, this crucial information is still missing. In these Cry toxins, biochemical data have shown that 7-helix domain I is involved in insertion in membranes, oligomerization and formation of a channel lined mainly by helix α4, whereas helices α1 to α3 seem to have a dynamic role during insertion. In the case of Cry1Aa, toxic against Manduca sexta larvae, a tetrameric oligomer seems to precede membrane insertion. Given the experimental difficulty in the elucidation of the membrane insertion steps, we used Alphafold-2 (AF2) to shed light on possible oligomeric structural intermediates in the membrane insertion of this toxin. AF2 very accurately (<1 Å RMSD) predicted the crystal monomeric and trimeric structures of Cry1Aa and Cry4Ba. The prediction of a tetramer of Cry1Aa, but not Cry4Ba, produced an 'extended model' where domain I helices α3 and α2b form a continuous helix and where hydrophobic helices α1 and α2 cluster at the tip of the bundle. We hypothesize that this represents an intermediate that binds the membrane and precedes α4/α5 hairpin insertion, together with helices α6 and α7. Another Cry1Aa tetrameric model was predicted after deleting helices α1 to α3, where domain I produced a central cavity consistent with an ion channel, lined by polar and charged residues in helix α4. We propose that this second model corresponds to the 'membrane-inserted' structure. AF2 also predicted larger α4/α5 hairpin n-mers (14 ≤n ≤ 17) with high confidence, which formed even larger (~5 nm) pores. The plausibility of these models is discussed in the context of available experimental data and current paradigms.

Distribution and Genetic Diversity of Genes from Brazilian Bacillus thuringiensis Strains Toxic to Agricultural Insect Pests Revealed by Real-Time PCR

Bacillus thuringiensis is a Gram-positive aerobic bacterium and the most used biopesticide worldwide. Given the importance of B. thuringiensis strain characterization for the development of new bioinsecticides or transgenic events and the identification and classification of new B. thuringiensis genes and strains to understand its distribution and diversity, this work is aimed at creating a gene identification system based on qPCR reactions utilizing core B. thuringiensis genes cry1, cry2, cry3, cry4, cry5, app6, cry7, cry8, cry9, cry10, cry11, vpb1, vpa2, vip3, cyt1, and cyt2 for the characterization of 257 strains of B. thuringiensis. This system was based on the Invertebrate Bacteria Collection from Embrapa Genetic Resources and Biotechnology and analyzed (a) the degree of correlation between the distribution of these strains and the origin of the substrate from which the strain was isolated and (b) between its distribution and geoclimatic conditions. This study made it possible to observe that the cry1, cry2, and vip3A/B genes occur homogeneously in the Brazilian territory, and some genes are found in specific regions. The biggest reservoir of variability is within B. thuringiensis strains in each region, and it is suggested that both geoclimatic conditions and regional crops interfere with the genetic diversity of the B. thuringiensis strains present in the region, and B. thuringiensis strains can constantly exchange genetic information.

Inheritance and fitness cost of laboratory-selected resistance to Vip3Aa in Helicoverpa zea (Lepidoptera: Noctuidae)

The polyphagous pest Helicoverpa zea (Lepidoptera: Noctuidae) has evolved practical resistance to transgenic corn and cotton producing Cry1 and Cry2 crystal proteins from Bacillus thuringiensis (Bt) in several regions of the United States. However, the Bt vegetative insecticidal protein Vip3Aa produced by Bt corn and cotton remains effective against this pest. To advance knowledge of resistance to Vip3Aa, we selected a strain of H. zea for resistance to Vip3Aa in the laboratory. After 28 generations of continuous selection, the resistance ratio was 267 for the selected strain (GA-R3) relative to a strain not selected with Vip3Aa (GA). Resistance was autosomal and almost completely recessive at a concentration killing all individuals from GA. Declines in resistance in heterogeneous strains containing a mixture of susceptible and resistant individuals reared in the absence of Vip3Aa indicate a fitness cost was associated with resistance. Previously reported cases of laboratory-selected resistance to Vip3Aa in lepidopteran pests often show partially or completely recessive resistance at high concentrations and fitness costs. Abundant refuges of non-Bt host plants can maximize the benefits of such costs for sustaining the efficacy of Vip3Aa against target pests.

Cry11Aa and Cyt1Aa exhibit different structural orders in crystal topography

Cry11Aa and Cyt1Aa are two pesticidal toxins produced by Bacillus thuringiensis subsp. israelensis. To improve our understanding of the nature of their oligomers in the toxic actions and synergistic effects, we performed the atomic force microscopy to probe the surfaces of their natively grown crystals, and used the L-weight filter to enhance the structural features. By L-weight filtering, molecular sizes of the Cry11Aa and Cyt1Aa monomers obtained are in excellent agreement with the three-dimensional structures determined by x-ray crystallography. Moreover, our results show that the layered feature of a structural element distinguishes the topographic characteristics of Cry11Aa and Cyt1Aa crystals, suggesting that the Cry11Aa toxin has a better chance than Cyt1Aa for multimerization and therefore cooperativeness of the toxic actions.

Group Selection as a Basis for Screening Mutagenized Libraries of Public Goods (Bacillus thuringiensis Cry Toxins)

The pesticidal toxins of Bacillus thuringiensis (Bt) supply the active proteins for genetically modified insect-resistant crops. There is therefore keen interest in finding new toxins, or improving known toxins, in order to increase the mortality of various targets. The production and screening of large libraries of mutagenized toxins are among the means of identifying improved toxins. Since Cry toxins are public goods, and do not confer advantages to producers in competition, conventional directed evolution approaches cannot be used here. Instead, thousands of individual mutants have to be sequenced and assayed individually, a costly and time-consuming process. In this study, we tested a group selection-based approach that could be used to screen an uncharacterized pool of Cry toxin mutants. This involved selecting for infectivity between subpopulations of Bt clones within metapopulations of infected insects in three rounds of passage. We also tested whether additional mutagenesis from exposure to ethyl methanesulfonate could increase infectivity or supply additional Cry toxin diversity during passage. Sequencing of pools of mutants at the end of selection showed that we could effectively screen out Cry toxin variants that had reduced toxicity with our group selection approach. The addition of extra mutagenesis during passage decreased the efficiency of selection for infectivity and did not produce any additional novel toxin diversity. Toxins with loss-of-function mutations tend to dominate mutagenized libraries, and so a process for screening out these mutants without time-consuming sequencing and characterization steps could be beneficial when applied to larger libraries. IMPORTANCE Insecticidal toxins from the bacterium Bacillus thuringiensis are widely exploited in genetically modified plants. This application creates a demand for novel insecticidal toxins that can be used to better manage resistant pests or control new or recalcitrant target species. An important means of producing novel toxins is via high-throughput mutagenesis and screening of existing toxins, a lengthy and resource-intensive process. This study describes the development and testing of an efficient means of screening a test library of mutagenized insecticidal toxins. Here, we showed that it is possible to screen out loss-of-function mutations with low infectivity within a pool without the need to characterize and sequence each mutant individually. This has the potential to improve the efficiency of processes used to identify novel proteins.

More than 10 years after commercialization, Vip3A-expressing MIR162 remains highly efficacious in controlling major Lepidopteran maize pests: laboratory resistance selection versus field reality

MIR162, a maize event that expresses Vip3Aa20 (Vip3A) approved for commercial cultivation around 2010, has been excellent for control of major Lepidopteran pests. However, development of fall armyworm (FAW) resistance to Vip3A is a serious concern. Resistant colonies selected in the laboratory can serve as valuable tools not only for better understanding of Vip3A's mode of action (MOA) and mechanism of resistance (MOR) but also for screening novel leads of new MOA that will help control FAW in case resistance to Vip3A in the field becomes a reality. We selected a Vip3A-resistant FAW strain, FAW_Vip3AR, by subjecting a FAW founder population containing field genetics to Vip3A exposure. FAW_Vip3AR had >9800-fold resistance to Vip3A by diet surface overlay bioassays and resistance was stable. Feeding bioassays using detached leaf tissues or whole plants indicated that FAW_Vip3AR larvae readily fed and completed the full life cycle on Vip3A-expressing MIR162 maize plants and leaf tissues that killed 100% of susceptible larvae. Yet, FAW_Vip3AR faced at least two challenges. First, FAW_Vip3AR suffered an apparent disadvantage (incomplete resistance) when feeding on MIR162 in comparison to FAW_Vip3AR feeding on Vip3A-free isoline AX5707 maize; and second, FAW_Vip3AR showed a fitness costs in comparison to a Vip3A-susceptible strain when both fed on AX5707. We also demonstrated that, >10 years after commercialization, MIR162 and Vip3A remain highly efficacious against field populations of three major Lepidopteran pests from different geographic locations and FAW strains resistant to other Bacillus thuringiensis (Bt) toxins that are currently on the market.

A Modified F2 Screen for Estimating Cry1Ac and Cry2Ab Resistance Allele Frequencies in Helicoverpa zea (Lepidoptera: Noctuidae)

Evaluating the frequency of resistance alleles is important for resistance management and sustainable use of transgenic crops that produce insecticidal proteins from Bacillus thuringiensis. Helicoverpa zea (Boddie) (Lepidoptera: Noctuidae) is a major crop pest in the United States that has evolved practical resistance to the crystalline (Cry) proteins in Bt corn and cotton. The standard F2 screen for estimating resistance allele frequency does not work well for H. zea because successful single-pair matings are rare. In this study, we developed and implemented a modified F2 screen for H. zea that generates F1 progeny by crossing three laboratory susceptible female moths with one feral male moth instead of single-pair crosses. During 2019-2020, we used this modified method to establish 192 F2 families from 623 matings between susceptible females and feral males from Arkansas, Louisiana, Mississippi, and Tennessee. From each F2 family, we screened 128 neonates against discriminating concentrations of Cry1Ac and Cry2Ab in diet overlay bioassays. Based on these discriminating concentration bioassays, families were considered positive for resistance if at least five larvae survived to second instar, including at least one to third instar. The percentage of positive families was 92.7% for Cry1Ac and 38.5% for Cry2Ab, which yields an estimated resistance allele frequency (with 95% confidence interval) of 0.722 (0.688-0.764) for Cry1Ac and 0.217 (0.179-0.261) for Cry2Ab. The modified F2 screen developed and implemented here may be useful for future resistance monitoring studies of H. zea and other pests.

Monitoring Insect Resistance to Bt Maize in the European Union: Update, Challenges, and Future Prospects

Transgenic maize producing the Cry1Ab toxin of Bacillus thuringiensis (Bt maize) was approved for cultivation in the European Union (EU) in 1998 to control the corn borers Sesamia nonagrioides (Lefèbvre) and Ostrinia nubilalis (Hübner). In the EU since then, Cry1Ab is the only Bt toxin produced by Bt maize and Spain is the only country where Bt maize has been planted every year. In 2021, about 100,000 hectares of Bt maize producing Cry1Ab were cultivated in the EU, with Spain accounting for 96% and Portugal 4% of this area. In both countries, Bt maize represented less than 25% of all maize planted in 2021, with a maximum regional adoption of 64% Bt maize in northeastern Spain. Insect resistance management based on the high-dose/refuge strategy has been implemented in the EU since 1998. This has been accompanied by monitoring to enable early detection of resistance. The monitoring data from laboratory bioassays show no decrease in susceptibility to Cry1Ab had occurred in either pest as of 2021. Also, control failures have not been reported, confirming that Bt maize producing Cry1Ab remains effective against both pests. Conditions in the EU preventing approval of new genetically modified crops, including maize producing two or more Bt toxins targeting corn borers, may limit the future effectiveness of resistance management strategies.

Evaluating response mechanisms of soil microbiomes and metabolomes to Bt toxin additions

The accumulation and persistence of Bt toxins in soils from Bt plants and Bt biopesticides may result in environmental hazards such as adverse impacts on soil microorganisms. However, the dynamic relationships among exogenous Bt toxins, soil characteristics, and soil microorganisms are not well understood. Cry1Ab is one of the most commonly used Bt toxins and was added to soils in this study to evaluate subsequent changes in soil physiochemical properties, microbial taxa, microbial functional genes, and metabolites profiles via 16S rRNA gene pyrosequencing, high-throughput qPCR, metagenomic shotgun sequencing, and untargeted metabolomics. Higher additions of Bt toxins led to higher concentrations of soil organic matter (SOM), ammonium (NH⁺₄-N), and nitrite (NO₂^--N) compared against controls without addition after 100 days of soil incubation. High-throughput qPCR analysis and shotgun metagenomic sequencing analysis revealed that the 500 ng/g Bt toxin addition significantly affected profiles of soil microbial functional genes involved in soil carbon (C), nitrogen (N), and phosphorus (P) cycling after 100 days of incubation. Furthermore, combined metagenomic and metabolomic analyses indicated that the 500 ng/g Bt toxin addition significantly altered low molecular weight metabolite profiles of soils. Importantly, some of these altered metabolites are involved in soil nutrient cycling, and robust associations were identified among differentially abundant metabolites and microorganisms due to Bt toxin addition treatments. Taken together, these results suggest that higher levels of Bt toxin addition can alter soil nutrients, probably by affecting the activities of Bt toxin-degrading microorganisms. These dynamics would then activate other microorganisms involved in nutrient cycling, finally leading to broad changes in metabolite profiles. Notably, the addition of Bt toxins did not cause the accumulation of potential microbial pathogens in soils, nor did it adversely affect the diversity and stability of microbial communities. This study provides new insights into the putative mechanistic associations among Bt toxins, soil characteristics, and microorganisms, providing new understanding into the ecological impacts of Bt toxins on soil ecosystems.

pH-Dependent Changes in Structural Stabilities of Bt Cry1Ac Toxin and Contrasting Model Proteins following Adsorption on Montmorillonite

The environmental fate of insecticidal Cry proteins, including time-dependent conservation of biological properties, results from their structural stability in soils. The complex cascade of reactions involved in biological action requires Cry proteins to be in solution. However, the pH-dependent changes in conformational stability and the adsorption-desorption mechanisms of Cry protein on soil minerals remain unclear. We used Derjaguin-Landau-Verwey-Overbeek (DLVO) calculation and differential scanning calorimetry to interpret the driving forces and structural stabilities of Cry1Ac and two contrasting model proteins adsorbed by montmorillonite. The structural stability of Cry1Ac is closer to that of the "hard" protein, α-chymotrypsin, than that of the "soft" bovine serum albumin (BSA). The pH-dependent adsorption of Cry1Ac and α-chymotrypsin could be explained by DLVO theory, whereas the BSA adsorption deviated from it. Patch-controlled electrostatic attraction, hydrophobic effects, and entropy changes following protein unfolding on a mineral surface could contribute to Cry1Ac adsorption. Cry1Ac, like chymotrypsin, was partly denatured on montmorillonite, and its structural stability decreased with an increase in pH. Moreover, small changes in the conformational heterogeneity of both Cry1Ac and chymotrypsin were observed following adsorption. Conversely, adsorbed BSA was completely denatured regardless of the solution pH. The moderate conformational rearrangement of adsorbed Cry1Ac may partially explain why the insecticidal activity of Bt toxin appears to be conserved in soils, albeit for a relatively short time period.

Novel miR-108 and miR-234 target juvenile hormone esterase to regulate the response of Plutella xylostella to Cry1Ac protoxin

Insect hormones, such as juvenile hormone (JH), precisely regulate insect life-history traits. The regulation of JH is tightly associated with the tolerance or resistance to Bacillus thuringiensis (Bt). JH esterase (JHE) is a primary JH-specific metabolic enzyme which plays a key role in regulating JH titer. Here, we characterized a JHE gene from Plutella xylostella (PxJHE), and found it was differentially expressed in the Bt Cry1Ac resistant and susceptible strains. Suppression of PxJHE expression with RNAi increased the tolerance of P. xylostella to Cry1Ac protoxin. To investigate the regulatory mechanism of PxJHE, two target site prediction algorithms were applied to predict the putative miRNAs targeting PxJHE, and the resulting putative miRNAs were subsequently verified for their function targeting PxJHE using luciferase reporter assay and RNA immunoprecipitation. MiR-108 or miR-234 agomir delivery dramatically reduced PxJHE expression in vivo, whilst only miR-108 overexpression consequently increased the tolerance of P. xylostella larvae to Cry1Ac protoxin. By contrast, reduction of miR-108 or miR-234 dramatically increased PxJHE expression, accompanied by the decreased tolerance to Cry1Ac protoxin. Furthermore, injection of miR-108 or miR-234 led to developmental defects in P. xylostella, whilst injection of antagomir did not cause any obvious abnormal phenotypes. Our results indicated that miR-108 or miR-234 can be applied as potential molecular targets to combat P. xylostella and perhaps other lepidopteran pests, providing novel insights into miRNA-based integrated pest management.

Relative fitness of susceptible and Cry1A.105/Cry2Ab2-single-/dual-protein-resistant Helicoverpa zea (Boddie) (Lepidoptera: Noctuidae) on non-Bt diet and a diet containing a low concentration of two proteins

Helicoverpa zea (Boddie) is a destructive agricultural pest species that is targeted by both Bacillus thuringiensis (Bt) maize and cotton in the United States. Cry1A.105 and Cry2Ab2 are two Bt proteins expressed in a widely planted maize event MON 89034. In this study, two tests (Test-I and Test-II) were conducted to evaluate the relative fitness of Bt-susceptible and -resistant H. zea on non-Bt diet (Test-I and Test-II) and a diet containing a mix of Cry1A.105 and Cry2Ab2 at a low concentration (Test-II only). Insect populations evaluated in Test-I were two Bt-susceptible strains and three Bt-resistant strains (a single-protein Cry1A.105-, a single-protein Cry2Ab2-, and a dual-protein Cry1A.105/Cry2Ab2-resistant strains). Test-II analyzed the same two susceptible strains, three backcrossed-and-reselected Cry1A.105/Cry2Ab2-single-/dual-protein-resistant strains, and three F₁ heterozygous strains. Measurements of life table parameters showed that neither the single- nor dual-protein Cry1A.105/Cry2Ab2 resistance in H. zea was associated with fitness costs under the test conditions. The single Cry protein resistances at a concentration of a mix of Cry1A.105 and Cry2Ab2 that resulted in a zero net reproductive rate for the two susceptible strains were functionally incomplete recessive or codominant, and the dual-protein resistance was completely dominant. The lack of fitness costs could be a factor contributing to the rapid revolution of resistance to the Cry proteins in this species. Data generated from this study should aid our understanding of Cry protein resistance evolution and help in refining IRM programs for H. zea.

Seven Years of Monitoring Susceptibility to Cry1Ab and Cry1F in Asian Corn Borer

Resistance monitoring in the Asian corn borer, Ostrinia furnacalis, is necessary to accommodate the commercial introduction and stewardship of Bt maize in China. The susceptibility of 56 O. furnacalis field populations, collected between 2015 and 2021 from the corn belt regions of China, to Cry1Ab and Cry1F toxins was determined. Neonate larvae (within 12 h after hatching) were placed on the surface of semi-artificial agar-free diet incorporating a series of concentrations of purified toxins, and mortality was evaluated after 7d. The median lethal concentration (LC₅₀) values of Cry1Ab and Cry1F were 0.05 to 0.37 µg/g (protein/diet) and 0.10 to 1.22 µg/g, respectively. Although interpopulation variation in susceptibility to the toxins was observed, the magnitude of the differences was 5.8-fold and 8.3-fold for Cry1Ab and Cry1F, respectively. These results suggested that the observed susceptibility differences reflect natural geographical variation in response and not variation caused by prior exposure to selection pressures. Therefore, the O. furnacalis populations were apparently still susceptible to Cry1Ab and Cry1F across their range within China. The monitoring data established here will serve as a comparative reference for early warning signs of field-evolved resistance after the cultivation of Bt maize in China.

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.

Bt maize can provide non-chemical pest control and enhance food safety in China

China is the world's second-largest maize producer and consumer. In recent years, the invasive fall armyworm Spodoptera frugiperda (J.E. Smith) has adversely affected maize productivity and compromised food security. To mitigate pest-inflicted food shortages, China's Government issued biosafety certificates for two genetically modified (GM) Bt maize hybrids, Bt-Cry1Ab DBN9936 and Bt-Cry1Ab/Cry2Aj Ruifeng 125, in 2019. Here, we quantitatively assess the impact of both Bt maize hybrids on pest feeding damage, crop yield and food safety throughout China's maize belt. Without a need to resort to synthetic insecticides, Bt maize could mitigate lepidopteran pest pressure by 61.9-97.3%, avoid yield loss by 16.4-21.3% (range -11.9-99.2%) and lower mycotoxin contamination by 85.5-95.5% as compared to the prevailing non-Bt hybrids. Yield loss avoidance varied considerably between experimental sites and years, as mediated by on-site infestation pressure and pest identity. For either seed mixtures or block refuge arrangements, pest pressure was kept below established thresholds at 90% Bt maize coverage in Yunnan (where S. frugiperda was the dominant species) and 70% Bt maize coverage in other sites dominated by Helicoverpa armigera (Hübner) and Ostrinia furnacalis (Guenée). Drawing on experiences from other crop/pest systems, Bt maize in se can provide area-wide pest management and thus, contribute to a progressive phase-down of chemical pesticide use. Hence, when consciously paired with agroecological and biodiversity-based measures, GM insecticidal crops can ensure food and nutrition security, contribute to the sustainable intensification of China's agriculture and reduce food systems' environmental footprint.

Segregation of Cry Genes in the Seeds Produced by F1 Bollgard® II Cotton Differs between Hybrids: Could This Be Linked to the Observed Field Resistance in the Pink Bollworm?

Indian populations of the Pink Bollworm (PBW) are resistant to Bt (Bacillus thuringiensis) cotton hybrids containing Cry1Ac and Cry2Ab genes. Segregation of these Cry genes in F₁ hybrids could subject PBW to sublethal concentrations. Moreover, planting hybrids with varying zygosities of Cry genes could produce diverse segregation patterns and expose PBW populations to highly variable toxin concentrations. This could potentially promote the rate of resistance development. Therefore, we studied the segregation patterns of Cry genes in different commercial Bt hybrids cultivated in India. Results showed that two hybrids segregated according to the Mendelian mono-hybrid ratio, three segregated according to the Mendelian di-hybrid ratio, and one showed a mixed segregation pattern. The assortment of seeds containing Cry genes varied between bolls of the same hybrid. In India, different Bt cotton hybrids are cultivated in small patches next to each other, exposing PBW populations to sublethal doses and wide variations in the occurrence of Cry genes. It is necessary to avoid segregation of Cry genes in the seeds produced by F₁ hybrids. This study recommends using Bt parents homozygous for Cry genes in commercial Bt cotton hybrid development. This breeding strategy could be effective for similar transgenic crop hybrids as well.

Potential Management Options for the Invasive Moth Spodoptera frugiperda in Europe

We here review and discuss management options that growers in Europe could take in response to the expected invasion of the fall armyworm, Spodoptera frugiperda (Lepidoptera: Noctuidae). The focus is put on maize but the information provided is also relevant for other crops potentially affected. A sound forecasting system for fall armyworm both on a regional as well as at local scale should be established to alert growers as early as possible. Whilst a number of cultural control methods are adopted by maize growers in different regions globally to fight fall armyworm, many of them may either not be highly effective, too laborious, or otherwise unfeasible within the mechanized crop production systems used in Europe. Potential is seen in the stimulation of natural enemies through conservation biocontrol approaches, e.g., the planting of flower strips or intermediate cover crops, reducing tillage intensity, and avoiding broad-spectrum insecticides. To manage fall armyworm infestations, several effective biologically-based products are available globally, and some in Europe, e.g., based on specific baculoviruses, certain Bacillus thuringiensis strains, few entomopathogenic nematodes, and a number of botanicals. These should be given priority to avoid a major influx of insecticides into the maize agro-ecosystem once the fall armyworm arrives and in case growers are not prepared. Plant protection companies, particularly biocontrol companies should act proactively in starting registration of ingredients and products against fall armyworm in Europe. European maize growers should be made aware, in time, of key features of this new invasive pest and appropriate control options.

Safety of Bacillus thuringiensis Cry1Ah and Vip3Aa toxins for the predatory stink bug Arma custos (Hemiptera: Pentatomidae)

The widespread adoption of Bt crops expressing insecticidal proteins derived from Bacillus thuringiensis has created a need to assess the potential effects of these toxins on non-target organisms, especially species such as Arma custos, a generalist predator that provides important biological control services in many field crops in Asia. Direct dietary exposure of A. custos to Cry1Ah and Vip3Aa proteins produced no adverse effects on life history traits, despite continuous exposure throughout development and early adult life to concentrations significantly higher than the Bt protein concentration likely encountered by A.custos in the field, even when feeding directly on Bt plants. Enzyme-linked immunosorbent assay confirmed the presence of Bt proteins in A. custos midguts, but quantitative real-time PCR analysis of 12 genes associated with detoxification, antioxidative responses, immune responses, and metabolism revealed no significant changes in expression in adult bugs. Indirect exposure to these toxins via consumption of intoxicated prey, larvae of Helicoverpa armigera (Hübner), likewise produced no negative impacts on survival, development, adult weight, or female fecundity in either the F0 (exposed) or F1 (unexposed) generation, but female fresh weight was reduced in the F0 generation by the Cry1Ah (50 μg/g) treatment. Finally, a competitive binding assay with labelled protein and a ligand blotting assay both demonstrated that the Cry1Ah protein could not bind to receptors on the midgut brush border membrane vesicles (BBMVs) of A. custos adults. Therefore, we conclude that Cry1Ah and Vip3Aa proteins are unlikely to have significant negative effects on A. custos populations if employed as plant-incorporated protectants in field crops.

Varied frequencies of resistance alleles to Cry1Ab and Cry1Ac among Brazilian populations of the sugarcane borer, Diatraea saccharalis (F.)

BACKGROUND

Brazil is the largest grower of the world's 26 million ha of sugarcane, Saccharum officinarum. Pest damage mainly by the sugarcane borer, Diatraea saccharalis (F.), is a great challenge to the sugarcane industry. To control D. saccharalis, Brazil launched the world's first commercial use of Bt sugarcane in 2017. As part of the resistance management programs for Bt sugarcane planting, 535 F₂ isoline families of D. saccharalis collected from three major sugarcane planting states (Goiás, Minas Gerais and São Paulo) in Brazil during 2019-2020 were screened for resistance to two Bt sugarcane varieties: CTC20BT expressing Cry1Ab and CTC9001BT expressing Cry1Ac. Here we report the results of the first study related to Bt resistance in a sugarcane cropping system.

RESULTS

Larval survivorships of these families in an F₂ screen on CTC20BT were highly correlated with their survival on CTC9001BT, whereas the Cry1Ac tissues exhibited greater insecticidal activities than Cry1Ab. Resistance allele frequencies (RAFs) for populations from Goiás and Minas Gerais were relatively low at 0.0034 for Cry1Ab and 0.0045 to Cry1Ac. By contrast, RAFs for São Paulo populations were considerably greater (0.0393 to Cry1Ab, 0.0245 to Cry1Ac).

CONCLUSIONS

RAFs to Cry1Ab and Cry1Ac varied among Brazilian D. saccharalis populations. Prior selection resulting from an intensive use of single-gene Bt maize under low compliance of refuge planting could be a main factor contributing to the high RAF in São Paulo. The results suggest that mitigation measures including sufficient non-Bt maize refuge planting, effective resistance monitoring, and use of pyramided Bt sugarcane traits should be implemented promptly to prevent further increase in the RAF to ensure the sustainable use of Bt sugarcane in Brazil.

MINI ABSTRACT

To control Diatraea saccharalis, Brazil launched the world's first commercial use of Bt sugarcane in 2017. As part of the resistance management programs for Bt sugarcane planting in Brazil, 535 F₂ isoline families of D. saccharalis collected from three major sugarcane planting states (Goiás, Minas Gerais and São Paulo) in Brazil during 2019-2020 were screened for resistance to Cry1Ab and Cry1Ac sugarcane plants Resistance allele frequencies (RAFs) for the populations from Goiás and Minas Gerais were relatively low at 0.0034 for Cry1Ab and 0.0045 to Cry1Ac. By contrast, RAFs for the São Paulo populations were considerably greater (0.0393 to Cry1Ab, 0.0245 to Cry1Ac). Prior selection resulting from an intensive use of single-gene Bt maize under low compliance of non-Bt maize refuge planting could be a main factor contributing to the high RAF in São Paulo. The results suggest that effective mitigation measures including sufficient non-Bt maize refuge planting, effective resistance monitoring and use of pyramided Bt sugarcane traits should be implemented promptly to prevent further increase in the RAF to ensure the sustainable use of Bt sugarcane in Brazil. © 2022 Society of Chemical Industry.

Toxicity of Bacillus thuringiensis Strains to Six Lepidopteran Pests of Brazilian Agricultural Landscape

The first step for biological insecticide production is the evaluating of the pathogenicity and strains virulence of entomopathogens against the target pests. Bacillus thuringiensis (Bt) is widespread in nature, lethal pathogen of a range of orders and the most widely used entomopathogenic biological control agent. Therefore, our objective was to study the susceptibility of six important lepidopteran pests in Brazilian agricultural production systems, Anticarsia gemmatalis (Hübner), Chrysodeixis includens (Walker) Helicoverpa armigera (Hübner), Spodoptera cosmioides (Walker), S. eridania (Cramer), and S. frugiperda (Smith), to Bt strains. In vitro bioassays were conducted with neonate larvae exposed to two Bt strains, 1680A and 775E, in seven concentrations of spores and protein crystals. After 7 days, the mortality was evaluated and indicated that lethal concentration (LC₅₀) for 775E strain ranged from 7.72 × 10⁴ to 2.45 × 10⁶ spores mL^-1 and for 1608A strain from 5.63 × 10³ to 1.21 × 10⁶ spores mL^-1. The strain 1608A was more toxic for A. gemmatalis and H. armigera than 775E strain, and then was further evaluated in greenhouse bioassays with maize and soybean plants infested separately with target insects. The strain 1608A showed an efficient control for most insect species studied in the greenhouse bioassays, promoting considerable reduction in leaf injury and demonstrating a high potential in biological control of important lepidopteran species in crop production systems in Brazil.

Practical resistance to Cry toxins and efficacy of Vip3Aa in Bt cotton against Helicoverpa zea

BACKGROUND

Crops genetically engineered to make insect-killing proteins from Bacillus thuringiensis (Bt) have revolutionized management of some pests. However, the benefits of such transgenic crops are reduced when pests evolve resistance to Bt toxins. We evaluated resistance to Bt toxins and Bt cotton plants using laboratory bioassays and complementary field trials focusing on Helicoverpa zea, one of the most economically important pests of cotton and other crops in the United States.

RESULTS

The data from 235 laboratory bioassays demonstrate resistance to Cry1Ac, Cry1Fa, and Cry2Ab occurred in most of the 95 strains of H. zea derived from Arkansas, Louisiana, Mississippi, Tennessee, and Texas during 2016 to 2021. Complementary field data show efficacy decreased for Bt cotton producing Cry1Ac + Cry1Fa or Cry1Ac + Cry2Ab, but not Cry1Ac + Cry1Fa + Vip3Aa. Moreover, analysis of data paired by field site and year shows higher survival in bioassays was generally associated with lower efficacy of Bt cotton.

CONCLUSIONS

The results confirm and extend previous evidence showing widespread practical resistance of H. zea in the United States to the Cry toxins produced by Bt cotton and corn, but not to Vip3Aa. Despite deployment in combination with Cry toxins in Bt crops, Vip3Aa effectively acts as a single toxin against H. zea larvae that are highly resistant to Cry toxins. Furthermore, Vip3Aa adoption is increasing and previous work provided an early warning of field-evolved resistance. Thus, rigorous resistance management measures are needed to preserve the efficacy of Vip3Aa against this highly adaptable pest. © 2022 Society of Chemical Industry.

Multifunctional Properties of a Bacillus thuringiensis Strain (BST-122): Beyond the Parasporal Crystal

Chemical products still represent the most common form of controlling crop pests and diseases. However, their extensive use has led to the selection of resistances. This makes the finding of new solutions paramount to countering the economic losses that pests and diseases represent in modern agriculture. Bacillus thuringiensis (Bt) is one of the most reliable alternatives to chemical-based solutions. In this study, we aimed to further expand the global applicability of Bt strains beyond their spores and crystals. To this end, we selected a new Bt strain (BST-122) with relevant toxicity factors and tested its activity against species belonging to different phyla. The spore and crystal mixture showed toxicity to coleopterans. Additionally, a novel Cry5-like protein proved active against the two-spotted spider mite. In vivo and plant assays revealed significant control of the parasitic nematode, Meloidogyne incognita. Surprisingly, our data indicated that the nematocidal determinants may be secreted. When evaluated against phytopathogenic fungi, the strain seemed to decelerate their growth. Overall, our research has highlighted the potential of Bt strains, expanding their use beyond the confinements of spores and crystals. However, further studies are required to pinpoint the factors responsible for the wide host range properties of the BST-122 strain.

Knockout of ABC transporter gene ABCA2 confers resistance to Bt toxin Cry2Ab in Helicoverpa zea

Evolution of pest resistance reduces the benefits of widely cultivated genetically engineered crops that produce insecticidal proteins derived from Bacillus thuringiensis (Bt). Better understanding of the genetic basis of pest resistance to Bt crops is needed to monitor, manage, and counter resistance. Previous work shows that in several lepidopterans, resistance to Bt toxin Cry2Ab is associated with mutations in the gene encoding the ATP-binding cassette protein ABCA2. The results here show that mutations introduced by CRISPR/Cas9 gene editing in the Helicoverpa zea (corn earworm or bollworm) gene encoding ABCA2 (HzABCA2) can cause resistance to Cry2Ab. Disruptive mutations in HzABCA2 facilitated the creation of two Cry2Ab-resistant strains. A multiple concentration bioassay with one of these strains revealed it had > 200-fold resistance to Cry2Ab relative to its parental susceptible strain. All Cry2Ab-resistant individuals tested had disruptive mutations in HzABCA2. We identified five disruptive mutations in HzABCA2 gDNA. The most common mutation was a 4-bp deletion in the expected Cas9 guide RNA target site. The results here indicate that HzABCA2 is a leading candidate for monitoring Cry2Ab resistance in field populations of H. zea.

Drop-off behaviour of Bt-resistant and Bt-susceptible Helicoverpa armigera (Hübner) (Lepidoptera: Noctuidae) larvae on Bt-cotton and non-Bt cotton plants

The highest natural mortality rate of larval Lepidoptera in field populations occurs in the first instar, but it is highly variable. The pattern and degree of survival is not easily predicted but depends on their ability to establish on host plants. Lepidopteran larval dispersal behaviour, known as 'drop-off', happens when the host is unsuitable for larvae to settle and begin feeding. Understanding drop-off behaviour of Helicoverpa armigera (Hübner) with and without physiological resistance to Bt toxins on Bt and non-Bt cotton plants is an important component for resistance management strategies for this insect. We examined the drop-off behaviour of H. armigera to determine: (1) whether they move the same way or differently in response to Bt and non-Bt, and (2) could H. armigera larvae detect Bt toxin levels in cotton plants or did they move independently of toxin levels? In this study, we assessed the drop-off behaviour of Bt-resistant and Bt-susceptible H. armigera neonates on artificial diets and cotton plants with and without Bt toxin during the first 12 h after hatching. Bt-resistant and Bt-susceptible H. armigera neonates behaved differently on Bt and non-Bt substrates. The percentages of Bt-resistant larvae that dropped off Bt and non-Bt cotton plants were not significantly different. In contrast, significantly more Bt-susceptible larvae dropped off Bt cotton than non-Bt cotton plants over time. Although Bt-susceptible larvae could not detect Bt toxin, they showed preference on non-Bt toxin substrates and were more likely to drop off substrates with Bt toxin.

Responses to Bt toxin Vip3Aa by pink bollworm larvae resistant or susceptible to Cry toxins

BACKGROUND

Transgenic crops that make insecticidal proteins from Bacillus thuringiensis (Bt) have revolutionized management of some pests. However, evolution of resistance to Bt toxins by pests diminishes the efficacy of Bt crops. Resistance to crystalline (Cry) Bt toxins has spurred adoption of crops genetically engineered to produce the Bt vegetative insecticidal protein Vip3Aa. Here we used laboratory diet bioassays to evaluate responses to Vip3Aa by pink bollworm (Pectinophora gossypiella), one of the world's most damaging pests of cotton.

RESULTS

Against pink bollworm larvae susceptible to Cry toxins, Vip3Aa was less potent than Cry1Ac or Cry2Ab. Conversely, Vip3Aa was more potent than Cry1Ac or Cry2Ab against laboratory strains highly resistant to those Cry toxins. Five Cry-susceptible field populations were less susceptible to Vip3Aa than a Cry-susceptible laboratory strain (APHIS-S). Relative to APHIS-S, significant resistance to Vip3Aa did not occur in strains selected in the laboratory for > 700-fold resistance to Cry1Ac or both Cry1Ac and Cry2Ab.

CONCLUSIONS

Resistance to Cry1Ac and Cry2Ab did not cause strong cross-resistance to Vip3Aa in pink bollworm, which is consistent with predictions based on the lack of shared midgut receptors between these toxins and previous results from other lepidopterans. Comparison of the Bt toxin concentration in plants relative to the median lethal concentration (LC₅₀ ) from bioassays may be useful for estimating efficacy. The moderate potency of Vip3Aa against Cry1Ac- and Cry2Ab-resistant and susceptible pink bollworm larvae suggests that Bt cotton producing this toxin together with novel Cry toxins might be useful as one component of integrated pest management. © 2022 Society of Chemical Industry.

Heterologous expression of cry1Ia12 insecticidal gene in cotton encodes resistance against pink bollworm, Pectinophora gossypiella (Lepidoptera: Gelechiidae); an alternate insecticidal gene for insect pest management

BACKGROUND

Cotton is continuously exposed to sucking and chewing insect pest pressure since emergence to harvesting. Pink bollworm (Pectinophora gossypiella) has become major chewing insect pest to reduce the cotton yield and results in bad lint quality even in transgenic crops. The efficiency of insecticidal genes has been compromised due to extensive utilization of transgenic crops.

METHODS AND RESULTS

The present study was conducted to evaluate the efficacy of an alternate cry1Ia12 insecticidal gene against pink bollworm (PBW) in cotton. Agrobacterium tumefaciens strain LBA4404 harboring pCAMBIA2300 expression vector containing cry1Ia12 gene under the control of 35S CaMV was used to transform a local cotton cultivar GS-01. The various molecular analyses revealed the transgene integration and expression in primary transformants. Among five selected transgenic plants, tcL-08 showed maximum (16.06-fold) mRNA expression of cry1Ia12 gene whereas tcL-03 showed minimum (2.33-fold) expression. Feeding bioassays of 2nd and 3rd instar pink bollworm (PBW) larvae on immature cotton bolls, flowers and cotton squares revealed up to 33.33% mortality on tcL-08 while lowest mortality (13.33%) was observed in tcL-03 and tcL-15. Furthermore, the average weight and size of survived larvae fed on transgenic plants was significantly lesser than the average weight of larvae survived on non-transgenic plants.

CONCLUSIONS

The present study suggests the cry1Ia12 gene as an alternate insecticidal gene for the resistance management of cotton bollworms, especially PBW.

Cry4Aa and Cry4Ba Mosquito-Active Toxins Utilize Different Domains in Binding to a Particular Culex ALP Isoform: A Functional Toxin Receptor Implicating Differential Actions on Target Larvae

The three-domain Cry4Aa toxin produced from Bacillus thuringiensis subsp. israelensis was previously shown to be much more toxic to Culex mosquito larvae than its closely related toxin—Cry4Ba. The interaction of these two individual toxins with target receptors on susceptible larval midgut cells is likely to be the critical determinant in their differential toxicity. Here, two full-length membrane-bound alkaline phosphatase (mALP) isoforms from Culex quinquefasciatus larvae, Cq-mALP1263and Cq-mALP1264, predicted to be GPI-linked was cloned and functionally expressed in Spodoptera frugiperda (Sf9) cells as 57- and 61-kDa membrane-bound proteins, respectively. Bioinformatics analysis disclosed that both Cq-mALP isoforms share significant sequence similarity to Aedes aegypti-mALP—a Cry4Ba toxin receptor. In cytotoxicity assays, Sf9 cells expressing Cq-mALP1264, but not Cq-mALP1263, showed remarkably greater susceptibility to Cry4Aa than Cry4Ba, while immunolocalization studies revealed that both toxins were capable of binding to each Cq-mALP expressed on the cell membrane surface. Molecular docking of the Cq-mALP1264-modeled structure with individual Cry4 toxins revealed that Cry4Aa could bind to Cq-mALP1264 primarily through particular residues on three surface-exposed loops in the receptor-binding domain—DII, including Thr⁵¹², Tyr⁵¹³ and Lys⁵¹⁴ in the β10-β11loop. Dissimilarly, Cry4Ba appeared to utilize only certain residues in its C-terminal domain—DIII to interact with such a Culex counterpart receptor. Ala-substitutions of selected β10-β11loop residues (T512A, Y513A and K514A) revealed that only the K514A mutant displayed a drastic decrease in biotoxicity against C. quinquefasciatus larvae. Further substitution of Lys⁵¹⁴ with Asp (K514D) revealed a further decrease in larval toxicity. Furthermore, in silico calculation of the binding affinity change (ΔΔG_bind) in Cry4Aa-Cq-mALP1264 interactions upon these single-substitutions revealed that the K514D mutation displayed the largest ΔΔG_bind value as compared to three other mutations, signifying an adverse impact of a negative charge at this critical receptor-binding position. Altogether, our present study has disclosed that these two related-Cry4 mosquito-active toxins conceivably exploited different domains in functional binding to the same Culex membrane-bound ALP isoform—Cq-mALP1264 for mediating differential toxicity against Culex target larvae.

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.

Critical Analysis of Multi-Omic Data from a Strain of Plutella xylostella Resistant to Bacillus thuringiensis Cry1Ac Toxin

Rapid evolution of resistance in crop pests to Bacillus thuringiensis (Bt) products threatens their widespread use, especially as pests appear to develop resistance through a range of different physiological adaptations. With such a diverse range of mechanisms reported, researchers have resorted to multi-omic approaches to understand the molecular basis of resistance. Such approaches generate a lot of data making it difficult to establish where causal links between physiological changes and resistance exist. In this study, a combination of RNA-Seq and iTRAQ was used with a strain of diamondback moth, Plutella xylostella (L.), whose resistance mechanism is well understood. While some of the causal molecular changes in the resistant strain were detected, other previously verified changes were not detected. We suggest that while multi-omic studies have use in validating a proposed resistance mechanism, they are of limited value in identifying such a mechanism in the first place.

Effect of genetically modified maize expressing the Cry1Ab and EPSPS proteins on growth, development, and gut bacterial diversity of the non-target arthropod Locusta migratoria

The widespread cultivation of genetically modified (GM) crops has raised concerns for their safety. Here, we evaluated the effects of a GM maize variety expressing the Cry1Ab (14.76 ± 0.87 μg/g FW) and EPSPS proteins (191.55 ± 15.69 μg/g FW) on the life-history traits and gut bacterial community of a non-target arthropod, Locusta migratoria, in the laboratory. We found that GM maize had no significant effect on the survival or body weight of different development stages of L. migratoria. The midgut and hindgut bacterial diversities and compositions were determined using high-throughput sequencing targeting the V3-V4 regions of the 16S rRNA. No significant changes were found in the species diversity or abundance between insects in the GM-fed treatment and the non-GM control. Furthermore, the concentration of Cry1Ab and EPSPS in the gut was determined after digestion of GM maize. Results showed that the contents of Cry1Ab/EPSPS rapidly decreased and were hard to detect after 72 h. Based on the parameters assessed, we can conclude that the GM maize variety examined has no significant adverse effect on L. migratoria.

Effect of Metabolic Changes in Aphis gossypii-Damaged Cotton Plants on Oviposition Preference and Larval Development of Subsequent Helicoverpa armigera

Aphis gossypii and Helicoverpa armigera are two important agricultural pests in cotton plants. However, whether early colonization of A. gossypii affects subsequent H. armigera is unknown. We implemented ecological experiments to reveal that A. gossypii-damaged cotton plants [Bacillus thuringiensis (Bt) and non-Bt] had a significant avoidance effect on the oviposition preference of H. armigera adults. However, A. gossypii-damaged cotton plants (non-Bt) increased the weight and pupation rate and reduced the mortality of H. armigera larvae. Transcriptomic and metabolomic analyses showed that 13 and 9 genes were significantly upregulated to be involved in salicylic acid (SA) and indole acetic acid (IAA) biosynthesis, and SA and IAA contents were significantly increased, respectively. However, 15 genes involved in jasmonic acid (JA) biosynthesis were significantly downregulated as a result of the antagonism of SA and JA. Moreover, there was significant upregulation in multiple genes involved in the biosynthesis of l-histidine, fructose, maltotetraose, melezitose, lecithin, stearidonic acid, and mannitol, in which metabolites were confirmed to promote the growth and development of H. armigera. Our study is a reference for investigating the evolutionary relationships and provides insights into implementing effective insect biocontrol between H. armigera and A. gossypii.

Food and feed safety of the Bacillus thuringiensis derived protein Vpb4Da2, a novel protein for control of western corn rootworm

Western corn rootworm (WCR), Diabrotica virgifera virgifera, LeConte, is an insect pest that poses a significant threat to the productivity of modern agriculture, causing significant economic and crop losses. The development of genetically modified (GM) crops expressing one or more proteins that confer tolerance to specific insect pests, such as WCR, was a historic breakthrough in agricultural biotechnology and continues to serve as an invaluable tool in pest management. Despite this, evolving resistance to existing insect control proteins expressed in current generation GM crops requires continued identification of new proteins with distinct modes of action while retaining targeted insecticidal efficacy. GM crops expressing insecticidal proteins must undergo extensive safety assessments prior to commercialization to ensure that they pose no increased risk to the health of humans or other animals relative to their non-GM conventional counterparts. As part of these safety evaluations, a weight of evidence approach is utilized to assess the safety of the expressed insecticidal proteins to evaluate any potential risk in the context of dietary exposure. This study describes the food and feed safety assessment of Vpb4Da2, a new Bacillus thuringiensis insecticidal protein that confers in planta tolerance to WCR. Vpb4Da2 exhibits structural and functional similarities to other insect control proteins expressed in commercialized GM crops. In addition, the lack of homology to known toxins or allergens, a lack of acute toxicity in mice, inactivation by conditions commonly experienced in the human gut or during cooking/food processing, and the extremely low expected dietary exposure to Vpb4Da2 provide a substantial weight of evidence to demonstrate that the Vpb4Da2 protein poses no indication of a risk to the health of humans or other animals.

MON 95379 Bt maize as a new tool to manage sugarcane borer (Diatraea saccharalis) in South America

BACKGROUND

The sugarcane borer (SCB), Diatraea saccharalis (Lepidoptera: Crambidae), is a key pest of maize in Argentina, and genetically modified maize, producing Bacillus thuringiensis (Bt) proteins, has revolutionized the management of this insect in South America. However, field-evolved resistance to some Bt technologies has been observed in SCB in Argentina. Here we assessed a new Bt technology, MON 95379, in the laboratory, greenhouse and field for efficacy against SCB.

RESULTS

In a laboratory leaf disc bioassay, both MON 95379 (producing Cry1B.868 and Cry1Da_7) and Cry1B.868_single maize (producing only Cry1B.868) resulted in 100% mortality of SCB. The level of Cry1B.868 in the Cry1B.868_single maize is comparable to that in MON 95379 maize. However, the Cry1Da_7 protein does not have high efficacy against SCB, as evidenced by < 20% mortality on Cry1Da_7_single leaf tissue. Total (100%) mortality of SCB in a Cry1B.868_single tissue dilution bioassay indicated that Cry1B.868_single maize meets the criteria to be classified as a high dose. Similar median lethal concentration (LC₅₀ ) values were observed for MON 89034-R and susceptible SCB strains exposed to Cry1B.868 protein. MON 95379 also controlled SCB strains resistant to MON 89034 (Cry1A.105/Cry2Ab2) and Cry1Ab. Under field conditions in Brazil and Argentina, MON 95379 maize plants were consistently protected from SCB damage.

CONCLUSION

MON 95379 maize will bring value to maize growers in South America by effectively managing SCB even in locations where resistance to other Bt-containing maize technologies has been reported. © 2022 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

MicroRNAs miR-14 and miR-2766 regulate tyrosine hydroxylase to control larval-pupal metamorphosis in Helicoverpa armigera

BACKGROUND

The cotton bollworm, Helicoverpa armigera, is a worldwide polyphagous pest, causing huge economic losses in vegetable, cotton and corn crops, among others. Owing to long-term exposure to Bacillus thuringiensis (Bt) toxins, evolution of resistance has been detected in this pest. As a conservative and effective neurotransmitter, dopamine (DA) has an important role in insect growth and development. In this study, we investigated the regulatory functions of DA and its associated non-coding RNA in metamorphosis in H. armigera.

RESULTS

Expression profiles indicated that DA and DA pathway genes were highly expressed during larval-pupal metamorphosis in H. armigera. RNA interference and pharmacological experiments confirmed that tyrosine hydroxylase (TH), dopa decarboxylase, vesicular amine transporter and DA receptor 2 are critical genes related to the development of H. armigera from larvae to pupae. We also found that miR-14 and miR-2766 targeted the 3' untranslated region to post-transcriptionally regulate HaTH function. Application of miR-2766 and miR-14 antagomirs significantly increased levels of HaTH transcripts and proteins, while injection of miR-2766 and miR-14 agomirs not only suppressed messenger RNA and protein levels of HaTH, but also resulted in defective pupation in H. armigera.

CONCLUSION

These results suggest that DA deficiency inhibits larval-pupal metamorphosis in H. armigera. Potentially, DA pathway genes and their microRNAs could be used as a novel target for H. armigera management. © 2022 Society of Chemical Industry.

Susceptibilities of the Invasive Fall Armyworm (Spodopterafrugiperda) to the Insecticidal Proteins of Bt maize in China

To control the fall armyworm (FAW), Spodoptera frugiperda (Smith), a serious threat to maize production in China, the Chinese government has issued biosafety certificates for transgenic insect-resistant maize expressing Bt (Bacillus thuringiensis) toxins including Bt-Cry1Ab maize (crop event DBN9936), Bt-Vip3Aa maize (event DBN9501), Bt-(Cry1Ab+Vip3Aa) maize with superimposed traits (event DBN9936 × DBN9501) and Bt-(Cry1Ab+Vip3Aa) maize with superimposed traits (event Bt11 × MIR162), but the susceptibility baselines of geographically distinct FAW populations to these events, which form the basis for managing resistance development in the pest to these events, are not clear. We used the diet-incorporated bioassays method to detect the susceptibilities of the seven FAW populations collected from Yunnan, Henan and Hubei provinces in China in 2021 to the insecticidal proteins of the four Bt maize events. The result showed that the susceptibilities of different geographical populations to Bt insecticidal proteins were significantly different. In the seven populations, the range in median lethal concentrations (LC₅₀) of Cry1Ab expressed in DBN9936 was 0.87-2.63 μg/g, 0.14-0.30 μg/g for Vip3Aa expressed in DBN9501, 0.78-1.86 μg/g for Cry1Ab+Vip3Aa expressed in DBN9936 × DBN9501, and 0.36-1.42 μg/g for CryAb+Vip3Aa expressed in Bt11 × MIR162. The growth inhibition responses also showed that the susceptibilities varied with the different median growth inhibitory concentration (GIC₅₀) ranges (0.38-1.22, 0.08-0.28, 0.28-0.87, and 0.24-0.78 μg/g, respectively). The variations in the ranges of the susceptibility baselines of the geographical populations of fall armyworm in China to the insecticidal proteins expressed in the four events provide a scientific basis for monitoring FAW population resistance to Bt maize and managing the populations using different Bt maize events.

Toxic effects of the combined cadmium and Cry1Ab protein exposure on the protective and transcriptomic responses of Pirata subpiraticus

Cadmium (Cd) pollution poses a serious threat to agricultural production and paddy field fauna. Crystalline proteins (e.g., Cry1Ab and Cry1Ac) are secreted by Bacillus thuringiensis, which can manage pests via a complicated toxic mechanism and have been widely used for pest control due to the commercialization of transgenic crops (e.g., cotton and rice) that expresses Bt insecticidal proteins. Nonetheless, studies on the effects of combined stress of Cd and Cry1Ab protein on field indicator species are limited. In the present study, we showed that spiders, Pirata subpiraticus, fed with Cd-containing flies+Cry1Ab had dramatically higher Cd accumulation than that in the spiders fed with Cd-containing flies (p < 0.05). In addition, the enrichment of Cd led to the activation of the protective mechanism by elevating the concentrations of glutathione peroxidase, glutathione S-transferase, and metallothionein in the spiders (p < 0.05). An in-depth transcriptome analysis revealed that the activities of ion metal binding proteins, transporters, and channels might play essential roles in the Cd accumulation process. More importantly, the higher Cd concentration in the combined Cd+Cry1Ab exposure prolonged developmental duration of P. subpiraticus, due to the down-regulated cuticle proteins (CPs) encoding genes involved in the molting process, which was regulated by a series of putative transcriptional factors such as ZBTB and zf-C2H2. Collectively, this integrated analysis illustrates that the combined Cd+Cry1Ab exposure increases the adverse effects of Cd stress on the growth, antioxidase, and CPs encoding genes of P. subpiraticus, thus providing a research basis and prospect for the rationality of transgenic Cry1Ab crops in the cultivation of heavy metal contaminated soil.

Molecular and Kinetic Models for Pore Formation of Bacillus thuringiensis Cry Toxin

Cry proteins from Bacillus thuringiensis (Bt) and other bacteria are pesticidal pore-forming toxins. Since 2010, when the ABC transporter C2 (ABCC2) was identified as a Cry1Ac protein resistant gene, our understanding of the mode of action of Cry protein has progressed substantially. ABCC2 mediates high Cry1A toxicity because of its high activity for helping pore formation. With the discovery of ABCC2, the classical killing model based on pore formation and osmotic lysis became nearly conclusive. Nevertheless, we are still far from a complete understanding of how Cry proteins form pores in the cell membrane through interactions with their host gut membrane proteins, known as receptors. Why does ABCC2 mediate pore formation with high efficiency unlike other Cry1A-binding proteins? Is the “prepore” formation indispensable for pore formation? What is the mechanism underlying the synergism between ABCC2 and the 12-cadherin domain protein? We examine potential mechanisms of pore formation via receptor interactions in this paper by merging findings from prior studies on the Cry mode of action before and after the discovery of ABC transporters as Cry protein receptors. We also attempt to explain Cry toxicity using Cry–receptor binding affinities, which successfully predicts actual Cry toxicity toward cultured cells coexpressing ABC transporters and cadherin.

RNA Interference-Mediated Knockdown of Bombyx mori Haemocyte-Specific Cathepsin L (Cat L)-Like Cysteine Protease Gene Increases Bacillus thuringiensis kurstaki Toxicity and Reproduction in Insect Cadavers

The silkworm’s Cat L-like gene, which encodes a lysosomal cathepsin L-like cysteine protease, is thought to be part of the insect’s innate immunity via an as-yet-undetermined mechanism. Assuming that the primary function of Cat L-like is microbial degradation in mature phagosomes, we hypothesise that the suppression of the Cat L-like gene expression would increase Bacillus thuringiensis (Bt) bacteraemia and toxicity in knockdown insects. Here, we performed a functional analysis of Cat L-like in larvae that were fed mulberry leaves contaminated with a commercial biopesticide formulation based on Bt kurstaki (Btk) (i.e., Dipel) to investigate its role in insect defence against a known entomopathogen. Exposure to sublethal doses of Dipel resulted in overexpression of the Cat L-like gene in insect haemolymph 24 and 48 h after exposure. RNA interference (RNAi)-mediated suppression of Cat L-like expression significantly increased the toxicity of Dipel to exposed larvae. Moreover, Btk replication was higher in RNAi insects, suggesting that Cat L-like cathepsin may be involved in a bacterial killing mechanism of haemocytes. Finally, our results confirm that Cat L-like protease is part of the antimicrobial defence of insects and suggest that it could be used as a target to increase the insecticidal efficacy of Bt-based biopesticides.

Downregulation of the CsABCC2 gene is associated with Cry1C resistance in the striped stem borer Chilo suppressalis

Chilo suppressalis is a major target pest of transgenic rice expressing the Bacillus thuringiensis (Bt) Cry1C toxin in China. The evolution of resistance of this pest is a major threat to Bt rice. Since Bt functions by binding to receptors in the midgut (MG) of target insects, identification of Bt functional receptors in C. suppressalis is crucial for evaluating potential resistance mechanisms and developing effective management strategies. ATP-binding cassette (ABC) transporters have been vastly reported to interact with Cry1A toxins, as receptors and their mutations cause insect Bt resistance. However, the role of ABC transporters in Cry1C resistance to C. suppressalis remains unknown. Here, we measured CsABCC2 expression in C. suppressalis Cry1C-resistant (Cry1C-R) and Cry1C-susceptible strains (selected in the laboratory) via quantitative real-time PCR (qRT-PCR); the transcript level of CsABCC2 in the Cry1C-R strain was significantly lower than that in the Cry1C-susceptible strain. Furthermore, silencing CsABCC2 in C. suppressalis via RNA interference (RNAi) significantly decreased Cry1C susceptibility. Overall, CsABCC2 participates in Cry1C mode of action, and reduced expression of CsABCC2 is functionally associated with Cry1C resistance in C. suppressalis.

Map-based cloning and functional analysis revealed ABCC2 is responsible for Cry1Ac toxin resistance in Bombyx mori

Bt toxins are parasporal crystals produced by Bacillus thuringiensis (Bt). They have specific killing activity against various insects and have been widely used to control agricultural pests. However, their widespread use has developed the resistance of many target insects. To maintain the sustainable use of Bt products, the resistance mechanism of insects to Bt toxins must be fully clarified. In this study, Bt-resistant and Bt-susceptible silkworm strains were used to construct genetic populations, and the genetic pattern of silkworm resistance to Cry1Ac toxin was determined. Sequence-tagged site molecular marker technology was used to finely map the resistance gene and to draw a molecular genetic linkage map, and the two closest markers were T1590 and T1581, indicating the resistance gene located in the 155 kb genetic region. After analyzing the sequence of the predicted gene in the genetic region, an ATP binding cassette transporter (ABCC2) was identified as the candidate gene. Molecular modeling and protein-protein docking result showed that a tyrosine insertion in the mutant ABCC2 might be responsible for the interaction between Cry1Ac and ABCC2. Moreover, CRISPR/Cas9-mediated genome editing technology was used to knockout ABCC2 gene. The homozygous mutant ABCC2 silkworm was resistant to Cry1Ac toxin, which indicated ABCC2 is the key gene that controls silkworm resistance to Cry1Ac toxin. The results have laid the foundation for elucidating the molecular resistance mechanism of silkworms to Cry1Ac toxin and could provide a theoretical basis for the biological control of lepidopteran pests.

[Advances in receptor-mediated resistance mechanisms of Lepidopteran insects to Bacillus thuringiensis toxin]

Bacillus thuringiensis is widely used as an insecticide which is safe and environmentally friendly to humans and animals. One of the important insecticidal mechanisms is the binding of Bt toxins to specific toxin receptors in insect midgut and forming a toxin perforation which eventually leads to insect death. The resistance of target pests to Bt toxins is an important factor hampering the long-term effective cultivation of Bt crops and the continuous use of Bt toxins. This review summarizes the mechanism of insect resistance to Bt toxins from the perspective of important Bt toxin receptors in midgut cells of Lepidopteran insects, which may facilitate the in-depth study of Bt resistance mechanism and pest control.

Role of Lectin in the Response of Aedes aegypti Against Bt Toxin

Aedes aegypti is one of the world’s most dangerous mosquitoes, and a vector of diseases such as dengue fever, chikungunya virus, yellow fever, and Zika virus disease. Currently, a major global challenge is the scarcity of antiviral medicine and vaccine for arboviruses. Bacillus thuringiensis var israelensis (Bti) toxins are used as biological mosquito control agents. Endotoxins, including Cry4Aa, Cry4Ba, Cry10Aa, Cry11Aa, and Cyt1Aa, are toxic to mosquitoes. Insect eradication by Cry toxin relies primarily on the interaction of cry toxins with key toxin receptors, such as aminopeptidase (APN), alkaline phosphatase (ALP), cadherin (CAD), and ATP-binding cassette transporters. The carbohydrate recognition domains (CRDs) of lectins and domains II and III of Cry toxins share similar structural folds, suggesting that midgut proteins, such as C-type lectins (CTLs), may interfere with interactions among Cry toxins and receptors by binding to both and alter Cry toxicity. In the present review, we summarize the functional role of C-type lectins in Ae. aegypti mosquitoes and the mechanism underlying the alteration of Cry toxin activity by CTLs. Furthermore, we outline future research directions on elucidating the Bti resistance mechanism. This study provides a basis for understanding Bti resistance, which can be used to develop novel insecticides.

The influence of Bt cotton cultivation on the structure and functions of the soil bacterial community by soil metagenomics

Bt cotton successfully controlled major devastating pests in cotton，such as Helicoverpa armigera and Spodoptera exigua, and led to a drastic decrease in insecticide use in cotton fields, and it has been grown commercially worldwide. However, Bt cotton cultivation left Bt toxin residues in the soil, resulting in a response by its microbiome that caused potential environmental risks. In this research, the metagenomics analysis was performed to investigate the structure and functions of the soil bacterial community in the Bt cotton field from the Binzhou, Shandong province of China, where the Bt cotton has been cultivated for over fifteen years. Analysis of the function genes proved that the receptors of Bt toxins were absent in the soil bacteria and Bt toxins failed to target the soil bacteria. The microbiome structure and function were highly influenced by Bt cotton cultivation, however, no significant change in the total abundance of the bacteria was observed. Proteobacteria was the largest taxonomic group in the soil bacterial (42-52%) and its abundance was significantly increased after Bt cotton cultivation. The increase of Proteobacteria abundance resulted in an increase in ABC transporters gene abundance, indicating the improved ability of detoxification metabolism over Bt cotton cultivation. Xanthomonadales could be a biomarker of the Bt cotton group, whose abundance was significantly increased to contribute to the increase of the genes abundance in ABC transporters. The abundance of apoptosis genes was significantly decreased, and it might be related to the increase of Proteobacteria abundance by Bt cotton cultivation. In addition, Myxococcales was responsible for carotenoid biosynthesis, whoes genes abundance was significantly decreased due to the decrease of Myxococcales abundance by Bt cotton cultivation. These changes in soil bacterial community structure and functions indicate the influence by Bt cotton cultivation, leading to an understanding of the bacteria colonization patterns due to successive years of Bt cotton cultivation. These research results should be significant for the rational risk assessment of Bt cotton cultivation.