Resistance to Fomesafen, Imazamox and Glyphosate in Euphorbia heterophylla from Brazil

: Euphorbia heterophylla is a species of weed that was previously controlled by fomesafen, imazamox and glyphosate, but continued use of these herbicides has selected resistant populations from the Rio Grande do Sul (Brazil). One resistant (R) strain and one susceptible (S) strain to fomesafen, imazamox and glyphosate were compared, the latter by recurrent selection. Dose-response tests showed multiple resistance to these herbicides. The required imazamox concentration to inhibit ALS by 50% was approximately 16 times greater in the R population than in the S population. Based on the EPSPS activity results, the R population was 10 fold less sensitive to glyphosate than the S counterpart. In addition, basal EPSPS activity from R plants was 3.3 fold higher than the level detected on S plants. The Proto IX assays showed high resistance to fomesafen in the R population that accumulated less Proto IX than the S population. Malathion assays showed the participation of CytP450 in fomesafen resistance, but a molecular mechanism could also be involved. To our knowledge, this is the ﬁrst characterisation of multiple resistance to these three groups of herbicides in E. heterophylla in the world.


Introduction
Euphorbia heterophylla (Wild poinsettia) is a summer annual plant with a short life cycle and with two or more generations per year. This species is a very common weed in South America that was later extended to Mexico and the Southern United States [1]. In soybean [Glycine max (L.) Merr.], E. heterophylla can cause a daily yield loss of up to 5.1 kg ha −1 , depending on the weed density [2]. Due to a high level of cross-pollination and elevated genetic recombination, there is high variability among populations of E. heterophylla which contributes to a rapid evolution of herbicide resistance [3,4].
Soybean is the most important crop in Brazil [5], and in 2018, occupied more than 30 million hectares corresponding to 45% of the cultivated area (~35 million ha) in the country [6]. In the 1980s, the control of E. heterophylla was achieved with protoporphyrinogen oxidase (PPO) inhibitor herbicides (fomesafen) and later in the mid-1980s, acetolactate synthase (ALS) inhibitor herbicides (mainly imidazolinones) began to be used [6]. However, in the early 2000s, biotypes of this species had already been selected for multiple resistance to fomesafen and imazamox [7]. On the other hand, the rapid adoption of glyphosate-resistant (GR) soybeans led to the exclusive use of glyphosate for the control E. heterophylla in the first years after the introduction of this technology. GR soybean varieties were officially introduced in Brazil in 2005; however, they began to be cultivated irregularly since 2000

Herbicides
The commercial herbicides used are shown in Table 1.

Plant Material
Euphorbia heterophylla seeds resistant (R) and susceptible (S) to glyphosate were supplied by Prof. Ribas Vidal [5]. The R population was collected in Porto Alegre, the Rio Grande do Sul (Brazil). This population was confirmed to be resistant to fomesafen, imazamox, and moderately resistant to glyphosate [5,12]. The S population was confirmed as glyphosate-susceptible, and their seeds were collected in places never treated with herbicides over the last 15 years [5]. S seeds were stored at 4 • C, and each year their germination level was texted, maintaining a percentage greater than 80%.

Evolution of Glyphosate-Resistance under Recurrent Selection
During the period from 2011 to 2019, the R population was subjected to nine cycles of recurrent glyphosate selection under field conditions. Seeds of the R population were germinated in Petri dishes containing filter paper moistened with distilled water. The Petri dishes were placed in a growth chamber at 28/18 • C (day/night) with a photoperiod of 16 h, 850 µmol m −2 s −1 of photosynthetic photon flux and 60% relative humidity. Four days later, 100 seedlings were transplanted into small plots (experimental unit of 2 × 6 m 2 with plants spaced at 35 cm (approximately) in a square planting pattern) in the field of experimentation at the University of Córdoba (Spain); the plants were watered three days by week during the experiment. Plants were treated at the 4-leaf growth stage using a Pulverex backpack sprayer with a T coupling for the wand equipped with four flat fan nozzles at a spraying pressure of 200 kPa and calibrated to give a 200 L ha −1 volume. Glyphosate was applied at 540, 720, 900, 1080 and 1800 g ae ha −1 ; after 5 months surviving plants finished their cycle and the seeds were harvested and used to repeat the selection process in successive years. The seeds collected from each dose of herbicide were cleaned and stored in a chamber at 4 • C until use. This selection process was repeated every year until 2019 in spring/summer, and during this last year a sequential application was made with fomesafen (plants with 3 leaves), imazamox (plants before tillering) and glyphosate (plants in tillering) at field dose ( Table 1). The percentage of surviving plants was determined each year to evaluate the evolution of resistance to glyphosate and finally multiple resistance to fomesafen and imazamox. The confirmation of the multiple levels of resistance as well as the putative resistance mechanisms were carried out with the seeds obtained in 2019. The experiment of recurrent selection was carried out employing three replicates for each glyphosate-dose and year.

