RNA-Seq transcriptome analysis to identify genes involved in metabolism-based diclofop resistance in Lolium rigidum

Mechanisms of evolved herbicide resistance

The widely successful use of synthetic herbicides over the past 70 years has imposed strong and widespread selection pressure, leading to the evolution of herbicide resistance in hundreds of weed species. Both target-site resistance (TSR) and nontarget-site resistance (NTSR) mechanisms have evolved to most herbicide classes. TSR often involves mutations in genes encoding the protein targets of herbicides, affecting the binding of the herbicide either at or near catalytic domains or in regions affecting access to them. Most of these mutations are nonsynonymous SNPs, but polymorphisms in more than one codon or entire codon deletions have also evolved. Some herbicides bind multiple proteins, making the evolution of TSR mechanisms more difficult. Increased amounts of protein target, by increased gene expression or by gene duplication, are an important, albeit less common, TSR mechanism. NTSR mechanisms include reduced absorption or translocation and increased sequestration or metabolic degradation. The mechanisms that can contribute to NTSR are complex and often involve genes that are members of large gene families. For example, enzymes involved in herbicide metabolism-based resistances include cytochromes P450, GSH S-transferases, glucosyl and other transferases, aryl acylamidase, and others. Both TSR and NTSR mechanisms can combine at the individual level to produce higher resistance levels. The vast array of herbicide-resistance mechanisms for generalist (NTSR) and specialist (TSR and some NTSR) adaptations that have evolved over a few decades illustrate the evolutionary resilience of weed populations to extreme selection pressures. These evolutionary processes drive herbicide and herbicide-resistant crop development and resistance management strategies.

Non-target-site resistance to PDS-inhibiting herbicides in a wild radish (Raphanus raphanistrum) population

BACKGROUND

Diflufenican resistance has been reported in wild radish populations since 1998, but the resistance mechanisms have not been investigated. Recently, we identified a wild radish population (H2/10) from the Western Australian grain belt that is resistant (R) to the phytoene desaturase (PDS)-inhibiting herbicide diflufenican.

RESULTS

Dose-response results showed this R population is 4.9-fold more resistant than the susceptible (S) population based on the LD₅₀ R/S ratio. In addition, the R population also exhibits cross-resistance to the PDS-inhibiting herbicide fluridone. The cytochrome P450 inhibitor malathion reversed diflufenican resistance and partially reversed fluridone resistance in the R population. The full coding sequences of the PDS gene were cloned from the S and R plants and there are natural variations in the PDS gene transcripts/alleles with no correlation to resistance. In addition, the R plants had a level of PDS gene expression that is not significantly different from the S plants.

CONCLUSION

These results demonstrated that diflufenican resistance in this R wild radish population is likely due to non-target-site based enhanced herbicide metabolism involving cytochrome P450s. © 2019 Society of Chemical Industry.

Evolution of resistance to HPPD-inhibiting herbicides in a wild radish population via enhanced herbicide metabolism

BACKGROUND

Relatively new herbicides that target 4-hydroxyphenylpyruvate dioxygenase (HPPD) are now available for use on the world's great grain crops (rice, wheat, corn and soybean) and for other uses. With widespread and persistent use of HPPD-inhibiting herbicides, the evolution of HPPD-inhibiting herbicide resistant weeds is inevitable. Currently, resistance to HPPD-inhibiting herbicides is known in two weed species, waterhemp and Palmer amaranth. Here, we report a HPPD-inhibiting herbicide resistant wild radish population from the Western Australia grain belt. This population was not selected with HPPD-inhibiting herbicides, rather it evolved resistance to earlier used herbicides with different modes of action and exhibits cross-resistance to HPPD-inhibiting herbicides.

RESULTS

Dose-response experiments showed the resistant (R) population exhibits 4 to 6.5-fold resistance to the HPPD-inhibiting herbicides mesotrione, tembotrione and isoxaflutole, compared to the susceptible (S) population. This resistance is not target-site based as cloning of full coding sequences of the HPPD genes from S and R plants did not reveal resistance-endowing single nucleotide polymorphisms. The HPPD gene expression levels are similar in S and R plants. In addition, no differences in [¹⁴ C]-mesotrione uptake and translocation were observed in the S and R plants. However, the time required for R plants to metabolise 50% [¹⁴ C]-mesotrione is 7.7-fold faster than for the S plants.

CONCLUSION

We confirm resistance to HPPD-inhibiting herbicides exists in a population of the economically damaging global weed wild radish. The resistance in this population is due to a non-target-site based enhanced rate of herbicide metabolism. © 2019 Society of Chemical Industry.

Comparative structural and functional analysis of phi class glutathione transferases involved in multiple-herbicide resistance of grass weeds and crops

Multiple-herbicide resistant (MHR) weeds are a global problem and a looming threat to weed control in crops. MHR weeds express a specific phi class glutathione transferase (MHR-GSTF) which seems to contribute to herbicide resistance. The present work aims to investigate the structure and catalytic properties of the MHR-GSTFs from different grass weeds and crops (Alopecurus myosuroides, Lolium rigidum, Hordeum vulgare, Triticum aestivum). Recombinant MHR-GSTFs were expressed in E. coli and purified by affinity chromatography. Kinetic analysis of substrate specificity using a range of thiol substrates and xenobiotic compounds suggested that all enzymes display a broad range of specificity and are capable of detoxifying major stress-induced toxic products. Notably, all tested enzymes exhibited high activity towards organic hydroperoxides. The crystal structure of MHR-GSTF from Alopecurus myosuroides (AmGSTF) was determined by molecular replacement at 1.33 Å resolution. The enzyme was resolved with bound glutathione sulfenic acid (GSOH) at the G-site and succinic acid at the H-site. The enzyme shows conserved structural features compared to other Phi class GSTs. However, some differences were observed at the C-terminal helix H9 that may affect substrate specificity. The structural and functional features of AmGSTF were compared with those of the homologue crop enzymes (HvGSTF and TaGSTF) and discussed in light of their contribution to the MHR mechanism.

Cyhalofop-butyl and Glyphosate Multiple-Herbicide Resistance Evolved in an Eleusine indica Population Collected in Chinese Direct-Seeding Rice

Eleusine indica is a typical xerophytic weed species with a cosmopolitan distribution. It is invasive and highly adaptable to diverse habitats and crops. Due to rice cropping-pattern changes, E indica has become one of the main dominant grass weeds infecting direct-seeding paddy fields. A Chinese E. indica population has evolved multiple-herbicide resistance to cyhalofop-butyl and glyphosate. In this study, the multiple-resistance profile of E. indica to these two different types of herbicides and their resistance mechanisms were investigated. Whole-plant dose-response assays indicated that the multiple-herbicide-resistant (MHR) population exhibited 10.8-fold resistance to cyhalofop-butyl and 3.1-fold resistance to glyphosate compared with the susceptible (S) population. ACCase sequencing revealed that the Asp-2078-Gly mutation was strongly associated with E. indica resistance to cyhalofop-butyl. The MHR plants accumulated less shikimic acid than S plants at 4, 6, and 8 days after glyphosate treatment. In addition, no amino acid substitution in the EPSPS gene was found in MHR plants. Further analysis revealed that the relative expression level of EPSPS in MHR plants was 6-10-fold higher than that in S plants following glyphosate treatment, indicating that EPSPS overexpression may contribute to the glyphosate resistance. Furthermore, the effectiveness of nine post-emergence herbicides against E. indica were evaluated, and one PPO inhibitor pyraclonil was identified as highly effective in controlling the S and MHR E. indica populations.

Herbicide drift exposure leads to reduced herbicide sensitivity in Amaranthus spp

While the introduction of herbicide tolerant crops provided growers new options to manage weeds, the widespread adoption of these herbicides increased the risk for herbicide spray drift to surrounding vegetation. The impact of herbicide drift in sensitive crops is extensively investigated, whereas scarce information is available on the consequences of herbicide drift in non-target plants. Weeds are often abundant in field margins and ditches surrounding agricultural landscapes. Repeated herbicide drift exposure to weeds could be detrimental to long-term management as numerous weeds evolved herbicide resistance following recurrent-selection with low herbicide rates. The objective of this study was to evaluate if glyphosate, 2,4-D, and dicamba spray drift could select Amaranthus spp. biotypes with reduced herbicide sensitivity. Palmer amaranth and waterhemp populations were recurrently exposed to herbicide drift in a wind tunnel study over two generations. Seeds from survival plants were used for the subsequent rounds of herbicide drift exposure. Progenies were subjected to herbicide dose-response studies following drift selection. Herbicide drift exposure rapidly selected for Amaranthus spp. biotypes with reduced herbicide sensitivity over two generations. Weed management programs should consider strategies to mitigate near-field spray drift and suppress the establishment of resistance-prone weeds on field borders and ditches in agricultural landscapes.

A transcriptomic analysis reveals soybean seed pre-harvest deterioration resistance pathways under high temperature and humidity stress

