Glyphosate resistance: state of knowledge.

Studies of mechanisms of resistance to glyphosate have increased current understanding of herbicide resistance mechanisms. Thus far, single-codon non-synonymous mutations of EPSPS (5-enolypyruvylshikimate-3-phosphate synthase) have been rare and, relative to other herbicide mode of action target-site mutations, unconventionally weak in magnitude for resistance to glyphosate. However, it is possible that weeds will emerge with non-synonymous mutations of two codons of EPSPS to produce an enzyme endowing greater resistance to glyphosate. Today, target-gene duplication is a common glyphosate resistance mechanism and could become a fundamental process for developing any resistance trait. Based on competition and substrate selectivity studies in several species, rapid vacuole sequestration of glyphosate occurs via a transporter mechanism. Conversely, as the chloroplast requires transporters for uptake of important metabolites, transporters associated with the two plastid membranes may separately, or together, successfully block glyphosate delivery. A model based on finite glyphosate dose and limiting time required for chloroplast loading sets the stage for understanding how uniquely different mechanisms can contribute to overall glyphosate resistance.

1 INTRODUCTION

Glyphosate stands alone in many categories.1,2 This N-phosphonomethyl-modified derivative of glycine is able to form a stabile dead-end complex with EPSPS (5-enolypyruvylshikimate-3-phosphate synthase, EC 2.5.1.19)3 as a transition-state analogue competitive with phosphoenolpyruvate (PEP).4 The inhibition of the shikimate pathway located in the chloroplast was attributed to glyphosate's mode of action (MOA) 8 years after commercialization.5,6 Some thought resistance to glyphosate was unlikely to evolve owing to its transition-state mimicry, lack of plant metabolism and absence of known active transporters for herbicides in plants,7 and, indeed, during the first 15 years of glyphosate use, primarily as a non-selective treatment before planting, glyphosate-resistant (GR) weeds were not found (1972–1997). However, in the last 15 years (1998–2013), glyphosate resistance in 24 species on six continents has been documented.8,9

The basic herbicide resistance mechanisms can be catalogued as follows: (1) target-site resistance, typically represented by amino acid substitutions that affect herbicide interactions at the target enzyme; (2) metabolism, a chemical modification of the herbicide either by conjugation or degradation; (3) exclusion of the herbicide from the target, either physically with enhanced cuticular and other structural barriers or physiologically with active transporters.10 A fourth mechanism could be added – avoidance, pertaining to biochemical ability to handle the toxic agent produced by the pesticide and thereby avoid a toxic result. An example of avoidance in herbicides is the overexpression of superoxide dismutase as a rescue for paraquat.11 Intriguingly, the known glyphosate resistance mechanisms now exceed those described for any other herbicide. Mechanisms include target-site mutation, target-site gene duplication, active vacuole sequestration, limited cellular uptake and a rapid necrosis response. Notably, glyphosate resistance is generally lower in magnitude and less frequently involves target-site and metabolic mechanisms than in the case of herbicides that inhibit acetolactate synthase (ALS) or acetyl-CoA carboxylase (ACCase).8

2 EXCLUSION MECHANISMS

The first investigations of glyphosate resistance in Lolium rigidum found no biochemical differences between GR and glyphosate-susceptible (GS) populations, indicating that quite subtle differences in physiology could result in glyphosate resistance.12,13 Early 14C glyphosate translocation studies revealed restricted translocation in GR L. rigidum.14 Similarly, 14C glyphosate showed restricted systemic delivery in Conyza canadensis when comparing GR and GS populations.15 The curiosity was that target-site resistance (a Pro106Ser substitution, see Section 4) afforded approximately twofold resistance in Eleusine indica,16 whereas the very slight differences (<2×) in restricted translocation in C. canadensis and L. rigidum afforded resistance magnitudes 6–8× the recommended field rates (1× = 860 g ha−1). Something substantial was happening to the glyphosate beyond the usual target-site and metabolic resistance mechanisms observed for other herbicides. The observations of restricted translocation quickly highlighted a resistance mechanism that became widely used in discussions of glyphosate resistance. But really, in the words of a famous bard, it is necessary to ‘Find out the cause of this effect, Or rather say, the cause of this defect, For this effect defective comes by cause.’ (Hamlet 2:2, Shakespeare).

2.1 Vacuolar sequestration

Experiments using 31P nuclear magnetic resonance (NMR) to observe glyphosate in vivo entering cells and cellular compartments have identified at least two new resistance mechanisms that result in restricted translocation correlated with the glyphosate resistance observed.17 The observation18 of rapid sequestration of glyphosate in C. canadensis by 31P NMR demonstrated the application of an old technique19–21 to the biochemistry of resistance in weedy species. The pH versus chemical shift (31P NMR signal) changes of the glyphosate phosphonate have been well documented.22 Hence, after dosing a plant with glyphosate prior to or during NMR observation, the pH changes make it possible to follow the passage of glyphosate from the apoplast buffer (pH 5) into the cell cytoplasm (pH 6.8) and further into the vacuole (pH 5.5). The overlay of the apoplast glyphosate chemical shift in the spectrum with the vacuole glyphosate signal was resolved using a pulse chase method, by first using the leaf-sample circulating stabilization buffer to challenge with glyphosate. Spectra were collected while observing initial cellular glyphosate uptake for 10 h, and then removing glyphosate from the circulating buffer solution to eliminate the apoplast glyphosate signal and reveal the underlying glyphosate vacuole signal. The difference in glyphosate vacuole sequestration rates between GR and GS C. canadensis was about 10×. The range of concentrations for the glyphosate challenge allowed the rates of cellular uptake and even vacuole sequestration to be measured.17 These NMR results on the Delaware GR C. canadensis population23 showed that the resistant trait is likely a singular biochemical activity (sequestration) and corroborated inheritance data showing a nuclear inherited single gene.24

