Skip to main content

Strong oviposition preference for Bt over non-Bt maize in Spodoptera frugiperda and its implications for the evolution of resistance



Transgenic crops expressing Bt toxins have substantial benefits for growers in terms of reduced synthetic insecticide inputs, area-wide pest management and yield. This valuable technology depends upon delaying the evolution of resistance. The ‘high dose/refuge strategy’, in which a refuge of non-Bt plants is planted in close proximity to the Bt crop, is the foundation of most existing resistance management. Most theoretical analyses of the high dose/refuge strategy assume random oviposition across refugia and Bt crops.


In this study we examined oviposition and survival of Spodoptera frugiperda across conventional and Bt maize and explored the impact of oviposition behavior on the evolution of resistance in simulation models. Over six growing seasons oviposition rates per plant were higher in Bt crops than in refugia. The Cry1F Bt maize variety retained largely undamaged leaves, and oviposition preference was correlated with the level of feeding damage in the refuge. In simulation models, damage-avoiding oviposition accelerated the evolution of resistance and either led to requirements for larger refugia or undermined resistance management altogether. Since larval densities affected oviposition preferences, pest population dynamics affected resistance evolution: larger refugia were weakly beneficial for resistance management if they increased pest population sizes and the concomitant degree of leaf damage.


Damaged host plants have reduced attractiveness to many insect pests, and crops expressing Bt toxins are generally less damaged than conventional counterparts. Resistance management strategies should take account of this behavior, as it has the potential to undermine the effectiveness of existing practice, especially in the tropics where many pests are polyvoltinous. Efforts to bring down total pest population sizes and/or increase the attractiveness of damaged conventional plants will have substantial benefits for slowing the evolution of resistance.


Transgenic plants expressing the Bacillus thuringiensis (Bt) Cry toxins (Bt crops), first commercialized in 1996, represent a technological revolution in agricultural pest control, and can confer improved pest management, reduced pesticide applications, increased yields and area-wide reductions in pest populations [13]. By 2012, Bt crops were planted on about 70 million hectares worldwide [4]. Currently, the market of Bt crops is dominated by maize and cotton cultivars engineered to target Lepidoptera and/or Coleopteran pests. The primary threat to the ongoing success of Bt crops is the evolution of resistance by insect pests [5] and resistance management strategies are a legal requirement of the licensing of transgenic crops in many countries. Bt resistance is typically controlled by single large effect genes, with resistance being recessive at high doses of Cry toxin [6]. Regulatory agencies commonly demand that Bt cultivars express Cry toxins at a dose that will impose near completely recessive resistance, that is, be high enough to kill nearly all heterozygous resistant insects. In conjunction with these high doses, farmers are required to plant a refuge area of a non-Bt variety of the same crop in close proximity to the Bt crop [7]. Under this ‘high dose/refuge strategy’, the more abundant Bt susceptible adults emerging from the refuge will tend to mate with any resistant homozygous moths emerging from the Bt crop field [7, 8]. The resulting heterozygous progeny should still be susceptible to the high dose of Bt toxins. So far, this strategy has successfully delayed resistance for Bt crops in many key target species [5, 7].

However, field-evolved resistance to crops expressing Bt toxins has occurred at least in five cases where both failed insect control and an increased frequency of resistance genes have been recorded [9, 10]. Notably, one of first documented cases of insect-resistance in Bt crops was in fall armyworm Spodoptera frugiperda (J.E. Smith) (Lepidoptera: Noctuidae) in Puerto Rico [11]. Resistance enabled S. frugiperda to survive and feed on Bt maize expressing Cry1F toxin in the field and led to the withdrawal of this product from Puerto Rico [11]. There are believed to be two main causes for resistance problems in general. First, growers may not strictly adhere to the requirements for planting refugia [12, 13]. Second, toxin doses have been too low or variable to consistently kill heterozygous resistant insects [5, 11, 12].

Several oviposition preference studies have found adult females cannot discriminate between transgenic and near isogenic conventional cultivars [1416]. However, these experiments only investigated the number of eggs on undamaged plants. This methodology, therefore, ignores any potential effects of subsequent larval feeding. Since many Bt crops kill newly hatched larvae soon after emergence, the feeding damage caused by larvae is likely to differ between conventional and Bt crops as the growing season progresses. This is significant as the volatiles released during larval feeding can deter oviposition in adult female Lepidoptera, a behavior that can reduce intra-specific competition or reduce detection by natural enemies that use plant volatiles as foraging cues [17, 18].

In the present study, we tested the hypothesis that oviposition rates of S. frugiperda would respond to differences in damage levels incurred by Bt and non-Bt maize under field conditions. The field study covered both dry and wet growing seasons over four years in the Caribbean island of Cuba, which has similar growing conditions to neighboring Puerto Rico, where resistance to Cry1F maize was first documented [11]. S. frugiperda is the primary pest of maize (Zea mays) in the tropical region of the western hemisphere, where frequent insecticide applications are needed to maintain the population below economic thresholds [19]. We recorded oviposition rates, crop damage and larval numbers of S. frugiperda in Bt maize and its conventional refuge in young maize plants. We used these data to estimate in-field survivorship from egg to third instar and to inform simulation models that explored the impact of non-random damage-avoiding oviposition on the evolution of resistance under a high dose/refuge management strategy.


Oviposition preference for Bt-maize in S. frugiperdacorrelates with the level of feeding damage in the refuge

Damage to the vegetative crop was assessed on a four-point visual scale (Figure 1A). Data from the first observation each season (20 to 25 days post-planting) showed that transgenic FR-Bt1 plants were almost entirely undamaged, while the conventional variety FR-28 suffered a wider range of damage levels (Figure 1B). Statistical comparisons of damage scores were made twice per season using Kruskal Wallis non-parametric tests. Every comparison showed highly significant differences [see Additional file 1: Table S1], although damage scores were slightly reduced at the later observation period.

