Investigating the mechanisms of glyphosate resistance in goosegrass (Eleusine indica (L.) Gaertn.) by RNA sequencing technology

Molecular responses reveal that two glutathione S-transferase CsGSTU8s contribute to detoxification of glyphosate in tea plants (Camellia sinensis)

Tea plant (Camellia sinensis) is an important economical crop that frequently suffers from various herbicides, especially glyphosate. However, the molecular responses and regulatory mechanisms of glyphosate stress in tea plants remain poorly understood. Here, we reported a transcriptome dataset and identified large number of differentially expressed genes (DEGs) under glyphosate exposure. Next, two glutathione S-transferase genes (CsGSTU8-1 and CsGSTU8-2) that upregulated significantly were screened as candidate genes. Tissue-specific expression patterns showed that both CsGSTU8-1 and CsGSTU8-2 had extremely high expression levels in the roots and were predominantly localized in the nucleus and plasma membrane based on subcellular localization. Both were significantly upregulated at different time points under various stressors, including drought, cold, salt, pathogen infections, and SA treatments. An enzymatic activity assay showed that CsGSTU8-1 catalyzes the conjugation of glutathione with 2,4-dinitrochlorobenzene (CDNB). Functional analysis in yeast verified that the two genes significantly contributed to the detoxification of glyphosate, and CsGSTU8-1 had a stronger role in detoxification than CsGSTU8-2. Taken together, these findings provide insights into the molecular responses of tea plants to glyphosate and the functions of CsGSTU8s in glyphosate detoxification, which can be used as a promising genetic resource for improving herbicide resistance in tea cultivars.

Herbicidal Sorbicillinoid Analogs Cause Lignin Accumulation in Aspergillus aculeatus TE-65L

Five new sorbicillinoid derivatives, including (±)-aspersorbicillin A [(±)-1], a pair of enantiomers at C-9, and aspersorbicillins B-D (2-4), together with two known analogs (5 and 6) were isolated from the endophytic fungus Aspergillus aculeatus TE-65L. Their structures including absolute configurations were determined by detailed spectroscopic analyses and electronic circular dichroism calculations. The herbicidal activity of sorbicillinoids on the germ and radicle elongation of various weed types was reported for the first time. Compound 1 displayed significant herbicidal activity against Eleusine indica germ elongation (IC50 = 28.8 μg/mL), while compound 6 inhibited radicle elongation (IC50 = 25.6 μg/mL). Both were stronger than those of glyphosate (66.2 and 30.9 μg/mL, respectively). Further transcriptomic and LC-MS/MS metabolomic analysis indicated that 6 induced the transcriptional expressions of genes related to the lignin biosynthetic pathway, resulting in lignin accumulation. Transmission electron microscopy confirmed the cell wall thickening of seeds treated with 6, suggesting weed growth inhibition. This study reveals new lead compounds for fabricating natural herbicides and expands the agricultural use of sorbicillinoid analogs.

Metabolic resistance mechanism to glufosinate in Eleusine indica

Eleusine indica is one of the most troublesome weeds in farmland worldwide, especially in Citrus Orchard of China. Glufosinate, as an efficient non-selective broad-spectrum herbicide, has been widely utilized for the control of E. indica in Citrus Orchard. The E. indica resistant population (R) was collected from a Citrus Orchard in Yichang City in Hubei province, China. Bioassay experiments showed that the R plants exhibited 3-fold resistance to glufosinate compared with the E. indica susceptible population (S). No known glutamine synthetase (GS) gene mutation associated with glufosinate resistance was found in R plants. And there was also no significant difference in GS activity between R and S plants. Those results indicated that the resistance to glufosinate in R did not involve target-site resistance. However, glutathione S-transferase (GST) inhibitor 4-chloro-7-nitrobenzoxadiazole (NBD-Cl) plus glufosinate gave a better control of R plants compared with glufosinate treatment alone. Moreover, both before and after glufosinate treatment, the GST activity in R plants was significantly higher than that in S plants. By RNA-seq, the expression of GSTU6 and GST4 up-regulated in R plants relative to S plants with or without glufosinate treatment. They were also significantly up-regulated expression in E. indica field resistant populations compared with S population. In summary, the study elucidated that R plants developed metabolic resistance to glufosinate involving GST. And GSTU6 and GST4 genes may play an important role in this glufosinate metabolic resistance. The research results provide a theoretical basis for a deeper understanding of resistance mechanism to glufosinate in E. indica.

