RNA-Seq analysis of rye-grass transcriptomic response to an herbicide inhibiting acetolactate-synthase identifies transcripts linked to non-target-site-based resistance

Fenoxaprop-P-ethyl resistance conferred by cytochrome P450s and target site mutation in Alopecurus japonicus

BACKGROUND

Alopecurus japonicus is a serious grass weed species in wheat fields in eastern Asia, and has evolved strong resistance to acetyl-CoA carboxylase (ACCase)-inhibiting herbicides. Although target-site resistance (TSR) to ACCase inhibitors in A. japonicus has been reported, non-target site resistance (NTSR) has not. This study investigated both TSR and NTSR in a fenoxaprop-P-ethyl-resistant A. japonicus population (AHFD-3), which was collected in Feidong County, Anhui Province, China.

RESULTS

We found that AHFD-3 exhibited high resistance to fenoxaprop-P-ethyl and low resistance to flucarbazone-sodium. The sensitivity of AHFD-3 to fenoxaprop-P-ethyl increased significantly after treatment with cytochrome P450 (P450) inhibitors; however, such synergies between P450 inhibitors and fenoxaprop-P-ethyl were not found in two control populations. Sequences of the entire carboxyltransferase domain of A. japonicus ACCase were obtained, and AHFD-3 plants showed an Asp-2078-Gly substitution in the ACCase. With the derived cleaved amplified polymorphic sequence (dCAPS) method, we found that 85.4% of the plants of AHFD-3 carried this mutation. The P450 content in AHFD-3 plants was significantly higher than those of the two control populations after treatment with fenoxaprop-P-ethyl. Ten partial sequences of P450 genes in A. japonicus were cloned. Three P450 genes were up-regulated 12 h after fenoxaprop-P-ethyl treatment, which were all from the P450 subfamily CYP72A. Moreover, a P450 gene from the P450 family CYP81 was up-regulated after fenoxaprop-P-ethyl treatment in all populations studied.

CONCLUSION

Fenoxaprop-P-ethyl resistance in AHFD-3 plants was conferred by up-regulation of cytochrome P450s in the CYP72A subfamily and target site mutation of the ACCase gene. © 2018 Society of Chemical Industry.

Microbe and plant assisted-remediation of organic xenobiotics and its enhancement by genetically modified organisms and recombinant technology: A review

Environmental problems such as the deterioration of groundwater quality, soil degradation and various threats to human, animal and ecosystem health are closely related to the presence of high concentrations of organic xenobiotics in the environment. Employing appropriate technologies to remediate contaminated soils is crucial due to the site-specificity of most remediation methods. The limitations of conventional remediation technologies include poor environmental compatibility, high cost of implementation and poor public acceptability. This raises the call to employ biological methods for remediation. Bioremediation and microbe-assisted bioremediation (phytoremediation) offer many ecological and cost-associated benefits. The overall efficiency and performance of bio- and phytoremediation approaches can be enhanced by genetically modified microbes and plants. Moreover, phytoremediation can also be stimulated by suitable plant-microbe partnerships, i.e. plant-endophytic or plant-rhizospheric associations. Synergistic interactions between recombinant bacteria and genetically modified plants can further enhance the restoration of environments impacted by organic pollutants. Nevertheless, releasing genetically modified microbes and plants into the environment does pose potential risks. These can be minimized by adopting environmental biotechnological techniques and guidelines provided by environmental protection agencies and other regulatory frameworks. The current contribution provides a comprehensive overview on enhanced bioremediation and phytoremediation approaches using transgenic plants and microbes. It also sheds light on the mitigation of associated environmental risks.

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.

Dinitroaniline herbicide resistance in a multiple-resistant Lolium rigidum population

BACKGROUND

The pre-emergence dinitroaniline herbicides (such as trifluralin and pendimethalin) are vital to Australian no-till farming systems. A Lolium rigidum population collected from the Western Australian grain belt with a 12-year trifluralin use history was characterised for resistance to dinitroaniline, acetyl CoA carboxylase (ACCase)- and acetolactate synthase (ALS)-inhibiting herbicides. Target-site resistance mechanisms were investigated.

RESULTS

