Limited fitness costs of herbicide-resistance traits in Amaranthus tuberculatus facilitate resistance evolution

Monitoring the temporal changes in herbicide-resistant Amaranthus tuberculatus: a landscape-scale probability-based estimation in Iowa

BACKGROUND

A landscape-scale probability-based sampling of Iowa soybean [Glycine max (L.) Merr.] fields was conducted in 2013 and 2019; Amaranthus tuberculatus [Moq.] J.D. Sauer seed was collected from 97 random geospatial selected fields. The objectives were to evaluate the prevalence and distribution of herbicide-resistant A. tuberculatus (waterhemp) in soybean fields and evaluate temporal changes over 6 years. Amaranthus tuberculatus seedlings were evaluated for resistance to imazethapyr, atrazine, glyphosate, lactofen and mesotrione at 1× and 4× label rates.

RESULTS

Resistance to imazethapyr, glyphosate, lactofen and mesotrione at the 1× rate increased significantly from 2013 to 2019 and was found in 99%, 97%, 16% and 15% of Iowa A. tuberculatus populations in 2019, respectively. Resistance to atrazine at the 4× rate increased over time; atrazine resistance was found in 68% of populations in 2019. Three-way multiple herbicide-resistant A. tuberculatus was the most frequent and increased significantly to 4× rates from 16% in 2013 to 43% of populations in 2019. All A. tuberculatus populations resistant to HPPD-inhibitor herbicides also were resistant to atrazine.

CONCLUSION

To the best of our knowledge, this is the first probability-based study that presented evolution of A. tuberculatus herbicide resistance over time. The results demonstrated that imazethapyr, atrazine and glyphosate resistance in Iowa A. tuberculatus populations was frequent whereas resistance to lactofen and mesotrione was less frequent. Most Iowa A. tuberculatus populations evolved resistance to multiple sites of action over time. The results of our study are widely applicable given the similarities in weed management practices throughout the Midwest United States. © 2023 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Rapid weed adaptation and range expansion in response to agriculture over the past two centuries

North America has experienced a massive increase in cropland use since 1800, accompanied more recently by the intensification of agricultural practices. Through genome analysis of present-day and historical samples spanning environments over the past two centuries, we studied the effect of these changes in farming on the extent and tempo of evolution across the native range of the common waterhemp (Amaranthus tuberculatus), a now pervasive agricultural weed. Modern agriculture has imposed strengths of selection rarely observed in the wild, with notable shifts in allele frequency trajectories since agricultural intensification in the 1960s. An evolutionary response to this extreme selection was facilitated by a concurrent human-mediated range shift. By reshaping genome-wide diversity across the landscape, agriculture has driven the success of this weed in the 21st century.

Dissecting weed adaptation: Fitness and trait correlations in herbicide-resistant Alopecurus myosuroides

BACKGROUND

Unravelling the genetic architecture of non-target-site resistance (NTSR) traits in weed populations can inform questions about the inheritance, trade-offs and fitness costs associated with these traits. Classical quantitative genetics approaches allow study of the genetic architecture of polygenic traits even where the genetic basis of adaptation remains unknown. These approaches have the potential to overcome some of the limitations of previous studies into the genetics and fitness of NTSR.

RESULTS

Using a quantitative genetic analysis of 400 pedigreed Alopecurus myosuroides seed families from nine field-collected populations, we found strong heritability for resistance to the acetolactate synthase and acetyl CoA carboxylase inhibitors (h²  = 0.731 and 0.938, respectively), and evidence for shared additive genetic variance for resistance to these two different herbicide modes of action, r_g  = 0.34 (survival), 0.38 (biomass). We find no evidence for genetic correlations between life-history traits and herbicide resistance, indicating that resistance to these two modes of action is not associated with large fitness costs in blackgrass. We do, however, demonstrate that phenotypic variation in plant flowering characteristics is heritable, h²  = 0.213 (flower height), 0.529 (flower head number), 0.449 (time to flowering) and 0.372 (time to seed shed), demonstrating the potential for adaptation to other nonchemical management practices (e.g. mowing of flowering heads) now being adopted for blackgrass control.

CONCLUSION