Pre-harvest soybean seeds in the field are susceptible to high temperature and humidity (HTH) stress, leading to pre-harvest seed deterioration, which will result in a reduction in grain quality, yield, and seed vigor. To understand the gene expression involved in seed deterioration response under HTH stress, in this study, we conducted an RNA-Seq analysis using two previously screened soybean cultivars with contrasting seed deterioration resistance. HTH stress induced 1081 and 357 differentially expressed genes (DEGs) in the sensitive cultivar Ningzhen No. 1 and resistant cultivar Xiangdou No. 3, respectively. The majority of DEGs in the resistant cultivar were up-regulated, while down-regulated DEGs were predominant in the sensitive cultivar. KEGG pathway analysis revealed that metabolic pathways, biosynthesis of secondary metabolites, and protein processing in endoplasmic reticulum were the predominant pathways in both cultivars during seed deterioration under HTH stress. The genes involved in photosynthesis, carbohydrate metabolism, lipid metabolism, and heat shock proteins pathways might contribute to the different response to seed deterioration under HTH treatment in the two soybean cultivars. Our study extends the knowledge of gene expression in soybean seed under HTH stress and further provides insight into the molecular mechanism of seed deterioration as well as new strategies for breeding soybean with improved seed deterioration resistance.

A Val-202-Phe α-tubulin mutation and enhanced metabolism confer dinitroaniline resistance in a single Lolium rigidum population

BACKGROUND

A Lolium rigidum population collected from Western Australia was previously reported as highly resistant to dinitroaniline herbicides mainly due to a Val-202-Phe substitution in the target site α-tubulin protein. To further determine the contribution of the 202 mutation to resistance, two sub-populations, respectively comprising the 202 mutant and wild-type (WT) individuals, were isolated from within the same resistant population and subject to dinitroaniline herbicide doses. A rice transgenic study was conducted to demonstrate whether the amino acid substitution at the 202 residue confers resistance. In addition, as indicated in the phenotyping and genotyping study, non-target enhanced trifluralin metabolism was further examined in the same population.

RESULTS

The 202 mutants were more resistant than the wild-type plants. Rice calli transformed with the L. rigidum mutant α-tubulin gene (Val-202-Phe) were more resistant to dinitroaniline herbicides relative to calli transformed with the wild-type gene. Also, enhanced trifluralin metabolism was detected in the 202 mutants in comparison to the susceptible seedlings.

CONLCUSION

Both target-site Val-202-Phe α-tubulin mutation and non-target-site enhanced trifluralin metabolism co-exist in this dinitroaniline-resistant L. rigidum population. © 2019 Society of Chemical Industry.

Herbicide resistant weeds: A call to integrate conventional agricultural practices, molecular biology knowledge and new technologies

Herbicide resistant (HR) weeds are of major concern in modern agriculture. This situation is exacerbated by the massive adoption of herbicide-based technologies along with the overuse of a few active ingredients to control weeds over vast areas year after year. Also, many other anthropological, biological, and environmental factors have defined a higher rate of herbicide resistance evolution in numerous weed species around the world. This review focuses on two central points: 1) how these factors have affected the resistance evolution process; and 2) which cultural practices and new approaches would help to achieve an effective integrated weed management. We claim that global climate change is an unnoticed factor that may be acting on the selection of HR weeds, especially those evolving into non-target-site resistance mechanisms. And we present several new tools -such as Gene Drive and RNAi technologies- that may be adopted to cope with herbicide resistance spread, as well as discuss their potential application at field level. This is the first review that integrates agronomic and molecular knowledge of herbicide resistance. It covers not only the genetic basis of the most relevant resistance mechanisms but also the strengths and weaknesses of traditional and forthcoming agricultural practices.

Cytochrome P450 Herbicide Metabolism as the Main Mechanism of Cross-Resistance to ACCase- and ALS-Inhibitors in Lolium spp. Populations From Argentina: A Molecular Approach in Characterization and Detection

Knowledge about the mechanisms of herbicide resistance provide valuable insights into evolving weed populations in response to selection pressure and should be used as a basis for designing management strategies for herbicide-resistant weeds. The selection pressure associated with reactive management against glyphosate-resistant Lolium spp. populations would have favored the herbicide resistance to ACCase- and ALS-inhibitors. This work was aimed to determine the sensitivity of 80 Argentinean Lolium spp. populations to ALS- and ACCase-inhibitor herbicides for use in wheat or barley and to study the mechanisms of resistance involved. Sensitivity to pinoxaden and iodosulfuron-mesosulfuron were positively correlated (r = 0.84), even though both affect different target sites. Inhibitors of cytochrome P450 monooxygenases (P450s) increased the sensitivity to pinoxaden and iodosulfuron-mesosulfuron in 94% of herbicide-resistant populations and target-site ACCase resistance mutations were detected only in two cases. Polymorphic variants were obtained with a pair primer designed on P450 sequences, cluster analysis discriminated around 80% of susceptible and P450-metabolic resistant plants sampled from a single population or different populations. Five markers corresponding to herbicide sensitivity were identified to be significantly associated with phenotypic variance in plants. Resistance to ALS- and ACCase-inhibitor herbicides were closely related, challenging the rotation of herbicides of both sites of action as a practice against resistance. In that sense, the use of pinoxaden and iodosulfuron-mesosulfuron would have provoked a selection on P450 genes that conduced a convergence of P450-metabolism based resistant Lolium spp. populations, which was detected by markers in a contribution to elucidate the molecular basis of this type of resistance.

Omics Potential in Herbicide-Resistant Weed Management

The rapid development of omics technologies has drastically altered the way biologists conduct research. Basic plant biology and genomics have incorporated these technologies, while some challenges remain for use in applied biology. Weed science, on the whole, is still learning how to integrate omics technologies into the discipline; however, omics techniques are more frequently being implemented in new and creative ways to address basic questions in weed biology as well as the more practical questions of improving weed management. This has been especially true in the subdiscipline of herbicide resistance where important questions are the evolution and genetic basis of herbicide resistance. This review examines the advantages, challenges, potential solutions, and outlook for omics technologies in the discipline of weed science, with examples of how omics technologies will impact herbicide resistance studies and ultimately improve management of herbicide-resistant populations.

Trp-1999-Ser mutation of acetyl-CoA carboxylase and cytochrome P450s-involved metabolism confer resistance to fenoxaprop-P-ethyl in Polypogon fugax

BACKGROUND

Asia minor bluegrass (Polypogon fugax Nees ex Steud.) is an invasive grass species severely infesting wheat and canola fields in China. In May 2017, a suspected resistant P. fugax population AHHY that survived fenoxaprop-P-ethyl applied at its field-recommended rate was collected from a wheat field in Huaiyuan County, Anhui Province, China. This study aimed to determine the resistance profile of AHHY to acetyl-CoA carboxylase (ACCase) inhibitors and to investigate its mechanisms of resistance to fenoxaprop.

RESULTS

Single-dose testing indicated that the AHHY population had evolved resistance to fenoxaprop. The partial carboxyltransferase domain of ACCase in P. fugax was amplified and compared. Four loci encoding plastidic ACCase were isolated from both the resistant and sensitive individuals. Combining gene sequencing with the derived cleaved amplified polymorphic sequence assay, we found that 100% of the plants of AHHY carried Trp-1999-Ser mutation in their ACCase1,1-2 allele. Whole-plant dose-response bioassay indicated that AHHY was highly resistant to fenoxaprop and pinoxaden (resistance index (RI) ≥ 10) with low resistance to clodinafop-propargyl, sethoxydim, and clethodim (2 ≤ RI < 5). Pre-treatment with piperonyl butoxide largely reduced (55%) the weed's resistance to fenoxaprop. Both basal and fenoxaprop-induced glutathione S-transferases activities toward 1-chloro-2, 4-dinitrobenzene were significantly higher in resistant plants than in susceptible plants.

CONCLUSION

This study revealed that P. fugax had multiple alleles encoding plastidic ACCase, and reported for the first time the occurrence of Trp-1999-Ser mutation and non-target-site resistance in this species. Fenoxaprop resistance in AHHY plants was conferred by target-site mutation and P450s-involved enhanced metabolism. © 2019 Society of Chemical Industry.

Non-Target-Site Resistance to Herbicides: Recent Developments

Non-target-site resistance (NTSR) to herbicides in weeds can be conferred as a result of the alteration of one or more physiological processes, including herbicide absorption, translocation, sequestration, and metabolism. The mechanisms of NTSR are generally more complex to decipher than target-site resistance (TSR) and can impart cross-resistance to herbicides with different modes of action. Metabolism-based NTSR has been reported in many agriculturally important weeds, although reduced translocation and sequestration of herbicides has also been found in some weeds. This review focuses on summarizing the recent advances in our understanding of the physiological, biochemical, and molecular basis of NTSR mechanisms found in weed species. Further, the importance of examining the co-existence of TSR and NTSR for the same herbicide in the same weed species and influence of environmental conditions in the altering and selection of NTSR is also discussed. Knowledge of the prevalence of NTSR mechanisms and co-existing TSR and NTSR in weeds is crucial for designing sustainable weed management strategies to discourage the further evolution and selection of herbicide resistance in weeds.

Response of Amaranthus spp. following exposure to sublethal herbicide rates via spray particle drift

