Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Sixteen Years of Bt Maize in the EU Hotspot: Why Has Resistance Not Evolved?

  • Pedro Castañera,

    Affiliation Department of Environmental Biology, Centro de Investigaciones Biológicas, CSIC, Ramiro de Maeztu 9, 28040 Madrid, Spain

  • Gema P. Farinós,

    Affiliation Department of Environmental Biology, Centro de Investigaciones Biológicas, CSIC, Ramiro de Maeztu 9, 28040 Madrid, Spain

  • Félix Ortego,

    Affiliation Department of Environmental Biology, Centro de Investigaciones Biológicas, CSIC, Ramiro de Maeztu 9, 28040 Madrid, Spain

  • David A. Andow

    dandow@umn.edu

    Affiliation Department of Entomology, University of Minnesota, 219 Hodson Hall, 1980 Folwell Ave., Saint Paul, Minnesota 55108, United States of America

Correction

25 Jul 2016: Castañera P, Farinós GP, Ortego F, Andow DA (2016) Correction: Sixteen Years of Bt Maize in the EU Hotspot: Why Has Resistance Not Evolved?. PLOS ONE 11(7): e0160294. https://doi.org/10.1371/journal.pone.0160294 View correction

Abstract

The majority of Bt maize production in the European Union (EU) is concentrated in northeast Spain, which is Europe’s only hotspot where resistance might evolve, and the main target pest, Sesamia nonagrioides, has been exposed to Cry1Ab maize continuously since 1998. The cropping system in northeast Spain has some similar characteristics to those that probably led to rapid resistance failures in two other target noctuid maize pests. These include repeated cultivation of Bt maize in the same fields, low use of refuges, recurring exposure of larvae to non-high dose concentrations of Cry1Ab toxin during the first years of cultivation, low migratory potential, and production concentrated in an irrigated region with few alternative hosts. Available data reveal no evidence of resistance in S. nonagrioides after 16 years of use. We explore the possible reasons for this resistance management success using evolutionary models to consider factors expected to accelerate resistance, and those expected to delay resistance. Low initial adoption rates and the EU policy decision to replace Event 176 with MON 810 Bt maize were key to delaying resistance evolution. Model results suggest that if refuge compliance continues at the present 90%, Bt maize might be used sustainably in northeast Spain for at least 20 more years before resistance might occur. However, obtaining good estimates of the present R allele frequency and level of local assortative mating are crucial to reduce uncertainty about the future success of resistance management.

Introduction

The approval in 1998 of Cry1Ab-expressing maize hybrids (Bt maize) in the European Union (EU) continues to be debated strongly today. Spain is the only EU country where Bt maize has been grown continuously at a large scale, accounting for 92% of the total EU area [1], with most of this in Ebro Valley in northeast Spain, to control recurring damage by the target pest, Mediterranean corn borer, Sesamia nonagrioides Lef. [Lepidoptera: Noctuidae] [2]. Adoption of Bt maize in northeast Spain has reached >75% (Fig 1), which makes this Europe’s only hotspot, as defined by the European Food Safety Authority (EFSA) [3], where resistance might be expected.

thumbnail
Fig 1. Adoption rate of Bt maize in the three main maize growing areas in Spain.

Percentage of Bt maize in northeast (Cataluña and Aragón), central (Madrid and Castilla-La Mancha) and southwest (Extremadura and Andalucía) Spain with respect to the total cultivated maize in each region (inner graph).

https://doi.org/10.1371/journal.pone.0154200.g001

Resistance in the target pests is the major threat to the long-term sustainability of insect resistant transgenic plants. Field-evolved resistance to Bt crops has been reported in only five insect pests during the last 17 years [4]. For Bt maize, resistance has been documented in two species of noctuid Lepidoptera: African stem borer, Busseola fusca, to Bt maize expressing Cry1Ab in South Africa [5, 6] and fall armyworm, Spodoptera frugiperda, to Cry1F-expressing maize in Puerto Rico and Brazil [7, 8, 9], which occurred after only 8 and 3 years of cultivation, respectively. Resistance can be ascribed to several causes, but factors common to South Africa and Brazil [5, 6, 8, 10, 11] were a) the high adoption rate and low refuge use, b) repeated cultivation of Bt maize in the same fields, c) with recurring exposure of larvae of the second and third generation on the crop to non-high dose levels of toxin due to reduced expression of the Bt toxin, d) low migratory potential, and e) production concentrated in an irrigated region where there were few wild or cultivated alternative hosts. In addition, for B. fusca in the Vaalharts area, there was increased selection pressure due to its preference for irrigated Bt maize rather than for adjacent non-irrigated non-Bt fields. For S. frugiperda, its preference for non-damaged maize [12] might have increased selection pressure as it might have preferentially oviposited on undamaged Bt maize instead of damaged non-Bt maize.