To characterize the specificity of glyphosate uptake, other solutes were added such as the glyphosate degradation product aminomethyl phosphonate (AMPA) and the structurally similar glycine. Both showed significant competitive inhibition of glyphosate uptake. This competitive nature of the glyphosate uptake is an expected property of an active transport system.17 Varying the glyphosate challenge concentration revealed cytoplasm concentrations accumulating with a linear dependence to >20 mM glyphosate, indicating that the cellular plasmalemma uptake transporter could not be saturated. During the chase period, the glyphosate vacuole accumulation observations indicated the rate of vacuole uptake. This vacuole uptake rate did saturate with increasing cytoplasmic glyphosate concentration. (Sammons RD, personal communication) Substrate saturation behavior (Km, Vmax) is a key characteristic of active transporters. Finally, another major observation in the pulse chase is that glyphosate does not leach out of the cytoplasm, nor out of the vacuole to the cytoplasm, unless there is evidence of cell death. That is, the transport of glyphosate across membranes is observed to be unidirectional in all species studied to date (23 in total).17 These observations also support the hypothesis that glyphosate cellular uptake and vacuole sequestration are mediated by active transport. Previously, glyphosate transport systems have been characterized insufficiently to explain observed uptake and translocation values.25,26 Given that these transporters are present, the realization that glyphosate may require assistance to cross any active bilayer further suggests that many transporters are working to allow glyphosate cell uptake, whole-plant translocation and subcellular distribution.

2.1.1 Temperature effects on vacuolar sequestration

C. canadensis survives freezing very well, which is consistent with its biennial lifecycle,27 making it a good candidate for low-temperature 31P NMR to determine how cold temperature affects the sequestration process. Pulse chase 31P NMR and whole-plant treatments of glyphosate in the cold (8 °C) demonstrated that vacuole sequestration was inhibited and that GR C. canadensis individuals were susceptible to glyphosate at cold temperatures.28 This provided an unequivocal demonstration that the resistance was due solely to the sequestration process and showed that the physiological properties of the resistance mechanism could be utilized to reverse the resistance to the herbicide. Interestingly, vacuolar sequestration has also been hypothesized as a mechanism of paraquat resistance,29 with a similar temperature dependence for paraquat resistance reported.30 Lipid bilayers are well known to respond to temperature differently, based on their lipid saturation levels, as these affect the membrane fluidity.31

The new exclusion mechanism of rapid vacuole sequestration was also demonstrated to be a major part of the glyphosate resistance in Lolium spp., with the sequestration rate correlating with the LD50.32 Glyphosate-resistant Lolium spp. from three continents were shown to sequester glyphosate in the vacuole, suggesting that the GR trait is present within diverse Lolium spp. and can be selected independently wherever Lolium is found. Glyphosate-resistant Sorghum halepense and Lolium spp. populations having reduced translocation were also made sensitive at suboptimal temperatures, further confirming that populations dependent on membrane mechanisms could be made sensitive in the cold.33 A comparison of GR L. multiflorum populations from Italy34 in cold acclimation and glyphosate treatments showed that the most resistant population, also shown to contain a EPSPS Pro106Ser mutation (see Section 4), was not made as sensitive by the cold acclimation. This supports the idea that multiple resistance mechanisms can contribute to the overall resistance magnitude, and not all mechanisms are temperature sensitive.

2.1.2 Variation among species in sensitivity to glyphosate

The demonstration of rapid vacuole sequestration then seemed to be the most reasonable explanation for the variation in species sensitivity to glyphosate. A hypothesis was developed that those species not using the vacuole would be more sensitive and those utilizing the vacuole to sequester glyphosate would be more tolerant and potential candidates for selection of resistance. Consequently, a survey of 23 species was carried out to compare glyphosate-sensitive [e.g. barley (Hordeum vulgare)] and hard-to-control weeds (e.g. Ipomoea lacunosa, Desmanthus illinoensis).17 The plants were all track sprayed with a 4× label rate of glyphosate and then examined 24 h later. No correlation of vacuole sequestration and glyphosate sensitivity was found, and in fact several different categories could be defined. Quite unexpectedly, it was observed that some species accumulated very little glyphosate in the cytoplasm. In one case, barley allowed very little glyphosate inside the cell, but a strong shikimate-3-phosphate (S3P) signal suggested a high degree of toxicity with the small amount of glyphosate, which is consistent with its sensitive categorization. However, Convolvulus arvensis also took up very little glyphosate, and there was no evidence of an S3P peak, which again is consistent with its hard-to-control classification. In another situation, Chenopodium album and I. lacunosa took up ∼10× more glyphosate than C. arvensis and did not sequester glyphosate in the vacuole, and yet no S3P was observed; therefore, glyphosate appeared to be less toxic, even though much more was entering the cytoplasm. Echinochloa colona uptake was similar to that of I. lacunosa, but 83% of the glyphosate was moved to the vacuole (the same as GR Conyza canadensis). Overall, 11 species took up very little glyphosate, including the following: Ambrosia artemisiifolia, E. colona, C. album, I. lacunosa, Xanthium strumarium, Camelina microcarpa, Setaria viridis, H. vulgare, Avena fatua, E. crus-galli and C. arvensis. Only E. colona sequestered glyphosate in the vacuole together with minimal uptake. A caveat would be that this lack of vacuole signal in these minimal uptake species was generally due to such low cytoplasm concentrations. The 13 species that took up moderate amounts of glyphosate included nine that used the vacuole, including Senna obtusifolia, Sida spinosa, Polygonum convolvulus, Ambrosia trifida, Sorghum bicolor, Urochloa platyphylla, L. rigidum and C. canadensis. P. convolvulus rivaled GR C. canadensis in its efficiency to sequester glyphosate.