The adverse consequences of herbicide drift towards sensitive crops have been extensively reported in the literature. However, little to no information is available on the consequences of herbicide drift onto weed species inhabiting boundaries of agricultural fields. Exposure to herbicide drift could be detrimental to long-term weed management as several weed species have evolved herbicide-resistance after recurrent selection with sublethal herbicide rates This study investigated the deposition of glyphosate, 2,4-D, and dicamba spray particle drift from applications with two different nozzles in a low speed wind tunnel, and their impact on growth and development of Amaranthus spp. Herbicide drift resulted in biomass reduction or complete plant mortality. Inflection points (distance to 50% biomass reduction) for Amaranthus tuberculatus were 7.7, 4.0, and 4.1 m downwind distance for glyphosate, 2,4-D, and dicamba applications with the flat-fan nozzle, respectively, whereas these values corresponded to 2.8, 2.5, and 1.9 m for applications with the air-inclusion nozzle. Inflection points for Amaranthus palmeri biomass reduction were 16.3, 10.9, and 11.5 m for glyphosate, 2,4-D, and dicamba applications with the flat-fan nozzle, respectively, whereas these values corresponded to 7.6, 5.4, and 5.4 m for applications with the air-inclusion nozzle. Plants were more sensitive to glyphosate at higher exposure rates than other herbicides, whereas plants were more sensitive to 2,4-D and dicamba at lower exposure rates compared to glyphosate. Applications with the flat-fan nozzle resulted in 32.3 and 11.5% drift of the applied rate at 1.0 and 3.0 m downwind, respectively, whereas the air-inclusion nozzle decreased the dose exposure in the same distances (11.4 and 2.7%, respectively). Herbicide drift towards field boundaries was influenced by nozzle design and exposed weeds to herbicide rates previously reported to select for herbicide-resistant biotypes.

Do plants pay a fitness cost to be resistant to glyphosate?

We reviewed the literature to understand the effects of glyphosate resistance on plant fitness at the molecular, biochemical and physiological levels. A number of correlations between enzyme characteristics and glyphosate resistance imply the existence of a plant fitness cost associated with resistance-conferring mutations in the glyphosate target enzyme, 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS). These biochemical changes result in a tradeoff between the glyphosate resistance of the EPSPS enzyme and its catalytic activity. Mutations that endow the highest resistance are more likely to decrease catalytic activity by reducing the affinity of EPSPS for its natural substrate, and/or slowing the velocity of the enzyme reaction, and are thus very likely to endow a substantial plant fitness cost. Prediction of fitness costs associated with EPSPS gene amplification and overexpression can be more problematic. The validity of cost prediction based on the theory of evolution of gene expression and resource allocation has been cast into doubt by contradictory experimental evidence. Further research providing insights into the role of the EPSPS cassette in weed adaptation, and estimations of the energy budget involved in EPSPS amplification and overexpression are required to understand and predict the biochemical and physiological bases of the fitness cost of glyphosate resistance.

Unravelling mesosulfuron-methyl phytotoxicity and metabolism-based herbicide resistance in Alopecurus aequalis: Insight into regulatory mechanisms using proteomics

Non-target-site based resistance (NTSR), a poorly understood multigenic trait, has evolved as the greatest threat to crop production worldwide, by endowing weed plants an unpredictable pattern of resistance to herbicides. Our recent work with multiple-herbicide-resistant shortawn foxtail (Alopecurus aequalis Sobol.) biotype has preliminary indicated that cytochrome P450s-involved enhanced rate of mesosulfuron-methyl metabolism may involve in the NTSR. Here by further determining the differences in glutathione S-transferase (GST) activity and uptake and metabolic rates of mesosulfuron between resistant (R) and susceptible (S) A. aequalis plants, and associating them with endogenous differently regulated proteins (DEPs) identified from combinational proteomics analyses, we provided direct evidences on the enhanced herbicide degradation in resistant plants. Subsequently, the physiological phenotypes of photosynthesis, chlorophyll fluorescence, and antioxidation were compared between R and S plants and linked with correlative DEPs, indicating a series of key pathways including solar energy capture, photosynthetic electron transport, redox homeostasis, carbon fixation, photorespiration, and reactive oxygen species scavenging in susceptible plants were broken or severely damaged by mesosulfuron stress. In comparison, resistant plants have evolved enhanced herbicide degradation to minimize the accumulation of mesosulfuron and protect the photosynthesis and ascorbate-glutathione cycle against the adverse effects of chemical injury, giving A. aequalis plants a NTSR phenotype. Additionally, three key proteins respectively annotated as esterase, GST, and glucosyltransferase were identified and enabled as potential transcriptional markers for quick diagnosing the metabolic mesosulfuron resistance in A. aequalis species.

Transcriptomic Analysis Identifies New Non-Target Site Glyphosate-Resistance Genes in Conyza bonariensis

Conyza bonariensis (hairy fleabane) is one of the most problematic and widespread glyphosate-resistant weeds in the world. This highly competitive weed species significantly interferes with crop growth and substantially decreases crop yield. Despite its agricultural importance, the molecular mechanisms of glyphosate resistance are still unknown. The present RNA-Seq study was performed with the goal of identifying differentially expressed candidate transcripts (genes) related to metabolism-based non-target site glyphosate resistance in C. bonariensis. The whole-transcriptome was de novo assembled from glyphosate-resistant and -sensitive biotypes of C. bonariensis from Southern Brazil. The RNA was extracted from untreated and glyphosate-treated plants at several timepoints up to 288 h after treatment in both biotypes. The transcriptome assembly produced 90,124 contigs with an average length of 777 bp and N50 of 1118 bp. In response to glyphosate treatment, differential gene expression analysis was performed on glyphosate-resistant and -sensitive biotypes. A total of 9622 genes were differentially expressed as a response to glyphosate treatment in both biotypes, 4297 (44.6%) being up- and 5325 (55.4%) down-regulated. The resistant biotype presented 1770 up- and 2333 down-regulated genes while the sensitive biotype had 2335 and 2800 up- and down-regulated genes, respectively. Among them, 974 up- and 1290 down-regulated genes were co-expressed in both biotypes. In the present work, we identified 41 new candidate target genes from five families related to herbicide transport and metabolism: 19 ABC transporters, 10 CYP450s, one glutathione S-transferase (GST), five glycosyltransferases (GT), and six genes related to antioxidant enzyme catalase (CAT), peroxidase (POD), and superoxide dismutase (SOD). The candidate genes may participate in metabolic-based glyphosate resistance via oxidation, conjugation, transport, and degradation, plus antioxidation. One or more of these genes might 'rescue' resistant plants from irreversible damage after glyphosate treatment. The 41 target genes we report in the present study may inform further functional genomics studies, including gene editing approaches to elucidate glyphosate-resistance mechanisms in C. bonariensis.

Defenses Against ROS in Crops and Weeds: The Effects of Interference and Herbicides

The antioxidant defense system acts to maintain the equilibrium between the production of reactive oxygen species (ROS) and the elimination of toxic levels of ROS in plants. Overproduction and accumulation of ROS results in metabolic disorders and can lead to the oxidative destruction of the cell. Several stress factors cause ROS overproduction and trigger oxidative stress in crops and weeds. Recently, the involvement of the antioxidant system in weed interference and herbicide treatment in crops and weeds has been the subject of investigation. In this review, we address ROS production and plant mechanisms of defense, alterations in the antioxidant system at transcriptional and enzymatic levels in crops induced by weed interference, and herbicide exposure in crops and weeds. We also describe the mechanisms of action in herbicides that lead to ROS generation in target plants. Lastly, we discuss the relations between antioxidant systems and weed biology and evolution, as well as the interactive effects of herbicide treatment on these factors.

Target-site and non-target-site-based resistance to tribenuron-methyl in multiply-resistant Myosoton aquaticum L

Myosoton aquaticum L., a widespread and competitive winter weed of wheat in China, has evolved resistance to many classes of herbicides. In one M. aquaticum population (AH03), collected from Anhui Province, where tribenuron-methyl and florasulam had been used to control this weed resistance to both herbicides had evolved. Compared with the sensitive population, HN03(S), the resistant (R) population, AH03, was highly resistant to tribenuron-methyl, flucarbazone-Na and pyroxsulam, moderately resistant to pyrithiobac‑sodium, and florasulam, and had low resistance to diflufenican. AH03 was still controlled by imazethapyr, 2,4-D butylate, fluroxypyr-meptyl, and isoproturon. Pretreatment with the P450 inhibitor malathion reduced the GR₅₀ value of tribenuron-methyl by 43% in the R population, and by 25% in the S population. This indicates that P450-mediated enhanced metabolism is one likely mechanism for tribenuron-methyl resistance in M. aquaticum. Glutathione-S-transferase (GST) activity could be induced by tribenuron-methyl in both the R and S populations. However, both the basal and induced GST activity of the R population was lower than that of the S population. The in vitro ALS assay confirmed that the ALS from the R plants showed a high resistance (52.93-fold) to tribenuron-methyl. ALS gene sequencing revealed a Pro197Ala substitution in the R plants. Based on the ALS gene sequence analysis, molecular markers were also developed to identify the specific Pro197Ala mutation. This population of M. aquaticum has multiple resistance and target-site (ALS Pro197Ala) and non-target-site resistance mechanisms contribute to tribenuron-methyl resistance.