Figure 1
figure 1

Feeding damage and oviposition preference of S. frugiperda in transgenic Bt plants and conventional refuge plants. (A) The damage scores of representative maize plants. (B) The variation in damage scores between transgenic Bt plants and conventional refugia plants. Box and whisker plots display medians (solid circles), lower and upper quartiles (boxes), the range (whiskers) and outliers (open circles). (C) Fall armyworm lay fewer egg masses on Bt plants (red bars) than on conventional plants (blue), than would be predicted by chance (Pearson’s χ2 = 92.9, df = 5, P < <0.0001), errors bars are ± SE. (D) The strength of preference for laying eggs in Bt maize increases with crop damage in the refuge within the growing season (overall slope = 1.4, SE = 0.38, t = 3.67, P = 0.0037). Data are from the wet growing seasons in 2009, 2011 and 2012. After taking the square roots of egg mass counts, preference is defined as egg masses per plant in Bt crop - egg masses per plant in refuge, while the damage differential is the difference in mean damage scores.

Scouting of egg masses allowed us to assess oviposition rates in Bt and the conventional parts of the crop. The number of egg masses per plant was significantly higher in Bt compared to non-Bt fields over six growing seasons (Figure 1C). Bias coefficients, which reflect the proportional reduction in egg counts in refugia relative to random expectation (see Methods), varied from 0.16 to 0.59, with a mean of 0.47. In other words, over the course of this study we recorded 145 egg masses in refugia, only 47% of the total of 309 masses predicted from a hypothesis of random oviposition. Crop damage levels varied within season, so we were also able to assess whether oviposition patterns responded to changing levels of damage in refugia. We analyzed data from the three seasons with the highest abundance of S. frugiperda and in which experimental plots were planted with 50% conventional refugia. The preference for Bt maize was significantly positively correlated with levels of damage in conventional refugia (Figure 1D, P = 0.0036, df = 1, likelihood ratio 8.49). This analysis also indicates that there were no intrinsic differences in preference for Bt maize over conventional maize, since when both crop types were undamaged oviposition rates were similar. This was statistically supported by our fitted model, in which the intercept was not significantly different from zero (value = 0.47, SE = 0.63, t = 0.75, P = 0.47) indicating random oviposition on undamaged conventional and Bt maize.

We inferred survivorship of S. frugiperda larvae from scouting observations and from the assumption that each egg mass contains, on average, 150 fertile eggs [20]. Survival on Bt maize ranged from 0.006 to 0.0003 [see Additional file 1: Table S2]. Moreover, comparisons between years showed that natural mortality in the conventional crop was strongly negatively density-dependent, with mortality following an inverse power law [see Additional file 2: Figure S1].

Biased oviposition has substantial effects on the evolution of resistance

We investigated the effects of damage-avoiding oviposition on the evolution of resistance using computer simulation models. These models were informed by our field data and by previous reports of the genetics of Bt resistance in S. frugiperda[11, 21]. Larval population densities in conventional refugia determined oviposition bias in these models. We considered how different pest management regimes with different forms of density dependence would affect the evolution of resistance, as these will impact larval densities, leaf damage and damage-avoiding oviposition behavior. We considered three scenarios: a spray threshold model in which application of synthetic pesticides in refugia was based on the number of egg masses laid per crop plant; a natural predation model with density-dependent mortality of larvae in refugia only (based upon the power law inferred from field data); and a predation model in which there was density-dependent mortality occurring throughout the Bt crop and the refuge in addition to the mortality imposed by the plant incorporated Bt toxin.

When oviposition was random, we generally recovered the typical result in which increasing refuge size provided better control of resistance (Figure 2). The one exception to this rule was that the evolution of fully recessive resistance was prevented at very low refuge sizes under the spray threshold scenario (Figure 2A). This occurred because simulated pests formed ‘population sinks’ with negative population growth at very low refuge sizes so that population sizes were never large enough (<10,000 adults) to produce any resistant homozygotes. Refugia of 10% to 30% of the total crop area were typically required to control the evolution of resistance when initial resistance was fully recessive (Figure 2A,C,E), although stable management of resistance was not possible under partially dominant resistance (Figure 2B,D,F). Damaging-avoiding oviposition accelerated resistance evolution, and undermined resistance management. Only under the most favorable conditions (recessive resistance, refugia >30% and density dependence throughout the crop) was stable resistance management possible with damage-avoiding oviposition (Figure 2E). When there was density dependence throughout the crop refuge size and the dominance of resistance appeared to interact: refuge size was a more important parameter as dominance decreased [see Additional file 3: Figure S2].Moreover, the problems caused by damage-avoiding oviposition could not typically be resolved by planting larger refugia, especially when resistance was partially dominant. This insensitivity to refuge size occurred because increasing refuge sizes increased adult population density and overall pest pressure. Higher population densities led to increased damage rates and more biased oviposition. With fully recessive resistance increasing refugia size had some benefits, especially when we imposed some density-dependent mortality throughout the crop (Figure 2C,E). The natural predation scenario could provide better resistance management at higher refuge size because strong density-dependent mortality ensured that oviposition bias did not increase with population size (Figure 2G). Conversely, under a spray threshold regime, the proportion of eggs laid in refugia decreased dramatically with population size, thereby offsetting the benefits of planting a larger refuge (Figure 2H). Note that under both scenarios egg counts in Bt plants were at most three times greater than expected by chance, a fairly conservative maximum given our field data.

Figure 2
figure 2