Some of these same factors are common to S. nonagrioides populations in northeast Spain, where nearly all maize is irrigated and concentrated in a small area and crop rotation is limited. Refuge compliance was low during the first eight years (<60% before 2006) of use of Cry1Ab-expressing Bt maize (see Table H in S1 File). Compa CB (Event 176) was the only variety grown until 2003, when it was gradually replaced by MON 810 hybrids and disallowed in 2006. Whereas MON 810 varieties express the Cry1Ab toxin in vegetative and reproductive stages at similar levels [13], the toxin titer for Compa CB decreased as the crop reached maturation [14], exposing larvae of the second and third generation S. nonagrioides to non-high dose concentrations of Cry1Ab toxin [15, 16], which would accelerate resistance by reducing selection against heterozygotes. Similarly, long distance migration in this species may be rare [17], which allows local selection to dominate the evolution of resistance. In addition, S. nonagrioides is functionally monophagous on maize in the region. Thus, there are some factors that would suggest that the resistance risk for S. nonagrioides could be high. At the same time, there are other factors that would suggest a reduced resistance risk for this species, including low initial adoption of Bt maize, and elimination of Event 176. See S1 File (Supporting Information) for more details on the biology of S. nonagrioides and maize cultivation in northeast Spain (Figs A-C in S1 File).

Here we establish that resistance has not been detected in S. nonagrioides in northeast Spain after 16 years of Bt maize cultivation, by measuring susceptibility in this hotspot and two other locations with historically low levels of adoption of Bt maize (Fig 1) and by comparing with information from EFSA [18]. Next we investigate possible reasons why resistance not been detected despite high exposure to a single Bt maize toxin. We considered factors expected to accelerate resistance (e.g., low refuge compliance, low female movement, and high local positive assortative mating), and those expected to delay resistance (e.g., low adoption, external colonization, and replacement of Event 176) using a resistance evolution model, parameterized for S. nonagrioides.

Materials and Methods

No specific permissions were required at any of the study locations for all of our activities. In Spain, permission for collections is not required because the species is a maize pest. Furthermore, the field sampling was part of the Spanish Post-Market Resistance Monitoring Programme for Bt-maize cultivation, which is required and coordinated by the Spanish Ministerio de Agricultura, Alimentación y Medio Ambiente (MAAMA). As a part of their responsibility, the MAAMA also informed farmers that samples were taken and explained the purpose of the study. This study did not involve endangered or protected species.

Bioassays

A unique batch of Cry1Ab crystals (81% purity), isolated from a bacterial culture of B. thuringiensis ssp. kurstaki HD1-9 strain and provided by Syngenta in 1999 [16], was used for all the bioassays. The Cry1Ab crystals were stored at -20°C and every year a toxin stock solution was prepared by suspending the crystals in a 0.1% solution of Triton X-100 to perform the bioassays programmed for that season.

Last instar pre-diapausing larvae of S. nonagrioides from the final generation of the season were collected from maize fields from northeast, central and southwest Spain (2–4 sampling sites from each geographical area, between 100–200 larvae from each site making a total of about 500 larvae, see S2 File for the geographic coordinates for the sampling sites) and the susceptibility to Cry1Ab was tested on the F1 progeny. In addition, a laboratory colony established in 1998, before Bt maize cultivation, and refreshed every two years with individuals collected from non-Bt fields (CIB, Madrid, Spain) was used as control. Bioassays were performed in “Bio-Ba-128” plastic trays (Color-Dec Italy, Capezzano Pianore, Italy) containing 128 wells, where ≈ 0.5 ml of rearing diet [16] were placed, flattened and let solidify. Serial dilutions (7 to 10 different concentrations) of Cry1Ab toxin that elicited insect responses ranging from near absent to near complete mortality, excluding dilutions with 0% and 100% responses, plus a control were tested. Fifty μl of each concentration was added in each well on the surface of the diet (1.767 cm2). The controls consisted of the buffer solution used to dilute the toxin. The treated trays were let dry in a laminar flow hood and shortly afterwards one neonate larva (<24 h old) was placed in each well and confined with a cover (Bio-Cv-16, Color-Dec Italy, Capezzano Pianore, Italy). Three replicates of 16–32 larvae at each concentration and 32–64 larvae for the control were performed. Neonates from different oviposition cages were used in each replicate. Trays were incubated in growth chambers at 25 ± 1°C, 70 ± 5% r.h. and constant dark. Larval mortality (larvae that did not show any reaction when prodded) was recorded after 7 days. The maximum mortality allowed in controls was 20%. Lethal concentrations (LC50 and LC90) values with a 95% confidence interval were estimated by probit analysis using POLO-PC (LeOra Software). Resistance Ratios (RR) with 95% confidence limits were calculated with respect to the LC50 or MIC50 (molt inhibition concentration) values obtained for the susceptible laboratory strain each year to normalize the data [19].

Mathematical Model