Weed biology and management in the multi-omics era: Progress and perspectives

Weeds pose a significant threat to crop production, resulting in substantial yield reduction. In addition, they possess robust weedy traits that enable them to survive in extreme environments and evade human control. In recent years, the application of multi-omics biotechnologies has helped to reveal the molecular mechanisms underlying these weedy traits. In this review, we systematically describe diverse applications of multi-omics platforms for characterizing key aspects of weed biology, including the origins of weed species, weed classification, and the underlying genetic and molecular bases of important weedy traits such as crop-weed interactions, adaptability to different environments, photoperiodic flowering responses, and herbicide resistance. In addition, we discuss limitations to the application of multi-omics techniques in weed science, particularly compared with their extensive use in model plants and crops. In this regard, we provide a forward-looking perspective on the future application of multi-omics technologies to weed science research. These powerful tools hold great promise for comprehensively and efficiently unraveling the intricate molecular genetic mechanisms that underlie weedy traits. The resulting advances will facilitate the development of sustainable and highly effective weed management strategies, promoting greener practices in agriculture.

Comprehensive insights into herbicide resistance mechanisms in weeds: a synergistic integration of transcriptomic and metabolomic analyses

Omics techniques, including genomics, transcriptomics, proteomics, and metabolomics have smoothed the researcher's ability to generate hypotheses and discover various agronomically relevant functions and mechanisms, as well as their implications and associations. With a significant increase in the number of cases with resistance to multiple herbicide modes of action, studies on herbicide resistance are currently one of the predominant areas of research within the field of weed science. High-throughput technologies have already started revolutionizing the current molecular weed biology studies. The evolution of herbicide resistance in weeds (particularly via non-target site resistance mechanism) is a perfect example of a complex, multi-pathway integration-induced response. To date, functional genomics, including transcriptomic and metabolomic studies have been used separately in herbicide resistance research, however there is a substantial lack of integrated approach. Hence, despite the ability of omics technologies to provide significant insights into the molecular functioning of weeds, using a single omics can sometimes be misleading. This mini-review will aim to discuss the current progress of transcriptome-based and metabolome-based approaches in herbicide resistance research, along with their systematic integration.

Cloning and functional characterization of a tau class glutathione transferase associated with haloxyfop-P-methyl resistance in Digitaria sanguinalis

BACKGROUND

Haloxyfop-P-methyl, an acetyl-CoA carboxylase (ACCase)-inhibiting herbicide, has been extensively used to control grass weeds. Widespread use of haloxyfop-P-methyl in cotton fields in China has led to the development of glutathione transferase (GST)-mediated resistance in Digitaria sanguinalis. An RNA-seq analysis identified DsGSTU1, a tau class glutathione transferase from the D. sanguinalis transcriptome as a potential candidate. Here, we cloned DsGSTU1 from D. sanguinalis young leaf tissues and subsequently characterized DsGSTU1 by a combination of sequence analysis, as well as functional heterologous expression in rice.

RESULTS

The full-length coding DNA sequence (CDS) of DsGSTU1 is 717 bp in length. Higher DsGSTU1 expression was observed in haloxyfop-P-methyl-resistant (HR) D. sanguinalis than in haloxyfop-P-methyl-susceptible (HS) plants. Overexpression of the DsGSTU1 gene was confirmed by transformation into the wild-type (WT) Nipponbare rice with pBWA(V)HS, a recombinant expression vector. GST activity in transgenic rice seedlings was 1.18-1.40-fold higher than the WT rice seedlings before and after haloxyfop-P-methyl treatment, respectively. Additionally, transgenic rice seedlings overexpressing DsGSTU1 were less sensitive to haloxyfop-P-methyl.

CONCLUSION

