Acessibilidade / Reportar erro

Transgenerational memory of drought stress and low rates of glyphosate reduce the sensitivity of Eragrostis plana to the herbicide

Abstract

Background

Eragrostis plana is the main invasive plant in the Pampa Biome of southern Brazil. This plant is highly competitive, tolerant to abiotic stresses, and very difficult to manage. Plants exposed to abiotic stresses and herbicides develop mechanisms that help alleviate or reduce damage caused by stressors and transmit this ability to the progenies.

Objective

The study’s objectives were to ascertain whether acclimatization to stress due to drought and sub-lethal doses of glyphosate may decrease plant sensitivity to glyphosate and investigate the possible memory mechanisms involved in this process.

Methods

A population of E. plana was submitted to drought, glyphosate, or a combination of drought plus glyphosate for two generations. The progenies were analyzed for sensitivity to the herbicide glyphosate and for biochemical, metabolites, and molecular responses.

Results

When testing sensitivity to glyphosate in the G 2 generations, the DRYxGLY (plants stressed with drought plus glyphosate) was the least sensitive (ED50 443.0 g a.e. ha-1), while CHK (without stresses) was the most sensitive (ED50 278.0 g a.e. ha-1). Evaluating the effects of glyphosate, the CHK population showed a greater accumulation of oxidative damage and shikimic acid. While DRYxGLY had greater antioxidant activity and higher expressions of the EPSPS and ABC-carrier MRP10 genes.

Conclusions

Recurrent selection with drought stress and sub-lethal rate of glyphosate (DRYxGLY) showed reduced sensitivity to glyphosate in the second generation ( G 2 ). The conjunction of factors, including the upregulation of EPSPS and the ABC MRP10 transporter, antioxidant enzymes seem to decrease the sensitivity of the DRYxGLY population to glyphosate.

South African lovegrass; Acclimatization; Shikimic acid; Stress memory

1.Introduction

Eragrostis plana (Ness), known as South African lovegrass or Tough lovegrass, is the primary invasive weed in over two million hectares of the Pampa Biome ( Medeiros, Focht, 2007Medeiros RB, Focht T. [Invasion, prevention, control and utilization of capim-annoni-2(Eragrostis plana Nees) in Rio Grande do Sul, Brazil]. Pesq Agrop Gaúcha. 2007;13(1-2):105-14. Portuguese. ) in southern Brazil. South African lovegrass is a slow-growing C4 grass with high seed production and dormancy, allelopathic activity, and tolerance to abiotic stresses ( Bastiani et al., 2021Bastiani MO, Roma-Burgos N, Langaro AC, Salas-Perez RA, Rouse CE, Fipke MV et al. Ammonium sulfate improves the efficacy of glyphosate on South African lovegrass ( Eragrostis plana ) under water stress. Weed Sci. 2021;69(2):167-76. Available from: https://doi.org/10.1017/wsc.2020.97
https://doi.org/10.1017/wsc.2020.97...
).

The most used herbicide to control E. plana in the grassland is glyphosate ( N -phosphonomethyl glycine). Glyphosate inhibits the enzyme 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS), reducing the synthesis of aromatic amino acids, which are intermediate compounds in the shikimate pathway (glycolytic and phosphate pathways), proteins, and secondary metabolism ( Maroli et al., 2015Maroli AS, Nandula VK, Dayan FE, Duke SO, Gerard P, Tharayil N. Metabolic profiling and enzyme analyses indicate a potential role of antioxidant systems in complementing glyphosate resistance in an Amaranthus palmeri biotype. J Agric Food Chem. 2015;63(41):9199-209. Available from: https://doi.org/10.1021/acs.jafc.5b04223
https://doi.org/10.1021/acs.jafc.5b04223...
). The shikimic acid pathway route expends around 30% of the fixed carbon, which is used by plants to synthesize intermediate compounds, auxins, lignins, secondary metabolism, and aromatic amino acids ( Maroli et al., 2015Maroli AS, Nandula VK, Dayan FE, Duke SO, Gerard P, Tharayil N. Metabolic profiling and enzyme analyses indicate a potential role of antioxidant systems in complementing glyphosate resistance in an Amaranthus palmeri biotype. J Agric Food Chem. 2015;63(41):9199-209. Available from: https://doi.org/10.1021/acs.jafc.5b04223
https://doi.org/10.1021/acs.jafc.5b04223...
). Therefore, interruption of the shikimate pathway by glyphosate triggers several metabolic and physiological processes, leading to plant death ( Gaines et al., 2020Gaines TA, Duke SO, Morran S, Rigon CAG, Tranel PJ, Küpper A et al. Mechanisms of evolved herbicide resistance. J Biol Chem. 2020;295(30):10307-30. Available from: https://doi.org/10.1074/jbc.REV120.013572
https://doi.org/10.1074/jbc.REV120.01357...
). In addition to the direct effects of EPSPS inhibition, glyphosate indirectly induces some ROS production, leading to lipid peroxidation, DNA, RNA, and protein oxidation, and slow cell death ( Gomes et al., 2017Gomes MP, Le Manach SG, Hénault-Ethier L, Labrecque M, Lucotte M, Juneau P. Glyphosate-dependent inhibition of photosynthesis in willow. Front Plant Sci. 2017;8:1-13. Available from: https://doi.org/10.3389/fpls.2017.00207
https://doi.org/10.3389/fpls.2017.00207...
; Maroli et al., 2015Maroli AS, Nandula VK, Dayan FE, Duke SO, Gerard P, Tharayil N. Metabolic profiling and enzyme analyses indicate a potential role of antioxidant systems in complementing glyphosate resistance in an Amaranthus palmeri biotype. J Agric Food Chem. 2015;63(41):9199-209. Available from: https://doi.org/10.1021/acs.jafc.5b04223
https://doi.org/10.1021/acs.jafc.5b04223...
).

Despite the efficient action of glyphosate, plants have evolved multiple mechanisms to mitigate the herbicide effects ( Gaines et al., 2020Gaines TA, Duke SO, Morran S, Rigon CAG, Tranel PJ, Küpper A et al. Mechanisms of evolved herbicide resistance. J Biol Chem. 2020;295(30):10307-30. Available from: https://doi.org/10.1074/jbc.REV120.013572
https://doi.org/10.1074/jbc.REV120.01357...
). The known resistance mechanisms are classified into target-site resistance (TSR) and non-target-site resistance (NTSR) mechanisms, being the latter the most complex, involving several genes and metabolic routes ( Gaines et al., 2020Gaines TA, Duke SO, Morran S, Rigon CAG, Tranel PJ, Küpper A et al. Mechanisms of evolved herbicide resistance. J Biol Chem. 2020;295(30):10307-30. Available from: https://doi.org/10.1074/jbc.REV120.013572
https://doi.org/10.1074/jbc.REV120.01357...
). TSR mechanisms include mutation and upregulation of the target gene ( Gaines et al., 2020Gaines TA, Duke SO, Morran S, Rigon CAG, Tranel PJ, Küpper A et al. Mechanisms of evolved herbicide resistance. J Biol Chem. 2020;295(30):10307-30. Available from: https://doi.org/10.1074/jbc.REV120.013572
https://doi.org/10.1074/jbc.REV120.01357...
). The main known NTSR mechanisms for glyphosate are reduction of absorption and translocation of the herbicide, metabolism of the molecule (forming mainly aminomethylphosphonic acid - AMPA), and exclusion of the herbicide (mainly in vacuole) ( Gaines et al., 2020Gaines TA, Duke SO, Morran S, Rigon CAG, Tranel PJ, Küpper A et al. Mechanisms of evolved herbicide resistance. J Biol Chem. 2020;295(30):10307-30. Available from: https://doi.org/10.1074/jbc.REV120.013572
https://doi.org/10.1074/jbc.REV120.01357...
).

Abiotic stress and herbicide sub-lethal doses lead to several plant responses, such as NTSR regulation and other pleiotropic effects ( Dyer, 2018Dyer WE. Stress-induced evolution of herbicide resistance and related pleiotropic effects. Pest Manag Sci. 2018;74(8):1759-68. Available from: https://doi.org/10.1002/ps.5043
https://doi.org/10.1002/ps.5043...
). These defense responses are similar and interconnected ( Burns et al., 2018Burns EE, Keith BK, Refai MY, Bothner B, Dyer WE. Constitutive redox and phosphoproteome changes in multiple herbicide resistant Avena fatua L. are similar to those of systemic acquired resistance and systemic acquired acclimation. J Plant Physiol. 2018;220:105-14. Available from: https://doi.org/10.1016/j.jplph.2017.11.004
https://doi.org/10.1016/j.jplph.2017.11....
). Therefore, it is expected that plants adapted to abiotic stresses may have greater tolerance to herbicide stresses and vice-versa ( Burns et al., 2018Burns EE, Keith BK, Refai MY, Bothner B, Dyer WE. Constitutive redox and phosphoproteome changes in multiple herbicide resistant Avena fatua L. are similar to those of systemic acquired resistance and systemic acquired acclimation. J Plant Physiol. 2018;220:105-14. Available from: https://doi.org/10.1016/j.jplph.2017.11.004
https://doi.org/10.1016/j.jplph.2017.11....
; Fipke et al., 2022Fipke MV, Feijó AR, Garcia NS, Heck T, Viana VE, Dayan FE et al. Transgenerational effect of drought stress and sub-lethal doses of quizalofop-p-ethyl: decreasing sensitivity to herbicide and biochemical adjustment in Eragrostis plana . Agriculture. 2022;12(3):396. Available from: https://doi.org/10.3390/agriculture12030396
https://doi.org/10.3390/agriculture12030...
). An increase in ROS-detoxifying enzymes promotes acclimation to drought and heat stress in Conyza bonariensis and Abutilon theophrasti , increasing the tolerance to paraquat and glyphosate ( Dyer, 2018Dyer WE. Stress-induced evolution of herbicide resistance and related pleiotropic effects. Pest Manag Sci. 2018;74(8):1759-68. Available from: https://doi.org/10.1002/ps.5043
https://doi.org/10.1002/ps.5043...
; Zhou et al., 2007Zhou J, Tao B, Messersmith CG, Nalewaja JD. Glyphosate efficacy on velvetleaf ( Abutilon theophrasti ) is affected by stress. Weed Sci. 2007;55(3):240-4. Available from: https://doi.org/10.1614/ws-06-173.1
https://doi.org/10.1614/ws-06-173.1...
).

