Mechanism(s) of Glyphosate Resistance in a Selected Plantago lanceolata (L.) R Biotype

: In 2003, a glyphosate-resistant plantago ( Plantago lanceolata L.) population located in the Robertson district of South Africa was subjected to different glyphosate dosages and the highest dosage (7200 g a.e. ha − 1 ) gave no acceptable levels of control. Here we reconﬁrm resistance and investigate the mechanism of glyphosate resistance. Dose-response curves indicated that the glyphosate dosage rate causing 50% survival (LD 50 ) for the resistant (R) biotype is 43 times greater than for the susceptible (S) biotype, i.e., 43-fold resistant to glyphosate. Investigation into the molecular mechanism of plantago showed shikimate accumulation of the R biotype was lower than that of the S biotype. The reported 31 P and 13 C nuclear magnetic resonance (NMR) spectra show rapid glyphosate translocation into the young untreated leaves of the S biotype. No glyphosate translocation was observed in the R biotype. A point mutation in the 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) gene, resulting in an amino acid substitution was also observed, indicating two unique glyphosate resistance mechanisms within the R biotype. The rapid evolution of glyphosate-resistant weeds threatens the usage of the world’s most important herbicide (glyphosate), which is essential in world food production and further limits grower options for weed control. New weed management strategies will be necessary to combat plantago R biotypes.

Glyphosate resistant weed species have become very common and they threaten glyphosate-based weed management strategies. This is because glyphosate is the world's most important herbicide, and is used worldwide to control a broad spectrum of weeds in various cropping systems. Consequently, high selection pressure from glyphosate abuse has led to the evolution of resistance to glyphosate in weeds [3][4][5][6][7]. Glyphosate resistance was first reported in 1996 in an apple orchard. Since the development of glyphosate, more than 71 different countries have reported glyphosate resistance in more than 52 weed species, and the number of glyphosate resistant weed species is expected to increase [3,7]. Although weeds can be destructive in various situations, the ones that have a major economic impact are those in glyphosate-resistant crops. Moreover, glyphosate is widely used to control weeds that are already resistant to other herbicides. Because of this, glyphosate resistant weeds pose a serious threat to sustainable weed management in agronomic crops.

Dose Response Trial
Seed from a susceptible (S) plantago biotype NW (S) and from a putatively resistant R2 (R) plantago biotype located in the Robertson district of South Africa was collected in 2018 and used in the current study ( Table 1). The growers in the Robertson area usually apply glyphosate and paraquat in rotation to control weeds. Plants from the resulting F2 generation were germinated in sand and also used for subsequent experiments. After germination, the seedlings were transplanted in small plastic pots filled with coarse gravel. The seedlings were placed in a glasshouse and irrigated with a standard Steiner nutrient solution [23]. Glasshouse was set at 20/25 • C, (12/12 h night/day) photoperiod, and 1500 mmol m −2 s −1 PAR. At the 2-3 leaf stage, the plants were subjected to glyphosate (Roundup ® Turbo, 450 g a.e. L −1 Bayer Crop Science (Pty) LtD, Isando, South Africa) applications in a dose response trial using a custom built, cabinet pneumatic sprayer equipped with a flat fan nozzle (Stellenbosch University, Department of Engineering) delivering herbicide at 108 L ha −1 at 200 kPa, with a speed of 1 m s −1 . Glyphosate was applied at dosage rates of 0, 270, 540, 1080, 2160, 4320, 8640, and 17,280 g a.e. ha −1 . Water delivery rate was 100 L ha −1 . The dose response trial was a 2 × 8 factorial arranged in a completely randomized block design (CRBD), replicated six times. Evaluation took place 42 days after herbicide application, when survival rate and dry mass were determined. Plants were recorded as alive or dead if they exhibited clear or lack of regrowth 42 days after herbicide application.