Our combined findings suggest that DsGSTU1 is involved in metabolic resistance to haloxyfop-P-methyl in D. sanguinalis. A better understanding of the major genes contributing to herbicide-resistant D. sanguinalis facilitates the development of resistance management strategies for this global invasive grass weed. © 2023 Society of Chemical Industry.

Phytochemical Cue for the Fitness Costs of Herbicide-Resistant Weeds

Despite increasing knowledge of the fitness costs of viability and fecundity involved in the herbicide-resistant weeds, relatively little is known about the linkage between herbicide resistance costs and phytochemical cues in weed species and biotypes. This study demonstrated relative fitness and phytochemical responses in six herbicide-resistant weeds and their susceptible counterparts. There were significant differences in the parameters of viability (growth and photosynthesis), fecundity fitness (flowering and seed biomass) and a ubiquitous phytochemical (-)-loliolide levels between herbicide-resistant weeds and their susceptible counterparts. Fitness costs occurred in herbicide-resistant Digitaria sanguinalis and Leptochloa chinensis but they were not observed in herbicide-resistant Alopecurus japonicas, Eleusine indica, Ammannia arenaria, and Echinochloa crus-galli. Correlation analysis indicated that the morphological characteristics of resistant and susceptible weeds were negatively correlated with (-)-loliolide concentration, but positively correlated with lipid peroxidation malondialdehyde and total phenol contents. Principal component analysis showed that the lower the (-)-loliolide concentration, the stronger the adaptability in E. crus-galli and E. indica. Therefore, not all herbicide-resistant weeds have fitness costs, but the findings showed several examples of resistance leading to improved fitness even in the absence of herbicides. In particular, (-)-loliolide may act as a phytochemical cue to explain the fitness cost of herbicide-resistant weeds by regulating vitality and fecundity.

Integrated physiological, metabolite and proteomic analysis reveal the glyphosate stress response mechanism in tea plant (Camellia sinensis)

Glyphosate residues can tremendously impact the physiological mechanisms of tea plants, thus threatening tea security and human health. Herein, integrated physiological, metabolite, and proteomic analyses were performed to reveal the glyphosate stress response mechanism in tea plant. After exposure to glyphosate (≥1.25 kg ae/ha), the leaf ultrastructure was damaged, and chlorophyll content and relative fluorescence intensity decreased significantly. The characteristic metabolites catechins and theanine decreased significantly, and the 18 volatile compounds content varied significantly under glyphosate treatments. Subsequently, tandem mass tags (TMT)-based quantitative proteomics was employed to identify the differentially expressed proteins (DEPs) and to validate their biological functions at the proteome level. A total of 6287 proteins were identified and 326 DEPs were screened. These DEPs were mainly catalytic, binding, transporter and antioxidant active proteins, involved in photosynthesis and chlorophyll biosynthesis, phenylpropanoid and flavonoid biosynthesis, sugar and energy metabolism, amino acid metabolism, and stress/defense/detoxification pathway, etc. A total of 22 DEPs were validated by parallel reaction monitoring (PRM), demonstrating that the protein abundances were consistent between TMT and PRM data. These findings contribute to our understanding of the damage of glyphosate to tea leaves and molecular mechanism underlying the response of tea plants to glyphosate.

EST-SSR Markers' Development Based on RNA-Sequencing and Their Application in Population Genetic Structure and Diversity Analysis of Eleusine indica in China

Goosegrass (Eleusine indica) is one of the worst agricultural weeds in China. Molecular markers were developed for genetic diversity and population structure analyses. In this study, we identified 8391 expressed sequence tag-simple sequence repeat (EST-SSR) markers from the de novo assembled unigenes of E. indica. Mononucleotides were the most abundant type of repeats (3591, 42.79%), followed by trinucleotides (3162, 37.68%). The most dominant mononucleotide and trinucleotide repeat motifs were A/T (3406, 40.59%) and AAT/ATT (103, 1.5%), respectively. Fourteen pairs of EST-SSR primers were verified and used to analyze the genetic diversity and population structure of 59 goosegrass populations. A total of 49 alleles were amplified, with the number of alleles (N_a) ranging from two to eleven per locus, and the effective number of alleles (N_e) ranged from 1.07 to 4.53. The average polymorphic information content (PIC) was 0.36. Genetic structure analysis (K = 2) and principal coordinate analysis divided 59 E. indica populations into two groups in a manner similar to the unweighted pair-group method (Dice genetic similarity coefficient = 0.700). This study developed a set of EST-SSR markers in E. indica and successfully analyzed the diversity and population genetic structures of 59 E. indica populations in China.