Evolution of resistance to a single Bt toxin crop under random oviposition (magenta lines) and damage-avoiding oviposition (blue lines). Models simulated three scenarios for pest management with different forms of density dependence under fully recessive (h = 0) or partially dominant resistance (h = 0.1). In A and B pest management in refugia uses a density-dependent spray threshold of 0.25 egg masses per plant; in C and D natural enemies cause density-dependent mortality in refugia while in E and F natural enemies cause density-dependent mortality throughout the refugia and the Bt crop. The response variable is the number of generations required for the resistance gene frequency (R) to equal or exceed 0.5, up to maximum of 200. Note that a value of 200 is returned if the frequency of resistance never exceeded 0.5. Fecundity was set at 300 eggs per female and the initial resistance gene frequency was 0.001. G and H show the relationship between population density and avoidance of refugia under the two modeling scenarios with h = 0.1. Parameter values in the last two panels are calculated before the evolution of full resistance, that is, from generation 1 until the time point at which the frequency R > = 0.5. The refuge avoidance factor is 1/b and represents the multiplicative reduction in eggs laid in the refuge relative to random expectation, so that 1 indicates random oviposition and 3 indicates a three-fold reduction in egg counts relative to random expectation.

Pest population dynamics strongly influence resistance evolution

We further explored the impact of population dynamics on the evolution of resistance in a series of sensitivity analyses. With random oviposition under the spray-threshold scenario, female fecundity affected evolution of resistance when it led to negative population growth under partially dominant (Figure 3) and fully recessive resistance [see Additional file 4: Figure S3]. However, under damage-avoiding oviposition, increased fecundity reduced the time to evolve resistance over a range of parameter values. Simulation results were, however, relatively insensitive to spray thresholds, so that decisions to minimize damage by spraying lightly infested refugia are not predicted to improve resistance management [see Additional file 5: Figure S4]. Spray thresholds were only important with large refugia under damage-avoiding oviposition; under those conditions higher spray thresholds (>0.6 egg masses per plant) tended to delay resistance [see Additional file 5: Figure S4].

Figure 3
figure 3

The effects of random and damage-avoiding oviposition behavior and female fecundity on evolution of resistance under the sprayed refuge scenario. In these simulations the spray action threshold for refugia was 0.25 egg masses per plant, the initial frequency of resistance = 0.001 and we assumed partially dominant resistance (h = 0.1).


Here, we have shown that S. frugiperda has a strong oviposition preference for undamaged Bt maize over damaged conventional maize planted as a refuge for resistance management. Oviposition by S. frugiperda affects the production of plant volatiles in maize [22]. Moreover, ovipositing females show a preference for the volatiles of undamaged maize plants [23]. Early reports have also described how the egg-laying behavior of this species changes when larval densities are high [24]. We controlled for the genetic background of the Bt and conventional crop species, so that differences in maize variety other than the presence or absence of Cry1F are unlikely to have affected insect behavior. Since oviposition preference within season also increased with the level of leaf damage in the conventional refuge the most plausible explanation for increased egg laying in the Bt crop field is that females are avoiding plants previously damaged by conspecifics. This behavior makes good adaptive sense for ovipositing females, as larval progeny on damaged plants may suffer increased competition for food and more attention from foraging parasitoids [25, 26]. One caveat of our observations is that it is uncertain whether differences in oviposition rates will also hold true when there are much greater distances between refugia and Bt crops. The ‘in-field’ refugia that were used in this study are relatively rare in some maize cropping areas [27]. That being said, all refugia must be planted close enough to Bt crops in order to facilitate random mating between adults emerging from transgenic and conventional plants, so that females on the wing should be able to readily sample volatiles from both crop types.

Modeling of the evolution of resistance to Bt under damage-avoiding oviposition has previously assumed that adults will prefer to lay eggs on conventional rather than Bt crops [28], perhaps because Bt and its toxins have been shown to deter larval feeding [29, 30]. This earlier theoretical work is consistent with our results, in that resistance management improves as more eggs are laid in refugia [28]. However, this earlier hypothesis of avoidance of Bt crops is unsupported [1416]. Damage on Bt crops is typically less than on their non-Bt counterparts. Given the many reports of avoidance of damaged plants or damage-induced volatiles by adult Lepidoptera [17, 18, 26, 3134], the opposite prediction is more likely to hold: reduced leaf damage in Bt crops will attract increased rates of oviposition. Since leaf damage is determined by larval density, our modeling showed that the population dynamics of pest species in conventional refugia have a strong impact on the efficacy of the high dose/refuge management strategy.

This conclusion is supported by the details of the evolution of resistance to Bt maize by S. frugiperda in Puerto Rico [11]. In that case, population densities of S. frugiperda were high prior to reports of control failure [11]. This population outbreak was compounded by a drought, which eliminated natural refugia in wild grasses, concentrating larval populations onto irrigated maize crops. This would have led to very high damage levels on conventional maize plants and a substantial weakening of the efficacy of any planted refugia. In Puerto Rico, the evolution of resistance occurred within 30 to 40 generations [11]. Simulation models that incorporate damage avoidance predict the evolution of resistance within a similar period of time (27 generations), when we make assumptions that fit with the Puerto Rican experience, that is, partial dominance and a sprayed 20% to 30% refuge (N. Storer, personal communication). In contrast, if we assume random oviposition, the predicted time range for the evolution of resistance is 39 to 65 generations, a range which is largely outside the Puerto Rican estimate. This correspondence between field data and modeling is reassuring and suggests damage avoidance could have accelerated the evolution of resistance to Cry1F in Puerto Rico. However, this match must be treated with caution as it was not our intention to explicitly simulate the Puerto Rican situation, and important parameters, such as the initial frequencies of resistance in Puerto Rico, were unknown.

