Characterization of the horseweed (Conyza canadensis) transcriptome using GS-FLX 454 pyrosequencing and its application for expression analysis of candidate non-target herbicide resistance genes

Physiological and biochemical characteristics and microbial responses of Medicago sativa (Fabales: Fabaceae) varieties with different resistance to atrazine stress

Atrazine, a commonly employed herbicide for corn production, can leave residues in soil, resulting in photosynthetic toxicity and impeding growth in subsequent alfalfa (Medicago sativa L.) crops within alfalfa-corn rotation systems. The molecular regulatory mechanisms by which atrazine affects alfalfa growth and development, particularly its impact on the microbial communities of the alfalfa rhizosphere, are not well understood. This study carried out field experiments to explore the influence of atrazine stress on the biomass, chlorophyll content, antioxidant system, and rhizosphere microbial communities of the atrazine-sensitive alfalfa variety WL-363 and the atrazine-resistant variety JN5010. The results revealed that atrazine significantly reduced WL-363 growth, decreasing plant height by 8.58 cm and root length by 5.42 cm (p < 0.05). Conversely, JN5010 showed minimal reductions, with decreases of 1.96 cm in height and 1.26 cm in root length. Chlorophyll content in WL-363 decreased by 35% under atrazine stress, while in JN5010, it was reduced by only 10%. Reactive oxygen species (ROS) accumulation increased by 60% in WL-363, compared to a 20% increase in JN5010 (p < 0.05 for both). Antioxidant enzyme activities, such as superoxide dismutase (SOD) and catalase (CAT), were significantly elevated in JN5010 (p < 0.05), suggesting a more robust defense mechanism. Although the predominant bacterial and fungal abundances in rhizosphere soils remained generally unchanged under atrazine stress, specific microbial groups exhibited variable responses. Notably, Promicromonospora abundance declined in WL-363 but increased in JN5010. FAPROTAX functional predictions indicated shifts in the abundance of microorganisms associated with pesticide degradation, resistance, and microbial structure reconstruction under atrazine stress, displaying different patterns between the two varieties. This study provides insights into how atrazine residues affect alfalfa rhizosphere microorganisms and identifies differential microbial responses to atrazine stress, offering valuable reference data for screening and identifying atrazine-degrading bacteria.

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.

Multiple Resistance Mechanisms Involved in Glyphosate Resistance in Eleusine indica

Glyphosate is a non-selective herbicide and is widely used for weed control in non-cultivated land in China. One susceptible (S) and five putative glyphosate-resistant (R1, R2, R3, R4, and R5) Eleusine indica biotypes were selected to investigate their resistance levels and the potential resistance mechanisms. Based on the dose-response assays, the R3 and R5 biotypes showed a low-level (2.4 to 3.5-fold) glyphosate resistance, and the R1, R2, and R4 biotypes exhibited a moderate- to high-level (8.6 to 19.2-fold) resistance, compared with the S biotype. The analysis of the target-site resistance (TSR) mechanism revealed that the P106A mutation and the heterozygous double T102I + P106S mutation were found in the R3 and R4 biotypes, respectively. In addition, the similar EPSPS gene overexpression was observed in the R1, R2, and R5 biotypes, suggesting that additional non-target-site resistance (NTSR) mechanisms may contribute to glyphosate resistance in R1 and R2 biotypes. Subsequently, an RNA-Seq analysis was performed to identify candidate genes involved in NTSR. In total, ten differentially expressed contigs between untreated S and R1 or R2 plants, and between glyphosate-treated S and R1 or R2 plants, were identified and further verified with RT-qPCR. One ATP-binding cassette (ABC) transporter gene, one aldo-keto reductases (AKRs) gene and one cytochrome P450 monooxygenase (CytP450) gene were up-regulated in R1 or R2 plants. These results indicated that EPSPS overexpression, single or double mutation was a common TSR mechanisms in E. indica. Additional NTSR mechanisms could play an essential role in glyphosate resistance. Three genes, ABCC4, AKR4C10, and CYP88, could serve as important candidate genes and deserve further functional studies.

DNA Methylation of ABC Transporters Differs in Native and Non-native Populations of Conyza canadensis L

While differences in the methylation patterns of ABC transporters under different environmental conditions and their role in plant growth, development, and response to biotic and abiotic stresses are well documented, less is known about the variation in the methylation patterns of ABC transporters in plant species in the native and non-native ranges. In this study, we present the results of differences in methylation of ABC transporters of Conyza canadensis L. in its native (North America) and non-native (Kashmir Himalaya) ranges. Our data show that ABC transporter genes have reduced DNA methylation in Kashmir Himalaya than in North America. Furthermore, in the non-native range of Kashmir Himalaya, we found that ABC transporter genes have enriched RNA Pol-II binding and reduced nucleosome occupancy, both hallmarks of transcriptional activity. Taken together, our study showed differential DNA methylation in the ABC transporter genes in the native range of North America and non-native range of Kashmir Himalaya in Conyza canadensis and that the reduced DNA methylation and increased RNA Pol-II binding is one of the possible mechanisms through which this species in the non-native range of Kashmir Himalaya may show greater gene expression of ABC transporter genes. This increased ABC transporter gene expression may help the plant to grow in different environmental conditions in the non-native range. Furthermore, this study could pave way for more studies to better explain the enigmatic plant invasions of C. canadensis in the non-native range of Kashmir Himalaya.

Multiple evolutionary origins of glyphosate resistance in Lolium multiflorum