The moderate oxidative stress induced by glyphosate is not detected in Amaranthus palmeri plants overexpressing EPSPS

The present study aimed to determine whether glyphosate-induced oxidative stress is directly related to the action mechanism of this herbicide (5-enolpyruvylshikimate-3-phosphate synthase or EPSPS inhibition) and analyse the role of oxidative stress in glyphosate toxicity of the weed Amaranthus palmeri S. Wats. Two kinds of populations were studied using EPSPS amplification: glyphosate-sensitive and glyphosate-resistant (by gene amplification). Plants were grown hydroponically and treated with different glyphosate doses, after which several oxidative stress markers were measured in the leaves. Untreated, sensitive and resistant plants showed similar values for the analysed parameters. Treated glyphosate-sensitive plants showed an increase in shikimate, superoxide and H2O2 contents and dose-dependent lipid peroxidation and antioxidant responses; however, none of these effects were observed in resistant plants, indicating that glyphosate-induced oxidative stress is related to EPSPS inhibition. Oxidative stress is associated with an increase in the activity of peroxidases due to EPSPS inhibition, although the link between both processes remains elusive. The fact that some glyphosate doses were lethal but did not induce major oxidative damage provides evidence that glyphosate toxicity is independent of oxidative stress.

Detection of the maximum resistance to the herbicides diuron and glyphosate, and evaluation of its phenotypic cost, in freshwater phytoplankton

One of the most important anthropogenic impacts on freshwater aquatic ecosystems close to intensive agriculture areas is the cumulative increase in herbicide concentrations. The threat is especially relevant for phytoplankton organisms because they have the same physiological targets as the plants for which herbicides have been designed. This led us to explore the evolutionary response of three phytoplanktonic species to increasing concentrations of two herbicides and its consequences in terms of growth and photosynthesis performance. Specifically, we used an experimental ratchet protocol to investigate the differential evolution and the limit of resistance of a cyanobacterium (Microcystis aeruginosa) and two chlorophyceans (Chlamydomonas reinhardtii and Dictyosphaerium chlorelloides) to two herbicides in worldwide use: glyphosate and diuron. Initially, the growth rate of M. aeruginosa and D. chlorelloides was completely inhibited when they were exposed to a dose of 0.23 ppm diuron or 40 ppm glyphosate, whereas a higher concentration of both herbicides (0.46 ppm diuron or 90 ppm glyphosate) was necessary to abolish C. reinhardtii growth. However, after running a ratchet protocol, the resistance of the three species to both herbicides increased by an adaptation process. M. aeruginosa and D. chlorelloides were able to grow at 1.84 ppm diuron and 80 ppm glyphosate and C. reinhardtii proliferated at twice these concentrations. Herbicide-resistant strains showed lower growth rates than their wild-type counterparts in the absence of herbicides, as well as changes on morphology and differences on photosynthetic pigment content. Besides, herbicide-resistant cells generally showed a lower photosynthetic performance than wild-type strains in the three species. These results indicate that the introduction of both herbicides in freshwater ecosystems could produce a diminution of primary production due to the selection of herbicide-resistant mutants, that would exhibit lower photosynthetic performance than wild-type populations.

Characterization of glyphosate and quizalofop-p-ethyl multiple resistance in Eleusine indica