The evolutionary origins of pesticide resistance

Durable crop protection is an essential component of current and future food security. However, the effectiveness of pesticides is threatened by the evolution of resistant pathogens, weeds and insect pests. Pesticides are mostly novel synthetic compounds, and yet target species are often able to evolve resistance soon after a new compound is introduced. Therefore, pesticide resistance provides an interesting case of rapid evolution under strong selective pressures, which can be used to address fundamental questions concerning the evolutionary origins of adaptations to novel conditions. We ask: (i) whether this adaptive potential originates mainly from de novo mutations or from standing variation; (ii) which pre-existing traits could form the basis of resistance adaptations; and (iii) whether recurrence of resistance mechanisms among species results from interbreeding and horizontal gene transfer or from independent parallel evolution. We compare and contrast the three major pesticide groups: insecticides, herbicides and fungicides. Whilst resistance to these three agrochemical classes is to some extent united by the common evolutionary forces at play, there are also important differences. Fungicide resistance appears to evolve, in most cases, by de novo point mutations in the target-site encoding genes; herbicide resistance often evolves through selection of polygenic metabolic resistance from standing variation; and insecticide resistance evolves through a combination of standing variation and de novo mutations in the target site or major metabolic resistance genes. This has practical implications for resistance risk assessment and management, and lessons learnt from pesticide resistance should be applied in the deployment of novel, non-chemical pest-control methods.

HRGPred: Prediction of herbicide resistant genes with k-mer nucleotide compositional features and support vector machine

