Evaluation of the Impact of Genetically Modified Cotton After 20 Years of Cultivation in Mexico

For more than 20 years cotton has been the most widely sown genetically modified (GM) crop in Mexico. Its cultivation has fulfilled all requirements and has gone through the different regulatory stages. During the last 20 years, both research-institutions and biotech-companies have generated scientific and technical information regarding GM cotton cultivation in Mexico. In this work, we collected data in order to analyze the environmental and agronomic effects of the use of GM cotton in Mexico. In 1996, the introduction of Bt cotton made it possible to reactivate this crop, which in previous years was greatly reduced due to pest problems, production costs and environmental concerns. Bt cotton is a widely accepted tool for cotton producers and has proven to be efficient for the control of lepidopteran pests. The economic benefits of its use are variable, and depend on factors such as the international cotton-prices and other costs associated with its inputs. So far, the management strategies used to prevent development of insect resistance to GM cotton has been successful, and there are no reports of insect resistance development to Bt cotton in Mexico. In addition, no effects have been observed on non-target organisms. For herbicide tolerant cotton, the prevention of herbicide resistance has also been successful since unlike other countries, the onset of resistance weeds is still slow, apparently due to cultural practices and rotation of different herbicides. Environmental benefits have been achieved with a reduction in chemical insecticide applications and the subsequent decrease in primary pest populations, so that the inclusion of other technologies—e.g., use of non-Bt cotton- can be explored. Nevertheless, control measures need to be implemented during transport of the bolls and fiber to prevent dispersal of volunteer plants and subsequent gene flow to wild relatives distributed outside the GM cotton growing areas. It is still necessary to implement national research programs, so that biotechnology and plant breeding advances can be used in the development of cotton varieties adapted to the Mexican particular environmental conditions and to control insect pests of regional importance.

Introduction

Cotton is one of the most important natural sources for fiber, oil, and seeds for livestock feeding. All the cotton produced in the world is obtained from four domesticated species of the Gossypium genus of the Malvaceae family. With 18 species, Central, and South America are the richest regions in Gossypium species globally, being Mexico one of the most diverse countries with 14 different species. The northeast of Africa and the southwest of Arabia also have 14 different species and Australia has 17 species (Cronquist, 1981; Fryxell, 1992; Percival et al., 1999).