Glyphosate and Acetyl-coenzyme A Carboxylase (ACCase) inhibitors are popular herbicides that control goosegrass. However, some populations are difficult to control due to resistance resulting from the increasing selection pressure. The objectives of this research were to detect the multiple resistance levels, resistance mechanisms, and fitness costs of two goosegrass populations collected in China. The resistance indices of two resistant populations (denominated as R1 and R2) to glyphosate were 3.8 and 2.3, respectively; and it was 18.0 and 14.2 to quizalofop-p-ethyl, respectively. Shikimate accumulation in R1 and R2 populations was only 8% of that of the susceptible population after glyphosate treatment. A Pro-106-Ala mutation in EPSPS and an Asp-2078-Gly mutation in ACCase were present in both resistant populations. Both the expression level of EPSPS and ACCase in resistant populations were similar to that of susceptible populations. The leaf area of the individuals in wild-type populations was more than three times of the leaf area in the resistant populations. Similarly, resistant plants were 45-49% shorter, had 70-76% less fresh shoot weight, and 67-69% fewer seeds than wild-type plants. Goosegrass populations have evolved multiple resistance to glyphosate and the ACCase inhibitor quizalofop-p-ethyl in China. The Pro-106-Ala mutation in the EPSPS and the Asp-2078-Gly mutation in the ACCase were responsible for this resistance. In addition, a fitness cost exists in the resistant populations, and more work should conduct to clear which mutation is responsible for the fitness penalty.

Pro-197-Ser Mutation in ALS and High-Level GST Activities: Multiple Resistance to ALS and ACCase Inhibitors in Beckmannia syzigachne

American sloughgrass (Beckmannia syzigachne Steud.) is one of the most troublesome weeds infesting wheat and canola fields in China. Some biotypes cannot be controlled, either by acetolactate synthase (ALS) or acetyl coenzyme A carboxylase (ACCase) inhibitors, which are the main herbicides for controlling this weed. However, very few studies have investigated multiple resistance mechanism in B. syzigachne. In this study, a B. syzigachne biotype with a high resistance to ALS inhibitors we have reported was also showed relatively lower resistance to ACCase inhibitors, with a resistance index around 7. RNA-seq analysis was used to investigate the factors responsible for multiple resistance, and 60,108 unigenes were assembled by de novo transcriptome assembly and then annotated across eight databases. A Pro-197-Ser mutation was identified in the ALS gene by SNPs analysis and validated by PCR, while no mutation was identified in the ACCase gene. Nineteen candidate metabolic genes were screened and their overexpression was confirmed by qPCR. The expression of GST-T3 and GST-U6 in resistant plants ranged from 7.5- to 109.4-folds than that in susceptible ones at different times after two kinds of herbicide treatment. In addition, GST activities in resistant plants were 3.0-5.0 times higher than that in susceptible plants. Other novel resistance factors also showed high correlation with multiple resistance which included four genes encoding disease resistance proteins, a transcription factor (MYC3), and one gene conferring blight resistance. In this research, a B. syzigachne biotype was confirmed to have evolved multiple resistance to ACCase and ALS inhibitors. The Pro-197-Ser mutation in ALS gene and high-level GST activities were confirmed responsible for the multiple resistance. Characterized disease-resistance proteins, transcription factor, and blight-resistance proteins may play an essential role in these multiple herbicide resistance.

Identification and expression of main genes involved in non-target site resistance mechanisms to fenoxaprop-p-ethyl in Beckmannia syzigachne

BACKGROUND

Non-target-site resistance (NTSR) to herbicides is a serious threat to global agriculture. Although metabolic resistance is the dominant mechanism of NTSR, the molecular mechanisms are not yet well-characterized. This study aimed to uncover the likely metabolism-related genes in Beckmannia syzigachne (American sloughgrass) resistant to fenoxaprop-p-ethyl.

RESULTS