Shikimic Acid Analysis
High shikimate accumulation in leaves is an indicator of susceptible EPSPS enzyme and the reverse is true for R biotypes [5]. To document this, resistant and susceptible plantago biotypes sprayed with glyphosate at 2160, 4320, 8640, and 17,280 g a.e. ha −1 were harvested after one week. Leaf tissue of treated and control plants were collected from each biotype. The plants were harvested and dried at 70 • C. This was due to very high shikimate accumulation after six [24] and seven days (168 h) [25], as reported in literature. The plants were then ground to a fine powder using a milling machine (Peter Rassloff Instruments and Services Pty (Ltd.), Cape Town, South Africa). One mL of 70% methanol (MeOH) was added to ca. 10 mg of each sample. The samples were vortexed and extracted in the oven at 60 • C for three hours. Subsequently, 250 µL of each sample was transferred into a two mL vial and completely dried under a gentle stream of nitrogen. The dried samples were derivatized with 100 µL of acetonitrile (ACN) and 50 µL N, O-Bis(trimethylsilyl)trifluoroacetamide (BSTFA) at 80 • C for one hour. The samples were vortexed and transferred to an insert (positioned in a vial) and injected with a 10:1 split ratio into a GC-MS instrument (Agilent Technologies Inc., Santa Clara, CA, USA) [26].

Gas Chromatograph Mass Spectrometry (GC-MS)
Separation was performed on a gas chromatograph (6890 N, Agilent Technologies Inc., Santa Clara, CA, USA). coupled to an Agilent Technologies Inc. inert XL EI/CI Mass Selective Detector (MSD) (5975, Agilent Technologies Inc., Santa Clara, CA, USA). The GC-MS system was coupled to a CTC Analytics PAL autosampler. Separation of sugars was performed on a ZB-Semivolatile (30 m, 0.25 mm ID, 0.25 µm film thickness) Zebron 7HG-G027-11-GGA capillary column. Helium was used as the carrier gas at a flow rate of 1 mL/min. The injector temperature was maintained at 250 • C. One µL of the sample was injected in 10:1 split ratio. The oven temperature was programmed as follows: 100 • C for two minutes and then increased to 320 • C at a flow rate of 20 • C per minute. The mass selective detector (MSD) was operated in a full scan mode and the source and quad temperatures were maintained at 240 and 150 • C, respectively. The transfer line temperature was maintained at 250 • C. The mass spectrometer was operated under electron impact (EI) mode at an ionization energy of 70 eV, scanning from 35 to 600 m/z. Shikimic acid determination using GC-MS was carried out at the Central Analytical Facilities (CAF) specialized laboratories as per [26], but the method was modified substantially.

31 P and 13 C NMR
Small concentrations of glyphosate in plant cells can be determined by using 31 P and 13 C NMR procedures (Agilent Technologies Inc., Santa Clara, CA, USA) [15,27]. Specifically, for this study, due to the easy sample preparation, the translocation of glyphosate in leaves was monitored by solid state NMR [27][28][29]. Comparison of glyphosate translocation in Agronomy 2021, 11, 884 4 of 12 resistant and susceptible plantago biotypes dried leaf tissue was evaluated using both 31 P CP MAS and 13 C CP MAS NMR after the application of 13 C 2 -labeled glyphosate at CAF specialized laboratories.
Twelve resistant and susceptible plantago biotypes grown in small plastic pots in coarse gravel placed were used. A drop of 13 C 2 enriched glyphosate made up in commercial formulation of four times the recommended rate (540 g a.e. ha −1 ) [15]. Roundup ® Turbo was applied at the middle of a mature fully expanded leaf [2]. After two days, plants were harvested and sectioned into exposed (treated) leaf and protected (untreated) leaves. The untreated leaves were milled using a milling machine, then vacuum-dried immediately before being subjected to solid state 31 P and 13 C NMR spectrometry [30]. The 13 C 2 -enriched glyphosate was chosen because of the low natural abundance and receptivity of the 13 C isotope. The spectra were acquired using an Agilent VNMRS (Agilent Technologies Inc., Santa Clara, CA, USA) 500 MHz two-channel solid state NMR spectrometer using 4 mm zirconia rotors spinning at 10,000 Hz, utilizing a 4 mm T3 HX MAS probe. All crosspolarization (CP) magic angle spinning (MAS) spectra were recorded with the latest VnmrJ 4.2 (Agilent Technologies Inc., Santa Clara, CA, USA) instrument software at ambient temperature with optimized high-power proton decoupling; a relaxation delay of 2 s and 35,587 scans were collected to obtain sufficient signal-to-noise. Both experiments employed were sensitivity-enhanced cross-polarization experiments where the magnetization is transferred from the most NMR sensitive 1 H nucleus to either the 31 P or 13 C nuclei utilizing a 2 us proton excitation pulse width and 2 ms contact time, and a 20 ms acquisition time. The probe was fine-tuned to 1 H, 31 P, or 31 C for each sample.