Dose-Response Assays
Plants from R and S populations were sprayed with the following increasing doses of glyphosate (0, X/16, X/8, X/4, X/2, X, 2X, 4X and 6X g ae ha −1 ), imazamox (0, X/4, X/2, X, 4X, 8X, 10X and 20X g ai ha −1 ) and fomesafen (0, X/4, X/2, X, 4X, 10X, 20X and 40X g ai ha −1 ), with X being the field doses of every herbicide (Table 1). Five pots with one plant per pot were treated for each dose of herbicide. Spraying was carried out in a laboratory chamber equipped with an 8001 flat fan nozzle delivering 200 L ha −1 at a constant pressure of 200 kPa. Twenty-one days after treatment (DAT), plant mortality (LD 50 ) and dry weight (g) per plant (GR 50 ) were evaluated [26]. Data were expressed as percentages relative to the untreated control. The experiment was repeated twice at different times.

ALS Enzyme Activity
The ALS activity in response to imazamox was determined in vitro following the protocol described by Rojano-Delgado et al. [15]. Three grams of young leaf tissue of each biotype were ground separately using liquid N 2 . Polyvinylpyrrolidone (0.5 g) was added to the fine powder as well as extraction buffer [1M K-phosphate buffer solution (pH 7.5), 10 mM sodium pyruvate, 5 mM MgCl 2 , 50 mM thiamine pyrophosphate, 100 µflavin adenine dinucleotide, 12 mM dithiothreitol and glycerol-water (1:9, v/v)] in a proportion of 1:3 tissue:buffer. Samples were agitated for 10 min at 4 • C, filtered through cheesecloth and centrifuged (20,000× g for 20 min). The obtained supernatant was immediately used for the enzyme assays. For the ALS enzyme activity, 90 µL of enzyme extract was used with 110 µL of freshly prepared assay buffer [0.08 M K-phosphate buffer solution (pH 7.5), 0.5 M sodium pyruvate, 0.1 M MgCl 2 , 0.5 mM thiamine pyrophosphate and 1000 µM flavin adenine dinucleotide]. The imazamox concentrations assayed were 0, 0.1, 1, 10, 100, 1000 and 10,000 µM. Aliquots of 250 µL of a solution 0.04 M K 2 HPO 4 at pH 7.0 were added. This mixture was incubated for 1 h at 37 • C. Afterwards, the reaction was stopped by adding 50 µL of H 2 SO 4 /water 1:10 (v/v). To decarboxylate acetolactate to acetoin, the tubes with the mixture were heated for 15 min at 60 • C. A red complex (λ 520 nm) formed after the addition of 250 µL of creatine (5 g L −1 freshly prepared in water) and 250 µL of 1-naphthol (50 g L −1 freshly prepared in 5 N NaOH) before incubation at 60 • C for 15 min. Total protein content was determined by the Bradford method [27]. The experiment was carried out twice and with three replicates for each concentration of the herbicide.