Herbicide resistance (HR) is a major concern for the agricultural producers as well as environmentalists. Resistance to commonly used herbicides are conferred due to mutation(s) in the genes encoding herbicide target sites/proteins (GETS). Identification of these genes through wet-lab experiments is time consuming and expensive. Thus, a supervised learning-based computational model has been proposed in this study, which is first of its kind for the prediction of seven classes of GETS. The cDNA sequences of the genes were initially transformed into numeric features based on the k-mer compositions and then supplied as input to the support vector machine. In the proposed SVM-based model, the prediction occurs in two stages, where a binary classifier in the first stage discriminates the genes involved in conferring the resistance to herbicides from other genes, followed by a multi-class classifier in the second stage that categorizes the predicted herbicide resistant genes in the first stage into any one of the seven resistant classes. Overall classification accuracies were observed to be ~89% and >97% for binary and multi-class classifications respectively. The proposed model confirmed higher accuracy than the homology-based algorithms viz., BLAST and Hidden Markov Model. Besides, the developed computational model achieved ~87% accuracy, while tested with an independent dataset. An online prediction server HRGPred (http://cabgrid.res.in:8080/hrgpred) has also been established to facilitate the prediction of GETS by the scientific community.

Target site mutations and cytochrome P450s confer resistance to fenoxaprop-P-ethyl and mesosulfuron-methyl in Alopecurus aequalis

BACKGROUND

Shortawn foxtail (Alopecurus aequalis Sobol.) is a competitive grass weed infesting winter wheat- and canola-growing fields in China. In May 2016, a suspected A. aequalis resistant population AHTC-06 that survived fenoxaprop-P-ethyl and mesosulfuron-methyl applied at their field-recommended rates was collected from a wheat field in Jinji County, Anhui Province, China. This study aimed to determine the resistance profile of this AHTC-06 population to ACCase- and ALS-inhibitors, and to investigate its mechanisms of resistance to fenoxaprop-P-ethyl and mesosulfuron-methyl.

RESULTS

Two mutations, a common Ile-2041-Asn (ACCase gene) and a very rare Pro-197-Tyr (ALS1 gene), were both identified in resistant individual plants. The homozygous subpopulation AHTC-06F1 for the two mutations was generated, and it showed broad-spectrum resistance to APPs, DENs, and ALS-inhibiting herbicides of all five chemical families tested, with resistance index (RI) values that ranged from 2.2 to 36.5. In vitro ALS activity assays showed the ALS from the resistant population was insensitive to all the tested ALS inhibitors, with RI values ranging from 3.10 to 22.51. Pre-treatment with piperonyl butoxide (PBO) and malathion significantly (P < 0.05) reversed the weed's resistance to fenoxaprop-P-ethyl and mesosulfuron-methyl, respectively. Two P450 genes, c21190_g1 and c43350_g3, were constitutively overexpressed and mesosulfuron-methyl-induced upregulated in resistant plants, for which c43350_g3 was also fenoxaprop-P-ethyl-induced upregulated.

CONCLUSION

This study confirms the first case of a grass weed featuring broad-spectrum resistance to ALS-inhibiting herbicides due to a Pro-197-Tyr mutation in the ALS gene. Fenoxaprop-P-ethyl and mesosulfuron-methyl resistances in AHTC-06 plants were conferred by target site mutations and P450s-based metabolism. © 2018 Society of Chemical Industry.

Genomic Clues for Crop-Weed Interactions and Evolution

Agronomically critical weeds that have evolved alongside crop species are characterized by rapid adaptation and invasiveness, which can result in an enormous reduction in annual crop yield worldwide. We discuss here recent genome-based research studies on agricultural weeds and crop-weed interactions that reveal several major evolutionary innovations such as de-domestication, interactions mediated by allelochemical secondary metabolites, and parasitic genetic elements that play crucial roles in enhancing weed invasiveness in agricultural settings. We believe that these key studies will guide future research into the evolution of crop-weed interactions, and further the development of practical applications in agricultural weed control and crop breeding.

Interspecific and intraspecific transference of metabolism-based mesotrione resistance in dioecious weedy Amaranthus

Pollen-mediated gene flow (PMGF) might play an important role in dispersing herbicide resistance alleles in dioecious weedy Amaranthus species. Field experiments in a concentric donor-receptor design were conducted to quantify two sets of PMGF studies, an interspecific (Amaranthus tuberculatus × Amaranthus palmeri) and an intraspecific (A. tuberculatus × A. tuberculatus). In both studies, PMGF was evaluated using a resistant A. tuberculatus phenotype with enhanced mesotrione detoxification via P450 enzymes as a source of resistance alleles. For interspecific hybridization, more than 104 000 putative hybrid seedlings were screened with three markers, one phenotypic and two molecular. The two molecular markers used, including 2-bp polymorphisms in the internal transcribed spacer region, distinguished A. palmeri, A. tuberculatus and their hybrids. Results showed that 0.1% hybridization between A. tuberculatus × A. palmeri occurred under field research conditions. For intraspecific hybridization, 22 582 seedlings were screened to assess the frequency of gene flow. The frequency of gene flow (F_GF ) varied with distance, direction and year of the study. The farthest distance for 90% reduction of F_GF was at 69 m in 2015 however, after averaging across directions it was 13.1 and 26.1 m in 2014 and 2015, respectively. This study highlights the transfer of metabolism-based mesotrione resistance from A. tuberculatus to A. palmeri under field research conditions. The results presented here might aid in the rapid detection of A. palmeri among other Amaranthus species and show that PMFG could be expediting the increase of herbicide resistance in A. palmeri and A. tuberculatus across US crop production areas.

Non-target site-based resistance to tribenuron-methyl and essential involved genes in Myosoton aquaticum (L.)

BACKGROUND

Water chickweed (Myosoton aquaticum (L.)) is a dicot broadleaf weed that is widespread in winter fields in China, and has evolved serious resistance to acetolactate synthase (ALS) inhibiting herbicides.

RESULTS

We identified a M. aquaticum population exhibiting moderate (6.15-fold) resistance to tribenuron-methyl (TM). Target-site ALS gene sequencing revealed no known resistance mutations in these plants, and the in vitro ALS activity assays showed no differences in enzyme sensitivity between susceptible and resistant populations; however, resistance was reversed by pretreatment with the cytochrome P450 (CYP) monooxygenase inhibitor malathion. An RNA sequencing transcriptome analysis was performed to identify candidate genes involved in metabolic resistance, and the unigenes obtained by de novo transcriptome assembly were annotated across seven databases. In total, 34 differentially expressed genes selected by digital gene expression analysis were validated by quantitative real-time (qRT)-PCR. Ten consistently overexpressed contigs, including four for CYP, four for ATP-binding cassette (ABC) transporter, and two for peroxidase were further validated by qRT-PCR using additional plants from resistant and susceptible populations. Three CYP genes (with homology to CYP734A1, CYP76C1, and CYP86B1) and one ABC transporter gene (with homology to ABCC10) were highly expressed in all resistant plants.

CONCLUSION

The mechanism of TM resistance in M. aquaticum is controlled by NTSR rather than TSR. Four genes, CYP734A1, CYP76C1, CYP86B1, and ABCC10 could play essential role in metabolic resistance to TM and justify further functional studies. To our knowledge, this is the first large-scale transcriptome analysis of genes associated with NTSR in M. aquaticum using the Illumina platform. Our data provide resource for M. aquaticum biology, and will facilitate the study of herbicide resistance mechanism at the molecular level in this species as well as in other weeds.

Optimizing RNA-seq studies to investigate herbicide resistance

Transcriptomic profiling, specifically via RNA sequencing (RNA-seq), is becoming one of the more commonly used methods for investigating non-target site resistance (NTSR) to herbicides due to its high throughput capabilities and utility in organisms with little to no previous sequence information. A review of the weed science RNA-seq literature revealed some basic principles behind generating quality data from these types of studies. First, studies that included more replicates per biotype and took steps to control for genetic background had significantly better control of false positives and, consequently, shorter lists of potential resistance genes to sift through. Pooling of biological replicates prior to sequencing was successful in some cases, but likely contributed to an overall increase in the false discovery rate. Although the inclusion of herbicide-treated samples was common across most of the studies, it ultimately introduced difficulties in interpretation of the final results due to challenges in capturing the right sampling window after treatment and to the induction of stress responses in the injured herbicide-sensitive plants. RNA-seq is an effective tool for NTSR gene discovery, but careful consideration should be given to finding the most powerful and cost-effective balance between replicate number, sequencing depth and treatment number. © 2017 Society of Chemical Industry.

Tembotrione detoxification in 4-hydroxyphenylpyruvate dioxygenase (HPPD) inhibitor-resistant Palmer amaranth (Amaranthus palmeri S. Wats.)

BACKGROUND

Resistance to the 4-hydroxyphenylpyruvate dioxygenase (HPPD)-inhibiting herbicide tembotrione in an Amaranthus palmeri population from Nebraska (NER) has previously been confirmed to be attributable to enhanced metabolism. The objective of this study was to identify and quantify the metabolites formed in Nebraska susceptible (NES) and resistant (NER) biotypes.

RESULTS

NER and NES formed the same metabolites. Tembotrione metabolism in NER differed from that in NES in that resistant plants showed faster 4-hydroxylation followed by glycosylation. The T₅₀ value (time for 50% production of the maximum 4-hydroxylation product) was 4.9 and 11.9 h for NER and NES, respectively. This process is typically catalyzed by cytochrome P450 enzymes. Metabolism differences between NER and NES were most prominent under 28 °C conditions and herbicide application at the four-leaf stage.

CONCLUSION

Further research with the aim of identifying the gene or genes responsible for conferring metabolic resistance to HPPD inhibitors should focus on cytochrome P450s. Such research is important because non-target-site-based resistance (NTSR) poses the threat of cross resistance to other chemical classes of HPPD inhibitors, other herbicide modes of action, or even unknown herbicides. © 2017 Society of Chemical Industry.

Metabolism of 2,4-dichlorophenoxyacetic acid contributes to resistance in a common waterhemp (Amaranthus tuberculatus) population

BACKGROUND

Synthetic auxins such as 2,4-dichlorophenoxyacetic acid (2,4-D) have been widely used for selective control of broadleaf weeds since the mid-1940s. In 2009, an Amaranthus tuberculatus (common waterhemp) population with 10-fold resistance to 2,4-D was found in Nebraska, USA. The 2,4-D resistance mechanism was examined by conducting [¹⁴ C] 2,4-D absorption, translocation and metabolism experiments.

RESULTS

No differences were found in 2,4-D absorption or translocation between resistant and susceptible A. tuberculatus plants. Resistant plants metabolized [¹⁴ C] 2,4-D more rapidly than did susceptible plants. The half-life of [¹⁴ C] 2,4-D in susceptible plants was 105 h, compared with 22 h in resistant plants. Pretreatment with the cytochrome P450 inhibitor malathion inhibited [¹⁴ C] 2,4-D metabolism in resistant plants and reduced the 2,4-D dose required for 50% growth inhibition (GR₅₀ ) of resistant plants by 7-fold to 27 g ha^-1 , similar to the GR₅₀ for susceptible plants in the absence of malathion.

CONCLUSION

Our results demonstrate that rapid 2,4-D metabolism is a contributing factor to resistance in A. tuberculatus, potentially mediated by cytochrome P450. Metabolism-based resistance to 2,4-D could pose a serious challenge for A. tuberculatus control because of the potential for cross-resistance to other herbicides. © 2017 Society of Chemical Industry.

Enhanced Herbicide Metabolism and Metabolic Resistance Genes Identified in Tribenuron-Methyl Resistant Myosoton aquaticum L

The evolved resistance of Myosoton aquaticum L. to acetolactate synthase (ALS) inhibitors is well established, but most research has focused on target-site resistance, while nontarget-site resistance remains neglected. Here, we investigated mechanisms of the latter. The pretreatment with the P450 inhibitor malathion significantly increased the sensitivity of resistant plants to tribenuron-methyl. The rapid P450-mediated tribenuron-methyl metabolism in resistant plants was confirmed by LC-MS/MS analysis. Besides, GST activity was higher among resistant than susceptible individuals. The next transcriptome analysis generated 544,102,236 clean reads from RNA sequencing libraries. De novo assembly yielded 102,529 unigenes with an average length of 866 bp, annotated across seven databases. Digital gene expression selected 25 differentially expressed genes, further validated with qRT-PCR. Three P450 genes, two GST genes, two glucosyltransferase genes, four ABC transporter genes, and four additional contigs were constitutively up-regulated in resistant individuals. Overall, our research confirmed that enhanced herbicide metabolism drives tribenuron-methyl resistance in M. aquaticum.

Stress-induced evolution of herbicide resistance and related pleiotropic effects

Herbicide-resistant weeds, especially those with resistance to multiple herbicides, represent a growing worldwide threat to agriculture and food security. Natural selection for resistant genotypes may act on standing genetic variation, or on a genetic and physiological background that is fundamentally altered because of stress responses to sublethal herbicide exposure. Stress-induced changes include DNA mutations, epigenetic alterations, transcriptional remodeling, and protein modifications, all of which can lead to herbicide resistance and a wide range of pleiotropic effects. Resistance selected in this manner is termed systemic acquired herbicide resistance, and the associated pleiotropic effects are manifested as a suite of constitutive transcriptional and post-translational changes related to biotic and abiotic stress adaptation, representing the evolutionary signature of selection. This phenotype is being investigated in two multiple herbicide-resistant populations of the hexaploid, self-pollinating weedy monocot Avena fatua that display such changes as well as constitutive reductions in certain heat shock proteins and their transcripts, which are well known as global regulators of diverse stress adaptation pathways. Herbicide-resistant populations of most weedy plant species exhibit pleiotropic effects, and their association with resistance genes presents a fertile area of investigation. This review proposes that more detailed studies of resistant A. fatua and other species through the lens of plant evolution under stress will inform improved resistant weed prevention and management strategies. © 2018 Society of Chemical Industry.

Pyroxasulfone resistance in Lolium rigidum is metabolism-based

The evolution of resistant weed populations in response to intensive herbicide selection pressure is a global issue. Resistance to post-emergence herbicides is widespread, whereas soil-applied pre-emergence herbicides can often remain effective. For example, in Australia pyroxasulfone is a new pre-emergence soil-applied herbicide which provides control of Lolium rigidum populations resistant to multiple post-emergence herbicide modes of action. A fundamental knowledge of the genetic basis of metabolic resistance in weeds is important for understanding plant evolution pathways under herbicide selection and sustaining long-term weed resistance management. In this study we define the mechanistic basis of resistance to pyroxasulfone in a L. rigidum population. TLC provides evidence that pyroxasulfone resistance is metabolism-based with approximately 88% of parental [¹⁴C]-labelled pyroxasulfone metabolized in resistant plants 24 h after the herbicide treatment. HPLC-MS allowed identification of several metabolites of pyroxasulfone formed via a glutathione (GSH) conjugation pathway in pyroxasulfone-resistant L. rigidum plants. However, the initial pyroxasulfone-glutathione conjugate was not found likely due to its labile nature. The observed constitutive over-expression from six to nine-fold of two putative resistance-endowing GST genes was associated with the pyroxasulfone resistance phenotype. The most logical conclusion, based on the data thus far available, is that rapid detoxification of pyroxasulfone mediates pyroxasulfone resistance in L. rigidum plants. Future research is warranted to confirm the hypothesis advanced by this study of rapid pyroxasulfone metabolism due to GSH conjugation mediated by GST over-expressed in pyroxasulfone-resistant plants which similarly leads to the production of distinctive GSH-pyroxasulfone metabolites in L. rigidum and wheat.

Metabolic Resistance to Acetolactate Synthase Inhibiting Herbicide Tribenuron-Methyl in Descurainia sophia L. Mediated by Cytochrome P450 Enzymes

Descurainia sophia is one of the most notorious broadleaf weeds in China and has evolved extremely high resistance to acetolactate synthase (ALS)-inhibiting herbicide tribenuron-methyl. The target-site resistance due to ALS gene mutations was well-known, while the non-target-site resistance is not yet well-characterized. Metabolic resistance, which is conferred by enhanced rates of herbicide metabolism, is the most important NTSR. To explore the mechanism of metabolic resistance underlying resistant (R) D. sophia plants, tribenuron-methyl uptake and metabolism levels, qPCR reference gene stability, and candidate P450 genes expression patterns were investigated. The results of liquid chromatography-mass spectrometry (LC-MS) analysis indicated that the metabolic rates of tribenuron-methyl in R plants was significantly faster than in susceptible (S) plants, and this metabolism differences can be eliminated by P450 inhibitor malathion. The genes for 18S rRNA and TIP41-like were identified as the most suitable reference genes using programs of BestKeeper, NormFinder, and geNorm. The P450 gene CYP96A146 constitutively overexpressed in R plants compared to S plants; this overexpression in R plants can be suppressed by malathion. Taken together, a higher expression level of P450 genes, leading to higher tribenuron-methyl metabolism, appears to be responsible for metabolic resistance to tribenuron-methyl in R D. sophia plants.

Changes in the proteome of the problem weed blackgrass correlating with multiple-herbicide resistance

Herbicide resistance in grass weeds is now one of the greatest threats to sustainable cereal production in Northern Europe. Multiple-herbicide resistance (MHR), a poorly understood multigenic and quantitative trait, is particularly problematic as it provides tolerance to most classes of chemistries currently used for post-emergence weed control. Using a combination of transcriptomics and proteomics, the evolution of MHR in populations of the weed blackgrass (Alopecurus myosuroides) has been investigated. While over 4500 genes showed perturbation in their expression in MHR versus herbicide sensitive (HS) plants, only a small group of proteins showed >2-fold changes in abundance, with a mere eight proteins consistently associated with this class of resistance. Of the eight, orthologues of three of these proteins are also known to be associated with multiple drug resistance (MDR) in humans, suggesting a cross-phyla conservation in evolved tolerance to chemical agents. Proteomics revealed that MHR could be classified into three sub-types based on the association with resistance to herbicides with differing modes of action (MoA), being either global, specific to diverse chemistries acting on one MoA, or herbicide specific. Furthermore, the proteome of MHR plants were distinct from that of HS plants exposed to a range of biotic (insect feeding, plant-microbe interaction) and abiotic (N-limitation, osmotic, heat, herbicide safening) challenges commonly encountered in the field. It was concluded that MHR in blackgrass is a uniquely evolving trait(s), associated with changes in the proteome that are distinct from responses to conventional plant stresses, but sharing common features with MDR in humans.

Protein identification before and after glyphosate exposure in Lolium multiflorum genotypes

BACKGROUND

Weeds reduce crop yields, and among the methods used to control these plants, the use of chemicals is preferred. However, the repeated application of herbicides with the same mechanism of action selects for resistant populations. The aim of this study was to evaluate glyphosate resistance in Lolium multiflorum (Lam.) and relate the resistance to protein expression in the absence and presence of the herbicide using a metabolic-proteomic approach.

RESULTS

Glyphosate resistance was confirmed, with a sevenfold difference in resistance between susceptible and resistant genotypes. Among the possible mechanisms affecting resistance, mutations in the enzyme 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS), herbicide differential translocation and overexpression of EPSPS are suggested. Susceptible plants had higher growth than did resistant plants in the absence of the herbicide, in addition to greater expression of protein groups related to photosynthesis and to tolerance to biotic and abiotic stresses. With application of glyphosate, resistant plants maintained their metabolism and began to express EPSPS and other candidate proteins related to herbicide resistance.