For the remaining four species that did not use the vacuole but had good cellular uptake, GR populations of three species were able to restrict cellular uptake when compared with GS populations, including Amaranthus palmeri, A. tuberculatus and S. halepense. One particular case demonstrates this observation. In a survey of GR A. palmeri35 for various EPSPS copy numbers, two female plants (R20 and R25) were found to have very similar genomic EPSPS copy numbers (∼35 by qRT-PCR, see Section 3). Upon vegetative cloning of R20 and R25 and glyphosate titration (five replicates, WeatherMax®, 1/2× to 16× by factors of 2), the R20 plants survived 16× the field rate (13.8 kg ha−1), whereas the R25 plants only survived 4× the rate. Owing to this surprising result, both plants were then vegetatively cloned and examined by 31P NMR pulse chase with a 10 mM glyphosate challenge (according to the methods described).17 The R25 cloned individuals had a glyphosate cellular uptake rate 4× faster than the R20 cloned individuals, so that after 10 h the ratio of glyphosate/internal standard in R25 was 0.4 and for R20 the ratio was 0.1. Therefore, the reduced cellular uptake of glyphosate in the R20 individuals was shown to contribute >4× the label rate in resistance magnitude. These data showing limitations to glyphosate cellular uptake support a resistance mechanism of restricted cytoplasm entry. The 31P NMR data also indicated that glyphosate is not pumped out of the cell after the challenge is removed, so the limitation to glyphosate uptake is due to restricted entry rather than a shifted equilibrium due to the increased activity of an export transporter. Previously, Hetherington et al.36 studied 14C glyphosate efflux in tissue culture using transgenic GR soybean (Glycine max) and described slow rates of efflux correlating with cytoplasmic and vacuole pools measurable over long times because GR plant tissue was used. The time course of A. palmeri NMR experiments was not long enough to measure such slow efflux rates. Export transporters are very commonly found in microbial systems, and so multidrug resistance pumps would have been a reasonable hypothesis.37 Similar restricted entry results were observed for GR A. tuberculatus (Sammons RD, unpublished) and for GR S. halepense from Arkansas.17 Restricted entry into the source-leaf cells has been proposed in studies examining L. rigidum with 14C glyphosate translocation and offered as the origin of the restricted translocation observed for glyphosate systemic movement in the resistant populations.14 The translocation of 14C glyphosate in S. halepense was also observed to be restricted,38,39 and so a similar explanation can be offered. Specifically, translocation of glyphosate requires source-leaf cell uptake to access the phloem for systemic translocation. The restricted cellular entry of the R20 A. palmeri line and the GR A. tuberculatus described above does not result in the reduced translocation phenotype (Sammons RD, unpublished).

In this last group of four species studied in the survey, D. illinoensis was very unique, taking up 100× more glyphosate than C. arvensis, not using the vacuole and showing no S3P signal. This latter case prompts a new hypothesis whereby glyphosate has restricted entry into the chloroplast, and this last physiological barrier could well be the most potent of all glyphosate resistance mechanisms. This hypothesis allows explanations for the reduced toxicity of glyphosate even while present in the cell.

2.1.3 Molecular characterization of vacuolar sequestration

Some molecular evidence suggests that proteins associated with the vacuolar membrane may have roles in the glyphosate vacuolar sequestration resistance mechanism in C. canadensis, including ATP-binding cassette (ABC) transporters and a tonoplast intrinsic protein (TIP).40–42 Using a transcriptome profiling technique called iGentifier, a TIP was identified that was generally higher in GR than GS individuals from several populations and in GR individuals from segregating F2 populations.42 The TIP was also highly expressed in one GS population. A limitation of the iGentifier approach was that only 72 sequence tags were linked to previously obtained EST sequences, and only 52 sequence tags were able to be annotated. Using an Arabidopsis heterologous microarray for hybridization of untreated and glyphosate-treated RNA samples from one GR population, several genes were induced by glyphosate treatment, including four ABC transporters.42 To date, genes identified from transcriptome analysis have not been shown to confer glyphosate resistance. Genetic marker evidence (SSRs) suggested that resistance had evolved several times independently in C. canadensis, and clustering of GR populations based on iGentifier expression patterns suggested that similar mechanisms had evolved in several GR populations.42

The successful application of next-generation sequencing technology to GR C. canadensis identified ABC-transporter unigenes M10 and M11.41 M11 is an ortholog of Arabidopsis gene At3g13080 (an ATP-dependent multidrug-resistant protein-like transporter able to transport glutathione conjugates), which was also identified in the microarray study with approximately 30-fold increased expression following glyphosate treatment in a TN GR C. canadensis population.42 The M11 Arabidopsis ortholog is associated with the vacuolar membrane, making this gene an especially interesting candidate given the vacuole sequestration mechanism hypothesis.18 Using qRT-PCR to measure the expression relative to an actin gene, both M10 and M11 were constitutively more expressed in GR C. canadensis than in GS; both were upregulated in response to glyphosate treatment, and expression in GR remained higher than in GS.41 Additional ABC-transporter unigenes M6, M7 and P3 were also constitutively higher in GR than in GS, were induced in GS by glyphosate treatment and were unaffected by glyphosate treatment in GR.

Additional confirmatory results suggest a possible role for M10 and M11 in C. canadensis glyphosate resistance, when gene expression was compared using qRT-PCR in two GR and one GS Mediterranean C. canadensis populations.40 Constitutive expression of both M10 and M11 was higher in the GS population than in both GR populations. However, both M10 and M11 were highly (11–23-fold) induced in both GR populations by glyphosate treatment, and significant induction did not occur in the GS population.