EPSPS Enzyme Activity
The extraction and activity of the EPSPS enzyme were carried out following the methodology described by Palma-Bautista et al. [26]. Five grams of young leaf tissue from the R and S populations of E. heterophylla were ground with liquid N 2 to a fine powder. Samples were transferred to tubes containing 25 mL of extraction buffer (100 mM MOPS, 5 mM EDTA, 10% glycerol, 50 mM KCl, and 0.5 mM benzamidine) with 70 µL of ß-mercaptoethanol and 1% in polyvinylpolypyrrolidone (PVPP) to extract the EPSPS. Samples were vortexed for 5 min, avoiding foaming, and then centrifuged (18,000× g, 30 min, 4 • C). The supernatant was filtered using a cheesecloth in a cold beaker and the (NH 4 ) 2 SO 4 was slowly added while under continuous stirring (30 min) and at the end centrifuged (15,000× g, 30 min, 4 • C). This step was repeated once more to obtain a 70% (NH 4 ) 2 SO 4 (w/v) solution to precipitate the fraction that contained the EPSPS activity. The pellets obtained in the successive centrifugations were mixed and re-suspended in 1 mL assay buffer (100 mM MOPS, 1 mM MgCl 2 , 10% glycerol (v/v), 2 mM sodium molybdate, 200 mM NaF). These samples were dialysed using Slide-A-Lyzer dialysis cassettes (1000-MWC, Thermo Scientific, Meridian, IL, USA) overnight in 2 L of dialysis buffer (100 mM MOPS and 5 mM EDTA) at 4 • C on a stir plate. The pH of the buffers was adjusted to 7.0.
The concentration of total protein soluble (TPS) was determined by Bradford assay [27]. The EPSPS enzyme activity was determined in a continuous assay quantifying the inorganic phosphate (Pi) released from shikimate-3-phosphate with the EnzCheck phosphate assay kit (Invitrogen, Carlsbad, CA, USA) following the manufacturer's instructions. The glyphosate concentrations used were 0, 0.01, 0.1, 1, 10, 100 and 1000 µM. The amount Pi released was measured at 360 nm for 10 min in a spectrophotometer (DU-640, Beckman Coulter Inc., Fullerton, CA, USA). The EPSPS activity was calculated by determining the amount of Pi (µmol) released in µg of TSP −1 min −1 . EPSPS enzyme activity was expressed as a percentage of enzyme activity in the presence of glyphosate relative to the control (basal activity without glyphosate). Three replicates for each concentration in a completely randomised design were tested and the experiment was conducted twice.

Without CytP450 Inhibitor
Determination of the Proto IX levels without a CytP450 inhibitor was carried out following the methodology described by Fernández-Moreno et al. [28].
Leaf disks (approx 0.1 g each) from each population were incubated individually in a 6 cm diameter Petri dish containing 6 mL of 2% (wt/v) sucrose, 1 mM of 2-(N-morpholine) ethanesulphonic acid, and different concentrations (100, 500, 1000, 2000 and 5000 µM) of fomesafen at technical grade during 20 h at 25 • C in darkness. The same experiment was conducted but without herbicide as control. All extractions of Proto IX from leaf discs were made under a dim, green light source. After incubation, the leaf disks were homogenised using 2 mL of a mix of methanol and NH 4 OH 1N (9:1) at maximum Agronomy 2020, 10, 1573 5 of 14 speed for 15 s and kept in the dark for 10 min. Homogenates were centrifuged at 3000× g for 3 min, and then supernatants were collected. Pellets were re-suspended in 1 mL of methanol in the dark for 10 min, collecting the supernatant. Supernatants from both centrifugations were combined and dried under an air stream. Concentrates were reconstituted in 1 mL of the methanol and filtered through a 0.2-µm nylon syringe filter to remove particles. The extracts were stored in opaque glass vials at −20 • C until their analysis by the HPLC method.
Fifty microliters of the reconstituted phase were injected into a Gold HPLC System Beckman Coulter 126 (Fullerton, CA, USA) equipped with a Jasco FP-1520 fluorescence detector (Easton, PA, USA), with excitation and emission wavelengths set at 400 and 630 nm, respectively, and a 25 cm × 4.6 mm, 5-micron particle size Phenomenex SphereClone ® ODS column (Torrance, CA, USA). Deionised water was used as mobile phase A, and pure LC grade methanol as mobile phase B. The elution programme started by 60% mobile phase B and followed the linear gradient until 100% methanol in 10 min. After 30 min, the solvent system was returned to the initial conditions. The constant flow rate and column temperature were 1.0 mL/min and 25 • C, respectively. Proto levels in the extracts were quantified using a calibration curve obtained with a commercially available Proto IX standard from Sigma Aldrich (St. Louis, MO, USA). Data were expressed as nmol g −1 of fresh weight. All treatments were carried out three times using three replicates with and without herbicide in a completely randomised design.