CONCLUSIONS

In the absence of glyphosate, the susceptible plants would replace the resistant plants over time, and abiotic or biotic stresses would accelerate this process. Resistance in plants resulted from a combination of target-site and non-target-site resistance mechanisms. We identified several candidate proteins that could be investigated in future studies on glyphosate resistance. © 2017 Society of Chemical Industry.

RNA-Seq transcriptome analysis of Amaranthus palmeri with differential tolerance to glufosinate herbicide

Amaranthus palmeri (Amaranthaceae) is a noxious weed in several agroecosystems and in some cases seriously threatens the sustainability of crop production in North America. Glyphosate-resistant Amaranthus species are widespread, prompting the use of alternatives to glyphosate such as glufosinate, in conjunction with glufosinate-resistant crop cultivars, to help control glyphosate-resistant weeds. An experiment was conducted to analyze the transcriptome of A. palmeri plants that survived exposure to 0.55 kg ha-1 glufosinate. Since there was no record of glufosinate use at the collection site, survival of plants within the population are likely due to genetic expression that pre-dates selection; in the formal parlance of weed science this is described as natural tolerance. Leaf tissues from glufosinate-treated and non-treated seedlings were harvested 24 h after treatment (HAT) for RNA-Seq analysis. Global gene expression was measured using Illumina DNA sequence reads from non-treated and treated surviving (presumably tolerant, T) and susceptible (S) plants. The same plants were used to determine the mechanisms conferring differential tolerance to glufosinate. The S plants accumulated twice as much ammonia as did the T plants, 24 HAT. The relative copy number of the glufosinate target gene GS2 did not differ between T and S plants, with 1 to 3 GS2 copies in both biotypes. A reference cDNA transcriptome consisting of 72,780 contigs was assembled, with 65,282 sequences putatively annotated. Sequences of GS2 from the transcriptome assembly did not have polymorphisms unique to the tolerant plants. Five hundred sixty-seven genes were differentially expressed between treated T and S plants. Of the upregulated genes in treated T plants, 210 were more highly induced than were the upregulated genes in the treated S plants. Glufosinate-tolerant plants had greater induction of ABC transporter, glutathione S-transferase (GST), NAC transcription factor, nitronate monooxygenase (NMO), chitin elicitor receptor kinase (CERK1), heat shock protein 83, ethylene transcription factor, heat stress transcription factor, NADH-ubiquinone oxidoreductase, ABA 8'-hydroxylase, and cytochrome P450 genes (CYP72A, CYP94A1). Seven candidate genes were selected for validation using quantitative real time-PCR. While GST was upregulated in treated tolerant plants in at least one population, CYP72A219 was consistently highly expressed in all treated tolerant biotypes. These genes are candidates for contributing tolerance to glufosinate. Taken together, these results show that differential induction of stress-protection genes in a population can enable some individuals to survive herbicide application. Elevated expression of detoxification-related genes can get fixed in a population with sustained selection pressure, leading to evolution of resistance. Alternatively, sustained selection pressure could select for mutation(s) in the GS2 gene with the same consequence.

Glyphosate Resistance and EPSPS Gene Duplication: Convergent Evolution in Multiple Plant Species

One of the increasingly widespread mechanisms of resistance to the herbicide glyphosate is copy number variation (CNV) of the 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) gene. EPSPS gene duplication has been reported in 8 weed species, ranging from 3 to 5 extra copies to more than 150 extra copies. In the case of Palmer amaranth (Amaranthus palmeri), a section of >300 kb containing EPSPS and many other genes has been replicated and inserted at new loci throughout the genome, resulting in significant increase in total genome size. The replicated sequence contains several classes of mobile genetic elements including helitrons, raising the intriguing possibility of extra-chromosomal replication of the EPSPS-containing sequence. In kochia (Kochia scoparia), from 3 to more than 10 extra EPSPS copies are arranged as a tandem gene duplication at one locus. In the remaining 6 weed species that exhibit EPSPS gene duplication, little is known about the underlying mechanisms of gene duplication or their entire sequence. There is mounting evidence that adaptive gene amplification is an important mode of evolution in the face of intense human-mediated selection pressure. The convergent evolution of CNVs for glyphosate resistance in weeds, through at least 2 different mechanisms, may be indicative of a more general importance for this mechanism of adaptation in plants. CNVs warrant further investigation across plant functional genomics for adaptation to biotic and abiotic stresses, particularly for adaptive evolution on rapid time scales.

Transcriptomic changes in Echinochloa colona in response to treatment with the herbicide imazamox

Presented here is the first Echinochloa colona leaf transcriptome. Analysis of gene expression before and after herbicide treatment reveals that E. colona mounts a stress response upon exposure to herbicide. Herbicides are the most frequently used means of controlling weeds. For many herbicides, the target site is known; however, it is considerably less clear how plant gene expression changes in response to herbicide exposure. In this study, changes in gene expression in response to herbicide exposure in imazamox-sensitive (S) and- resistant (R) junglerice (Echinochloa colona L.) biotypes was examined. As no reference genome is available for this weed, a reference leaf transcriptome was generated. Messenger RNA was isolated from imazamox-treated- and untreated R and S plants and the resulting cDNA libraries were sequenced on an Illumina HiSeq2000. The transcriptome was assembled, annotated, and differential gene expression analysis was performed to identify transcripts that were upregulated or downregulated in response to herbicide exposure for both biotypes. Differentially expressed transcripts included transcription factors, protein-modifying enzymes, and enzymes involved in metabolism and signaling. A literature search revealed that members of the families represented in this analysis were known to be involved in abiotic stress response in other plants, suggesting that imazamox exposure induced a stress response. A time course study examining a subset of transcripts showed that expression peaked within 4-12 h and then returned to untreated levels within 48 h of exposure. Testing of plants from two additional biotypes showed a similar change in gene expression 4 h after herbicide exposure compared to the resistant and sensitive biotypes. This study shows that within 48 h junglerice mounts a stress response to imazamox exposure.

Constitutive redox and phosphoproteome changes in multiple herbicide resistant Avena fatua L. are similar to those of systemic acquired resistance and systemic acquired acclimation