These results highlight that quantitative genetics can provide important insight into the inheritance and genetic architecture of NTSR, and can be used alongside emerging molecular techniques to better understand the evolutionary and fitness landscape of herbicide resistance. © 2022 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Perspective: Scientific rigor or ritual? Statistical significance in pest management science

P values, binary hypothesis tests, and statistical significance are too often overused or used incorrectly in pest management reports. These statistical results are not the point of an analysis, they are just a by-product that can sometimes be informative. The point of the analysis is always a biological interpretation of the data focusing on whether something new has been learned and then quantifying what was learned in terms of biological effects or supported hypotheses. The term 'statistical significance' should be retired if only because of the common misunderstanding that it implies biological relevance; these terms are not equivalent. © 2021 Society of Chemical Industry.

Can new herbicide discovery allow weed management to outpace resistance evolution?

While herbicides are the most effective and widely adopted weed management approach, the evolution of multiple herbicide resistance in damaging weed species threatens the yield and profitability of many crops. Weeds accumulate multiple resistance mechanisms through sequential selection and/or gene flow, with long-range and international transport of herbicide-resistant weeds proving to be a serious issue. Metabolic resistance mechanisms can confer resistance across multiple sites of action and even to herbicides not yet discovered. When a new site of action herbicide is introduced to control a key driver weed, it likely will be one of very few effective available herbicide options for that weed in a specific crop due to the continuous use of herbicides over the years and the resulting accumulation of resistance mechanisms, placing it at even higher risk to be rapidly lost to resistance due to the high selection pressure it will experience. The number of available, effective herbicides for certain driver weeds is decreasing over time because the rate of resistance evolution is faster than the rate of new herbicide discovery. Effective monitoring for species movement and diagnostics for resistance should be deployed to rapidly identify emerging resistance to any new site of action. While innovation in herbicide discovery is urgently needed to combat the pressing issue of resistance in weeds, the rate of selection for herbicide resistance in weeds must be slowed through changes in the patterns of how herbicides are used. © 2021 Society of Chemical Industry. © 2021 Society of Chemical Industry.

Fitness Cost Associated With Enhanced EPSPS Gene Copy Number and Glyphosate Resistance in an Amaranthus tuberculatus Population

The evolution of resistance to pesticides in agricultural systems provides an opportunity to study the fitness costs and benefits of novel adaptive traits. Here, we studied a population of Amaranthus tuberculatus (common waterhemp), which has evolved resistance to glyphosate. The growth and fitness of seed families with contrasting levels of glyphosate resistance was assessed in the absence of glyphosate to determine their ability to compete for resources under intra- and interspecific competition. We identified a positive correlation between the level of glyphosate resistance and gene copy number for the 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) glyphosate target, thus identifying gene amplification as the mechanism of resistance within the population. Resistant A. tuberculatus plants were found to have a lower competitive response when compared to the susceptible phenotypes with 2.76 glyphosate resistant plants being required to have an equal competitive effect as a single susceptible plant. A growth trade-off was associated with the gene amplification mechanism under intra-phenotypic competition where 20 extra gene copies were associated with a 26.5 % reduction in dry biomass. Interestingly, this growth trade-off was mitigated when assessed under interspecific competition from maize.

Investigating the presence of compensatory evolution in dicamba resistant IAA16 mutated kochia (Bassia scoparia)†

BACKGROUND

Lack of fitness costs has been reported for multiple herbicide resistance traits, but the underlying evolutionary mechanisms are not well understood. Compensatory evolution that ameliorates resistance costs, has been documented in bacteria and insects but rarely studied in weeds. Dicamba resistant IAA16 (G73N) mutated kochia was previously found to have high fecundity in the absence of competition, regardless of significant vegetative growth defects. To understand if costs of dicamba resistance can be compensated through traits promoting reproductive success in kochia, we thoroughly characterized the reproductive growth and development of different G73N kochia biotypes. Flowering phenology, seed production and reproductive allocation were quantified through greenhouse studies, floral (stigma-anthers distance) and seed morphology, as well as resulting mating and seed dispersal systems were studied through time-course microcopy images.

RESULTS

G73N covaried with multiple phenological, morphological and ecological traits that improve reproductive fitness: (i) 16-60% higher reproductive allocation; (ii) longer reproduction phase through early flowering (2-7 days); (iii) smaller stigma-anthers separation (up to 60% reduction of herkogamy and dichogamy) that can potentially promote selfing and reproductive assurance; (iv) 'winged' seeds with 30-70% longer sepals that facilitate long-distance seed dispersal.