Based on results available to date, it is not clear whether M10, M11 or other ABC transporters have a direct role in the C. canadensis vacuolar sequestration glyphosate resistance mechanism. Observed ABC-transporter expression has been generally constitutively higher in GR populations, and often induced to higher levels by glyphosate treatment.40–42 Additional research is necessary to establish the possible role(s) of ABC transporters in glyphosate resistance, including M10 and M11.

2.2 Other reduced translocation mechanisms

Some GR A. trifida populations have shown an interesting response to glyphosate application. Within 12 h after treatment, mature leaves show rapid necrosis and cell death.43 Glyphosate translocation is substantially reduced in GR individuals, presumably owing to the rapid necrosis and the inability of glyphosate to be exported from mature leaves. The molecular and genetic basis of the rapid necrosis response remains unknown, but future research will very likely provide important insights into cell signaling responses that occur after glyphosate treatment.

3 GENE DUPLICATION

Gene duplication is defined as the heritable replication of a DNA segment, resulting in one or more additional gene copies within the genome of an organism.44 Gene amplification is sometimes used synonymously with gene duplication, but gene amplification can also have a contrasting definition, such as the non-heritable replication of a DNA segment.44 In other fields, gene amplification is defined as the heritable replication of multiple copies of a DNA segment, while gene duplication indicates the replication of only a single copy.45 Gene duplication is vital for generating genomic diversity and is a common process in the evolutionary history of plants.46 Duplication can occur through unequal crossing-over events, chromosome duplication or transposable element activity.47 When a duplicated gene sequence contains transcription signals including promoters, or when a duplicated gene is inserted behind a promoter, the immediate results are increased expression of mRNA and protein. The extra gene copies can accumulate mutations without potentially negative consequences (e.g. loss of function), thereby increasing opportunities to select more advantageous mutations. Given the fitness advantages of improved mutations, duplicated gene copies can be deleted over time while retaining the improved gene copies, an important evolutionary process. Of particular relevance when considering glyphosate resistance mechanisms, examples have been reported of positive selection for gene dosage, leading to adaptive gene duplication in response to stressful environmental conditions.48

3.1 EPSPS gene duplication

Cell culture selection studies using glyphosate provided the first examples of EPSPS gene duplication in selected GR cell lines.49–53 In several cases, duplicated EPSPS copies were not maintained when glyphosate selection was removed, and not inherited when plantlets were cultured from selected cells, suggesting that the duplicated EPSPS genes were unstable or carried a fitness penalty.54

The first reported occurrence of EPSPS gene duplication in a field-evolved weed species was a GR A. palmeri population from Georgia, United States.55 Both EPSPS and ALS were measured using quantitative real-time PCR on genomic DNA, and EPSPS relative genomic copy number was calculated using ALS as an internal standard. Using this technique, the EPSPS gene was duplicated from two- to over 100-fold in comparison with a GS population. Both EPSPS mRNA and EPSPS protein expression correlated positively with genomic EPSPS copy number in the Georgia GR population55 and in a different GR population from Mississippi.56 In the Georgia A. palmeri population, the EPSPS gene was present on all chromosomes. Subsequently, a North Carolina GR A. palmeri population was reported with EPSPS gene duplication from 22- to 63-fold.57 A New Mexico GR A. palmeri population was sevenfold more resistant than a GS population and had EPSPS gene duplication ranging from two- to tenfold, and individuals with EPSPS gene duplication were able to survive label-rate glyphosate applications (0.82 kg ha−1).58,59

Additional species have recently been reported with EPSPS gene duplication in GR populations, including L. multiflorum,60 Kochia scoparia,61 A. tuberculatus62 and A. spinosus63 (summarized in Table 1). The reported number of duplicated EPSPS gene copies in K. scoparia and A. tuberculatus populations is lower (from four- to tenfold) than the Georgia, North Carolina and Mississippi55–57 A. palmeri populations (Table 1), and similar to the New Mexico A. palmeri population.58 EPSPS mRNA expression and EPSPS protein levels were found to have a linear correlation with EPSPS genomic copy number in L. multiflorum populations,60 as well as in K. scoparia populations.61

Table 1

EPSPS gene duplication reported in glyphosate-resistant weed species

Species	Population origin	EPSPS relative genomic copy number range	Reference
Amaranthus palmeri	Georgia	40–100	55
Amaranthus palmeri	North Carolina	20–60	57
Amaranthus palmeri	New Mexico	2–10	58
Amaranthus palmeri	Mississippi	33–59	56
Amaranthus tuberculatus	Missouri, Illinois	4	62
Lolium multiflorum	Arkansas	15–25	60
Kochia scoparia	Kansas, Colorado	3–9	61
Amaranthus spinosus	Mississippi	26–37	63

A. spinosus and A. palmeri have been shown to be closely related using both molecular markers and gene sequence phylogeny.64,65 As glyphosate resistance can transfer from A. palmeri to A. spinosus through interspecific hybridization and the hybrid progeny are self-fertile,66 an important question is whether A. spinosus with EPSPS gene duplication63 has evolved independently, or whether the trait has transferred via interspecific hybridization.

The glyphosate resistance level appears to increase with higher EPSPS genomic copy number in several examples. In L. multiflorum, the dose required to achieve 50% reduction in plant growth increased in a linear relationship with EPSPS protein expression.60 Higher EPSPS copy number correlated with higher-level resistance in A. palmeri populations from both Georgia and Mississippi.56,67 Shikimate accumulation following glyphosate treatment decreased as EPSPS cDNA expression level increased in a New Mexico A. palmeri population.58 An open question is whether EPSPS genomic copy number in species such as K. scoparia and A. tuberculatus may respond to continuing glyphosate selection pressure and increase to higher levels than previously observed. The number of species using gene duplication begs the question as to how common the mechanism may be and why it has not been observed for other herbicide mechanisms of action? Gene duplication could be a possible explanation for the increased gene expression associated with rapid vacuole sequestration, or any trait that could be conferred by increasing gene expression through increased gene copy number.