The multitude of herbicide resistance patterns that have evolved in different weed species is a remarkable example of the rapid adaptation to anthropogenic-driven disturbance. Recently, resistance to glyphosate was identified in multiple populations of Lolium multiflorum in Oregon. We used phenotypic approaches, as well as population genomic and gene expression analyses, to determine whether known mechanisms were responsible for glyphosate resistance and whether resistance phenotypes evolved independently in different populations, and to identify potential loci contributing to resistance. We found no evidence of genetic alterations or expression changes at known target and non-target-site resistance mechanisms of glyphosate. Population genomic analyses indicated that resistant populations tended to have largely distinct ancestry from one another, suggesting that glyphosate resistance did not spread among populations by gene flow. Rather, resistance appears to have evolved independently on different genetic backgrounds. We also detected potential loci associated with the resistance phenotype, some of which encode proteins with potential effects on herbicide metabolism. Our results suggest that Oregon populations of L. multiflorum evolved resistance to glyphosate due to a novel mechanism. Future studies that characterize the gene or genes involved in resistance will be necessary to confirm this conclusion.

Resistance to mesosulfuron-methyl in Beckmannia syzigachne may involve ROS burst and non-target-site resistance mechanisms

Herbicide resistance to chemical herbicide is a global issue that presents an ongoing threat to grain production. Though it has been frequently implicated that the production of detoxification enzymes increased in resistance development, the mechanisms for overexpression of these genes employed by herbicide-resistant weeds remain complicated. In this study, a mesosulfuron-methyl resistant Beckmannia syzigachne population (R) was found to be cross-resistant to another herbicide pyriminobac-methyl. No known target-site mutations were detected in the R population. In contrast, the decreased uptake and enhanced metabolic rates of mesosulfuron-methyl were detected in the R than the susceptible (S) population. Two candidate ATP-binding cassette (ABC) transporter genes (ABCB25 and ABCC14) that were constitutively up-regulated in the R population were identified by RNA-sequencing and validated by RT-qPCR. Alteration of antioxidant enzyme activities and gene expressions implied that mesosulfuron-methyl-induced antioxidant defenses provoked reactive oxygen species (ROS) burst. ROS scavenger assay showed that ROS induces ABCB25 and ABCC14 expression. This study reported for the first time that ABC transporters mediated non-target-site resistance contributes to mesosulfuron-methyl resistance in a B. syzigachne population, and implicated that ROS burst might be involved in the overexpression of ABC transporter genes in weeds.

Novel Candidate Genes Differentially Expressed in Glyphosate-Treated Horseweed (Conyza canadensis)

The evolution of herbicide-resistant weed species is a serious threat for weed control. Therefore, we need an improved understanding of how gene regulation confers herbicide resistance in order to slow the evolution of resistance. The present study analyzed differentially expressed genes after glyphosate treatment on a glyphosate-resistant Tennessee ecotype (TNR) of horseweed (Conyza canadensis), compared to a susceptible biotype (TNS). A read size of 100.2 M was sequenced on the Illumina platform and subjected to de novo assembly, resulting in 77,072 gene-level contigs, of which 32,493 were uniquely annotated by a BlastX alignment of protein sequence similarity. The most differentially expressed genes were enriched in the gene ontology (GO) term of the transmembrane transport protein. In addition, fifteen upregulated genes were identified in TNR after glyphosate treatment but were not detected in TNS. Ten of these upregulated genes were transmembrane transporter or kinase receptor proteins. Therefore, a combination of changes in gene expression among transmembrane receptor and kinase receptor proteins may be important for endowing non-target-site glyphosate-resistant C. canadensis.

Molecular mechanisms underlying glyphosate resistance in bacteria

Glyphosate is a nonselective herbicide that kills weeds and other plants competing with crops. Glyphosate specifically inhibits the 5-enolpyruvyl-shikimate-3-phosphate (EPSP) synthase, thereby depleting the cell of EPSP serving as a precursor for biosynthesis of aromatic amino acids. Glyphosate is considered to be toxicologically safe for animals and humans. Therefore, it became the most-important herbicide in agriculture. However, its intensive application in agriculture is a serious environmental issue because it may negatively affect the biodiversity. A few years after the discovery of the mode of action of glyphosate, it has been observed that bacteria evolve glyphosate resistance by acquiring mutations in the EPSP synthase gene, rendering the encoded enzyme less sensitive to the herbicide. The identification of glyphosate-resistant EPSP synthase variants paved the way for engineering crops tolerating increased amounts of the herbicide. This review intends to summarize the molecular mechanisms underlying glyphosate resistance in bacteria. Bacteria can evolve glyphosate resistance by (i) reducing glyphosate sensitivity or elevating production of the EPSP synthase, by (ii) degrading or (iii) detoxifying glyphosate and by (iv) decreasing the uptake or increasing the export of the herbicide. The variety of glyphosate resistance mechanisms illustrates the adaptability of bacteria to anthropogenic substances due to genomic alterations.

An ABCC-type transporter endowing glyphosate resistance in plants

Glyphosate is the most widely used herbicide in world agriculture and for general vegetation control in a wide range of situations. Global and often intensive glyphosate selection of very large weedy plant populations has resulted in widespread glyphosate resistance evolution in populations of many weed species. Here, working with a glyphosate-resistant (GR) Echinochloa colona population that evolved in a Western Australia agricultural field, we identified an ATP-binding cassette (ABC) transporter (EcABCC8) that is consistently up-regulated in GR plants. When expressed in transgenic rice, this EcABCC8 transporter endowed glyphosate resistance. Equally, rice, maize, and soybean overexpressing the EcABCC8 ortholog genes were made resistant to glyphosate. Conversely, CRISPR/Cas9-mediated knockout of the EcABCC8 ortholog gene OsABCC8 increased rice susceptibility to glyphosate. Subcellular localization analysis and quantification of glyphosate cellular levels in treated ABCC8 transgenic rice plants and isolated leaf protoplasts as well as structural modeling support that EcABCC8 is likely a plasma membrane-localized transporter extruding cytoplasmic glyphosate to the apoplast, lowering the cellular glyphosate level. This is a report of a membrane transporter effluxing glyphosate in a GR plant species, and its function is likely conserved in crop plant species.

The Genetic Architecture of Emerging Fungicide Resistance in Populations of a Global Wheat Pathogen