Ultra-performance liquid chromatography - tandem mass spectrometry experiments showed that the resistant American sloughgrass biotype (R, SD-04-SS) showed enhanced degradation of this herbicide compared to the susceptible biotype (S, SD-12). R and S biotype were harvested at 24 h after fenoxaprop-p-ethyl treatment to conduct RNA sequencing (RNA-Seq) analysis to investigate the likely fenoxaprop-p-ethyl metabolic genes. The RNA-Seq libraries yield 417 969 980 clean reads. The de novo assembly generated 115 112 unigenes, of which 57 906 unigenes were annotated. Finally, we identified 273 cytochrome P450s, 178 oxidases, 47 glutathione S-transferases (GSTs), 166 glucosyltransferases (GTs) and 180 ABC transporter genes to determine the likely fenoxaprop-p-ethyl metabolism-related genes in R biotype. Twelve overlapping up-regulated genes in the R biotype (fenoxaprop-p-ethyl-treated R/non-treated R, fenoxaprop-p-ethyl-treated R/fenoxaprop-p-ethyl-treated S) were identified by RNA-Seq and the results were validated using qRT-PCR. Ten were identified as fenoxaprop-p-ethyl metabolism-related genes, including three P450s (homologous to CYP71D7, CYP99A2 and CYP71D10), one GST (homologous to GSTF1), two GTs (homologous to UGT90A1 and UGT83A1) and four oxidase genes.

CONCLUSION

This work demonstrates that the NTSR mechanism by means of enhanced detoxification of fenoxaprop-p-ethyl in American sloughgrass is very likely driven by herbicide metabolism related genes. The RNA-Seq data presented here provide a valuable resource for understanding the molecular mechanism of NTSR in American sloughgrass. © 2020 Society of Chemical Industry.

Glutathione S-transferase: a versatile protein family

Glutathione-S transferase (GST) is a most ancient protein superfamily of multipurpose roles and evolved principally from gene duplication of an ancestral GSH binding protein. They have implemented in diverse plant functions such as detoxification of xenobiotic, secondary metabolism, growth and development, and majorly against biotic and abiotic stresses. The vital structural features of GSTs like highly divergent functional topographies, conserved integrated architecture with separate binding pockets for substrates and ligand, the stringent structural fidelity with high Tm values (50º-60º), and stress-responsive cis-regulatory elements in the promoter region offer this protein as most flexible plant protein for plant breeding approaches, biotechnological applications, etc. This review article summarizes the recent information of GST evolution, and their distribution and structural features with emphasis on the assorted roles of Ser and Cys GSTs with the signature motifs in their active sites, alongside their recent biotechnological application in the area of agriculture, environment, and nanotechnology have been highlighted.

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.

Target-Site and Metabolic Resistance Mechanisms to Penoxsulam in Barnyardgrass (Echinochloa crus-galli (L.) P. Beauv)

Herbicide resistance identification is essential for effective chemical weed control. In this study, we quantified the differences in growth response between penoxsulam resistant (R) and sensitive (S) Echinochloa crus-galli populations, explored the changes in ALS, and performed genetic analyses to identify metabolic genes that are up-regulated by the application of penoxsulam and other common herbicides. The R population showed a 26.0-fold higher resistance to penoxsulam and varied resistance to most tested herbicides with indices ranging from 4.9 to 145.9. A Trp-574-Arg amino acid mutation in ALS and low penoxsulam ALS sensitivity were the main mechanisms underlying herbicide resistance. The penoxsulam resistance can be significantly reversed by two P450s inhibitors and one GST inhibitor. By RNA-Seq, thirty-six highly expressed contigs were selected, and 30 of them were up-regulated in the R population treated by penoxsulam. Many of these genes were significantly expressed when treated with pyroxsulam, metamifop, and quinclorac. These upregulated genes appear to be complementary for plant resistance to penoxsulam and other common herbicides.

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.

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.

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.

Gene expression in response to glyphosate treatment in fleabane (Conyza bonariensis) - glyphosate death response and candidate resistance genes

BACKGROUND

This study takes a whole-transcriptome approach to assess gene expression changes in response to glyphosate treatment in glyphosate-resistant fleabane. We assessed gene expression changes in both susceptible and resistant lines so that the glyphosate death response could be quantified, and constitutively expressed candidate resistance genes identified. There are three copies of the glyphosate target site (5-enolpyruvylshikimate-3-phosphate; EPSPS) gene in Conyza and because Conyza bonariensis is allohexaploid, there is a baseline nine copies of the gene in any individual.

RESULTS

Many genes were differentially expressed in response to glyphosate treatment. Known resistance mutations are present in EPSPS2 but they are present in a glyphosate-susceptible line as well as resistant lines and therefore not sufficient to confer resistance. EPSPS1 is expressed four times more than EPSPS2, further reducing the overall contribution of these mutations.