Our results have several implications for resistance management in Bt crops. One clear recommendation is that there are substantial gains to be made by improving area-wide pest population control, as lower population sizes will decrease the difference in damage levels between Bt crops and refugia. A second implication is that resistance management strategies that are less dependent on refugia will be less vulnerable to non-random oviposition. One such possible strategy is the large-scale release of transgenic female-specific RIDL (release of insects with dominant lethality) insects for population control, as mating with large numbers of Bt susceptible released males is the mechanism that lowers the frequency of resistance in the field [35]. Insect population dynamics and density dependence were important for resistance evolution, but had the most impact when we assumed recessive resistance. Large, unsprayed conventional refugia with control of pests by natural enemies could provide good resistance management under damage-avoiding oviposition, if we assumed there was also some density-dependent mortality in the Bt portion of the crop. One caveat here is that our endemic population of S. frugiperda was relatively well controlled by density-dependent larval mortality. This natural regulation can be disrupted by the over-use of broad-spectrum insecticides or by the migration of S. frugiperda, since migrating populations can escape natural enemy control [24]. Nevertheless, when natural enemy regulation is poor, high spraying rates in refugia did not benefit resistance management in our simulations.

The presented simulation results depend upon the naïve assumption that pest population dynamics are entirely driven by population processes in crops, an assumption that will be true to a greater or lesser extent in different pest species and at different times. For example, while S. frugiperda can feed on a range of wild grass species in the wet season [24], during the dry season this species is more dependent on irrigated crops [11]. Pests with a narrow host range are clearly more dependent on agro-ecosystems; however, agricultural practice can have substantial impacts on the populations of polyphagous pests [1, 2]. In the event that agricultural practice has no impact on pest population dynamics, it would still be beneficial to plant the largest non-Bt refuge possible, as larger refugia would not affect oviposition bias.

While we found that damage-avoiding oviposition has the potential to undermine current resistance management strategies, there are several possibilities that might mitigate this effect. Maize plants, for example, show considerable genetic variation in the extent and make-up of volatile release after damage [36]. It may, therefore, be possible to plant non-Bt refugia with crop plants that remain attractive to ovipositing females after being damaged by larvae. Refugia are more effective when plants are less damaged and insect populations reduced. Efforts to bring down regional population sizes in an integrated pest management program (IPM) be they via pheromone trapping, habitat management (control of alternative hosts) or mechanical ‘pupae busting’ , will all yield benefits for diverse insect pests in the face of damage averse oviposition [37]. We expected that lowering the density of insects in refugia by imposing tighter spray thresholds might also mitigate the impact of damage-avoiding oviposition. However, simulation results were largely insensitive to refugia spray thresholds, possibly because increasing spray rates also decrease insect population sizes in refugia [see Additional file 3: Figure S2] or because density-dependent mortality in refugia can itself accelerate the evolution of resistance [38, 39].


Modern multiple toxin Bt crop varieties are expected to improve resistance management [40, 41]. Nevertheless, these crops are not expected to be immune to the evolution of resistance [10, 42] and non-Bt refugia are still required for effective resistance management [40, 41]. In both single and multiple toxins systems the efficacy of any planted refuge will be reduced by oviposition bias, especially where pests are polyvoltinous or have overlapping generations and are, therefore, capable of responding to the damage caused by larvae feeding earlier in the season. A number of insect species are deterred from laying eggs on plants damaged by conspecific larvae [17, 18, 26, 3134]. These species include prominent target pests for a wide range of Bt crops, such as Heliothis virescens, S. frugiperda, Spodoptera littoralis and Ostrinia nubilalis[17, 18, 23, 34]. For H. virescens and O. nubilalis current resistance management practices for slowing the evolution of resistance to Bt crops appear to be effective [1, 43], despite possible issues with damage-avoiding oviposition. This could be because area-wide management of these pests has been effective enough to prevent high levels of damage in refugia [1, 2], because the results here which are based on a small-scale field trial may not scale up to larger refugia or because bivoltine dynamics limit the opportunities for damage-avoiding oviposition to a single generation per year [44]. Despite these caveats it is likely that the impact of oviposition behavior on the evolution of resistance to Bt crops has been underestimated and needs careful re-evaluation.



The Bt (Cry1F) maize variety FR-Bt1 (referred to as Bt maize) along with the Cuban open-pollinated variety FR-28 (referred to as non-Bt maize) [45] were used in all the experiments. FR-Bt1 was produced in the CIGB by crossing plants of a commercial Bt maize bearing the TC1507 event into FR-28 plants and back-crossing four additional times into this conventional maize variety. Before each backcross, FR-Bt1 was selected for Bt Cry1F expression by ELISA immunoassay (Agdia Inc., Elkhart, IN, USA). Based on the number of backcrosses, the resultant FR-Bt1 plants contain 96.75% of the genetic material of the FR-28 genotype.

Study area and experimental conditions

The experiments were conducted at the CIGB experimental field station from 2009 to 2012. Bt and non-Bt maize cultivars were planted during the second week of January (for the dry season) and first week of May (for the wet season). Plants were arranged in a randomized complete block design with four replicates. Each field plot had three rows, and row size ranged from 4 to 10 m long. The distance between plots of Bt corn and non-Bt corn was five meters. All studies followed standard agronomic practices for maize. Seeds used for testing were treated only with a commercial fungicide (Poncho 250, Bayer CropSciences, Monheim, Germany) and no insecticides were applied to any part of the plots. Natural rainfall was supplemented by irrigating weekly with 7 cm of water as needed, particularly during the dry season.