Containing fungal diseases often depends on the application of fungicidal compounds. Fungicides can rapidly lose effectiveness due to the rise of resistant individuals in populations. However, the lack of knowledge about resistance mutations beyond known target genes challenges investigations into pathways to resistance. We used whole-genome sequencing data and association mapping to reveal the multilocus genetic architecture of fungicide resistance in a global panel of 159 isolates of Parastagonospora nodorum, an important fungal pathogen of wheat. We found significant differences in azole resistance among global field populations. The populations evolved distinctive combinations of resistance alleles which can interact when co-occurring in the same genetic background. We identified 34 significantly associated single nucleotide polymorphisms located in close proximity to genes associated with fungicide resistance in other fungi, including a major facilitator superfamily transporter. Using fungal colony growth rates and melanin production at different temperatures as fitness proxies, we found no evidence that resistance was constrained by genetic trade-offs. Our study demonstrates how genome-wide association studies of a global collection of pathogen strains can recapitulate the emergence of fungicide resistance. The distinct complement of resistance mutations found among populations illustrates how the evolutionary trajectory of fungicide adaptation can be complex and challenging to predict.

Mechanisms of evolved herbicide resistance

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

A chromosome-scale draft sequence of the Canada fleabane genome

BACKGROUND

Due to the accessibility of underlying technologies the 'Omics', in particular genomics, are becoming commonplace in several fields of research, including the study of agricultural pests. The weed community is starting to embrace these approaches; genome sequences have been made available in the past years, with several other sequencing projects underway, as promoted by the International Weed Genome Consortium. Chromosome-scale sequences are essential to fully exploit the power of genetics and genomics.

RESULTS

We report such an assembly for Conyza canadensis, an important agricultural weed. Third-generation sequencing technology was used to create a genome assembly of 426 megabases, of which nine chromosome-scale scaffolds cover more than 98% of the entire assembled sequence. As this weed was the first to be identified with glyphosate resistance, and since we do not have a firm handle on the genetic mechanisms responsible for several herbicide resistances in the species, the genome sequence was annotated with genes known to be associated with herbicide resistance. A high number of ABC-type transporters, cytochrome P450 and glycosyltransferases (159, 352 and 181, respectively) were identified among the list of ab initio predicted genes.

CONCLUSION

As C. canadensis has a small genome that is syntenic with other Asteraceaes, has a short life cycle and is relatively easy to cross, it has the potential to become a model weed species and, with the chromosome-scale genome sequence, contribute to a paradigm shift in the way non-target site resistance is studied. © 2020 Her Majesty the Queen in Right of CanadaPest Management Science © 2020 Society of Chemical Industry.

Target-site EPSPS Pro-106-Ser mutation in Conyza canadensis biotypes with extreme resistance to glyphosate in Ohio and Iowa, USA

Documenting the diversity of mechanisms for herbicide resistance in agricultural weeds is helpful for understanding evolutionary processes that contribute to weed management problems. More than 40 species have evolved resistance to glyphosate, and at least 13 species have a target-site mutation at position 106 of EPSPS. In horseweed (Conyza canadensis), this p106 mutation has only been reported in Canada. Here, we sampled seeds from one plant (= biotype) at 24 sites in Ohio and 20 in Iowa, screened these biotypes for levels of resistance, and sequenced their DNA to detect the p106 mutation. Resistance categories were based on 80% survival at five glyphosate doses: S (0×), R1 (1×), R2 (8×), R3 (20×), or R4 (40×). The p106 mutation was not found in the19 biotypes scored as S, R1, or R2, while all 25 biotypes scored as R3 or R4 had the same proline-to-serine substitution at p106. These findings represent the first documented case of target-site mediated glyphosate resistance in horseweed in the United States, and the first to show that this mutation was associated with very strong resistance. We hypothesize that the p106 mutation has occurred multiple times in horseweed and may be spreading rapidly, further complicating weed management efforts.

Parallel and nonparallel genomic responses contribute to herbicide resistance in Ipomoea purpurea, a common agricultural weed

The repeated evolution of herbicide resistance has been cited as an example of genetic parallelism, wherein separate species or genetic lineages utilize the same genetic solution in response to selection. However, most studies that investigate the genetic basis of herbicide resistance examine the potential for changes in the protein targeted by the herbicide rather than considering genome-wide changes. We used a population genomics screen and targeted exome re-sequencing to uncover the potential genetic basis of glyphosate resistance in the common morning glory, Ipomoea purpurea, and to determine if genetic parallelism underlies the repeated evolution of resistance across replicate resistant populations. We found no evidence for changes in 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS), glyphosate's target protein, that were associated with resistance, and instead identified five genomic regions that showed evidence of selection. Within these regions, genes involved in herbicide detoxification-cytochrome P450s, ABC transporters, and glycosyltransferases-are enriched and exhibit signs of selective sweeps. One region under selection shows parallel changes across all assayed resistant populations whereas other regions exhibit signs of divergence. Thus, while it appears that the physiological mechanism of resistance in this species is likely the same among resistant populations, we find patterns of both similar and divergent selection across separate resistant populations at particular loci.

Non-Target-Site Resistance to Herbicides: Recent Developments

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

Population Genomic Approaches for Weed Science

Genomic approaches are opening avenues for understanding all aspects of biological life, especially as they begin to be applied to multiple individuals and populations. However, these approaches typically depend on the availability of a sequenced genome for the species of interest. While the number of genomes being sequenced is exploding, one group that has lagged behind are weeds. Although the power of genomic approaches for weed science has been recognized, what is needed to implement these approaches is unfamiliar to many weed scientists. In this review we attempt to address this problem by providing a primer on genome sequencing and provide examples of how genomics can help answer key questions in weed science such as: (1) Where do agricultural weeds come from; (2) what genes underlie herbicide resistance; and, more speculatively, (3) can we alter weed populations to make them easier to control? This review is intended as an introduction to orient weed scientists who are thinking about initiating genome sequencing projects to better understand weed populations, to highlight recent publications that illustrate the potential for these methods, and to provide direction to key tools and literature that will facilitate the development and execution of weed genomic projects.

