The fall armyworm, Spodoptera frugiperda(J. E. Smith) (Lepidoptera: Noctuidae), is currently the most important maize pest in Mexico. Its control is mainly based on the use of conventional insecticides. Additionally, Bt-maize expressing Cry1F protein represents an alternative to control this pest. We estimated the baseline susceptibility in Mexican populations of S. frugiperda to Cry1F protein.Twenty-eight geographical populations were field collected from Baja California Sur, Chihuahua, Coahuila, Durango, Sinaloa, Sonora, and Tamaulipas states.The F1 neonate larvae of each population were subjected to diet-overlay bioassay. After 7 d of Cry1F exposure, the percent mortality and the percent growth inhibition with respect to the untreated control were recorded (S-LAB).The LC50 ranged from 14.4 (6.3–24.0) (Cajeme 1, Sonora) to 161.8 ng/cm2 (92.0–320) (Ahumada 2, Chihuahua), while the LC95 was between 207.1 (145–363) (Obregón, Sonora) and 1,217 ng/cm2 (510.8–7,390.0) (Río Bravo 2,Tamaulipas). The sensitivity ratios at 50% mortality, (LC50 field/LC50 S-Lab) and 95% mortality were ≤6.45 and ≤5.05-fold, respectively.The 50% growth inhibition (GI50) ranged from 2.8 (0.008–9.3) (Obregón, Sonora) to 42.4 ng/cm2 (3.6–147.0) (Cajeme 1, Sonora).The GI95 was between 75.4 (San Luis Río Colorado, Sonora) to 1,198 ng/cm2 (Cajeme 1, Sonora).The relative inhibition at 50% of the growth, (RI50 = GI50 field /GI50 S-LAB) was ≤3.5 and at 95% (RI95) was ≤1.91-fold.These results indicated susceptibility to Cry1F protein in the evaluated populations of S. frugiperda.
The fall armyworm, Spodoptera frugiperda (J. E. Smith) (Lepidoptera: Noctuidae), is one of the most important pests of maize (Torres and Cotes 2005). This species originated in the Americas and spread to the tropical and subtropical regions of the world (Andrews 1988, Clavijo and Perez Greiner 2000). It has currently dispersed in the Western hemisphere (Storer et al. 2012; Velez et al. 2013, 2016; Bernardi et al. 2015) and more recently in Africa (Bateman et al. 2018), India (CGIAR 2018), and China (Jing et al. 2019, USDA-FAS 2019). This species has been observed feeding on 186 hosts located in 42 botanical families (Pashley 1986). Moreover, S. frugiperda can cause damage to late-planted crops in subtropical regions and throughout the cropping season in the tropical regions (Chandrasena et al. 2018).
In Mexico, S. frugiperda is widely distributed (Pacheco 1985) and considered the most economically important pest of maize. It can reduce yield from 45% (Hruska and Gladstone 1988) to 100% (Cruz et al. 1999). The application of organo-synthetic insecticides has been the most-used strategy to control this pest (Pacheco 1994). However, insecticide resistance in S. frugiperda is currently reported in Americas to 29 active ingredients belonging to six insecticide groups (Mota-Sanchez and Wise 2017). It is estimated that up to 1,152 tons of pyrethroid are used in México against S. frugiperda (Terán-Vargas 2008) and annually, around 3,000 tons of organo-synthetic insecticides are sprayed against this pest (Blanco et al. 2014). Despite the abundant use of chemical compounds, the fall armyworm still causes significant yield loss. As an alternative approach, developers of Bt technology have been generating transgenic maize plants that express genes from the bacterium Bacillus thuringiensis Berliner that confer resistance to economically important pests such as the fall armyworm.
Dow AgroSciences LLC and Pioneer Hi-Bred International, Inc. developed the GM-maize event TC1507 (Herculex I Insect Protection), which expresses the Cry1F protein derived from Bacillus thuringiensis var. aizawai (Bt). TC1507 provides control of different foliar, ear-feeding and stalk boring lepidopteran including S. frugiperda and has been adopted in major maize production countries such as the United States, Argentina, and Brazil (Storer et al. 2012, Chandrasena et al. 2018).
Bt maize has potential to effectively control the fall armyworm and significantly reduce the use of insecticides in almost 4 million acres in Mexico. Since October 2013, a moratorium on the use of Bt maize was established limiting its cultivation, including field experiments. However, farmers anticipate a positive resolution to this matter, given the lack of scientific evidence related to adverse effects derived from the use of this technology (NAS 2017). Meanwhile, the efforts have been centered on conducting the studies supporting commercialization such as establishing baselines for susceptibility of the target pest populations, to enable future development of resistance management programs in Mexico.
To date, insect pests have demonstrated the ability to develop resistance to Bt crops around the globe (Tabashnik et al. 2003). For example, the level of resistance to the Cry1F protein was so high in Puerto Rico that no fall armyworm larvae exposed to 10,000 ng/ cm2 were killed. Field resistance is influenced heavily by the strength of the implemented resistance management program (Siegfried et al. 2007), importance of nonstructured refuge (Shelton and Zhao 2009), number of maize growing seasons per year (Storer 2010), and even the risk of resistance allele immigration mediated by wind movement across wide areas (Drake and Farrow 1988) or commercial trade of infested agricultural products (Nagoshi et al. 2018).
Deployment of a sound resistance management program is of importance to extend the durability of this technology. Before a Bt crop is deployed in the field, it is helpful to determine the baseline susceptibility of primary pests to estimate the natural geographical variation in susceptibility to specific Bt proteins expressed by the genetically modified crop. This information constitutes a reference that enables the measurement of potential changes in population susceptibility once these technologies are used in the field and selection for resistance occurs. Therefore, the objective of this study was to determine the baseline response to the Cry1F protein derived from B. thuringiensis in neonate larvae of S. frugiperda populations from the main maize-producing Mexican States of Baja California Sur, Chihuahua, Coahuila, Durango, Sonora, and Tamaulipas.
Materials and Methods
Populations
In total, 28 geographical populations of S. frugiperda were field-collected from non-Bt commercial field maize in Baja California Sur, Chihuahua, Coahuila, Durango, Sinaloa, Sonora, and Tamaulipas (Fig. 1) between 2013 and 2017 (Table 1). For each population, five maize fields with a distance of approximately ≥15 km between them were selected. At least five sampling points were randomly selected per field, and the maize plants were visually inspected to detect those with new damage. Then, a single more than or equal to third instar was field-collected per infested plant. In total, 280–560 larvae were field-collected per population and placed individually in cups containing 8.0 ml of meridic diet (Southland products, Lake Village, AR). A population which has been maintained approximately for 10 yr under pesticide-free laboratory conditions, with occasional input of larvae collected from maize in Montecillo, State of Mexico, where this species has not been subjected to insecticide applications, was used as a reference for comparison.
Insect Rearing
Once in the laboratory, the field-collected larvae were quarantined for 1 wk. Then, the dead, diseased, and parasitized larvae were discarded. The healthy larvae were individually transferred to a glass vial containing 10 ml of artificial diet prepared according to manufacturer's instructions (Insect Media for Spodoptera frugiperda, Southland Products, Inc., Lake Village, AR). In all cases, at least 90% of the field-collected larvae were healthy and reared to the adult stage to obtain at least 1,000 F1 healthy neonate larvae for each bioassay. Once the adults emerged, they were visually inspected to confirm species identity (Capinera 2000). Groups of 20 pairs of adults were placed in a 16-size brown paper bag, and 15% sucrose water was provided as a food source. Females laid their eggs on the inner surface of the bag, which were collected every other day and incubated to obtain neonate larvae. The rearing of this species was at 27 ± 1°C with a photoperiod of 14:10 (L:D) h and a relative humidity of 70–80%.
Protein
The insecticidal Cry1F protein of B. thuringiensis (90% a.i, Lot TSN304065) was provided by Dow AgroSciences de México S.A. de C.V. The lyophilized protein was stored under refrigeration at 4°C in a vacuum-sealed container that contained a desiccant. CAPS (N-cyclohexyl-3-aminopropanesulfonic acid) (Sigma–Aldrich, St. Louis, MO) buffer was used to dilute the protein to prepare the required concentrations.
Bioassay
In each well of the bioassay trays (Bio-Ba-128, C-D International, Pitman, NJ), 1.5 ml of meridic diet (Southland Products, Lake Village, AR) was deposited. After 2 h, once the diet solidified, 50 µl of the required concentration of the Cry1F protein was deposited using overlay method in each well to form a layer on the diet surface of 1.5 cm2. Two hours later, each well was infested with an active neonate larva (F1, 0 to 24 h of age and without previous exposure to meridic diet). At least seven concentrations that covered the range from 0 to 100% mortality determined by preliminary assays were used. For each concentration, 16 larvae were used per dose, with at least five replications. Each replication included an untreated control to which only 50 µl of the CAPS buffer was applied.
After 7 d of exposure to this protein, the following variables were determined: percentage mortality, percent inhibition of growth measured by comparing weight of surviving larvae in a given treatment with the untreated control. Those that did not respond to the stimulus when prodded with a dissecting needle or failed to develop to second instar (L2 stage) were counted as dead. The maximum accepted mortality in the control was 20% and treatment mortality was corrected for control mortality using Abbott's formula (Abbott 1925). For calculating the percentage of growth inhibition, as indicated by weight, all the surviving larvae in the control, as well as those in the respective concentrations, were separately grouped and weighed. Then, the average weight of the larvae was measured. Based on the comparison of the average weights, the relative weight of the larvae surviving on each concentration compared to the weight of the larvae surviving on control treatment was used to calculate the growth inhibition. Dead insects and larvae that did not develop to second instar (L2) were combined for calculating percent effective mortality with the argument that any instar failing to reach second instar by the duration indicate that they are unfit to survive further.
Statistical Analysis
The variable, percentage of effective mortality was analyzed using the Probit model (Finney 1971) to estimate LC50 and LC95 values, as well as 95% Fiducial confidence intervals (CI) by using the Proc Probit procedure of the SAS statistical package (SAS Institute 2008). The sensitivity ratios (SR) at the level of LC50 (SR50) and LC95 (SR95) were calculated by dividing the LC50 (LC95) values of the respective field population by the LC50 (LC95) of the susceptible population (S-LAB) as suggested by Storer et al. 2010. Growth inhibition was calculated as percent weight reduction between the treatment and control. GI50 and GI95 values were estimated by the Proc Probit procedure using the SAS statistical package (SAS Institute 2008). The relative inhibition ratios (RI) at the level of GI50 (RI50) and GI95 (RI95) were calculated by dividing the GI50 (GI95) values of the respective field population by the GI50 (GI95) of the susceptible population (S-LAB).
Results
The results by each state from which S. frugiperda were collected are presented below. The range within a parenthesis indicate 95% Fiducial CIs for the corresponding LC or EC value.
We detected differences, both at LC50 and LC95 level, between some field-collected populations and the susceptible laboratory population (S-LAB) (Tables 2 and 3). These differences were more common at the LC50 level, where 18 out of 28 (64.2%) evaluated field populations required higher dosage to reach the same level of mortality (50%). At the CL95, only two field population (Canachi-Sinaloa and Navojoa-Sonora) displayed differences with S-LAB; these two population showed differences with the susceptible population at both LC50 and LC95, with SR95 values of 3.23 and 7.1×, respectively (Table 2). We observed fewer differences regarding the growth inhibition values; only one field-collected population (Ahumada 1- Chihuahua) displayed differences at the GI50 level (2.42×) with the S-LAB reference (Table 3).
Table 1.
Origin of the populations of Spodoptera frugiperda field-collected in four states of Mexico
Baja California Sur
The LC50 and the LC95 were 28.9 (23.9–34.3) and 245.4 (180–370.2) ng/cm2, respectively (Table 2). The estimated concentration that inhibited 50% of growth in comparison to the untreated control (GI50) was 8.6 (0.02–23.0) ng/cm2, and the GI95 was 187.5 (66.8–132,841.0) ng/cm2 (Table 3).
Chihuahua
The LC50 ranged from 27.3 (21.0–34.1; Cuauhtémoc) to 161.8 (92.0–320.0; Ahumada 2) ng/cm2 and the LC95 was between 390.9 (317.9–507.4; Camargo) and 1,065 (273.4–158,954,188; Ahumada 1) ng/cm2 (Table 2) as shown with wide Fiducial limits. When compared with the S-LAB, these field-collected populations showed wide variation in the sensitive ratio at 50% mortality (SR50) from 1.09- to 6.45-fold, and the SR95 was between 2.09- and 4.42-fold (Table 2). The GI50 ranged from 12.0 (2.8–21.3; Cuauhtémoc) to 29.3 (23.1–36.2; Ahumada 2) ng/cm2, and the GI95 from 98.9 (49.6–984.4; Cuauhtémoc) to 658.6 (448.3–1,097.0; Ahumada 2) ng/cm2 (Table 3).
Coahuila
Only a single population was tested from this region. The LC50 for Torreón population was 39.9 (23.6–68.2) ng/cm2, and the LC95 was 214.6 (107.3–1,560.0) ng/cm2 (Table 2). When compared with the susceptible population, the SR50 of the population was 1.59-fold and the SR95 was 0.89-fold (Table 2). The GI50 and the GI95 values were 30.4 (3.9–77.4) and 661.9 (184.8–550,300.0) ng/cm2, respectively (Table 3). The Torreón population had no survivors beyond 300 ng/cm2.
Durango
A single population (Gómez Palacio; collected February 2016) was assayed from this region. The LC50 of Gómez Palacio population was 50.3 (44.0–57.0) ng/cm2 and the LC95 was 209.6 (168.0–276.0) ng/ cm2 (Table 2). The GI50 was 4.6 (0.0007–16.2) and the GI95, 376.3 (110.8–1,804,182) ng/cm2 (Table 3). At 600 ng/cm2, none of the treated larvae survived.
Sinaloa
Ten field-collected populations were evaluated from the state of Sinaloa. The LC50 ranged from 29.4 (22.9–37.0; Culiacán 2) to 107.2 (74.0–157.0; Culiacán 1) ng/cm2. The LC95 was between 261.5 (117.3–3,757.0; Los Mochis) and 702 (396.0–1,934.0; Culiacán 1) ng/cm2 (Table 2). The lowest GI50 value was 1.5 (0.2–3.8; Chinitos) ng/cm2 and the highest GI50 was 41.4 (16.0–85.5; Culiacán 2) ng/cm2 (Table 3). The GI95 ranged from 205.5 (82.7–10,380.0; Los Mochis) to 788.2 (221.8–166,846.0; Culiacán 1) ng/cm2 (Culiacán 1). All the populations evaluated had >92% mortality at the highest evaluated dose (600 ng/cm2).
Sonora
The LC50 ranged from 14.4 (6.3–24.0; Cajeme 1) to 128.1 (110.0–149.0; Navojoa) ng/cm2 and the LC95 was between 157.1 (85.3–571.0; Cajeme 1) and 778 (589–1,119.0; Navojoa) ng/cm2 (Table 2). The GI50 was between 2.8 (0.008–9.3; Obregón) to 42.4 (3.6–147.0; Cajeme 1) ng/cm2, and the GI95 varied from 75.4 (52.4–139.5; San Luis Río Colorado) to 1,198 (256.8–195,915,619.0; Cajeme 1) ng/ cm-2 (Table 3). At 600 ng/cm2, only three populations had live larvae (Cajeme 1, Obregón, and Navojoa). The SR50 varied from 0.23 (San Obregón population, Sonora) to 3.5 (Cajeme 1 population, Sonora). The SR95 ranged from 0.12 (San Luis Río Colorado, Sonora) to 1.91 (Cajeme 1, Sonora). In both cases, the state of Sonora had the populations with the highest variation in RI50(95).
Tamaulipas
The LC50 ranged from 37.8 (16.0–96.0; Díaz Ordaz 2) to 79.7(48.3–134.2; Rio Bravo 2) ng/cm2 and the LC95 was between 228.3 (183.0–307.0; Díaz Ordaz 1) and 1,217 (510.8–7,390.0; Rio Bravo 2) ng/ cm2 (Table 2). The GI50 was between 1.4 (0.2–3.7; Rivereña) and 7.0 (0.3–16.7; Rio Bravo 1) ng/cm2 (Table 3). The GI95 ranged from 57.7 (39.9–116.2; Díaz Ordaz 1) to 656.8 (364.9–1,650.0; Rio Bravo 2) ng/cm2 (Table 3).
At the (SR50–SR95) level, the highest variation in all of the field populations respect to the susceptible one (S-LAB) were as follows: Baja California (1.15–1.02), Chihuahua (6.45–4.42), Coahuila (1.59–0.89), Durango (2.0–0.87), Sinaloa (4.27–2.91), Sonora (5.1–3.23), and Tamaulipas (3.18–5.05) (Table 2)
Discussion
Susceptibility to Cry1F across the various regions was indicated by mortality and growth parameters as discussed. To date, the established baseline is the first record of regional susceptibility of S. frugiperda to Cry1F in Mexico. RR value ≥10 could be indicating insect resistant development to insecticides as defined by Tabashnik et al. (2009, 2014) and Mota-Sanchez et al. (2002). Therefore, these populations appeared susceptible to the Cry1F protein because none of them reached that status.
Table 2.
Mortality of neonate larvae of Spodoptera frugiperda from several Mexican States, exposed for seven days to the Cry1F protein of Bacillus thuringiensis
In several regions of Latin America such as Argentina, Brazil, Colombia, Puerto Rico, Uruguay and Mexico, S. frugiperda is still a devastating pest in conventional and Bt maize (Blanco et al. 2014). Resistance to Cry1F (TC1507) in this pest was first documented in Puerto Rico during 2006 (Storer et al. 2010). Recently, Cry1F resistance and reduced susceptibility in S. frugiperda have been identified in Florida (Huang et al. 2014) and North Carolina (Li et al. 2016). Storer et al. (2012) found that resistant fall armyworm population from Puerto Rico showed no mortality when exposed to the Cry1F protein, even at the highest dosage evaluated (10,000 ng/cm2). In some areas of South America, resistance to Cry1F was detected after several years of commercialization, such as Brazil (>5,000-fold; Farias et al. 2014a,b) and Argentina (survival of 16.6–97% to the diagnostic concentration of 2,000 ng/cm2; Chandrasena et al. 2018). In Puerto Rico, resistance was characterized as autosomal and recessive (Storer et al. 2010), similar to the inheritance of resistance found in Florida (Huang et al. 2014), Brazil (Farias et al. 2014a,b), and Argentina (Chandrasena et al. 2018). Although flights are unlikely to happen between more distant geographies (Nagoshi et al. 2007a,b; 2018), migratory nature of this pest may also enable spread of resistant individuals to Mexico from nearby geographies, such as Brazil and Argentina, where resistance has been confirmed and where environmental conditions are not a limitation to their dispersal.
The highest LC95 obtained from all six provinces in Mexico ranged between 209 and 1217 ng/cm2 indicating S. frugiperda susceptibility to Cry1F in the geography. Additionally, a diagnostic concentration (LC99) derived from such a field-collected population can be considered as a field-relevant concentration to detect field evolved resistance. Brazil and Argentina reported a concentration of 2,000 ng/cm2 as a discriminating dose (based on upper confidence limit of the highest LC99) to distinguish S. frugiperda populations developing resistance to Cry1F (Farias et al. 2014a,b; Chandrasena et al. 2018) in the field. Despite the observed variability, our results did not indicate resistance to Cry1F protein in Mexico as several nearby geographies.
Table 3.
Percentage of weight inhibition compared to the untreated control of Spodoptera frugiperda neonate larvae from several states of Mexico, exposed for 7 d to the Cry1F protein of Bacillus thuringiensis
The sustainability of a Bt crop depends in part on the effectiveness of the implemented insecticide resistance management program (Tabashnik et al. 2003) and involves baseline susceptibility estimation, implementation of resistance monitoring plan, monitoring, and remediation. The characterization of susceptibility in a population is a complex issue given the fact that the variation in response to the Bt protein depends on the complex interaction of the genes with its environment. In this case, we considered that larval response is not be mediated by the expression of the alleles for resistance since they are assumed to be at very low frequency, and generally are recessive. Also, the bioassays are carried out with a specific meridic diet to provide uniform conditions, thus minimizing variation in susceptibility imposed by factors such as differences in nutritional source where the target pest feeds under field conditions (Carrie et al. 2016, Deans et al. 2016). Other factors such as bioassay method, mortality criteria, type of parameter evaluated (mortality or growth inhibition) and characteristics of toxicological lot(s) should remain similar in order to minimize the variation in susceptibility due to external factors (Saeglitz et al. 2006, Bird and Akhurst 2007, Blanco et al. 2008b, Storer et al. 2010). Immigration of resistant alleles might be another source of variation in the response to Bt proteins since the fall armyworm displays long distance flights assisted by wind (Nagoshi et al. (2017). Under laboratory conditions, baseline susceptibility is determined by phenotypic bioassays (Blanco et al. 2008b) or by estimating the frequency of alleles for resistance (Blanco et al. 2008a, Farias et al. 2016). The phenotypic bioassays consist of either depositing a layer of Bt protein (diet-overlay) on the surface of the diet (Farias et al. 2014a,b) or mixing the protein with the diet (diet-incorporated). The overlay method requires less amount of protein, nevertheless irregularities in the layer deposition on the diet could result in inconsistent results derived from nonuniform exposure of the larvae (Siegfried et al. 2007) and this method may not be appropriate when dealing with larvae that typically burrow into the diet to feed, such as the pink bollworm, Pectinophora gossypiella Saunders (Lepidoptera: Gelechiidae) (Patin et al. 1999). Phenotypic bioassays are not expected to detect small changes in the allele frequency (Andow and Alstad 1998) but are cheaper, use less protein, is less time-consuming and can provide information meaningful to assess changes in population susceptibility. The methodologies such as the F2 screen to estimate the frequency of alleles for resistance provide reliable information about the proportion of this type of allele in comparison to the susceptible counterpart, even though its partially recessive (Andow and Ives 2002), but the use of this option is expensive and time-consuming. However, regardless of the methodology deployed, laboratory research itself is not sufficient to detect resistance (Natália et al. 2016) but generates alerts that must be confirmed with field observations. It is recommended to conduct susceptibility evaluations at least annually in both, the laboratory and field (Shelton and Zhao 2009), in order to detect resistant individuals as early as possible (Dennehy 1987). Under field conditions, the performance of the pest in both, Bt and the conventional crop, should be evaluated in several locations before large-scale deployment. Phenotypic resistance is adequately detected in the field, rather than in the laboratory, when the larva survives on Bt plants, reaches the adult stage and transmits its genes to the next generation (Andow and Ives 2002). In order to provide a clearer concept, Tabashnik and Carrière (2017) defined field-evolved resistance as a ‘genetically based decreased in susceptibility of an insect population to a Bt toxin caused by selection in the field'.
Before Bt crop deployment in the field, it is considered that the alleles conferring resistance to those Bt genes are rare in the population, with a frequency in the order of 10–3 (Tabashnik 1994, Gould et al. 1997, Burd et al. 2003) making the homozygous resistant individuals rare and recessive with a high proportion of homozygous and heterozygous susceptible insects. Even though Bt maize is not grown in Mexico, we consider that the frequency of alleles for resistance is above the indicated value; wide-area transgenic crops in Texas are grown and the dispersal of individuals to host crops in the studied areas is possible. Additionally, illegal cultivation of Bt maize in Mexico is not discarded and may contribute to selection pressure.
Several elements of resistance management have to be implemented in addition to monitoring susceptibility to Cry1F in target pest populations with bioassays, to address the complexity of introducing Bt crops in countries such as Mexico. Resistance must be managed by a combination of tactics such as application of Integrated Pest Management practices, continuing to promote the use of refuge by outreach education programs for growers, introduction of pyramided Bt crops, and also involving proactive engagement with stakeholders, such as the academia, industry, regulatory authorities, and grower organizations.
Conclusion
The susceptibility to Cry1F across various regions were indicative of mortality and growth parameters discussed. The results indicated that the evaluated populations of S. frugiperda are susceptible to the Cry1F protein. The established baseline is the first record of regional susceptibility of S. frugiperda to Cry1F in Mexico to date. This multi-region assessment completes the investigation of baseline susceptibility to Cry1F in S. frugiperda in Mexico indicating potential of Cry1F containing maize to be an effective tool for the control of this pest.
Acknowledgements
We like to extend our sincere thanks to academic research and technical laboratory staff who participated in various activities; Chad Boeckman (Corteva Agriscience) and Mike Connelly (Pioneer Hi-Bred) for their support to coordinate this project.