Next Article in Journal
A Functional ClpXP Protease is Required for Induction of the Accessory Toxin Genes, tst, sed, and sec
Previous Article in Journal
Are Cyanotoxins the Only Toxic Compound Potentially Present in Microalgae Supplements? Results from a Study of Ecological and Non-Ecological Products
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxins
Volume 12
Issue 9
10.3390/toxins12090551
toxins-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Demographic Performance of Helicoverpa zea Populations on Dual and Triple-Gene Bt Cotton

by
Marcelo M. Rabelo
1,2,3,*,
Silvana V. Paula-Moraes
2,
Eliseu Jose G. Pereira
1,3 and
Blair D. Siegfried
4
1
Department of Entomology, Universidade Federal de Viçosa, Viçosa 36570-900, MG, Brazil
2
Department of Entomology and Nematology, West Florida Research and Education Center, University of Florida, Jay, FL 32565, USA
3
National Institute of Science and Technology in Plant-Pest Interactions, Universidade Federal de Viçosa, Viçosa 36570-900, MG, Brazil
4
Department of Entomology and Nematology, University of Florida, Gainesville, FL 32608, USA
*
Author to whom correspondence should be addressed.
Toxins 2020, 12(9), 551; https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12090551
Submission received: 2 July 2020 / Revised: 10 August 2020 / Accepted: 26 August 2020 / Published: 28 August 2020
(This article belongs to the Section Bacterial Toxins)
Browse Figure
Versions Notes

Abstract

:
Insecticidal toxins from Bacillus thuringiensis (Bt) are valuable tools for pest management worldwide, contributing to the management of human disease insect vectors and phytophagous insect pests of agriculture and forestry. Here, we report the effects of dual and triple Bt toxins expressed in transgenic cotton cultivars on the fitness and demographic performance of Helicoverpa zea (Boddie)—a noctuid pest, known as cotton bollworm and corn earworm. Life-history traits were determined for individuals of three field populations from a region where H. zea overwintering is likely. Triple-gene Bt cotton cultivars that express Cry and Vip3Aa toxins killed 100% of the larvae in all populations tested. In contrast, dual-gene Bt cotton that express Cry1Ac+Cry1F and Cry1Ac+Cry2Ab allowed population growth with the intrinsic rate of population growth (rm) 38% lower than on non-Bt cotton. The insects feeding on Bt cotton plants that express Cry1Ac+Cry2Ab, Cry1Ac+Cry1F, or Cry1Ab+Cry2Ae exhibited reduced larval weight, survival rate, and increased development time. Additionally, fitness parameters varied significantly among the insect populations, even on non-Bt cotton plants, likely because of their different genetic background and/or previous Bt toxin exposure. This is the first report of the comparative fitness of H. zea field populations on dual-gene Bt cotton after the recent reports of field resistance to certain Bt toxins. These results document the population growth rates of H. zea from an agricultural landscape with 100% Bt cotton cultivars. Our results will contribute to the development and validation of resistance management recommendations.
Keywords:
fitness; life table; cotton bollworm; corn earworm; toxin; resistance management
Key Contribution: Triple-gene Bt cotton cultivars expressing Cry and Vip3Aa toxins killed 100% of the larvae in all three populations tested. In contrast; dual-gene Bt cotton expressing Cry1Ac+Cry1F and Cry1Ac+Cry2Ab allowed population growth with the intrinsic rate of population growth (rm) 38% lower than on non-Bt cotton.

1. Introduction

Transgenic crops that express insecticidal toxins from the bacterium Bacillus thuringiensis (Berliner) (Bt) provide valuable pest management options for pests of field crops worldwide [1,2,3]. Positive socio-economic and environmental impacts of Bt crop adoption have been reported since commercial release in 1996 [4,5]. In the United States (U.S.), the Bt technology provides control of the major cotton pests Chloridea virescens (Fabricius) (Lepidoptera: Noctuidae), Pectinophora gossypiella (Saunders) (Lepidoptera: Gelichiidae), and Helicoverpa zea (Boddie) (Lepidoptera: Noctuidae) [6,7]. Bt cotton also improves the management of other lepidopteran pests, such as Spodoptera exigua (Hübner), Trichoplusia ni (Hübner), Spodoptera frugiperda (J.E. Smith), S. eridania (Stoll), and Chrysodeixis includens (Walker) [8,9,10,11]. From 1996–2003, commercial Bt cotton in the U.S. was limited to events that expressed the Cry1Ac toxin [12]. Second generation dual-gene Bt cotton that expressed the toxins Cry1Ac+Cry2Ab, Cry1Ac+Cry1F, and Cry1Ab+Cry2Ae became available in 2003 and were widely adopted [7,13]. The goal of these second-generation cotton events is to reduce the risk of resistance evolution by targeting unique and independent target sites. Since 2014 the third generation of cotton events with triple Bt traits became available, including those expressing Cry1Ac+Cry1F+Vip3Aa19, Cry1Ac+Cry2Ab+Vip3Aa19, and Cry1Ab+Cry2Ae+Vip3Aa19. The Cry and Vip families are produced during different stages of the B. thuringiensis life cycle [14]. Vip3A shares no sequence homology with any known Bt Cry toxins. The amount of shared sequence homology between Bt toxins is an important indicator of the risk that the two Bt toxins will share binding sites on the midgut of the insect, predisposing the toxins to cross-resistance. Although the two toxin classes are thought to have a similar mode of action against the target insects, they have different receptors in the insect midgut [14,15].
The tobacco budworm, C. virescens, is highly susceptible to most commercial Cry Bt toxins, and Bt cotton usually reaches a high-dose condition, killing almost all heterozygotes for Bt resistance [8,16], even for single toxin events. However, the cotton bollworm, H. zea is less susceptible to Cry toxins expressed in cotton and corn, which, therefore, do not satisfy high-dose criteria [17]. Toxicological bioassays performed with populations of H. zea from the southeastern U.S., including populations from the Florida Panhandle, have indicated a decrease in susceptibility of H. zea populations to Cry1Ab, Cry1Ac, Cry1A.105, and Cry2Ab, but not to Vip3Aa [5,18,19,20]. However, a major resistance allele conferring high levels of Vip3Aa resistance in a field-derived strain of H. zea in Texas has been recently reported [21].
Field-evolved resistance in target pests is a threat to the success of the Bt technology, which may lead to control failures [22] and the need for traditional insecticides for supplemental control [23]. Helicoverpa zea exhibits a sequence of host crop utilization based on the temporal dynamics of the southeastern U.S. agricultural landscape. In general, the first generation of this pest feeds on Bt field corn, cultivated during the spring to the beginning of summer. Corn is a major host plant of H. zea. However, this pest has a low impact on the yield of field corn when it is planted early in the season [24]. Later in the summer, subsequent generations of H. zea disperse from corn to cotton, which becomes the prevalent Bt crop in the southeastern U.S. agricultural landscape until the end of the crop season. Thus, corn serves as a source of H. zea populations, and if these source populations develop on Bt corn, a dispersion of Bt pre-exposed survivors from corn to cotton is likely. Larval feeding in both crops producing the very same or similar Bt toxin exerts continuous selection pressure and raises concerns about the selection of resistant populations [13].
The United States Environmental Protection Agency (EPA) has listed the knowledge of pest biology and ecology as key elements in formulating a Bt insect resistance management (IRM) programs [25]. While complete studies documenting the reproductive potential of H. zea populations feeding on dual- or triple-toxin Bt cotton are scarce, Cry1Ab Bt corn, for which is only moderately toxic, might reduce H. zea growth potential [26,27]. Developing life tables for H. zea on dual and triple-gene Bt cotton events allow the documentation of key fitness parameters, such as survival, development time, fertility, and population growth [28]. Fitness components and rates of population growth (i.e., demographic performance), if available, contribute to developing models to predict the rate of resistance evolution in target pests or to comparatively assess different resistance management practices, such as the use of structured and/or natural refuges [22,29,30,31].
Helicoverpa zea populations from the Florida Panhandle represent valuable resources for documenting the fitness components and demographic performance of lepidopteran pests targeted by Bt toxins in cotton. The region is in the Gulf Coastal Plain of the southeastern U.S., an ecological transition zone between temperate and subtropical climates. In the region, H. zea populations can overwinter and disperse throughout the growing season [32,33]. The overwintering survival of H. zea could be a carry-over source of Bt resistance alleles for other regions and between seasons [13]. Furthermore, the region has a distinctive regional landscape consisting of natural vegetation, forests, and field crops. Cotton is cultivated in large areas (approximately 50,000 hectares), with 100% adoption of Bt cultivars [34,35]. In this study, we report individual and population fitness of representative H. zea larvae challenged with dual- and triple-toxin Bt cotton technologies, information that contributes to the development and validation of resistance management recommendations.

2. Results

2.1. Life-History Traits

The interaction between cotton cultivar and insect population was significant (Table 1, p < 0.05) for larval weight, larval development time, larval survival, and pupal development time. Pupal weight, pre-pupa time, and egg viability varied only with the main effects of either cultivar or population or both (Table 1), and the pupal viability did not significantly vary (p > 0.05).
The dual-gene Bt cotton cultivars significantly (p < 0.05) reduced larval and pupal weights in all populations tested (Table 2). The population from Escambia County exhibited the lowest larval weight on Cry1Ab+Cry2Ae and Cry1Ac+Cry2Ab relative to non-Bt cotton, while Cry1Ac+Cry1F had the least negative impact on larval and pupal weights relative to the other cultivars. All three dual-gene Bt cotton cultivars reduced larval weight equally in the Santa Rosa population compared to the non-Bt cotton. The population from Jackson had the lowest larval weight on Cry1Ac+Cry2Ab and both lowest larval and pupal weights on Cry1Ab+Cry2Ae, while Cry1Ac+Cry1F cotton did not impact larval or pupal weights compared to non-Bt cotton. Among populations, H. zea from Jackson had the lowest larval and pupal weights, even when feeding on non-Bt cotton. Regarding larval survival rates, the triple-gene Bt cotton cultivars caused 100% mortality of all populations tested, and therefore, were not included in further analysis of life-history traits (Table 3). Helicoverpa zea from Escambia and Jackson had similar larval survival on non-Bt and Cry1Ac+Cry1F cotton. However, larval survival was reduced on Cry1Ac+Cry2Ab and Cry1Ab+Cry2Ae. Helicoverpa zea from Santa Rosa had larval survival reduced by all dual-gene Bt cotton, with Cry1Ab+Cry2Ae resulting in the most severe reduction. Among dual-gene Bt cotton cultivars, Cry1Ac+Cry1F and Cry1Ab+Cry2Ae allowed the highest and lowest larval survival, respectively. The survivorship of pupa (i.e., pupal viability) ranged from 91 to 100% and did not vary among cotton cultivars or insect populations (Table 3). The larvae developed more slowly on dual-gene Bt cotton than on non-Bt cotton, except for H. zea from Escambia and Santa Rosa feeding on Cry1Ab+Cry2Ae (Table 4). The insects from Jackson county had longer larval development time when feeding on non-Bt and Cry1Ab+Cry2Ae than the other populations. The Santa Rosa insects had longer pre-pupa development time on Cry1Ac+Cry2Ab, but this trait was not affected by the other cultivars, populations, or their interaction (Table 4). The duration of the pupal stage was shorter for insects feeding on Cry1Ac+Cry1F than on non-Bt or the other Bt cotton cultivars (Table 4). The egg viability was similar in all-cotton cultivars, but was higher for the Escambia population (Table 5).

2.2. Life Table Parameters

The demographic performance of H. zea feeding on non-Bt, Cry1Ac+Cry2Ab, and Cry1Ac+Cry1F varied among the cotton cultivars and insect populations (Table 6, Figure 1). Reproductive capacity on the other cultivars was not determined, due to low survival. The net reproductive rate (R0) of insects reared on non-Bt cotton was approximately 50% higher than those on Cry1Ac+Cry1F and Cry1Ac+Cry2Ab, except those from Jackson County, which exhibited the same R0 value on non-Bt and Cry1Ac+Cry1F. The intrinsic rate of population increase (rm) of the insects reared on non-Bt cotton was 30% greater than on Cry1Ac+Cry1F and Cry1Ac+Cry2Ab, except those from Jackson County, which exhibited the same rm value on non-Bt and Cry1Ac+Cry1F. The generation time (T) was nearly ten days shorter for insects reared on non-Bt cotton compared to those on the other cultivars, except for Jackson insects, which exhibited the same generation time on non-Bt and Cry1Ac+Cry1F. In contrast, the Jackson population had a higher fitness (higher R0, rm, and lower T) on Cry1Ac+Cry1F compared to the others.

3. Discussion