The resistance evolution model is based on Comins (1977) [20,21,22, 23], and is a single-locus, two-allele, deterministic, three-patch, discrete generation model with stochastic weather that determines the diapause rate and survival of the third generation (Fig D in S1 File). Each generation is modelled by tracking population density and allele frequencies through adult emergence, pre-mating movement, mating, post-mating movement, reproduction, selection, density-dependent larval survival, and adult emergence. Each generation has a unique parameterization based on the ecology of S. nonagrioides. The full model is a system of six coupled non-linear difference equations characterizing resistance allele frequency and population size in three different habitat patches (see SI for mathematical details). The three habitat patches are Mon810 Bt maize, Event 176 Bt maize, and non-Bt maize. For each generation there are 24 parameters, and most of the parameters are allowed to vary over the three generations per year. Nearly all of the refuge plants in northeast Spain are non-Bt maize. Resistance is assumed to be conferred by a single resistance allele that gives resistance to both Mon810 and Event 176. Selection occurs on larvae prior to density-dependent mortality. Adult moths disperse and mate. Dispersal is allowed to vary among males, virgin females and mated females. Mating can be locally assortative or disassortative.

Parameter Estimation

All parameters were estimated independently, based on the best information available on the ecology and genetics of S. nonagrioides, and the maize cropping system in northeast Spain, and the set of parameter values is referred to as the Best Parameter Values (BPV). For example, the proportion of non-Bt maize refuge is taken from historical records of planting of Bt maize in northeast Spain (Fig 1). Future projections assume that adoption of Bt maize rises up to the limit allowed by the refuge policy in Spain and the European Union (80% of the area planted to maize in northeast Spain). The initial resistance allele frequency in S. nonagrioides was estimated by using an F2 screen and a Bayesian estimator [24], resulting in a best available estimate of initial resistance allele frequency of 1.5 x 10−3. Selection will be different in Event 176 and MON 810, because toxin concentration in Event 176 declines after maize anthesis [14]. Selection parameters are equivalent to a high dose for the first generation on Event 176 and for all generations on Mon810, and heterozygosity increases for each successive generation on Event 176. The best parameter values (BPV) were used to simulate the most likely evolutionary trajectory in northeast Spain, and a range of parameters were used to investigate possible alternative trajectories (Table 1) [2538]. Additional details on these and all other parameter values are in the S1 File (Tables B-I and Figs E-G in S1 File).

thumbnail
Table 1. Parameter values used in resistance evolution model.

https://doi.org/10.1371/journal.pone.0154200.t001

Model Experiments

Although the model structure and parameterization were designed to provide the best predictions of the time to resistance, models are more useful for understanding how the relative time to resistance changes as parameter values change. Because of their heuristic nature, these kinds of investigations of model behaviour can be considered to be thought experiments, answering questions like “What if the parameters were different?” We varied the following parameter values comparing the simulation results to the BPV model: initial R allele frequency, female dispersal, Bt maize adoption rate, spatial structure of the refuge (and degree of assortative mating), external colonization and interactions among the parameters (see SI for additional details).

All simulations began with a “warm-up” period to remove artifacts from initial density transients. First, populations of S. nonagrioides were simulated on landscapes that had no Bt maize with the initial resistance allele frequency until the population densities reached equilibrium (no selection). This also generated an equilibrium distribution of the R allele (equal to the initial frequency in all patch types). After the populations reached equilibrium densities, Bt maize was introduced using the varieties and rates actually experienced in northeast Spain during the past 17 years. In addition, based on the actual weather data from northeast Spain, we estimated a diapause rate for the second generation and an end-of-the-season mortality rate for the third generation for the past 17 years (Table F in S1 File). Simulations continued until the average resistance allele frequency across all patch types exceeded 0.5, which was recorded as the number of years to resistance. Because success of the third generation was stochastic, we ran 1000 replications of each parameter combination and report the mean and standard deviations of these simulations.

Results and Discussion

Sustained efficacy of Bt maize for 16 years

Available data revealed no shift in susceptibility to Cry1Ab in Spanish field populations of S. nonagrioides after 16 years (32–48 generations) of cultivation of Bt maize in this hotspot. A post-market resistance monitoring program was initiated in 1998, allowing us to establish a baseline and to measure the natural variations in susceptibility of field populations of S. nonagrioides [39]. Thereafter, susceptibility to Cry1Ab has been monitored every one or two years in the most representative Spanish Bt maize production areas. The results obtained with dose-probit analysis (see Table A in S1 File) demonstrate that the LC50s of field populations and the laboratory susceptible strain have remained within a narrow range during the period 1999–2011. Resistance ratios in northeast Spain showed no clear trends, being always below 4 except for an occurrence of 7 in 2001, and were similar to those found in regions with much lower adoption rates (central and southwest Spain) (Fig 2). A resistance monitoring program, using a different batch of Cry1Ab toxin was undertaken by Monsanto Europe S.A. for field populations from northeast Spain in the period 2007–2013 [18]. Interestingly, when comparing MIC50 values of the northeast population with respect to the laboratory susceptible population over time (Fig 2), no shifts in the susceptibility of field populations were reported. Moreover, an early alert system, based on farmer complaints due to unexpected plant damage, was included in this monitoring programme. No control failures against corn borers have been yet reported from any of the Spanish Bt maize areas.