EPSPS cDNA Sequencing
RNA was extracted using the NucleoSpin ® RNA Plant Kit (Macherey-Nagel, Düren, Germany). The RNA samples were quantified and assessed for quality using the Bioanalyzer RNA 6000 Nano Assay (Agilent Technologies, Santa Clara, CA, USA) and cDNA was synthesized from 200 ng RNA for each of the samples and extracted using the SuperScript™ VILO™ cDNA synthesis kit (ThermoFisher Scientific, Waltham, MA, USA). The EPSPS gene was amplified with 20 ng first-strand synthesized cDNA, using 10 pmol of each of the sense 5 -GATGTCGAGGAAGACAAAGC-3 and antisense 5 -CCTTCCAGAAGAGTTACGTG-3 primers in a final reaction volume of 25 uL, containing 200 µM of each dNTP and one-unit Platinum™ II Taq Hot-Start DNA Polymerase (ThermoFisher Scientific). The cDNA for the samples were amplified using standard polymerase chain reaction (PCR) conditions at an annealing temperature of 58 • C to yield a fragment of 795 bp. The resulting PCR fragment was purified with 1.8× volume Agencourt AMPure XP Reagent (Beckman Coulter, Brea, CA, USA). Direct cycle sequencing was performed on the purified PCR fragment using standard protocols (DNA sequencing Unit, Central Analytical Facilities, Stellenbosch University, Stellenbosch). DNA sequences from each of the plants were aligned to the partial plantago 5-enol-pyruvylshikimate-phosphate synthase (EPSPS) coding sequence (Genbank accession number: AY545665.1). For the multiple sequence alignment, the reference and alternative alleles were separated using the online analysis software Poly Peak Parser (http://yosttools.genetics.utah.edu/PolyPeakParser/) (accessed on 7 June 2019) [31]. Six independent RNA extractions from each biotype were analyzed and the sequence results were aligned and compared. Sequencing was conducted at CAF specialized laboratories.

Statistical Analysis
All of the experiments had six individual replicates, but each experiment was repeated twice. Nonsignificant differences were detected between repeated experiments; therefore, data from the two experiments were pooled and used for subsequent analysis. For the NMR studies, Mestrenova version 11.02 was employed (Mestrelab research 2020, Santiago de Compostela, Spain). The herbicide dosage rate causing 50% reduction in survival rate % (SR) (LD 50 ) or growth reduction (GR 50 ) was calculated for each biotype by a log-logistic model using STATISTICA 2019 software version 13.6 (TIBCO SOFTWARE Inc., 2019, Palo Alto, CA, USA).
Data sets were fitted to the log-logistic model: where C = lower limit, D = upper limit, b = slope, and I 50 = dosage rate yielding 50% response [32].

Dose Response Trial
The dose response trial established 0% survival of a known plantago S biotype at less than half the recommended glyphosate dosage rate (130 g a.e. ha −1 ). On the contrary, putative R biotype was markedly less affected; 5775 g a.e. ha −1 was needed to achieve LD 50 ( Table 2 and Figure 1). The resistance index (RI) of 43 for the R biotype in response to glyphosate applications was noted ( Table 2). The putative resistant R2 (R) biotype was therefore confirmed to be resistant to glyphosate. Glyphosate rate causing 50% reduction of growth for the R biotype was found to be 784 g a.e. ha −1 , fivefold greater than the S biotype (134 g a.e. ha −1 ) ( Table 2 and Figure 2). STATISTICA 2019 software version 13.6 (TIBCO SOFTWARE Inc., 2019, Palo Al USA). Data sets were fitted to the log-logistic model: where C = lower limit, D = upper limit, b = slope, and I50 = dosage rate yielding sponse [32].

Dose Response Trial
The dose response trial established 0% survival of a known plantago S biotyp than half the recommended glyphosate dosage rate (130 g a.e. ha −1 ). On the contra tative R biotype was markedly less affected; 5775 g a.e. ha −1 was needed to achie ( Table 2 and Figure 1). The resistance index (RI) of 43 for the R biotype in resp glyphosate applications was noted ( Table 2). The putative resistant R2 (R) bioty therefore confirmed to be resistant to glyphosate. Glyphosate rate causing 50% re of growth for the R biotype was found to be 784 g a.e. ha −1 , fivefold greater tha biotype (134 g a.e. ha −1 ) ( Table 2 and Figure 2).

Shikimic Acid
The highest glyphosate dosage resulted in 93% S biotype shikimate accumulation compared to 36% accumulation by the R biotype. Accumulation of shikimate occurs at higher concentrations in S biotypes. Much larger glyphosate dosages are required before the same accumulation is seen in R biotypes ( Figure 3).

Shikimic Acid
The highest glyphosate dosage resulted in 93% S biotype shikimate accumulation compared to 36% accumulation by the R biotype. Accumulation of shikimate occurs at higher concentrations in S biotypes. Much larger glyphosate dosages are required before the same accumulation is seen in R biotypes (Figure 3).

Shikimic Acid
The highest glyphosate dosage resulted in 93% S biotype shikimate accumulation compared to 36% accumulation by the R biotype. Accumulation of shikimate occurs at higher concentrations in S biotypes. Much larger glyphosate dosages are required before the same accumulation is seen in R biotypes (Figure 3).

Relationship between Shikimate, Survival Rate, and Dry Mass
There was a significant and negative correlation between shikimate and the dose response results for the R biotype, specifically, survival rate and dry mass (Table 3). Therefore, Agronomy 2021, 11, 884 7 of 12 lower shikimate implies higher survival rate and dry mass in the R biotype. For the S biotype, due to low or lack of variance, correlation was not feasible.

Glyphosate Translocation Using 31 P and 13 C NMR
The 31 P CP MAS NMR spectra from the untreated leaf of the S biotype gave a glyphosate signal at 13 ppm, at the same chemical shift as the (standard) C 2 -labeled glyphosate and both relative to a H 3 PO 4 (aq) reference standard (Figure 4). Similarly, for 13 C NMR, the untreated leaf of the R biotype gave no glyphosate signal. A small glyphosate signal was observed for the S biotype, which was established at 50 ppm by using a 13 C 2 -labeled glyphosate standard ( Figure 5).

Relationship between Shikimate, Survival Rate, and Dry Mass
There was a significant and negative correlation between shikimate and the dose response results for the R biotype, specifically, survival rate and dry mass (Table 3). Therefore, lower shikimate implies higher survival rate and dry mass in the R biotype. For the S biotype, due to low or lack of variance, correlation was not feasible.

Glyphosate Translocation Using 31 P and 13 C NMR
The 31 P CP MAS NMR spectra from the untreated leaf of the S biotype gave a glyphosate signal at 13 ppm, at the same chemical shift as the (standard) C2-labeled glyphosate and both relative to a H3PO4 (aq) reference standard (Figure 4). Similarly, for 13 C NMR, the untreated leaf of the R biotype gave no glyphosate signal. A small glyphosate signal was observed for the S biotype, which was established at 50 ppm by using a 13 C2-labeled glyphosate standard ( Figure 5).

EPSPS cDNA Sequencing
To identify DNA variants in resistant and susceptible plantago biotypes, the partial coding sequence of the EPSPS gene was PCR amplified and sequenced (Figures 6 and 7). A nucleotide variant, resulting in a predicted Proline (Pro) to Serine (Ser) substitution at position 106 of the amino acid sequence of EPSPS [31], was identified in the R biotype ( Figure 6). The reference and alternative alleles were included in the multiple sequence alignment against the plantago EPSPS gene reference sequence (Figure 7). A Cytosine to Adenine transition was observed in the resistant biotype ( Figure 6) and in the susceptible alternative allele (Figure 7). Cytosine to Thymine transition DNA variant, predicted to result in a substitution of a Proline at position 106 with a Serine, is denoted by the red arrow (Figures 6 and 7).
A partial alignment of the EPSPS amino acid sequence was performed for eleven Eudicots against the partial plantago reference sequence (GenBank accession number: AAT45242.1) (Figure 8). The alignment shows conservation of the Proline amino acid at position 106 (Figure 8). Included in this alignment is the Arabidopsis thaliana and Eleusine indica EPSPS sequences presented by [31]. The only amino acid change shown in the multiple sequence alignment for EPSPS at position 106 of the amino acid sequence was that of the plantago R biotype (Figure 8).

EPSPS cDNA Sequencing
To identify DNA variants in resistant and susceptible plantago biotypes, the partial coding sequence of the EPSPS gene was PCR amplified and sequenced (Figures 6 and 7). A nucleotide variant, resulting in a predicted Proline (Pro) to Serine (Ser) substitution at position 106 of the amino acid sequence of EPSPS [31], was identified in the R biotype ( Figure 6). The reference and alternative alleles were included in the multiple sequence alignment against the plantago EPSPS gene reference sequence (Figure 7). A Cytosine to Adenine transition was observed in the resistant biotype ( Figure 6) and in the susceptible alternative allele (Figure 7). Cytosine to Thymine transition DNA variant, predicted to result in a substitution of a Proline at position 106 with a Serine, is denoted by the red arrow (Figures 6 and 7).    Eleusine indica amino acid sequences, as described by [31], against that of plantago (Plantago lanceolata L.) EPSPS resistant and susceptible biotypes. Highlighted amino acids indicate similarities against the plantago reference sequence. The predicted nonsynonymous variant in the resistant plantago biotype is shown with a red arrow and results in the substitution of a Proline to a Serine at position 106 (Pro-106-Ser), as described by [31].

Discussion
The R2 biotype exhibited a high level of resistance, viz., R/S ratios of greater than 40 to glyphosate. These results are similar to Tehranchian et al. [33], but greater than those reported by Yu et al. [2], who reported an R/S ratio of 14 in a ryegrass population after glyphosate application, which is three times less resistant than plantago. These shows that glyphosate resistant weeds are increasing exponentially and becoming virulent [3]. Unfortunately, no recent studies have reported on plantago resistance to glyphosate for a direct comparison. There are currently more than 52 weed species resistant to glyphosate in more than 73 countries. This is likely to be an undercount since there may have been some reports that have not yet been reported [3,7].
The shikimate accumulation observed in the NW (S) biotype after glyphosate application serves as an indicator for susceptible EPSPS enzyme because the shikimate pathway is blocked [34]. Shikimate accumulation occurs in higher concentrations in S biotypes, relative to the R biotypes. This was also reported by several other authors [5,11,25,33,34]. The differences in shikimate accumulation is in line with dose response curve results, where the R biotype showed 43-fold less sensitivity to glyphosate compared to the S biotype. In addition, the survival rate and dry mass of the R biotype was found to be negatively correlated to shikimic accumulation. Accumulation of shikimic acid has been reported to be related to glyphosate dosages [35]. Lack of shikimate in the R biotypes could be because glyphosate does not reach the target sites due to sequestration or due to the (observed) mutations on the EPSPS enzyme [5,8,15,34].
Glyphosate moves through the phloem from high concentration (treated leaves) to Eleusine indica amino acid sequences, as described by [31], against that of plantago (Plantago lanceolata L.) EPSPS resistant and susceptible biotypes. Highlighted amino acids indicate similarities against the plantago reference sequence. The predicted nonsynonymous variant in the resistant plantago biotype is shown with a red arrow and results in the substitution of a Proline to a Serine at position 106 (Pro-106-Ser), as described by [31].

Discussion
The R2 biotype exhibited a high level of resistance, viz., R/S ratios of greater than 40 to glyphosate. These results are similar to Tehranchian et al. [33], but greater than those reported by Yu et al. [2], who reported an R/S ratio of 14 in a ryegrass population after glyphosate application, which is three times less resistant than plantago. These shows that glyphosate resistant weeds are increasing exponentially and becoming virulent [3]. Unfortunately, no recent studies have reported on plantago resistance to glyphosate for a direct comparison. There are currently more than 52 weed species resistant to glyphosate in more than 73 countries. This is likely to be an undercount since there may have been some reports that have not yet been reported [3,7].
The shikimate accumulation observed in the NW (S) biotype after glyphosate application serves as an indicator for susceptible EPSPS enzyme because the shikimate pathway is blocked [34]. Shikimate accumulation occurs in higher concentrations in S biotypes, relative to the R biotypes. This was also reported by several other authors [5,11,25,33,34]. The differences in shikimate accumulation is in line with dose response curve results, where the R biotype showed 43-fold less sensitivity to glyphosate compared to the S biotype. In addition, the survival rate and dry mass of the R biotype was found to be negatively correlated to shikimic accumulation. Accumulation of shikimic acid has been reported to be related to glyphosate dosages [35]. Lack of shikimate in the R biotypes could be because glyphosate does not reach the target sites due to sequestration or due to the (observed) mutations on the EPSPS enzyme [5,8,15,34].
Glyphosate moves through the phloem from high concentration (treated leaves) to lower concentration (untreated leaves) [36]. The observed 31 P CP MAS NMR spectra from the S biotype were corroborated by other authors who reported glyphosate signal at approximately 12 ppm [28]. Slight shift differences may be attributed to electron shielding differences between the crystalline and amorphous glyphosate used as the standard. Previous studies using 14 C-labeled glyphosate have shown that just after one day of application, about 35% of the glyphosate was translocated from the exposed leaf to the protected leaf [37]. This movement appears to be very quick in S biotypes compared to R biotypes, or perhaps absent in R biotypes [2] since no glyphosate signal was observed for the untreated leaf of the R biotype. This means glyphosate remained at the site of application in the R biotype, even after 48 h following glyphosate application. If there was any translocation, it was below the level of 31 P CP MAS NMR detection. Similarly, the observed 13 C NMR low glyphosate signal chemical shift was corroborated by [29]. The reason for a small signal is due to 13 C low absolute receptivity/natural abundance relative to 31 P, based on NMR frequency table. The R biotype gave no glyphosate signal. Reduced glyphosate translocation allows R biotypes to survive glyphosate applications by reducing the amount of glyphosate that moves from the treated leaves to the meristems [34]; this was certainly the case in plantago R biotypes.
In this study, a substitution of GCC to GCA was observed in the resistant biotype and in the susceptible alternative allele. However, the substitution did not result in amino acid changes. Silent mutations are plentiful in literature [2,5]. In contrast, Pro substitutes Ser; these point mutations very likely grants the resistance to glyphosate, together with the reduced translocation mechanism. This target-site resistance is usually accompanied by a nontarget site resistance mechanism [38]. Conversely, no other studies have reported on the mechanism of glyphosate resistance in plantago.
It was expected, however, that mutations of Pro-106-Ser would also occur in plantago, since it has been documented in other glyphosate resistant weed species [39]. Proline to Serine substitution at position 106 has been shown to cause glyphosate resistance EPSPS in bacteria (Escherichia coli) and weed species (Eleusine indica), among others [2,31]. Some other weed species have been found to possess TIPS, which is the double mutation of Thr-102-Ile and Pro-106-Ser, and grants a very high level of resistance [5]. Most recently, the so-called TAP-IVS triple mutation (Thr-102 + Ala-103 + Pro-106) was documented in Amaranthus hybridus [40].
In contrast to the aforementioned weed species, windmill grass (Chloris elata) R biotypes displayed no Pro-106-Ser substitution at any of the sequence position or any other amino acid substitutions (target site resistance), but showed differences in herbicide translocation (nontarget site resistance). This was the case in [12]. In this study, there were no mutations at the Pro106 codon in the gene encoding (EPSPS); however, the biotypes showed differences in translocation with the S biotype, showing higher glyphosate absorption into the leaf compared to the R biotype. In our study, mutation in the EPSPS enzyme of plantago was observed, pointing to target site resistance. However, the 31 P and 13 C NMR spectra showed reduced glyphosate translocation in the R biotype, implying a nontarget site resistance (presumably due to glyphosate sequestration). Worldwide, many species resistant to glyphosate are known [7]. However, this is the first report on the mechanism of glyphosate resistance in plantago.

Conclusions
Conclusively, 31 P NMR remains the preferred method to detect glyphosate even when 13 C enriched glyphosate is employed. The results unambiguously show the first evidence of glyphosate resistance mechanisms in plantago. Glyphosate resistance in this R biotype is due to a Proline to Serine substitution at amino acid position 106 of the EPSPS gene, and reduced glyphosate translocation to young leaves. This R biotype displayed a total of two mechanisms endowing glyphosate resistance, suggesting the presence of multiple mechanisms that confer glyphosate resistance in these R biotypes. This is very disturbing as it threatens the world's most important weed control resource (glyphosate). To preserve this essential herbicide resource, growers must be weaned of their total dependence on glyphosate. Herbicide rotation should be encouraged and other urgent alternative nonchemical control management strategies will also be necessary to combat resistant plantago. These results provide a basis for growers and industry to immediately amend how they perceive and control weeds.  Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.
Data Availability Statement: The datasets used and/or analyzed in the current study are available from the corresponding author on reasonable request.