In this study, the life-history traits and demographic performance of H. zea from the Florida Panhandle varied among cotton cultivars and field populations, indicating differences among the cultivars in the efficacy against H. zea and the current population susceptibility to the Bt toxins. Gassmann et al. (2009) suggest that survival, developmental time, and body weight are key individual fitness components [30]. Here, the effects on immature insect fitness components associated with both population and cotton cultivar translated to negative effects on the growth potential of H. zea. Although the non-Bt cotton cultivar used as control is not isoline of the Bt cultivars, the differences in the effects on H. zea life-history between the Bt and non-Bt cotton documented here are likely associated with the expression of Bt toxins in each cotton cultivar tissue. Cotton plants are rich in terpenoid compounds, which may function as a barrier against herbivores impairing growth/development and/or behavioral traits. However, the cotton plant has been modified during domestication and breeding for high yield and quality, including low gossypol oil in cottonseeds, which may have lessened the content of anti-herbivory secondary metabolites (such as gossypol) [36,37,38]. Information on secondary compounds in the cotton cultivars used in the present study was not available. However, other studies comparing non-Bt cotton cultivars have shown low or no change in noctuid life-history [39].
The triple-gene Bt cotton that expressed Vip3Aa (Cry1Ac+Cry1F+Vip3Aa19, Cry1Ac+Cry2Ab+Vip3Aa19, or Cry1Ab+Cry2Ae+Vip3Aa19) caused 100% larval mortality in all H. zea populations, which reinforces the high efficacy of this toxin for H. zea control [18,22]. These results suggest a low frequency of resistant alleles to Vip3A in H. zea populations tested. Data from laboratory and field in the U.S. consistently indicate high efficacy of the Vip3A against H. zea [6,40,41], and the debate has been focused on whether the trait meets the high dose definition [42]. Due to the relatively recent adoption of Vip3A toxins in commercial cultivars and limited insect sample size (40-130 individuals) in the present study, it would be unlikely to detect resistance to Vip3A at its current low frequency in the field. Vip toxins show limited amino-acid sequence homology with Cry toxins and cause pore formation with unique properties, thus, having a low risk for cross-resistance between them [14,43,44]. Cry toxins co-expressed in some Bt cotton cultivars have reportedly low impact in some H. zea populations [19,45], which compromises the pyramid of Bt genes. Our study demonstrates that life-history traits (body weight, survival, and development time) of insects from Escambia, Santa Rosa, and Jackson populations were negatively affected by the dual-gene Bt cotton cultivars that expressed Cry toxins.
Life-history traits were more affected by Cry1Ab+Cry2Ae and Cry1Ac+Cry2Ab, while Cry1Ac+Cry1F caused fewer negative impacts given the relative lack of toxicity that Cry1F has on the H. zea larvae and the widespread Cry1Ac resistance [19]. Significant mortality from Cry1Ac+Cry1F was observed only in the Santa Rosa population. Cry toxins have been expressed in Bt cotton cultivars since its first commercial release, and the first report of H. zea Cry1Ac resistance was documented 15 years later in the U.S., and recently the widespread resistance to Cry2Ab [6,20,46]. These may be a contributing factor to the considerable rates of larval survival on Cry1Ac+Cry2Ab and Cry1Ab+Cry2Ae cultivars. It also confirms the Cry1A and possibly Cry2A resistance alleles occurrence at high frequencies in the H. zea populations tested, which does not mean that the dual gene Bt cotton cultivars lost the benefit on H. zea management completely. Overall, our data on life-history traits (survival rates, body weight, developmental time) are consistent with previous reports that Cry1Ac+Cry1F affects H. zea larvae less than Cry1Ac+Cry2Ab, Cry1Ab+Cry2Ae, and cultivars expressing Vip3Aa [5,7,19].
Sublethal effects of Bt toxins on H. zea, as indicated by reduced body weight and the prolonged larval development, may have implications for pest management. Delayed larval development and low body weight are expected to increase the likelihood of exposure to other mortality factors. For example, early-instar larvae are unable to bore into the cotton bolls [47]. Consequently, they may be more exposed to insecticide applications and vulnerable to natural enemies [48]. Moreover, slow larval growth tends to increase the intervals for insecticide applications, which should target the most vulnerable stage of smaller larvae (about 1 cm) [49].
Cotton is the last summer crop to be planted in the Florida Panhandle region, remaining for a longer period than other crops in the agricultural landscape before the fallow season. The longer larval development time of H. zea when feeding on Cry1Ac+Cry2Ab, associated with infestations during mid- and late season, could expose larvae to shorter days and decreasing temperatures, factors that regulate insect diapause [50,51]. The Florida Panhandle is considered a “hybrid zone” of populations of noctuids, such as S. frugiperda, which flies from south Florida and Texas to the northern U.S. [52,53]. Diapausing and migration of H. zea populations from the Florida Panhandle may contribute to infestations in corn and cotton, North to 40 N latitude, where H. zea cannot permanently survive [32].
The pupal viability of H. zea was similar among different cotton cultivars and populations. In contrast, H. zea pupal weight varied when feeding on different cotton cultivars. Pupal weight is often correlated with fecundity [54], although this correlation might be affected by several other factors [55]. In our study, the heavier pupal weight was linked with higher fecundity, which agrees with reports for other noctuids, such as H. armigera [56].
Cumulative effects on specific life-history traits of H. zea (i.e., larval survival and development time) impact the population growth potential on Bt cotton cultivars. The life table parameters indicated that H. zea populations tested are expected to grow when feeding on Cry1Ac+Cry1F or Cry1Ac+Cry2Ab Bt cotton, but with reduced growth rates (i.e., R0, rm). Overall, insects feeding on Cry1Ac+Cry1F and Cry1Ac+Cry2Ab are expected to generate 30% and 42% fewer individuals per day compared to non-Bt cotton, respectively. The growth potential of H. zea from Jackson county population was similar when feeding on non-Bt or on Cry1Ac+Cry1F, which could result in a higher number of exposed offspring [12]. However, the insects of the Jackson county had lower growth rates on non-Bt cotton compared with Escambia and Santa Rosa populations, indicating the presence of fitness costs [57] when they do not feed on Cry1Ac+Cry1F cotton. In a theoretical scenario where only Cry1Ac+Cry1F cotton is cultivated in the Florida Panhandle, the H. zea population from Jackson county is expected to produce in one generation 10–27% more females per female than the populations from Escambia and Santa Rosa. These differences between H. zea populations reinforce that resistance may develop because of local selection [58,59].
During a period of over 23 years, in which commercialized Bt crops have been used in the U.S., IRM programs have relied on models to predict how quickly resistance to Bt may occur in different scenarios [60,61]. Although our study was performed in the laboratory and conclusions about field-evolved resistance are limited, this is one of the most complete life table studies of H. zea in Bt cotton available in the literature. The information provided in this study (life-history traits and life table parameters) can contribute to the refinement of predictive models and delayed resistance to important Bt toxins, such as Vip3Aa [22]. Our results reinforce the need for region-specific knowledge of target pests of Bt technology when designing IRM programs [60]. Helicoverpa zea has a high dispersal capacity and reproductive biology, which leads to extensive gene flow [60,61]. However, fitness components and their variability across environments should be taken into consideration in the simulation of predictive models [60].
In conclusion, this paper has quantified the dual and triple-gene Bt cotton effect on the life-history and demographic performance of three populations of H. zea from the Florida Panhandle. Triple-gene Bt cotton caused 100% larval mortality in all populations tested, indicating the value of Vip3Aa toxin on H. zea management in the region. Despite resistance, dual-gene Bt cotton containing Cry1A and Cry2A toxins significantly affected the fitness and demographic growth of the three populations of H. zea evaluated. However, the magnitude of the effect on the life-history, and consequently on the life table parameters of H. zea in a landscape containing 100% Bt cotton varied. Interaction between H. zea populations (Escambia, Santa Rosa, and Jackson counties) and cotton cultivars (Bt and non-Bt) was detected. These findings improve our understanding of how data on demographic growth rates of target pests to Bt technology matters and fill a gap by providing region-specific information when developing IRM programs. The results of this study also provide valuable parameters for the refinement of models to better predict the risk of resistance evolution and validate resistance management strategies, including refuge recommendations.

4. Materials and Methods

4.1. Cotton Plants

This study was conducted during 2018 at the West Florida Research and Education Center (WFREC), the University of Florida at Jay, FL, USA. The cotton cultivars utilized are adapted to the region and described in Table 7. The cultivars were planted in a Randomized Complete Block Design with four replications. Each cotton cultivar was planted on 5 m wide × 8 m long plots containing eight rows. The agronomic practices were based on standard recommendations for the region [62], and no applications of insecticides were performed in the experimental plots. Fully expanded cotton leaves were collected from the upper part of the plant canopy in each plot during the first bloom to open boll plant stages, placed in a ziplock bag (Johnson, Racine, WI, USA), and held in Styrofoam (ULINE, Chicago, FL, USA) boxes with an ice pack. In the laboratory, the cotton leaves were tested using Envirologix GMO quick Stix to confirm Cry1Ac, Cry2Ab, Cry1F, and Vip3Aa expression (EnviroLogix Kit, Portland, ME, USA) among the different events. Similar procedures were used to collect blooms, squares, and bolls at the early-middle stages of development also during the first bloom to open boll plant stages.

4.2. Insect Populations