thumbnail
Fig 2. Susceptibility to Cry1Ab of Spanish field populations of Sesamia nonagrioides.

Resistance Ratios (RR) with respect to a susceptible laboratory strain were calculated from: LC50 data of field populations from northeast, central and southwest Spain in the period 1999–2011 (see Table A in S1 File); and MIC50 data reported in the resistance monitoring programme undertaken by Monsanto Europe S.A. for field populations from northeast Spain in the period 2007–2013 (in EFSA, 2015 [18]). Error bars correspond to confidence limits at 95% for RR values calculated from LC50 values.

https://doi.org/10.1371/journal.pone.0154200.g002

Factors Affecting the Rate of Resistance Evolution

Overall, the factors that could delay resistance evolution (time to an R allele frequency that exceeds 0.5) dominated those that could accelerate it. We investigated the role of poor refuge compliance by using the actual estimated use of Bt maize and refuges in northeast Spain (Fig 1 and Table H in S1 File). For two different initial allele frequencies (northeast Spain: 0.0029; EU: 0.0015, see SI for additional details), we found that even with the observed poor refuge compliance, resistance evolution would take 26–53 years (Fig 3A, BPV female movement).

thumbnail
Fig 3. Predicted rate of resistance evolution for Sesamia nonagrioides in northeast Spain.

Mean years to resistance (R allele frequency ≥ 0.5) for the interactions a) between amount of female movement and initial R allele frequency (BVP indicates the best parameter values, including the level of non-compliance to the refuge recommendation; Maximum: all females leave natal site before mating and leave mating site before ovipositing; Half Maximum: female movement parameters are halfway between BPV and Maximum); b) between local assortative mating in Bt fields and maximum female movement (initial R allele frequency = 0.0015); c) between level of external colonization and maximum female movement (legend as in 6b) (initial R allele frequency = 0.0015) d) between adoption rate of different Bt maize events and local assortative mating in Bt fields (Full: 80% adoption from year 1; Partial: observed adoption rates to 2014 and full adoption thereafter; Both 176/810: sum of adoption of Event 176 and MON 810 is 80% from year 1, and the proportion of each equals the observed proportion to 2014 and only MON 810 thereafter; initial R allele frequency = 0.0015); error bars are standard deviations of 1000 random simulations.

https://doi.org/10.1371/journal.pone.0154200.g003

Female S. nonagrioides have limited pre-mating and post-mating dispersal [25]. Because limited female dispersal is predicted to greatly reduce the rate of resistance evolution [21], we modelled the effect of increased levels of female S. nonagrioides movement on resistance evolution up to the maximum movement possible. Increased movement reduced the time to resistance, but quantitatively, the effect was small, and even at the maximum possible movement, the projected time to resistance was still 43 years, only 10 years faster than BPV (Fig 3A). At the higher initial R allele frequency, the effect of female movement virtually disappeared (Fig 3A). Limited female movement reduced the rate of resistance evolution, but was not a primary factor influencing the durability of Bt maize in northeast Spain.

In addition to the amount of refuge, the location of refuges is considered critical to delay the evolution of resistance [20, 40]. By ensuring that refuges are located near Bt fields, there can be sufficient mixing of the population so that resistant moths are likely to mate with susceptible moths, thereby delaying resistance evolution. Resistant females emerging from a Bt maize field will mate with males that come from other fields, whereas resistant males will fly to adjacent fields to mate with sedentary susceptible females [25]. If refuges are too far away, resistant females will tend to mate with resistant males, with local positive assortative mating. To evaluate the role of the spatial structure of refuge, we allowed small amounts of assortative mating only in Bt maize fields (see SI for justification). With even 1% positive assortative mating, the projected time to resistance drops by about 30% to 34 years (Fig 3B). With 4% assortative mating, this drops to only 22 years. Thus properly located refuges, whether by happenstance or plan, have been critical for maintaining the durability of Bt maize for S. nonagrioides.

In Australian cotton, Helicoverpa armigera evolved resistance to all major classes of insecticides [41], while it never occurred in the congener, H. punctigera. The main reason accounting for this difference was that there was a large influx of unselected H. punctigera each year from another geographic region followed by mating between colonists and residents, while this did not occur for H. armigera [42, 43]. It is unlikely that this was happening in northeast Spain, because S. nonagrioides is a sedentary species [17], so long-distance migration is probably uncommon. However, even if it were to occur, colonists would likely originate from the south, either southern Spain or northwestern Africa, and ride the prevailing winds to northeast Spain. Consequently, colonists would probably arrive earlier than emergence of the overwintering population in northeast Spain, because emergence at the source locations would happen much earlier. So there might be little mating between colonists and overwintered residents, and colonists would contribute only their offspring to the first generation with genetic mixing occurring in the second and third generations. In this case, the model showed that even a large influx of colonists would have no discernible effect on the rate of resistance evolution (Fig 3C).