Evolutionary and ecological insights from herbicide-resistant weeds: what have we learned about plant adaptation, and what is left to uncover?

The evolution of herbicide resistance in crop weeds presents one of the greatest challenges to agriculture and the production of food. Herbicide resistance has been studied for more than 60 yr, in the large part by researchers seeking to design effective weed control programs. As an outcome of this work, various unique questions in plant adaptation have been addressed. Here, I collate recent research on the herbicide-resistant problem in light of key questions and themes in evolution and ecology. I highlight discoveries made on herbicide-resistant weeds in three broad areas - the genetic basis of adaptation, evolutionary constraints, experimental evolution - and similarly discuss questions left to be answered. I then develop how one would use herbicide-resistance evolution as a model for studying eco-evolutionary dynamics within a community context. My overall goals are to highlight important findings in the weed science literature that are relevant to themes in plant adaptation and to stimulate the use of herbicide-resistant plants as models for addressing key questions within ecology and evolution.

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

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

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

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

EPSPS duplication and mutation involved in glyphosate resistance in the allotetraploid weed species Poa annua L

BACKGROUND

Poa annua is a widespread winter annual weed species in California. Recently, poor control of this species with glyphosate was reported by growers in an almond orchard in California with a history of repetitive glyphosate use. The objectives of this research were to evaluate the level of glyphosate resistance in a developed S4 P. annua line (R) and identify the mechanisms of resistance involved.

RESULTS

Whole-plant dose-response experiments confirmed glyphosate resistance in R, which required 18-fold more glyphosate to achieve a 50% growth reduction compared with a susceptible line (S), results that were supported by the lower shikimate accumulation observed in R compared with S. No differences in glyphosate absorption, translocation, or metabolism were observed, suggesting that non-target-site mechanisms of resistance are not involved in the resistance phenotype. A missense single nucleotide polymorphism was observed in EPSPS coding position 106 in R, resulting in a leucine to proline substitution. This polymorphism was observed exclusively in P. supina EPSPS homeologs. A seven-fold increase in the number of copies of EPSPS alleles was observed in R compared with S.

CONCLUSIONS

We report the first case of glyphosate resistance associated with both EPSPS duplication and target-site mutation at position 106, leading to high levels of glyphosate resistance in the allotetraploid weed species Poa annua L. Data obtained in this research will be useful for the development of diagnostic tools for rapid glyphosate resistance identification, monitoring and containment. © 2018 Society of Chemical Industry.

A novel triple amino acid substitution in the EPSPS found in a high-level glyphosate-resistant Amaranthus hybridus population from Argentina

BACKGROUND

The evolution of herbicide-resistant weeds is one of the most important concerns of global agriculture. Amaranthus hybridus L. is a competitive weed for summer crops in South America. In this article, we intend to unravel the molecular mechanisms by which an A. hybridus population from Argentina has become resistant to extraordinarily high levels of glyphosate.

RESULTS

The glyphosate-resistant population (A) exhibited particularly high parameters of resistance (GR50  = 20 900 g ai ha-1 , Rf = 314), with all plants completing a normal life cycle even after 32X dose application. No shikimic acid accumulation was detected in the resistant plants at any of the glyphosate concentrations tested. Molecular and genetic analyses revealed a novel triple substitution (TAP-IVS: T102I, A103V, and P106S) in the 5-enol-pyruvylshikimate-3-phosphate synthase (EPSPS) enzyme of population A and an incipient increase on the epsps relative copy number but without effects on the epsps transcription levels. The novel mechanism was prevalent, with 48% and 52% of the individuals being homozygous and heterozygous for the triple substitution, respectively. In silico conformational studies revealed that TAP-IVS triple substitution would generate an EPSPS with a functional active site but with an increased restriction to glyphosate binding.

CONCLUSION

The prevalence of the TAP-IVS triple substitution as the sole mechanism detected in the highly glyphosate resistant population suggests the evolution of a new glyphosate resistance mechanism arising in A. hybridus. This is the first report of a naturally occurring EPSPS triple substitution and the first glyphosate target-site resistance mechanism described in A. hybridus. © 2018 Society of Chemical Industry.

Sequencing the Plastid Genome of Giant Ragweed (Ambrosia trifida, Asteraceae) From a Herbarium Specimen

We report the first plastome sequence of giant ragweed (Ambrosia trifida); with this new genome information, we assessed the phylogeny of Asteraceae and the transcriptional profiling against glyphosate resistance in giant ragweed. Assembly and genic features show a normal angiosperm quadripartite plastome structure with no signatures of deviation in gene directionality. Comparative analysis revealed large inversions across the plastome of giant ragweed and the previously sequenced members of the plant family. Asteraceae plastid genomes contain two inversions of 22.8 and 3.3 kb; the former is located between trnS-GCU and trnG-UCC genes, and the latter between trnE-UUC and trnT-GGU genes. The plastid genome sequences of A. trifida and the related species, Ambrosia artemisiifolia, are identical in gene content and arrangement, but they differ in length. The phylogeny is well-resolved and congruent with previous hypotheses about the phylogenetic relationship of Asteraceae. Transcriptomic analysis revealed divergence in the relative expressions at the exonic and intronic levels, providing hints toward the ecological adaptation of the genus. Giant ragweed shows various levels of glyphosate resistance, with introns displaying higher expression patterns at resistant time points after the assumed herbicide treatment.

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.

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.

High Levels of Glyphosate Resistance in Conyza canadensis from Agricultural and Non-Agricultural Sites in Ohio and Iowa