As transcriptome, proteome, and plant metabolism changes caused by herbicide stress are similar to abiotic stresses, plants adapted to abiotic stresses are likely more tolerant to herbicides ( Yuan et al., 2010Yuan JS, Abercrombie LLG, Cao Y, Halfhill MD, Zhou X, Peng Y et al. Functional genomics analysis of horseweed ( Conyza canadensis ) with special reference to the evolution of non–target-site glyphosate resistance. Weed Sci. 2010;58(2):109-17. Available from: https://doi.org/10.1614/ws-d-09-00037.1
https://doi.org/10.1614/ws-d-09-00037.1...
). In this context, the study’s objectives were to ascertain whether the transgenerational effect of acclimatization to stress due to water deficit and sub-lethal doses of glyphosate may decrease herbicide sensitivity and investigate whether antioxidant enzymes and other resistance mechanisms are involved in this process.

2.Material and Methods

2.1 Establishment of the populations

Seeds of E. plana were collected in an experimental field (bulk of seed) in the Embrapa located in Bagé (RS Brazil) and denominated G 0 . In the greenhouse facility at UFPel, the populations G 1 and G 2 were generated in 2017 and 2018 ( Figure 1 ). In order to build the populations, the initial G 0 population seeds were divided into four groups (each group with eight plants) and exposed to the treatments: untreated Check (CHK - without stress); drought treatment (DRY - water deficit until 35% of stomatal conductance in relation to unstressed plants); glyphosate treatment (GLY - 120 g a.e. ha-1 of glyphosate); DRY×GLY (water deficit until 35% of stomatal conductance and after recovery, sprayed with 120 g a.e. ha-1 of glyphosate). The dose of glyphosate (dose to caused 30% to injury, data are not shown) or level of drought were defined in the pre-experiment (data are not shown), aiming at plant stress. The drought stress was performed according to the methodology described by Bastiani et al. (2021)Bastiani MO, Roma-Burgos N, Langaro AC, Salas-Perez RA, Rouse CE, Fipke MV et al. Ammonium sulfate improves the efficacy of glyphosate on South African lovegrass ( Eragrostis plana ) under water stress. Weed Sci. 2021;69(2):167-76. Available from: https://doi.org/10.1017/wsc.2020.97
https://doi.org/10.1017/wsc.2020.97...
. In each treatment group, eight plants were used. As E. plana is probably predominantly allogamous, the plants from each treatment were isolated (in different greenhouse, avoiding any crossing among treatments) and subsequently seeds were collected (in bulk), generating G 1 population. G 1 populations were planted and exposed to the same treatments to obtain G 2 generation.

Figure 1
Design of CHK (check), DRY (drought), GLY (glyphosate), and DRY×GLY (drought plus glyphosate) populations in two generations ( G 1 and G 2 ) from this G 0 E. plana plants

2.2 Dose-response curves with glyphosate

After obtaining all the population’s offspring, a growth chamber experiment was conducted, with experimental units consisting of 1 L pots filled with sieved soil (Haplosol) organized in a randomized block design with six replications. Factor A included the different populations G 0 , G 1 (CHK, DRY, GLY, and DRY×GLY) and G 2 (CHK, DRY, GLY, and DRY×GLY), and factor B was of glyphosate doses (0, 90, 120, 180, 250, 360, 540 and 720 a.e. g ha-1). All treatments were performed when plants reached an average of eight tillers, using glyphosate (Monsanto, Roundup Original™) applied with a CO2 sprayer calibrated to deliver 150 L ha-1 of spray solution. At 35 days after application (DAA), the visual injury (%) and shoot dry weight (SDW) ratings were analyzed. The SDW was transformed as the percentage compared to the plants without herbicide (%) of each population and generation. The experiment was repeated twice.

2.3 Physiological, biochemical, and profiling analysis of G2 populations

The experiments were carried out in a growth chamber (temperature 32/28 °C and photoperiod of 12 h). The experiments had three replications and consisted of 8 L pots filled with the same soil described in the 2.2 section. Factor A consisted of populations of G 2 generation: CHK and DRY×GLY (more and less sensitive glyphosate, respectively). Factor B: doses of 0 and 360 g a.e. ha-1 of glyphosate. The glyphosate dose used was the one that previously demonstrated plant injury without plant death, allowing the study of the physiological state of responses to the herbicide. The herbicide was applied when plants reached an average of 10 tillers. In experiment 1, plants were collected at 24, 48, 96, and 144 h after glyphosate application (HAA), and in experiment 2 were collected at 48, 96, and 192 HAA. Samples were kept in an ultra-freezer (-80 °C). Leaf samples (second and third fully expanded leaves) were collected and immediately frozen in liquid nitrogen and stored in an ultra-freezer (-80 °C).

2.3.1. Measurements

Stomatal conductance was measured with an LI-1600 steady-state porometer (Li-Cor Biosciences™ - Lincoln, NE USA). All plants were evaluated, taking measures in the adaxial leaf face.

H2O2content and lipid peroxidation measurement were performed using leaf tissues (0.250 g) to determine H2O2 content according to Fipke et al. (2022)Fipke MV, Feijó AR, Garcia NS, Heck T, Viana VE, Dayan FE et al. Transgenerational effect of drought stress and sub-lethal doses of quizalofop-p-ethyl: decreasing sensitivity to herbicide and biochemical adjustment in Eragrostis plana . Agriculture. 2022;12(3):396. Available from: https://doi.org/10.3390/agriculture12030396
https://doi.org/10.3390/agriculture12030...
; and malondialdehyde (MDA) determination was performed according to ( Fipke et al., 2022Fipke MV, Feijó AR, Garcia NS, Heck T, Viana VE, Dayan FE et al. Transgenerational effect of drought stress and sub-lethal doses of quizalofop-p-ethyl: decreasing sensitivity to herbicide and biochemical adjustment in Eragrostis plana . Agriculture. 2022;12(3):396. Available from: https://doi.org/10.3390/agriculture12030396
https://doi.org/10.3390/agriculture12030...
; Velikova et al., 2000Velikova V, Yordanov I, Edreva A. Oxidative stress and some antioxidant systems in acid rain-treated bean plants protective role of exogenous polyamines. Plant Sci. 2000;151(1):59-66. Available from: https://doi.org/10.1016/S0168-9452(99)00197-1
https://doi.org/10.1016/S0168-9452(99)00...
).

Antioxidant enzymes: were performed with leaf tissues (0.250 g) to extract and determine SOD (EC 1.15.1.1), CAT (EC 1.11.1.6), APX (EC 1.11.1.11) activity adapted by Fipke et al. (2022)Fipke MV, Feijó AR, Garcia NS, Heck T, Viana VE, Dayan FE et al. Transgenerational effect of drought stress and sub-lethal doses of quizalofop-p-ethyl: decreasing sensitivity to herbicide and biochemical adjustment in Eragrostis plana . Agriculture. 2022;12(3):396. Available from: https://doi.org/10.3390/agriculture12030396
https://doi.org/10.3390/agriculture12030...
.

2.3.2. Quantification of glyphosate, AMPA, and aromatic amino acids

The concentrations of glyphosate, AMPA, shikimic acid, phenylalanine, and tyrosine were determined in duplicate using High-Performance Liquid Chromatography (HPLC) coupled with Tandem Mass Spectrometry (MS/MS). First, the treated leaves were immersed in a 10-mL wash solution (water in pH 2.5 with phosphoric acid 6 M) and the unabsorbed glyphosate was removed from the leaf surface. Herbicide extraction was performed according to the methodology described by Gomes et al. (2015)Gomes GLGC, Carbonari CA, Velini ED, Trindade MLB, Silva JRM. Extraction and simultaneous determination of glyphosate, AMPA and compounds of the shikimic acid pathway in plants. Planta Daninha. 2015;33(2):295-304. Available from: https://doi.org/10.1590/0100-83582015000200015
https://doi.org/10.1590/0100-83582015000...
.

To quantify glyphosate and metabolites, HPLC-MS/MS system equipped with a mass spectrometer (Quadrupole-Orbitrap, Q-Exactive Focus, Thermo ScientificTM). For separation, a Hypercarb (5 µm 250 Å, 50 mm x 3 mm) column was used, and as mobile phase: (A) 5 mM of ammonium acetate in water (pH 7.0) and (B) 5 mM of ammonium acetate in methanol (pH 7.0), in gradient condition. Calibration curves were prepared in water employing analytical standards with purity ≥99% (Sigma-Aldrich, St. Louis, MO, USA). Method linear range (R2 ≥ 0.99) was 1 to 500 µg L-1 for all compounds, and samples were diluted when necessary.

2.4 Gene expression analysis

In this study we evaluated the expression of the target genes 5-Enolpyruvylshikimate-3-phosphate synthase gene ( EcEPSPS ) Aldo-keto reductase ( EcAKR ), ABC-C Family MRP10 ( EcM10 ), and ABC-C Family MRP8 ( EcM11 ). EPSPS was chosen because it is directly involved in synthesizing the enzyme inhibited by glyphosate ( Gaines et al., 2010Gaines TA, Zhang W, Wang D, Bukun B, Chisholm ST, Shaner DL et al. Gene amplification confers glyphosate resistance in Amaranthus palmeri . Proc Natl Acad Sci USA. 2010;107(3):1029-34. Available from: https://doi.org/10.1073/pnas.0906649107
https://doi.org/10.1073/pnas.0906649107...
). The AKR , M10 , M11 genes are known to be involved in the process of glyphosate metabolism ( Pan et al., 2019Pan L, Yu Q, Han H, Mao L, Nyporko A, Fan L et al. Aldo-keto reductase metabolizes glyphosate and confers glyphosate resistance in Echinochloa colona . Plant Physiol. 2019;181(4):1519-34. Available from: https://doi.org/10.1104/pp.19.00979
https://doi.org/10.1104/pp.19.00979...
; Piasecki et al., 2019aPiasecki C, Carvalho IR, Cechin J, Goulart FAP, Maia LC, Agostinetto D et al. Oxidative stress and differential antioxidant enzyme activity in glyphosate-resistant and-sensitive hairy fleabane in response to glyphosate treatment. Bragantia 2019a;78(3):379-96. Available from: https://doi.org/10.1590/1678-4499.20180289
https://doi.org/10.1590/1678-4499.201802...
).

Oligonucleotides for the target genes were obtained based on the Eragrostis curvula’s genome available in the National Center for Biotechnology Information – NCBI (National Center for Biotechnology Information, 2020) and were designed using Primer3Plus software ( Untergasser et al., 2007Untergasser A, Nijveen H, Rao X, Bisseling T, Geurts R, Leunissen JAM. Primer3Plus, an enhanced web interface to Primer3. Nucleic Acids Res. 2007;35(Web Server):71-4. Available from: https://doi.org/10.1093/nar/gkm306
https://doi.org/10.1093/nar/gkm306...
). In addition, rice ( Oryza sativa L.) reference genes were used ( Table 1 ). The GenBank IDs are EcAKR (TVU36522.1), EcEPSPS (AP014962.1), EcMRP10 (TVU07241.,1), and EcMRP8 (TVU00075.1).

Table 1
Oligonucleotides for the reference and target genes used to study the Eragrostis plana gene expression in RT-qPCR.

The total RNA was extracted from leaves (three biological replicates) using PureLink™ (Plant RNA Reagent-InvitrogenTM, Carlsbas, USA) following the manufacturer’s instructions. RNA concentration, quality and integrity were assessed using NanoVue™ (GE HealthcareTM, Buckinghamshire, UK) and agarose gel electrophoresis, respectively. RNA (1 µg) of each sample was treated with DNase I (Invitrogen) according to the manufacturer’s recommendations and converted into cDNA using oligo(dT) and the SuperScript™ III First-Strand Synthesis System kit (Invitrogen).

RT-qPCR experiments were performed in LightCycler™ 480 thermocycler with three biological and three technical replicates using oligonucleotides for the target and reference genes ( Table 1 ), according to the MIQE Guidelines ( Bustin et al., 2009Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M et al. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009;55(4):611-22. Available from: https://doi.org/10.1373/clinchem.2008.112797
https://doi.org/10.1373/clinchem.2008.11...
). The amplification efficiency and specificity of each oligonucleotide were determined in validation experiments using four dilutions of cDNA. The reactions were performed following the protocol described by Fipke et al. (2022)Fipke MV, Feijó AR, Garcia NS, Heck T, Viana VE, Dayan FE et al. Transgenerational effect of drought stress and sub-lethal doses of quizalofop-p-ethyl: decreasing sensitivity to herbicide and biochemical adjustment in Eragrostis plana . Agriculture. 2022;12(3):396. Available from: https://doi.org/10.3390/agriculture12030396
https://doi.org/10.3390/agriculture12030...
. Gene expression was calculated following the 2^(-delta-delta CT) method (Livak, Schmittgen, 2001), using CHK population without glyphosate (0 g a.e. ha-1) at each time of collection as a baseline and normalized with OsACT1 , Os18S, and OsEF1 α reference genes.

2.5 Statistics analysis

Data were tested for homogeneity of variance and normality. Visual injury and SDW data were fitted with a nonlinear sigmoid y=a/1+expxED50/b and log-logistic y=a/1+x/GR50b regressions model. The sensitivity reduction rate was calculated with GR50(DRY,GLYorDRY×GLY)/ED50 or GR50 (CHK) in the same generation.

The stomatal conductance, H2O2 content, lipid peroxidation, antioxidant enzyme activity, glyphosate, AMPA, shikimate acid, L-tyrosine, and L-phenylalanine concentration were subjected to analysis of variance, calculated confidence intervals (CI; 95%). For plants of the same population with and without herbicide, a t-test was performed. The statistical analysis was conducted using SAS University EditionTM (SAS Institute, Inc.TM, Cary, NC, USA) statistical program.

3.Results and Discussion

3.1 Dose-response curves for populations to check for sensitivity shift

In G 0 plants ( Figure 2A ), the dose of glyphosate required to cause 50% plant injury (ED50) was 173 g a.e. ha-1 ( Table 2 ). Regarding G 1 and G 2 generation, CHK populations showed an increased dose of glyphosate compared to G 0 , with ED50 229 and 278 g a.e. ha-1, respectively. In G 1 there were no differences in the ED50 of the DRY, GLY, and DRY×GLY populations compared to the CHK population. When analyzing the G 2 generation, the DRY×GLY population presented the highest ED50, 443 g a.e. ha-1, among populations. In G 2 generation, the sensitivity reduction rate ( Table 2 ) showed 1.59 in the DRY×GLY population.

Figure 2
Visual injury ( A , B and C ) and shoot dry weight ( D , E and F ) of Eragrostis plana populations from generations G 0 ( A and D ), G 1 ( B and E ), and G 2 ( C and F ) at 35 DAA of glyphosate doses. Parameters estimated for the curves are presented in Table 1. Vertical error bars represent a 95% confidence interval

Table 2
Parameters estimates for the curves for visual injury and shoot dry weight of Eragrostis plana populations of glyphosate doses.

In the dose of glyphosate required to cause 50% of shoot dry weight reduction (GR50), there were no differences between the CHK populations of the G 0 , G 1 , and G 2 generations, with 122, 199, and 152 g a.e. ha-1, respectively. In G 1 generations, there were no differences in GR50 between populations. However, in the G 2 generation, the DRY, GLY, and DRY×GLY populations presented higher GR50 (284, 308, and 436 g a.e. ha-1, respectively) than the CHK population. In G 2 generation, the sensitivity reduction rate demonstrated 1.85, 2.02, and 2.85 for DRY, GLY, and DRY×GLY populations, respectively.

These results indicated that the offspring of plants exposed to water stress decreased sensitivity to glyphosate, mainly in G 2 . Plants from the DRY×GLY population, which received both stresses (drought followed by a sub-lethal dose of glyphosate), were less sensitive to glyphosate. It is well known that plants exposed to previous stresses may have improved tolerance to other stresses, a phenomenon related to memory mechanisms in plants ( Galviz et al., 2020Galviz YCF, Ribeiro RV, Souza GM. Yes, plants do have memory. Theor Exp Plant Physiol. 2020;32:195-202. Available from: https://doi.org/10.1007/s40626-020-00181-y
https://doi.org/10.1007/s40626-020-00181...
). For instance, in another study, E. plana submitted to drought, quizalofop, and drought plus quizalofop had a lesser sensitivity to quizalofop than non-stressed plants ( Fipke et al., 2022Fipke MV, Feijó AR, Garcia NS, Heck T, Viana VE, Dayan FE et al. Transgenerational effect of drought stress and sub-lethal doses of quizalofop-p-ethyl: decreasing sensitivity to herbicide and biochemical adjustment in Eragrostis plana . Agriculture. 2022;12(3):396. Available from: https://doi.org/10.3390/agriculture12030396
https://doi.org/10.3390/agriculture12030...
). Consequently, the current study also sought to determine the possible biochemical mechanisms involved in priming for glyphosate reduced sensitivity in E. plana .

3.2 Physiological, biochemical, and metabolites analysis of G2 populations

3.2.1. Stomatal conductance

At 24 and 48 HAA ( Figures 3A and 3B ), all plants treated with glyphosate had lower stomatal conductance than untreated plants (without glyphosate). DRY×GLY population glyphosate-treated at 24, 48, and 144 HAA ( Figures 3A and 3B ) also presented reduced stomatal conductance compared with untreated plants.

Figure 3
Stomatal conductance (mmol m-2 s-1) of Eragrostis plana populations CHK ( A , no stress) and DRY×GLY ( B , drought plus glyphosate) with and without herbicide at 24, 48, 96, and 144 h after application of glyphosate. Error bars represent 95% confidence intervals. *t-test (α=0.05); **t-test (α=0.01)