The initial low adoption of Bt maize in northeast Spain (Fig 1) reduced selection for resistance, increasing the time to resistance. To evaluate the effect of low adoption, we modelled the evolutionary process with full adoption (80% Bt maize in northeast Spain) from the first years of commercial use under the assumption that Event 176 and MON 810 split the market at their observed shares of the Bt maize market. Event 176 would allow heterozygote survival during the 2nd and 3rd generations (Table E in S1 File), and the switch to MON 810 would maintain high-dose selection on all three generations (see SI for additional details). We also considered the possibilities that only Event 176 or only MON 810 was available and fully adopted. Low initial adoption, which allowed the timely transition from Event 176 to MON 810, was critical for the lack of detectable resistance in this EU hotspot (Fig 3D). The BPV projection indicates that resistance would evolve in 53 years. If there had been full adoption during the first years of commercialization, resistance would have occurred in only 24 years (Fig 3D), so low adoption provided 29 more years before resistance would be expected to occur under BPV conditions. If in addition, Event 176 had not been replaced by MON 810, resistance would be predicted to be widespread in 18 years (2016). Instead of accounting for the “paradox” of no observed resistance, we would be accounting for sporadic resistance failures! If only MON 810 had been available and fully adopted, resistance was projected in 70 years, significantly longer than the 53 years projected under the present conditions. Clearly, the EU policy to take Event 176 off the market was an important factor that delayed the evolution of resistance.

Conclusions

Our results suggest that if two conditions hold, 1) the initial R allele frequency was not higher than 0.0015 and mating has been and continues to be locally random, then Cry1Ab Bt maize is likely to last an additional 20+ years in northeast Spain. Although some may consider this a long time, it is important to confirm these key assumptions and to implement strategies that postpone resistance further into the future. For example, if refuges had been located so that only 2% assortative mating occurred in Bt fields, the expected time to resistance would be reduced to only 10 years from now (Fig 3B). Clearly, obtaining more accurate estimates of R allele frequency by improving rearing methods to enable routine use of F2 screens [44], and determining that mating is locally random are high priorities to improve these projections of resistance evolution.