Glyphosate is an important herbicide worldwide, but its efficacy has been compromised where weed species have evolved glyphosate resistance. To better understand evolutionary outcomes of continued and strong selection from glyphosate exposure, we characterized variation in resistance in self-pollinating Conyza canadensis (horseweed) in Ohio and Iowa, where glyphosate resistance was first reported in 2002 and 2011, respectively. In 2015, we collected seeds from a total of 74 maternal plants (biotypes) from no-till soybean fields vs. non-agricultural sites in each state, using one representative plant per site. Young plants from each biotype were sprayed with glyphosate rates of 0x, 1x (840 g ae ha-1), 8x, 20x, or 40x. Resistant biotypes with at least 80% survival at each dosage were designated as R1 (1x), R2 (8x), R3 (20x), or R4 (40x). Nearly all Ohio agricultural biotypes were R4, as were 62% of biotypes from the non-agricultural sites. In Iowa, R4 biotypes were clustered in the southeastern soybean fields, where no-till agriculture is more common, and 45% of non-agricultural biotypes were R1-R4. Our results show that resistance levels to glyphosate can be very high (at least 40x) in both states, and that non-agricultural sites likely serve as a refuge for glyphosate-resistant biotypes.

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

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

Investigation of glyphosate resistance levels and target-site based resistance (TSR) mechanisms in Conyza canadensis (L.) from apple orchards around areas of Bohai seas and Loess Plateau in China

The resistance levels to glyphosate and target-site based resistance mechanisms in susceptible (S) and resistant (R) Conyza canadensis (L.) populations, which were collected from apple orchards around areas of Bohai seas and Loess Plateau in China, were investigated. Among forty C. canadensis populations, eighteen populations (45%) were still susceptible; fourteen populations (35%) evolved low resistance levels resistance to glyphosate with resistance index (RI) of 2.02 to 3.90. In contrast, eight populations (20%) evolved medium resistance levels with RI of 4.35 to 8.38. The shikimic acid concentrations in R populations were highly negative relative with the glyphosate resistance levels in C. canadensis, the Pearson correlation coefficient was -0.82 treated by glyphosate at 1.8mg/L. Three 5-enoylpyruvylshikimate 3'-phosphate synthase genes (EPSPS1, EPSPS2 and EPSPS3) were cloned in all S and glyphosate-resistant C. canadensis populations. No amino acid substitution was identified at site of 102 and 106 in three EPSPS genes, which were reported to confer glyphosate resistance in other weed species. The relative expression level of EPSPS mRNA in R populations (SD07, LN05, SHX06 and SD09) was 4.5 to 13.2 times higher than in S biotype. The Pearson correlation coefficient between EPSPS expression levels and RI was 0.79, which indicated the over expression of EPSPS mRNA may cause these R populations evolve higher resistance level to glyphosate.

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

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

Foliar Glyphosate Treatment Alters Transcript and Hormone Profiles in Crown Buds of Leafy Spurge and Induces Dwarfed and Bushy Phenotypes throughout its Perennial Lifecycle

Leafy spurge ( L.) is an invasive weed of North America and its perennial nature attributed to underground adventitious buds (UABs) that undergo seasonal cycles of para-, endo-, and ecodormancy. Recommended rates of glyphosate (∼1 kg ha) destroy aboveground shoots but plants still regenerate vegetatively; therefore, it is considered glyphosate-tolerant. However, foliar application of glyphosate at higher rates (2.2-6.7 kg ha) causes sublethal effects that induce UABs to produce stunted, bushy phenotypes. We investigated the effects of glyphosate treatment (±2.24 kg ha) on vegetative growth, phytohormone, and transcript profiles in UABs under controlled environments during one simulated seasonal cycle. Because shoots derived from UABs of foliar glyphosate-treated plants produced stunted, bushy phenotypes, we could not directly determine if these UABs transitioned through seasonally induced endo- and ecodormancy. However, transcript abundance for leafy spurge dormancy marker genes and principal component analyses suggested that UABs of foliar glyphosate-treated plants transitioned through endo- and ecodormancy. Glyphosate treatment increased shikimate abundance in UABs 7 d after treatment; however, the abundance of shikimate gradually decreased as UABs transitioned through endo- and ecodormancy. The dissipation of shikimate over time suggests that glyphosate's target site was no longer affected, but these changes did not reverse the altered phenotypes observed from UABs of foliar glyphosate-treated leafy spurge. Transcript profiles further indicated that foliar glyphosate treatment significantly affected phytohormone biosynthesis and signaling, particularly auxin transport; gibberellic acid, abscisic acid and jasmonic acid biosynthesis; ethylene responses; and detoxification and cell cycle processes in UABs. These results correlated well with the available phytohormone profiles and altered phenotypes.

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.

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

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

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.

Comparative proteomic analysis of horseweed (Conyza canadensis) biotypes identifies candidate proteins for glyphosate resistance

Emergence of glyphosate-resistant horseweed (Conyza canadensis) biotypes is an example of how unrelenting use of a single mode of action herbicide in agricultural weed control drives genetic adaptation in targeted species. While in other weeds glyphosate resistance arose from target site mutation or target gene amplification, the resistance mechanism in horseweed uses neither of these, being instead linked to reduced herbicide uptake and/or translocation. The molecular components underpinning horseweed glyphosate-resistance remain unknown. Here, we used an in vitro leaf disc system for comparative analysis of proteins extracted from control and glyphosate-treated tissues of glyphosate-resistant and glyphosate-susceptible biotypes. Analysis of shikimic acid accumulation, ABC-transporter gene expression, and cell death were used to select a suitable glyphosate concentration and sampling time for enriching proteins pivotal to glyphosate resistance. Protein gel analysis and mass spectrometry identified mainly chloroplast proteins differentially expressed between the biotypes before and after glyphosate treatment. Chloroplasts are the organelles in which the shikimate pathway, which is targeted by glyphosate, is located. Calvin cycle enzymes and proteins of unknown function were among the proteins identified. Our study provides candidate proteins that could be pivotal in engendering resistance and implicates chloroplasts as the primary sites driving glyphosate-resistance in horseweed.

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.

Nucleotide Diversity at Site 106 of EPSPS in Lolium perenne L. ssp. multiflorum from California Indicates Multiple Evolutionary Origins of Herbicide Resistance