This L. rigidum population exhibited 32-fold resistance to trifluralin, as compared with the susceptible population. It also displayed 12- to 30-fold cross-resistance to other dinitroaniline herbicides (pendimethalin, ethalfluralin and oryzalin). In addition, this population showed multiple resistance to commonly used post-emergence ACCase- and ALS-inhibiting herbicides. Two target-site α-tubulin gene mutations (Val-202-Phe and Thr-239-Ile) previously documented in other dinitroaniline-resistant weed species were identified, and some known target-site mutations in ACCase (Ile-1781-Leu, Asp-2078-Gly and Cys-2088-Arg) and ALS (Pro-197-Gln/Ser) were found in the same population. An agar-based Petri dish screening method was established for the rapid diagnosis of resistance to dinitroaniline herbicides.

CONCLUSION

Evolution of target-site resistance to both pre- and post-emergence herbicides was confirmed in a single L. rigidum population. The α-tubulin mutations Val-202-Phe and Thr-239-Ile, documented here for the first time in L. rigidum, are likely to be responsible for dinitroaniline resistance in this population. Early detection of dinitroaniline herbicide resistance and integrated weed management strategies are needed to maintain the effectiveness of dinitroaniline herbicides. © 2017 Society of Chemical Industry.

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.

Novel α-Tubulin Mutations Conferring Resistance to Dinitroaniline Herbicides in Lolium rigidum

The dinitroaniline herbicides (particularly trifluralin) have been globally used in many crops for selective grass weed control. Consequently, trifluralin resistance has been documented in several important crop weed species and has recently reached a level of concern in Australian Lolium rigidum populations. Here, we report novel mutations in the L. rigidum α-tubulin gene which confer resistance to trifluralin and other dinitroaniline herbicides. Nucleotide mutations at the highly conserved codon Arg-243 resulted in amino acid substitutions of Met or Lys. Rice calli transformed with the mutant 243-Met or 243-Lys α-tubulin genes were 4- to 8-fold more resistant to trifluralin and other dinitroaniline herbicides (e.g., ethalfluralin and pendimethalin) compared to calli transformed with the wild type α-tubulin gene from L. rigidum. Comprehensive modeling of molecular docking predicts that Arg-243 is close to the trifluralin binding site on the α-tubulin surface and that replacement of Arg-243 by Met/Lys-243 results in a spatial shift of the trifluralin binding domain, reduction of trifluralin-tubulin contacts, and unfavorable interactions. The major effect of these substitutions is a significant rise of free interaction energy between α-tubulin and trifluralin, as well as between trifluralin and its whole molecular environment. These results demonstrate that the Arg-243 residue in α-tubulin is a determinant for trifluralin sensitivity, and the novel Arg-243-Met/Lys mutations may confer trifluralin resistance in L. rigidum.

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.

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.

Fate of Flumioxazin in Aquatic Plants: Two Algae (Pseudokirchneriella subcapitata, Synechococcus sp.), Duckweed (Lemna sp.), and Water Milfoil (Myriophyllum elatinoides)

Flumioxazin separately ¹⁴C-labeled at 1,2-positions of the tetrahydrophthalimide moiety or uniformly labeled at the phenyl ring was exposed to two algae and duckweed via the water layer and water milfoil via the water layer or bottom sediment for 14 days to investigate uptake and metabolic profiles in these aquatic plants. While ¹⁴C-flumioxazin received immediate hydrolysis through maleimide ring opening and amide bond cleavage with its hydrolytic half-life of <1 day in both water and sediment, the ¹⁴C-plant uptake was ≤4.7% of the applied radioactivity (%AR) with water exposure for all plants and 0.9%AR with sediment exposure for water milfoil. No ¹⁴C-translocation between shoot/leaves and roots occurred in water milfoil. The components of ¹⁴C residues in plants were common among the species, which were the above hydrolysates and their transformation products, that is, dicarboxylic acid derivative metabolized via hydroxylation at the double bond of the cyclohexene ring followed by sugar conjugation with its counterpart amine derivative via acid conjugations.

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.

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.

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.