An outstanding feature of cotton domestication is that it occurred simultaneously in different continents from local cotton wild ancestors. This process of parallel and convergent domestication occurred for the species Gossypium hirsutum in Mexico, G. barbadense in Peru, G. arboreum in Sudan and G. herbaceum in Pakistan. In each of these four cases, the unique properties of cotton fiber useful to make ropes and textiles were noticed thousands of years ago. From these four species, G. hirsutum, commonly referred to as Mexican cotton or highland cotton, is the most widely planted, accounting for 90% of world production. This is relevant, since Mexico is an important center of origin and domestication of many other cultivated crops, such as corn, squash, pumpkin, bean, and chilies. Currently in Mexico several native cotton species are present, including the wild relatives of G. hirsutum. The highest concentrations of wild cotton relatives are located in the southeast region of the country, the only place where G. hirsutum is found as a common species in the native flora (Coppens d'Eeckenbrugge and Lacape, 2014).

Before the deployment of GM technology, cotton production was associated to high environmental, economic, and sanitary costs due to the necessity of large amounts of pesticide applications. A different strategy was necessary to improve yields, thus technology involving GM cotton cultivars with inserted genes that confers resistance to lepidopteran pests and to herbicides was adopted by the growers (Deguine et al., 2008; Benbrook, 2012).

In Mexico as in other parts of the world, the cultivation of cotton was characterized by the application of large quantities of chemical insecticides. For example, in the 1970s decade, cotton cultivation required almost 20 applications of chemical insecticides from the plant emergence to harvest, since cotton plants must be protected from insect attack when the plant emerges, until the profitable bolls open (a period that lasts ~20 weeks). In the middle of the Twentieth century, at the peak of cotton production in Mexico, the cotton area that was planted reached 900,000 hectares with 2 million bales produced per year (the term “white gold” was used at that time to describe cotton). Years later, the increasing pest pressure and high doses of pesticides resulted in the evolution of insect resistance to chemical insecticides. In addition, reductions in international prices of the fiber resulted in a production decline due to unsustainable operating costs (Martínez-Carrillo and Díaz-López, 2005; Martínez-Carrillo, 2015).

In 1996, GM cotton was for the first time commercially planted in Mexico as well as in five other countries (James, 2016), due to the impossibility of cultivating conventional cotton in areas of severe pest pressure (Terán-Vargas et al., 2005). Since then, a total of 15 countries have commercialized GM cotton (Argentina, Australia, Burma, Brazil, Burkina Faso, China, Colombia, Costa Rica, United States, India, Mexico, Paraguay, Pakistan, South Africa, and Sudan). In Mexico, the increase in GM cotton adoption was gradual (Martínez-Carrillo, 2005), and since 2008 the 96% of the area cultivated with cotton was GM cotton (Purcell et al., 2008).

Nevertheless, the area planted with GM cotton in Mexico has fluctuated, depending on international fiber prices, input costs and the prevalence pests, weeds, and diseases. The main cotton production areas of Mexico are located in the northern region of the country. This region has an arid climate and growers used irrigation systems. These areas of cotton production are not in close proximity to areas containing wild relatives of cotton, as stated in the Mexican law (CIBIOGEM, 2018).

The transformation events or transgenes that have been authorized in Mexico since 1996 confer two main traits, one is the tolerance to herbicides and the other is the resistance to lepidopteran pests. In the first case, plants are tolerant to herbicides such as glyphosate (Nida et al., 1996), ammonium glufosinate (Blair-Kerth et al., 2001) and dicamba (Cahoon et al., 2015) that are used to combat weeds. In the second, resistance to lepidopteran pests is due to the insertion of cry genes from the bacterium Bacillus thuringiensis (Bt) that confers resistance to larval stages of different lepidopteran pest such as Pectinophora gossypiella, Helicoverpa zea, Heliothis virescens (Benedict et al., 1993), and Spodoptera exigua (Wilson et al., 1992; James, 2016).

In Mexico, the “Biosafety Law of Genetically Modified Organisms” regulates the cultivation of GM cotton and other biotech crops in a step-by-step and case-by-case basis. The different steps refer to the different stages of release: experimental, pilot and commercial plantings. Prior to the commercial release, the authorities evaluate the results of the experimental and pilot (semi-commercial) scale releases, carrying out risk assessment studies and examining the experimental results, as well as the compliance and effectiveness of the biosafety measures (DOF, 2005). Academic institutions must endorse the research carried out in Mexico. A total of 15 GM cotton release events were requested from 2005 to 2015, in 342 dossiers [Figure 1; (CIBIOGEM, 2018)].

FIGURE 1

Figure 1. Timeline for GM cotton transformation events released in Mexico. Black dots indicate the first and last year that each release was requested. The figure shows the different GM cotton released applications that have received a permit in Mexico, from 2005 until 2015 (CIBIOGEM, 2018).

The environmental risk assessment studies aim to identify potential damage to the environment where the level of risk is estimated, the potential negative effects are identified, and actions needed to reduce environmental risks are determined (EPA, 1998). In the case of the environmental risks associated with the release of agricultural GMOs, it is important to compare them with the risks associated to the agricultural practices used on conventional crops. This is why a “case by case” analysis should be performed, that is, to consider the modified organism, the intended use, and the likely environment and environmental conditions in which it will be grown. The risk assessment studies for the release of GM cotton in the case of Mexico included an evaluation of the risks of gene flow to wild relatives, the possible effect on non-target organisms, the risks of selection of resistant weeds to herbicides and the evolution of resistance to Cry proteins by the insect pests (SEMARNAT, 2018).

In this work, we present an updated analysis of the data available since the release of GM cotton in 1996. Two main hypothesis were questioned: the first hypothesis is if there is potential risk in gene flow to native species, while the second is if the use of GM-cotton in Mexico would result in a reduction of pesticides use and in higher yields.

Methods

Analysis of Wild Cotton Species Distribution

For the analysis of the wild cotton species distribution, we used the CONABIO database where 2,238 records were cured and verified (including 16 cotton species: G. thurberi, G. armourianum, G. harknessii, G. davidsonii, G. aridum, G. raimondii, G. gossypioides, G. lobatum, G. laxum, G. trilobum, G. turneri, G. schwendimanii, G. lanceolatum, G. hirsutum, and G. barbadense; CONABIO, 2018).

In order to assess the likelihood of gene flow, the cotton growing regions were characterized and a distribution model of wild G. hirsutum was constructed. The environmental characteristics of these cotton growing regions were identified by a classification tree, using as covariates of 19 current bioclimatic layers (Hijmans et al., 2005), 12 solar radiation layers (WorldClim), terrain slopes and ruggedness index.

Development of an Ecological Niche Models

To elaborate ecological niche models (ENM) of two different scenarios of cultivated cotton (without volunteers and with volunteers), we used a database constructed with 259 unique presence records of GM cotton plots and 17 records of cotton volunteers reported by several volunteer monitoring campaigns carried out in the cotton growing regions. Records from plots in the Northeast region (Tamaulipas) were not available and were not included in the analysis.

Nineteen current bioclimatic layers were downloaded from the WorldClim 1.4 data set (Hijmans et al., 2005) and six topographical layers from the HYDRO1k Elevation Derivative Database (available at: http://lta.cr.usgs.gov/HYDRO1K), using a resolution of 30 arc-s (ca. 1 km).

Maxent 3.3.3e (Phillips et al., 2006) runs were performed, one for each scenario. Each run included 30 replicates using the logistic model, and 20% random test by bootstrap. All the distribution models were evaluated using AUC scores (0.98 with and without volunteers). The models were transformed into binomial data, with a total presence value as the cut-off for each scenario (0.01 without volunteers and 0.15 with volunteers).

Surveys of Cotton Farmers

In Mexico, cotton farming is commonly managed by the owner of the land or the farmer that uses it, and a “technical advisor,” that is a professional pest control crop advisor.

In order to determine the perception of the Mexican farmers on the impacts of planting of GM cotton, a survey was designed and applied to 167 farmers in 20 municipalities of the main cotton-producing states. The objectives of the survey were to identify factors associated with the use of GM cotton in Mexico, to know the willingness of farmers to use this biotechnology and the perception of benefits or problems that they have observed, to identify changes in yields, production costs, control of pests, handling, and use of pesticides from the transition from conventional to GM cotton and to evaluate the indirect effects of the use of this technology on the environment and in human health. The survey was designed according to the methodology of agricultural surveys with multiple sampling frames and the sample design for the study of rural organizations in Mexico (Kish, 1990; González-Villalobos and Wallace, 1998). The margin of error of this survey was ±7.46% with a total estimated population of 5,000 cotton farmers and a confidence level of 95% (Survey System, 2018).

Surveys of Technical Advisors

A survey was applied to 165 technical advisors specialized in cotton management. This survey was based on Shaw et al. (2009), to assess the impact of GM-crops with tolerance to glyphosate. Questions related to the pest management were also added.

The technician advisors' sample size was: Mexicali (n = 46); Chihuahua (n = 39); and La Laguna (n = 80) (Figure 2). The margin of error of this survey was ± 7.5% and a confidence level of 95% (Survey System, 2018).

FIGURE 2

Figure 2. Climatic suitability model for geographical space projection regions for wild cotton G. hirsutum (green dots) and cultivated GM cotton regions (red dots). Regions do not overlap, but show proximity in the area known as “La Laguna” (black circle).

Results and Discussion

Exploring the First Hypothesis: Gene Flow From Cultivated Cotton (Conventional or Transgenic) to Wild Relatives

Since Mexico is a center of the origin and diversification of G. hirsutum, one of the main environmental concerns for the release of GM cotton was the possibility of transgene flow to native cotton populations (Ellstrand, 2002, 2012; Ellstrand et al., 2013).

In Mexico there is a continuum of G. hirsutum cotton varieties that range from wild, feral and locally domesticated to improved varieties, therefore the potential for gene flow among them exists if they coexist in the same area. To assess such risk, it is necessary to know the geographic distribution patterns of the different varieties, and also the dispersal mechanisms of the species. The geographical distribution of wild populations and cultivated cotton was taken into account in the risk assessment evaluation and the geographical separation constitutes one of the conditions in México for the release of GM cotton into the environment and before sowing field visits were done to identify the possible presence of wild cotton relatives (BCH, 2018; SAGARPA, 2018).

The geographical overlap between native species distribution and the region in which GM cotton is currently planted is minimal, according to the records of the “National Commission for Knowledge and Use of Biodiversity” (CONABIO, 2018). The delimited GM cotton growing regions correspond to semi-arid regions (Figure 2, red dots) that do not geographically overlap with the area of climatic suitable zones of wild G. hirsutum. However, they are close to the La Laguna region (Figure 2). Few GM cotton regions were not included in the analysis either due to security issues or restrictions in technical support (i.e., North of Tamaulipas, Valleys of Yaqui, and Mayo, and Planicie Huasteca), all of them coincided with the climatic suitability zones of G. hirsutum. Nevertheless, according to the National Statistics (INEGI, 2012), Tamaulipas is the state with less cotton production in the country and the Yaqui valley as well as the Planicie Huasteca are not even in the statistics of cotton production.

For gene flow through pollen to occur, it is not only required that the plants coexist in the same area and that they are compatible, but also that the pollen containing transgenes is dispersed via pollinators. In the case of cotton, the rate of cross-pollination (the probability that a plant is pollinated with pollen from other plant) is 10% or less, since 90% of the plants resulted from self-pollination (Meredith and Bridge, 1973; Llewellyn and Fitt, 1996; Sen et al., 2004; Van Deynze et al., 2005; Zhang et al., 2005). It was also reported that, in cases where cross-pollination by bees occurs, it significantly decreases with the distance between plants. High cross-pollination probability occurs only when plants are located in close proximity (Umbeck et al., 1991; Yan et al., 2015). Moreover, the cross-pollination rate depends, to a large extent, on the climatic and ecological condition that determine, for example, the patterns of activity and abundance of insect species carrying out pollination and pollen flow (Llewellyn et al., 2007).

However, in our study we observed that the most imminent risk of gene flow is not by pollen, but by seeds spilled during transportation. Cotton-seeds can be efficiently dispersed by either wind or water. During several field visits to the cotton productions areas, it was observed that there is a very strict control and biosafety measures during the movement of the GM cotton-seeds from the seed-companies to the fields. The GM seeds arrive in closed packages and closed vehicles. However, after the harvest, such controls relaxed, and the seeds are transported to the gins in open vehicles that spill seeds in the roads. Volunteer plants can grow from spilled seeds and have been observed in the edge of roads. Sanity authorities and seed companies are in charge of removing the volunteer plants, but unnoticeable escapes are always possible.

From the two scenarios of cultivated cotton (without volunteers and with volunteers), we further elaborate an ENM as described in Methods. Figure 3 shows the Principal Component Analysis (PCA) of environmental conditions of the analyzed cotton records (wild, GM, and volunteer). It can be seen that the conditions in which GM cotton is planted (blue dots) are very restrictive and conditions are clearly differentiable from the rest of the cotton species (wild in black, gray, and colors). However, the presence of GM volunteers (red dots) in environments other than GM growing regions demonstrates the environmental plasticity of GM cotton, and broadens the environmental component of the GM cotton niche toward the environmental space occupied by wild species. In Figure 4 we show the potential distribution of GM and wild cotton. According to the models describing the two possible scenarios (without and with volunteers), this figure shows that the presence of volunteers significantly expands the niche of GM cotton in its geographic component (Figure 4).

FIGURE 3

Figure 3. Principal Component Analysis (PCA) of environmental conditions of the analyzed cotton records (wild, GM, and volunteer). Conditions in which GM cotton is planted (blue dots) are very restrictive and differentiable from the rest of the cotton species (wild in black, gray, and colors). GM volunteers are represented with red dots.

FIGURE 4

Figure 4. GM cotton distribution models adding the presence of volunteers. In green: potential distribution of GM cotton. In yellow: potential distribution of wild cotton. Black dots indicate the records of GM cotton plots while red dots indicate the records of volunteer plants used for the elaboration of the models.

It is important to mention that Wegier et al. (2011) reported the existence of gene flow at long distances between cultivated and wild populations of G. hirsutum, by the identification of recombinant proteins in wild populations of cotton. These authors proposed that the gene flow may be possible through the dispersion of seeds (Wegier et al., 2011). Hence, it is necessary to follow up the monitoring of hybrid populations and implement sensitive methods such as RT-PCR and digital-PCR to evaluate in detail the changes in transgene frequencies in these populations (Holst-Jensen, 2009; Fraiture et al., 2015; Randhawa et al., 2016).

What Do People That Work With the GM Cotton in Mexico Think

Surveys of Cotton Farmers

Overall, farmers pointed out that the use of GM cotton resulted in better pest control and easier pest management. Also, higher yields of GM cotton were generally mentioned. The reasons for stopping the planting of non-Bt conventional seed include difficulty for controlling pests and high costs of insecticides. According to the opinion of the farmers, GM cotton showed higher yields and required less use of insecticides and crop management. Nevertheless, according to farmers' opinions GM cotton-seeds are expensive and the use of herbicides is higher. In addition, farmers agreed that the highest yields of GM cotton are due to better seed quality and favorable weather conditions.

Cotton is planted in the arid areas of northern Mexico, where adverse weather conditions are prevalent, including the lack of water, extreme temperatures, drought, and frost. Inputs such as special planting equipment, irrigation, and fertilizers result in high production costs. In addition, an increase in seed prices, machinery, and fuels in recent years exacerbated the production costs.

The high operation costs as well as fluctuations in international fiber prices, led to a large fluctuation in the total cotton area planted. For instance, in 2016 the total cotton area in Mexico was reduced to 104,000 ha, due to the decrease in international prices and the increase in input costs. However, the cotton area was doubled to 210,000 ha in 2017 due to an increase in international fiber prices. The decrease in grain prices could be another important factor that favors cotton growing for some farmers.

Despite the cost of production, 80% of the farmers are highly satisfied with the use of the GM varieties, since the lepidopteran pests are controlled and excellent weed control is obtained. The remaining 11% of farmers are moderately satisfied, and 9% are not satisfied. Ten percent of the farmers considered that GM cotton is not profitable.

Interestingly, 40% of the farmers would be willing to plant conventional seeds if available in Mexico (conventional seeds are not produced now in Mexico), because it is assumed by these farmers that those seeds would cost less. Furthermore, due to current pest populations observed for the past few years, they considered that current pests are not necessarily controlled by GM varieties.

From the point of view of the effects on human health, farmers have a positive perception about the adoption of GM cotton. They believe that the intoxication cases due to chemical pesticide exposure have been reduced with the adoption of GM cotton. They reported less intoxication cases due to a lower use of chemical insecticides (Nava-Camberos et al., unpublished results).

Surveys of the Technical Advisors

In order to analyze changes in pest and weed management after the adoption of GM cotton a survey was applied to 165 technical advisors specialized in cotton management.

With respect to the management of weeds and herbicides, the responses of the technicians indicated that glyphosate is practically applied to the entire cotton growing area in Mexico at least once during the production cycle. The main weed species associated with cotton are field bindweed Convolvulus arvensis L., annual morning glories Ipomoea hederacea Jacq. and Ipomoea purpurea (L.) Roth, palmer amaranth Amaranthus palmeri S. Wats, johnsongrass Sorghum halepense (L.) Pers. and various annual grasses, mainly barnyardgrass Echinochloa colona (L.) Link.

According to these surveys, weed management in cotton in Mexico generally consists of the application of glyphosate that is complemented by deep tillage for soil preparation and in-row cultivation in more than 90% of the cotton area. The application of other herbicides in addition to glyphosate in pre-planting and pre-emergence is done in about 21% of the area where trifluralin represents the most used herbicide in these early applications.

Technicians indicated that problems associated with weed management were reduced in Mexicali and La Laguna, but they were increased in the state of Chihuahua, where the control of weeds with glyphosate was qualified as low. Sixty two percent of the technicians indicated that they have observed changes in the response of weeds to glyphosate. This response of the weeds implies the need of a dose increase of herbicides in order to have an effective control in the most difficult weeds. Nevertheless, 85% of the technicians are currently carrying out management practices to prevent the selection of glyphosate-resistant weeds, focusing mainly to in-row cultivation, hand weeding and crop rotation.

Before the use of Bt cotton, the Lepidoptera complex (P. gossypiella, H. zea, H. virescens, and S. exigua) comprised the majority (ca. 60%) of the total reported pests; followed by sucking insects (whitefly, Chlorochroa ligata and Lygus; ca. 20%). The reported insects list is presented in Table 1, where it is observed this drastic drop in lepidopteran counts, while other insects such as aphids, mites, weevils, thrips, and whiteflies increased in counts by the technicians. The technicians consider that the pressure of the Lepidoptera complex was very high before the use of GM-cotton and now it has effectively been reduced.

TABLE 1

Table 1. Insects reported by technicians as important pests in all regions.

After 20 years of using Bt-cotton, the interviewed technical advisors have observed drastic changes in the composition of insect pest species. Currently, the most important are Anthonomus grandis, C. ligata, Bemisia tabaci, several species of sucking insect pests, and thrips. The Lepidoptera complex represented only up to 5% of the reported pests (mentioned by 0, 0, 0, and 5% of the technical advisors in Mexicali, Sonora, La Laguna, and Chihuahua, respectively) while the sucking insect pests comprised around 73% (60, 60, 80, and 95% of the survey in Sonora, La Laguna, Chihuahua, and Mexicali, respectively). Due to environmental differences in the cotton growing regions of Mexico, it is difficult to rank the overall importance of pests. For example, whiteflies are of primary importance in Mexicali, Sonora, and La Laguna, but in Chihuahua, they are considered a secondary pest. Conchuela (C. ligata) is still considered the primary pest in La Laguna and Chihuahua, but it is not a concern in Mexicali and Sonora. A. grandis once a menacing pest throughout Mexico, currently is only important in La Laguna and Sonora, but in Mexicali and Chihuahua this pest is eradicated. This eradication is due to the joint A. grandis eradication Mexico-USA program. After using Bt-cotton, P. gossypiella, H. virescens, and Bucculatrix thurberiella now have very low population levels in the different cotton regions. H. zea and S. exigua are currently considered pests of secondary importance in all cotton areas (Table 1).

Regarding the number of total insecticide applications, the technicians reported a significant decrease due to the use of GM cotton. Due to the effectiveness of Bt-cotton, and its high rate of adoption in most of the growing areas, in Chihuahua and La Laguna the synthetic insecticides sprays have been reduced to 3.5 and 5.0 applications, respectively, from the previous ~12 applications used in a crop season. Nevertheless, in other regions such as Mexicali and Sonora that showed high pressure of pests that are not targeted by Bt-cotton (whiteflies, Lygus bugs, and boll weevils) the insecticide sprays are still high.

Exploring the Second the Hypothesis: Effects and Impacts of GM Cotton Cultivation in Mexico

Different lines of evidence indicated that the use of GM cotton has contributed to reducing the number of insecticide applications necessary to achieve adequate control of lepidopteran pests in the cotton regions of Mexico. Cotton is one of the crops in which the greatest amount of pesticides is applied in the world, so the alternative of using Bt-cotton represents an advantage from the environmental point of view (Abedullah et al., 2015). It is known that the use of pesticides can have negative impacts on the quality of water and soil, human health, aquatic species, and beneficial insects and other organisms (Boatman et al., 2004; Arias-Estevez et al., 2008; Athukorala et al., 2012).

According to most farmers, GM cotton in Mexico, despite its costs, is still economically profitable and is one of the main income sources in the municipalities where it is planted. In those places, GM cotton seems to ensure production, and prevent losses by lepidopteran insect pests, while reducing costs and labor activities as well as the use of vehicles to spray pesticides (Skevas et al., 2013). The impact on crop yield has also been significant since in Chihuahua, La Laguna and Mexicali the yield increments are 1.8, 2.4, and 3.7 bales per ha, respectively, which is equivalent to increases of $8,700, $11,500, and $17,700 Mexican pesos per ha.

It is difficult to illustrate the agronomic advances that the cotton industry has experienced in recent decades without also involving factors such as the improvement of seeds, the better use of water and fertilizers. Great effects are the result of better training of the agricultural technicians and government campaigns for crop health. Pest eradication is an additional benefit of this technology. For example, since 2007 it has not been necessary to apply insecticides against P. gossypiella in Chihuahua. It is calculated that the P. gossypiella-eradication program resulted in 1.7 million less liters of chemicals saving of more than 207 million Mexican pesos for cotton producers (CESAVECH, 2015).

Few studies have analyzed the effect on human health and the environment of GM cotton. Adoption of Bt-cotton reduced acute pesticide poisoning in farmers in China and India (Hossain et al., 2004; Kouser and Qaim, 2011). The compounds present in the pesticides used in conventional crops tend to accumulate in human tissues and are very dangerous for workers if the appropriate safety equipment is not used.

As mentioned before, different data and our surveys indicate that the intensity with which pesticides were used before GM cotton was very high. The intense use of broad-spectrum insecticides in conventional cotton was highly toxic, since those compounds affect many kinds of animals, including humans, and usually have high permanence in the field, affecting food chains of predators, parasitoids, and pollinator insects.

Ecological and Evolutionary Aspects of GM Cotton

Effect of Bt-cotton on Non-target Insects

Annual crops such as cotton require a field season comprised of 6–7 months and involve the intensive management of both weeds and insect pests. The Cry toxins produced by Bt that are expressed in different cotton events (Bt cotton) are specific to insects of the order Lepidoptera. These toxins are active against common cotton pests such as P. gossypiella, H. zea, H. virescens, and S. exigua. Thus, the control of other pests of different insect orders that attack cotton such as the coleopteran A. grandis, or the hemipteran B. tabaci or other insect pests still require applications of synthetic insecticides.

It is important to note that formulated insecticides based on Bt are used in integrated pest management (IPM) and organic agriculture because of their high specificity. Bt is also integrated into pest management, due to its biodegradable nature and ability to control specific pests, lacking impact on non-target organisms such as bees, parasitoid wasps, earthworms, beneficial true bugs, or predatory beetles, which do not possess an active target site (or receptor) where the Bt protein can interact (Pardo-López et al., 2013). The results of numerous studies with Bt toxins show that when non-target organisms are exposed to Bt toxins in similar amounts, or higher than those produced by the Bt-crops, they are not affected (Zwahlen et al., 2003; Ferry et al., 2005; Lu et al., 2010; Schuler et al., 2013). Among the most detailed studies are those in which a pest (e.g., an aphid, mite or worm) is fed on Bt-cotton and is consequently consumed or parasitized by a predator/parasitoid without any effect on the non-target insect (Zwahlen et al., 2003; Ferry et al., 2005; Lu et al., 2010; Schuler et al., 2013).

Due to the high effectiveness of Bt cotton against the most important lepidopteran pests, the damage induced by these Lepidoptera complex in Bt cotton is substantially smaller, or non-existent, when compared with the damage that they produced on conventional cotton if they were not controlled by chemical insecticides. However, the reduction of lepidopteran pests in Bt cotton may result in an increase of other cotton pests that are not controlled by Bt cotton. This phenomenon has been observed worldwide (Wang et al., 2006, 2009; Zhao et al., 2011) suggesting that secondary pests can occupy the resources previously used by lepidopteran insects. However, it was also reported that the lower use of chemical insecticides promotes the increase of natural enemies than can decrease populations of other non-target pests (Tian et al., 2015).

This increase of secondary pests apparently has been erroneously interpreted as an undesired effect of Bt cotton (Wang et al., 2008; Li et al., 2011; Zhao et al., 2011). Nevertheless, farmers generally control outbreaks of secondary pests with broad-spectrum insecticides. This practice, although effective against the target insects, also kills beneficial organisms.

It has also been shown that populations of non-target organisms may fluctuate in conventional cotton fields compared to those of Bt cotton, since the density of a pest may have consequences on the abundance of predators and parasitoids (Romeis et al., 2006). The reduced applications of the broad-spectrum pesticides may favor the increase of beneficial insect populations. However, a lower number of lepidopteran eggs and larvae in Bt cotton can affect the availability of food and hosts of natural enemies. Since the vast majority of these biological control agents have broad diets, the decrease in eggs, and larvae of lepidopteran insects affects their populations only temporarily (Theiling and Croft, 1988; Bradbury and Coats, 1989; Pisa et al., 2015).

Considering the ongoing controversy regarding the environmental impact of Bt cotton and particularly the scarce information on its effects on the diversity of the non-target insects under Mexican conditions, a study was carried out comparing arthropod populations in non-Bt and Bt cotton in the states of Durango and Coahuila (known as “La Laguna;” Nava-Camberos et al., unpublished results). Key target pests H. zea and S, exigua were only abundant in non Bt-cotton, while no differences were found in overall arthropod species composition and abundance between conventional and Bt-cotton areas. Among them, insects of three orders (Hemiptera, Thysanoptera, and Diptera) and three families (Aleyrodidae, Anthocoridae, and Thripidae) were the most abundant. At the trophic level, the total number of entomophagous and phytophagous insects was similar in both types of cotton. However, the non-Bt cotton presented a reduced diversity index, after several applications of insecticides (Nava-Camberos et al., unpublished results).

Evolution of Resistance in Insects

One of the most important economic risks of genetically modified crops is the evolution of resistance to Cry proteins by insects (Tabashnik et al., 2008) and to herbicides by weeds (Powles, 2008; Heap, 2018). In the case of Bt crops, the evolution of resistance to these crops has already been reported in different parts of the world in Bt corn and Bt cotton that express a single Cry protein (Tabashnik et al., 2008, 2013) or two Cry proteins (Jurat-Fuentes et al., 2003), although there has been no report of such resistance in Mexico (Tamez, 2010; Aguilar-Medel et al., 2017; Mota-Sanchez and Wise, 2018).

One strategy to delay the evolution of resistance is the deployment of “refuges,” which consist of plots with non-Bt plants near GM crops (Georghiou and Taylor, 1977; Gould, 1998). For the refuge strategy to be effective, insect resistance should be recessive (Carrière et al., 2010). This means that a resistant insect must carry two copies of the resistant allele. Heterozygous individuals with just one copy of the recessive allele are sensitive to a Cry toxin present in Bt cotton, and only homozygous individuals carrying two copies of the resistant alleles survive on the Bt plants. Therefore, the refuge has the purpose of maintaining a healthy population of susceptible insects. The idea is that when homozygous susceptible insects from the refuge mate with the resistant from the Bt crop field, their progeny will be heterozygous, meaning that they will have one susceptible allele, and one resistant allele. If this occurs effectively in the fields, the pest population will remain sensitive to the Cry toxin expressed in the Bt crop (Andow and Alstad, 1998).

Nevertheless, if two heterozygous insects mate, $\frac{1}{4}$ of their progeny will be resistant. For this reason it was suggested that in addition to the refuge strategy, the stacking of two or more cry genes that have different modes of action has been widely used to delay the evolution of resistance to Bt crops. For example, the stacked MON-88913-8 X MON-15985-7 event expresses the Cry1Ac and Cry2Ab toxins, which have been shown to have a different mode of actions, as they recognize distinct protein receptors in the guts of the same sensitive larvae (Caccia et al., 2010). The consequence that Cry1Ac and Cry2Ab recognize different proteins in the target pests, greatly reduces the probability of having a pest with double mutation (Caccia et al., 2010).

The eradication program of P. gossypiella implemented in the United States and Mexico since 2002 established the use of Cry toxins in conjunction with other control strategies. The adoption of Bt cotton with dual toxins by the local farmers resulted in the dramatic decline of this insect and its practically eradication in the Northern region of Mexico (SAGARPA, 2012, 2016; Martínez-Carrillo, 2015).

However, the secondary lepidopteran pest, S. exigua shows low susceptibility to Cry1A and Cry2A toxins, and recently it is causing significant damages to the Bt cotton crop in Mexico. To overcome this issue a new stacked event containing the vip3Aa gene plus cry1A, and cry2Ab (Kurtz et al., 2007; Carrière et al., 2015) might be deployed. Vip3A is a highly effective Bt protein that exhibits high toxicity against S. exigua, and it has a different mechanism of action than Cry proteins, thus these new pyramided events expressing also Vip3A could effectively control S. exigua (Lee et al., 2003; Chakroun et al., 2016). Therefore, this new stacked-Bt cotton variety has a wider spectrum of control than the previous ones, and it will be very helpful in insecticide resistance management. However, due to the high usage of Bt cotton in the American continent, the eventual evolution of resistance, even to the newly stacked events, cannot be ruled out. Therefore, it is necessary to continue searching for novel insecticidal proteins with different modes of action and high efficacy against different cotton pests.

Evolution of Resistance in Weeds

Regarding resistance to herbicides, the first cotton events used in Mexico and elsewhere contained glyphosate resistance genes, which caused an intense use of this herbicide in fields of GM cotton in very large areas of the planet, with the consequence of the evolution of resistance to glyphosate by a diversity of weeds (Powles, 2008). Currently, there are 40 weed species already resistant to glyphosate (Heap, 2018). For this reason, it is recommended the use GM cotton resistant to alternative herbicides with different mechanism of action and other integrated weed management practices that would allow an effective control of weeds, avoiding the evolution of herbicide resistance.

It is interesting that in Mexico there are no reports of weed resistance to the herbicides used in GM cotton (SENASICA, 2016; Heap, 2018). This may be due to the fact that Mexican cotton farmers commonly use conventional tillage and in-row cultivation. Adoption of no-tillage systems in herbicide-resistant GM crops seems to be part of the problem of evolution of herbicide-resistant weeds in countries such as USA, Brazil, and Argentina (Powles, 2008). To cope with this problem there has been a worldwide request to release events that have more than one gene of resistance to different herbicides such as ammonium glufosinate and glyphosate or glyphosate and dicamba.

Thus, in Mexico, the deep tillage along with manual removal of early weeds, in-row cultivation, and crop rotation have apparently delayed the appearance of glyphosate-resistant weeds despite the fact that GM cotton technology has been adopted for more than 15 years (CIBIOGEM, 2018). In contrast, in the United States the first case of Palmer amaranth A. palmeri resistant to glyphosate was reported in 2005 (Culpepper et al., 2006), only 8 years after this technology was adopted (Norswhorty et al., 2016).

Conclusions, Recommendations, and Perspectives

G. hirsutum is a native species in Mexico, from which several highly efficient GM cultivars have been developed for the production of cotton worldwide, and some of them are now used in the north region of Mexico.

The tetraploid cotton G. hirsutum has a relatively large genome and diverged from its diploid ancestors several million years ago (Shan et al., 2016). Due to the distribution and chromosomal composition of this species, it is expected that there is low risk of introgression or mixing with other diploid wild species of Mexico by pollen flow, but seeds represent an important risk. However, it is still possible that the mixing of GM cotton with wild populations of the same species or another tetraploid specie occurs. It is also possible that the effect of this introgression may be diluted in the wild by processes like meiotic drive or by the lack of selective pressure to maintain the GM genes in complex communities and if the GM genes represent a cost to carry and to express them. Nevertheless, direct experiments will be required to follow the introgressed plants for several generations in the field. Also, given the possibility of introgression is a potential risk, careful monitoring programs for transgenes should be maintained, in particular focusing on the fate and dispersal of the seeds due to spills that occur during transportation from the fields to the gins.

We need detailed socioeconomic studies, as well as epidemiological studies on the health of Mexican cotton farmers, as nowadays there is not enough data to conclude on those aspects.

So far no cases of weed resistance to glyphosate associated with cotton have been reported in Mexico (Heap, 2018). However, it is strongly recommended to encourage the use of appropriate management practices and alternative herbicides with different mechanism of action to delay the evolution of resistance to glyphosate (Devine et al., 1992). In cases of resistance, the use of GM-glyphosate resistant seeds should be avoided since there is a greater danger of increasing the populations of glyphosate-resistant weeds species. This has already occurred in the United States, where weeds such as Palmer amaranth, Johnson grass and barnyard grass are now resistant.

It is known that the use of herbicides with two or more modes of action significantly delays the evolution of herbicide resistance (Neve et al., 2011). Besides, it is necessary to continue integrating the use of herbicides with other management practices, such as deep tillage, in row cultivation, and crop rotation to diversify weed management and decrease selection pressure for herbicide resistance.

The impact of Bt cotton on the use of chemical insecticides has been significant. Since its introduction 20 years ago, there has been a decrease in the use of chemical insecticides, but the data varies between regions due to differences in the ecological and management conditions, different composition of pests and other non-target pests. The evolution of resistance in target-pests cannot be ruled out, even despite the proper use of refuges.

The reduction in the number of applications ranges from one application in Sonora and Mexicali, to almost five applications of chemical insecticide per crop cycle in La Laguna. Also, it is important that the chemical insecticides that are currently used to control the pest complex have, in average, a lower environmental impact than the ones used a couple of decades ago.

Despite the relative good news, it is necessary that farmers and cotton technicians continue to get involved in the detection of a possible loss of efficacy of Bt cotton against the target pests. It is very important also to maintain the active participation of farmers and technicians for the prevention of the evolution of resistance, particularly in the adequate implementation of refuge areas.

In the future, the integration of various pest management tactics will be important, such as cultural control through the destruction of crop residues and biological control through the use of natural enemies (entomopathogens, predators, and parasitoids). The monitoring of insect resistance to Cry toxins expressed by the approved cultivars and those that are envisaged for their introduction in the Mexican market should continue. Federal support for cotton producers is considered crucial to continue with the Binational (Mexico-USA) Program for the eradication of P. gossypiella and A. grandis, in several regions to declare more free zones in a short term.

The change in the composition of primary insect pests and the increasing possibility of the development of glyphosate-resistant weeds, suggest the urgent need of developing new biotechnological tools to meet national needs. Policies directed toward federal funding for scientific research in Mexico, as well as a national program of seed production should be also strongly encouraged. Mexico has now the human and scientific capabilities and consistent funding of long-term goals directed to a more sustainable agriculture is needed. This is particularly important due to the lack of possibilities for producers, since there is no national policy for seed production, which puts at risk not only cotton, but also the national food security. Today Mexico depends totally on seeds from the large international companies for its cotton production.

Mexico has been careful in observing the principles of the Cartagena Protocol and the national regulation is highly demanding and expensive to meet. However, in many cases these regulations can only be met by the large companies; as a result, researchers and national institutions with low budgets find impossible to comply with all the requirements established in the biosafety law.

Finally we strongly recommended the agricultural and scientific authorities of Mexico to support a healthy long-term program of national research in order to meet the new needs of agriculture, conventional or GM, for the next 20 years.

Ethics Statement

Either ethics approval or written consent are not necessary in our study because it is not a clinical study, but instead shows a collection of technical opinions of a group of experts and did not involve disclosure of sensitive personal data. According to the Declaration of Helsinki, medical research is subject to ethical standards that promote and ensure respect for all human subjects and protect their health and rights and it is addressed primarily to physicians. In this work the information obtained from technicians and farmers is related to the crop cultivation practices, and it did not involve any medical study.

Author Contributions

MR-M, VS, and AB coordinated the research. LE and MS planned the obtainment of the data. EN, SC, ES collected data and apply surveys. DM-S, ER-R, UN-C, JM-C, CB, MR-M, EN, and SC designed the technician's survey. SC and ES performed the geographical modeling. All authors contributed with the data analyses, discussion of the results, and writing of the paper.

Funding

Funding was provided by the Interministerial Commission of Biosafety of Genetically Modified Organisms in Mexico (CIBIOGEM-1000/655/2015).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The reviewer JI and handling Editor declared their shared affiliation.

Acknowledgments

This work was supported by the interministerial commission of Biosafety, (CIBIOGEM), Project 1000/655/2015 and in part by program PASPA-DGAPA, UNAM to LE and VS. We thank the help of Sol Ortiz, Laura Tovar, Pedro Macías, Jesús García-Feria, Aurora Ávila, Luis Omar Jimenez, Ricardo Mora, Dr. Ramón Cinco, and Ing. Verduzco and all the personnel that helped to data collection. To all the farmers and technicians that agreed to answer the surveys. To Victor Gutierrez, Jorge Martinez, Rolando Rios, and Gerardo Montejano for applying the farmers survey. To the cotton farmers from Sistema Producto Algodón, that provided valuable information. To the memory of Jorge Medina-Medina and Dr. Cándido Márquez, who showed great enthusiasm to collaborate in the project.

References

Abedullah, Kouser, S., and Qaim, M. (2015). Bt Cotton, pesticide use and environmental efficiency in Pakistan. J. Agric. Econ. 66, 66–86. doi: 10.1111/1477-9552.12072

CrossRef Full Text | Google Scholar

Aguilar-Medel, S., Rodríguez, J. C., Martínez-Carrillo, J. L., and Silva-Aguayo, G. (2017). Susceptibilities of geographic populations of Helicoverpa zea (Lepidoptera: Noctuidae) in Mexico to Bt ∂-endotoxins Cry1Ac and Cry2Ab: an 18-yr study. J. Econ. Entomol. 110, 2207–2216. doi: 10.1093/jee/tox203

PubMed Abstract | CrossRef Full Text | Google Scholar

Andow, D. A., and Alstad, D. N. (1998). F2 screen for rare resistant alleles. J. Econ. Entomol. 91, 572–578. doi: 10.1093/jee/91.3.572

CrossRef Full Text | Google Scholar

Arias-Estevez, M., Lopez-Periago, E., Martinez-Carballo, E., Simal-Gandara, J., Mejuto, J. C., and Garcia-Rio, L. (2008). The mobility and degradation of pesticides in soils and the pollution of groundwater resources. Agric. Ecosyst. Environ. 123, 247–260. doi: 10.1016/j.agee.2007.07.011

CrossRef Full Text | Google Scholar

Athukorala, W., Wilson, C., and Robinson, T. (2012). Determinants of health costs due to farmers' exposure to pesticides: an empirical analysis. J. Agric. Econ. 63, 158–174. doi: 10.1111/j.1477-9552.2011.00326.x

CrossRef Full Text | Google Scholar

BCH (2018). Biosafety Clearing House. Country Profile. Mexico. Risk Assessment. Available Online at: https://bch.cbd.int/about/countryprofile.shtml?country=mx (Accessed January 3, 2018).

Benbrook, C. (2012). Impacts of genetically engineered crops on pesticide use in the U.S. - the first sixteen years. Environ. Sci. Eur. 24:24. doi: 10.1186/2190-4715-24-24

CrossRef Full Text | Google Scholar

Benedict, J. H., Sachs, E. S., Altman, D. W., Ring, D. R., Stone, T. B., and Sims, S. R. (1993). Impact of δ-endotoxin-producing transgenic cotton on insect–plant interactions with Heliothis virescens and Helicoverpa zea (Lepidoptera: Noctuidae). Environ. Entomol. 22, 1–9. doi: 10.1093/ee/22.1.1

CrossRef Full Text | Google Scholar

Blair-Kerth, L. K., Dotray, P. A., Keeling, J. W., Gannaway, J. R., Oliver, M. J., and Quisenberry, J. E. (2001). Tolerance of transformed cotton to glufosinate. Weed Sci. 49, 375–380. doi: 10.1614/0043-1745(2001)049[0375:TOTCTG]2.0.CO;2

CrossRef Full Text | Google Scholar

Boatman, N. D., Brickle, N. W., Hart, J. D., Milsom, T. P., Morris, A. J., Murray, A. W. A., et al. (2004). Evidence for the indirect effects of pesticides on farmland birds. Ibis 146, 131–143. doi: 10.1111/j.1474-919X.2004.00347.x

CrossRef Full Text | Google Scholar

Bradbury, S. P., and Coats, J. R. (1989). Comparative toxicology of the pyrethroid insecticides. Rev. Environ. Contam. Toxicol. 108, 133–177. doi: 10.1007/978-1-4613-8850-0_4

PubMed Abstract | CrossRef Full Text | Google Scholar

Caccia, S., Hernández-Rodríguez, C. S., Mahon, R. J., Downes, S., James, W., Bautsoens, N., et al. (2010). Binding site alteration is responsible for field-isolated resistance to Bacillus thuringiensis Cry2A insecticidal proteins in two Helicoverpa species. PLoS ONE 5:e9975. doi: 10.1371/journal.pone.0009975

PubMed Abstract | CrossRef Full Text | Google Scholar

Cahoon, C. W., York, A. C., Jordan, D. L., Everman, W. J., Seagroves, R. W., Culpepper, A. S., et al. (2015). Palmer amaranth (Amaranthus palmeri) management in dicamba-resistant cotton. Weed Technol. 29, 758–770. doi: 10.1614/WT-D-15-00041.1

CrossRef Full Text | Google Scholar

Carrière, Y., Crickmore, N., and Tabashnik, B. E. (2015). Optimizing pyramided transgenic Bt crops for sustainable pest management. Nat. Biotechnol. 33, 161–168. doi: 10.1038/nbt.3099

PubMed Abstract | CrossRef Full Text | Google Scholar

Carrière, Y., Crowder, D. W., and Tabashnik, B. E. (2010). Evolutionary ecology of adaptation to Bt crops. Evol. Appl. 3, 561–573. doi: 10.1111/j.1752-4571.2010.00129.x

PubMed Abstract | CrossRef Full Text | Google Scholar

CESAVECH (2015). Comité Estatal de Sanidad del Estado de Chihuahua. Historia de la Erradicación del Gusano Rosado en el Estado. Programa de Supresión/Erradicación de Picudo y Gusano Rosado del Algodonero. Chihuahua: SAGARPA.

Chakroun, M., Banyuls, N., Bel, Y., Escriche, B., and Ferre, J. (2016). Bacterial vegetative insecticidal proteins (Vip) from entomopathogenic bacteria. Microbiol. Mol. Biol. Rev. 80, 329–350. doi: 10.1128/MMBR.00060-15

PubMed Abstract | CrossRef Full Text | Google Scholar

CONABIO (2018). Comisión Nacional para el Conocimiento y Uso de la Biodiversidad. Available Online at: http://www.conabio.gob.mx/remib_ingles/doctos/remibnodosdb.html (Accessed January 3, 2018).

CIBIOGEM (2018). Registro Nacional de Bioseguridad de OGMs. Available Online at: https://www.conacyt.gob.mx/cibiogem/index.php/sistema-nacional-de-informacion/registro-nacional-bioseguridad-ogms (Accessed January 29, 2018).

Coppens d'Eeckenbrugge, G., and Lacape, J. M. (2014). Distribution and differentiation of wild, feral, and cultivated populations of perennial upland cotton (Gossypium hirsutum L.) in Mesoamerica and the Caribbean. PLoS ONE 9:e107458. doi: 10.1371/journal.pone.0107458

PubMed Abstract | CrossRef Full Text | Google Scholar

Cronquist, A. (1981). An Integrated System of Classification of Flowering Plants. New York, NY: Columbia University Press.

Culpepper, A. S., Grey, T. L., Vencill, W. K., Kichler, J. M., Webster, T. M., Brown, S. M., et al. (2006). Glyphosate-resistant Palmer amaranth (Amaranthus palmeri) confirmed in Georgia. Weed Sci. 54, 620–626. doi: 10.1614/WS-06-001R.1

CrossRef Full Text | Google Scholar

Deguine, J. P., Ferron, P., and Russell, D. (2008). Sustainable pest management for cotton production. A review. Agron. Sust. Dev. 28, 113–137. doi: 10.1051/agro:2007042

CrossRef Full Text | Google Scholar

Devine, M., Duke, S. O., and Fedtke, C. (1992). Physiology of Herbicide Action. Englewood Cliffs, NJ: PTR Prentice Hall.

Google Scholar

DOF (2005). (Diario Oficial de la Federación) [Official Gazette of the Federation]. (2005, March 18). Ley de Bioseguridad de Organismos Genéticamente Modificados [Biosafety law for genetically modified organisms]. México. available online at: http://www.diputados.gob.mx/LeyesBiblio/pdf/LBOGM.pdf and https://www.conacyt.gob.mx/cibiogem/images/cibiogem/eng/Docs/Ing_LBOGM_P.pdf

Ellstrand, N. C. (2002). “Gene flow from transgenic crops to wild relatives: what have we learned, what do we know, what do we need to know?,” in Gene Flow Workshop (Columbus, OH: The Ohio State University).

Google Scholar

Ellstrand, N. C. (2012). “Over a decade of crop transgenes out-of-place,” in Regulation of Agricultural Biotechnology: The United States and Canada, eds C. A. Wozniak and A. McHughen (Dordrecht: Springer), 123–135.

Google Scholar

Ellstrand, N. C., Meirmans, P., Rong, J., Bartsch, D., Ghosh, A., de Jong, T. J., et al. (2013). Introgression of crop alleles into wild or weedy populations. Annu. Rev. Ecol. Evol. Syst. 44, 325–345. doi: 10.1146/annurev-ecolsys-110512-135840

CrossRef Full Text | Google Scholar

EPA (1998). Guidelines for ecological risk assessment. Fed. Regist. 63, 26846–26924.

Ferry, N., Mulligan, E. A., Majerus, M. E., and Gatehouse, A. M. (2005). Bitrophic and tritrophic effects of Bt Cry3A transgenic potato on beneficial, non-target, beetles. Transgenic Res. 16, 795–812. doi: 10.1007/s11248-007-9088-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Fraiture, M. A., Herman, P., Taverniers, I., De Loose, M., Deforce, D., and Roosens, N. H. (2015). Current and new approaches in GMO detection: challenges and solutions. Biom. Res. Int. 2015:392872. doi: 10.1155/2015/392872

PubMed Abstract | CrossRef Full Text | Google Scholar

Fryxell, P. A. (1992). A revised taxonomic interpretation of Gossypium, L. (Malvaceae). Rheedea 2, 108-165.

Google Scholar

Georghiou, G. P., and Taylor, C. E. (1977). Operational influences in the evolution of insecticide resistance. J. Econ. Entomol. 70, 653–658. doi: 10.1093/jee/70.5.653

PubMed Abstract | CrossRef Full Text | Google Scholar

González-Villalobos, A., and Wallace, M. A. (1998). Encuestas Agrícolas Con Múltiples Marcos de Muestreo. Rome: Food and Agriculture Organization of the United Nations.

Google Scholar

Gould, F. (1998). Sustainability of transgenic insecticidal cultivars: integrating pest genetics and ecology. Annu. Rev. Entomol. 43, 701–726. doi: 10.1146/annurev.ento.43.1.701

PubMed Abstract | CrossRef Full Text | Google Scholar

Heap, I. M. (2018). International Survey of Herbicide Resistant Weeds. Available online at: http://www.weedscience.org/In.asp (Accessed January 3, 2018).

Hijmans, R. J., Cameron, S. E., Parra, J. L., Jones, P. G., and Jarvis, A. (2005). Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 25, 1965–1978. doi: 10.1002/joc.1276

CrossRef Full Text | Google Scholar

Holst-Jensen, A. (2009). Testing for genetically modified organisms (GMOs): past, present and future perspectives. Biotechnol. Adv. 27, 1071–1082. doi: 10.1016/j.biotechadv.2009.05.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Hossain, M. D., Shitomi, Y., Moriyama, K., Higuchi, M., Hayakawa, T., Mitsui, T., et al. (2004). Characterization of a novel plasma membrane protein, expressed in the midgut epithelia of bombyx mori, that binds to Cry1A toxins. Appl. Environ. Microbiol. 70, 4604–4612. doi: 10.1128/AEM.70.8.4604-4612.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

INEGI (2012). Encuesta Nacional Agropecuaria. Instituto Nacional de Estadística y Geografía. Available online at: http://www.inegi.org.mx/est/contenidos/proyectos/encuestas/agropecuarias/ena/ena2012/ (Accessed May 17, 2018).

James, C. (2016). Global Status of Commercialized Biotech/GM Crops: 2016. Ithaca, NY: ISAAA Brief No. 52. ISAAA.

Google Scholar

Jurat-Fuentes, J. L., Gould, F. L., and Adang, M. J. (2003). Dual resistance to Bacillus thuringiensis Cry1Ac and Cry2Aa toxins in Heliothis virescens suggests multiple mechanisms of resistance. Appl. Environ. Microbiol. 69, 5898–5906. doi: 10.1128/AEM.69.10.5898-5906.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Kish, L. (1990). Métodos de Muestreo Para Las Encuestas Agricolas. Rome: Food and Agriculture Organization of the United Nations.

Google Scholar

Kouser, S., and Qaim, M. (2011). Impact of Bt cotton on pesticide poisoning in smallholder agriculture: a panel data analysis. Ecol. Econ. 70, 2105–2113. doi: 10.1016/j.ecolecon.2011.06.008

CrossRef Full Text | Google Scholar

Kurtz, R. W., McCaffery, A., and O'Reilly, D. (2007). Insect resistance management for Syngenta's VipCot (TM) transgenic cotton. J. Invertebr. Pathol. 95, 227–230. doi: 10.1016/j.jip.2007.03.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, M. K., Walters, F. S., Hart, H., Palekar, N., and Chen, J. S. (2003). Mode of action of the Bacillus thuringiensis vegetative insecticidal protein Vip3A differs from that of Cry1Ab delta-endotoxin. Appl. Environ. Microbiol. 69, 4648–4657. doi: 10.1128/AEM.69.8.4648-4657.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, G. P., Feng, H. Q., McNeil, J. N., Liu, B., Chen, P. Y., and Qiu, F. (2011). Impacts of transgenic Bt cotton on a non-target pest, Apolygus lucorum (Meyer-Dur) (Hemiptera: Miridae), in northern China. Crop Protect. 30, 1573–1578. doi: 10.1016/j.cropro.2011.08.015

CrossRef Full Text | Google Scholar

Llewellyn, D., and Fitt, G. (1996). Pollen dispersal from two field trials of transgenic cotton in the Namoi valley, Australia. Mol. Breed. 2, 157–166. doi: 10.1007/BF00441430

CrossRef Full Text | Google Scholar

Llewellyn, D. J., Tyson, C., Constable, G. A., Duggan, B., Beale, S., and Steel, P. (2007). Containment of regulated genetically modified cotton in the field. Agric. Ecosyst. Environ. 121, 419–429. doi: 10.1016/j.agee.2006.11.019

CrossRef Full Text | Google Scholar

Lu, Y., Wu, K., Jiang, Y., Xia, B., Li, P., Feng, H., et al. (2010). Mirid bug outbreaks in multiple crops correlated with wide-scale adoption of Bt cotton in China. Science 328, 1151–1154. doi: 10.1126/science.1187881

PubMed Abstract | CrossRef Full Text | Google Scholar

Martínez-Carrillo, J. L. (2005). “Evolución del algodón transgénico en México,” in VII Congreso Internacional de Ciencia Agrícolas (Mexicali: Universidad Autónoma de Baja California), 1–5.

Martínez-Carrillo, J. L. (2015). “El algodón GM en México,” in VIII Congreso Internacional en Ciencias Agrícolas (Mexicali: Universidad Autónoma de Baja California), 1166–1170.

Martínez-Carrillo, J. L., and Díaz-López, N. (2005). “Nine years of transgenic cotton in México, adoption and resistance management results,” in Beltwide Cotton Conferences (Memphis, TN: National Cotton Council of America), 1368–1372.

Meredith, W. R., and Bridge, R. R. (1973). Natural crossing in cotton (Gossypium hirsutum L.) in the Delta of Mississippi. Crop Sci. 13, 551–552. doi: 10.2135/cropsci1973.0011183X001300050016x

CrossRef Full Text | Google Scholar

Mota-Sanchez, D., and Wise, J. C. (2018). The Arthropod Pesticide Resistance Database. Michigan State University. Available online at: http://www.pesticideresistance.org

Neve, P., Norsworthy, J. K., Smith, K. L., and Zelaya, I. A. (2011). Modeling Glyphosate Resistance management strategies for Palmer amaranth (Amaranthus palmeri) in Cotton. Weed Technol. 25, 335–343. doi: 10.1614/WT-D-10-00171.1

CrossRef Full Text | Google Scholar

Nida, D. L., Kolacz, K. H., Buehler, R. E., Deaton, W. R., Schuler, W. R., Armstrong, T. A., et al. (1996). Glyphosate-tolerant cotton: genetic characterization and protein expression. J. Agric. Food Chem. 44, 1960–1966. doi: 10.1021/jf9505640

CrossRef Full Text | Google Scholar

Norswhorty, J. K., Schwartz, L. M., and Barber, L. M. (2016). The incidence and ramifications of glyphosate resistance in cotton. Outlooks Pest Manage. 27, 31–35. doi: 10.1564/v27_feb_07

CrossRef Full Text | Google Scholar

Pardo-López, L., Soberón, M., and Bravo, A. (2013). Bacillus thuringiensis insecticidal toxins: mode of action, insect resistance and consequences for crop protection. FEMS Microbiol. Rev. 37, 3–22. doi: 10.1111/j.1574-6976.2012.00341.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Percival, A. E., Wendel, J. F., and Stewart, J. M. (1999). “Taxonomy and germplasm resources” in Cotton Origin, History, Technology and Production, eds C. W. Smith and J. T. Cothren (New York, NY: Wiley Series in Crop Science), 33–63.

Google Scholar

Phillips, S. J., Anderson, R. P., and Schapire, R. E. (2006). Maximum entropy modeling of species geographic distributions. Ecol. Modell. 190, 231–259. doi: 10.1016/j.ecolmodel.2005.03.026

CrossRef Full Text | Google Scholar

Pisa, L. W., Amaral-Rogers, V., Belzunces, L. P., Bonmatin, J. M., Downs, C. A., Goulson, D., et al. (2015). Effects of neonicotinoids and fipronil on non-target invertebrates. Environ. Sci. Pollut. Res. 22, 68–102. doi: 10.1007/s11356-014-3471-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Powles, S. B. (2008). Evolved glyphosate-resistant weeds around the world: lessons to be learnt. Pests Manage. Sci. 64, 360–365. doi: 10.1002/ps.1525

PubMed Abstract | CrossRef Full Text | Google Scholar

Purcell, J. P., Rubinstein, C., and Vaughn, T. (2008). “El futuro de la biotecnología en latinoamérica,” in Cultivos Transgénicos Para la Agricultura Latinoamericana, ed C. Blanco (Mexico: Fondo de Cultura económica, Colección La ciencia para todos), 131–162.

Randhawa, G., Singh, M., and Sood, P. (2016). DNA-based methods for detection of genetically modified events in food and supply chain. Curr. Sci. 110, 1000–1009. doi: 10.18520/cs/v110/i6/1000-1009

CrossRef Full Text | Google Scholar

Romeis, J., Meissle, M., and Bigler, F. (2006). Transgenic crops expressing Bacillus thuringiensis toxins and biological control. Nat. Biotechnol. 24, 63–71. doi: 10.1038/nbt1180

PubMed Abstract | CrossRef Full Text | Google Scholar

SAGARPA (2012). “Acuerdo por el que se Declara Zona Libre de Gusano Rosado (Pectinophora gossypiella) y Picudo Del Algodonero (Anthonomus grandis) a Los Municipios de Juárez, Praxedis, G. Guerrero, Guadalupe, Ahumada, Janos, Ascensión, Nuevo Casas Grandes, Casas Grandes, Galeana y Buenaventura, en el Estado de Chihuahua”.

SAGARPA (2016). Acuerdo Por el que se Declara Como Zona Libre de Gusano Rosado (Pectinophora gossypiella) a los Estados de Baja California y Sonora.

SAGARPA (2018). Secretaria de Agricultura, Ganaderia, Desarrollo Rural, Pesca y Alimentacion. Available online at: https://www.gob.mx/sagarpa

SEMARNAT (2018). Secretaria de Medio Ambiente y Recursos Naturales. Available online at: http://www.semarnat.gob.mx/temas/gestion-ambiental/impacto-ambiental-y-tipos/estudios-de-riesgo-e-impacto-ambiental

SENASICA (2016). Servicio Nacional de Sanidad, Inocuidad y Calidad Agroalimentaria. Available online at: https://www.gob.mx/senasica/acciones-y-programas/acciones-estrategicas-de-sanidad-vegetal

Schuler, T. H., Potting, R. P. J., Denholm, I., Clark, S. J., Clark, A. J., Stewart, C. N., et al. (2013). Tritrophic choice experiments with Bt plants, the diamondback moth (Plutella xylostella) and the parasitoid Cotesia plutellae. Trans. Res. 12, 351–361. doi: 10.1023/A:1023342027192

PubMed Abstract | CrossRef Full Text | Google Scholar

Sen, I., Oglakci, M., Bolek, Y., Cicek, B., Kisakurek, N., and Aydin, S. (2004). Assessing the Out-crossing ratio, isolation distance and pollinator insects in cotton (Gossypium hirsutum L.). Asian J. Plant Sci. 3, 724–727. doi: 10.3923/ajps.2004.724.727

CrossRef Full Text | Google Scholar

Shan, W., Jiang, Y., Han, J., and Wang, K. (2016). Comprehensive cytological characterization of the Gossypium hirsutum genome based on the development of a set of chromosome cytological markers. Crop J. 4, 256–265. doi: 10.1016/j.cj.2016.04.001

CrossRef Full Text | Google Scholar

Shaw, D. R., Givens, W. A., Farno, L. A., Gerard, P. D., Jordan, D., Johnson, W. G., et al. (2009). Using a grower survey to assess the benefits and challenges of glyphosate-resistant cropping systems for weed management in U.S. corn, cotton, and soybean. Weed Technol. 23, 134–149. doi: 10.1614/WT-08-042.1

CrossRef Full Text | Google Scholar

Skevas, T., Stefanou, S. E., and Lansink, A. O. (2013). Do farmers internalise environmental spillovers of pesticides in production? J. Agric. Econ. 64, 624–640. doi: 10.1111/1477-9552.12007

CrossRef Full Text | Google Scholar

Survey System (2018). Sample Size Calculator. Available online at: https://www.surveysystem.com/sscalc.htm (Accessed January 29, 2018).

Tabashnik, B. E., Brévault, T., and Carrière, Y. (2013). Insect resistance to Bt crops: lessons from the first billion acres. Nat. Biotechnol. 31, 510–521. doi: 10.1038/nbt.2597

PubMed Abstract | CrossRef Full Text | Google Scholar

Tabashnik, B. E., Gassmann, J. A., Crowder, W. D., and Carrière, Y. (2008). Insect resistance to Bt crops: evidence versus theory. Nat. Biotechnol. 26, 199–202. doi: 10.1038/nbt1382

PubMed Abstract | CrossRef Full Text | Google Scholar

Tamez, P. (2010). A review of U.S. and Mexican cooperation to develop insect resistance management and monitoring methods for surveying transgenic crops expressing Bacillus thuringiensis proteins: 2003 to 2010. Southwest. Entomol. 35, 373–382. doi: 10.3958/059.035.0318

CrossRef Full Text | Google Scholar

Terán-Vargas, A. P., Rodríguez, J. C., Blanco, C. A., Martínez-Carrillo, J. L., Cibrián-Tovar, J., Sánchez-Arroyo, H., et al. (2005). Bollgard cotton and resistance of tobacco budworm (Lepidoptera: Noctuidae) to conventional insecticides in southern Tamaulipas, Mexico. J. Econ. Entomol. 98, 2203–2209. doi: 10.1093/jee/98.6.2203

PubMed Abstract | CrossRef Full Text | Google Scholar

Theiling, K. M., and Croft, B. A. (1988). Pesticide side-effects on arthropod natural enemies: a database summary. Agr. Ecosys. Environ. 21, 191–218. doi: 10.1016/0167-8809(88)90088-6

CrossRef Full Text | Google Scholar

Tian, J. C., Yao, J., Long, L. P., Romeis, J., and Shelton, A. M. (2015). Bt crops benefit natural enemies to control non-target pests. Sci. Rep. 5:16636. doi: 10.1038/srep16636

PubMed Abstract | CrossRef Full Text | Google Scholar

Umbeck, P. F., Barton, K. A., Nordheim, E. V., McCarty, J. C., Parrot, W. L., and Jenkins, J. N. (1991). Degree of pollen dispersal by insects from a field test of genetically engineered cotton. J. Econ. Entomol. 84, 1943–1950. doi: 10.1093/jee/84.6.1943

CrossRef Full Text | Google Scholar

Van Deynze, A., Sundstrom, F., and Bradford, K. (2005). Pollen-mediated gene flow in California cotton depends on pollinator activity. Crop Sci. 45, 1565–1570. doi: 10.2135/cropsci2004.0463

CrossRef Full Text | Google Scholar

Wang, S., Just, D. R., and Pinstrup-Andersen, P. (2006). Damage from secondary pests and the need for refuge in China. Nat. Res. Manag. Pol. 30, 625–637. doi: 10.1007/978-0-387-36953-2_28

CrossRef Full Text | Google Scholar

Wang, S., Just, D. R., and Pinstrup-Andersen, P. (2008). Bt- cotton and secondary pests. Int. J. Biotechnol. 10, 113–121. doi: 10.1504/IJBT.2008.018348

CrossRef Full Text | Google Scholar

Wang, Z. J., Lin, H., Huang, J. K., Hu, R., Rozelle, S., and Pray, C. (2009). Bt cotton in China: are secondary insect infestations offsetting the benefits in farmer fields? Agr. Sci. China 8, 83–90. doi: 10.1016/S1671-2927(09)60012-2

CrossRef Full Text | Google Scholar

Wegier, A., Piñeyro-Nelson, A., Alarcón, J., Gálvez-Mariscal, A., Álvarez-Buylla, E. R., and Piñero, D. (2011). Recent long-distance transgene flow into wild populations conforms to historical patterns of gene flow in cotton (Gossypium hirsutum) at its centre of origin. Mol. Ecol. 20, 4182–4194. doi: 10.1111/j.1365-294X.2011.05258.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilson, D. F., Flint, H. M., Deaton, R. W., Fischhoff, D. A., Perlak, F. J., Armstrong, T. A., et al. (1992). Resistance of cotton lines containing a Bacillus thuringiensis toxin to pink bollworm (Lepidoptera: Gelechiidae) and other insects. J. Econ. Entomol. 85, 1516–1521. doi: 10.1093/jee/85.4.1516

CrossRef Full Text | Google Scholar

Yan, S., Zhu, J., Zhu, W., Li, Z., Shelton, A. M., Luo, J., et al. (2015). Pollen-mediated gene flow from transgenic cotton under greenhouse conditions is dependent on different pollinators. Sci. Rep. 5:15917. doi: 10.1038/srep15917

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, B. H., Pan, X. P., Guo, T. L., Wang, Q. L., and Anderson, T. A. (2005). Measuring gene flow in the cultivation of transgenic cotton (Gossypium hirsutum L.). Mol. Biotechnol. 31, 11–20. doi: 10.1385/MB:31:1:011

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, J. H., Ho, P., and Azadi, H. (2011). Benefits of Bt cotton counterbalanced by secondary pests? Perceptions of ecological change in China. Environ. Monitor. Assess. 173, 985–994. doi: 10.1007/s10661-010-1439-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Zwahlen, C., Nentwig, W., Bigler, F., and Hilbeck, A. (2003). Tritrophic interactions of transgenic Bacillus thuringiensis corn, Anaphothrips obscurus (Thysanoptera: Thripidae), and the predator Orius majusculus (Heteroptera: Anthocoridae). Environ. Entomol. 29, 846–850. doi: 10.1603/0046-225X-29.4.846

CrossRef Full Text | Google Scholar