The repeated evolution of herbicide resistance in weeds is an ongoing problem in agricultural regions across the world, and presents a unique system in which to study the origins and spread of adaptive traits across heterogeneous landscapes. Lolium perenne ssp. multiflorum (Lam.) (Italian ryegrass) is a widespread grass weed of agricultural crops that has repeatedly evolved resistance to herbicides across the world. In California, resistance to glyphosate has become increasingly common. To identify the mechanisms conferring glyphosate resistance in California populations of L. perenne and to gain insights into the evolutionary origins and spread of resistance in the region, we investigated the frequency of target-site mutations conferring resistance to glyphosate combined with the frequency of resistant individuals in 14 populations. A region of 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) was sequenced in 401 individuals to assay for target site mutations. Seven unique alleles were detected at codon site 106, four of which have been previously shown to confer target-site-based resistance to glyphosate. Four different resistance alleles were detected, indicating that resistance to glyphosate has evolved multiple times in the region. In two populations, no EPSPS mutations were detected despite the presence of resistant plants, strongly suggesting that non-target-site-based mechanisms confer resistance to glyphosate in these populations. It is likely that resistance to glyphosate in these 14 California populations of L. perenne derives from at least five evolutionary origins, indicating that adaptive traits can evolve repeatedly over agricultural landscapes.

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.

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.

Diversity of ABC transporter genes across the plant kingdom and their potential utility in biotechnology

Background

The ATP-binding cassette (ABC) transporter gene superfamily is ubiquitous among extant organisms and prominently represented in plants. ABC transporters act to transport compounds across cellular membranes and are involved in a diverse range of biological processes. Thus, the applicability to biotechnology is vast, including cancer resistance in humans, drug resistance among vertebrates, and herbicide and other xenobiotic resistance in plants. In addition, plants appear to harbor the highest diversity of ABC transporter genes compared with any other group of organisms. This study applied transcriptome analysis to survey the kingdom-wide ABC transporter diversity in plants and suggest biotechnology applications of this diversity.

Results

We utilized sequence similarity-based informatics techniques to infer the identity of ABC transporter gene candidates from 1295 phylogenetically-diverse plant transcriptomes. A total of 97,149 putative (approximately 25 % were full-length) ABC transporter gene members were identified; each RNA-Seq library (plant sample) had 88 ± 30 gene members. As expected, simpler organisms, such as algae, had fewer unique members than vascular land plants. Differences were also noted in the richness of certain ABC transporter subfamilies. Land plants had more unique ABCB, ABCC, and ABCG transporter gene members on average (p < 0.005), and green algae, red algae, and bryophytes had significantly more ABCF transporter gene members (p < 0.005). Ferns had significantly fewer ABCA transporter gene members than all other plant groups (p < 0.005).

Conclusions

We present a transcriptomic overview of ABC transporter gene members across all major plant groups. An increase in the number of gene family members present in the ABCB, ABCC, and ABCD transporter subfamilies may indicate an expansion of the ABC transporter superfamily among green land plants, which include all crop species. The striking difference between the number of ABCA subfamily transporter gene members between ferns and other plant taxa is surprising and merits further investigation. Discussed is the potential exploitation of ABC transporters in plant biotechnology, with an emphasis on crops.

Electronic supplementary material

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

Environmental Conditions Influence Induction of Key ABC-Transporter Genes Affecting Glyphosate Resistance Mechanism in Conyza canadensis

Conyza canadensis has been reported to be the most frequent weed species that evolved resistance to glyphosate in various parts of the world. The objective of the present study was to investigate the effect of environmental conditions (temperature and light) on the expression levels of the EPSPS gene and two major ABC-transporter genes (M10 and M11) on glyphosate susceptible (GS) and glyphosate resistant (GR) horseweed populations, collected from several regions across Greece. Real-time PCR was conducted to determine the expression level of the aforementioned genes when glyphosate was applied at normal (1×; 533 g·a.e.·ha(-1)) and high rates (4×, 8×), measured at an early one day after treatment (DAT) and a later stage (four DAT) of expression. Plants were exposed to light or dark conditions, at three temperature regimes (8, 25, 35 °C). GR plants were made sensitive when exposed to 8 °C with light; those sensitized plants behaved biochemically (shikimate accumulation) and molecularly (expression of EPSPS and ABC-genes) like the GS plants. Results from the current study show the direct link between the environmental conditions and the induction level of the above key genes that likely affect the efficiency of the proposed mechanism of glyphosate resistance.

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.

Molecular Basis for Resistance Against Phosphonate Antibiotics and Herbicides

Research in recent years have illuminated data on the mechanisms and targets of phosphonic acid antibiotics and herbicides, including fosfomycin, glyphosate, fosmidomycin and FR900098. Here we review the current state of knowledge of the structural and biochemical characterization of resistance mechanisms against these bioactive natural products. Advances in the understanding of these resistance determinants have spurred knowledge-based campaigns aimed towards the design of derivatives that retain biological activity but are less prone to tolerance.

454 Pyrosequencing of Olive (Olea europaea L.) Transcriptome in Response to Salinity

Olive (Olea europaea L.) is one of the most important crops in the Mediterranean region. The expansion of cultivation in areas irrigated with low quality and saline water has negative effects on growth and productivity however the investigation of the molecular basis of salt tolerance in olive trees has been only recently initiated. To this end, we investigated the molecular response of cultivar Kalamon to salinity stress using next-generation sequencing technology to explore the transcriptome profile of olive leaves and roots and identify differentially expressed genes that are related to salt tolerance response. Out of 291,958 obtained trimmed reads, 28,270 unique transcripts were identified of which 35% are annotated, a percentage that is comparable to similar reports on non-model plants. Among the 1,624 clusters in roots that comprise more than one read, 24 were differentially expressed comprising 9 down- and 15 up-regulated genes. Respectively, inleaves, among the 2,642 clusters, 70 were identified as differentially expressed, with 14 down- and 56 up-regulated genes. Using next-generation sequencing technology we were able to identify salt-response-related transcripts. Furthermore we provide an annotated transcriptome of olive as well as expression data, which are both significant tools for further molecular studies in olive.