Acetolactate Synthase-Inhibiting Gametocide Amidosulfuron Causes Chloroplast Destruction, Tissue Autophagy, and Elevation of Ethylene Release in Rapeseed

Background: Acetolactate synthase (ALS)-inhibiting herbicides amidosulfuron (Hoestar) is an efficient gametocide that can induce male sterility in rapeseed (Brassica napus L.). We conducted an integrated study of cytological, transcriptomic, and physiological analysis to decipher the gametocidal effect of amidosulfuron. Results: In the first several days after exposure to amidosulfuron at a gametocidal dose of ca. 1 μg per plant, the plants showed the earliest symptoms including short retard of raceme elongation, slight chlorosis on leaf, and decrease of photosynthesis rate. Chloroplasts in leaf and anther epidermis, and tapetal plastids were deformed. Both tapetal cell and uni-nucleate microspore showed autophagic vacuoles and degenerated quickly. The amidosulfuron treatment caused reduction of photosynthetic rate and the contents of leaf chlorophyll, soluble sugar and pyruvate, as well as content alteration of several free amino acids in the treated plants. A comparison of transcriptomic profiling data of the young flower buds of the treated plants with the control identified 142 up-regulated and 201 down-regulated differential expression transcripts with functional annotations. Down-regulation of several interesting genes encoding PAIR1, SDS, PPD2, HFM1, CSTF77, A6, ALA6, UGE1, FLA20, A9, bHLH91, and putative cell wall protein LOC106368794, and up-regulation of autophagy-related protein ATG8A indicated functional abnormalities about cell cycle, cell wall formation, chloroplast structure, and tissue autophagy. Ethylene-responsive transcription factor RAP2-11-like was up-regulated in the flower buds and ethylene release rate was also elevated. The transcriptional regulation in the amidosulfuron-treated plants was in line with the cytological and physiological changes. Conclusions: The results suggested that metabolic decrease related to photosynthesis and energy supply are associated with male sterility induced by amidosulfuron. The results provide insights into the molecular mechanisms of gametocide-induced male sterility and expand the knowledge on the transcriptomic complexity of the plants exposure to sulfonylurea herbicide.

Herbicide Safeners Decrease Sensitivity to Herbicides Inhibiting Acetolactate-Synthase and Likely Activate Non-Target-Site-Based Resistance Pathways in the Major Grass Weed Lolium sp. (Rye-Grass)

Herbicides are currently pivotal to control weeds and sustain food security. Herbicides must efficiently kill weeds while being as harmless as possible for crops, even crops taxonomically close to weeds. To increase their selectivity toward crops, some herbicides are sprayed in association with safeners that are bioactive compounds exacerbating herbicide-degrading pathways reputedly specifically in crops. However, exacerbated herbicide metabolism is also a key mechanism underlying evolved non-target-site-based resistance to herbicides (NTSR) in weeds. This raised the issue of a possible role of safeners on NTSR evolution in weeds. We investigated a possible effect of the respective field rates of the two broadly used safeners cloquintocet-mexyl and mefenpyr-diethyl on the sensitivity of the troublesome global weed Lolium sp. (rye-grass) to the major herbicides inhibiting acetolactate-synthase (ALS) pyroxsulam and iodosulfuron + mesosulfuron, respectively. Three Lolium sp. populations were studied in three series of experiments. The first experiment series compared the frequencies of plants surviving application of each herbicide alone or in association with its safener. Safener co-application caused a net increase ranging from 5.0 to 46.5% in the frequency of plants surviving the field rate of their associated herbicide. In a second series of experiments, safener effect was assessed on individual plant sensitivity using vegetative propagation. A reduction in sensitivity to pyroxsulam and to iodosulfuron + mesosulfuron was observed for 44.4 and 11.1% of the plants in co-treatment with cloquintocet-mexyl and mefenpyr-diethyl, respectively. A third series of experiments investigated safener effect on the expression level of 19 Lolium sp. NTSR marker genes. Safeners showed an enhancing effect on the expression level of 10 genes. Overall, we demonstrated that cloquintocet-mexyl and mefenpyr-diethyl both reduced the sensitivity of Lolium sp. to their associated ALS-inhibiting herbicide and most likely exacerbated herbicide-degrading secondary metabolism pathways. This suggests that genetic variation for safener response is present in Lolium sp. Thus, a possible, uninvestigated way to NTSR evolution could be selection for increased responsiveness to safener action. Delivering safeners exclusively to the crop could mitigate the risk for NTSR evolution in weeds.

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

