INTRODUCTION
Agricultural weeds are mostly wild plant species, which thrive in agricultural ecosystems, competing for water, light and space with crops, and some are even intermediate hosts of diseases and pests. American slough grass (Beckmannia syzigachne (Steud.( Fernald) is adapted to the northern United States, southern Canada and part of eastern Asia (GBIF, 2017). In China, it is a troublesome weed which widely distributed in the wheat (Triticum aestivum L.) or oilseed rape (Brassica rapa L. subsp. oleifera (DC.) Metzg.) fields rotated with rice (Oryza sativa L.), mainly in the Yangtze Delta and the southwestern region (Li et al., 2013).
Herbicides are a major tool for controlling weeds to protect crop yields. Fenoxaprop-P-ethyl, belonging to the acetyl coenzyme A carboxylase (ACCase; EC 6.4.1.2) inhibitors, has long been used to selectively control grass weeds in wheat fields in China. However, the susceptibility of B. syzigachne to fenoxaprop-P-ethyl has decreased after many years use, and the resistance of B. syzigachne to fenoxaprop-P-ethyl was first reported in 2008 (Liu and Zhang, 2008). Consequently, mesosulfuron-methyl, a widely used acetolactate synthase (ALS; EC 4.1.3.18) inhibitor, rapidly became the most important herbicide for controlling ACCase inhibitor-resistant grass weeds (Bi et al., 2013). Mesosulfuron-methyl is well-known for its high efficiency, low animal toxicity worldwide and a broad herbicidal spectrum (Powles and Yu, 2010). However, resistance to ALS-inhibiting herbicides in B. syzigachne started to appear after several years of successful control (Li et al., 2015). Therefore, B. syzigachne became the most economically important crop weed, and the widespread evolution of multiple-herbicide resistance in this species makes its control more difficult.
There are generally two primary mechanisms in resistant weeds (Délye, 2005): target-site resistance (TSR) and non-target-site resistance (NTSR). TSR is due to changes in herbicide-binding sites or differences in the expression of target enzymes and has been well described (Délye, 2005; Powles and Yu, 2010). For example, the ACCase-inhibiting herbicides resistance is mainly caused by alterations of the ACCase enzyme. To date, 12 amino acid substitutions (AASs) at seven positions (i.e. 1781, 1999, 2027, 2041, 2078, 2088, and 2096) in the CT domain are implicated in conferring resistance to ACCase-inhibiting herbicides (Powles and Yu, 2010; Beckie and Tardif, 2012; Kaundun, 2014). Different ACCase mutations may confer distinct cross-resistance patterns to ACCase-inhibitors (Beckie and Tardif, 2012), and the homo/heterozygous status of the plants for the mutation can also influence the resistance levels (Délye, 2005; Scarabel et al., 2011). The substitutions at Leu1781 can confer resistance to aryloxyphenoxypropionates (APPs), cyclohexanediones (CHDs) and pinoxaden (Délye et al., 2002), 2027Cys and 2041Asn alleles can cause resistance to APPs and pinoxaden, but not to CHDs (Délye et al., 2003; Petit et al., 2010; Scarabel et al., 2011); 2096Ala allele confers resistance mainly to APPs (Délye, 2005). For the target site ALS gene mutations, at least 28 resistance endowing AASs at eight positions of the ALS gene have been identified over the past two decades (i.e. 122, 197, 205, 376, 377, 574, 645 and 653) (Yu and Powles, 2014a; Heap, 2017). Generally, substitutions at Ala122 confer resistance to Imidazolinones (IMIs) herbicides, but not to Sulfonylureas (SUs) herbicides; substitutions at Pro197 confer resistance to SU herbicides, but not to IMI herbicides; and substitutions at Trp574 provide resistance to both SU and IMI herbicides (Yu and Powles, 2014a).
NTSR is a complicated phenomenon compared with the TSR. It is thought to be caused by constitutive or herbicide-induced mechanisms and therefore mostly be driven by inheritable differences in the expression of a set of genes (Dinelli et al., 2006; Délye, 2013). To date, there are four well-established gene families participating in NTSR: the cytochrome P450 monooxygenases (P450s), glutathione S-transferases (GSTs), glycosyltransferases (GTs) and ATP-binding cassette (ABC) transporters (Yuan et al., 2007). The P450s constitute the largest family of enzymatic proteins in higher plants and protect plants from harmful exogenous chemicals, including pesticides and industrial pollutants, making them less phytotoxic. P450s are reported to play a primary role in phase I herbicide metabolism in weeds and crops (Werck-Reichhart et al., 2000; Siminszky, 2006), and has often been detected using specific P450 inhibitors, such as the organophosphate insecticide malathion, 1-aminobenzotriazole and piperonyl butoxide (PBO) (Werck-Reichhart et al., 2000). There were reports revealed that P450 inhibitors could enhance the rate of metabolism of several herbicides and had been confirmed in some herbicide-resistant weeds (Preston, 2004; Yun et al., 2005; Yasuor et al., 2009).
However, there were no reports about metabolic resistance in fenoxaprop-P-ethyl and mesosulfuron-methyl resistant B. syzigachne without any AASs in the target enzymes genes. The objective of this research was (1) to clone and compare the carboxyl transferase (CT) domain of the ACCase gene and the ALS gene in fenoxaprop-P-ethyl-resistant or sensitive and mesosulfuron-methyl-resistant or sensitive B. syzigachne, (2) to obtain the offspring of resistant population (RF1) and offspring of sensitive population (SF1) which sharing homogenized genetic backgrounds, (3) to investigate the effects of P450 inhibitors on the dose response to fenoxaprop-P-ethyl and mesosulfuron-methyl in the specific resistant B. syzigachne, (4) to establish more accurate cross resistance patterns at plant level including six ACCase inhibitors, three ALS inhibitors and one photosystem II inhibitor in the RF1.
MATERIALS AND METHODS
Plant material and herbicides
A suspected resistant population (R) was collected from wheat fields under a rice-wheat rotation field where fenoxaprop-P-ethyl and mesosulfuron-methyl had been used annually for several years in Anhui province, China. A suspected sensitive population (S) was collected from uncultivated land in Shandong province, China. To ensure the variability in field collection, each seed sample was collected from approximately 30 mature plants. Seed samples were air-dried and stored in sealed paper bags at 4 °C until used.
The herbicides used for dose response tests were: fenoxaprop-P-ethyl ((R)-2-[4-(6-chlorobenzoxazol-2-yloxy)phenoxy]propionic acid; 69 g L-1 EW, Bayer, Hangzhou, China); quizalofop-P-ethyl ((R)-2-[4-(6-chloroquinoxalin-2-yloxy)phenoxy]propionic acid; 50 g L-1 EC, LIER-Chemical, Mianyang, China); clodinafop-propargyl (prop-2-ynyl (R)-2-[4-(5-chloro-3-fluoro-2-pyridyloxy)phenoxy]propionate; 15% WP, Jinr, Shandong, China); sethoxydim ((5RS)-2-[(EZ)-1-(ethoxyimino)butyl]-5-[(2RS)-2-(ethylthio)propyl]-3-hydroxycyclohex-2-en-1-one 12.5% EC, Cynda, Jinan, China); clethodim ((5RS)-2-{(1EZ)-1-[(2E)-3-chloroallyloxyimino]propyl}-5-[(2RS)-2-(ethylthio)propyl]-3-hydroxycyclohex-2-en-1-one; 120 g L-1 EC, Cynda, Jinan, China); pinoxaden (8-(2,6-diethyl-p-tolyl)-1,2,4,5-tetrahydro-7-oxo-7H-pyrazolo[1,2-d][1,4,5]oxadiazepin-9-yl 2,2-dimethylpropionate; 5% EC, Syngenta, Suzhou, China); mesosulfuron-methyl (methyl 2-[(4,6-dimethoxypyrimidin-2-ylcarbamoyl)sulfamoyl]-α-(methanesulfonamido)-p-toluate; 30 g L-1 OF, Bayer, Hangzhou, China); flucarbazone-sodium (sodium [(4,5-dihydro-3-methoxy-4-methyl-5-oxo-1H-1,2,4-triazol-1-yl)carbonyl]{[2-(trifluoromethoxy)phenyl]sulfonyl}azanide; 70%, WG, Arysta, Shanghai, China); pyroxsulam (N-(5,7-dimethoxy[1,2,4]triazolo[1,5-a]pyrimidin-2-yl)-2-methoxy-4-(trifluoromethyl)pyridine-3-sulfonamide; 7.5% WG, Dow AgroSciences, Shanghai, China); isoproturon (3-p-cumenyl-1,1-dimethylurea; 50% WP, Bianjing, Suzhou, China), piperonyl butoxide (PBO) (2-(2-butoxyethoxy)ethyl 6-propylpiperonyl ether; 97%, Aladdin, Shanghai, China) and malathion (S-1,2-bis(ethoxycarbonyl)ethyl O,O-dimethyl phosphorodithioate; 95%, Ningbo Sunjoy Cropscience, Ningbo, China).
Acquisition of purified R F1 and S F1
Preliminary screening test. Prior to planting, seeds were first stored at -20 °C for 7 d and then deposited in Petri dishes containing two layers of filter paper soaked with about 5 mL distilled water. All of the Petri dishes were placed in chamber at 22/15 °C (12 h day/12 h night temperatures) under 75% RH. After germination, the seeds were sown in plastic pots and transferred to a glasshouse (temperature maintained at approximately 15-25 °C, with 63%-75% RH and natural sunlight) (Li et al., 2013). Seedlings were thinned to six evenly sized plants per pot and watered as needed. At the three to four leaf stage, seedlings of suspected R population were sprayed with fenoxaprop-P-ethyl at 0, 62.1 (1×), 186.3 (3×) and 558.9 (9×) g ai ha-1; mesosulfuron-methyl at 0, 5.25 (1/3×), 15.75 (1×) and 47.25 (3×) g ai ha-1. While the seedlings of suspected S population were sprayed with fenoxaprop-P-ethyl at 0, 20.7 (1/3×), 62.1 (1×) and 186.3 (3×) g ai ha-1; mesosulfuron-methyl at 0, 1.75 (1/9×), 5.25 (1/3×) and 15.75 (1×) g ai ha-1. Spraying was performed using a compressed air, moving nozzle cabinet sprayer equipped with one flat-fan nozzle (9503EVS, TeeJet, Wheaton, Illinois, USA) at a height of 50 cm above the foliage, calibrated to deliver 450 L ha-1 at a pressure of 0.28 MPa. The treated plants were returned to the greenhouse and randomly placed on raised beds. At 21 d after treatment (DAT), the result was checked visually.
Amplification, cloning and sequencing of the target enzyme genes. Seeds of both populations were treated as described above, with three seedlings per pot, and each population containing four pots. Forty milligrams of shoot tissue from individual plants in each pot was mixed and frozen in liquid nitrogen for DNA extraction. Total DNA was extracted using the cetyltrimethylammonium bromide (CTAB) method (Doyle, 1987). As this harvest only included a small part of the plant and the remaining plant tissue was left and transplanted into larger plastic pots to produce seeds that would later be used to produce the RF1 and SF1 generation. When they were ripe, seeds from the two populations (RF1 and SF1) were harvested and stored at room temperature for 2-mo to enable after-ripening.
To identify the molecular mechanism of resistance to fenoxaprop-P-ethyl in B. syzigachne, the primers and procedure in Li et al. (2013) were used to amplify the carboxyltransferase (CT) domain of the ACCase gene of B. syzigachne by polymerase chain reaction (PCR). In addition, the primers and procedure in Li et al. (2015) were used, which amplified a fragment including all known eight ALS inhibitor-resistant positions. The amplified DNA fragments were purified using the TIANgel Midi Purification Kit and cloned directly into pEASY-T1 vector (TransGen Biotech, Beijing, China). The sequences of the inserts were determined using the ABI PRISM 3730 DNA sequencer (Sangon Biotech, Shanghai, China). Twelve individual plants from the R population and six plants from the S population were sequenced and at least four positive clones from each biological replicate were sent for sequencing. Each fragment was sequenced in the forward and reverse directions by Shanghai Sangon Biological Engineering Technology & Services Co. (Shanghai, China). The sequencing data from the two B. syzigachne populations were compared to determine whether there were AASs associated with resistance using DNAMAN version 5.2.2 software (Lynnon Biosoft, Quebec, Canada).
To ensure the background of RF1 and SF1, we randomly selected 60 plants from the RF1 and 30 plants from the SF1 sequenced as described above.
Sensitivity to fenoxaprop- P -ethyl and mesosulfuron-methyl following the P450s inhibitors
The seeds of RF1 and SF1 were pretreated as described above and then six seedlings were sown per pot. Seedlings at the three to four leaf stage were sprayed with PBO, fenoxaprop-P-ethyl, PBO plus fenoxaprop-P-ethyl, malathion, mesosulfuron-methyl and malathion plus mesosulfuron-methyl, in addition, we set up three concentration gradients of each P450 inhibitors (Table 1). There were a total 35 different treatments of each P450 inhibitors, including the blank control (spraying with the same volume of water). Spraying was performed as described above. At 21 d after treatment (DAT), the aboveground materials were harvested and the fresh-weight data were recorded.
Table 1: Piperonyl butoxide (PBO) and malathion treatment applied for P450 inhibitors.

Underlined values are the recommended doses.
RF1: Offspring of resistant Beckmannia syzigachne population; SF1: offspring of sensitive B. syzigachne population.
PBO was applied in two applications of 2100 g ai ha-1 in 97 L ha-1 water, for total application of 4200 g ai ha-1 in 194 L ha-1 water. Malathion and PBO were formulated in a mixture of emulsifier (Tween-80, 1 mL L-1 water) and acetone and were applied 1 h prior to herbicide application (Preston et al., 1996). Reference plants were treated with a mixture of emulsifier and acetone. All treatments were replicated three times, and the experiment was conducted twice.
Sensitivity to other herbicides
Seedlings of RF1 and SF1 at the three to four leaf stage were sprayed with 10 herbicides at different rates (Table 2). Three weeks later, the results were checked as described above. All treatments were replicated three times, and the experiment was conducted twice. Resistance to six ACCase-inhibiting herbicides, including three aryloxyphenoxypropionates (APPs), two cyclohexanediones (CHDs), and one phenylpyraxoline (DEN), at different doses was determined. Three ALS inhibitor and one photosystem II inhibitor were also tested in this experiment.
Statistical analyses
The results were expressed as percentage inhibition [(fresh weight of treated pots/fresh weight of control pots) × 100]. ANOVA was selected to determine whether results were significantly different between the data of the two runs of the experiment for each herbicide and was performed on all data using SPSS software (Version 20.0; IBM Corporation, Armonk, New York, USA). When the ANOVA showed nonsignificant difference between the two runs of the experiment, thus the data were pooled and used for subsequent analysis.
The herbicide rate that inhibited growth by 50% (GR50) was estimated based on a four-parameter non-linear logistic model (sigmoidal logistic, four parameters) using SigmaPlot software (SigmaPlot v.12.0; Systat Software, San Jose, California, USA) (Hochberg et al., 2009):
Y = C + (D - C)/[1 + (x/GR50)b]
In this model, C is the lower limit, D is the upper limit, GR50 is the effective dose causing 50% growth reduction, and b describes the slope of the GR50 curve. The fitted equation was used to estimate the GR50 value. The resistance index (RI) value was calculated by dividing the GR50 value of the resistant population by that of the sensitive population.
RESULTS
Preliminary screening test
The suspected R population could grow normally even at three times of the recommend dose of fenoxaprop-P-ethyl, while all the suspect S population dead at the recommend dose. For the mesosulfuron-methyl test, the suspect R individuals’ growth was inhibited at the recommend dose, but all are survived; while for the suspect S population all dead at the recommend dose. The result of preliminary screening test showed that suspect R population did resistance to fenoxaprop-P-ethyl and mesosulfuron-methyl.
Sequencing and alignment of the target genes
A 1437 bp fragment including all seven ACCase inhibitor-resistant positions within the CT domain of ACCase gene was amplified, sequenced and blasted between the two populations and also between the RF1 and SF1. A 1729 bp fragment encompassing all eight ALS inhibitor-resistant positions within the partial of the ALS gene was amplified to compare sequences between the two populations and also between the RF1 and SF1. The blasted results indicated that there were not any reported AASs about resistance in the target enzymes genes in the R population as well as the RF1 (Figures 1 and 2).

ACCase: Acetyl coenzyme A carboxylase; ACCase-R: CT domain of resistant population’s ACCase gene; ACCase-S: CT domain of sensitive population’s ACCase gene; RF1: offspring of resistant Beckmannia syzigachne population; SF1: offspring of sensitive B. syzigachne population.
Figure 1: Alignment of Beckmannia syzigachne ACCase with Alopecurus myosuroides. Boxes indicate amino acid positions whose substitutions are known to confer ACCase inhibitor resistance.

ALS: Acetolactate synthase; ALS-R: acetolactate synthase of resistant population; ALS-S: acetolactate synthase of sensitive population; RF1: offspring of resistant Beckmannia syzigachne population; SF1: offspring of sensitive B. syzigachne population.
Figure 2: Alignment of Beckmannia syzigachne ALS with Arabidopsis thaliana. Boxes indicate amino acid positions whose substitutions are known to confer ALS inhibitor resistance.
Sensitivity to fenoxaprop- P -ethyl and mesosulfuron-methyl following the P450s inhibitors
In the absence of fenoxaprop-P-ethyl, treatment with 4200 g ai ha-1 PBO had no obvious effect on the RF1 and SF1, nor did malathion (2000 g ai ha-1). Plants from the RF1 treated with fenoxaprop-P-ethyl plus PBO exhibited significantly increased mortality compared with fenoxaprop-P-ethyl treatment alone. Pretreatment with the maximum concentration of PBO significantly inhibited growth of the RF1, decreased the GR50 value by 60.2% compared with fenoxaprop-P-ethyl-only treatment. The malathion combined with mesosulfuron-methyl increased toxicity in plants, and pretreatment with malathion decreased the GR50 value by 65.6% compared with the mesosulfuron-methyl-only treatment in the RF1. Moreover, we observed that mesosulfuron-methyl plus malathion or fenoxaprop-P-ethyl plus PBO also reduced the GR50 of the SF1 partly compared with herbicide treatment alone (Table 3 and Figure 3).
Table 3: Response to herbicides following different concentrations of P450 inhibitors.

GR50: Herbicide ratio required to decrease plant fresh weight by 50% compared to the untreated control. Each value represents the mean ± standard error; RI: resistance index was calculated as the ratio between the GR50 of the resistant population and the GR50 of the susceptible population; RF1: offspring of Beckmannia syzigachne resistant population; SF1: offspring of sensitive B. syzigachne population.

RF1: Offspring of resistant Beckmannia syzigachne population; SF1: offspring of sensitive B. syzigachne population; PBO: piperonyl butoxide.
Figure 3: Dose-response curve for above ground fresh weights of RF1 and SF1 plants treated with fenoxaprop-P-ethyl and mesosulfuron-methyl following the P450s inhibitors. Each point represents the mean of two experiments, each containing three replicates.
Levels of resistance to herbicides in R F1
The results of whole-plant dose-response experiments showed that the RF1 produced high resistance to fenoxaprop-P-ethyl, based on the RI of GR50 which RI was as high as 31.2 (Table 4 and Figure 4) (Beckie and Tardif, 2012). The results showed that the RF1 exhibited low resistance to quizalofop-P-ethyl, clodinafop-propargyl and sethoxydim, with RI values of 2.8, 4.2 and 2.7, respectively. The RF1 displayed no resistance to clethodim and pinoxaden, with corresponding RIs below 2. Moreover, the RF1 showed multiple resistance to ALS inhibitors, but not including photosystem II inhibitors. The RF1 showed low resistance to mesosulfuron-methyl, intermediate resistance to flucarbazone-sodium, and high resistance to pyroxsulam, with RI values of 2.6, 7.7 and 23.6, respectively. The GR50 of RF1 to isoproturon was 63.3 g ai ha-1, which value was much lower than the recommended rate (1050 g ai ha-1) showing its high efficiency in controlling the multiple-resistant B. syzigachne.
Table 4: The sensitivities of resistant and sensitive populations to other herbicides.

GR50: Herbicide ratio required to decrease plant fresh weight by 50% compared to the untreated control. Each value represents the mean ± standard error; RI: resistance index was calculated as the ratio between the GR50 of the resistant population and the GR50 of the susceptible population; RF1: offspring of resistant Beckmannia syzigachne population; SF1: offspring of sensitive B. syzigachne population.

RF1: Offspring of resistant Beckmannia syzigachne population; SF1: offspring of sensitive B. syzigachne population.
Figure 4: Dose-response curve for above ground fresh weights of RF1 and SF1 plants treated with 10 herbicides. Each point represents the mean of two experiments, each containing three replicates.
DISCUSSION
Metabolic resistance in weeds is much more complicated compared with TSR and due to the threat of NTSR to agriculture, extensive investigation of the underlying mechanism is necessary (Yu and Powles, 2014b). The population examined was shown to be resistant to fenoxaprop-P-ethyl and mesosulfuron-methyl though the preliminary screening test. The result of sequencing and alignment of the target genes showed that there were not any endow resistance AASs in the resistant population, neither in the RF1. From above phenomenon, we suspect that NTSR is present in this specific resistant population. Therefore, we added an experiment using P450 inhibitors to verify this hypothesis and did a preliminary exploration about this kind of resistant B. syzigachne.
The result in our research also suggested that synergism may occur and there was an effect (particularly for the RF1) when the PBO and malathion were added, which in accordance with a previous report (Pan et al., 2015). Moreover, we observed that mesosulfuron-methyl plus malathion or fenoxaprop-P-ethyl plus PBO also reduced the GR50 of the SF1 partly compared with herbicide treatment alone, probably at rates not sufficient to confer resistance though SF1 metabolizes some mesosulfuron-methyl or fenoxaprop-P-ethyl (Gaines et al., 2014). Based on above results, P450-based metabolism is most likely involved in the RF1. There were several inhibitors of P450s had been identified, including 1-aminobenzo-triazole, tetcyclacis, PBO, carbon monoxide and tridiphane (McFadden et al., 1989; Zimmerlin and Durst, 1990; Siminszky, 2006; Pan et al., 2015). These inhibitors have provided evidence of the involvement of P450s in herbicide resistance and tolerance. Application of PBO to Poa annua in vivo showed that PBO exerted a synergistic effect with fenoxaprop-P-ethyl (Wang et al., 2013). Malathion has shown to present synergistic activity with the ALS-inhibiting herbicide chlorsulfuron against the resistant Lolium rigidum biotype (Preston et al., 1996). Similar effects of malathion on the metabolism of the ALS inhibitor propoxycarbazone-sodium have been reported for Bromus tectorum (Park et al., 2004).
P450-mediated mechanism can trigger resistance to herbicides with similar mechanisms or different mechanisms of action, including herbicides that have never been used. To obtain a more comprehensive understanding of the background of the RF1, we characterized the cross and multiple resistance patterns to 10 herbicides at plant level. The evolution of multiple resistance to ACCase- and ALS-inhibitors in the RF1 indicates a further threat to alternative herbicides for controlling the weed. Fortunately, there was no resistance to isoproturon, clethodim emulsifiable and pinoxaden, and then our results indicate that these four herbicides should be used. When managing this kind of resistant B. syzigachne with metabolic resistance, an integrated system should be used that includes rotation of crops and tillage systems as well as using a fully diversified herbicide program.
The results in this research were consistent with previous work showing that NTSR generally confers resistance to group of herbicides, or general for some herbicide groups (Yu and Powles, 2014b). The RF1 in this study was resistant to fenoxaprop-P-ethyl and mesosulfuron-methyl but without finding any AASs in the target enzymes. Therefore, we speculated that this might result from NTSR and/or other unknown mechanism and this inconclusive resistance mechanism may confer resistance in the RF1 and further analysis is in progress.
CONCLUSION
In conclusion, the two target enzymes were cloned and compared between the two Beckmannia syzigachne populations, resistant (RF1) and sensitive (SF1). The effects of P450 inhibitors on the dose response to fenoxaprop-P-ethyl and mesosulfuron-methyl in the RF1 were explored, as well as the comprehensive cross-resistance patterns. Further study is needed to determine whether resistance is associated with the expression level of the target enzymes. It is also possible that mutations in the upstream regulatory regions of the target genes influence resistance. Accordingly, these possibilities will be investigated in future studies.