Three H. zea populations were collected during the 2018 crop season from commercial fields located in the main cotton-producing counties in the Florida Panhandle: Santa Rosa, Escambia, and Jackson, USA. Cotton fields were located in areas where a peanut/cotton rotation is prevalent, with corn planted on a smaller scale. The populations from Santa Rosa (n = 100) and Jackson (n = 130) were collected from ears of Bt corn (Cry1A.105+Cry2Ab). The Escambia population (n = 40) was collected from blooms and bolls of Bt cotton (Cry1Ac+Cry2Ab). Information on the collections, including location, and the number of generations in the laboratory are shown in Table 8.
Collected larvae were identified based on their morphology and validated after adult emergence [63]. The larvae were removed from the plant individually and placed in plastic cups containing a multispecies lepidopteran diet (Southland Products, Lake Village, AR, USA). The cups were held in Styrofoam boxes with an ice pack during transport to the laboratory where they were maintained at 25 ± 2 °C, 70 ± 10% relative humidity and 14L: 10D photoperiod. Pupae were transferred to Petri dishes and covered with vermiculite moistened with water and placed in rearing cages (22 × 30 × 2.5 cm) for adult emergence. The adults were fed a solution of 10% honey which was replaced every two days. Paper towels (Great Value, Bentonville, AR, USA) were used to cover the internal walls of the cages as an oviposition substrate. The eggs were collected and transferred to ziplock bags until hatching. Neonates were transferred to a multispecies lepidopteran diet (Southland Products, Lake Village, AR, USA) in rearing containers (Southland Products, Lake Village, AR, USA) and maintained individually until pupation.

4.3. Life-History Traits and Life Table Parameters

One hundred H. zea neonates of each population (Santa Rosa, Escambia, and Jackson (Table 8)) were transferred in groups of five to 473 mL polypropylene containers (Fabri-Kal Corp. Kalamazoo, MI, USA) and fed with cotton leaves, blooms, squares, and bolls of the cultivars described in Table 7. After 5 d, the larvae were placed in individual containers to avoid cannibalism, as previously described [64]. The plant tissues were replaced every four days until pupation. When the larvae reached fourth instar, wet vermiculite was added to the bottom of the rearing containers as a substrate for pupation and to avoid desiccation. The pupae were left in the containers until adult emergence. Larvae weight was determined after seven days. Once larval development was completed, and within 24 h after pupation, each pupa was weighed, and the sex was determined. Other life-history components were recorded, including survival rate (neonate to pupa) and development time of larvae, pre-pupae, pupae, and adults. The experiment was arranged in a completely randomized design with 100 larvae per cotton cultivar (1 larva per replication) for each population.
One male and female from each cultivar that emerged within two days of one another were confined in mating cages (30 cm high × 20 cm diameter polymerized vinyl chloride tube). The cages were covered with a waxed brown paper (Roberts Consolidated Industries Inc., Boca Raton, FL, USA) as an oviposition substrate, and supplied with a 10% aqueous honey solution, replaced every day. Adult survival and the number of eggs were recorded daily. The brown paper containing the eggs was transferred to ziplock bags until hatching. An additional egg viability estimation was performed based on daily evaluation of the presence of the neonates in each ziplock bag. The sex ratio, number of eggs (fecundity), survival, and age of females at the onset of egg-laying were determined to estimate the life table parameters. These included the net reproductive rate (R0), which represents the multiplication rate per generation, the intrinsic rate of population increase (rm), which reflects the ability of one female to generate another female per unit of time, and generation time (T), the mean time between two successive generations. The life table experiment was conducted in a completely randomized design, with 16 to 18 replications (couples) per cotton cultivar.

4.4. Statistical Analyses

Differences in the survival rate, body weight, development time, and egg viability of the three populations reared on the seven cotton cultivars were compared using a two-way analysis of variance in R software (version 3.5.1) [65]. The fixed effects tested were the H. zea population, cotton cultivar, and their interaction. Pairwise comparisons were made using Tukey’s HSD (honestly significant difference) post hoc test using a level of significance of 0.05. The population growth parameters (R0, rm, T) were determined using the SAS programming developed by Maia (2000) [66], and the variances associated with the estimates were obtained by the Jackknife method. This procedure constructs confidence intervals for the estimated parameters in addition to comparisons by the t-test.

Author Contributions

Conceptualization, M.M.R. and S.V.P.-M.; investigation, M.M.R.; resources, S.V.P.-M.; data curation, M.M.R.; writing—original draft preparation, M.M.R.; writing—review and editing, M.M.R., S.V.P.-M., E.J.G.P., B.D.S.; visualization, S.V.P.-M., E.J.G.P., B.D.S.; supervision, S.V.P.-M. and E.J.G.P.; project administration, S.V.P.-M.; funding acquisition, S.V.P.-M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by USDA National Institute of Food and Agriculture, Hatch project 005643, and CPPM EIP FLA-ENY-005649.

Acknowledgments