Glyphosate is an important non-selective herbicide that is in common use worldwide. However, evolved glyphosate-resistant (GR) weeds significantly affect crop yields. Unfortunately, the mechanisms underlying resistance in GR weeds, such as goosegrass (Eleusine indica (L.) Gaertn.), an annual weed found worldwide, have not been fully elucidated. In this study, transcriptome analysis was conducted to further assess the potential mechanisms of glyphosate resistance in goosegrass. The RNA sequencing libraries generated 24 597 462 clean reads. De novo assembly analysis produced 48 852 UniGenes with an average length of 847 bp. All UniGenes were annotated using seven databases. Sixteen candidate differentially expressed genes selected by digital gene expression analysis were validated by quantitative real-time PCR (qRT-PCR). Among these UniGenes, the EPSPS and PFK genes were constitutively up-regulated in resistant (R) individuals and showed a higher copy number than that in susceptible (S) individuals. The expressions of four UniGenes relevant to photosynthesis were inhibited by glyphosate in S individuals, and this toxic response was confirmed by gas exchange analysis. Two UniGenes annotated as glutathione transferase (GST) were constitutively up-regulated in R individuals, and were induced by glyphosate both in R and S. In addition, the GST activities in R individuals were higher than in S. Our research confirmed that two UniGenes (PFK, EPSPS) were strongly associated with target resistance, and two GST-annotated UniGenes may play a role in metabolic glyphosate resistance in goosegrass.

Fitness of ALS-Inhibitors Herbicide Resistant Population of Loose Silky Bentgrass (Apera spica-venti)

Herbicide resistance is an example of plant evolution caused by an increased reliance on herbicides with few sites of action to manage weed populations. This micro-evolutionary process depends on fitness, therefore the assessment of fitness differences between susceptible and resistant populations are pivotal to establish management strategies. Loose silky bentgrass (Apera spica-venti) is a serious weed in Eastern, Northern, and Central Europe with an increasing number of herbicide resistant populations. This study examined the fitness and growth characteristics of an ALS resistant biotype. Fitness and growth characteristics were estimated by comparing seed germination, biomass, seed yield and time to key growth stages at four crop densities of winter wheat (0, 48, 96, and 192 plants m^-2) in a target-neighborhood design. The resistant population germinated 9-20 growing degree days (GDD) earlier than the susceptible population at 10, 16, and 22°C. No differences were observed between resistant and susceptible populations in tiller number, biomass, time to stem elongation, time to first visible inflorescence and seed production. The resistant population reached the inflorescence emergence and flowering stages in less time by 383 and 196 GDD, respectively, at a crop density of 96 winter wheat plants m^-2 with no differences registered at other densities. This study did not observe a fitness cost to herbicide resistance, as often hypothesized. Inversely, a correlation between non-target site resistance (NTSR), earlier germination and earlier flowering time which could be interpreted as fitness benefits as these plant characteristics could be exploited by modifying the timing and site of action of herbicide application to better control ALS NTSR populations of A. spica-venti.

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.

Herbicide-related signaling in plants reveals novel insights for herbicide use strategies, environmental risk assessment and global change assessment challenges

Herbicide impact is usually assessed as the result of a unilinear mode of action on a specific biochemical target with a typical dose-response dynamics. Recent developments in plant molecular signaling and crosstalk between nutritional, hormonal and environmental stress cues are however revealing a more complex picture of inclusive toxicity. Herbicides induce large-scale metabolic and gene-expression effects that go far beyond the expected consequences of unilinear herbicide-target-damage mechanisms. Moreover, groundbreaking studies have revealed that herbicide action and responses strongly interact with hormone signaling pathways, with numerous regulatory protein-kinases and -phosphatases, with metabolic and circadian clock regulators and with oxidative stress signaling pathways. These interactions are likely to result in mechanisms of adjustment that can determine the level of sensitivity or tolerance to a given herbicide or to a mixture of herbicides depending on the environmental and developmental status of the plant. Such regulations can be described as rheostatic and their importance is discussed in relation with herbicide use strategies, environmental risk assessment and global change assessment challenges.