CONCLUSION

The current study demonstrates that costs of herbicide resistance can be ameliorated through coevolution of other fitness penalty alleviating traits. As illustrated in a hypothetical model, the evolution of herbicide resistance is an ongoing fitness maximization process, which poses challenges to contain the spread of resistance. © 2020 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

A dicamba resistance-endowing IAA16 mutation leads to significant vegetative growth defects and impaired competitiveness in kochia (Bassia scoparia)†

BACKGROUND

Precise quantification of the fitness cost of synthetic auxin resistance has been impeded by lack of knowledge about the genetic basis of resistance in weeds. Recent elucidation of a resistance-endowing IAA16 mutation (G73N) in the key weed species kochia (Bassia scoparia), allows detailed characterization of the contribution of resistance alleles to weed fitness, both in the presence and absence of herbicides. Different G73N genotypes from a segregating resistant parental line (9425) were characterized for cross-resistance to dicamba, 2,4-d and fluroxypyr, and changes on stem/leaf morphology and plant architecture. Plant competitiveness and dominance of the fitness effects was quantified through measuring biomass and seed production of three F₂ lines in two runs of glasshouse replacement series studies.

RESULTS

G73N confers robust resistance to dicamba but only moderate to weak resistance to 2,4-D and fluroxypyr. G73N mutant plants displayed significant vegetative growth defects: (i) they were 30-50% shorter, with a more tumbling style plant architecture, and (ii) they had thicker and more ovate (versus lanceolate and linear) leaf blades with lower photosynthesis efficiency, and 40-60% smaller stems with less-developed vascular bundle systems. F₂ mutant plants had impaired plant competitiveness, which can lead to 80-90% less biomass and seed production in the replacement series study. The pleiotropic effects of G73N were mostly semidominant (0.5) and fluctuated with the environments and traits measured.

CONCLUSION

G73N is associated with significant vegetative growth defects and reduced competitiveness in synthetic auxin-resistant kochia. Management practices should target resistant kochia's high vulnerability to competition in order to effectively contain the spread of resistance.

Herbicide resistance in Amaranthus tuberculatus†

Amaranthus tuberculatus is the major weed species in many midwestern US row-crop production fields, and it is among the most problematic weeds in the world in terms of its ability to evolve herbicide resistance. It has now evolved resistance to herbicides spanning seven unique sites of action, with populations and even individual plants often possessing resistance to several herbicides/herbicide groups. Historically, herbicide target-site changes accounted for most of the known resistance mechanisms in this weed; however, over the last few years, non-target-site mechanisms, particularly enhanced herbicide detoxification, have become extremely common in A. tuberculatus. Unravelling the genetics and molecular details of non-target-site resistance mechanisms, understanding the extent to which they confer cross resistance to other herbicides, and understanding how they evolve remain as critical research endeavors. Transcriptomic and genomics approaches are already facilitating such studies, the results of which hopefully will inform better resistance-mitigation strategies. The largely unprecedented level of herbicide resistance in A. tuberculatus is not only a fascinating example of evolution in action, but it is a serious and growing threat to the sustainability of midwestern US cropping systems. © 2020 Society of Chemical Industry.

Agricultural adaptation in the native North American weed waterhemp, Amaranthus tuberculatus (Amaranthaceae)

There is increasing interest in documenting adaptation of weedy plant species to agricultural ecosystems, beyond the evolution of herbicide resistance. Waterhemp (Amaranthus tuberculatus) is a native plant of the Midwestern U.S. that began infesting agricultural fields in the 20th century within the central portion of its range. We hypothesized that the vegetative growth and reproductive traits of waterhemp from this heavily infested central region provide differential fitness benefits in agricultural environments. We collected seeds from across the species' native range, representing regions with varying degrees of waterhemp infestation, and planted them together in common garden soybean plots. A 2010 common garden experiment was conducted within the range of agriculturally weedy waterhemp (in Missouri), and a 2011 common garden experiment was conducted outside of this range (in Ohio). Days to flowering and flowering plant height, mature plant size data (height, number of branches, and length of the longest branch), and above-ground biomass were measured to estimate relative fitness. In both common garden locations, plants from regions where waterhemp occurs as an agricultural weed - including those from the heavily infested Mississippi Valley region (Iowa, Illinois, and Missouri) and the less severely infested Plains region (Nebraska, Kansas, and Oklahoma) - had higher relative performance in almost all fitness-related measures than plants from the Northeast region (Ohio, Michigan, and Ontario), which had little to no agriculturally weedy waterhemp at the time of our study. Further analysis revealed that fewer days to flowering in the Northeast populations can be largely accounted for by latitude of origin, suggesting a strong genetic influence on this reproductive trait. These findings suggest intraspecific variation in agricultural adaptation in a native U.S. weed, and support the use of agricultural weeds to study adaptation.