Plants are routinely confronted with numerous biotic and abiotic stressors, and in response have evolved highly effective strategies of systemic acquired resistance (SAR) and systemic acquired acclimation (SAA), respectively. A much more evolutionarily recent abiotic stress is the application of herbicides to control weedy plants, and their intensive use has selected for resistant weed populations that cause substantial crop yield losses and increase production costs. Non-target site resistance (NTSR) to herbicides is rapidly increasing worldwide and is associated with alterations in generalized stress defense networks. This work investigated protein post-translational modifications associated with NTSR in multiple herbicide resistant (MHR) Avena fatua, and their commonalities with those of SAR and SAA. We used proteomic, biochemical, and immunological approaches to compare constitutive protein profiles in MHR and herbicide susceptible (HS) A. fatua populations. Phosphoproteome and redox proteome surveys showed that post-translational modifications of proteins with functions in core cellular processes were reduced in MHR plants, while those involved in xenobiotic and stress response, reactive oxygen species detoxification and redox maintenance, heat shock response, and intracellular signaling were elevated in MHR as compared to HS plants. More specifically, MHR plants contained constitutively elevated levels of three protein kinases including the lectin S-receptor-like serine/threonine-protein kinase LecRK2, a well-characterized component of SAR. Analyses of superoxide dismutase enzyme activity and protein levels did not reveal constitutive differences between MHR and HS plants. The overall results support the idea that herbicide stress is perceived similarly to other abiotic stresses, and that A. fatua NTSR shares analogous features with SAR and SAA. We speculate that MHR A. fatua's previous exposure to sublethal herbicide doses, as well as earlier evolution under a diversity of abiotic and biotic stressors, has led to a heightened state of stress preparedness that includes NTSR to a number of unrelated herbicides.

Multiple mechanisms are involved in new imazamox-resistant varieties of durum and soft wheat

Weed control in wheat is one of the major goals of farmers in their efforts toward obtaining the highest crop yields for human foods. Several studies (dose-response, enzyme activity, absorption-translocation and metabolism) were conducted to characterize the resistance level of two new wheat cultivars called Rafalín (Triticum aestivum) and Antoñín (T. durum) that were obtained by conventional breeding based on Clearfield® technology; they are resistant (R) to imazamox compared to their sensitive (S) counterparts (Gazul and Simeto, respectively). The R-cultivars were 93.7-fold (Rafalín) and 43.7-fold (Antoñín) more resistant than their respective S-cultivars. The acetolactate synthase (ALS) enzyme activity revealed high resistance to imidazolinone (IMI) herbicides in R-cultivars, but no cross-resistance to other ALS herbicides was found. The Ser653Asn mutation that confers resistance to IMI herbicides was identified in the imi1 and imi2 genes of Rafalín and only in the imi1 gene of Antoñín. The ¹⁴C-imazamox absorption did not differ between the R- and S-cultivars. Imazamox was metabolized by Cyt-P450 into imazamox-hydroxyl and imazamox-glucoside in the R-cultivars, altering their translocation patterns. The differential sensitivity to imazamox between R-cultivars was due to the number of resistance genes that carry each genotype. The R-cultivars Rafalín and Antoñín could be excellent weed control tools.

Intensive herbicide use has selected for constitutively elevated levels of stress-responsive mRNAs and proteins in multiple herbicide-resistant Avena fatua L

BACKGROUND

Intensive use of herbicides has led to the evolution of two multiple herbicide-resistant (MHR) Avena fatua (wild oat) populations in Montana that are resistant to members of all selective herbicide families available for A. fatua control in US small grain crops. We used transcriptome and proteome surveys to compare constitutive changes in MHR and herbicide-susceptible (HS) plants associated with non-target site resistance.

RESULTS

Compared to HS plants, MHR plants contained constitutively elevated levels of differentially expressed genes (DEGs) with functions in xenobiotic catabolism, stress response, redox maintenance and transcriptional regulation that are similar to abiotic stress-tolerant phenotypes. Proteome comparisons identified similarly elevated proteins including biosynthetic and multifunctional enzymes in MHR plants. Of 25 DEGs validated by RT-qPCR assay, differential regulation of 21 co-segregated with flucarbazone-sodium herbicide resistance in F₃ families, and a subset of 10 of these were induced or repressed in herbicide-treated HS plants.

CONCLUSION

Although the individual and collective contributions of these DEGs and proteins to MHR remain to be determined, our results support the idea that intensive herbicide use has selected for MHR populations with altered, constitutively regulated patterns of gene expression that are similar to those in abiotic stress-tolerant plants. © 2017 Society of Chemical Industry.

Biochemical characterization of metabolism-based atrazine resistance in Amaranthus tuberculatus and identification of an expressed GST associated with resistance

Rapid detoxification of atrazine in naturally tolerant crops such as maize (Zea mays) and grain sorghum (Sorghum bicolor) results from glutathione S-transferase (GST) activity. In previous research, two atrazine-resistant waterhemp (Amaranthus tuberculatus) populations from Illinois, U.S.A. (designated ACR and MCR), displayed rapid formation of atrazine-glutathione (GSH) conjugates, implicating elevated rates of metabolism as the resistance mechanism. Our main objective was to utilize protein purification combined with qualitative proteomics to investigate the hypothesis that enhanced atrazine detoxification, catalysed by distinct GSTs, confers resistance in ACR and MCR. Additionally, candidate AtuGST expression was analysed in an F₂ population segregating for atrazine resistance. ACR and MCR showed higher specific activities towards atrazine in partially purified ammonium sulphate and GSH affinity-purified fractions compared to an atrazine-sensitive population (WCS). One-dimensional electrophoresis of these fractions displayed an approximate 26-kDa band, typical of GST subunits. Several phi- and tau-class GSTs were identified by LC-MS/MS from each population, based on peptide similarity with GSTs from Arabidopsis. Elevated constitutive expression of one phi-class GST, named AtuGSTF2, correlated strongly with atrazine resistance in ACR and MCR and segregating F₂ population. These results indicate that AtuGSTF2 may be linked to a metabolic mechanism that confers atrazine resistance in ACR and MCR.

Perspectives on non-target site mechanisms of herbicide resistance in weedy plant species using evolutionary physiology

Evolutionary physiology merges the disciplines of evolution and physiology, and it is a research approach that has not received much attention for studying the development of herbicide resistance. This paper makes a case for using evolutionary physiology more frequently when studying herbicide resistance, and illustrates this using three areas where more work would be useful: (i) the interaction among major and minor alleles over many generations during the evolution of physiological responses that lead to specific mechanisms of resistance; (ii) the role of epigenetic factors, especially at an early stage of evolution, on the physiological modifications that result in phenotypes that become insensitive to herbicides; and (iii) the interaction between fitness and physiological performance over time, with emphasis on understanding mechanisms that improve the fitness of herbicide-resistant phenotypes during selection.

Transcription of putative tonoplast transporters in response to glyphosate and paraquat stress in Conyza bonariensis and Conyza canadensis and selection of reference genes for qRT-PCR

Herbicide resistance is a challenge for modern agriculture further complicated by cases of resistance to multiple herbicides. Conyza bonariensis and Conyza canadensis are invasive weeds of field crops, orchards, and non-cropped areas in many parts of the world. In California, USA, Conyza populations resistant to the herbicides glyphosate and paraquat have recently been described. Although the mechanism conferring resistance to glyphosate and paraquat in these species was not elucidated, reduced translocation of these herbicides was observed under experimental conditions in both species. Glyphosate and paraquat resistance associated with reduced translocation are hypothesized to be a result of sequestration of herbicides into the vacuole, with the possible involvement of over-expression of genes encoding tonoplast transporters of ABC-transporter families in cases of glyphosate resistance or cationic amino acid transporters (CAT) in cases of paraquat resistance. However, gene expression in response to herbicide treatment has not been studied in glyphosate and paraquat resistant populations. In the current study, we evaluated the transcript levels of genes possibly involved in resistance using real-time PCR. First, we evaluated eight candidate reference genes following herbicide treatment and selected three genes that exhibited stable expression profiles; ACTIN, HEAT-SHOCK-PROTEIN-70, and CYCLOPHILIN. The reference genes identified here can be used for further studies related to plant-herbicide interactions. We used these reference genes to assay the transcript levels of EPSPS, ABC transporters, and CAT in response to herbicide treatment in susceptible and resistant Conyza spp. lines. No transcription changes were observed in EPSPS or CAT genes after glyphosate or paraquat treatment, suggesting that these genes are not involved in the resistance mechanism. Transcription of the two ABC transporter genes increased following glyphosate treatment in all Conyza spp. lines. Transcription of ABC transporters also increased after paraquat treatment in all three lines of C. bonariensis. However, in C. canadensis, paraquat treatment increased transcription of only one ABC transporter gene in the susceptible line. The increase in transcription of ABC transporters after herbicide treatment is likely a stress response based on similar response observed across all Conyza lines regardless of resistance or sensitivity to glyphosate or paraquat, thus these genes do not appear to be directly involved in the mechanism of resistance in Conyza spp.

Genetic diversity and structure of Lolium perenne ssp. multiflorum in California vineyards and orchards indicate potential for spread of herbicide resistance via gene flow