In our study, the reduction in stomatal conductance in all populations was a typical response to glyphosate ( Figure 3 ). However, this response is not necessarily involved in harmful effects caused by the herbicide in plants. Lolium perenne resistant to glyphosate had a stomatal closure regulation; however, it was linked to carbon fixation optimization according to assimilated demands ( Yanniccari et al., 2012Yanniccari M, Tambussi E, Istilart C, Castro AM. Glyphosate effects on gas exchange and chlorophyll fluorescence responses of two Lolium perenne L. biotypes with differential herbicide sensitivity. Plant Physiol Biochem. 2012;57:210-7. Available from: https://doi.org/10.1016/j.plaphy.2012.05.027
https://doi.org/10.1016/j.plaphy.2012.05...
). Stomatal conductance can be used as a biomarker in plants sensitive to glyphosate ( Yanniccari et al., 2012Yanniccari M, Tambussi E, Istilart C, Castro AM. Glyphosate effects on gas exchange and chlorophyll fluorescence responses of two Lolium perenne L. biotypes with differential herbicide sensitivity. Plant Physiol Biochem. 2012;57:210-7. Available from: https://doi.org/10.1016/j.plaphy.2012.05.027
https://doi.org/10.1016/j.plaphy.2012.05...
).

3.2.2. H2O2 content and lipid peroxidation

Glyphosate-treated CHK population produced 82 and 80% higher hydrogen peroxide than untreated plants at 96 and 144 HAA ( Figure 4a ). At 144 HAA, the glyphosate-treated DRY×GLY population ( Figure 4b ) also showed higher hydrogen peroxide than untreated plants. Depending on the efficiency of redox homeostasis, it can be related to apparent antagonistic functions such as, on the one hand, signaling to prevent stressful situations and, on the other hand, causing oxidative stresses ( Moretti et al., 2017Moretti ML, Alárcon-Reverte R, Pearce S, Morran S, Hanson BD. Transcription of putative tonoplast transporters in response to glyphosate and paraquat stress in Conyza bonariensis and Conyza canadensis and selection of reference genes for qRT-PCR. Plos One. 2017;12(7):1-16. Available from: https://doi.org/10.1371/journal.pone.0180794
https://doi.org/10.1371/journal.pone.018...
).

Figure 4
Hydrogen peroxide content ( A and B ) and thiobarbituric acid-reactive substances ( C and D ) at 24, 48, 96, and 144 h of Eragrostis plana populations CHK ( A and C , no stress) and DRY×GLY ( B and D , drought plus glyphosate) with and without herbicide hours after application of glyphosate. Error bars represent a 95% confidence interval. *t-test (α=0.05); **t-test (α=0.01)

Inhibition of EPSPS by glyphosate indirectly induces ROS production ( Maroli et al., 2015Maroli AS, Nandula VK, Dayan FE, Duke SO, Gerard P, Tharayil N. Metabolic profiling and enzyme analyses indicate a potential role of antioxidant systems in complementing glyphosate resistance in an Amaranthus palmeri biotype. J Agric Food Chem. 2015;63(41):9199-209. Available from: https://doi.org/10.1021/acs.jafc.5b04223
https://doi.org/10.1021/acs.jafc.5b04223...
). This unregulated accumulation of H2O2 and other ROS causes lipid peroxidation and loss of membrane integrity, and, ultimately, cell death (Foyer, Noctor, 2005). The glyphosate-treated CHK population showed the highest lipid peroxidation levels ( Figure 4c ), with an increase of 49, 75, and 70% (48, 96, and 144 HAA, respectively). On the other hand, primed DRY×GLY plants ( Figure 4d ) had no ROS increases in the evaluated periods.

Accumulation of H2O2 and the concomitant increase of lipid peroxidation ( Figure 4 ) could explain the susceptibility of the CHK population to glyphosate. At 96 and 144 HAA, the high production of H2O2 likely triggered an increase in lipid peroxidation. At 48 HAA, although H2O2 apparently did not induce lipid peroxidation, it could be reached by other types of ROS responsible for this process. In C. bonariensis , susceptible plants showed 5.2 and 3.2-fold higher in H2O2 and lipid peroxidation than untreated plants (Piasecki et al., 2019a). Salix miyabeana plants subjected to glyphosate had an increase in the production of H2O2 after 6 HAA, reaching a peak at 72 HAA. Likely, from 6 HAA, a gradual increase in lipid peroxidation has occurred ( Gomes et al., 2017Gomes MP, Le Manach SG, Hénault-Ethier L, Labrecque M, Lucotte M, Juneau P. Glyphosate-dependent inhibition of photosynthesis in willow. Front Plant Sci. 2017;8:1-13. Available from: https://doi.org/10.3389/fpls.2017.00207
https://doi.org/10.3389/fpls.2017.00207...
).

3.2.3. Antioxidant enzymes

SOD activity differed only at 24 and 96 HAA in the glyphosate-treated DRY×GLY population compared to the untreated plants ( Figure 5b ). In offspring of plants DRY×GLY population, it was observed to have a higher baseline of SOD activity than untreated CHK population. All plants presented increased SOD activity following glyphosate application, with 16 and 56% for CHK and DRY×GLY, respectively, demonstrating that plants submitted to stresses in previous generations abet stimulated SOD activity. SOD is a crucial enzyme in regulating oxidative stress. For instance, C. bonariensis GLY resistant had 1.6-fold SOD activity than the susceptible biotype ( Piasecki et al., 2019aPiasecki C, Carvalho IR, Cechin J, Goulart FAP, Maia LC, Agostinetto D et al. Oxidative stress and differential antioxidant enzyme activity in glyphosate-resistant and-sensitive hairy fleabane in response to glyphosate treatment. Bragantia 2019a;78(3):379-96. Available from: https://doi.org/10.1590/1678-4499.20180289
https://doi.org/10.1590/1678-4499.201802...
).

CHK population exhibited differences in CAT activity at 24, 48, and 144 HAA ( Figure 5c ), showing increases of 81, 60, and 49% in glyphosate-treated plants than the untreated. However, this increase in CAT activity in the CHK population was like or less than the glyphosate-treated DRY×GLY population. In 144 HAA ( Figure 5d ), the glyphosate-treated DRY×GLY population presented the highest CAT activity. APX activity ( Figure 5f ) exhibited differences for the DRY×GLY population (24 and 96 HAA). CHK population did not differ from APX activity with or without quizalofop-treated ( Figure 5e ).

SOD, CAT, and APX enzymes play an essential role in ROS detoxification. In C. bonariensis , a glyphosate-resistant biotype demonstrated higher CAT and APX activity than a susceptible biotype when the herbicide was applied ( Piasecki et al., 2019aPiasecki C, Carvalho IR, Cechin J, Goulart FAP, Maia LC, Agostinetto D et al. Oxidative stress and differential antioxidant enzyme activity in glyphosate-resistant and-sensitive hairy fleabane in response to glyphosate treatment. Bragantia 2019a;78(3):379-96. Available from: https://doi.org/10.1590/1678-4499.20180289
https://doi.org/10.1590/1678-4499.201802...
). Application of glyphosate in S. miyabeana stimulated ROS production and, consequently, increased enzymatic activities of CAT and APX ( Gomes et al., 2017Gomes MP, Le Manach SG, Hénault-Ethier L, Labrecque M, Lucotte M, Juneau P. Glyphosate-dependent inhibition of photosynthesis in willow. Front Plant Sci. 2017;8:1-13. Available from: https://doi.org/10.3389/fpls.2017.00207
https://doi.org/10.3389/fpls.2017.00207...
). Antioxidant enzymes SOD, CAT, and APX are responsible for reducing oxidative stress and complement important resistance mechanisms enabling plants to tolerate herbicide glyphosate ( Maroli et al., 2015Maroli AS, Nandula VK, Dayan FE, Duke SO, Gerard P, Tharayil N. Metabolic profiling and enzyme analyses indicate a potential role of antioxidant systems in complementing glyphosate resistance in an Amaranthus palmeri biotype. J Agric Food Chem. 2015;63(41):9199-209. Available from: https://doi.org/10.1021/acs.jafc.5b04223
https://doi.org/10.1021/acs.jafc.5b04223...
).

3.2.4. Glyphosate, AMPA, shikimic acid, and amino acids concentration in plants

Glyphosate concentration in CHK and DRYxGLY populations ( Figure 6a ) reveals that in the CHK population, there was 72% higher glyphosate concentration at 48 HAA compared to the DRYxGLY population. There were no differences in glyphosate concentration regarding sampling times in CHK and DRYxGLY populations. In CHK and DRYxGLY populations, glyphosate concentration remained similar in the leaf tissues until 144 HAA, indicating low translocation and/or metabolization in the leaf tissues.

Figure 6
Glyphosate ( A ) and AMPA ( B ) concentration (µg g-1) of Eragrostis plana populations (CHK – no stress and DRYxGLY – drought plus glyphosate) at 24, 48, 96, and 144 h after applying glyphosate. Error bars represent a 95% confidence interval. *t-test (α=0.05); **t-test (α=0.01)

Glyphosate degradation in plants generates metabolites AMPA and glyoxylate ( Rojano-Delgado et al., 2012Rojano-Delgado AM, Cruz-Hipolito H, Prado R, Castro MDL, Franco AR. Limited uptake, translocation and enhanced metabolic degradation contribute to glyphosate tolerance in Mucuna pruriens var. utilis plants. Phytochemistry. 2012;73:34-41. Available from: https://doi.org/10.1016/j.phytochem.2011.09.007
https://doi.org/10.1016/j.phytochem.2011...
). Herein, AMPA metabolite was found in both populations CHK and DRYxGLY ( Figure 6b ). In the CHK population, AMPA was detected only at 24 and 48 HAA, with a concentration of 4.4 and 2.5 µg g-1. In the DRYxGLY population, AMPA was found in all evaluation periods, with the highest concentration at 24 HAA and the lowest at 144 HAA (3.1 and 0.8 µg g-1, respectively).