3.2 Molecular mechanism and inheritance of EPSPS gene duplication

The A. palmeri genome likely contains two original EPSPS loci, and one of the loci appears to have been duplicated in GR populations.68,69 Genomic sequencing suggested that the gene duplication proceeded via a DNA-mediated mechanism, and the duplicated section of DNA including the 10 kb EPSPS gene is at least 30 kb long.68 In that study, various mobile genetic elements were identified within the duplicated sequence bordering EPSPS gene copies, including an unclassified DNA transposon, an Activator transposase and sequences with homology to miniature inverted-repeat transposable elements (MITEs). However, it remains unclear whether the gene duplication mechanism involves MITEs or any other mobile genetic elements.

EPSPS gene duplication inheritance appears to be more complex than inheritance of resistance-endowing single-gene mutations. The EPSPS gene duplication in several A. palmeri populations is inherited on the nuclear genome.56–58,67 Inheritance in a North Carolina A. palmeri population was generally consistent with polygenic inheritance and did not fit a single-gene model.57 EPSPS gene duplication inheritance in a New Mexico A. palmeri population was likewise not consistent with a single-gene model.58 Inheritance of high EPSPS copy number into F1 and F2 generations from a Georgia population was quantitative and not always stable, including transgressive segregation and non-Mendelian inheritance in at least one F2 family.69 This variable gene duplication transmission has parallels with the instability of duplicated EPSPS copies from glyphosate-selected cell cultures.54 EPSPS gene duplication inheritance in a Mississippi A. palmeri population showed maternal influence, as a higher copy number was inherited in progeny with a GR maternal parent than in progeny from a GS maternal parent.56 Evidence for apomixis was shown by the production of seeds on maternal A. palmeri individuals in the complete absence of pollination.56 Such apomixis may increase the frequency of the maternal genotype in the progeny, in part explaining the maternal influence on EPSPS copy number observed in the inheritance studies. Apomixis could be very advantageous for the evolution of gene duplication, ensuring that the selected trait is maintained in the population. Additional factors that may contribute to the observed unstable inheritance include ongoing activity of the duplication mechanism, protective genomic mechanisms that remove duplicated EPSPS gene copies and/or other unknown mechanisms. Notably, a lack of fitness cost due to EPSPS gene duplication has been reported, accounting for variations in genetic background and measuring multiple phenotypic parameters.70,71 Clearly, the inheritance and population genetics of EPSPS gene duplication are more complex than most known herbicide resistance mechanisms.

4 EPSPS TARGET-SITE MUTATIONS

The EPSPS gene in plants contains coding sequence for a transit peptide used to direct the protein to the chloroplast.72,73 The transit peptide, variable in sequence and length among different species, is cleaved upon chloroplast delivery, producing the mature EPSPS protein. For consistency and ease of comparison, most publications number the amino acid positions relative to the start codon of the mature plant EPSPS protein, without the chloroplast transit peptide leader sequence.

Reported mutations conferring glyphosate resistance in weeds change the hydrophobic Pro106 amino acid of EPSPS to the hydrophobic amino acids Ala or Leu, or the hydrophilic amino acids Ser and Thr (summarized in Table 2). The Pro106 codon (CCx) can mutate to the commonly reported Ser, Ala or Thr through substitutions at the first base of the codon (TCx, GCx and ACx respectively). Substitutions at the second base of the codon can produce Leu (CTx, the most infrequently reported mutation), Arg (CGx), Gln (CAA/G) or His (CAC/T). Pro106Arg, Gln and His mutations have not been reported to date. While Gln is hydrophilic, both Arg and His are positively charged and may be disruptive to the active site. EPSPS target-site mutations have been reported in six species, most frequently in the genus Lolium, including L. rigidum and L. multiflorum (Table 2). The first dicotyledonous species with a Pro106 mutation is A. tuberculatus, with two populations recently reported.74,75

Table 2

Reported Pro106 target-site mutations in EPSPS endowing glyphosate resistance in weed species

Species	Reference	Pro106a to	Fold resistance	Other mechanisms detected?
Eleusine indica	16	Ser	2–4	No
Amaranthus tuberculatus	74	Ser	5	Yes, reduced translocation
Amaranthus tuberculatus	75	Ser	5	No, Ser mutation did not fully account for resistance
Echinochloa colona	108	Ser	6.6	No
Lolium multiflorum	109	Ser	2–5	No
		Ala	5–15
Eleusine indica	110, 111	Ser	3	No
		Thr	3
Digitaria insularis	88	Thr	4	Yes, reduced absorption and reduced translocation
Eleusine indica	112	Ser	2	No
Lolium rigidum	107	Leu	1.7	Yes, unknown mechanism
Lolium rigidum	106	Ser	6–8	Yes, reduced translocation
		Thr	8–11	Yes, reduced translocation
Lolium multiflorum	113	Ser	5	Yes, reduced translocation
Lolium multiflorum	114	Ser	5	No, reduced translocation detected in different population
Lolium rigidum	115	Thr	2–3	No
Lolium rigidum	34	Ser & Leu	16–21	No, but other mechanisms suspected
Lolium rigidum	116	Ser	—	No
Lolium rigidum	117	Ala	14	Yes, reduced translocation

Numbered relative to the start of the mature EPSPS enzyme in plants; for example, the first ten amino acids of the mature petunia EPSPS (GenBank M21084.1) are KPSEIVLQPI, and for the Arabidopsis EPSPS gene AT2G45300 the first ten are KASEIVLQPI.