The number of egg masses, the rank of leaf-feeding injury and the larval development stage were scored from the middle row of the four independent plots repeatedly over the growing season. Observations began at the V3 (three true-leaves) stage and were taken weekly until no new egg masses were found for two consecutive weeks. To prevent duplicate counting, each new egg mass was circled on the leaf with a waterproof marker. The total plants sampled in the Bt crop and conventional refuge respectively were: wet season 2009, N = 376 and N = 236; dry season 2009, N = 175 and N = 170; wet season 2010, N = 670 and N = 80; dry season 2010, N = 530 and N = 53; for all remaining seasons N = 100 in both sections of the crop. The damage from S. frugiperda herbivory was assessed at around 20 days after planting and on the last day any significant damage was detected. Leaf damage was estimated using a four-point damage visual-rating scale (1 to 4), where a score of 1 indicates no visible damage, while a score of 4 represents >90% destruction of whorl and furl leaves (Figure 1A). For the estimates of larval survival, a mean of 150 eggs per egg mass was assumed. Larval scouting began at the V3 stage and was maintained until one week before initiation of flowering. Instars were determined by measuring larval size and head capsule width [46]. The larvae were considered dead if they did not move when touched with a needle.

Data analysis

Data on oviposition rates were analyzed using Pearson’s Chi-squared test, using mean egg masses per plant in the whole crop to calculate expected frequencies. All statistical tests were carried out in R v 2.6.2 [47]. Mixed models of oviposition preference within season used mean damage in the refuge as a fixed factor and year of planting as a random effect. Models with different fixed effect structures were compared using likelihood ratio tests following maximum likelihood model fitting. Final model fitting and model assumptions were checked with graphical analyses [48].

Simulation models

The deterministic two-patch population models followed a published framework [49] and were carried out in R v 2.6.2. A single gene determined resistance. Storer et al. (2010) estimated the dominance of resistance in terms of survival on Cry1F plants as between 0.07 and 0.14 [11]. All models considered at least two values for dominance: completely recessive (h = 0) and partially dominant (h = 0.1). Further sensitivity analyses on the effects of dominance were conducted as appropriate [see Additional file 3: Figure S2]. Fitness costs of resistance were recessive and fixed at a 25% reduction in larval survival. We assumed random mating and used data from previous reports of CryIF resistance [11, 21] and our field observations to determine life history parameters [see Additional file 1: Table S3]. Pest populations in the model were closed to immigration and emigration. We assumed non-overlapping generations and imposed a minimum population size of 4,000 adults and a carrying capacity of 400,000 adults at the end of each generation.

The oviposition bias coefficient, b, was calculated as 1 - (number of adults in refuge/number of plants in refuge) and was bounded by 0.2 and 1.0, after carrying capacity had been imposed on adult populations. The proportion of total oviposition in refugia was calculated according to b X p where p is the proportional size of the refuge. Oviposition in the Bt crop was calculated as 1 - (b X p). Thus, a bias coefficient of 1.0 was equivalent to random oviposition, while at a bias of 0.2, oviposition in refugia was reduced five-fold relative to random expectation. Field observations showed that oviposition rates in the summer growing seasons were reduced by two- to five-fold relative to random expectation.