Management of agroecosystems with herbicides imposes strong selection pressures on weedy plants leading to the evolution of resistance against those herbicides. Resistance to glyphosate in populations of Lolium perenne L. ssp. multiflorum is increasingly common in California, USA, causing economic losses and the loss of effective management tools. To gain insights into the recent evolution of glyphosate resistance in L. perenne in perennial cropping systems of northwest California and to inform management, we investigated the frequency of glyphosate resistance and the genetic diversity and structure of 14 populations. The sampled populations contained frequencies of resistant plants ranging from 10% to 89%. Analyses of neutral genetic variation using microsatellite markers indicated very high genetic diversity within all populations regardless of resistance frequency. Genetic variation was distributed predominantly among individuals within populations rather than among populations or sampled counties, as would be expected for a wide-ranging outcrossing weed species. Bayesian clustering analysis provided evidence of population structuring with extensive admixture between two genetic clusters or gene pools. High genetic diversity and admixture, and low differentiation between populations, strongly suggest the potential for spread of resistance through gene flow and the need for management that limits seed and pollen dispersal in L. perenne.

Unraveling the Transcriptional Basis of Temperature-Dependent Pinoxaden Resistance in Brachypodium hybridum

Climate change endangers food security and our ability to feed the ever-increasing human population. Weeds are the most important biotic stress, reducing crop-plant productivity worldwide. Chemical control, the main approach for weed management, can be strongly affected by temperature. Previously, we have shown that temperature-dependent non-target site (NTS) resistance of Brachypodium hybridum is due to enhanced detoxification of acetyl-CoA carboxylase inhibitors. Here, we explored the transcriptional basis of this phenomenon. Plants were characterized for the transcriptional response to herbicide application, high-temperature and their combination, in an attempt to uncover the genetic basis of temperature-dependent pinoxaden resistance. Even though most of the variance among treatments was due to pinoxaden application (61%), plants were able to survive pinoxaden application only when grown under high-temperatures. Biological pathways and expression patterns of members of specific gene families, previously shown to be involved in NTS metabolic resistance to different herbicides, were examined. Cytochrome P450, glucosyl transferase and glutathione-S-transferase genes were found to be up-regulated in response to pinoxaden application under both control and high-temperature conditions. However, biological pathways related to oxidation and glucose conjugation were found to be significantly enriched only under the combination of pinoxaden application and high-temperature. Analysis of reactive oxygen species (ROS) was conducted at several time points after treatment using a probe detecting H2O2/peroxides. Comparison of ROS accumulation among treatments revealed a significant reduction in ROS quantities 24 h after pinoxaden application only under high-temperature conditions. These results may indicate significant activity of enzymatic ROS scavengers that can be correlated with the activation of herbicide-resistance mechanisms. This study shows that up-regulation of genes related to metabolic resistance is not sufficient to explain temperature-dependent pinoxaden resistance. We suggest that elevated activity of enzymatic processes at high-temperature may induce rapid and efficient pinoxaden metabolism leading to temperature-dependent herbicide resistance.

Selection of relatively exact reference genes for gene expression studies in goosegrass (Eleusine indica) under herbicide stress

Goosegrass (Eleusine indica) is one of the most serious annual grassy weeds worldwide, and its evolved herbicide-resistant populations are more difficult to control. Quantitative real-time PCR (qPCR) is a common technique for investigating the resistance mechanism; however, there is as yet no report on the systematic selection of stable reference genes for goosegrass. This study proposed to test the expression stability of 9 candidate reference genes in goosegrass in different tissues and developmental stages and under stress from three types of herbicide. The results show that for different developmental stages and organs (control), eukaryotic initiation factor 4 A (eIF-4) is the most stable reference gene. Chloroplast acetolactate synthase (ALS) is the most stable reference gene under glyphosate stress. Under glufosinate stress, eIF-4 is the best reference gene. Ubiquitin-conjugating enzyme (UCE) is the most stable reference gene under quizalofop-p-ethyl stress. The gene eIF-4 is the recommended reference gene for goosegrass under the stress of all three herbicides. Moreover, pairwise analysis showed that seven reference genes were sufficient to normalize the gene expression data under three herbicides treatment. This study provides a list of reliable reference genes for transcript normalization in goosegrass, which will facilitate resistance mechanism studies in this weed species.

Physiological and Molecular Characterization of Hydroxyphenylpyruvate Dioxygenase (HPPD)-inhibitor Resistance in Palmer Amaranth (Amaranthus palmeri S.Wats.)

Herbicides that inhibit hydroxyphenylpyruvate dioxygenase (HPPD) such as mesotrione are widely used to control a broad spectrum of weeds in agriculture. Amaranthus palmeri is an economically troublesome weed throughout the United States. The first case of evolution of resistance to HPPD-inhibiting herbicides in A. palmeri was documented in Kansas (KS) and later in Nebraska (NE). The objective of this study was to investigate the mechansim of HPPD-inhibitor (mesotrione) resistance in A. palmeri. Dose response analysis revealed that this population (KSR) was 10-18 times more resistant than their sensitive counterparts (MSS or KSS). Absorbtion and translocation analysis of [¹⁴C] mesotrione suggested that these mechanisms were not involved in the resistance in A. palmeri. Importantly, mesotrione (>90%) was detoxified markedly faster in the resistant populations (KSR and NER), within 24 hours after treatment (HAT) compared to sensitive plants (MSS, KSS, or NER). However, at 48 HAT all populations metabolized the mesotrione, suggesting additional factors may contribute to this resistance. Further evaluation of mesotrione-resistant A. palmeri did not reveal any specific resistance-conferring mutations nor amplification of HPPD gene, the molecular target of mesotrione. However, the resistant populations showed 4- to 12-fold increase in HPPD gene expression. This increase in HPPD transcript levels was accompanied by increased HPPD protein expression. The significant aspects of this research include: the mesotrione resistance in A. palmeri is conferred primarily by rapid detoxification (non-target-site based) of mesotrione; additionally, increased HPPD gene expression (target-site based) also contributes to the resistance mechanism in the evolution of herbicide resistance in this naturally occurring weed species.

Transcriptional markers enable identification of rye-grass (Lolium sp.) plants with non-target-site-based resistance to herbicides inhibiting acetolactate-synthase

Molecular detection of herbicide non-target-site-based resistance (NTSR) classically requires extensively validated NTSR genes. We assessed the feasibility of predicting NTSR phenotypes using expression data of NTSR transcriptional markers, i.e., transcripts which expression levels are statistically correlated to NTSR. Markers were sought by comparative RNA-Seq analysis of untreated NTSR or sensitive plants from four rye-grass populations followed by expression quantification in 299 individual plants with characterised sensitivity to two acetolactate-synthase-inhibiting herbicides. Multivariate analyses were implemented to predict NTSR using combined marker expression data. Nineteen markers (four cytochromes P450, four glutathione-S-transferases, three glycosyltransferases, two ABC transporters, two hydrolases, one aldolase, one peptidase, one transferase and one esterase) expressed significantly higher in NTSR plants were identified. Expression was highest in the most resistant plants. Some markers appeared co-regulated. Combined marker expression data enabled prediction of NTSR phenotypes in individual plants or of resistant plant frequencies in populations. Thus, NTSR detection based on transcriptional markers proved feasible. Accuracy can be improved by identifying additional markers, especially markers associated to NTSR regulation. Additionally, our data suggest that NTSR mechanisms emerged in different populations via redundant evolution, and that NTSR can evolve by selection for higher constitutive expression of whole herbicide-response pathways.

De novo transcriptome assembly analysis of weed Apera spica-venti from seven tissues and growth stages

BACKGROUND

Loose silky bentgrass (Apera spica-venti) is an important weed in Europe with a recent increase in herbicide resistance cases. The lack of genetic information about this noxious weed limits its biological understanding such as growth, reproduction, genetic variation, molecular ecology and metabolic herbicide resistance. This study produced a reference transcriptome for A. spica-venti from different tissues (leaf, root, stem) and various growth stages (seed at phenological stages 05, 07, 08, 09). The de novo assembly was performed on individual and combined dataset followed by functional annotations. Individual transcripts and gene families involved in metabolic based herbicide resistance were identified.

RESULTS

Eight separate transcriptome assemblies were performed and compared. The combined transcriptome assembly consists of 83,349 contigs with an N50 and average contig length of 762 and 658 bp, respectively. This dataset contains 74,724 transcripts consisting of total 54,846,111 bp. Among them 94% had a homologue to UniProtKB, 73% retrieved a GO mapping, and 50% were functionally annotated. Compared with other grass species, A. spica-venti has 26% proteins in common to Brachypodium distachyon, and 41% to Lolium spp. Glycosyltransferases had the highest number of transcripts in each tissue followed by the cytochrome P450s. The GSTF1 and CYP89A2 transcripts were recovered from the majority of tissues and aligned at a maximum of 66 and 30% to proven herbicide resistant allele from Alopecurus myosuroides and Lolium rigidum, respectively.

CONCLUSIONS

De novo transcriptome assembly enabled the generation of the first reference transcriptome of A. spica-venti. This can serve as stepping stone for understanding the metabolic herbicide resistance as well as the general biology of this problematic weed. Furthermore, this large-scale sequence data is a valuable scientific resource for comparative transcriptome analysis for Poaceae grasses.