The Pro106 is not directly involved in molecular interactions with either glyphosate or the substrate PEP, but changing the Pro106 to a different amino acid causes a structural change in the active site and shifts other amino acids towards the inhibitor, reducing the available space.76 The Pro106Ser mutation has previously been shown to increase the inhibitory constant Ki(app) for glyphosate from 0.4 μM in the wild-type petunia (Petunia × hybrida) enzyme to 3 μM in the mutant, but also to increase the Km(app) for PEP from 5 to 44 μM.72 The ratio Ki/Km, considered to be a selectivity factor for PEP over glyphosate binding,72 showed little change due to the Pro106Ser mutation, from 0.08 in wild type to 0.07 in the mutant. Subsequent research using a Pro106Ser mutant cloned from E. indica showed a wild-type Ki/Km ratio of 0.013 and a Pro106Ser mutant ratio of 0.12, a 9.2-fold increase.16 The expected loss of substrate affinity was less severe than reported in petunia,72 suggesting that in the E. indica EPSPS the Pro106Ser mutation confers glyphosate resistance sufficient to confer survival in the field while maintaining adequate PEP affinity.16

Additional mutations in the EPSPS gene able to confer glyphosate resistance have been identified using directed mutagenesis and expression in E. coli and transgenic plants, but have not yet been reported in GR weeds (Table 3). These mutations should be considered as candidate glyphosate resistance mechanisms when considering EPSPS target-site resistance by sequencing the entire EPSPS gene. Several reports involve double mutations that provide high-level glyphosate resistance while maintaining sufficient affinity for PEP (Table 3). The double mutants Gly101Ala + Gly144Asp77 and Gly101Ala + Ala192Thr78 have been shown to confer glyphosate resistance in transgenic maize (Zea mays) and canola (Brassica napus) respectively. The Thr102Ile + Pro106Ser (TIPS) double mutation has been described in transgenic crops (GA21 event used in first-generation Roundup Ready maize) as providing high-level glyphosate resistance,79,80 and the enzyme kinetics shows improved glyphosate resistance and relatively acceptable PEP affinity in comparison with the Pro106Ser mutation alone.81 Recently, Jalaludin et al.82 reported a GR E. indica population that had the Thr102Ile mutation in combination with the Pro106Ser mutation, which had previously evolved. This first report of a double EPSPS mutation in a field-evolved GR weed population demonstrates the evolutionary process by which mutations can accumulate to confer increasingly efficient resistance when selection pressure is persistent. Consistent with empirical and theoretical predictions for evolution of fitter proteins,83,84 the Thr102Ile mutation would be unlikely to occur first or independently owing to its significant reduction in PEP Km (Table 3). However, as a second mutation, Thr102Ile can improve the Pro106Ser mutation that has evolved first in many GR species.

Table 3

Additional target-site mutations in EPSPS shown to confer glyphosate resistance (increased Ki for glyphosate) with variable effects on PEP affinity (Km) when expressed in E. coli. Numbering system relative to the start of the mature EPSPS enzyme in plants, including petunia and Arabidopsis. Positions Gly101, Thr102, Pro106, Gly144 and Ala192 in plant mature EPSPS consensus correspond to Gly96, Thr97, Pro101, Gly137 and Ala183 in E. coli

Mutation(s)	Ki(app) for glyphosate (μM)	Km (app) for PEP (μM)	Ki/Km	Reference
Petunia wild-type EPSPS	0.4	5.2	0.08	118
Gly101Ala	2000	200	10.0	118
Gly101Ala, Ala192Thr	683	54	12.6	118
Gly101Ala, Gly144Asp	348	40	7.7	119
Gly101Ala, Gly144Asn	960	91	10.5	119
Maize wild type EPSPS	0.5	27	0.02	120
Maize Thr102Ile	148.6	233.0	0.6	120
Maize Pro106Ser	1.0	17.1	0.06	120
Maize Thr102Ile, Pro106Ser	58.0	10.6	5.5	120
Maize Thr102Ile, Pro106Thr	101.3	11.2	9.0	120
Maize Thr102Ile, Pro106Gly	38.6	23.0	1.68	120
Maize Thr102Ile, Pro106Cys	818.2	47.0	17.4	120
Maize Thr102Ile, Pro106Ala	148.3	10.2	14.5	120
Maize Thr102Ile, Pro106Ile	2500	60.3	41.5	120
Maize Thr102Ile, Pro106Val	1600	109.3	14.6	120
Maize Thr102Ile, Pro106Met	37 200	143.3	260	120
Maize Thr102Ile, Pro106Leu	2100	99.5	21.1	120
Maize Pro106Thr	4.0	24.6	0.16	120
Maize Pro106Leu	28.6	86.7	0.33	120
Agrobacterium spp. CP4	5100	14.4	354	120

5 OTHER MECHANISMS

Because very high (e.g. >50× the label rate) resistance levels have not been reported in GR weeds as has been observed in some cases for herbicides that inhibit ALS or ACCase, it is hypothesized that low dose selection has been effective for glyphosate resistance mechanisms. Low dose recurrent selection is a major scientific methodology for investigating pesticide resistance evolution, including herbicides,85–87 primarily because weak resistance traits are successfully selected with a low dose. Cross-pollination among survivors of a low dose followed by recurrent selection combines additive contributions from genes of small effect to evolve commercially significant resistance.86 Such traits involving any mechanism that contributes to improved survival under glyphosate selection have appeared in the literature. Reduced glyphosate absorption has been reported in some GR populations.88,89 In the case of Chilean L. multiflorum, the leaf angle in the GR population was different from that in the GS population, enabling reduced foliar spray interception.89 In addition, the GR population absorbed less glyphosate and had lower glyphosate translocation from treated leaves.