A Trp574Leu Target-Site Mutation Confers Imazamox Resistance in Multiple Herbicide-Resistant Wild Poinsettia Populations from Brazil

Wild poinsettia (Euphorbia heterophylla L.) is an important weed species in southern Brazil, especially due to the evolution of multiple herbicide resistance (e.g., acetolactate synthase (ALS)- inhibitors, protoporphyrinogen oxidase inhibitors, and glyphosate). The mechanism of resistance to imazamox was investigated in two wild poinsettia populations (R1 and R2) from southern Brazil and compared to a known susceptible (S) population. Imazamox dose-response experiments revealed high levels of resistance: 45-fold and 224.5-fold based on dry biomass reduction, for R1 and R2, respectively. Extremely high concentrations of imazamox (20,000 µM) were not sufficient to provide 50% inhibition of ALS enzyme activity (I50) for R1 or R2. Hence, resistance levels were estimated to be greater than 123-fold for both populations based on in vitro ALS assays. The ALS gene from all R1 and R2 plants had a Trp574Leu mutation. A genotyping assay was developed to discriminate resistant and susceptible alleles based on the Trp574Leu mutation.

Herbicide resistant weeds: A call to integrate conventional agricultural practices, molecular biology knowledge and new technologies

Herbicide resistant (HR) weeds are of major concern in modern agriculture. This situation is exacerbated by the massive adoption of herbicide-based technologies along with the overuse of a few active ingredients to control weeds over vast areas year after year. Also, many other anthropological, biological, and environmental factors have defined a higher rate of herbicide resistance evolution in numerous weed species around the world. This review focuses on two central points: 1) how these factors have affected the resistance evolution process; and 2) which cultural practices and new approaches would help to achieve an effective integrated weed management. We claim that global climate change is an unnoticed factor that may be acting on the selection of HR weeds, especially those evolving into non-target-site resistance mechanisms. And we present several new tools -such as Gene Drive and RNAi technologies- that may be adopted to cope with herbicide resistance spread, as well as discuss their potential application at field level. This is the first review that integrates agronomic and molecular knowledge of herbicide resistance. It covers not only the genetic basis of the most relevant resistance mechanisms but also the strengths and weaknesses of traditional and forthcoming agricultural practices.

No evidence for early fitness penalty in glyphosate-resistant biotypes of Conyza canadensis: Common garden experiments in the absence of glyphosate

Strong selection from herbicides has led to the rapid evolution of herbicide-resistant weeds, greatly complicating weed management efforts worldwide. In particular, overreliance on glyphosate, the active ingredient in RoundUp®, has spurred the evolution of resistance to this herbicide in ≥40 species. Previously, we reported that Conyza canadensis (horseweed) has evolved extreme resistance to glyphosate, surviving at 40× the original 1× effective dosage. Here, we tested for underlying fitness effects of glyphosate resistance to better understand whether resistance could persist indefinitely in this self-pollinating, annual weed. We sampled seeds from a single maternal plant ("biotype") at each of 26 horseweed populations in Iowa, representing nine susceptible biotypes (S), eight with low-level resistance (LR), and nine with extreme resistance (ER). In 2016 and 2017, we compared early growth rates and bolting dates of these biotypes in common garden experiments at two sites near Ames, Iowa. Nested ANOVAs showed that, as a group, ER biotypes attained similar or larger rosette size after 6 weeks compared to S or LR biotypes, which were similar to each other in size. Also, ER biotypes bolted 1-2 weeks earlier than S or LR biotypes. These fitness-related traits also varied among biotypes within the same resistance category, and time to bolting was inversely correlated with rosette size across all biotypes. Disease symptoms affected 40% of all plants in 2016 and 78% in 2017, so we did not attempt to measure lifetime fecundity. In both years, the frequency of disease symptoms was greatest in S biotypes and similar in LR versus ER biotypes. Overall, our findings indicate there are no early growth penalty and possibly no lifetime fitness penalty associated with glyphosate resistance, including extremely strong resistance. We conclude that glyphosate resistance is likely to persist in horseweed populations, with or without continued selection pressure from exposure to glyphosate.