With CytP450 Inhibitor
For this assay, we used whole plants at the four to the six-leaf stage, that were sprayed without and with malathion (2000 g ai ha −1 ). After approximately 120 min, the fomesafen was applied at 375 g ai ha −1 to the populations. The doses were based on previous studies [29,30]. Non-fomesafen incubation was used because the metabolism is a slower mechanism.
Plants not treated with fomesafen were used as controls. Ten plants from each population and fomesafen and malathion dose were harvested at 96 HAT using leaf discs (approximately 0.1 g each).
All Proto IX extractions from the leaf discs were carried out following the methodology described by Fernández-Moreno et al. [28]. Proto IX levels in the extracts were also quantified using a calibration curve obtained with a commercially available Proto IX standard from Sigma Aldrich (St. Louis, MO, USA). All treatments were carried out three times using three replicates with and without herbicide in a completely randomised design.

Statistical Analysis
The herbicide concentration causing a 50% reduction in plant growth (GR 50 ), the dose required to kill 50% of the population (LD 50 ) and the herbicide concentration causing a 50% inhibition of enzyme activity (I 50 ) were calculated by a non-linear regression analysis adjusted to a three-parameter model using the R package drc (R Core Team) (Equation (1)): where d is the coefficient corresponding to the upper asymptote, the coefficient b is the slope at the inflection point, and the herbicide concentration required to inhibit shoot growth or enzyme activity by 50% (i.e., GR 50 , LD 50 or I 50 , respectively), and x is the herbicide dose. The resistance index (RI = R/S) was computed as R-to-S GR 50 , GR 50 or I 50 ratios. Data on the field evolution of glyphosate, basal EPSPS activity and the proto IX were subjected to ANOVA using Statistix 10.0 (Analytical Software, Tallahassee, FA, USA). Percentage data were transformed (arcsine of the square root) to meet model assumptions of normality of the error distribution and variance homogeneity. Model assumptions were graphically inspected. When needed, differences between means were separated using Tukey's honestly significant difference test. The replicates of the experiments were pooled due to the lack of statistical differences between them.

Evolution of Glyphosate Resistance under Recurrent Selection
The selection pressure exerted by glyphosate on a population of E. heterophylla was rapid, and after 9 years of recurrent selection, the maximum dose in which around 90% of plants survived to the herbicide increased from 540 to 1800 g ae ha −1 ( Table 2). The original R2011 population survived 540 g ae ha −1 , while only 20.67% of plants survived from that population treated with 720 glyphosate (field doses in Brazil), but in 2019 a survival rate of 98.00% was obtained as a product of the recurrent selection (Table 2). Similarly, plants treated independently with 900, 1080 and 1800 g ae ha −1 of glyphosate had survival rates of 94, 96 and 89% after 7, 5 and 3 years of application, respectively ( Table 2).

Dose-Response Assays
The LD 50 values were far less for the S population, than the ones that survived to higher doses used in soybean cultivation in southern Brazil (Table 3 and Figure 1). Based on the GR 50 values of the R population, we calculated resistance factors (RF: GR 50 of R/GR 50 of S) from 13.6, 19.1 and 68.1 for glyphosate, imazamox and fomesafen (Table 3 and Figure 1). Table 3. Parameters of the log-logistic equation used (Equation (1)) to calculate the fomesafen, imazamox and glyphosate rates required for 50% reduction dry weight (GR 50 ), and % survival (LD 50 ) expressed as a percentage of the mean untreated control of the R and S E. heterophylla populations.

EPSPS Enzyme Activity
Based on the EPSPS activity data, a log-logistic model was fitted for each population (Figure 2A). From there, the required glyphosate concentration to inhibit EPSPS by 50% (I 50 ) was estimated, which was approximately 10 times greater in the R population than in the S population (Table 4). Table 4. Parameters of the log-logistic model (Equation (1)) of the inhibition of EPSPS enzymatic activity in response to glyphosate concentration for R and S E. heterophylla populations.

ALS Enzyme Activity
The quantity of imazamox needed to inhibit the ALS activity by 50% (I50) in R plants was 651.13 μM, while that in the S population was 40.64 μM (Table 5 and Figure 3A). In other words, the R population was 1602 fold more resistant to imazamox than the S counterpart (Table 5). No differences were found in the basal enzyme activities ( Figure 3B). Table 5. Parameters of the log-logistic model (Equation (1)) of the inhibition of ALS enzymatic activity in response to imazamox concentration for R and S E. heterophylla populations.   Apparent differences were found between S and R population in the enzyme basal activities that were 56.20 ± 2.47 and 189.51 ± 4.48 nmol Pi µg TSP −1 min −1 , respectively ( Figure 2B).