Trends and Challenges in Pesticide Resistance Detection

Pesticide resistance is a crucial factor to be considered when developing strategies for the minimal use of pesticides while maintaining pesticide efficacy. This goal requires monitoring the emergence and development of resistance to pesticides in crop pests. To this end, various methods for resistance diagnosis have been developed for different groups of pests. This review provides an overview of biological, biochemical, and molecular methods that are currently used to detect and quantify pesticide resistance. The agronomic, technical, and economic advantages and drawbacks of each method are considered. Emerging technologies are also described, with their associated challenges and their potential for the detection of resistance mechanisms likely to be selected by current and future plant protection methods.

Utility of a Phylogenetic Perspective in Structural Analysis of CYP72A Enzymes from Flowering Plants

Plant adaptation to external pressures depends on functional diversity in cytochrome P450 (CYP) enzymes. CYPs contain structural domains necessary for the characteristic P450 fold that allows monooxygenation, but they also have great variation in substrate binding affinity. Plant genomes typically contain hundreds of CYPs that contribute to essential functions and species-specific metabolism. The CYP72A subfamily is conserved in angiosperms but its contribution to physiological functions is largely unknown. With genomic information available for many plants, a focused analysis of CYP subfamily diversity is important to understand the contributions of these enzymes to plant evolution. This study examines the extent to which independent gene duplication and evolution have contributed to structural diversification of CYP72A enzymes in different plant lineages. CYP72A genes are prevalent across angiosperms, but the number of genes within each genome varies greatly. The prevalence of CYP72As suggest that the last common ancestor of flowering plants contained a CYP72A sequence, but gene duplication and retention has varied greatly for this CYP subfamily. Sequence comparisons show that CYP72As are involved in species-specific metabolic functions in some plants while there is likely functional conservation between closely related species. Analysis of structural and functional domains within groups of CYP72As reveals clade-specific residues that contribute to functional constraints within subsets of CYP72As. This study provides a phylogenetic framework that allows comparisons of structural features within subsets of the CYP72A subfamily. We examined a large number of sequences from a broad collection of plant species to detect patterns of functional conservation across the subfamily. The evolutionary relationships between CYPs in plant genomes are an important component in understanding the evolution of biochemical diversity in plants.

Multiple Herbicide Resistance in Lolium multiflorum and Identification of Conserved Regulatory Elements of Herbicide Resistance Genes

Herbicide resistance is a ubiquitous challenge to herbicide sustainability and a looming threat to control weeds in crops. Recently four genes were found constituently over-expressed in herbicide resistant individuals of Lolium rigidum, a close relative of Lolium multiflorum. These include two cytochrome P450s, one nitronate monooxygenase and one glycosyl-transferase. Higher expressions of these four herbicide metabolism related (HMR) genes were also observed after herbicides exposure in the gene expression databases, indicating them as reliable markers. In order to get an overview of herbicidal resistance status of L. multiflorum L, 19 field populations were collected. Among these populations, four populations were found to be resistant to acetolactate synthase (ALS) inhibitors while three exhibited resistance to acetyl-CoA carboxylase (ACCase) inhibitors in our initial screening and dose response study. The genotyping showed the presence of mutations Trp-574-Leu and Ile-2041-Asn in ALS and ACCase, respectively, and qPCR experiments revealed the enhanced expression of HMR genes in individuals of certain resistant populations. Moreover, co-expression networks and promoter analyses of HMR genes in O. sativa and A. thaliana resulted in the identification of a cis-regulatory motif and zinc finger transcription factors. The identified transcription factors were highly expressed similar to HMR genes in response to xenobiotics whereas the identified motif is known to play a vital role in coping with environmental stresses and maintaining genome stability. Overall, our findings provide an important step forward toward a better understanding of metabolism-based herbicide resistance that can be utilized to devise novel strategies of weed management.