We simulated a range of management scenarios that imposed different forms of density-dependent mortality on insect larvae. Density independent mortality followed the constant survival of 14% [see Additional file 1: Table S2]. Spray action thresholds for control of pests in refugia were calculated according to the number of egg masses per plant, assuming each female laid two egg masses of 150 eggs with a maximum of one egg mass per plant. These were informed by agriculture extension service recommendations for S. frugiperda across North America on the proportion of infested field corn plants (typically 25% to 75% of plants with damage) [5052]. There was clear evidence for density-dependent larval mortality in experimental refugia from field observations of mortality from egg to third instar [see Additional file 2: Figure S1]. Mortality imposed by natural enemies in the refugia was based on a power law fitted to these data. We halved the overall survivorship relative to field conditions to account for the death of later instars, pupae and adults so that survival = 0.01x-1.274 where x is the total density of adults per plant in the whole crop, calculated after the imposition of carrying capacity. The field data could not reliably indicate whether overall mortality in the Bt crop was entirely density independent or incorporated a density-dependent component. We, therefore, also simulated a scenario in which density-dependent mortality occurred throughout the crop, with the survival of susceptibles calculated as = 0.00025x-1.274 [see Additional file 2: Figure S1] on the Bt portion of the crop. After the imposition of fitness costs, the mortality of resistant larvae was unaffected by toxin, as indicated by previous work [21].


  1. Hutchison WD, Burkness EC, Mitchell PD, Moon RD, Leslie TW, Fleischer SJ, Abrahamson M, Hamilton KL, Steffey KL, Gray ME, Hellmich RL, Kaster LV, Hunt TE, Wright RJ, Pecinovsky K, Rabaey TL, Flood BR, Raun ES: Areawide suppression of European corn borer with Bt maize reaps savings to non-Bt maize growers. Science. 2010, 330: 222-225. 10.1126/science.1190242.

    Article  CAS  PubMed  Google Scholar 

  2. Carriere Y, Ellers-Kirk C, Sisterson M, Antilla L, Whitlow M, Dennehy TJ, Tabashnik BE: Long-term regional suppression of pink bollworm by Bacillus thuringiensis cotton. Proc Natl Acad Sci U S A. 2003, 100: 1519-1523. 10.1073/pnas.0436708100.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  3. Christou P, Capell T, Kohli A, Gatehouse JA, Gatehouse AM: Recent developments and future prospects in insect pest control in transgenic crops. Trends Plant Sci. 2006, 11: 302-308. 10.1016/j.tplants.2006.04.001.

    Article  CAS  PubMed  Google Scholar 

  4. James C: Global status of commercialized Biotech/GM crops: 2012. Int Serv Aquisition Agri-Biotech Appl. 2012, Brief 44–2012.

    Google Scholar 

  5. Tabashnik BE, Gassmann AJ, Crowder DW, Carriere Y: Insect resistance to Bt crops: evidence versus theory. Nat Biotechnol. 2008, 26: 199-202. 10.1038/nbt1382.

    Article  CAS  PubMed  Google Scholar 

  6. Carrière Y, Crowder DW, Tabashnik BE: Evolutionary ecology of insect adaptation to Bt crops. Evol Appl. 2010, 3: 561-573. 10.1111/j.1752-4571.2010.00129.x.

    Article  PubMed Central  PubMed  Google Scholar 

  7. Huang F, Andow DA, Buschman LL: Success of the high-dose/refuge resistance management strategy after 15 years of Bt crop use in North America. Entomol Exp Appl. 2011, 140: 1-16. 10.1111/j.1570-7458.2011.01138.x.

    Article  Google Scholar 

  8. Gould F: Sustainability of transgenic insecticidal cultivars: integrating pest genetics and ecology. Ann Rev Entomol. 1998, 43: 701-726. 10.1146/annurev.ento.43.1.701.

    Article  CAS  Google Scholar 

  9. Tabashnik BE, Gassmann AJ, Crowder DW, Carriere Y: Field-evolved resistance to Bt toxins. Nat Biotechnol. 2008, 26: 1074-1076.

    Article  CAS  Google Scholar 

  10. Downes S, Parker T, Mahon R: Incipient resistance of Helicoverpa punctigera to the Cry2Ab Bt toxin in Bollgard II (R) cotton. PLoS One. 2010, 5: 9-

    Google Scholar 

  11. Storer NP, Babcock JM, Schlenz M, Meade T, Thompson GD, Bing JW, Huckaba RM: Discovery and characterization of field resistance to Bt maize: Spodoptera frugiperda (Lepidoptera: Noctuidae) in Puerto Rico. J Econ Entomol. 2010, 103: 1031-1038. 10.1603/EC10040.

    Article  PubMed  Google Scholar 

  12. Gassmann AJ, Petzold-Maxwell JL, Keweshan RS, Dunbar MW: Field-evolved resistance to Bt maize by western corn rootworm. PLoS One. 2011, 6: e22629-10.1371/journal.pone.0022629.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  13. Kruger M, Van Rensburg J, Van den Berg J: Resistance to Bt maize in Busseola fusca (Lepidoptera: Noctuidae) from Vaalharts. S Afr Environ Entomol. 2011, 40: 477-483. 10.1603/EN09220.

    Article  CAS  Google Scholar 

  14. Hellmich R, Higgins L, Witkowski J, Campbell J, Lewis L: Oviposition by European corn borer (Lepidoptera: Crambidae) in response to various transgenic corn events. J Econ Entomol. 1999, 92: 1014-1020.

    Article  Google Scholar 

  15. Kjaer C, Damgaard C, Lauritzen A: Assessment of effect of Bt-oilseed rape on large white butterfly (Pieris brassicae) in natural habitats. Entomol Exp Appl. 2010, 134: 304-311. 10.1111/j.1570-7458.2009.00958.x.

    Article  Google Scholar 

  16. Lei Z, Liu T, Greenberg S: Feeding, oviposition and survival of Liriomyza trifolii (Diptera: Agromyzidae) on Bt and non-Bt cottons. Bull Entomol Res. 2009, 99: 3-

    Article  Google Scholar 

  17. De Moraes C, Mescher M, Tumlinson J: Caterpillar-induced nocturnal plant volatiles repel conspecific females. Nature. 2001, 410: 577-580. 10.1038/35069058.

    Article  CAS  PubMed  Google Scholar 

  18. Harmon JP, White JA, Andow DA: Oviposition behavior of Ostrinia nubilalis (Lepidoptera: Crambidae) in response to potential intra- and interspecific interactions. Environ Entomol. 2003, 32: 334-339. 10.1603/0046-225X-32.2.334.

    Article  Google Scholar 

  19. Andrews K: The whorlworm, Spodoptera frugiperda, in Central America and neighboring areas. Fla Entomol. 1980, 63: 456-467. 10.2307/3494530.

    Article  Google Scholar 

  20. Hardke JT, Leonard BR, Temple JH: Fall armyworm oviposition on cotton plants expressing wide strike (TM) Bollgard (R), and Bollgard II (R) Cry proteins. Southwest Entomol. 2012, 37: 295-303. 10.3958/059.037.0305.

    Article  Google Scholar 

  21. Storer N, Kubiszak M, King JJ, Thompson G, Santos A: Status of resistance to Bt maize in Spodoptera frugiperda: lessons from Puerto Rico. J Invert Pathol. 2012, 110: 294-300. 10.1016/j.jip.2012.04.007.

    Article  Google Scholar 

  22. Penaflor MF, Erb M, Robert CA, Miranda LA, Werneburg AG, Alda Dossi FC, Turlings TCJ, Bento JM: Oviposition by a moth suppresses constitutive and herbivore-induced plant volatiles in maize. Planta. 2011, 234: 207-215. 10.1007/s00425-011-1409-9.

    Article  CAS  PubMed  Google Scholar 

  23. Signoretti AG, Penaflor MF, Bento JM: Fall Armyworm, Spodoptera frugiperda (JE Smith) (Lepidoptera: Noctuidae), female moths respond to herbivore-induced corn volatiles. Neotrop Entomol. 2012, 41: 22-26. 10.1007/s13744-011-0003-y.

    Article  CAS  PubMed  Google Scholar 

  24. Sparks AN: Review of the biology of the fall armyworm (Leipodopera, Noctuidae). Fla Entomol. 1979, 62: 82-87. 10.2307/3494083.

    Article  Google Scholar 

  25. Penaflor MF, Erb M, Miranda LA, Werneburg AG, Bento JM: Herbivore-induced plant volatiles can serve as host location cues for a generalist and a specialist egg parasitoid. J Chem Ecol. 2011, 37: 1304-1313. 10.1007/s10886-011-0047-9.

    Article  CAS  PubMed  Google Scholar 

  26. Fatouros NE, Lucas-Barbosa D, Weldegergis BT, Pashalidou FG, van Loon JJ, Dicke M, Harvey JA, Gols R, Huigens ME: Plant volatiles induced by herbivore egg deposition affect insects of different trophic levels. PLoS One. 2012, 7: 8-

    Article  Google Scholar 

  27. Andow DA, Farrell SL, Hu Y: Planting patterns of in-field refuges observed for Bt maize in Minnesota. J Econ Entomol. 2010, 103: 1394-1399. 10.1603/EC09201.

    Article  CAS  PubMed  Google Scholar 

  28. Onstad DW, Buschman LL: Evaluation of oviposition deterrence in management of resistance to transgenic corn by European corn borer (Lepidoptera: Crambidae). J Econ Entomol. 2006, 99: 2100-2109. 10.1603/0022-0493-99.6.2100.

    Article  CAS  PubMed  Google Scholar 

  29. Polanczyk RA, Alves SB: Biological parameters of Spodoptera frugiperda (J.E. Smith) (Lepidoptera: Noctuidae) assayed with Bacillus thuringiensis berliner. Sci Agr. 2005, 62: 464-468. 10.1590/S0103-90162005000500009.

    Article  Google Scholar 

  30. Berdegue M, Trumble JT, Moar WJ: Effect of CryIC toxin from Bacillus thuringiensis on larval feeding behavior of Spodoptera exigua. Entomol Exp Appl. 1996, 80: 389-401. 10.1111/j.1570-7458.1996.tb00951.x.

    Article  CAS  Google Scholar 

  31. Huang CH, Yan FM, Byers JA, Wang RJ, Xu CR: Volatiles induced by the larvae of the Asian corn borer (Ostrinia furnacalis) in maize plants affect behavior of conspecific larvae and female adults. Insect Sci. 2009, 16: 311-320. 10.1111/j.1744-7917.2009.01257.x.

    Article  CAS  Google Scholar 

  32. Martins CB, Zarbin PH: Volatile organic compounds of conspecific-damaged Eucalyptus benthamii influence responses of mated females of Thaumastocoris peregrinus. J Chem Ecol. 2013, 39: 602-611. 10.1007/s10886-013-0287-y.

    Article  CAS  PubMed  Google Scholar 

  33. Reisenman CE, Riffell JA, Duffy K, Pesque A, Mikles D, Goodwin B: Species-specific effects of herbivory on the oviposition behavior of the moth Manduca sexta. J Chem Ecol. 2013, 39: 76-89. 10.1007/s10886-012-0228-1.

    Article  CAS  PubMed  Google Scholar 

  34. Zakir A, Sadek MM, Bengtsson M, Hansson BS, Witzgall P, Anderson P: Herbivore-induced plant volatiles provide associational resistance against an ovipositing herbivore. J Ecol. 2013, 101: 410-417. 10.1111/1365-2745.12041.

    Article  Google Scholar 

  35. Alphey N, Coleman PG, Donnelly CA, Alphey L: Managing insecticide resistance by mass release of engineered insects. J Econ Entomol. 2007, 100: 1642-1649. 10.1603/0022-0493(2007)100[1642:MIRBMR]2.0.CO;2.

    Article  CAS  PubMed  Google Scholar 

  36. Degen T, Dillmann C, Marion-Poll F, Turlings TC: High genetic variability of herbivore-induced volatile emission within a broad range of maize inbred lines. Plant Physiol. 2004, 135: 1928-1938. 10.1104/pp.104.039891.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  37. Kos M, van Loon J, Dicke M, Vet L: Transgenic plants as vital components of integrated pest management. Trends Biotechnol. 2009, 27: 621-627. 10.1016/j.tibtech.2009.08.002.

    Article  CAS  PubMed  Google Scholar 

  38. Carrière Y, Tabashnik BE: Reversing insect adaptation to transgenic insecticidal plants. Proc Biol Sci. 2001, 268: 1475-1480. 10.1098/rspb.2001.1689.

    Article  PubMed Central  PubMed  Google Scholar 

  39. Ives AR, Glaum PR, Ziebarth NL, Andow DA: The evolution of resistance to two-toxin pyramid transgenic crops. Ecol Appl. 2011, 21: 503-515. 10.1890/09-1869.1.

    Article  PubMed  Google Scholar 

  40. Ghimire MN, Huang F, Leonard R, Head GP, Yang Y: Susceptibility of Cry1Ab-susceptible and -resistant sugarcane borer to transgenic corn plants containing single or pyramided Bacillus thuringiensis genes. Crop Prot. 2011, 30: 74-81. 10.1016/j.cropro.2010.09.008.

    Article  CAS  Google Scholar 

  41. Roush RT: Two-toxin strategies for management of insecticidal transgenic crops: can pyramiding succeed where pesticide mixtures have not?. Philos Trans R Soc Lond B Biol Sci. 1998, 353: 1777-1786. 10.1098/rstb.1998.0330.

    Article  CAS  PubMed Central  Google Scholar 

  42. Brévault T, Heuberger S, Zhang M, Ellers-Kirk C, Ni X, Masson L, Li X, Tabashnik BE, Carrière Y: Potential shortfall of pyramided transgenic cotton for insect resistance management. Proc Natl Acad Sci U S A. 2013, 110: 5806-5811. 10.1073/pnas.1216719110.

    Article  PubMed Central  PubMed  Google Scholar 

  43. Blanco CA, Andow DA, Abel CA, Sumerford DV, Hernandez G, López JD, Adams L, Groot A, Leonard R, Parker R, Payne G, Perera OP, Terán-Vargas AP, Azuara-Domínguez A: Bacillus thuringiensis Cry1Ac resistance frequency in tobacco budworm (Lepidoptera: Noctuidae). J Econ Entomol. 2009, 102: 381-387. 10.1603/029.102.0149.

    Article  CAS  PubMed  Google Scholar 

  44. Sappington T: First-flight adult European Corn Borer (Lepidoptera: Crambidae) distribution in roadside vegetation relative to cropping patterns and corn phenology. Environ Entomol. 2005, 34: 1541-1548. 10.1603/0046-225X-34.6.1541.

    Article  Google Scholar 

  45. Rabí O, Permuy N, Garcia E: Nueva variedad de maíz FR-28 con alto potencial de rendimiento. Agrotecnia. 1997, 27: 43-44.

    Google Scholar 

  46. Ashley T, Wiseman B, Davis F, Andrews K: The fall armyworm: a bibliography. Fla Entomol. 1989, 72: 152-202. 10.2307/3494982.

    Article  Google Scholar 

  47. The R Project for Statistical Computing.,

  48. Crawley MJ: The R Book. 2007, Chichester: John Wiley

    Book  Google Scholar 

  49. Raymond B, Sayyed AH, Hails RS, Wright DJ: Exploiting pathogens and their impact on fitness costs to manage the evolution of resistance to Bacillus thuringiensis. J Appl Ecol. 2007, 44: 768-780. 10.1111/j.1365-2664.2007.01285.x.

    Article  Google Scholar 

  50. Bessin R: Fall armyworm in corn.,

  51. Insects and Pests of Field Crops: Corn Insects and Pests. []

  52. Boyd ML, Bailey WC: Management of the armyworm complex in Missouri field crops. Available at:

Download references


This work was supported by Grant 3031–243 from BioCubaFarma and by an NERC fellowship (to BR). This study is a contribution from the Imperial College Grand Challenges in Ecosystems and the Environment initiative. Bt (Cry1F) maize was grown under a National Centre for Biosafety license CH47-L (68)08. The authors thank Bruce Tabashnik, Nick Storer and two anonymous reviewers for useful suggestions.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Camilo Ayra-Pardo.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

PTR obtained the Bt-maize variety, prepared the documentation for the biosafety license to grow a Bt-maize in open field conditions, participated in the design of the study and performed the field experiments. BR analyzed the field data, conducted the simulation modeling and wrote the paper. IMB and LRC participated in obtaining the Bt-maize variety, and the biosafety license for experiments with Bt-maize in open field conditions and assisted with the field experiment. CGB participated in the study design and coordination of the field experiments. DJW participated in the study design and helped to draft the manuscript. CAP conceived and designed the study, and drafted the manuscript. All authors read and approved the final manuscript.

Pilar Téllez-Rodríguez, Ben Raymond contributed equally to this work.

Electronic supplementary material


Additional file 1: Table S1: Comparisons of crop damage; data are means from a four-point damage scale (0- no damage, 4- heavily damaged). Kruskal Wallis tests were conducted for data from no more than two time points per season: at around 20 days after planting and also on the last day any significant damage was detected. Key: **** P <0.0001, *** P <0.001. Table S2 . Estimated survivorship from egg to third instar larva based on field observations over four independent plantings of Bt maize and conventional refuge. Survivorship is based on females laying an average of 150 eggs per egg mass. Table S3. Summary of simulation model parameters. Note: survival is given before effects of fitness costs are calculated. An oviposition bias of 1 is equivalent to random egg-laying behavior. (DOC 64 KB)