The 31P NMR survey of species (Section 2.1.1) in part chose D. illinoensis and S. obtusifolia because of the report of approximately 50% metabolism of glyphosate to AMPA in 7 days.90 AMPA also has a phosphonate, making it visible by 31P NMR, and the signal is easily distinguished from glyphosate. However, in the short time of these 31P NMR experiments, no AMPA could be observed (Sammons RD, unpublished). The metabolism of glyphosate in plants at least is rare and is only documented in a few cases, such as Equisetum arvense91 and soybeans.92 Examination of glyphosate metabolism in multiple species found no relationship between conversion to AMPA and glyphosate resistance level.90,93 Some recent reports employing a new methodology claim that two metabolic pathways for glyphosate are present in C. canadensis and Digitaria insularis.88,94 The detection of both glyoxylate and sarcosine in a glyphosate-dependent manner, when these compounds can also be readily metabolized in plants, requires confirmation by other methods. The unequivocal observations that weedy legumes90 can oxidize glyphosate to AMPA will spur on others to reinvestigate this long-overlooked mechanism and suggest that glyphosate degradation mechanisms may contribute to resistance eventually in selected species, particularly legumes.

6 INTEGRATIVE MECHANISTIC MODEL FOR GLYPHOSATE RESISTANCE

The concept starts with the realization that a finite dose of glyphosate is delivered to a plant, more particularly to the source leaves on the plant. Previous research has shown that sink leaves do not take up appreciable glyphosate and are incapable of transporting it systemically.95 The idea that the dose is finite builds on the loading that eventually ends. Many different formulations were evaluated,96 and while the kinetics of uptake varied among formulations, no formulation continued to complete glyphosate loading, indicating that there is a time period for loading and it does end. The model consequently has each source cell presented with a finite dose of glyphosate, and the fate of the glyphosate in the cell would depend on the pathways available once inside. Figure 1A shows theoretical glyphosate uptake and fate within a normal source-leaf cell. The y-axis is a measure of the glyphosate concentration in a relative sense, and the x-axis is time limited to 1 day. The accumulation of glyphosate in the cytoplasm presents glyphosate to organelles and most particularly the chloroplast, which must accumulate glyphosate. Chloroplast accumulation rates depend on the concentration in the cytoplasm (whether uptake is passive or active transport). Because most source cells supply assimilate, usually sucrose, to the phloem by way of cell-to-cell connectivity via plasmodesmata, glyphosate phloem mobilization (translocation) is dependent on the glyphosate concentration in the cytoplasm.97 The model needs some leniency from the reader for the ‘mobilized’ curve, as it really represents the rate of glyphosate being translocated out of the cell rather than a concentration in the phloem. The applied glyphosate dose is finite, limited by the concentration-volume of the application solution and the loading competency of the surfactant system.96 The amount of glyphosate entering the source cell slows and stops at about 10 h (Fig. 1A). The cell then continues to disperse the cytoplasmic glyphosate to the organelles and the plasmodesmata-phloem until some equilibrium is reached or the cytoplasm is drained of glyphosate. The concentration of glyphosate in the chloroplast determines the herbicidal activity, as partial EPSPS inhibition results in varying degrees of stunting and sublethal effects on the cell and the leaf. A concentration higher than the minimum inhibitory concentration sustained for an adequate period of time inhibits the EPSPS and blocks the shikimate pathway sufficiently to destroy the chloroplast (Fig. 1A). Because chloroplast destruction is fast98 and appears as swollen plastids, the idea that chloroplasts have lost their osmotic integrity seems likely. The potential phytotoxic action of rapid chloroplast destruction differs from aromatic amino acid starvation often attributed to the glyphosate mechanism of action (slow death due to lack of aromatic amino acids). The actual cellular toxic mechanism is not clear, although it could be hyperaccumulation of shikimate faster than it can be exported99 or some other indirectly linked process interrupted.100 However, it should be pointed out that killing a cell or a leaf, or a root tip, etc., is not sufficient to kill a plant, as this is only part of the consequential biochemical activity, and so the model does not mean to explain overall plant death.

Figure 1

A theoretical cellular-level model of glyphosate uptake and distribution for (A) a normal, glyphosate-susceptible source cell and (B) a glyphosate-resistant source cell using the vacuole to sequester glyphosate. The units for glyphosate are relative concentration, with a theoretical chloroplast minimum inhibitory concentration (25) indicated by a red line, and a chloroplast glyphosate concentration (35) consistent with saturated inhibition indicated by a dashed blue line.

Figure 1B represents a GR cell where, in addition, the vacuole is used to sequester glyphosate. The concentration of glyphosate in the cytoplasm of the GR source-leaf cell (Fig. 1B) is somewhat less at any given time than in the cytoplasm of the GS source-leaf cell (Fig. 1A) because there is now an extra pathway competing for glyphosate. The cytoplasm might be considered to be a conduit for the other destinations, so the total flux (loading and dispersal) is the same. The vacuole captures most of the glyphosate entering the cell, thereby decreasing the amount available to the chloroplast and phloem. The three pathways together decrease the amount of glyphosate in the GR cytoplasm faster than in the GS case, resulting in glyphosate concentrations lower than the minimum inhibitory concentration in the chloroplast and less glyphosate being mobilized (Fig. 1B). Resistance results, then, from a lower cytoplasm glyphosate concentration and less time to load the chloroplast. Therefore, Fig. 1B demonstrates the case where sequestration would provide a shield to the plant by decreasing the net available dose of glyphosate. The glyphosate trapped in the vacuole might be considered to be a Trojan horse, but this never seems to be realized in surviving herbaceous plants.