Evolutionary and ecological insights from herbicide-resistant weeds: what have we learned about plant adaptation, and what is left to uncover?

The evolution of herbicide resistance in crop weeds presents one of the greatest challenges to agriculture and the production of food. Herbicide resistance has been studied for more than 60 yr, in the large part by researchers seeking to design effective weed control programs. As an outcome of this work, various unique questions in plant adaptation have been addressed. Here, I collate recent research on the herbicide-resistant problem in light of key questions and themes in evolution and ecology. I highlight discoveries made on herbicide-resistant weeds in three broad areas - the genetic basis of adaptation, evolutionary constraints, experimental evolution - and similarly discuss questions left to be answered. I then develop how one would use herbicide-resistance evolution as a model for studying eco-evolutionary dynamics within a community context. My overall goals are to highlight important findings in the weed science literature that are relevant to themes in plant adaptation and to stimulate the use of herbicide-resistant plants as models for addressing key questions within ecology and evolution.

Do plants pay a fitness cost to be resistant to glyphosate?

We reviewed the literature to understand the effects of glyphosate resistance on plant fitness at the molecular, biochemical and physiological levels. A number of correlations between enzyme characteristics and glyphosate resistance imply the existence of a plant fitness cost associated with resistance-conferring mutations in the glyphosate target enzyme, 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS). These biochemical changes result in a tradeoff between the glyphosate resistance of the EPSPS enzyme and its catalytic activity. Mutations that endow the highest resistance are more likely to decrease catalytic activity by reducing the affinity of EPSPS for its natural substrate, and/or slowing the velocity of the enzyme reaction, and are thus very likely to endow a substantial plant fitness cost. Prediction of fitness costs associated with EPSPS gene amplification and overexpression can be more problematic. The validity of cost prediction based on the theory of evolution of gene expression and resource allocation has been cast into doubt by contradictory experimental evidence. Further research providing insights into the role of the EPSPS cassette in weed adaptation, and estimations of the energy budget involved in EPSPS amplification and overexpression are required to understand and predict the biochemical and physiological bases of the fitness cost of glyphosate resistance.

Molecular and physiological characterization of six-way resistance in an Amaranthus tuberculatus var. rudis biotype from Missouri

BACKGROUND

Previous research reported the first case of six-way herbicide resistance in a common waterhemp (Amaranthus tuberculatus var. rudis) biotype from Missouri, USA designated MO-Ren. This study investigated the mechanisms of multiple-resistance in the MO-Ren biotype to herbicides from six site-of-action (SOA) groups, i.e. synthetic auxins, 5-enolypyruvyl-shikimate-3-phosphate synthase (EPSPS)-, protoporphyrinogen oxidase (PPO)-, acetolactate synthase (ALS)-, photosystem II (PSII)-, and 4-hydroxyphenyl-pyruvate-dioxygenase (HPPD)-inhibitors.

RESULTS

Genomic DNA sequencing confirmed the presence of known mutations associated with ALS- or PPO-inhibiting herbicide resistance: the Trp-574-Leu amino acid substitution in the ALS enzyme and the codon deletion corresponding to the ΔG210 in the PPX2 enzyme. No target-site point mutations associated with resistance to PSII- and EPSPS-inhibitors were detected. Quantitative polymerase chain reaction (qPCR) indicated that MO-Ren plants contained five-fold more copies of the EPSPS gene than susceptible plants. Malathion in combination with 2,4-D (2,4-dichlorophenoxyacetic acid), mesotrione, and chlorimuron POST enhanced the activity of these herbicides indicating that metabolism due to cytochrome P450 monooxygenase activity was involved in herbicide resistance. 4-Chloro-7-nitrobenzofurazan (NBD-Cl), a glutathione-S-transferase (GST)-inhibitor, in combination with atrazine did not reduce the biomass accumulation. Reduced absorption or translocation of 2,4-D did not contribute to resistance. However, the resistant biotype metabolized 2,4-D, seven- to nine-fold faster than the susceptible.