Pyramiding multiple toxins may allow a reduction in the size of the refuge [45, 46] and/or postpone resistance further into the future if there is no cross-resistance between the toxins [47] and the refuge is large enough to sustain a viable population [48]. A number of single and stacked maize events expressing different toxins targeting corn borers are available in the USA, but only those expressing Cry1F are currently under consideration for cultivation in EU (http://www.gmo-compass.org/eng/gmo/db/). Maize hybrids expressing Cry1F toxin provide good control of European populations of S. nonagrioides [49] and Ostrinia nubilalis [50], another major maize pest in Spain. However, cross-resistance between Cry1Ab and Cry1F toxins cannot be excluded [51]. Competition experiments indicated that, in S. nonagrioides, Cry1F and Cry1Ab have different high-affinity binding sites, though Cry1Fa binds to Cry1Ab binding site with very low affinity and vice versa [52]. In addition, though low levels of cross-resistance between Cry1Ab and Cry1F have been reported for O. nubilalis [53], partial cross-recognition of the binding sites for these two toxins was found in this species [54]. Thus, substitution of single events with stacks should be carefully managed to fully profit from them [55]. Pyramided maize may also contribute to better control of other secondary lepidopteran pests, such as Mythimna unipuncta, whose potential to develop resistance to Cry1Ab-expressing maize cultivars has been recently demonstrated [56, 57]. Although reducing the refuge with pyramids might improve refuge compliance, current compliance in northeast Spain is over 90% (Table H in S1 File), and research in the US suggested that a 20% refuge was acceptable to farmers [58], so reducing the size of the refuge would probably result in only marginal improvements in compliance in Spain and would reduce the durability of the pyramid. Thus pyramiding Cry1Ab and Cry1F while maintaining the present 20% refuge requirement might provide better control of some primary and secondary pests and help to achieve long-term durability of Bt maize.

Supporting Information

S1 File. Includes 9 Tables and 7 Figures.

Study System in Spain, Formulation and Parameterization of Model, and Simulation Experiments.

https://doi.org/10.1371/journal.pone.0154200.s001

(PDF)

S2 File. Geographic coordinates for the sampling sites.

https://doi.org/10.1371/journal.pone.0154200.s002

(PDF)

Acknowledgments

We are grateful to C. Martín (Monsanto) and F. Tribó (Pioneer), and M. Eizaguirre (UdL), for providing useful information on agronomic and biological issues, respectively. We thank Syngenta Company for providing the Cry1Ab toxin.

Author Contributions

Conceived and designed the experiments: PC FO GPF DAA. Performed the experiments: PC FO GPF DAA. Analyzed the data: PC FO GPF DAA. Contributed reagents/materials/analysis tools: PC FO GPF DAA. Wrote the paper: PC FO GPF DAA.

References

  1. 1. James C. Global status of commercialized biotech/GM crops: 2014. ISAAA Brief 49. ISAAA, Ithaca, NY. 2014
  2. 2. Gómez-Barbero M, Berbel J, Rodríguez-Cerezo E. Bt corn in Spain-the performance of the EU´s first GM crop. Nat Biotechnol. 2008; 26: 384–386. pmid:18392015
  3. 3. EFSA. Scientific Opinion on the annual Post-Market Environmental Monitoring (PMEM) report from Monsanto Europe S.A. on the cultivation of genetically modified maize MON 810 in 2010. EFSA Journal 2012; 10(4): 2610.
  4. 4. Tabashnik BE, Brevault T, Carriere Y. Insect resistance to Bt crops: lessons from the first billion acres. Nat Biotechnol. 2013; 31: 510–521. pmid:23752438
  5. 5. Van den Berg J, Hilbeck A, Bøhn T. Pest resistance to Cry1Ab Bt maize: Field resistance, contributing factors and lessons from South Africa. Crop Prot. 2013; 54: 154–160.
  6. 6. Van Rensburg JBJ. First report of field resistance by the stem borer, Busseola fusca (Fuller) to Bt-transgenic maize. S Afr J Plant Soil. 2007; 24: 147–151.
  7. 7. Storer NP, Babcock JM, Schlenz M, Meade T, Thompson GD, Bing JW, et al. Discovery and characterization of field resistance to Bt Maize: Spodoptera frugiperda (Lepidoptera: Noctuidae) in Puerto Rico. J Econ Entomol. 2010; 103: 1031–1038. pmid:20857709
  8. 8. Farias J, Andow DA, Horikoshi RJ, Sorgatto PF, Dos Santos AC, Omoto C. Field-evolved resistance to Cry1F maize by Spodoptera frugiperda (Lepidoptera: Noctuidae) in Brazil. Crop Prot. 2014, 64: 150–158.
  9. 9. Omoto C, Bernardi O, Salmeron E, Sorgatto RJ, Dourado PM, Crivellari A, Carvalho RA, Willse A, Martinelli S, Head GP. Field-evolved resistance to Cry1Ab maize by Spodoptera frugiperda in Brazil. Pest Manag Sci. 2016. on-line first,
  10. 10. Van Rensburg JBJ, Walters MC, Giliomee JH. Geographical variation in the seasonal moth flight activity of the maize stalk borer, Busseola fusca (Fuller), in South Africa. S Afr J Plant Soil. 1985; 2: 123–126.
  11. 11. Farias J, Andow DA, Horikoshi RJ, Sorgatto PF, Dos Santos AC, Omoto C. Dominance of a Cry1F resistance in Spodoptera frugiperda (Lepidoptera: Noctuidae) on TC1507 Bt maize in Brazil. Pest Manag Sci. 2015; on-line.
  12. 12. Tellez-Rodriguez P, Raymond B, Morán-Bertot I, Rodríguez-Cabrera L, Wright DJ, Borroto CG, et al. Strong oviposition preference for Bt over non-Bt maize in Spodoptera frugiperda and its implications for the evolution of resistance. BMC Biology. 2014; 12: 48. pmid:24935031
  13. 13. Székács A, Lauber É, Juracsek J, Darvas B. Cry1Ab toxin production of MON 810 transgenic maize. Environ Toxicol Chem. 2010; 29: 182–190. pmid:20821434
  14. 14. Fearing PL, Brown D, Vlachos D, Meghji M, Privalle L. Quantitative analysis of CryIA(b) expression in Bt maize plants, tissues, and silage and stability of expression over successive generations. Mol Breed. 1997; 3: 169–176.
  15. 15. Eizaguirre M, Tort S, López C, Albajes R. Effects of sublethal concentrations of Bacillus thuringiensis on larval development of Sesamia nonagrioides. J Econ Entomol. 2005; 98: 464–470. pmid:15889739
  16. 16. Farinós GP, de la Poza M, Hernández-Crespo P, Ortego F, Castañera P. Resistance monitoring of field populations of the corn borers Sesamia nonagrioides and Ostrinia nubilalis after 5 years of Bt maize cultivation in Spain. Entomol Exp Appl. 2004; 110: 23–30.
  17. 17. De la Poza M, Farinós GP, Beroiz B, Ortego F, Hernández-Crespo P, Castañera P. Bt maize resistance management: genetic structure of Sesamia nonagrioides (Lefebvre) populations in the Mediterranean area. Environ Entomol. 2008; 37: 1354–1360.
  18. 18. EFSA. Scientific Opinion on the annual post-market environmental monitoring (PMEM) report from Monsanto Europe S.A. on the cultivation of genetically modified maize MON 810 in 2013. EFSA Journal 2015;13(3):4039.
  19. 19. Robertson JL, Russell RM, Preisler HK, Savin NE. Pesticide Bioassays with arthropods. CRC, Boca Raton, FL. 2007.
  20. 20. Comins HN. The development of insecticide resistance in the presence of migration. J Theor Biol. 1977; 64: 177–197. pmid:556789
  21. 21. Ives AR, Andow DA. Evolution of resistance to Bt crops: Directional selection in structured environments. Ecol Lett. 2002; 5: 792–801.
  22. 22. Alstad DN, Andow DA. Managing the evolution of insect resistance to transgenic plants. Science. 1995; 268: 1894–1896. pmid:17797533
  23. 23. Andow DA, Ives AR. Monitoring and adaptive resistance management. Ecol Appl. 2002; 12: 1378–1390.
  24. 24. Andreadis SS, Álvarez-Alfageme F, Sánchez-Ramos I, Stodola TJ, Andow DA, Milonas PG, et al. Frequency of resistance to Bacillus thuringiensis toxin Cry1Ab in Greek and Spanish population of Sesamia nonagrioides (Lepidoptera: Noctuidae). J Econ Entomol. 2007; 100: 195–201. pmid:17370828
  25. 25. Eizaguirre M, López C, Albajes R. Dispersal capacity in the Mediterranean corn borer, Sesamia nonagrioides (Lepidoptera: Noctuidae). Entomol Exp Appl. 2004; 113: 25–34.
  26. 26. López C, Sans A, Eizaguirre M. Influencia de la planta de maíz en el apareamiento de Sesamia nonagrioides Lefèbvre (Lepidoptera: Noctuidae). Inv Agr: Prod Prot Veg. 1999; 14: 415–422.
  27. 27. Larue P. La Sésamie du maïs (Sesamia nonagrioides Lef.). Dégâts et actualisation de lutte. Défense Végétaux. 1984; 227: 163–181.
  28. 28. López C, Eizaguirre M, Albajes R. Courtship and mating behaviour of the Mediterranean corn borer, Sesamia nonagrioides (Lepidoptera: Noctuidae). Spanish J Agric Res. 2003; 1: 43–51.
  29. 29. Albajes R, Eras J, López C, Ferran X, Vigata J, Eizaguirre M. Testing rubidium marking for measuring adult dispersal of the corn borer Sesamia nonagrioides: first results. IOBC/WPRS Bull. 2004; 27 (3): 15–22.
  30. 30. Eizaguirre M, Albajes R, López C, Eras J, Lumbierres B, Pons X. Six years after the commercial introduction of Bt maize in Spain: field evaluation, impact and future prospects. Transg Res. 2006; 15: 1–12.
  31. 31. Fantinou AA, Perdikis DC, Zota KF. Reproductive responses to photoperiod and temperature by diapausing and nondiapausing populations of Sesamia nonagrioides Lef. (Lepidoptera–Noctuidae). Physiol Entomol. 2004; 29: 169–175.
  32. 32. Eizaguirre M, Albajes R. Diapause induction in the stem corn borer, Sesamia nonagrioides (Lepidoptera: Noctuidae). Entomol Gen. 1992; 17: 277–283.
  33. 33. Gillyboeuf N, Anglade R, Lavenseau L, Peypelut L. Cold hardiness and overwintering strategy of the pink maize stalk borer, Sesamia nonagrioides Lef. (Lepidoptera, Noctuidae). Oecologia. 1994; 99: 366–373.
  34. 34. Dutton A, Klein H, Romeis J, Bigler F. Uptake of Bt-toxin by herbivores feeding on transgenic maize and consequences for the predator Chrysoperla carnea. Ecol Entomol. 2002; 27: 441–447.
  35. 35. Obrist LB, Dutton A, Albajes R, Bigler F. Exposure of arthropod predators to Cry1Ab toxin in Bt maize fields. Ecol Entomol. 2006; 31: 143–154.
  36. 36. Zoerb AC, Spencer T, Hellmich RL, Wright RJ, Siegfried DB. Larval distribution and survival of second generation European corn borer, Ostrinia nubilalis (Hübner) (Lepidoptera: Crambidae) on Event 176 Bt Corn. Crop Prot. 2003; 22: 179–184.
  37. 37. Eizaguirre M, López C, Albajes R. Sublethal effects of Bacillus thuringiensis on larval development in Sesamia nonagrioides J Econ Entomol. 2005; 98: 237–247.
  38. 38. Gasmann AJ, Carrière Y, Tabashnik BE. Fitness costs of insect resistance to Bacillus thuringiensis. Annu Rev Entomol. 2009; 54: 147–163. pmid:19067630
  39. 39. González-Núñez M, Ortego F, Castañera P. Susceptibility of Spanish populations of the corn borers Sesamia nonagrioides (Lepidoptera: Noctuidae) and Ostrinia nubilalis (Lepidoptera: Crambidae) to a Bacillus thuringiensis endotoxin. J Econ Entomol. 2000; 93: 459–463. pmid:10826200
  40. 40. Tabashnik BE, Croft BA. Managing pesticide resistance in crop-arthropod complexes: Interactions between biological and operational factors. Environ Entomol. 1982; 11: 1137–1144.
  41. 41. Forrester NW, Cahill M, Bird LJ, Layland JK. Management of pyrethroid and endosulfan resistance in Helicoverpa armigera (Hübner) in Australia. Bull Entomol Res. 1993; Supplement Series No. 1, 132 pp.
  42. 42. Fitt GP, Daly JC. Abundance of overwintering pupae and the spring generation of Helicoverpa spp. (Lepidoptera Noctuidae) in New South Wales, Australia: consequences for pest management. J Econ Entomol. 1990; 83: 1827–1836.
  43. 43. Gregg PC, Fitt GP, Zalucki MP, Murray DAH. Insect migration in an arid continent II. Helicoverpa spp. in Australia. In Drake VA, Gatehouse AG, editors. Insect Migration: Physical Factors and Physiological Mechanisms. Cambridge University Press, Cambridge, UK. 1995. pp. 151–172.
  44. 44. Andow DA, Alstad DN. The F2 screen for rare resistance alleles. J Econ Entomol. 1998. 91: 572–578.
  45. 45. Roush RT. Two-toxin strategies for management of insecticidal transgenic crops: Can pyramiding succeed where pesticide mixtures have not? Phil Trans R Soc London. 1998; 353: 1777–1786.
  46. 46. Alyokhin A. Scant evidence supports EPA’s pyramided Bt corn refuge size of 5%. Nat Biotechnol. 2011; 29: 577–578. pmid:21747379
  47. 47. Curtis CF. Theoretical models of the use of insecticide mixtures for the management of resistance. Bull Entomol Res. 1985; 75: 259–65.
  48. 48. Ives AR, Glaum PR, Ziebarth NL, Andow DA. The evolution of resistance to two-toxin pyramid transgenic crops. Ecol Appl. 2011; 21: 503–515. pmid:21563580
  49. 49. Farinós GP, de la Poza M, Ortego F, Castañera P. Susceptibility to the Cry1F toxin of field populations of Sesamia nonagrioides (Lepidoptera: Noctuidae) in Mediterranean maize cultivation regions. J Econ Entomol. 2012; 105: 214–221. pmid:22420274
  50. 50. Gaspers C, Siegfried BD, Spencer T, Alves AP, Storer NP, Schuphan I, Eber S. Susceptibility of European and North American populations of the European corn borer to the Cry1F insecticidal protein. J Appl Entomol 2011; 135: 7–16.
  51. 51. Jakka SRK, Gong L, Hasler J, Banerjee R, Sheets JJ, Narva K, Blanco CA, Jurat-Fuentes JL. Field-evolved Mode 1 fall armyworm resistance to Bt corn associated with reduced Cry1Fa toxin binding and midgut alkaline phosphatase expression. Appl Environ Microbiol 2016.
  52. 52. González-Cabrera J, Farinós GP, Caccia S, Díaz-Mendoza M, Castañera P, Leonardi MG, et al. Toxicity and mode of action of Bacillus thuringiensis Cry proteins in the Mediterranean corn borer, Sesamia nonagrioides (Lefebvre). Appl Environ Microbiol. 2006; 72: 2594–2600. pmid:16597962
  53. 53. Crespo LB, Rodrigo-Simón A, Siqueira HAA, Pereira EJG, Ferré J, Siegfried BD. Cross-resistance and mechanism of resistance to Cry1Ab toxin from Bacillus thuringiensis in a field-derived strain of European corn borer, Ostrinia nubilalis. J Invertebrate Pathol. 2011; 107: 185–192.
  54. 54. Hua G, Masson L, Jurat-Fuentes JL, Schwab G, Adang MJ. Binding analyses of Bacillus thuringiensis Cry δ-endotoxins using brush border membrane vesicles of Ostrinia nubilalis. Appl Environ Microbiol. 2001; 67: 872–879. pmid:11157257
  55. 55. Carrière Y, Crickmore N, Tabashnik BE. Optimizing pyramided transgenic Bt crops for sustainable pest management. Nature Biotech 2015; 33: 161–168.
  56. 56. Gonzalez-Cabrera J, Garcia M, Hernandez-Crespo P, Farinos GP, Ortego F, Castañera P. Resistance to Bt maize in Mythimna unipuncta (Lepidoptera: Noctuidae) is mediated by alteration in Cry1Ab protein activation. Insect Biochem Mol Biol. 2013, 43: 635–643. pmid:23603093
  57. 57. García M, Ortego F, Hernández-Crespo P, Farinós GP, Castañera P. Inheritance, fitness costs, incomplete resistance and feeding preferences in a laboratory-selected MON810-resistant strain of the true armyworm Mythimna unipuncta. Pest Manag Sci. 2015; 71:1631–1639. pmid:25565147
  58. 58. Hurley TH, Babcock BA, Hellmich RA. Bt Corn and insect resistance: An economic assessment of refuges. J Agric Res Economics. 2001; 26: 176–194.