CONCLUSION

We demonstrate that glyphosate resistance in C. bonariensis is not the result of EPSPS mutations or overexpression, but due to a non-target-site mechanism. A large number of genes are affected by glyphosate treatment. We present a list of candidate non-target-site-resistance (NTSR) genes in fleabane for future studies into these mechanisms. © 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.

Overview of glyphosate-resistant weeds worldwide

Glyphosate is the most widely used and successful herbicide discovered to date, but its utility is now threatened by the occurrence of several glyphosate-resistant weed species. Glyphosate resistance first appeared in Lolium rigidum in an apple orchard in Australia in 1996, ironically the year that the first glyphosate-resistant crop (soybean) was introduced in the USA. Thirty-eight weed species have now evolved resistance to glyphosate, distributed across 37 countries and in 34 different crops and six non-crop situations. Although glyphosate-resistant weeds have been identified in orchards, vineyards, plantations, cereals, fallow and non-crop situations, it is the glyphosate-resistant weeds in glyphosate-resistant crop systems that dominate the area infested and growing economic impact. Glyphosate-resistant weeds present the greatest threat to sustained weed control in major agronomic crops because this herbicide is used to control weeds with resistance to herbicides with other sites of action, and no new herbicide sites of action have been introduced for over 30 years. Industry has responded by developing herbicide resistance traits in major crops that allow existing herbicides to be used in a new way. However, over reliance on these traits will result in multiple-resistance in weeds. Weed control in major crops is at a precarious point, where we must maintain the utility of the herbicides we have until we can transition to new weed management technologies. © 2017 Society of Chemical Industry.

Transcriptomics analysis of the flowering regulatory genes involved in the herbicide resistance of Asia minor bluegrass (Polypogon fugax)

Background

Asia minor bluegrass (Polypogon fugax, P. fugax), a weed that is both distributed across China and associated with winter crops, has evolved resistance to acetyl-CoA carboxylase (ACCase) herbicides, but the resistance mechanism remains unclear. The goal of this study was to analyze the transcriptome between resistant and sensitive populations of P. fugax at the flowering stage.

Results

Populations resistant and susceptible to clodinafop-propargyl showed distinct transcriptome profiles. A total of 206,041 unigenes were identified; 165,901 unique sequences were annotated using BLASTX alignment databases. Among them, 5904 unigenes were classified into 58 transcription factor families. Nine families were related to the regulation of plant growth and development and to stress responses. Twelve unigenes were differentially expressed between the clodinafop-propargyl-sensitive and clodinafop-propargyl-resistant populations at the early flowering stage; among those unigenes, three belonged to the ABI3VP1, BHLH, and GRAS families, while the remaining nine belonged to the MADS family. Compared with the clodinafop-propargyl-sensitive plants, the resistant plants exhibited different expression pattern of these 12 unigenes.

Conclusion

This study identified differentially expressed unigenes related to ACCase-resistant P. fugax and thus provides a genomic resource for understanding the molecular basis of early flowering.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4324-z) contains supplementary material, which is available to authorized users.

Characterization of Eleusine indica with gene mutation or amplification in EPSPS to glyphosate

The evolution of weed-resistant species threatens the sustainable use of glyphosate, which is the most important herbicide widely used in agriculture worldwide. Moreover, the high glyphosate resistance (>180-fold based on LD50) of Eleusine indica found in Malaysia, which carries a double mutation in its 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS), made the control of this species more difficult. By contrast, the same species carrying the same double mutation in EPSPS (T102I+P106S) but found in China only shows a resistance level of not more than 14-fold based on GR50. The resistance level of this population is four times higher than that of the population carrying a single mutation (P106L). Although the members of this population survive under a high glyphosate dosage of 10,080gaeha-1, their growth was significantly inhibited by glyphosate under the recommend dose (840gaeha-1), where in the fresh weight was 85.4% of the control. EPSPS expression, relative copy number, and EPSPS activity in this population were similar to those of the susceptible population. In addition, the expression of two glutathione transferase (GST) genes (GST-U8 and GST-23) and the enzyme activity of the GST in this population did not significantly differ from those of the susceptible population. This finding is important in elucidating the resistance of the naturally evolved glyphosate-resistant (GR) weed species carrying a double mutation in EPSPS to glyphosate.

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.

