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.

Gene mutations conferring herbicide resistance are hypothesized to have negative pleiotropic effects on plant growth and fitness, which may in turn determine the evolutionary dynamics of herbicide resistance alleles. We used the widespread, annual, diploid grass weed Alopecurus aequalis as a model species to investigate the effect of two resistance mutations—the rare Pro-197-Tyr mutation and the most common mutation, Trp-574-Leu—on acetolactate synthase (ALS) functionality and plant growth. We characterized the enzyme kinetics of ALS from two purified A. aequalis populations, each homozygous for the resistance mutation 197-Tyr or 574-Leu, and assessed the pleiotropic effects of these mutations on plant growth. Both mutations reduced sensitivity of ALS to ALS-inhibiting herbicides without significant changes in extractable ALS activity. The 197-Tyr mutation slightly decreased the substrate affinity (corresponding to an increased K_m for pyruvate) and maximum reaction velocity (V_max) of ALS, whereas the 574-Leu mutation significantly increased these kinetics. Significant decrease or increase in plant growth associated, respectively, with the 197-Tyr and 574-Leu resistance mutations was highly correlated with their impact on ALS kinetics, suggesting more likely persistence of the 574-Leu mutation than the 197-Tyr mutation if herbicide application is discontinued.