These results indicated that there was detoxification of glyphosate forming AMPA and that it was present in both more sensitive and less sensitive populations to glyphosate (CHK and DRYxGLY, respectively). However, AMPA levels represent 13.6 and 4.9% of glyphosate concentration at 24 and 48 HAA in the CHK population, while in the DRYxGLY population, AMPA levels were equivalent to 10.2, 8.8, 6.9, and 1.6% of the glyphosate concentration at 24, 48, 96 and 144 HAA, respectively.

EPSPS inhibition by glyphosate leads to rapid accumulation of shikimic acid and reduced synthesis of aromatic amino acids ( Rojano-Delgado et al., 2012Rojano-Delgado AM, Cruz-Hipolito H, Prado R, Castro MDL, Franco AR. Limited uptake, translocation and enhanced metabolic degradation contribute to glyphosate tolerance in Mucuna pruriens var. utilis plants. Phytochemistry. 2012;73:34-41. Available from: https://doi.org/10.1016/j.phytochem.2011.09.007
https://doi.org/10.1016/j.phytochem.2011...
). All glyphosate-treated plants had a greater accumulation of shikimic acid than plants without glyphosate ( Figure 7a ). At all evaluated periods, the CHK population showed a higher concentration of shikimic acid than DRYxGLY populations, exhibiting an increase of 291, 97, 35, and 195% at 24, 48, 96, and 144 HAA, respectively. These results explain the lower sensitivity of DRYxGLY compared to CHK population. Comparing both populations, when the plants were not treated with glyphosate, there were no differences in shikimic acid concentration.

Figure 7
Shikimate acid ( A ), L-tyrosine ( B ), and L-phenylalanine ( C ) concentration (μg g-1) of Eragrostis plana populations (CHK – no stress and DRYxGLY – drought plus glyphosate) at 24, 48, 96 and 144 h after applying glyphosate. Error bars represent a 95% confidence interval. t-test was performed between populations in same glyphosate dose; *t-test (α=0.05); **t-test (α=0.01)

In the glyphosate-treated CHK population, an increase in shikimic acid was observed after 24 HAA and later remained constant at 48 and 144 HAA. At 96 HAA, a reduction in shikimic acid was detected, consistent with the lower glyphosate concentration ( Figure 6a ) in the CHK population. In the glyphosate-treated DRYxGLY population, shikimic acid accumulation initiates at 24 HAA, reaching the highest concentration at 48 and 96 HAA and reducing to 144 HAA. Plants with faster and higher shikimic acid accumulation over time demonstrate substantial inhibition of EPSPS and, therefore, greater sensitivity to glyphosate ( Tani et al., 2015Tani E, Chachalis D, Travlos IS. A glyphosate resistance mechanism in Conyza canadensis involves synchronization of EPSPS and ABC-transporter genes. Plant Mol Biol Rep. 2015;33:1721-30. Available from: https://doi.org/10.1007/s11105-015-0868-8
https://doi.org/10.1007/s11105-015-0868-...
). Conyza canadensis plants susceptible to glyphosate showed higher shikimic acid accumulation than resistant plants ( Tani et al., 2015Tani E, Chachalis D, Travlos IS. A glyphosate resistance mechanism in Conyza canadensis involves synchronization of EPSPS and ABC-transporter genes. Plant Mol Biol Rep. 2015;33:1721-30. Available from: https://doi.org/10.1007/s11105-015-0868-8
https://doi.org/10.1007/s11105-015-0868-...
).

Tyrosine and phenylalanine are aromatic amino acids from the shikimic acid pathway ( Fernández-Escalada et al., 2017Fernández-Escalada M, Zulet-González A, Gil-Monreal M, Zabalza A, Ravet K, Gaines T, et al. Effects of EPSPS copy number variation (CNV) and glyphosate application on the aromatic and branched chain amino acid synthesis pathways in Amaranthus palmeri . Front Plant Sci. 2017;8:1-11. Available from: https://doi.org/10.3389/fpls.2017.01970
https://doi.org/10.3389/fpls.2017.01970...
). Between CHK and DRYxGLY, there was a difference in L-tyrosine concentration only for untreated plants ( Figure 7b ). In glyphosate-treated CHK population had increased L-tyrosine concentration compared to untreated plants at 48 and 96 HAA. However, at 144 HAA, there was an inversion, and glyphosate-treated plants presented a reduction of 11%. Regarding the DRYxGLY population, contrasts were found at 24 and 96 HAA, with a difference of -25% and 16% in glyphosate-treated plants compared to untreated plants.

The L-phenylalanine concentration ( Figure 7c ) varied between populations at different times and glyphosate application. In the glyphosate-treated CHK population, an increased L-phenylalanine concentration was observed from the initial evaluation at 24 HAA until 96 HAA and a reduction at 144 HAA ( Figure 7C ). From the glyphosate-treated DRYxGLY population, differences of -9, -42, -18, and 92% occurred compared to untreated plants (24, 48, 96, and 144 HAA, respectively).

Amino acid concentrations usually vary among plants; however, the increase in phenylalanine and tyrosine concentrations may be due to increased protein turnover following glyphosate treatment ( Fernández-Escalada et al., 2017Fernández-Escalada M, Zulet-González A, Gil-Monreal M, Zabalza A, Ravet K, Gaines T, et al. Effects of EPSPS copy number variation (CNV) and glyphosate application on the aromatic and branched chain amino acid synthesis pathways in Amaranthus palmeri . Front Plant Sci. 2017;8:1-11. Available from: https://doi.org/10.3389/fpls.2017.01970
https://doi.org/10.3389/fpls.2017.01970...
). In Amaranthus palmeri, susceptible to glyphosate, an increase in phenylalanine and tyrosine concentration at 72 HAA was reported. In contrast, in an A. palmeri population resistant to glyphosate, differences in concentration of these amino acids were not detected ( Fernández-Escalada et al., 2017Fernández-Escalada M, Zulet-González A, Gil-Monreal M, Zabalza A, Ravet K, Gaines T, et al. Effects of EPSPS copy number variation (CNV) and glyphosate application on the aromatic and branched chain amino acid synthesis pathways in Amaranthus palmeri . Front Plant Sci. 2017;8:1-11. Available from: https://doi.org/10.3389/fpls.2017.01970
https://doi.org/10.3389/fpls.2017.01970...
). Also, Lolium multiflorum plants susceptible to glyphosate showed a higher concentration of tyrosine than plants resistant to 72 HAA, while no changes in phenylalanine concentrations were observed ( Barroso et al., 2018Barroso AAM, Costa MGS, Neto NJ, Santos JI, Balbuena TS, Carbonari CA et al. Protein identification before and after glyphosate exposure in Lolium multiflorum genotypes. Pest Manag Sci. 2018;74(5):1125-33. Available from: https://doi.org/10.1002/ps.4831
https://doi.org/10.1002/ps.4831...
).

3.3 Gene expression

In this study, EPSPS expression was relative to the untreated CHK population ( Figure 8A ). At 48 HAA, EPSPS was downregulated in glyphosate-treated CHK and DRYxGLY populations. However, EPSPS was upregulated relative to the untreated CHK population in all other treatments and evaluation periods. At 96 HAA, the glyphosate-treated DRYxGLY population demonstrated the highest EPSPS expression, with 4-fold more transcripts than the untreated CHK population and 1.92-fold more than the untreated DRYxGLY population. At 192 HAA, glyphosate-treated CHK and untreated DRYxGLY populations showed the higher EPSPS expressions, with 2.3-fold and 2.6-fold, respectively, compared to the untreated CHK population.

Figure 8
Relative expression levels in log2 of EPSPS (A), AKR (B), MRP8 (C), and MRP10 (D) gene in leaves of Eragrostis plana populations CHK – with 360 g a.e. ha-1 of glyphosate (CHK- 360 g ha-1— dark bars), DRYxGLY – without herbicide (DRYxGLY- 0 g ha-1—dark gray bars) and DRYxGLY – 360 g a.e. ha-1 of glyphosate (DRYxGLY-360 g ha-1—gray bars) at 48, 96 and 192 h after application of glyphosate. Vertical error bars represent 95% confidence intervals

