agronomy Evolving Multiple Resistance to EPSPS, GS, ALS, PSI, PPO, and Synthetic Auxin Herbicides in Dominican Republic Parthenium hysterophorus Populations. A Physiological and Biochemical Study

: Two Parthenium hysterophorus populations resistant (R) and susceptible (S) harvested in banana crop from the Dominican Republic were studied. All S plants died when the herbicides were applied at ﬁeld dose, except with paraquat. For the R population, the order of plant survival was as follows: glyphosate and paraquat > ﬂazasulfuron > glufosinate > fomesafen > 2,4-D. The resistance factors obtained in the dose–response assays showed a high resistance to glyphosate, ﬂazasulfuron, and fomesafen, medium resistance to glufosinate and 2,4-D, and a natural tolerance to paraquat (resistance factor (RF) = 1.0). The I 50 values obtained in the 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS), acetolactate synthase (ALS), and glutamine synthetase (GS) activity studies with glyphosate, ﬂazasulfuron, and glufosinate, respectively, were greater in R than in S. The e ﬀ ect of fomesafen was measured by the Proto IX levels, obtaining ﬁve times more Proto IX in the S than in the R population. The resistance to 2,4-D in the R was determined by the lower accumulation of ethylene compared to the S population. The studies with 14 C-paraquat conclude that the lower absorption and translocation in both the R and S populations would explain the natural tolerance of P. hysterophorus . This is the ﬁrst case of multiple resistance to herbicides with di ﬀ erent mechanisms of action conﬁrmed in P. hysterophorus.


Introduction
Parthenium hysterophorus is a weed that is widely distributed in Africa, Asia, and Oceania and is native to tropical America. For example, in Cuba, it is considered as one of the most noxious species [1,2]. The International Union for Conservation of Nature (IUCN) considers it to be one of the 100 most invasive species in the world [3].
This weed has a high seed production capacity (130,000-200,000 seeds m −2 ), as well as persisting in the soil and germinating in a wide range of temperatures at any time of the year. These characteristics have contributed to the propagation via flowing water, the movement of vehicles, animals, and machinery, or it can be blown by wind (presence of achenes), facilitating its dissemination

Plant Material
Seeds of P. hysterophorus from the provinces of Monte Cristi and Azua (Dominican Republic) were used in this study. The seeds of populations considered to be resistant (R) were collected during 2016 from banana orchards located in Monte Cristi, where the technicians and farmers of the area reported difficulty in controlling outbreaks with different herbicides and MOAs (EPSPS, GS, ALS, PSI, PPO, and synthetic auxin). The seeds of the population considered to be susceptible (S) were collected from organic banana orchards in a region of Azua that has never been exposed to glyphosate or other herbicide treatment. The seeds were germinated in plastic containers (10 × 10 cm x 6.3 cm) that were filled with moistened peat and covered with parafilm. The plastic containers were kept in a growth chamber (28/18 • C (day/night), 16 h photoperiod, 850 µmol m −2 s −1 light density and 60% relative humidity) until germination. Seedlings were transplanted into 250 cm 3 pots (1 plant/pot) filled with sand/peat (1:1 v/v). Throughout the experiments, pots were placed in a room with controlled conditions and were watered daily. Table 1 shows the commercial herbicides used in this study to determine the physiological effects on the populations of P. hysterophorus. To perform the absorption and translocation test in P. hysterophorus, we used 14 C-paraquat (specific activity 6.2789 MBq/mg, American Radiolabeled Chemicals, Inc., Saint Louis, MO, USA). The commercial formulation of paraquat is described in Table 1.

Fast Screening Assays
To ascertain whether these populations survive the application of herbicides by a mechanism of resistance and not because of incorrect application, greenhouse fast screening assays were carried out. The characteristics of the applied herbicides are shown in Table 1. All herbicides were applied in post-emergence plants with four to six true leaves. The experiments were arranged in a completely randomized design using ten plants per herbicide, applying the field doses (Table 1) in a treatment chamber (De Vries Manufacturing, Hollandale, MN, USA) equipped with a Tee Jet 8002 EVS flat fan nozzle calibrated to deliver 200 L ha −1 at 200 kPa at a height of 50 cm. Ten untreated plants were used as a control for all herbicides. This was done for both the resistant (R) and susceptible (S) P. hysterophorus populations. Three weeks after application (WAA), the number of surviving plants was counted. The experiments were conducted twice at different times (spring and fall) and the average of the two was obtained.