ALS Enzyme Activity
The quantity of imazamox needed to inhibit the ALS activity by 50% (I 50 ) in R plants was 651.13 µM, while that in the S population was 40.64 µM (Table 5 and Figure 3A). In other words, the R population was 1602 fold more resistant to imazamox than the S counterpart (Table 5). No differences were found in the basal enzyme activities ( Figure 3B). Table 5. Parameters of the log-logistic model (Equation (1)) of the inhibition of ALS enzymatic activity in response to imazamox concentration for R and S E. heterophylla populations.

Without CytP450 Inhibitor
The Proto IX accumulation was significantly different at 500 μM of fomesafen, in that R population accumulated significantly less Proto IX than the S population (Table 6 and Figure 4A). In R plants, the Proto IX accumulation seems to be constant from 1000 to 5000 μM.

Without CytP450 Inhibitor
The Proto IX accumulation was significantly different at 500 µM of fomesafen, in that R population accumulated significantly less Proto IX than the S population (Table 6 and Figure 4A). In R plants, the Proto IX accumulation seems to be constant from 1000 to 5000 µM. Table 6. Proto IX accumulation (expressed as nmol per gram of fresh weight ± standard error of the mean n = 9) in both R and S E. heterophylla populations treated with different fomesafen concentrations.

Without CytP450 Inhibitor
The Proto IX accumulation was significantly different at 500 μM of fomesafen, in that R population accumulated significantly less Proto IX than the S population (Table 6 and Figure 4A). In R plants, the Proto IX accumulation seems to be constant from 1000 to 5000 μM. Table 6. Proto IX accumulation (expressed as nmol per gram of fresh weight ± standard error of the mean n = 9) in both R and S E. heterophylla populations treated with different fomesafen concentrations.

With CytP450 Inhibitor
The R E. heterophylla population accumulated significantly less Proto IX than the S population with and without malathion (Table 7 and Figure 4B). However, while in the S population the accumulation was approximately the same with and without the CytP450 inhibitor, the application of malathion on R plants produced an increase in the Proto IX levels, approximately five times higher compared to the treatment of fomesafen without the inhibitor (Table 7 and Figure 4B). Table 7. Proto IX accumulation (expressed in nmol per gram of fresh weight ± standard error of the mean n = 9) in R and S E. heterophylla populations treated with a field dose of fomesafen without and with malathion (CytP450 inhibitor).

Physiological Assays
The weed survival after repeated applications of glyphosate (as the unique control tool) for more than five consecutive years triggers evolutionary adaptations that result in resistance to glyphosate [25,31,32]. In our case, the evolution of the R population (2011) was very fast, and with five years of repeated annual glyphosate applications of 720 g ae ha −1 (field dose in glyphosate-resistant soybean), the survival percentage reached 92% (Table 2). At the same time, the resistance level continued increasing exponentially, and 89% of surviving plants were recorded at the maximum evaluated rate (1800 g ae ha −1 ) after 3 years of recurrent selection with that dose ( Table 2). In this recurrent selection, high levels of resistance were achieved after nine years (one application per year), however, under field conditions, the time to select for resistance could be less, because the areas used for GR soybeans production in the summer cycle (safra) are generally occupied with GR maize in the winter cycle (safrinha) [6]. In addition, often more than one application of glyphosate is made per crop season, summer or winter, if GR varieties are grown [6].
The studies of the dose-response carried out in S and R E. heterophylla populations show high resistance to fomesafen, imazamox and glyphosate (Table 2 and Figure 1). In the case of dicotyledonous weeds, only four cases with multiple resistances (three sites of action) to these herbicides have been confirmed worldwide. The first was a population of Amaranthus palmeri detected in 2016 in Arkansas (USA) [33], the second was a population of A. tuberculatus detected in 2017 in Ontario (Canada) [34], and both species were found in glyphosate-resistant cotton and soybeans fields [33]. The third population was Epilobium ciliatum confirmed in 2018 in Lolol (Chile) in an olive orchard [35] and the fourth a Parthenium hystherophorus population confirmed in 2020 in Monte Cristi (Dominican Republic) in a field of banana orchards [26].