CONCLUSION

Target-site point mutations, gene amplification, and elevated rates of metabolism contribute to six-way resistance in the MO-Ren biotype, suggesting both target site and non-target site mechanisms contribute to multiple herbicide resistance in this Amaranthus tuberculatus biotype. © 2018 Society of Chemical Industry.

An ecological perspective on managing weeds during the great selection for herbicide resistance

More than 70 years after modern agriculture declared a 'war on weeds', they continue to thrive and suppress crop yields. Viewing weeds as an enemy that can be defeated if only a powerful enough technology can be deployed is a losing proposition. The latest evidence for the inadequacy of this approach, rampant evolution of multiple herbicide-resistant weed genotypes and dwindling options for chemical control in many production systems, should be seen as an urgent message to all those concerned with the science and practice of weed management: we need another way of thinking about the weed resistance issue. Fortunately, the theoretical and practical foundation of this alternative approach, ecological weed management, has been in development for decades. Here, we use Barry Commoner's laws of ecology as a conceptual frame for a review of some of the fundamental concepts of ecological weed management. © 2018 Society of Chemical Industry.

Fenoxaprop-P-ethyl and mesosulfuron-methyl resistance status of shortawn foxtail (Alopecurus aequalis Sobol.) in eastern China

Resistance to the acetyl-coenzyme A carboxylase (ACCase)- and acetolactate synthase (ALS)- inhibiting herbicides in shortawn foxtail (Alopecurus aequalis) has been reported in wheat fields of eastern China. To better understand the distribution of the resistant populations and the occurrence of the target-site mutations, 74 populations collected from Anhui, Jiangsu and Shandong province were surveyed, and the ACCase and ALS gene fragments, encompassing all the documented mutant codon positions, were amplified and sequenced. Plants from 37 and 34 populations survived fenoxaprop-P-ethyl and mesosulfuron-methyl treatment at 62.1 g a.i. ha^-1 and 9 g a.i. ha^-1 respectively, with different survival rates. Twenty-seven populations exhibited multiple resistance to fenoxaprop-P-ethyl and mesosulfuron-methyl. Whole-plant dose-response experiments showed that the resistance index ranged from 6.2 to 167.8 for fenoxaprop-P-ethyl, and from 7.8 to 139.5 for mesosulfuron-methyl. Four ACCase (I1781L, I2041N, I2041T and D2078G) and four ALS (P197R, P197S, P197T and W574 L) resistance mutations were detected respectively. Individuals containing two amino acid substitutions were also found. D2078G and W574 L were predominant ACCase and ALS gene mutations respectively. This study has shown that fenoxaprop-P-ethyl and mesosulfuron-methyl resistance was prevalent in A. aequalis in eastern China, and target site mutations in the ACCase and ALS gene were one of the most common mechanisms.

Transgene escape and persistence in an agroecosystem: the case of glyphosate-resistant Brassica rapa L. in central Argentina

Brassica rapa L. is an annual Brassicaceae species cultivated for oil and food production, whose wild form is a weed of crops worldwide. In temperate regions of South America and especially in the Argentine Pampas region, this species is widely distributed. During 2014, wild B. rapa populations that escaped control with glyphosate applications by farmers were found in this area. These plants were characterized by morphology and seed acidic profile, and all the characters agreed with B. rapa description. The dose-response assays showed that the biotypes were highly resistant to glyphosate. It was also shown that they had multiple resistance to AHAS-inhibiting herbicides. The transgenic origin of the glyphosate resistance in B. rapa biotypes was verified by an immunological test which confirmed the presence of the CP4 EPSPS protein and by an event-specific GT73 molecular marker. The persistence of the transgene in nature was confirmed for at least 4 years, in ruderal and agrestal habitats. This finding suggests that glyphosate resistance might come from GM oilseed rape crops illegally cultivated in the country or as a seed contaminant, and it implies gene flow and introgression between feral populations of GM B. napus and wild B. rapa. The persistence and spread of the resistance in agricultural environments was promoted by the high selection pressure imposed by intensive herbicide usage in the prevalent no-till farming systems.