A mechanism of resistance to glyphosate that has been widely studied is the upregulation of the EPSPS gene. Increased copies of the EPSPS gene allows for a higher synthesis of the EPSPS enzyme, as there is more enzyme than glyphosate molecules for inhibition, making the plant resistant ( Gaines et al., 2010Gaines TA, Zhang W, Wang D, Bukun B, Chisholm ST, Shaner DL et al. Gene amplification confers glyphosate resistance in Amaranthus palmeri . Proc Natl Acad Sci USA. 2010;107(3):1029-34. Available from: https://doi.org/10.1073/pnas.0906649107
https://doi.org/10.1073/pnas.0906649107...
). Plants C. bonariensis and C. canadensis, both glyphosate-resistant, EPSPS expression was observed two times compared to susceptible plants. However, the authors report other NTSR mechanisms involved ( Dinelli et al., 2008Dinelli G, Marotti I, Bonetti A, Catizone P, Urbano JM, Barnes J. Physiological and molecular bases of glyphosate resistance in Conyza bonariensis biotypes from Spain. Weed Res. 2008;48(3):257-65. Available from: https://doi.org/10.1111/j.1365-3180.2008.00623.x
https://doi.org/10.1111/j.1365-3180.2008...
). In A. palmeri , it was found that increased expression of EPSPS was responsible for the resistance ( Gaines et al., 2010Gaines TA, Zhang W, Wang D, Bukun B, Chisholm ST, Shaner DL et al. Gene amplification confers glyphosate resistance in Amaranthus palmeri . Proc Natl Acad Sci USA. 2010;107(3):1029-34. Available from: https://doi.org/10.1073/pnas.0906649107
https://doi.org/10.1073/pnas.0906649107...
). The increase in EPSPS expression was explained by the amplification of genes, where several copies of the gene are present in the genome ( Gaines et al., 2010Gaines TA, Zhang W, Wang D, Bukun B, Chisholm ST, Shaner DL et al. Gene amplification confers glyphosate resistance in Amaranthus palmeri . Proc Natl Acad Sci USA. 2010;107(3):1029-34. Available from: https://doi.org/10.1073/pnas.0906649107
https://doi.org/10.1073/pnas.0906649107...
).

Studies have recently shown that plants that have submitted stress, such as herbicides, transmit pre-regulated expression patterns to the progeny (transgenerational) through epigenetic processes ( Délye et al., 2013Délye C, Jasieniuk M, Le Corre V. Deciphering the evolution of herbicide resistance in weeds. Trends Genet. 2013;29(11):649-58. Available from: https://doi.org/10.1016/j.tig.2013.06.001
https://doi.org/10.1016/j.tig.2013.06.00...
; Margaritopoulou et al., 2018Margaritopoulou T, Tani E, Chachalis D, Travlos I. Involvement of epigenetic mechanisms in herbicide resistance: the case of Conyza canadensis . Agriculture. 2018;8(1):1-9. Available from: https://doi.org/10.3390/agriculture8010017
https://doi.org/10.3390/agriculture80100...
). In C. canadensis resistant to glyphosate, there was an upregulation mechanism of EPSPS , and it was reported that this upregulation was due to the DNA methylation process affecting regions of EPSP-synthase 1 ( Margaritopoulou et al., 2018Margaritopoulou T, Tani E, Chachalis D, Travlos I. Involvement of epigenetic mechanisms in herbicide resistance: the case of Conyza canadensis . Agriculture. 2018;8(1):1-9. Available from: https://doi.org/10.3390/agriculture8010017
https://doi.org/10.3390/agriculture80100...
).

The AKR gene at 48 HAA was downregulated in all populations compared to the CHK population without herbicide ( Figure 8B ). The CHK population with herbicide did not differ in expression levels at any time. In contrast, the DRYxGLY population showed a relative increase of 1.1 and 0.6-fold at 96 HAA and 2.6 and 1.6-fold at 192 HAA compared to the CHK population without herbicide. DRYxGLY population had no difference in glyphosate application at any evaluation time. Aldo-keto reductase is a superfamily of enzymes responsible for catalyzing the reduction of NAD(P) H-dependent aldehydes and ketones ( Pan et al., 2019Pan L, Yu Q, Han H, Mao L, Nyporko A, Fan L et al. Aldo-keto reductase metabolizes glyphosate and confers glyphosate resistance in Echinochloa colona . Plant Physiol. 2019;181(4):1519-34. Available from: https://doi.org/10.1104/pp.19.00979
https://doi.org/10.1104/pp.19.00979...
). However, this superfamily is also involved in the metabolism of xenobiotics, secondary metabolism, and protection of osmolytes. It was recently reported that EcAKR4-1 is involved in the metabolization of glyphosate ( Pan et al., 2019Pan L, Yu Q, Han H, Mao L, Nyporko A, Fan L et al. Aldo-keto reductase metabolizes glyphosate and confers glyphosate resistance in Echinochloa colona . Plant Physiol. 2019;181(4):1519-34. Available from: https://doi.org/10.1104/pp.19.00979
https://doi.org/10.1104/pp.19.00979...
). However, it is reported that Echinochloa colona had resistance to glyphosate by Pro-106-Thr EPSPS target-site mutation ( McElroy, Hall, 2020McElroy JS, Hall ND. Echinochloa colona with reported resistance to glyphosate conferred by aldo-keto reductase also contains a pro-106-Thr EPSPS target site mutation. Plant Physiol. 2020;183(2):447-50. Available from: https://doi.org/10.1104/pp.20.00064
https://doi.org/10.1104/pp.20.00064...
).

These data suggest that the DRYxGLY population that has undergone two generations of drought stress and a sub-lethal dose of glyphosate was able to positively regulate AKR gene expression compared to the CHK population. A previous report demonstrates that E. colona could detoxify glyphosate, forming AMPA, by the AKR enzyme. The EcAKR4-1 gene was upregulated in this species ( Pan et al., 2019Pan L, Yu Q, Han H, Mao L, Nyporko A, Fan L et al. Aldo-keto reductase metabolizes glyphosate and confers glyphosate resistance in Echinochloa colona . Plant Physiol. 2019;181(4):1519-34. Available from: https://doi.org/10.1104/pp.19.00979
https://doi.org/10.1104/pp.19.00979...
).

ABC-transporter MRP8 ( M11 ) was differentially expressed among treatments ( Fig 8C ). The glyphosate-treated CHK and DRYxGLY populations were not very responsive to M11 , except for CHK at 96 HAA, where an M11 was upregulated. At 192 HAA, there was 1.7-fold M11 expression in the DRYxGLY population without herbicide. These results demonstrate that M11 was not responsive to glyphosate in E. plana and suggest that other ABC-transporter codifying gene may be involved in glyphosate responses in this species.

ABC-transporter MRP10 ( M10 ) was upregulated in almost all populations and times analyzed ( Figure 8D ). In the glyphosate-treated CHK population, M10 was upregulated in all periods analyzed, with 4.8-fold, 1.4-fold, and 2.8-fold at 48, 96, and 192 HAA, respectively. M10 was upregulated in the glyphosate-treated DRYxGLY population at 48 and 96 HAA and downregulation at 192 HAA. The higher difference between populations was detected at 96 HAA, where glyphosate-treated DRYxGLY showed the highest expression of M10 , 6.3-fold, thus being upregulated in this period. The untreated DRYxGLY population also showed upregulation of the M10 gene, demonstrating a greater baseline expression than the untreated CHK population.

M10 and M11 genes are ABC-transporters located in the tonoplast ( Moretti et al., 2017Moretti ML, Alárcon-Reverte R, Pearce S, Morran S, Hanson BD. Transcription of putative tonoplast transporters in response to glyphosate and paraquat stress in Conyza bonariensis and Conyza canadensis and selection of reference genes for qRT-PCR. Plos One. 2017;12(7):1-16. Available from: https://doi.org/10.1371/journal.pone.0180794
https://doi.org/10.1371/journal.pone.018...
). These transporters are related to weed resistance to glyphosate, acting on the active transport of this molecule to the vacuole ( Moretti et al., 2017Moretti ML, Alárcon-Reverte R, Pearce S, Morran S, Hanson BD. Transcription of putative tonoplast transporters in response to glyphosate and paraquat stress in Conyza bonariensis and Conyza canadensis and selection of reference genes for qRT-PCR. Plos One. 2017;12(7):1-16. Available from: https://doi.org/10.1371/journal.pone.0180794
https://doi.org/10.1371/journal.pone.018...
; Piasecki et al., 2019bPiasecki C, Yang Y, Benemann DP, Kremer FS, Galli V, Millwood RJ et al. Transcriptomic analysis identifies new non-target site glyphosate-resistance genes in Conyza bonariensis . Plants. 2019b;8(6):1-26. Available from: https://doi.org/10.3390/plants8060157
https://doi.org/10.3390/plants8060157...
). In C. canadensis and C. bonariensis, the ABC- transporters M10 and M11 are involved in resistance to glyphosate ( Piasecki et al., 2019bPiasecki C, Yang Y, Benemann DP, Kremer FS, Galli V, Millwood RJ et al. Transcriptomic analysis identifies new non-target site glyphosate-resistance genes in Conyza bonariensis . Plants. 2019b;8(6):1-26. Available from: https://doi.org/10.3390/plants8060157
https://doi.org/10.3390/plants8060157...
; Tani et al., 2015Tani E, Chachalis D, Travlos IS. A glyphosate resistance mechanism in Conyza canadensis involves synchronization of EPSPS and ABC-transporter genes. Plant Mol Biol Rep. 2015;33:1721-30. Available from: https://doi.org/10.1007/s11105-015-0868-8
https://doi.org/10.1007/s11105-015-0868-...
). In resistant C. canadensis , an upregulation of M10 and M11 was reported upon glyphosate treatment ( Nol et al., 2012Nol N, Tsikou D, Eid M, Livieratos IC, Giannopolitis CN. Shikimate leaf disc assay for early detection of glyphosate resistance in Conyza canadensis and relative transcript levels of EPSPS and ABC transporter genes. Weed Res. 2012;52(3):233-41. Available from: https://doi.org/10.1111/j.1365-3180.2012.00911.x
https://doi.org/10.1111/j.1365-3180.2012...
). The ABC-transporter M10 was 2.3-fold overexpressed in the resistant C. canadensis than the sensitive biotype ( Tani et al., 2015Tani E, Chachalis D, Travlos IS. A glyphosate resistance mechanism in Conyza canadensis involves synchronization of EPSPS and ABC-transporter genes. Plant Mol Biol Rep. 2015;33:1721-30. Available from: https://doi.org/10.1007/s11105-015-0868-8
https://doi.org/10.1007/s11105-015-0868-...
). The authors suggest that the ABC-transporters resistance in glyphosate-treated plants is mainly due to the exclusion of glyphosate in the vacuole and reduced herbicide translocation to the other plant tissues ( Nol et al., 2012Nol N, Tsikou D, Eid M, Livieratos IC, Giannopolitis CN. Shikimate leaf disc assay for early detection of glyphosate resistance in Conyza canadensis and relative transcript levels of EPSPS and ABC transporter genes. Weed Res. 2012;52(3):233-41. Available from: https://doi.org/10.1111/j.1365-3180.2012.00911.x
https://doi.org/10.1111/j.1365-3180.2012...
; Tani et al., 2015Tani E, Chachalis D, Travlos IS. A glyphosate resistance mechanism in Conyza canadensis involves synchronization of EPSPS and ABC-transporter genes. Plant Mol Biol Rep. 2015;33:1721-30. Available from: https://doi.org/10.1007/s11105-015-0868-8
https://doi.org/10.1007/s11105-015-0868-...
).