Target-site and non-target-site based resistance to the herbicide tribenuron-methyl in flixweed (Descurainia sophia L.)

Background

Flixweed (Descurainia sophia L.) is a troublesome and widespread broadleaf weed in winter fields in China, and has evolved high level resistance to acetolactate synthase (ALS)-inhibiting sulfonylurea herbicide tribenuron-methyl.

Results

We identified a resistant flixweed population (N11) exhibiting 116.3-fold resistance to tribenuron-methyl relative to the susceptible population (SD8). Target-site ALS gene mutation Pro-197-Thr was identified in resistant plants. Moreover, the resistance can be reversed to 28.7-fold by the cytochrome P450 inhibitor malathion. The RNA-Sequencing was employed to identify candidate genes involved in non-target-site metabolic resistance in this population. Total 26 differentially expressed contigs were identified and eight of them (four P450s, one ABC transporter, three glycosyltransferase) verified by qRT-PCR. Consistent over-expression of the two contigs homology to CYP96A13 and ABCC1 transporter, respectively, were further qRT-PCR validated using additional plants from the resistant and susceptible populations.

Conclusions

Tribenuron-methyl resistance in flixweed is controlled by target-site ALS mutation and non-target-site based mechanisms. Two genes, CYP96A13 and ABCC1 transporter, could play an important role in metabolic resistance to tribenuron-methyl in the resistant flixweed population and justify further functional studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2915-8) contains supplementary material, which is available to authorized users.

High-throughput sequencing reveals differential regulation of miRNAs in fenoxaprop-P-ethyl-resistant Beckmannia syzigachne

Non-target site resistance (NTSR) to herbicides is an increasing concern for weed control. The majority of previous studies have focused on metabolic resistance mechanisms of NTSR, but no research exists on gene regulation mechanisms behind herbicide resistance, such as microRNA (miRNA). Here, we identified 3 American sloughgrass (Beckmannia syzigachne Steud.) populations containing fenoxaprop-P-ethyl-resistant plants. We then constructed small RNA libraries and subjected them to deep sequencing and bioinformatics analyses. Forty known and 36 potentially novel, predicted miRNAs were successfully identified. Of these, we identified 3 conserved, predicted candidate NTSR-determinant miRNAs and their potential corresponding target genes, as well as 4 novel potential miRNAs with high count. Target gene prediction and annotation indicated that these 7 differentially expressed miRNAs potentially play a role in regulating specific stress-responsive genes, very likely related to herbicide resistance. Expression profiles were determined with quantitative real-time PCR. The present study is a novel, large-scale characterization of weed miRNAs. The results should further our understanding of miRNA expression profiles associated with herbicide resistance, allowing for the development of more effective weed management strategies.

Establishing a herbicide-metabolizing enzyme library in Beckmannia syzigachne to identify genes associated with metabolic resistance

Non-target site resistance (NTSR) to herbicides is an increasing concern for weed control. Metabolic herbicide resistance is an important mechanism for NTSR. However, little is known about metabolic resistance at the genetic level. In this study, we have identified three fenoxaprop-P-ethyl-resistant American sloughgrass (Beckmannia syzigachne Steud.) populations, in which the molecular basis for NTSR remains unclear. To reveal the mechanisms of metabolic resistance, the genes likely to be involved in herbicide metabolism (e.g. for cytochrome P450s, esterases, hydrolases, oxidases, peroxidases, glutathione S-transferases, glycosyltransferases, and transporter proteins) were isolated using transcriptome sequencing, in combination with RT-PCR (reverse transcription-PCR) and RACE (rapid amplification of cDNA ends). Consequently, we established a herbicide-metabolizing enzyme library containing at least 332 genes, and each of these genes was cloned and the sequence and the expression level compared between the fenoxaprop-P-ethyl-resistant and susceptible populations. Fifteen metabolic enzyme genes were found to be possibly involved in fenoxaprop-P-ethyl resistance. In addition, we found five metabolizing enzyme genes that have a different gene sequence in plants of susceptible versus resistant B. syzigachne populations. These genes may be major candidates for herbicide metabolic resistance. This established metabolic enzyme library represents an important step forward towards a better understanding of herbicide metabolism and metabolic resistance in this and possibly other closely related weed species. This new information may help to understand weed metabolic resistance and to develop novel strategies of weed management.