Biochemical Assays
The activity of the ALS enzyme showed that there were no differences in the basal activity, but differences were found between the populations in the amount of imazamox necessary to inhibit the ALS activity by 50%. These results could imply the presence of a TSR mechanism for imazamox (IMI) and cross-resistance to other ALS chemical families cannot be discarded, as has been demonstrated for another population collected in Nova Boa Vista Li Perau, Brazil [15]. This R population has also shown that resistance to imazamox is due to two resistance mechanisms, the main one was the Ser653Asn mutation and a secondary mechanism due to the exudation of the 14 C-imazamox by roots [15]. However, Mendes et al. [16] have studied resistance to imazamox in two populations collected from soybean fields in Mamborê and Kaloré (State of Paraná, Brazil) that showed a Trp574Leu mutation.
Significant differences in the EPSPS basal activity and response to glyphosate concentration were detected between S and R populations, concluding that resistance levels depend on the TSR mechanism [10,32]. Our results showed a resistance index of 9.98, which indicates that one or more mutations could be involved in resistance to glyphosate. In addition, the R populations presented a basal enzyme activity 3.37 times greater than that detected in the S population, which implies that other molecular mechanisms such as copy number and/or overexpression of EPSPS could also be involved in the resistance of E heterophylla to glyphosate (Table 5 and Figure 3B). Glyphosate resistance in dicot weeds is due to different mechanisms grouped and commonly known as non-target site resistance (NTSR) and target site resistance (TSR) mechanisms [32,36]. The TSR mechanisms are those related to the EPSPS, either by a loss of affinity between the linking enzyme and glyphosate caused by mutations and/or by the EPSPS overexpression as have been concluded for species such as Amaranthus spp, Bidens pilosa and B. subalternans, among others [37][38][39][40][41].
Differences in poor control and resistance factors between S and R populations confirm the high resistance of E. heterophylla to fomesafen (Table 3 and Figure 1). The accumulation of Protoporphyrin IX (Proto IX) due to the application of PPO inhibitor herbicides has been a widely used tool for the determination of resistance in weeds, where susceptible weeds accumulate much more of this metabolite [10,26,42]. Our results show that the accumulation of Proto IX in the S population was 7.18 times greater than in the R population treated with 2000 µM of fomesafen (Table 6 and Figure 4A), which is in agreement with the results obtained by other authors on dicot weeds such as E. ciliatum [42], P. hysterophorus [26] and recently on a population of E. heterophylla [11]. The use of malathion as an inhibitor of CytP450 would interfere with metabolisation of fomesafen in non-toxic products and therefore increased the accumulation of Proto IX in R plants, helping us to be able to interpret the possible mechanisms involved in resistance to PPO (Table 7 and Figure 4B). The application of fomesafen + malathion did not affect the accumulation of Proto IX in the S population; however, this treatment increased by 4.66 times the accumulation of Proto IX on the R population compared to the treatment done under the same conditions but only with fomesafen (Table 7 and Figure 4B). Nevertheless, the accumulation of Proto IX detected on R plants treated with fomesafen and malathion represented 22% of the Proto IX measured on S plants (Table 7). This result could be due to two mechanisms; due to the enhanced metabolism of fomesafen (NTSR), which is in the fixation stage since it is not yet able to fully confer resistance to fomesafen, or due to a low fomesafen sensitivity of PPO putatively associated with a TSR. In contrast, a trial conducted in a different population of E. heterophylla treated with lactofen (PPO inhibitor) and lactofen + malathion concluded that herbicide metabolism is not involved in resistance and only a mutation in PPO can explain the high level of resistance to lactofen [8]. In addition, molecular studies have shown that a mutation in the PPO gene at position 98 is the most likely resistance mechanism involved in Ambrosia artemisiifolia [43] and A. palmeri [44].

Conclusions
The first case of multiple resistances to ALS, EPSPS and PPO in Euphorbia heterophylla collected from a soybean field in Southern Brazil has been confirmed, which was obtained by recurrent selection from 2011 to 2019. Biochemical analyses suggest that resistance to imazamox and glyphosate is based on TSR, while fomesafen resistance would be based on TSR and NTSR mechanisms. The condition of multiple resistance to three MOAs (ALS, EPSPS and PPO) could occur soon in Brazilian soybean fields, as demonstrated by the recurrent selection, if alternatives to weed management, other than herbicides, are not included.