Utilizing next-generation sequencing to study homeologous polymorphisms and herbicide-resistance-endowing mutations in Poa annua acetolactate synthase genes

BACKGROUND

Detection of single nucleotide polymorphisms (SNPs) related to herbicide resistance in non-model polyploid weed species is fraught with difficulty owing to the gene duplication and lack of reference sequences. Our research seeks to overcome these obstacles by Illumina HiSeq read mapping, SNP calling and allele frequency determinations. Our focus is on the acetolactate synthase (ALS) gene, the target site of ALS-inhibiting herbicides, in Poa annua, an allotetraploid weed species originating from two diploid parents, P. supina and P. infirma.

RESULTS

ALS contigs with complete coding regions of P. supina, P. infirma and P. annua were assembled and compared with ALS genes from other plant species. The ALS infirma-homeolog of P. annua showed higher levels of nucleotide sequence variability than the supina-homeolog. Comparisons of read mappings of P. annua and a simulated P. supina × P. infirma hybrid showed high resemblance. Two homeolog-specific primer pairs were designed and used to amplify a 1860 bp region covering all resistance-conferring codons in the ALS gene. Four P. annua populations, GN, RB, GW and LG, showed high resistance to two ALS inhibitors, bispyribac-sodium and foramsulfuron, and two populations, HD and RS, showed lower resistance in the rate-response trial. Mutations conferring Trp-574-Leu substitution were observed in the infirma-homeolog of GN and RB and in the supina-homeolog of GW and LG, but no resistance-conferring mutation was observed in the two populations of lower resistance, HD and RS.

CONCLUSION

In this study we have demonstrated the use of NGS data to study homeologous polymorphisms, parentage and herbicide resistance in an allotetraploid weed species, P. annua. Complete coding sequences of the ALS gene were assembled for P. infirma, P. supina, infirma-homeolog and supina-homeolog in P. annua. A pipeline consisting of read mapping, SNP calling and allele frequency calculation was developed to study the parentage of P. annua, which provided a new perspective on this topic besides the views of morphology, karyotype and phylogeny. Our two homeolog-specific primer pairs can be utilized in future research to separate the homeologs of the ALS gene in P. annua and cover all the codons that have been reported to confer herbicide resistance.

Concerted action of target-site mutations and high EPSPS activity in glyphosate-resistant junglerice (Echinochloa colona) from California

BACKGROUND

Echinochloa colona is an annual weed affecting field crops and orchards in California. An E. colona population carrying a mutation in the EPSPS gene endowing resistance to glyphosate, the most widely used non-selective herbicide, was recently identified in the Northern Sacramento Valley of California. Plants from this population, from a suspected glyphosate-resistant (GR) population, and from one susceptible (S) population collected in the Northern Sacramento Valley of California, were used to generate three GR and one S selfed lines to study possible mechanisms involved in glyphosate resistance.

RESULTS

Based on the amount of glyphosate required to kill 50% of the plants (LD50 ), GR lines were 4-9-fold more resistant than S plants and accumulated less shikimate after glyphosate treatment. GR and S lines did not differ in glyphosate absorption, translocation or metabolism. A different target-site mutation was found in each of two of the GR lines corresponding to Pro106Thr and Pro106Ser substitutions; the mutations were found in different homoeologous EPSPS genes. No mutation was found in the third GR line, which exhibited 1.4-fold higher basal EPSPS activity and a fivefold greater LD50 than S plants. Quantitative RT-PCR revealed that GR lines had similar or lower EPSPS expression than S plants.

CONCLUSION

It is demonstrated that individuals with different glyphosate resistance mechanisms can coexist in the same population, individuals from different populations may carry different resistance mechanisms and different mechanisms can act in concert within single E. colona plants. However, other plant factors or resistance mechanisms appear to modulate plant expression of EPSPS sensitivity to glyphosate.

Glyphosate's impact on vegetative growth in leafy spurge identifies molecular processes and hormone cross-talk associated with increased branching

BACKGROUND

Leafy spurge (Euphorbia esula) is a perennial weed that is considered glyphosate tolerant, which is partially attributed to escape through establishment of new vegetative shoots from an abundance of underground adventitious buds. Leafy spurge plants treated with sub-lethal concentrations of foliar-applied glyphosate produce new vegetative shoots with reduced main stem elongation and increased branching. Processes associated with the glyphosate-induced phenotype were determined by RNAseq using aerial shoots derived from crown buds of glyphosate-treated and -untreated plants. Comparison between transcript abundance and accumulation of shikimate or phytohormones (abscisic acid, auxin, cytokinins, and gibberellins) from these same samples was also done to reveal correlations.

RESULTS

Transcriptome assembly and analyses confirmed differential abundance among 12,918 transcripts (FDR ≤ 0.05) and highlighted numerous processes associated with shoot apical meristem maintenance and stem growth, which is consistent with the increased number of actively growing meristems in response to glyphosate. Foliar applied glyphosate increased shikimate abundance in crown buds prior to decapitation of aboveground shoots, which induces growth from these buds, indicating that 5-enolpyruvylshikimate 3-phosphate (EPSPS) the target site of glyphosate was inhibited. However, abundance of shikimate was similar in a subsequent generation of aerial shoots derived from crown buds of treated and untreated plants, suggesting EPSPS is no longer inhibited or abundance of shikimate initially observed in crown buds dissipated over time. Overall, auxins, gibberellins (precursors and catabolites of bioactive gibberellins), and cytokinins (precursors and bioactive cytokinins) were more abundant in the aboveground shoots derived from glyphosate-treated plants.

CONCLUSION

Based on the overall data, we propose that the glyphosate-induced phenotype resulted from complex interactions involving shoot apical meristem maintenance, hormone biosynthesis and signaling (auxin, cytokinins, gibberellins, and strigolactones), cellular transport, and detoxification mechanisms.

Glyphosate's impact on vegetative growth in leafy spurge identifies molecular processes and hormone cross-talk associated with increased branching

BACKGROUND

Leafy spurge (Euphorbia esula) is a perennial weed that is considered glyphosate tolerant, which is partially attributed to escape through establishment of new vegetative shoots from an abundance of underground adventitious buds. Leafy spurge plants treated with sub-lethal concentrations of foliar-applied glyphosate produce new vegetative shoots with reduced main stem elongation and increased branching. Processes associated with the glyphosate-induced phenotype were determined by RNAseq using aerial shoots derived from crown buds of glyphosate-treated and -untreated plants. Comparison between transcript abundance and accumulation of shikimate or phytohormones (abscisic acid, auxin, cytokinins, and gibberellins) from these same samples was also done to reveal correlations.

RESULTS

Transcriptome assembly and analyses confirmed differential abundance among 12,918 transcripts (FDR ≤ 0.05) and highlighted numerous processes associated with shoot apical meristem maintenance and stem growth, which is consistent with the increased number of actively growing meristems in response to glyphosate. Foliar applied glyphosate increased shikimate abundance in crown buds prior to decapitation of aboveground shoots, which induces growth from these buds, indicating that 5-enolpyruvylshikimate 3-phosphate (EPSPS) the target site of glyphosate was inhibited. However, abundance of shikimate was similar in a subsequent generation of aerial shoots derived from crown buds of treated and untreated plants, suggesting EPSPS is no longer inhibited or abundance of shikimate initially observed in crown buds dissipated over time. Overall, auxins, gibberellins (precursors and catabolites of bioactive gibberellins), and cytokinins (precursors and bioactive cytokinins) were more abundant in the aboveground shoots derived from glyphosate-treated plants.

CONCLUSION

Based on the overall data, we propose that the glyphosate-induced phenotype resulted from complex interactions involving shoot apical meristem maintenance, hormone biosynthesis and signaling (auxin, cytokinins, gibberellins, and strigolactones), cellular transport, and detoxification mechanisms.

Protein kinase GCN2 mediates responses to glyphosate in Arabidopsis

BACKGROUND

The increased selection pressure of the herbicide glyphosate has played a role in the evolution of glyphosate-resistance in weedy species, an issue that is becoming a threat to global agriculture. The molecular components involved in the cellular toxicity response to this herbicide at the expression level are still unidentified.

RESULTS

In this study, we identify the protein kinase GCN2 as a cellular component that fosters the action of glyphosate in the model plant Arabidopsis thaliana. Comparative studies using wild-type and gcn2 knock-out mutant seedlings show that the molecular programme that the plant deploys after the treatment with the herbicide, is compromised in gcn2. Moreover, gcn2 adult plants show a lower inhibition of photosynthesis, and both seedlings and adult gcn2 plants accumulate less shikimic acid than wild-type after treatment with glyphosate.

CONCLUSIONS

These results points to an unknown GCN2-dependent factor involved in the cascade of events triggered by glyphosate in plants. Data suggest either that the herbicide does not equally reach the target-enzyme in a gcn2 background, or that a decreased flux in the shikimate pathway in a gcn2 plants minimize the impact of enzyme inhibition.

De novo genome assembly of the economically important weed horseweed using integrated data from multiple sequencing platforms

Horseweed (Conyza canadensis), a member of the Compositae (Asteraceae) family, was the first broadleaf weed to evolve resistance to glyphosate. Horseweed, one of the most problematic weeds in the world, is a true diploid (2n = 2x = 18), with the smallest genome of any known agricultural weed (335 Mb). Thus, it is an appropriate candidate to help us understand the genetic and genomic bases of weediness. We undertook a draft de novo genome assembly of horseweed by combining data from multiple sequencing platforms (454 GS-FLX, Illumina HiSeq 2000, and PacBio RS) using various libraries with different insertion sizes (approximately 350 bp, 600 bp, 3 kb, and 10 kb) of a Tennessee-accessed, glyphosate-resistant horseweed biotype. From 116.3 Gb (approximately 350× coverage) of data, the genome was assembled into 13,966 scaffolds with 50% of the assembly = 33,561 bp. The assembly covered 92.3% of the genome, including the complete chloroplast genome (approximately 153 kb) and a nearly complete mitochondrial genome (approximately 450 kb in 120 scaffolds). The nuclear genome is composed of 44,592 protein-coding genes. Genome resequencing of seven additional horseweed biotypes was performed. These sequence data were assembled and used to analyze genome variation. Simple sequence repeat and single-nucleotide polymorphisms were surveyed. Genomic patterns were detected that associated with glyphosate-resistant or -susceptible biotypes. The draft genome will be useful to better understand weediness and the evolution of herbicide resistance and to devise new management strategies. The genome will also be useful as another reference genome in the Compositae. To our knowledge, this article represents the first published draft genome of an agricultural weed.

Glyphosate resistance: state of knowledge

Studies of mechanisms of resistance to glyphosate have increased current understanding of herbicide resistance mechanisms. Thus far, single-codon non-synonymous mutations of EPSPS (5-enolypyruvylshikimate-3-phosphate synthase) have been rare and, relative to other herbicide mode of action target-site mutations, unconventionally weak in magnitude for resistance to glyphosate. However, it is possible that weeds will emerge with non-synonymous mutations of two codons of EPSPS to produce an enzyme endowing greater resistance to glyphosate. Today, target-gene duplication is a common glyphosate resistance mechanism and could become a fundamental process for developing any resistance trait. Based on competition and substrate selectivity studies in several species, rapid vacuole sequestration of glyphosate occurs via a transporter mechanism. Conversely, as the chloroplast requires transporters for uptake of important metabolites, transporters associated with the two plastid membranes may separately, or together, successfully block glyphosate delivery. A model based on finite glyphosate dose and limiting time required for chloroplast loading sets the stage for understanding how uniquely different mechanisms can contribute to overall glyphosate resistance.