Pro-106-Ser mutation and EPSPS overexpression acting together simultaneously in glyphosate-resistant goosegrass (Eleusine indica)

Glyphosate has been used for more than 15 years for weed management in citrus groves in the Gulf of Mexico, at up to 3-4 applications per year. Goosegrass (Eleusine indica (L.) Gaertn.) control has sometimes failed. In this research, the mechanisms governing three goosegrass biotypes (Ein-Or from an orange grove, and Ein-Pl1 and Ein-Pl2 from Persian lime groves) with suspected resistance to glyphosate were characterized and compared to a susceptible biotype (Ein-S). Dose-response and shikimate accumulation assays confirmed resistance of the resistant (R) biotypes. There were no differences in glyphosate absorption, but the R biotypes retained up to 62-78% of the herbicide in the treated leaf at 96 h after treatment (HAT), in comparison to the Ein-S biotype (36%). The 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) activity in the Ein-Or and Ein-S biotypes was over 100-fold lower than the Ein-Pl1 and Ein-Pl2 ones. The latter showed a high EPSPS-basal activity, a mutation at Pro-106-Ser position in the EPSPS gene, and EPSPS overexpression. The EPSPS basal and EPSPS overexpression were positively correlated. The R goosegrass biotypes displayed poor glyphosate translocation. Furthermore, this grassweed showed, for the first time, two mechanisms at the target-site level (Pro-106-Ser mutation + EPSPS overexpression) acting together simultaneously against glyphosate.

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.

Transcriptome Profiling to Identify Genes Involved in Mesosulfuron-Methyl Resistance in Alopecurus aequalis

Non-target-site resistance (NTSR) to herbicides is a worldwide concern for weed control. However, as the dominant NTSR mechanism in weeds, metabolic resistance is not yet well-characterized at the genetic level. For this study, we have identified a shortawn foxtail (Alopecurus aequalis Sobol.) population displaying both TSR and NTSR to mesosulfuron-methyl and fenoxaprop-P-ethyl, yet the molecular basis for this NTSR remains unclear. To investigate the mechanisms of metabolic resistance, an RNA-Seq transcriptome analysis was used to find candidate genes that may confer metabolic resistance to the herbicide mesosulfuron-methyl in this plant population. The RNA-Seq libraries generated 831,846,736 clean reads. The de novo transcriptome assembly yielded 95,479 unigenes (averaging 944 bp in length) that were assigned putative annotations. Among these, a total of 29,889 unigenes were assigned to 67 GO terms that contained three main categories, and 14,246 unigenes assigned to 32 predicted KEGG metabolic pathways. Global gene expression was measured using the reads generated from the untreated control (CK), water-only control (WCK), and mesosulfuron-methyl treatment (T) of R and susceptible (S). Contigs that showed expression differences between mesosulfuron-methyl-treated R and S biotypes, and between mesosulfuron-methyl-treated, water-treated and untreated R plants were selected for further quantitative real-time PCR (qRT-PCR) validation analyses. Seventeen contigs were consistently highly expressed in the resistant A. aequalis plants, including four cytochrome P450 monooxygenase (CytP450) genes, two glutathione S-transferase (GST) genes, two glucosyltransferase (GT) genes, two ATP-binding cassette (ABC) transporter genes, and seven additional contigs with functional annotations related to oxidation, hydrolysis, and plant stress physiology. These 17 contigs could serve as major candidate genes for contributing to metabolic mesosulfuron-methyl resistance; hence they deserve further functional study. This is the first large-scale transcriptome-sequencing study to identify NTSR genes in A. aequalis that uses the Illumina platform. This work demonstrates that NTSR is likely driven by the differences in the expression patterns of a set of genes. The assembled transcriptome data presented here provide a valuable resource for A. aequalis biology, and should facilitate the study of herbicide resistance at the molecular level in this and other weed species.