Dose-Response Assays
After the fast screening assays, we created a dose-response curve for each populationof P. hysterophorus (R and S) for the previously described herbicides (Table 1). Seven doses were used (0, X/16, X/8, X/4, X/2, X, 2X), where X is the field dose for each herbicide and zero (0) indicates untreated plants as control. The treatments were applied in the laboratory chamber mentioned above using the same application calibration. The experiment was arranged in a completely randomized design using five plants that were selected for each dose, for a total of 35 plants for each herbicide, for each population, for both the resistant population (R) and the susceptible population (S). After the treatments were completed, the plants were taken to a greenhouse and held under a temperature regime of 26/18 • C day/night, with water irrigation controlled at appropriate levels. Three WAA, the plants were weighed, and they were dried at 60 • C for four days in a climatic chamber.

EPSPS Enzyme Activity
The extraction of the EPSPS enzyme was carried out according to the methodology described by Bracamonte et al. [34]. Five grams of young leaf tissue from the R and S populations of P. hysterophorus were ground to a fine powder using liquid nitrogen in a cold mortar for EPSPS extraction and total protein content measurement [37]. The specific activity of EPSPS was tested in plants of both populations in the presence of different concentrations of glyphosate (0, 0.1, 1, 10, 50, 100, 200, 500, and 1000 µM) and in the absence (basal activity) of glyphosate. The EPSPS enzyme activity was expressed as percentage of enzymatic activity in the presence of glyphosate with respect to the control (without glyphosate). The EPSPS activity was determined by the amount of phosphate (µmol) released per µg of total soluble protein (TSP) per min. Three replicates for each concentration were tested in a completely randomized design and the experiment was conducted twice.

GS Enzyme Activity
The response of the glutamine synthetase enzyme (S) to glufosinate was determined using extracts of crude protein isolated from the R and S populations of P. hysterophorus following the methodology described by Tahmasebi et al. [38]. The following glufosinate concentrations were applied: 0, 15.6, 31.25, 62.5, 125, 250, 500, 1000, 2000, and 4000 µM. The concentration of glufosinate that reduced the S activity by 50% (I 50 ) was used to calculate the resistance factor (RF) values. The total protein was also measured following the Bradford method [37]. This experiment was conducted twice using three replicates per herbicide concentration in a completely randomized design.

ALS Enzyme Activity
The acetolactate synthase (ALS) activity was determined following the methodology described by Hatami et al. [39] and Tahmasebi et al. [38], using 3 g of young leaf tissue per each population. The crude extract of the ALS enzyme was used to measure the enzyme activity at increasing concentrations (0, 0.1, 1, 10, 100, 1000, and 10,000 µM) of flazasulfuron. The acetoin absorbance, obtained from the decarboxylation of acetolactate, was measured with a spectrophotometer (Beckman DU-640, Fullerton, CA, USA) at 520 nm. The concentration of flazasulfuron that reduced the ALS activity by 50% (I 50 ) was used to calculate the RF values. The total content of protein in the crude extract was also measured by the Bradford colorimetric method [37]. The experiment was carried out twice with three replicates per herbicide concentration and per population, following a completely randomized design.

14 C-Paraquat Absorption and Translocation
14 C-paraquat absorption and translocation were evaluated at 6, 12, 24, 48, and 96 h after treatment (HAT). Once the R and S plants had four leaves, they were treated with a solution of 14 C-paraquat plus a commercial paraquat formulation. The solution applied contained 14 C-paraquat, providing 0.834 kBq µL −1 at a final concentration of 100 g ia ha −1 of paraquat in 300 L. One microliter per plant of solution was applied on the adaxial surface of the second youngest fully expanded leaf. After treatment, the plants were maintained in the growth chamber for 12 h in the dark before being exposed to light [38,40].
To determine the absorption, the 14 C-paraquat-treated plants were harvested, and the treated leaves were washed three times with 1 mL of distilled water to recover the non-absorbed 14 C-paraquat. The washing solution was mixed with 2 mL of scintillation liquid (Ultima Gold, Perkin-Elmer, BV BioScience Packard Groningen, Netherlands, MA). The samples were reserved for radioactivity analysis. To determine translocation, the whole plants were carefully removed from the pot and washed, mainly the roots. The plants were individually divided into treated leaf, the root system, and the remainder of the plant. The samples were stored in flexible combustion cones (Perkin-Elmer, BV BioScience Packard) and dried in an oven at 60 • C for 48 h. Subsequently, the samples were combusted in a Packard Tri Carb 307 biological oxidizer (Packard Instrument Co., Downers Grove, IL, USA). The CO 2 released from the combustion was captured in 18 mL of a mix of Carbo-Sorb E and Permafluor (1:1 v/v) (Perkin-Elmer, BV BioScience Packard, Groningen, Netherlands). The radioactivity in the washing solution and combustion samples was quantified by liquid scintillation spectrometry in an LS 6500 scintillation counter (Beckman Coulter Inc., Fullerton, CA, USA) for 5 min per sample, and the measurement was repeated again 24 h later. The radioactive values were used to calculate the recovery percentage as % 14 C-paraquat recovered = (kBq in treated leaf + kBq remaining in plant + kBq in root system + kBq from washes/total kBq applied) x 100. The experiments were arranged in a completely random design with four replicates per population at each evaluated time point, and the experiment was conducted twice.

Proto IX Levels
Determination of the Proto IX levels was carried out following the methodology described by Dayan et al. [41]. Leaf discs (approximately 0.1 g each) from the two populations were incubated individuality in a 6 cm diameter Petri dish containing 6 mL of 2% (wt/v) sucrose, 1 mM 2-(N-morpholin) ethanesulfonic acid, and 1000 µM fomesafen (technical grade) for 20 h at 25 • C in the dark. The same experiment was conducted but without herbicide for the control. All Proto IX extractions from the leaf discs were carried out under a dim, green light source. The extracts were stored in light tight glass vials at −20 • C until analysis. The analysis was carried out by LC using a 126 Gold HPLC system from Beckman Coulter (Fullerton, California, USA) equipped with a Jasco FP−1520 fluorescence detector (Easton, Pennsylvania, USA, excitation and emission wavelengths at 400 and 630 nm, respectively) and a SphereClone ® ODS column (25 cm x 4.6 mm, 5 micron particle size) from Phenomenex (Torrance, California, CA, USA).
Proto IX levels in the extracts were quantified using a calibration curve obtained with a commercially available Proto IX standard from Sigma Aldrich (St. Louis, Missouri, USA). The data are expressed as nmol per g of fresh weight. All treatments were carried out three times using three replicates with and without herbicide in a completely randomized design.

2,4-D and Ethylene Production
R and S plants with 6 leaves were treated with 2,4-D at 0, 125, 250, 500, 1000, and 2000 g ai ha −1 following the methodology described above for the dose-response assays. After 24 h of treatment, 400 g of leaf tissue was placed in a 10 mL syringe with 1 mL of distilled water and then sealed [42]. The syringes were placed in a dark incubator at 27 • C for 4 h. The ethylene (C 2 H 4 ) in 1 mL of the gas at the top of the syringe was analyzed by gas chromatography [38,43]. The C 2 H 4 content was expressed as nanoliters per gram of fresh weight per hour. There were five replicates per treatment, and the experiment was conducted twice.

Statistical Analysis
The herbicide concentration causing 50% growth reduction (GR 50 ) of the plants and the herbicide concentration causing 50% inhibition of enzyme activity (I 50 ) were calculated by non-linear regression analysis using the R package drc (R Core Team) [44]. A three-parameter log-logistic model was selected for modeling growth reduction or enzyme inhibition (Equation (1)): where Y is the dry weight of the harvested plants or enzyme activity, d is the coefficient corresponding to the upper asymptote, c the lower limit (fixed at 0), the coefficient b is the slope at the inflection point, e the herbicide concentration required to inhibit shoot growth or enzyme activity by 50% (i.e., GR 50 or I 50 , respectively), and X the herbicide dose. The resistance factors (RF = R/S) were computed as R-to-S GR 50 or I 50 ratios. The effects of the population, time (HAT), and 14 C-paraquat absorption and translocation were tested using ANOVA. The population was treated as a fixed factor, and time was considered as a random factor. The means and standard errors of 14 C-paraquat absorption and translocation were computed for all plant parts, and the means were tested for group differences and then compared. For each analysis, the assumptions of equality of variance and normal error distribution were evaluated and they met after angular transformation of the response variable (arcsine of the square root of fractions). The Tukey HSD test at 5% probability was used to separate the means. The ANOVAs were performed using the Statistix software (version 10.0) from Analytical Software.

Fast Screening Assays
The fast screening test at field doses showed that the S populations were 100% controlled by glyphosate, glufosinate, fomesafen, flazasulfuron, and 2,4-D herbicides ( Table 2). The herbicides 2,4-D and fomesafen showed the lowest values of plant survival with respect to the control when they were applied at post-emergence, while the other herbicides (glufosinate and flazasulfuron) had higher percentages of plants that survived at field doses. Glyphosate had zero effectiveness on the R population. Paraquat had no effect on the R and S P. hysterophorus populations, which leads us to suspect that this species has a natural tolerance to paraquat.

Dose-Response Assays
Dose-response assays showed the existence of the first case of multiple resistance to herbicides with different mechanisms of action in P. hysterophorus collected from banana orchards in the Dominican Republic (Table 3). The GR 50 values obtained for the different herbicides (except paraquat) showed a high efficiency in reducing the 50% dry weight of the S population of P. hysterophorus, which is in good agreement with the results obtained in the fast screening test (Tables 2 and 3). Based on the GR 50 values of the R population, we calculated resistance factors (RF: GR 50 of R/GR 50 of S) from 4.9 to 27.1, which indicated the following order of resistance to the tested herbicides: glyphosate > flazasulfuron > fomesafen > glufosinate > 2,4-D (Table 3).
; where Y is the dry weight of the harvested plants or enzyme activity, d is the coefficient corresponding to the upper asymptote, c the lower limit (fixed at 0), the coefficient b is the slope at the inflection point, e the herbicide concentration required to inhibit shoot growth or enzyme activity by 50% (GR 50 ), and X the herbicide dose. Resistance factors (RF = R/S) are the ratio of R-to-S GR 50 . The p-value shows the level of significance of the non-linear model. b g ae ha −1.
Paraquat is a contact herbicide that shows very rapid efficacy in the appearance of chlorotic leaves after 24 HAA and necrosis one week later. However, in tolerant species, the plant begins to regrow in meristematic areas. In our study, the GR 50 results for the R and S populations showed very high and similar values (Table 3). These results, together with those obtained in the fast screening test, allow us to conclude that there is a natural tolerance to paraquat in this species (Tables 2 and 3).

EPSPS Enzyme Activity
The EPSPS was inhibited in situ by increasing glyphosate concentrations. The required glyphosate concentration to inhibit EPSPS by 50% (I 50 ) was approximately 24 times greater in the R population than in the S population (Table 4). However, no differences were found between R and S in the enzyme basal activity (128.46 and 119.11 nmol Pi/ µg TSP/min, respectively).
; where Y is the dry weight of the harvested plants or enzyme activity, d is the coefficient corresponding to the upper asymptote, c the lower limit (fixed at 0), the coefficient b is the slope at the inflection point, e the herbicide concentration required to inhibit shoot growth or enzyme activity by 50% (I 50 ), and X the herbicide dose. Resistance factors (RF = R/S) are the ratio of R-to-S I 50 . The p-value shows the level of significance of the non-linear model.

GS Enzyme Activity
The glufosinate doses required to reduce GS activity by 50% (I 50 ) were 946.6 and 32.2 µM for the R and S populations, respectively. The RF value for R was 29.3 times greater that of S (Table 4). No differences were found in the basal enzyme activity (226.29 and 240.01 nmol of glutamine/mg TSP/h, respectively).

ALS Enzyme Activity
The quantity of flazasulfuron needed to inhibit the ALS activity by 50% (I 50 ) in the R population was 624.4 µM, while that in the S population was 46.3 µM. In other words, the R population was 13.5 fold more resistant to flazasulfuron than the S population (Table 4). No differences were found in the basal enzyme activity (332.34 and 317.20 nmol of acetoine/mg TSP/h, respectively).

14 C-Paraquat Absorption and Translocation
The average total recovery of the applied 14 C-paraquat was >90% in both the R and S P. hysterophorus populations. The 14 C-paraquat absorption was not significantly different between the R and S plants throughout the time periods and increased over time up to 70% at 96 HAT (Table 5). In addition, the translocation of 14 C-paraquat from the treated leaf to the rest of the shoot and roots was minimal and similar in both populations. There was a minimal translocation of 1.2% to 4.6% to the rest of the shoot plant, and a negligible amount of 14 C-paraquat translocation to the root system during experimental time (Table 5).

Proto IX Levels
After the application of 1000 µM fomesafen, the P. hysterophorus R population accumulated significantly less Proto IX than the S population. The accumulation values were 0.50 ± 0.03 and 2.75 ± 0.14 nmol g −1 fresh weight in the R and S populations, respectively.

2,4-D and Ethylene Production
The ethylene accumulation in the R and S P. hysterophorus plants increased along with the dose of herbicide applied (Figure 1). The R and S control plants produced similar amounts of ethylene (data not shown). The ethylene accumulation at 1000 g ai ha −1 of 2,4-D was 5.0 times smaller in the R population than in the S population.

Discussion
Our results demonstrated different levels of resistance in P. hysterophorus to the tested herbicides. The highest resistance factor was for glyphosate, followed by flazasulfuron, fomesafen, and glufosinate. The lower values of GR50 for the susceptible population showed that all of the herbicides except paraquat were useful to control the susceptible population of this weed in banana crops [34]. The natural tolerance to paraquat in this species was confirmed in this study; it was reported for the first time in 1991 in coffee crops [45].
Compared with previous glyphosate research, the level of P. hysterophorus resistance was four times higher than reported in Colombia and 1.5 times lower than in the United States [35,46]. According to a study conducted in Cuba and the Dominican Republic, the values are similar if the greatest resistance factor reported in that study is taken into account [34]. The difference in resistance values can be attributed to the doses used in each country and crop, environmental factors, soil, tillage practices, and management practices [47][48][49].
Regarding 2,4-D, the resistance factor was similar to that reported in citrus orchards (also in the Dominican Republic) [36]. Finally, for ALS inhibitors, Gazziero et al. [50] reported RF values ranging between 17 and 21 for the different active ingredients evaluated (chloransulam-methyl, chlorimuronethyl, imazethapyr, iodosulfuron + foramsulfuron), and although flazasulfuron is not included, the values are similar.
When there is no significant difference in the EPSPS activity and basal activity, it demonstrates the absence of Target-Site Resistance (TSR) as a resistance mechanism in P. hysterophorus, according to the conclusions made by Palma-Bautista et al. [19]. However, in the study conducted in the Caribbean Islands, accessions showed different resistance mechanisms. Some of them showed significant differences in the EPSPS activity and basal activity between susceptible and resistant populations, concluding that resistance levels depend on the TSR mechanism. Furthermore, in other accessions, the results indicated that Non-Target-Site Resistance (NTSR) mechanisms that accumulated in this adhesion were also involved [34].
A similar GS (glutamine synthetase) basal activity was found in this study in both the R and S populations. However, the enzyme activity inhibition (I50) was 29-fold greater in the susceptible

Discussion
Our results demonstrated different levels of resistance in P. hysterophorus to the tested herbicides. The highest resistance factor was for glyphosate, followed by flazasulfuron, fomesafen, and glufosinate. The lower values of GR 50 for the susceptible population showed that all of the herbicides except paraquat were useful to control the susceptible population of this weed in banana crops [34]. The natural tolerance to paraquat in this species was confirmed in this study; it was reported for the first time in 1991 in coffee crops [45].
Compared with previous glyphosate research, the level of P. hysterophorus resistance was four times higher than reported in Colombia and 1.5 times lower than in the United States [35,46]. According to a study conducted in Cuba and the Dominican Republic, the values are similar if the greatest resistance factor reported in that study is taken into account [34]. The difference in resistance values can be attributed to the doses used in each country and crop, environmental factors, soil, tillage practices, and management practices [47][48][49].
Regarding 2,4-D, the resistance factor was similar to that reported in citrus orchards (also in the Dominican Republic) [36]. Finally, for ALS inhibitors, Gazziero et al. [50] reported RF values ranging between 17 and 21 for the different active ingredients evaluated (chloransulam-methyl, chlorimuron-ethyl, imazethapyr, iodosulfuron + foramsulfuron), and although flazasulfuron is not included, the values are similar.
When there is no significant difference in the EPSPS activity and basal activity, it demonstrates the absence of Target-Site Resistance (TSR) as a resistance mechanism in P. hysterophorus, according to the conclusions made by Palma-Bautista et al. [19]. However, in the study conducted in the Caribbean Islands, accessions showed different resistance mechanisms. Some of them showed significant differences in the EPSPS activity and basal activity between susceptible and resistant populations, concluding that resistance levels depend on the TSR mechanism. Furthermore, in other accessions, the results indicated that Non-Target-Site Resistance (NTSR) mechanisms that accumulated in this adhesion were also involved [34].
A similar GS (glutamine synthetase) basal activity was found in this study in both the R and S populations. However, the enzyme activity inhibition (I 50 ) was 29-fold greater in the susceptible population than in the resistance population. This result indicates that glufosinate resistance in P. hysterophorus is likely to be target-site based, although this hypothesis needs to be investigated. Similar results were found for Epilobium ciliatum in Chilean olive orchards [38]. Avila-Garcia and Mallory-Smith [59] and Jalaludin et al. [60] noted a similar sensitivity of the GS enzyme between susceptible and resistant populations, indicating that glufosinate resistance is not conferred by the target site.
Similar results were determined for the ALS enzyme activity, where there were no differences in the basal activity but differences between the populations were found in the quantity and needed to inhibit the ALS activity by 50%. These results demonstrate the presence of a TSR mechanism for flazasulfuron. Similar studies were conducted in Rapistrum rugosum for different ALS herbicides, in Euphorbia heterophylla and soft wheat for imazamox, and in Conyza canadensis and Epilobium ciliatum for flazasulfuron [38,39,[61][62][63].
The poor control, resistance factors, and lack of differences in the absorption and translocation of paraquat between susceptible and resistant populations confirm the natural tolerance of P. hysterophorus to paraquat, explained by the low absorption and translocation of the herbicide [38,40,64]. This natural tolerance was reported in Kenya for the first time [45]. By contrast, in perennial ryegrass (Lolium perenne L.), the tolerance was explained due to the higher activities of catalase and superoxide dismutase, along with elevated peroxidase activity, in paraquat-tolerant lines [65,66]. In young leaves of squash cultivars, the tolerance to paraquat is due to oxidative and abiotic stresses, such as the change in catalase, peroxidase, ascorbate peroxidase, and glutathione reductase activities [67]. In Chenopodium rubrum, the tolerance is explained by the activity of superoxide dismutase and peroxidase [68].
A useful methodology to determine the effectiveness of PPO herbicides is to evaluate the accumulation of protoporphyrin IX, where susceptible plants will accumulate more of this enzyme [41]. The lower accumulation of Proto IX in the R population confirmed the resistant to fomesafen. Similar results were found by Tahmasebi et al. [38], which supported the evolution of resistance to pyraflufen-ethyl by recurrent selection in E. ciliatum.
In the case of 2,4-D, the accumulation of ethylene was evidenced in the R and S plants treated with 2,4-D and has been reported by several authors [38,42,[69][70][71]. The lower accumulation of ethylene in R plants may be due to the herbicide not reaching its nuclear protein receptor complex to suppress the genes responsible for ethylene production, which respond to auxins [71][72][73]. Susceptible plants exposed to 2,4-D die due to unregulated auxin activity and the accumulation of reactive oxygen species (ROS), ABA, and ethylene [72,74]. In previous P. hysterophorus studies, resistance to 2,4-D was found to be due to two simultaneous mechanisms of NTSR-less translocation of the herbicide and enhanced metabolism-in which cytochorme-P450 may be involved [36]. The mechanism of resistance to synthetic auxins in Papaver rhoeas is due to the reduced translocation of 2,4-D, which leads to a lower production and accumulation of ethylene and the survival of resistant plants [73]. However, in a more recent study in the same species, an additional mechanism was found to be due to enhanced metabolism mediated by cytochrome P450 [42].

Conclusions
The first cases of multiple resistance to EPSPS, GS, ALS, PPO, synthetic auxin, and natural tolerance to PSI in Parthenium hysterophorus harvested from banana orchards in the Dominican Republic has been confirmed. Biochemical analyses suggest that resistance to glyphosate, glufosinate, flazasulfuron, and fomesafen is based on TSR and 2,4-D on NTSR, while tolerance to paraquat is also based on NTSR. Due to the complex evolution of multiple resistance to five MOAs and naturally tolerant to paraquat in this species, farmers should implement integrated weed management with non-chemical alternatives for the management of these weed populations.