3.4 Enhanced mechanisms of response to glyphosate in E. plana

The ability of plants to adapt to environmental stresses allows regulating their metabolism to increase their tolerance ( Délye et al., 2013Délye C, Jasieniuk M, Le Corre V. Deciphering the evolution of herbicide resistance in weeds. Trends Genet. 2013;29(11):649-58. Available from: https://doi.org/10.1016/j.tig.2013.06.001
https://doi.org/10.1016/j.tig.2013.06.00...
). Mechanisms of resistance to herbicides that do not involve a site of action are the same as those used by plants in responding to abiotic stresses ( Délye et al., 2013Délye C, Jasieniuk M, Le Corre V. Deciphering the evolution of herbicide resistance in weeds. Trends Genet. 2013;29(11):649-58. Available from: https://doi.org/10.1016/j.tig.2013.06.001
https://doi.org/10.1016/j.tig.2013.06.00...
; Fipke et al., 2022Fipke MV, Feijó AR, Garcia NS, Heck T, Viana VE, Dayan FE et al. Transgenerational effect of drought stress and sub-lethal doses of quizalofop-p-ethyl: decreasing sensitivity to herbicide and biochemical adjustment in Eragrostis plana . Agriculture. 2022;12(3):396. Available from: https://doi.org/10.3390/agriculture12030396
https://doi.org/10.3390/agriculture12030...
). According to our results, E. plana G 1 and G 2 generations, when exposed to combined stress (drought plus a sub-lethal dose of glyphosate), allowed a decrease in sensitivity compared with plants that did not undergo any stress. This transgenerational effect can be caused by epigenetic processes, allowing differentiated regulation of the expression of stress tolerance mechanisms, such as antioxidant machinery ( Galviz et al., 2020Galviz YCF, Ribeiro RV, Souza GM. Yes, plants do have memory. Theor Exp Plant Physiol. 2020;32:195-202. Available from: https://doi.org/10.1007/s40626-020-00181-y
https://doi.org/10.1007/s40626-020-00181...
). However, there is a need for further studies to explain the epigenetic role in this observed effect.

In addition to the antioxidant processes presented, it is suggested that the reduced sensitivity of the DRYxGLY population to glyphosate is due to multiple mechanisms. It is important to note that a high accumulation of shikimic acid was detected in both populations. However, the DRYxGLY population had a lower shikimate acid concentration than the CHK population in all evaluation periods ( Figure 7a ). The increased expression of EPSPS in DRYxGLY at 96 HAA ( Figure 8a ) may be one factor in reducing shikimic acid concentration at 144 HAA ( Figure 7a ). The tyrosine and phenylalanine amino acids pool showed contrasting responses ( Figures 7b and 7c ).

Regarding the involved NTSR mechanisms, it was found that the concentration of glyphosate in the leaf tissue remained constant over time, which may demonstrate a reduction in translocation to other plant tissues ( Figure 6A ). Also, AMPA was detected due to glyphosate detoxification ( Figure 6B ), together with the increase in AKR expression ( Figure 6B ). However, AMPA concentration was relatively low to explain the difference in sensitivity to glyphosate. Also, an increased expression of M10 in DRYxGLY at 96 HAA ( Figure 8D ) may help exclude glyphosate in the vacuole and, consequently, reduce its translocation to other plant tissues.

4.Conclusions

Eragrostis plana exposed to combined stress (drought followed by a sub-lethal dose of glyphosate) showed induced transgenerational memory allowing acclimation to stressful conditions and decreasing sensitivity to glyphosate. The increased activity of the enzimes SOD, CAT, and APX following homeostatic adjustment in the DRYxGLY population, relieves glyphosate’s oxidative stress. Additionally, the set of mechanisms such as the possible metabolization of glyphosate in AMPA and the upregulation of AKR, and the upregulation of MRP10 (exclusion of glyphosate in the vacuole) seem to be involved in decreasing glyphosate sensitivity in Eragrostis plana DRYxGLY population.

Figure 5
Superoxide dismutase activity ( A and B ), catalase activity ( C and D ), and ascorbate peroxidase activity ( E and F ) at 24, 48, 96 and 144 h of Eragrostis plana populations CHK ( A , C and E , no stress) and DRY×GLY ( B , D and F drought plus glyphosate) with and without herbicide hours after application of quizalofop. Error bars represent a 95% confidence interval. *t-test (α=0.05); **t-test (α=0.01)