Additional file 2: Figure S1: (A) Density-dependent survivorship of fall armyworm in the conventional maize refuge. Estimates of survivorship (diamonds) are based on scouting of larvae and on the assumption that each egg mass contains 150 eggs. (B) Estimates of larval survival in the experimental Bt crop (asterisks) showing the fitted density-independent and density-dependent mortality used in simulations. (PDF 91 KB)


Additional file 3: Figure S2: The impact of dominance of resistance on the evolution of resistance in computer simulations when there is density-dependent mortality throughout the refuge and the Bt crop. (A) Under random oviposition increasing refuge size and decreasing dominance slow the evolution of resistance. (B) Under damage-avoiding oviposition these factors appear to interact more: increasing refuge size had a greater impact on delaying resistance as the dominance of resistance decreased. (PDF 538 KB)


Additional file 4: Figure S3: The impact of population dynamics and biased oviposition on the evolution of resistance in computer simulations when resistance to Bt toxins is fully recessive. Here, we explored the relationship between female fecundity and refuge size under random (A) and damage-avoiding oviposition (B). As in previous models, simulation artefacts meant that the evolution of resistance was prevented in sink populations with small refugia. In general, larger refugia did not effectively prevent evolution of resistance when females were avoiding damaged plants. Simulation results were largely insensitive to female fecundity under random oviposition, while female fecundity was an important parameter under damage-avoiding oviposition. The spray action threshold was set at 0.25 egg masses per plant in all simulations. (PDF 533 KB)


Additional file 5: Figure S4: The impact of population dynamics and biased oviposition on the evolution of resistance in computer simulations when resistance to Bt toxins is under partially dominant (h = 0.1, A, B) and fully recessive (h = 0) resistance (C,D). Here, we explored the relationship between the density dependent action threshold for spraying refugia (in egg masses per plant) and refuge size under random and damage-avoiding oviposition . As in previous models, larger refugia had a limited impact on the evolution of resistance when females were avoiding damaged plants. This result was relatively insensitive to spray thresholds, except when refugia were very large. Note that sprays are rarely applied with spray thresholds of 0.8 or above. Once again, small refugia could lead to the simulated populations becoming sinks, which also prevented the evolution of fully recessive resistance. Fecundity was set at 300 eggs per female in all models. (PDF 698 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( ) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Téllez-Rodríguez, P., Raymond, B., Morán-Bertot, I. et al. Strong oviposition preference for Bt over non-Bt maize in Spodoptera frugiperda and its implications for the evolution of resistance. BMC Biol 12, 48 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: