Plant Materials

In September 2017, ripe seeds (F₀ generation) of a suspected resistant (R, HN-02) A. retroflexus population were collected from at least 100 plants in a soybean field in Nenjiang County (49.11°N, 124.94°E), Heilongjiang Province, where fomesafen and nicosulfuron have failed to control A. retroflexus since 2014. The sampled field was approximately 0.5 ha and had been under a repeated soybean–corn rotation since 1980. A herbicide-susceptible population (S, ST-1) was collected from Mount Tai (36.05°N, 117.03°E), Shandong Province, and had never been exposed to herbicides. After being cleaned and dried, the seeds from individual plants were collected and pooled into a single composite sample for each field and then stored at 4 C in paper bags for later use.

Single-dose Herbicide Resistance Testing

Before planting, both populations of seeds (F₀ generation, approximately 110) were germinated in petri dishes that contained two layers of filter paper (Whatman No. 1) moistened with 6 ml water and incubated in a growth chamber (33/23 C, 12/12 h day/night). The germinated seedlings were then transplanted into plastic pots (15-cm diameter, 12-cm height) containing loam; the final quantity of plants that were treated with herbicide was 10 per pot. The plants were grown to the 3- to 4-leaf stage in a controlled greenhouse environment (30/20 C, day/night) and watered every other day to maintain soil moisture. For each population, 50 plants were treated with fomesafen at 225 g ai ha⁻¹ (the field-recommended rate), another 50 plants were treated with nicosulfuron at 48 g ai ha⁻¹ (the field-recommended rate), and 10 plants were maintained as untreated controls. All herbicides were sprayed using a moving-nozzle cabinet sprayer equipped with one TeeJet^® 9503EVS (Greenman Machinery Company, Beijing, China) flat-fan nozzle (450 L ha⁻¹ at 280 kPa). The percentage of surviving plants was visually measured at 21 d after treatment. The plants were considered dead if their leaf tissue was brittle and necrotic and if no new growth point had developed. Afterward, approximately 100 mg of young leaf tissue of each individual surviving plant was removed and used for total RNA extraction via TRIzol Reagent (Thermo Fisher Scientific, Waltham, MA, USA).

Gene Sequencing and Acquisition of Homozygous TSR Resistant Seeds

The isolated RNA was immersed in 70% ethanol and stored at −80 C. RNA was treated with TransScript One-Step gDNA Removal and cDNA Synthesis SuperMix (TransGen Biotech, Beijing, China) to remove any gDNA contamination. A pair of primers (forward primer [ALS-F]: 5′-TGCGATGTTCTCGTTGAAGC-3′; reverse primer [ALS-R]: 5′-CCCTTCTTCCATCACCCT-3′; 1,711 bp) previously reported by Liu et al. (Reference Liu, Yuan, Du, Guo, Li, Bi and Wang2015) was used to amplify the ALS sequence. To amplify the PPX2 sequence, primers (PPX2-F: 5′-GCCATCGCCATTGTCAGTTTCA-3′; PPX2-R: 5′-ATGGATGAGAAGACCGCGTAA-3′; 1,450 bp) were designed (Primer Premier 5.0, Biosoft, Palo Alto, CA, USA) according to GenBank accession no. KY882137.1 and a study by Patzoldt et al. (Patzoldt et al. Reference Patzoldt, Hager, McCormick and Tranel2006). Both primer pairs were designed to amplify all reported mutation sites with PPO and ALS resistance, respectively. The PPX2 sequence for A. retroflexus has been published and is available in the NCBI database (MK716317-MK716318). Polymerase chain reaction (PCR) was performed in a final reaction volume of 25 µl, which consisted of 1 μl of cDNA, 1 μl of each primer (synthesized by Sangon Biotech, Shanghai, China), 12.5 μl of TransTaq HiFi PCR SuperMix II (TransGen Biotech, Beijing, China), and 9.5 μl of nuclease-free water. PCR was performed on a T100 thermal cycler (Bio-Rad, USA) with the following amplification conditions of the target gene: 94 C; 35 cycles of 30 s at 94 C, 30 s at X C (56 C and 58 C for the PPX2 and ALS genes, respectively), and 1 min 30 s at 72 C; and then 10 min at 72 C. The PCR products were separated by a 1.0% agarose gel run in 1X TAE buffer at 120 V and then purified by an EasyPure Quick Gel Extraction Kit (TransGen Biotech, Beijing, China). The purified products were subsequently sequenced in both the forward and reverse directions by a commercial sequencing company (Sangon Biotech, Shanghai, China).

DNAMAN v. 5.2.2 software (Lynnon Biosoft, Quebec, Canada) was used to align the nucleotide sequence positions of the PPX2 and ALS genes with those of the PPX2 gene of Palmer amaranth (Amaranthus palmeri S. Watson) (GenBank accession no. KY882137.1) and the ALS gene of A. retroflexus (GenBank accession no. AF363369.1), respectively. On the basis of results of the sequencing analysis, plants homozygous for both the ALS and PPX2 gene mutations were identified and grown in isolation in a growth chamber (33/23 C, 12/12 h day/night) to ensure that no other pollen was introduced. Subsequently, whole-plant herbicide-sensitivity bioassays and in vitro ALS activity assays were performed only for seeds collected from this separated R population, which was designated HN-02F1-1. For the HN-02F1-1 resistance test, 20 seeds were randomly selected for single-dose herbicide-resistance tests (following methods described earlier), and five plants were randomly selected for gene sequencing.

Whole-Plant Herbicide Sensitivity Bioassays

The plant germination and growth conditions of the R (HN-02F1-1) and S (ST-1) populations followed the methods described earlier. When the seedlings reached the 3- to 4-leaf stage, herbicides (four PPO inhibitors, six ALS inhibitors, and two synthetic auxins) were sprayed at a series of rates (Table 1) in accordance with the methods described earlier. The aboveground A. retroflexus materials were harvested at 21 d after herbicide treatment and then oven-dried for 72 h at 75 C, after which their dry weights were recorded. Each treatment was repeated three times, and the entire experiment was carried out twice.

Table 1.

Herbicides and herbicide application rates used in the whole-plant experiments.^a

^a Abbreviations: R, resistant; S, susceptible; ALS, acetolactate synthase; PPO, protoporphyrinogen oxidase; AS, aqueous solution; ME, microemulsion; WP, water power; OF, oil miscible suspension; WG, water dispersible granule; SC, suspension concentrate; EC, emulsifiable concentrate.

^b Fomesafen (250 g L⁻¹ AS, BrightMart CropScience, Foshan, China), fluoroglycofen-ethyl (10% ME, Huifeng Biological Agriculture, Jiangsu, China), acifluorfen (21.4% AS, Hanshen Biotechnology, Qindao, China), flumioxazin (50% WP, Sumitomo Chemical Corp, Japan), nicosulfuron (40 g L⁻¹ OF, Zhongshan Chemical Group, Huzhou, China), thifensulfuron-methyl (75% WG, FMC Corporation, Shanghai, China), cloransulam-methyl (84% WG, Dow AgroSciences, Beijing, China), flumetsulam (80% WG, Dow AgroSciences, Beijing, China), imazethapyr (5% AS, Lvbang Crop Science, Jinan, China), imazamox (4% AS, Jiangsu Agricultural Hormone Engineering Technology Research Center, Nanjing, China), fluroxypyr-meptyl (50% SC, Qingyuan Crown Resistant Weed Control, Qingdao, China), 2,4-D butylate (57% EC, Songliao Chemical, Liaoning, China).

^c The bold values represent the field-recommended rates.

In Vitro ALS Activity Assays

Seedlings of R (HN-02F1-1) and S (ST-1) populations were germinated and grown under the same conditions as those described earlier. At the 3- to 4-leaf stage, young leaf tissue was harvested and stored at –80 C. In accordance with the methods of Yu et al. (Reference Yu, Han, Vila-Aiub and Powles2010) and Han et al. (Reference Han, Yu, Purba, Li, Walsh, Friesen and Powles2012), ALS was extracted, and an in vitro assay of inhibitory activity was performed. The ALS activity was determined via colorimetry (530 nm) with a UV spectrophotometer (Thermo Fisher Scientific) by measuring the production of acetoin. The herbicide concentrations used for the activity assays are listed in Table 2. Three replications per assay were carried out, and the entire experiment was conducted twice, using independent ALS extraction.

Table 2.

Herbicides and herbicide application rates used in in vitro ALS activity assays.^a

^a Abbreviations: ALS, acetolactate synthase; R, resistant; S, susceptible.

Statistical Analyses

The data from both the whole-plant bioassays and the in vitro ALS activity assays were analyzed via ANOVA (SPSS v. 19.0, IBM, Armonk, NY, USA). The results indicated that there was no statistically significant (P > 0.05) interaction between repeats of either experiment, and thus, the combined data were analyzed using a four-parameter log-logistic curve via SigmaPlot (v. 12.5; SigmaPlot Software) as follows:

([1]) $$y = c + (d - c)/\{ 1 + \exp [b(\log x - {\rm{ }}\log G{R_{50}})]\} $$

where c represents the lower limit, d represents the upper limit, x represents the herbicide application rate, b represents the slope of the curve for the dose causing a 50% dry weight reduction in the aboveground growth (at GR₅₀), and y represents the growth response (percentage of the untreated controls). The resistance index (RI) was calculated as the ratio of the GR₅₀ value or the I₅₀ value (herbicide concentration required for 50% ALS activity inhibition) of the resistant and susceptible populations.

Single-Dose Testing, Gene Sequencing, and Resistant Population Purification

With respect to the single-dose testing, as expected, all plants of the HN-02 (R) population survived fomesafen at 225 g ai ha⁻¹, and 46 plants survived the nicosulfuron treatment; however, all plants of the ST-1 (S) population died.

Generally, all known target point mutations (ΔGly210, Arg-128-Leu, Arg-128-Gly/Met, Gly-399-Ala) causing PPO-inhibitor resistance in weed species are reported to occur in PPX2 instead of in PPX1 (Giacomini et al. Reference Giacomini, Umphres, Nie, Mueller, Steckel, Young, Scott and Tranel2017; Patzoldt et al. Reference Patzoldt, Hager, McCormick and Tranel2006; Rangani et al. Reference Rangani, Salas-Perez, Aponte, Knapp, Craig, Meitzner and Burgos2019; Rousonelos et al. Reference Rousonelos, Lee, Moreira, VanGessel and Tranel2012). Thus, a 1,437-bp fragment that encompasses all three known mutation sites within the A. retroflexus PPX2 gene was amplified from 40 individual plants (randomly selected from 96 surviving plants) of HN-02 and from the 10 control plants of ST-1. The results indicated that the similarity of the PPX2 gene sequences between HN-02 and ST-1 was approximately 99.6%. The change of a single nucleotide from AGG to GGG resulting in an Arg-128-Gly substitution was revealed when the PPX2 gene fragments were compared between the R and S samples. Furthermore, the corresponding sequence chromatograms revealed a single peak at all amino acid 128 mutation sites (Figure 1A), which indicated that only one allele is expressed, and the analyses of cDNA indicate that it is homozygous.

Figure 1.

cDNA sequencing results indicating (A) the Arg-128-Gly mutation in the PPX2 gene, (B) the Ala-205-Val (frequency of 65%) and Trp-574-Leu mutations (frequency of 20%) in the ALS gene in resistant Amaranthus retroflexus plants (HN-02, bottom), and (C) no Trp-574-Leu mutation in the ALS gene in the subpopulation HN-02F1-1 (bottom), compared with the susceptible plants (ST-1, top).

Furthermore, ALS gene fragments that include all eight known mutation sites (Ala-122, Pro-197, Ala-205, Asp-376, Arg-377, Trp-574, Ser-653, and Gly-654) (Tranel et al. Reference Tranel, Wright and Heap2019) were also amplified from the abovementioned 40 A. retroflexus plants. The gene-sequencing results of the 40 individual plants from HN-02 revealed a Trp-574-Leu substitution produced by the TGG-574-TTG mutation in 8 plants and an Ala-205-Val substitution produced by the GCT-205-GTT mutation in 26 plants, while 6 plants had no ALS-resistance mutations. However, no plants had the two abovementioned mutations coexisting in the ALS gene. The corresponding sequence chromatograms revealed a single peak in all 205 or 574 mutation sites (Figure 1B), indicating that all 205 and 574 mutants were homozygous. This phenomenon may have occurred because A. retroflexus is diploid and is mainly self-fertilized (Mandák et al. Reference Mandák, Zákravský, Dostál and Plačková2011; Yu and Powles Reference Yu and Powles2014b). Six plants had no ALS-resistance mutations but were resistant to nicosulfuron, indicating that NTSR is also present within the HN-02 population. All sequences of the 10 ST-1 individuals were identical, and no known mutations providing resistance to ALS or PPO inhibitors were found.

On the basis of the TSR gene-sequencing analysis, all 26 plants homozygous for both the PPX2 gene Arg-128-Gly mutation and the ALS gene Ala-205-Val mutation were crossed, producing an HN-02F1-1 subpopulation. To verify the purity of the HN-02F1-1 subpopulation, 20 seeds were randomly selected for cultivation, and five plants for sequencing. The results showed that all plants survived the field-recommended rates of fomesafen and nicosulfuron. The five individuals selected for sequencing were found to have identical sequences and to be homozygous for the Ala-205-Val mutation in the ALS gene and Arg-128-Gly mutation in the PPX2 gene (data not shown). Neither the Trp-574-Leu mutation (Figure 1C) nor other herbicide-resistance mutations were detected in the ALS gene or PPX2 gene.

In the past 20 yr, several PPO inhibitors, including oxidiazoles and diphenyl ethers, and ALS inhibitors, including SU and IMI, have been continuously used to control A. retroflexus in northeast China. Owing to unique planting patterns (soybean–corn rotations) and the complex history of herbicide treatment (the long-term use of PPO and ALS inhibitors) in this region, the specific crop-planting system and herbicide application history have played an important role in driving the evolution of multiple herbicide resistance in the HN-02 populations collected in the area. Therefore, investigation of the resistance mechanism against fomesafen and nicosulfuron in this specific A. retroflexus population will facilitate the understanding of the evolution of the multiple herbicide resistance in this weed species.

Molecular analysis of the surviving plants confirmed that three amino acid substitutions, Arg-128-Gly in the PPX gene and both Ala-205-Val and Trp-574-Leu in the ALS gene, were present in the plants from the HN-02 population. The most common amino acid substitution sites in ALS genes are Ala-205-Val and Trp-574-Leu, which perhaps reflect herbicide use patterns (SU + IMI select mostly for Trp-574-Leu, whereas IMI herbicides select mostly for Ala-205), selection pressure, and fitness costs correlated with ALS genetic mutations (Yu and Powles Reference Yu and Powles2014b). Imazethapyr is one of the most widely used herbicides in soybean fields (Li et al. Reference Li, Fan and Luan2008), while nicosulfuron is the main herbicide applied in cornfields in northeast China. These facts may explain why there are mainly Ala-205-Val mutations and relatively few Trp-574-Leu mutations in the HN-02 population. The Val mutation, in which Val is the most common substituted amino acid at the Ala-205 position in the ALS gene, has been documented in many weed species, including A. retroflexus (Ashigh and Tardif Reference Ashigh and Tardif2009; Matzrafi et al. Reference Matzrafi, Lazar, Sibony and Rubin2015; McNaughton et al. Reference McNaughton, Letarte, Lee and Tardif2005). The evolution of PPO-inhibitor resistance has been relatively rarely documented; to date, only 13 weed species have evolved resistance to PPO inhibitors (Heap Reference Heap2019). Moreover, the Arg-128 substitution has rarely been found in common ragweed (Ambrosia artemisiifolia L.) and A. palmeri (Dayan et al. Reference Dayan, Barker and Tranel2018). On the basis of our phenotype and sequencing results, the detection of the Arg-128-Gly mutation in the PPX2 gene and the Ala-205-Val mutation in the ALS gene may be related to resistance to fomesafen and nicosulfuron, respectively, in HN-02F1-1 plants. Many studies have shown that the frequency of herbicide-resistance mutations is correlated with specific weed species, geographic environment, herbicide selection pressure, and suitability (Délye et al. Reference Délye, Michel, Bérard, Chauvel, Brunel, Guillemin, Dessaint and Le Corre2010). Therefore, additional studies on the difference in mutation frequencies between amino acids 205 and 574 in the HN-02 population will be important to understand the multiple resistance of A. retroflexus.

Whole-Plant Herbicide Bioassays

The resistance profile of A. retroflexus to 12 herbicides was characterized in this study. The results indicated that the S population (ST-1) was sensitive (RI < 2.0) to all tested herbicides and that the R population (HN-02F1-1) displayed different levels of resistance to the herbicides tested. On the basis of RI values (Table 3), compared with the ST-1 population, the HN-02F1-1 population was 53-fold resistant to fomesafen, 15-fold resistant to fluoroglycofen-ethyl, 17-fold resistant to acifluorfen, and 1.9-fold resistant to flumioxazin. Therefore, the HN-02F1-1 population has evolved a high level of resistance to fomesafen, fluoroglycofen-ethyl, and acifluorfen while still being sensitive to flumioxazin. Whole-plant bioassays also revealed that HN-02F1-1 was highly resistant (RI > 10) to cloransulam-methyl (TP), nicosulfuron (SU), imazethapyr (IMI), and imazamox (IMI) and moderately resistant (RI, 5.0~10) to flumetsulam (TP) and thifensulfuron-methyl (SU). On the basis of the GR₅₀ values and the RI values (Table 3), compared with ST-1, HN-02F1-1 was 8.3- to 113-fold more resistant to all tested ALS inhibitors. Moreover, compared with those of the S group, the RI values of the R group for fluroxypyr-meptyl and 2,4-D butylate were lower than 2, indicating that the HN-02F1-1 population was still sensitive to those herbicides.

Table 3.

Parameter values of the four-parameter log-logistic equationfor calculating the GR₅₀values of Amaranthus retroflexus for multiple herbicides with standard errors in parentheses.^a

^a y = c + (d − c)/{1 + exp [b(logx − logED₅₀)]}, where b represents the relative slope, c represents the lower limit, d represents the upper limit, x represents the herbicide dose, and y represents the growth response (percentage of the untreated control); GR₅₀, dose required to reduce the plant dry weight by 50%.

^b R, resistant population HN-02F1-1; S, susceptible population ST-1.

^c Resistance index (RI) = GR₅₀ (HN-02F1-1)/GR₅₀ (ST-1).

Resistance to PPO inhibitors has been well represented in waterhemp [Amaranthus tuberculatus (Moq.) J. D. Sauer], the first weed species identified to have evolved resistance to these herbicides. The ΔGly210 mutation in the PPX2 gene provides A. tuberculatus with broad cross-resistance to most PPO inhibitors (Wuerffel et al. Reference Wuerffel, Young, Matthews and Young2015). The Arg-128 mutation was subsequently shown to confer cross-resistance against PPO inhibitors in both A. artemisiifolia and A. palmeri (Giacomini et al. Reference Giacomini, Umphres, Nie, Mueller, Steckel, Young, Scott and Tranel2017; Rousonelos et al. Reference Rousonelos, Lee, Moreira, VanGessel and Tranel2012). The Ala-205-Val substitution mutation is commonly considered to confer high IMI resistance to plants rather than resistance to SU herbicides (Ashigh and Tardif Reference Ashigh and Tardif2007). However, cross-resistance has also been reported to depend on specific mutations, on chemical groups of ALS inhibitors, and on specific herbicides within a particular group, and sometimes on weed species (Yu and Powles Reference Yu and Powles2014b). For instance, the Ala-205-Phe mutation in the ALS gene provides broad-spectrum resistance to ALS inhibitors in annual blugrass (Poa annua L.) (Brosnan et al. Reference Brosnan, Vargas, Breeden, Grier, Aponte, Tresch and Laforest2016). NTSR also plays an important role in ALS-inhibitor resistance. Délye et al. (Reference Délye, Pernin and Scarabel2011) reported that different corn poppy (Papaver rhoeas L.) individuals with identical ALS genotypes exhibited different levels of herbicide resistance, which was the result of NTSR. Therefore, in this study, the possibility of a metabolic resistance mechanism also cannot be ruled out in HN-02F1-1 with respect to the patterns of cross-resistance to PPO inhibitors and ALS inhibitors. The HN-02F1-1 population was still sensitive to flumioxazin, fluroxypyr-meptyl, and 2,4-D butylate. Therefore, flumioxazin, fluroxypyr-meptyl, and 2,4-D butylate can be used to prevent and manage the resistance of A. retroflexus, while crop rotation and/or mixtures of herbicides with different mechanisms of action as well as integrated pest management strategies are encouraged.

In Vitro Assays of ALS Activity

To confirm the ALS-inhibitor resistance related to the Ala-205-Val substitution of target resistance observed in the R population, ALS enzymes from the HN-02F1-1 and ST-1 populations were extracted and assayed. In the absence of ALS inhibitors, the HN-02F1-1 and ST-1 populations presented similar (RI = 1.1) activity of extracted ALS (17 nmol acetoin mg⁻¹ protein min⁻¹ and 15 nmol acetoin mg⁻¹ protein min⁻¹, respectively). Moreover, some of the same ALS inhibitors used for the whole-plant studies were also used for in vitro ALS activity assays in both R and S plants. As shown in Table 4, the I₅₀ of imazethapyr (IMI) for ST-1 was 0.86 μM, while it was 27 μM for HN-02F1-1. In addition, the I₅₀ of flumetsulam (TP) for ST-1 and HN-02F1-1 was 0.051 μM and 0.74 μM, respectively. The I₅₀ of nicosulfuron (SU) I₅₀ for the ST-1 population was 0.076 μM, while it was 2.0 μM for HN-02F1-1. On the basis of the RI value at I₅₀, compared with the ST-1 population, the HN-02F1-1 population was highly resistant to imazethapyr (320-fold), flumetsulam (15-fold), and nicosulfuron (28-fold). These results indicate that the Ala-205-Val substitution within ALS resulting in insensitivity provides cross-resistance to at least three types of ALS inhibitors (TP, SU, and IMI).

Table 4.

Parameter values of the four-parameter log-logistic equationfor calculating the I₅₀values of the in vitro ALS activity at different herbicide doses.^a

^a y = c + (d − c)/1 + exp [b (logx − logED₅₀)]}, where b represents the slope of the curve at I₅₀ (herbicide dose required to inhibit ALS activity by 50% compared with that of the untreated control), c represents the lower limit, d represents the upper limit, x represents the herbicide dose, and y represents the in vitro ALS activity response (as a percentage of the untreated control). The total acetolactate synthase (ALS) activity for HN-02F1-1 and ST-1 was 16.78 ± 0.53 nmol acetoin mg⁻¹ protein min⁻¹ and 15.11 ± 0.48 nmol acetoin mg⁻¹ protein min⁻¹, respectively.

^b R, resistant population HN-02F1-1; S, susceptible population ST-1.

^c Resistance index (RI) = I₅₀ (HN-02F1-1)/ I₅₀ (ST-1).

The results of the in vitro ALS activity inhibition assays of HN-02F1-1 plants were consistent with the results of the whole-plant herbicide bioassays. This finding further confirmed that the Ala-205-Val mutation plays an essential role in SU and IMI herbicide resistance, which is similar to the conclusions of other reports (Ashigh and Tardif Reference Ashigh and Tardif2006; Chen et al. Reference Chen, Huang, Zhang, Huang, Wei, Chen and Wang2015). Moreover, the HN-02F1-1 plants were highly resistant to TPs; this is the first reported case in which the Ala-205-Val mutation provides resistance to TP herbicides (Tranel et al. Reference Tranel, Wright and Heap2019). The reduced susceptibility of herbicide target enzymes in vitro typically translates into herbicide resistance in whole plants (Li et al. Reference Li, Li, Yu, Wang and Cui2017). Hence, the results of this study are consistent with those of previous research on SU and IMI resistance due to the Ala-205-Val mutation in horseweed (Erigeron canadensis L.) (Matzrafi et al. Reference Matzrafi, Lazar, Sibony and Rubin2015).

In most cases, amino acid substitutions in ALS provide resistance to ALS inhibitors (Shaner Reference Shaner1999). Here, we did not evaluate the role of the NTSR mechanism; thus, the possibility of NTSR cannot be excluded. In addition, the RIs did not completely correlate with the GR₅₀ and IC₅₀, which may due to a variety of reasons. NTSR may contribute to resistance in some HN-02F1-1 individuals, and the overexpression of the ALS gene was not ruled out. NTSR has recently been confirmed to be responsible for fomesafen resistance in A. palmeri populations from Arkansas (Varanasi et al. Reference Varanasi, Brabham and Norsworthy2018). Given that many weed species develop NTSR to ALS inhibitors (Zhao et al. Reference Zhao, Li, Bai, Guo, Yuan, Wang, Liu and Wang2017), it is reasonable to deduce that NTSR is present in A. retroflexus. Further investigations should focus on this aspect.

Amaranthus retroflexus is one of the main troublesome weed species in both soybean and cornfields (Wang et al. Reference Wang, Guo, Zhang, Zhao, Ge, Lv, Liu and Wang2017). In China, most farmers prefer to use one herbicide repeatedly as a single weed control method in a crop field. Additionally, farmers tend to increase the application rates of herbicides to ensure control efficacy. Given these factors, it is not surprising that the HN-02 population evolved multiple resistance to PPO inhibitors and ALS inhibitors. Evolution of severe herbicide resistance to photosystem II inhibitors, ALS inhibitors, PPO inhibitors, and glyphosate in A. tuberculatus populations in the U.S. Midwest has also been reported (Tranel et al. Reference Tranel, Riggins, Bell and Hager2010). Therefore, increased attention should be paid to investigating the resistance mechanism fully along with its fitness cost, biology, and ecology to control this notorious weed species.