References

  • Barroso AAM, Costa MGS, Neto NJ, Santos JI, Balbuena TS, Carbonari CA et al. Protein identification before and after glyphosate exposure in Lolium multiflorum genotypes. Pest Manag Sci. 2018;74(5):1125-33. Available from: https://doi.org/10.1002/ps.4831
    » https://doi.org/10.1002/ps.4831
  • Bastiani MO, Roma-Burgos N, Langaro AC, Salas-Perez RA, Rouse CE, Fipke MV et al. Ammonium sulfate improves the efficacy of glyphosate on South African lovegrass ( Eragrostis plana ) under water stress. Weed Sci. 2021;69(2):167-76. Available from: https://doi.org/10.1017/wsc.2020.97
    » https://doi.org/10.1017/wsc.2020.97
  • Burns EE, Keith BK, Refai MY, Bothner B, Dyer WE. Constitutive redox and phosphoproteome changes in multiple herbicide resistant Avena fatua L. are similar to those of systemic acquired resistance and systemic acquired acclimation. J Plant Physiol. 2018;220:105-14. Available from: https://doi.org/10.1016/j.jplph.2017.11.004
    » https://doi.org/10.1016/j.jplph.2017.11.004
  • Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M et al. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009;55(4):611-22. Available from: https://doi.org/10.1373/clinchem.2008.112797
    » https://doi.org/10.1373/clinchem.2008.112797
  • Délye C, Jasieniuk M, Le Corre V. Deciphering the evolution of herbicide resistance in weeds. Trends Genet. 2013;29(11):649-58. Available from: https://doi.org/10.1016/j.tig.2013.06.001
    » https://doi.org/10.1016/j.tig.2013.06.001
  • Dinelli G, Marotti I, Bonetti A, Catizone P, Urbano JM, Barnes J. Physiological and molecular bases of glyphosate resistance in Conyza bonariensis biotypes from Spain. Weed Res. 2008;48(3):257-65. Available from: https://doi.org/10.1111/j.1365-3180.2008.00623.x
    » https://doi.org/10.1111/j.1365-3180.2008.00623.x
  • Dyer WE. Stress-induced evolution of herbicide resistance and related pleiotropic effects. Pest Manag Sci. 2018;74(8):1759-68. Available from: https://doi.org/10.1002/ps.5043
    » https://doi.org/10.1002/ps.5043
  • Fernández-Escalada M, Zulet-González A, Gil-Monreal M, Zabalza A, Ravet K, Gaines T, et al. Effects of EPSPS copy number variation (CNV) and glyphosate application on the aromatic and branched chain amino acid synthesis pathways in Amaranthus palmeri . Front Plant Sci. 2017;8:1-11. Available from: https://doi.org/10.3389/fpls.2017.01970
    » https://doi.org/10.3389/fpls.2017.01970
  • Fipke MV, Feijó AR, Garcia NS, Heck T, Viana VE, Dayan FE et al. Transgenerational effect of drought stress and sub-lethal doses of quizalofop-p-ethyl: decreasing sensitivity to herbicide and biochemical adjustment in Eragrostis plana . Agriculture. 2022;12(3):396. Available from: https://doi.org/10.3390/agriculture12030396
    » https://doi.org/10.3390/agriculture12030396
  • Foyer CH, Noctor G. Redox homeostasis and antioxidant signaling: a metabolic interface between stress perception and physiological responses. Plant Cell. 2005;17(7):1866-75. Available from: https://doi.org/10.1105/tpc.105.033589
    » https://doi.org/10.1105/tpc.105.033589
  • Gaines TA, Duke SO, Morran S, Rigon CAG, Tranel PJ, Küpper A et al. Mechanisms of evolved herbicide resistance. J Biol Chem. 2020;295(30):10307-30. Available from: https://doi.org/10.1074/jbc.REV120.013572
    » https://doi.org/10.1074/jbc.REV120.013572
  • Gaines TA, Zhang W, Wang D, Bukun B, Chisholm ST, Shaner DL et al. Gene amplification confers glyphosate resistance in Amaranthus palmeri . Proc Natl Acad Sci USA. 2010;107(3):1029-34. Available from: https://doi.org/10.1073/pnas.0906649107
    » https://doi.org/10.1073/pnas.0906649107
  • Galviz YCF, Ribeiro RV, Souza GM. Yes, plants do have memory. Theor Exp Plant Physiol. 2020;32:195-202. Available from: https://doi.org/10.1007/s40626-020-00181-y
    » https://doi.org/10.1007/s40626-020-00181-y
  • Gomes GLGC, Carbonari CA, Velini ED, Trindade MLB, Silva JRM. Extraction and simultaneous determination of glyphosate, AMPA and compounds of the shikimic acid pathway in plants. Planta Daninha. 2015;33(2):295-304. Available from: https://doi.org/10.1590/0100-83582015000200015
    » https://doi.org/10.1590/0100-83582015000200015
  • Gomes MP, Le Manach SG, Hénault-Ethier L, Labrecque M, Lucotte M, Juneau P. Glyphosate-dependent inhibition of photosynthesis in willow. Front Plant Sci. 2017;8:1-13. Available from: https://doi.org/10.3389/fpls.2017.00207
    » https://doi.org/10.3389/fpls.2017.00207
  • Jain M, Nijhawan A, Tyagi AK, Khurana JP. Validation of housekeeping genes as internal control for studying gene expression in rice by quantitative real-time PCR. Biochem Biophys Res Comm. 2006;345(2):646-51. Available from: https://doi.org/10.1016/j.bbrc.2006.04.140
    » https://doi.org/10.1016/j.bbrc.2006.04.140
  • Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods. 2001;25(4):402-8. Available from: https://doi.org/10.1006/meth.2001.1262
    » https://doi.org/10.1006/meth.2001.1262
  • Margaritopoulou T, Tani E, Chachalis D, Travlos I. Involvement of epigenetic mechanisms in herbicide resistance: the case of Conyza canadensis . Agriculture. 2018;8(1):1-9. Available from: https://doi.org/10.3390/agriculture8010017
    » https://doi.org/10.3390/agriculture8010017
  • Maroli AS, Nandula VK, Dayan FE, Duke SO, Gerard P, Tharayil N. Metabolic profiling and enzyme analyses indicate a potential role of antioxidant systems in complementing glyphosate resistance in an Amaranthus palmeri biotype. J Agric Food Chem. 2015;63(41):9199-209. Available from: https://doi.org/10.1021/acs.jafc.5b04223
    » https://doi.org/10.1021/acs.jafc.5b04223
  • McElroy JS, Hall ND. Echinochloa colona with reported resistance to glyphosate conferred by aldo-keto reductase also contains a pro-106-Thr EPSPS target site mutation. Plant Physiol. 2020;183(2):447-50. Available from: https://doi.org/10.1104/pp.20.00064
    » https://doi.org/10.1104/pp.20.00064
  • Medeiros RB, Focht T. [Invasion, prevention, control and utilization of capim-annoni-2(Eragrostis plana Nees) in Rio Grande do Sul, Brazil]. Pesq Agrop Gaúcha. 2007;13(1-2):105-14. Portuguese.
  • Moretti ML, Alárcon-Reverte R, Pearce S, Morran S, Hanson BD. Transcription of putative tonoplast transporters in response to glyphosate and paraquat stress in Conyza bonariensis and Conyza canadensis and selection of reference genes for qRT-PCR. Plos One. 2017;12(7):1-16. Available from: https://doi.org/10.1371/journal.pone.0180794
    » https://doi.org/10.1371/journal.pone.0180794
  • National Center for Biotechnology Information – NCBI. Welcome to NCBI. Washington: National Center for Biotechnology Information; 2020[access Jan 10, 2020]. Available from: www.ncbi.nlm.nih.gov/
    » www.ncbi.nlm.nih.gov/
  • Nol N, Tsikou D, Eid M, Livieratos IC, Giannopolitis CN. Shikimate leaf disc assay for early detection of glyphosate resistance in Conyza canadensis and relative transcript levels of EPSPS and ABC transporter genes. Weed Res. 2012;52(3):233-41. Available from: https://doi.org/10.1111/j.1365-3180.2012.00911.x
    » https://doi.org/10.1111/j.1365-3180.2012.00911.x
  • Pan L, Yu Q, Han H, Mao L, Nyporko A, Fan L et al. Aldo-keto reductase metabolizes glyphosate and confers glyphosate resistance in Echinochloa colona . Plant Physiol. 2019;181(4):1519-34. Available from: https://doi.org/10.1104/pp.19.00979
    » https://doi.org/10.1104/pp.19.00979
  • Piasecki C, Carvalho IR, Cechin J, Goulart FAP, Maia LC, Agostinetto D et al. Oxidative stress and differential antioxidant enzyme activity in glyphosate-resistant and-sensitive hairy fleabane in response to glyphosate treatment. Bragantia 2019a;78(3):379-96. Available from: https://doi.org/10.1590/1678-4499.20180289
    » https://doi.org/10.1590/1678-4499.20180289
  • Piasecki C, Yang Y, Benemann DP, Kremer FS, Galli V, Millwood RJ et al. Transcriptomic analysis identifies new non-target site glyphosate-resistance genes in Conyza bonariensis . Plants. 2019b;8(6):1-26. Available from: https://doi.org/10.3390/plants8060157
    » https://doi.org/10.3390/plants8060157
  • Rojano-Delgado AM, Cruz-Hipolito H, Prado R, Castro MDL, Franco AR. Limited uptake, translocation and enhanced metabolic degradation contribute to glyphosate tolerance in Mucuna pruriens var. utilis plants. Phytochemistry. 2012;73:34-41. Available from: https://doi.org/10.1016/j.phytochem.2011.09.007
    » https://doi.org/10.1016/j.phytochem.2011.09.007
  • Tani E, Chachalis D, Travlos IS. A glyphosate resistance mechanism in Conyza canadensis involves synchronization of EPSPS and ABC-transporter genes. Plant Mol Biol Rep. 2015;33:1721-30. Available from: https://doi.org/10.1007/s11105-015-0868-8
    » https://doi.org/10.1007/s11105-015-0868-8
  • Untergasser A, Nijveen H, Rao X, Bisseling T, Geurts R, Leunissen JAM. Primer3Plus, an enhanced web interface to Primer3. Nucleic Acids Res. 2007;35(Web Server):71-4. Available from: https://doi.org/10.1093/nar/gkm306
    » https://doi.org/10.1093/nar/gkm306
  • Velikova V, Yordanov I, Edreva A. Oxidative stress and some antioxidant systems in acid rain-treated bean plants protective role of exogenous polyamines. Plant Sci. 2000;151(1):59-66. Available from: https://doi.org/10.1016/S0168-9452(99)00197-1
    » https://doi.org/10.1016/S0168-9452(99)00197-1
  • Yanniccari M, Tambussi E, Istilart C, Castro AM. Glyphosate effects on gas exchange and chlorophyll fluorescence responses of two Lolium perenne L. biotypes with differential herbicide sensitivity. Plant Physiol Biochem. 2012;57:210-7. Available from: https://doi.org/10.1016/j.plaphy.2012.05.027
    » https://doi.org/10.1016/j.plaphy.2012.05.027
  • Yuan JS, Abercrombie LLG, Cao Y, Halfhill MD, Zhou X, Peng Y et al. Functional genomics analysis of horseweed ( Conyza canadensis ) with special reference to the evolution of non–target-site glyphosate resistance. Weed Sci. 2010;58(2):109-17. Available from: https://doi.org/10.1614/ws-d-09-00037.1
    » https://doi.org/10.1614/ws-d-09-00037.1
  • Zhou J, Tao B, Messersmith CG, Nalewaja JD. Glyphosate efficacy on velvetleaf ( Abutilon theophrasti ) is affected by stress. Weed Sci. 2007;55(3):240-4. Available from: https://doi.org/10.1614/ws-06-173.1
    » https://doi.org/10.1614/ws-06-173.1
  • Zhou Y, Lu D, Li C, Luo J, Zhu BF, Zhu J et al. Genetic control of seed shattering in rice by the APETALA2 transcription factor Shattering Abortion1 . Plant Cell. 2012;24(3):1034-48. Available from: https://doi.org/10.1105/tpc.111.094383
    » https://doi.org/10.1105/tpc.111.094383
  • Funding: This research was funded by CAPES - financing code 001 for Marcus Vinícius Fipke and Vívian Ebeling Viana by CAPES/PNPD; by CNPq for the Research Fellowship of Luis Antonio de Avila/N.Proc. 310830/2019-2 and for the Universal grant Universal 2018 - CNPq N.Proc. 426714/2018-0; by FAPERGS - Program for the internationalization of Graduate Programs - FAPERGS/CAPES 06/2018: 19/2551-0000715-9.

Edited by

Approved by:
Editor in Chief: Carlos Eduardo Schaedler
Associate Editor: Lauren Lazaro

Publication Dates

  • Publication in this collection
    04 Nov 2022
  • Date of issue
    2022

History

  • Received
    1 June 2022
  • Accepted
    15 Sept 2022
Sociedade Brasileira da Ciência das Plantas Daninhas - SBCPD Rua Santa Catarina, 50, sala 1302 , 86010-470 - Londrina - Paraná / Brasil , +55 (51) 3308-6006 - Londrina - PR - Brazil
E-mail: sbcpd@sbcpd.org