We thank the extension agent Ethan Carter, the Entomology team at WFREC for the help collecting the insects, and Barry Brecke for the English proofreading. We also acknowledge the farmers from the Florida Panhandle who provided access to their commercial fields and the companies Bayer, BASF, and Corteva for providing the cotton seeds. MMR received a scholarship from CAPES (Brazilian Ministry of Education, Finance code 001).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Naranjo, S.E. Impacts of Bt transgenic cotton on integrated pest management. J. Agric. Food Chem. 2011, 59, 5842–5851. [Google Scholar] [CrossRef]
  2. Bates, S.L.; Zhao, J.Z.; Roush, R.T.; Shelton, A.M. Insect resistance management in GM crops: Past, present and future. Nat. Biotechnol. 2005, 23, 57–62. [Google Scholar] [CrossRef]
  3. Koch, M.S.; Ward, J.M.; Levine, S.L.; Baum, J.A.; Vicini, J.L.; Hammond, B.G. The food and environmental safety of Bt crops. Front. Plant Sci. 2015, 6, 283. [Google Scholar] [CrossRef]
  4. ISAAA. Global Status of Commercialized Biotech/GM Crops in 2017: Biotech Crop Adoption Surges as Economic Benefits Accumulate in 22 Years; ISAAA: Kenya, PH, USA, 2017; pp. 1–153. [Google Scholar]
  5. Fleming, D.; Musser, F.; Reisig, D.; Greene, J.; Taylor, S.; Parajulee, M.; Lorenz, G.; Catchot, A.; Gore, J.; Kerns, D.; et al. Effects of transgenic Bacillus thuringiensis cotton on insecticide use, heliothine counts, plant damage, and cotton yield: A meta-analysis, 1996–2015. PLoS ONE 2018, 13, e0200131. [Google Scholar] [CrossRef]
  6. Braswell, L.R.; Reisig, D.D.; Sorenson, C.E.; Collins, G.D. Development and dispersal of Helicoverpa zea (Lepidoptera: Noctuidae) on non-Bt and Bt pyramided cotton. Environ. Entomol. 2019, 48, 465–477. [Google Scholar] [CrossRef]
  7. Luttrell, R.G.; Jackson, R.E. Helicoverpa zea and Bt cotton in the United States. GM Crops Food 2012, 3, 213–227. [Google Scholar] [CrossRef] [Green Version]
  8. Fox, J.L. Bt cotton infestations renew resistance concerns. Nat. Biotechnol. 1996, 14, 1070. [Google Scholar] [CrossRef]
  9. Rabelo, M.M.; Matos, J.M.L.; Santos-Amaya, O.F.; França, J.C.; Gonçalves, J.; Paula-Moraes, S.V.; Guedes, R.N.C.; Pereira, E.J.G. Bt-toxin susceptibility and hormesis-like response in the invasive southern armyworm (Spodoptera eridania). Crop Prot. 2020, 132, 105129. [Google Scholar] [CrossRef]
  10. Stewart, S.D.; Adamczyk, J.J.; Knighten, K.S.; Davis, F.M. Impact of bt cottons expressing one or two insecticidal proteins of Bacillus thuringiensis Berliner on growth and survival of noctuid (Lepidoptera) larvae. J. Econ. Entomol. 2009, 94, 752–760. [Google Scholar] [CrossRef]
  11. Sorgatto, R.J.; Bernardi, O.; Omoto, C. Survival and development of Spodoptera frugiperda and Chrysodeixis Includens (Lepidoptera: Noctuidae) on Bt cotton and implications for resistance management strategies in Brazil. Environ. Entomol. 2015, 44, 189–192. [Google Scholar] [CrossRef]
  12. Tabashnik, B.E.; Gassmann, A.J.; Crowder, D.W.; Carriére, Y. Insect resistance to Bt crops: Evidence versus theory. Nat. Biotechnol. 2008, 26, 199–202. [Google Scholar] [CrossRef]
  13. US EPA Insect Resistance Management for Bt Plant-Incorporated Protectants. Available online: https://www.epa.gov/regulation-biotechnology-under-tsca-and-fifra/insect-resistance-management-bt-plant-incorporated (accessed on 5 April 2020).
  14. Palma, L.; Muñoz, D.; Berry, C.; Murillo, J.; Caballero, P.; Caballero, P. Bacillus thuringiensis toxins: An overview of their biocidal activity. Toxins 2014, 6, 3296–3325. [Google Scholar] [CrossRef] [Green Version]
  15. Pardo-López, L.; Soberón, M.; Bravo, A. Bacillus thuringiensis insecticidal three-domain Cry toxins: Mode of action, insect resistance and consequences for crop protection. FEMS Microbiol. Rev. 2013, 37, 3–22. [Google Scholar] [CrossRef] [Green Version]
  16. Blanco, C.A. Heliothis virescens and Bt cotton in the United States. GM Crops Food 2012, 3, 201–212. [Google Scholar] [CrossRef] [Green Version]
  17. Stone, T.B.; Sims, S.R. Geographic susceptibility of Heliothis virescens and Helicoverpa zea (Lepidoptera: Noctuidae) to Bacillus thuringiensis. J. Econ. Entomol. 1993, 86, 989–994. [Google Scholar] [CrossRef]
  18. Leite, N.A.; Pereira, R.M.; Durigan, M.R.; Amado, D.; Fatoretto, J.; Medeiros, F.C.L.; Omoto, C. Susceptibility of Brazilian populations of Helicoverpa armigera and Helicoverpa zea (Lepidoptera: Noctuidae) to Vip3Aa20. J. Econ. Entomol. 2018, 111, 399–404. [Google Scholar] [CrossRef]
  19. Reisig, D.D.; Huseth, A.S.; Bacheler, J.S.; Aghaee, M.A.; Braswell, L.; Burrack, H.J.; Flanders, K.; Greene, J.K.; Herbert, D.A.; Jacobson, A.; et al. Long-term empirical and observational evidence of practical Helicoverpa zea resistance to cotton with pyramided bt toxins. J. Econ. Entomol. 2018, 111, 1824–1833. [Google Scholar] [CrossRef]
  20. Kaur, G.; Guo, J.; Brown, S.; Head, G.P.; Price, P.A.; Paula-Moraes, S.; Ni, X.; Dimase, M.; Huang, F. Field-evolved resistance of Helicoverpa zea (Boddie) to transgenic maize expressing pyramided Cry1A.105/Cry2Ab2 proteins in northeast Louisiana, the United States. J. Invertebr. Pathol. 2019, 163, 11–20. [Google Scholar] [CrossRef]
  21. Yang, F.; González, J.C.S.; Little, N.; Reisig, D.; Payne, G.; Dos Santos, R.F.; Jurat-Fuentes, J.L.; Kurtz, R.; Kerns, D.L. First documentation of major Vip3Aa resistance alleles in field populations of Helicoverpa zea (Boddie) (Lepidoptera: Noctuidae) in Texas, USA. Sci. Rep. 2020, 10, 1–8. [Google Scholar]
  22. Reisig, D.D.; Kurtz, R. Bt resistance implications for Helicoverpa zea (Lepidoptera: Noctuidae) insecticide resistance management in the United States. Environ. Entomol. 2018, 47, 1357–1364. [Google Scholar] [CrossRef] [Green Version]
  23. Reisig, D.; Kerns, D.; Gore, J.; Musser, F. Managing pyrethroid-and Bt-resistant bollworm in southern U.S. cotton. Crop. Soils 2019, 52, 30–35. [Google Scholar] [CrossRef] [Green Version]
  24. Bibb, J.L.; Cook, D.; Catchot, A.; Musser, F.; Stewart, S.D.; Leonard, B.R.; Buntin, G.D.; Kerns, D.; Allen, T.W.; Gore, J.; et al. Impact of corn earworm (Lepidoptera: Noctuidae) on field corn (Poales: Poaceae) yield and grain quality. J. Econ. Entomol. 2018, 111, 1249–1255. [Google Scholar] [CrossRef] [PubMed]
  25. US EPA. Bt plant-incorporated protectants October 15, 2001 biopesticides registration action document. Biopestic. Regist. Action Doc. 2001. Available online: https://www3.epa.gov/pesticides/chem_search/reg_actions/pip/bt_brad.htm (accessed on 22 August 2020).
  26. Storer, N.P.; van Duyn, J.W.; Kennedy, G.G. Life history traits of Helicoverpa zea (Lepidoptera: Noctuidae) on non-Bt and Bt transgenic corn hybrids in Eastern North Carolina. J. Econ. Entomol. 2001, 94, 1268–1279. [Google Scholar] [CrossRef]
  27. Reisig, D.D.; Reay-Jones, F.P.F. Inhibition of Helicoverpa zea (Lepidoptera: Noctuidae) growth by transgenic corn expressing bt toxins and development of resistance to Cry1Ab. Environ. Entomol. 2015, 44, 1275–1285. [Google Scholar] [CrossRef]
  28. Birch, L.C. The intrinsic rate of natural increase of an insect population. J. Anim. Ecol. 1948, 17, 15. [Google Scholar] [CrossRef]
  29. Hackett, S.C.; Bonsall, M.B. Type of fitness cost influences the rate of evolution of resistance to transgenic Bt crops. J. Appl. Ecol. 2016, 53, 1391–1401. [Google Scholar] [CrossRef] [Green Version]
  30. Gassmann, A.J.; Carrière, Y.; Tabashnik, B.E. Fitness costs of insect resistance to Bacillus thuringiensis. Annu. Rev. Entomol. 2009, 54, 147–163. [Google Scholar] [CrossRef]
  31. Pittendrigh, B.R.; Gaffney, P.J.; Huesing, J.E.; Onstad, D.W.; Roush, R.T.; Murdock, L.L. Active refuges can inhibit the evolution of resistance in insects towards transgenic insect-resistant plants. J. Theor. Biol. 2004, 231, 461–474. [Google Scholar] [CrossRef]
  32. Morey, A.C.; Hutchison, W.D.; Venette, R.C.; Burkness, E.C. Cold hardiness of Helicoverpa zea (Lepidoptera: Noctuidae) Pupae. Environ. Entomol. 2012, 41, 172–179. [Google Scholar] [CrossRef]
  33. Raulston, J.R.; Pair, S.D.; Pedraza Martinez, F.A.; Westbrook, J.; Sparks, A.N.; Sanchez Valdez, V.M. Ecological studies indicating the migration of Heliothis zea, Spodoptera frugiperda, and Heliothis virescens from Northeastern Mexico and Texas. Proc. Life Sci. 1986, 204–220. [Google Scholar] [CrossRef]
  34. Cook, D.R. National cotton council of America. In Proceedings of the 2017 Cotton Insect Loss Estimates; 2018 Beltwide Cotton Conferences, San Antonio, TX, USA, 5 January 2018; pp. 721–780. [Google Scholar]
  35. USDA/NASS National Agricultural Statistics Service-USDA’s Florida Field Office. Available online: https://www.nass.usda.gov/Statistics_by_State/Florida/Publications/County_Estimates/ (accessed on 5 April 2020).
  36. Ti, X.; Zhang, Q. Advances in research of induced resistance to insects in cotton. Front. Biol. China 2009, 4, 289–297. [Google Scholar] [CrossRef]
  37. Chen, Y.H.; Gols, R.; Benrey, B. Crop domestication and its impact on naturally selected trophic interactions. Annu. Rev. Entomol. 2015, 60, 35–58. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Bernal, J.S.; Medina, R.F. Agriculture sows pests: How crop domestication, host shifts, and agricultural intensification can create insect pests from herbivores. Curr. Opin. Insect Sci. 2018, 26, 76–81. [Google Scholar] [CrossRef]
  39. Rabelo, M.M.; Matos, J.M.L.; Orozco-Restrepo, S.M.; Paula-Moraes, S.V.; Pereira, E.J.G. Like parents, like offspring? Susceptibility to Bt toxins, development on dual-gene Bt cotton, and parental effect of Cry1Ac on a non-target lepidopteran pest. J. Econ. Entomol. 2020, 113, 1234–1242. [Google Scholar] [CrossRef]
  40. Bilbo, T.R.; Reay-Jones, F.P.F.; Reisig, D.D.; Greene, J.K.; Turnbull, M.W. Development, survival, and feeding behavior of Helicoverpa zea (Lepidoptera: Noctuidae) relative to Bt protein concentrations in corn ear tissues. PLoS ONE 2019, 14, e0221343. [Google Scholar] [CrossRef] [Green Version]
  41. Ali, M.I.; Luttrell, R.G. Susceptibility of Helicoverpa zea and Heliothis virescens (Lepidoptera: Noctuidae) to Vip3A insecticidal protein expressed in VipCotTM cotton. J. Invertebr. Pathol. 2011, 108, 76–84. [Google Scholar] [CrossRef] [PubMed]
  42. US EPA. Resistance of lepidopteran pests to Bacillus thuringiensis (Bt) plant incorporated protectants (PIPs) in The United States. FIFRA Sci. 2018. Available online: https://cfpub.epa.gov/si/si_public_record_report.cfm?Lab=OPP&dirEntryID=342052 (accessed on 15 October 2020).
  43. Dourado, P.M.; Bacalhau, F.B.; Amado, D.; Carvalho, R.A.; Martinelli, S.; Head, G.P.; Omoto, C. High susceptibility to Cry1Ac and low resistance allele frequency reduce the risk of resistance of Helicoverpa armigers to Bt soybean in Brazil. PLoS ONE 2016, 11, e0161388. [Google Scholar] [CrossRef] [Green Version]
  44. Jurat-Fuentes, J.L.; Crickmore, N. Specificity determinants for cry insecticidal proteins: Insights from their mode of action. J. Invertebr. Pathol. 2017, 142, 5–10. [Google Scholar] [CrossRef] [Green Version]
  45. Dively, G.P.; Venugopal, P.D.; Finkenbinder, C. Field-evolved resistance in corn earworm to cry proteins expressed by transgenic sweet corn. PLoS ONE 2016, 11, e0169115. [Google Scholar] [CrossRef]
  46. Sivasupramaniam, S.; Moar, W.J.; Ruschke, L.G.; Osborn, J.A.; Jiang, C.; Sebaugh, J.L.; Brown, G.R.; Shappley, Z.W.; Oppenhuizen, M.E.; Mullins, J.W.; et al. Toxicity and characterization of cotton expressing Bacillus thuringiensis Cry1Ac and Cry2Ab2 proteins for control of lepidopteran pests. J. Econ. Entomol. 2008, 101, 546–554. [Google Scholar] [CrossRef]
  47. Ellsworth, P.C.; Bradley, J.R. Comparative damage potential and feeding dynamics of the European corn borer (Lepidoptera: Pyralidae) and cotton bollworm (Lepidoptera: Noctuidae) on cotton bolls. J. Econ. Entomol. 1992, 85, 402–410. [Google Scholar] [CrossRef]
  48. Liu, X.; Chen, M.; Collins, H.L.; Onstad, D.W.; Roush, R.T.; Zhang, Q.; Earle, E.D.; Shelton, A.M. Natural enemies delay insect resistance to Bt crops. PLoS ONE 2014, 9, e90366. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Pedigo, L.P. Entomology and pest management. In Entomology and Pest Management; Waveland Press: Long Grove, IL, USA, 1996; pp. 505–526. ISBN 0130195677. [Google Scholar]
  50. Meola, R.; Gray, R. Temperature-sensitive mechanism regulating diapause in Heliothis zea. J. Insect Physiol. 1984, 30, 743–749. [Google Scholar] [CrossRef]
  51. Clemmensen, S.F.; Hahn, D.A. Dormancy cues alter insect temperature–size relationships. Oecologia 2015, 177, 113–121. [Google Scholar] [CrossRef] [Green Version]
  52. Nagoshi, R.N.; Meagher, R.L.; Hay-Roe, M. Inferring the annual migration patterns of fall armyworm (Lepidoptera: Noctuidae) in the United States from mitochondrial haplotypes. Ecol. Evol. 2012, 2, 1458–1467. [Google Scholar] [CrossRef]
  53. Nagoshi, R.N.; Meagher, R.L.; Hay-Roe, M. Assessing the resolution of haplotype distributions to delineate fall armyworm (Lepidoptera: Noctuidae) migratory behaviors. J. Econ. Entomol. 2014, 107, 1462–1470. [Google Scholar] [CrossRef] [Green Version]
  54. Leuck, D.B.; Perkins, W.D. A method of estimating fall armyworm progeny reduction when evaluating control achieved by host-plant resistance12. J. Econ. Entomol. 1972, 65, 482–483. [Google Scholar] [CrossRef]
  55. Bilbo, T.R.; Reay-Jones, F.P.F.; Reisig, D.D.; Musser, F.R.; Greene, J.K. Effects of Bt corn on the development and fecundity of corn earworm (Lepidoptera: Noctuidae). J. Econ. Entomol. 2018, 111, 2233–2241. [Google Scholar] [CrossRef]
  56. Reigada, C.; Guimarães, K.F.; Parra, J.R.P. Relative fitness of Helicoverpa armigera (Lepidoptera: Noctuidae) on seven host plants: A perspective for IPM in Brazil. J. Insect Sci. 2016, 16, 3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Jakka, S.R.K.; Knight, V.R.; Jurat-Fuentes, J.L. Fitness costs associated with field-evolved resistance to Bt maize in Spodoptera frugiperda (Lepidoptera: Noctuidae). J. Econ. Entomol. 2014, 107, 342–351. [Google Scholar] [CrossRef]
  58. Andow, D.A. The risk of resistance evolution in insects to transgenic insecticidal crops. Collect. Biosaf. Rev. 2008, 4, 142–199. [Google Scholar]
  59. Endersby, N.M.; Ridland, P.M.; Hoffmann, A.A. The effects of local selection versus dispersal on insecticide resistance patterns: Longitudinal evidence from diamondback moth (Plutella xylostella (Lepidoptera: Plutellidae)) in Australia evolving resistance to pyrethroids. Bull. Entomol. Res. 2008, 98, 145. [Google Scholar] [CrossRef] [PubMed]
  60. Kennedy, G.G.; Storer, N.P. Life systems of polyphagous arthropod pests in temporally unstable cropping systems. Annu. Rev. Entomol. 2000, 45, 467–493. [Google Scholar] [CrossRef] [PubMed]
  61. Caprio, M.A. Evaluating resistance management strategies for multiple toxins in the presence of external refuges. J. Econ. Entomol. 1998, 91, 1021–1031. [Google Scholar] [CrossRef]
  62. Wright, D.; Marois, J.; Rich, J. Cotton cultural practices and fertility management. Agron. Dep. Florida Coop. Ext. Serv. Inst. Food Agric. Sci. Univ. Florida. P 2011, 20, 1–5. [Google Scholar]
  63. Hardwick, D.F. The corn earworm complex. Mem. Entomol. Soc. Can. 1965, 97, 5–247. [Google Scholar] [CrossRef]
  64. Capinera, J.L. Corn Earworm Scientific Name: Helicoverpa zea (Boddie) (Insecta: Lepidoptera: Noctuidae). July 2000. Reviewed: June 2020. Future Creature. University of Florida. Available online: http://entnemdept.ufl.edu/creatures/veg/corn_earworm.htm (accessed on 2 May 2020).
  65. R Core Team. R: A Language and Environment For Statistical Computing; Vienna, Austria, 2013; Available online: https://www.R-project.org/ (accessed on 23 April 2020).
  66. Maia, A.d.H.N.; Luiz, A.J.B.; Campanhola, C. Statistical inference on associated fertility life table parameters using jackknife technique: Computational aspects. J. Econ. Entomol. 2000, 93, 511–518. [Google Scholar] [CrossRef]