The model also works for sink tissue. In this case, glyphosate is delivered to the cytoplasm symplastically via the phloem, as long as the phloem is delivering assimilate.101 With decreased vacuole size in sink tissue, the volume of the cytoplasm is a much larger relative proportion, holding more glyphosate and allowing for an increased dose delivered to the chloroplasts. Sink tissues seem to be more reliant on their chloroplasts, as each unit of glyphosate is more toxic in sink tissues,98,102 but this could be explained by less organelle competition for glyphosate, resulting in much higher concentrations in the developing sink-leaf chloroplasts. One argument could be that the increased sensitivity of apical tissue would stop phloem delivery due to toxicity (blocked respiration of sucrose), leaving a high concentration of glyphosate in the source-leaf cytoplasm that could further load the chloroplasts there with increased time. This seems possible, but still it remains a fact that source leaves are generally the last green tissues in a lethally dosed glyphosate plant.

The model works very well for explaining the data observed by 31P NMR, 14C glyphosate and shikimate determinations, so the consequence of other potential resistance mechanisms can be examined. Firstly, an assumption illustrated in the model needs to be highlighted because of its critical impact on the predictions. The 31P NMR experiments demonstrate that glyphosate only crosses a lipid bilayer in one direction in the time course of 2–3 days (for living cells).17 This means that, once entering the cell, glyphosate cannot diffuse back to the apoplast, and similarly, once in the vacuole, glyphosate does not go back to the cytoplasm. For this reason it has been assumed that glyphosate entering the chloroplast also cannot go back to the cytoplasm. This makes the modeling somewhat easier, as a reversible equilibrium need not be included, but it can be accommodated if necessary when data are produced that support the efflux.

Why is the very weak mutation Pro106 to Ser or Thr found so often in so many species? The mutation itself on purified enzyme makes very slight changes to the Km of PEP and Ki for glyphosate. The experimentally determined EPSPS I50 values will be used primarily because the interactions of Ki and Km for mutants is beyond the scope of the present discussion.103 If the glyphosate concentration in the chloroplast is well below or well above the I50, then the mutant EPSPS may be of very little consequence, but small differences leading to life or death could be observed at chloroplast glyphosate concentrations near the EPSPS I50.104 Mutant EPSPS would effectively raise the minimum inhibitory concentration shown in Fig. 1. Because the Pro106 mutation is consistently selected in field situations, chloroplast glyphosate concentrations must be near the EPSPS I50. While the Pro106 mutation is weak, the addition of any other mechanism that lowers the glyphosate cytoplasm concentration, and therefore reduces the chloroplast glyphosate concentration below the mutant EPSPS I50, would have an additive interaction beneficial for survival. A good example described is when this weak Pro106Ser mutation was combined with sequestration as observed in L. multiflorum, R336.32,34 Such Pro106 mutations are frequently reported, which may suggest the presence of other contributing mechanisms. Obviously, a potent mutation that significantly raises the Ki will easily be selected, as the chloroplast is unlikely to accumulate such very high levels of glyphosate in the time allowed (Powles S, private communication). Higher EPSPS concentrations in the chloroplast would increase the total units of enzyme activity, requiring more glyphosate to obtain sufficient EPSPS inhibition. The I50 remains the same, but enough EPSPS activity remains for plant survival even when EPSPS is 95% inhibited. The EPSPS gene duplication system recently characterized in several species provides more enzyme, which requires more glyphosate in the chloroplast for sufficient inhibition of EPSPS.55

Other barriers in the model can easily be envisioned. If glyphosate has reduced uptake by a change in plasmalemma transport, then the cytoplasm concentration is decreased, limiting access to the chloroplast as measured for S. halepense.17,38 Similarly, if glyphosate uptake into the chloroplast is reduced, then the effective concentration of glyphosate could be drastically reduced as proposed for D. illinoensis.17 The combination of reduced cellular uptake and reduced chloroplast uptake together with vacuole sequestration would create a series of impressive barriers to delivering glyphosate.

This model has two unique perspectives built in. Firstly, there is a finite dose that must be dealt with by the plant at the cellular level. Secondly, there is a time component, as the limited amount of glyphosate available must reach the minimum inhibitory concentration (shown in Fig. 1) within a particular timeframe. Very many sequential barriers could decrease the glyphosate available to reach the chloroplast, and even physical cuticle barriers could play a role.89 Therefore, a number of specific physical and biochemical barriers can reduce the effective glyphosate available to inhibit EPSPS. The NMR data for a survey of cellular uptake suggests that there are three barriers used in very different ways by plants.17 Some reduce the cellular uptake, some sequester glyphosate using the vacuole and others seem to be able to restrict the entry of glyphosate into the chloroplast.

7 SUMMARY

The study of resistance mechanisms in GR weeds has provided fascinating insights into the evolutionary processes utilized by genetically diverse weed species in response to intense selection pressure. Previously unknown resistance mechanisms have been discovered in GR weed species, including vacuolar glyphosate sequestration,18 gene duplication55 and rapid mature leaf necrosis resulting in reduced translocation.43 Target-site mutations for glyphosate resistance have been reported at one codon, are apparently rare within populations and endow lower levels of resistance than target-site mutations for some other herbicide mechanisms of action. A common trend is for glyphosate resistance mechanisms to combine within populations and within individuals.105 Several examples have been reported of Pro106 target-site mutations acting additively with reduced translocation, or other unknown mechanisms, to provide higher-level resistance than either mechanism alone.74,106,107 In other cases, gene duplication or target-site mutations were detected but did not seem fully to explain the higher than expected levels of resistance.34,75 Accumulation of multiple mechanisms is a common occurrence, particularly in cross-pollinated species, producing higher resistance levels in response to continuing glyphosate selection pressure. The recruitment of a variety of biochemical processes to establish resistance to glyphosate in so many species suggests that weed management will have to diversify and use additional methods to control glyphosate-resistant weeds.