Genome-Wide Transcriptional Profiling and Metabolic Analysis Uncover Multiple Molecular Responses of the Grass Species Lolium perenne Under Low-Intensity Xenobiotic Stress

Lolium perenne, which is a major component of pastures, lawns, and grass strips, can be exposed to xenobiotic stresses due to diffuse and residual contaminations of soil. L. perenne was recently shown to undergo metabolic adjustments in response to sub-toxic levels of xenobiotics. To gain insight in such chemical stress responses, a de novo transcriptome analysis was carried out on leaves from plants subjected at the root level to low levels of xenobiotics, glyphosate, tebuconazole, and a combination of the two, leading to no adverse physiological effect. Chemical treatments influenced significantly the relative proportions of functional categories and of transcripts related to carbohydrate processes, to signaling, to protein-kinase cascades, such as Serine/Threonine-protein kinases, to transcriptional regulations, to responses to abiotic or biotic stimuli and to responses to phytohormones. Transcriptomics-based expressions of genes encoding different types of SNF1 (sucrose non-fermenting 1)-related kinases involved in sugar and stress signaling or encoding key metabolic enzymes were in line with specific qRT-PCR analysis or with the important metabolic and regulatory changes revealed by metabolomic analysis. The effects of pesticide treatments on metabolites and gene expression strongly suggest that pesticides at low levels, as single molecule or as mixture, affect cell signaling and functioning even in the absence of major physiological impact. This global analysis of L. perenne therefore highlighted the interactions between molecular regulation of responses to xenobiotics, and also carbohydrate dynamics, energy dysfunction, phytohormones and calcium signaling.

ALOMYbase, a resource to investigate non-target-site-based resistance to herbicides inhibiting acetolactate-synthase (ALS) in the major grass weed Alopecurus myosuroides (black-grass)

Background

Herbicide resistance in agrestal weeds is a global problem threatening food security. Non-target-site resistance (NTSR) endowed by mechanisms neutralising the herbicide or compensating for its action is considered the most agronomically noxious type of resistance. Contrary to target-site resistance, NTSR mechanisms are far from being fully elucidated. A part of weed response to herbicide stress, NTSR is considered to be largely driven by gene regulation. Our purpose was to establish a transcriptome resource allowing investigation of the transcriptomic bases of NTSR in the major grass weed Alopecurus myosuroides L. (Poaceae) for which almost no genomic or transcriptomic data was available.

Results

RNA-Seq was performed from plants in one F2 population that were sensitive or expressing NTSR to herbicides inhibiting acetolactate-synthase. Cloned plants were sampled over seven time-points ranging from before until 73 h after herbicide application. Assembly of over 159M high-quality Illumina reads generated a transcriptomic resource (ALOMYbase) containing 65,558 potentially active contigs (N50 = 1240 nucleotides) predicted to encode 32,138 peptides with 74 % GO annotation, of which 2017 were assigned to protein families presumably involved in NTSR. Comparison with the fully sequenced grass genomes indicated good coverage and correct representation of A. myosuroides transcriptome in ALOMYbase. The part of the herbicide transcriptomic response common to the resistant and the sensitive plants was consistent with the expected effects of acetolactate-synthase inhibition, with striking similarities observed with published Arabidopsis thaliana data. A. myosuroides plants with NTSR were first affected by herbicide action like sensitive plants, but ultimately overcame it. Analysis of differences in transcriptomic herbicide response between resistant and sensitive plants did not allow identification of processes directly explaining NTSR. Five contigs associated to NTSR in the F2 population studied were tentatively identified. They were predicted to encode three cytochromes P450 (CYP71A, CYP71B and CYP81D), one peroxidase and one disease resistance protein.

Conclusions

Our data confirmed that gene regulation is at the root of herbicide response and of NTSR. ALOMYbase proved to be a relevant resource to support NTSR transcriptomic studies, and constitutes a valuable tool for future research aiming at elucidating gene regulations involved in NTSR in A. myosuroides.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1804-x) contains supplementary material, which is available to authorized users.