Toxins 12 00551 g001 550
Figure 1. Reproductive schedule of H. zea feeding on non-Bt and Bt cotton cultivars as represented by fecundity (number of eggs per day) and female longevity. Each line represents an average of 16 H. zea females mated in pairs in mating cages. Panels A, B, and C represent the populations from Escambia, Santa Rosa, and Jackson county, respectively. The black line refers to insects feeding on non-Bt cotton, while the red and green are for insects feeding on Bt cotton Cry1Ac+Cry1F and Cry1Ac+Cry2Ab, respectively.
Figure 1. Reproductive schedule of H. zea feeding on non-Bt and Bt cotton cultivars as represented by fecundity (number of eggs per day) and female longevity. Each line represents an average of 16 H. zea females mated in pairs in mating cages. Panels A, B, and C represent the populations from Escambia, Santa Rosa, and Jackson county, respectively. The black line refers to insects feeding on non-Bt cotton, while the red and green are for insects feeding on Bt cotton Cry1Ac+Cry1F and Cry1Ac+Cry2Ab, respectively.
Toxins 12 00551 g001
Table
Table 1. Two-way ANOVA for life-history traits of Helicoverpa zea populations feeding on cotton cultivars.
Table 1. Two-way ANOVA for life-history traits of Helicoverpa zea populations feeding on cotton cultivars.
VariableSource of VariationFp
Larval weightPopulation8.8670.0003
Cultivar75.853<0.0001
Population × Cultivar8.471<0.0001
Larval development timePopulation12.959<0.0001
Cultivar174.10<0.0001
Population × Cultivar4.9270.0002
Larval survivalPopulation7.1360.0010
Cultivar182.51<0.0001
Population × Cultivar4.529<0.0001
Pre-pupal development timePopulation0.990.3717
Cultivar5.180.0016
Population × Cultivar1.7150.1302
Pupal weightPopulation9.430.0001
Cultivar15.70<0.0001
Population × Cultivar1.580.1647
Pupal development timePopulation0.9330.3943
Cultivar31.435<0.0001
Population × Cultivar4.846<0.0001
Pupal survivalPopulation0.2710.764
Cultivar0.6910.561
Population × Cultivar0.5200.760
Egg viabilityPopulation3.7120.0275
Cultivar0.1340.8751
Population × Cultivar0.1290.9427
p values of 0.05 or lower were considered significant as calculated using two-way ANOVA in R software (version 3.5.1).
Table
Table 2. Larval and pupal weight (mg) of H. zea reared on Bt and non-Bt cotton cultivars.
Table 2. Larval and pupal weight (mg) of H. zea reared on Bt and non-Bt cotton cultivars.
StagePopulationCotton Cultivar
Non-BtCry1Ac+Cry1Ac+Cry1Ab+Cry1Ac+Cry1Ac+Cry1Ab+
Cry1FCry2AbCry2AeCry1F+Cry2Ab+Cry2Ae+
Vip3AaVip3AaVip3Aa
LarvaEscambia59.6 ± 24.4 Aa25.5 ± 6.18 Ba2.02 ± 2.14 Ca0.29 ± 0.14 Ca***
Santa Rosa56.8 ± 16.0 Aa6.79 ± 1.96 Bb2.31 ± 1.91 Ba0.34 ± 0.17 Ba***
Jackson20.1 ± 7.98 Ab22.6 ± 6.03 Aa0.26 ± 0.02 Bb1.64 ± 1.66 Ba***
PupaEscambia429.0 ± 69.3 Aa363.0 ± 49.7 Ca393.0 ± 49.0 Ba401.0 ± 00.0 Ba***
Santa Rosa427.0 ± 67.7 Aa386.0 ± 67.9 Ba409.0 ± 40.4 Aa324.0 ± 00.0 Ba***
Jackson391.0 ± 70.7 Ab373.0 ± 53.2 Aa*302.0 ± 80.8 Ba***
Means (± SE) followed by the same capital letter within lines or the same lowercase latter within columns for each parameter do not significantly differ (p > 0.05; Tukey HSD). * not determined due to the high larval mortality.
Table
Table 3. Larva and pupal survival rates (%) of H. zea reared on Bt and non-Bt cotton cultivars.
Table 3. Larva and pupal survival rates (%) of H. zea reared on Bt and non-Bt cotton cultivars.
StagePopulationCotton Cultivar
non-BtCry1Ac+Cry1Ac+Cry1Ab+Cry1Ac+Cry1Ac+Cry1Ab+
Cry1FCry2AbCry2AeCry1F+Cry2Ab+Cry2Ae+
Vip3AaVip3AaVip3Aa
LarvaEscambia93 ± 4.83 Aa80 ± 0.7 Aa30.0 ± 28.3 Ba2.0 ± 0.4 Ca0.0 ± 0.0 Ca0.0 ± 0.0 Ca0.0 ± 0.0 Ca
Santa Rosa83 ± 10.6 Aab40 ± 1.0 Bb18.0 ± 16.2 Bab8.0 ± 0.8 Ca0.0 ± 0.0 Ca0.0 ± 0.0 Ca0.0 ± 0.0 Ca
Jackson77 ± 16.4 Ab80 ± 1.5 Aa4.22± 2.00 Bb13 ± 16.4 Ba0.0 ± 0.0 Ca0.0 ± 0.0 Ca0.0 ± 0.0 Ca
PupaEscambia96.9 ± 6.5 Aa93.0 ± 11.4 Aa95.2 ± 12.6 Aa100 ± 0.0 Aa***
Santa Rosa100 ± 0.0 Aa96.3 ± 11.1 Aa91.7 ± 20.0 Aa100 ± 0.0 Aa***
Jackson100 ± 0.0 Aa95.9 ± 10.8 Aa*100 ± 0.0 Aa***
Means (± SE) followed by the same capital letter within lines or the same lowercase latter within columns for each parameter do not significantly differ (p > 0.05; Tukey HSD). * not determined due to the high larval mortality.
Table
Table 4. Development time (days) of H. zea reared on Bt and non-Bt cotton cultivars.
Table 4. Development time (days) of H. zea reared on Bt and non-Bt cotton cultivars.
StagePopulationCotton Cultivar
Non-BtCry1Ac+Cry1Ac+Cry1Ab+Cry1Ac+Cry1Ac+Cry1Ab+
Cry1FCry2AbCry2AeCry1F+Cry2Ab+Cry2Ae+
Vip3AaVip3AaVip3Aa
LarvaEscambia21.7 ± 1.43 Cb28.3 ± 3.22 Ba32.5 ± 4.79 Aa21.00 ± 0.00 Cb***
Santa Rosa21.9 ± 2.02 Cb28.7 ± 2.89 Ba32.9 ± 3.51 Aa20.00 ± 3.50 Cb***
Jackson24.3 ± 3.63 Ca28.0 ± 3.18 Ba*37.2 ± 3.50 Aa***
Pre-PupaEscambia3.38 ± 0.71 Aa3.45 ± 1.03 Aa3.75 ± 1.08 Aa3.00 ± 0.00 Aa***
Santa Rosa3.52 ± 0.89 Ba3.27 ± 0.84 Ba4.36 ± 0.80 Aa3.00 ± 0.00 Ba***
Jackson3.48 ± 1.02 Ba3.75 ± 0.91 Ba*2.75 ± 0.50 Aa***
PupaEscambia19.0 ± 1.05 Aa17.4 ± 2.00 Ba19.2 ± 1.18 Aa21.00 ± 0.00 Aa***
Santa Rosa18.4 ± 1.44 Ba17.2 ± 2.78 Ca20.9 ± 1.14 Aa20.00 ± 1.81 Aa***
Jackson19.3 ± 2.04 Aa17.4 ± 2.00 Ba*17.00 ± 0.00 Ab***
Means (± SE) followed by the same capital letter within lines or the same lowercase latter within columns for each parameter do not significantly differ (p > 0.05; Tukey HSD). * not determined due to the high larval mortality.
Table
Table 5. Egg viability (%) of H. zea reared on Bt and non-Bt cotton cultivars.
Table 5. Egg viability (%) of H. zea reared on Bt and non-Bt cotton cultivars.
PopulationCotton Cultivar
on-BtCry1Ac+Cry1Ac+Cry1Ab+Cry1Ac+Cry1Ac+Cry1Ab+
Cry1FCry2AbCry2AeCry1F+Cry2Ab+Cry2Ae+
Vip3AaVip3AaVip3Aa
Escambia65.2 ± 37.6 Aa59.0 ± 44.3 Aa69.2 ± 34.2 Aa****
Santa Rosa44.3 ± 39.5 Ab44.5 ± 33.9 Ab43.8 ± 39.8 Ab****
Jackson43.3 ± 32.6 Ab46.7 ± 36.6 Ab43.3 ± 32.6 Ab****
Means (± SE) followed by the same lowercase latter within columns do not significantly differ (p > 0.05; Tukey HSD). * not determined due to the high larval mortality.
Table
Table 6. Life table of H. zea populations from different counties in the Florida Panhandle reared on Bt and non-Bt cotton cultivars.
Table 6. Life table of H. zea populations from different counties in the Florida Panhandle reared on Bt and non-Bt cotton cultivars.
ParameterPopulationCotton Cultivar
Non-BtCry1Ac+Cry1Ac+Cry1Ab+Cry1Ac+Cry1Ac+Cry1Ab+
Cry1FCry2AbCry2AeCry1F+Cry2Ab+Cry2Ae+
Vip3AaVip3AaVip3Aa
R0Escambia320.62 ± 61.74 Aa139.32 ± 18.99 Bb133.54 ± 16.94 Ba****
Santa Rosa406.80 ± 67.24 Aa87.26 ± 12.71 Bc80.57 ± 20.34 Ba****
Jackson289.59 ± 53.32 Aa270.09 ± 35.50 Aa*****
rmEscambia0.13 ± 0.006 Aa0.10 ± 0.003 Bb0.08 ± 0.003 Ca****
Santa Rosa0.13 ± 0.006 Aa0.08 ± 0.003 Bc0.07 ± 0.005 Bb****
Jackson0.11 ± 0.05 Ab0.11 ± 0.03 Aa*****
TEscambia41.94 ± 0.67 Ab48.88 ± 0.65 Bb55.37 ± 0.65 Cb****
Santa Rosa45.08 ± 0.95 Aa52.99 ± 0.87 Ba60.34 ± 1.73 Ca****
Jackson47.62 ± 1.20 Aa49.57 ± 0.75 Ab*****
Means (± SE) followed by the same capital letter within lines or the same lowercase latter within columns for each parameter do not significantly differ (p >  0.05) through pairwise comparisons using two-tailed t-tests after the jackknife method to estimate variance. R0—intrinsic rate of population increase (females per female per generation); rm—net reproductive rate (females per female per day); T—generation time (days). * not determined due to the high larval mortality.
Table
Table 7. Non-Bt and Bt cotton cultivars, expressing Bacillus thuringiensis toxins, that were used in this study.
Table 7. Non-Bt and Bt cotton cultivars, expressing Bacillus thuringiensis toxins, that were used in this study.
Trade NameCultivarFamilyBt Event NameYear LaunchedBt Toxin
Non-BtDP 1822XFDeltapine---
Bollgard IIDP 1646B2XFDeltapineMON159852003Cry1Ac, Cry2Ab
WideStrikePHY 444WRFPhytogen3006-210-23, 281-24-2362005Cry1Ac, Cry1F
TwinLinkST 5122GLTStonevilleT304-40, GHB1192014Cry1Ab, Cry2Ae
Bollgard IIIDP 1851B3XFDeltapineMON15985, COT1022014Cry1Ac, Cry2Ab, Vip3Aa19
WideStrike 3PHY480W3FEPhytogen3006-210-23, 281-24-236, COT1022015Cry1Ac, Cry1F, Vip3Aa
TwinLink PlusST 5471GLTPStonevilleT304-40 × GHB119 × COT1022017Cry1Ab, Cry2Ae, Vip3Aa19
Table
Table 8. Helicoverpa zea populations from the Florida Panhandle, 2018 crop season.
Table 8. Helicoverpa zea populations from the Florida Panhandle, 2018 crop season.
CountyGeospatial CoordinateNumber of Insects
Collected		Generation Tested
LatitudeLongitude
Santa Rosa30.7757−87.1432100F3
Escambia30.9842−87.469640F2
Jackson30.8041−85.0805130F3

Share and Cite

MDPI and ACS Style

Rabelo, M.M.; Paula-Moraes, S.V.; Pereira, E.J.G.; Siegfried, B.D. Demographic Performance of Helicoverpa zea Populations on Dual and Triple-Gene Bt Cotton. Toxins 2020, 12, 551. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12090551

AMA Style

Rabelo MM, Paula-Moraes SV, Pereira EJG, Siegfried BD. Demographic Performance of Helicoverpa zea Populations on Dual and Triple-Gene Bt Cotton. Toxins. 2020; 12(9):551. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12090551

Chicago/Turabian Style

Rabelo, Marcelo M., Silvana V. Paula-Moraes, Eliseu Jose G. Pereira, and Blair D. Siegfried. 2020. "Demographic Performance of Helicoverpa zea Populations on Dual and Triple-Gene Bt Cotton" Toxins 12, no. 9: 551. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12090551

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop