Population Genetic Structure in Glyphosate-Resistant and -Susceptible Palmer Amaranth (Amaranthus palmeri) Populations Using Genotyping-by-sequencing (GBS)

Genomic characterisation and dissection of the onset of resistance to acetyl CoA carboxylase-inhibiting herbicides in a large collection of Digitaria insularis from Brazil

An in-depth genotypic characterisation of a diverse collection of Digitaria insularis was undertaken to explore the neutral genetic variation across the natural expansion range of this weed species in Brazil. With the exception of Minas Gerais, populations from all other states showed high estimates of expected heterozygosity (HE > 0.60) and genetic diversity. There was a lack of population structure based on geographic origin and a low population differentiation between populations across the landscape as evidenced by average Fst value of 0.02. On combining haloxyfop [acetyl CoA carboxylase (ACCase)-inhibiting herbicide] efficacy data with neutral genetic variation, we found evidence of presence of two scenarios of resistance evolution in this weed species. Whilst populations originating from north-eastern region demonstrated an active role of gene flow, populations from the mid-western region displayed multiple, independent resistance evolution as the major evolutionary mechanism. A target-site mutation (Trp2027Cys) in the ACCase gene, observed in less than 1% of resistant populations, could not explain the reduced sensitivity of 15% of the populations to haloxyfop. The genetic architecture of resistance to ACCase-inhibiting herbicides was dissected using a genome wide association study (GWAS) approach. GWAS revealed association of three SNPs with reduced sensitivity to haloxyfop and clethodim. In silico analysis of these SNPs revealed important non-target site genes belonging to families involved in herbicide detoxification, including UDPGT91C1 and GT2, and genes involved in vacuolar sequestration-based degradation pathway. Exploration of five genomic prediction models revealed that the highest prediction power (≥0.80) was achieved with the models Bayes A and RKHS, incorporating SNPs with additive effects and epistatic interactions, respectively.

Repeated evolution of herbicide resistance in Lolium multiflorum revealed by haplotype-resolved analysis of acetyl-CoA carboxylase

Herbicide resistance in weeds is one of the greatest challenges in modern food production. The grass species Lolium multiflorum is an excellent model species to investigate evolution under similar selection pressure because populations have repeatedly evolved resistance to many herbicides, utilizing a multitude of mechanisms to neutralize herbicide damage. In this work, we investigated the gene that encodes acetyl-CoA carboxylase (ACCase), the target site of the most successful herbicide group available for grass weed control. We sampled L. multiflorum populations from agricultural fields with history of intense herbicide use, and studied their response to three ACCase-inhibiting herbicides. To elucidate the mechanisms of herbicide resistance and the genetic relationship among populations, we resolved the haplotypes of 97 resistant and susceptible individuals by sequencing ACCase amplicons using long-read DNA sequencing technologies. Our dose-response data indicated the existence of many, often unpredictable, resistance patterns to ACCase-inhibiting herbicides, where populations exhibited as much as 37-fold reduction in herbicide response. The majority of the populations exhibited resistance to all three herbicides studied. Phylogenetic and molecular genetic analyses revealed multiple evolutionary origins of resistance-endowing ACCase haplotypes, as well as widespread admixture in the region regardless of cropping system. The amplicons generated were diverse, with haplotypes exhibiting 26-110 polymorphisms. Polymorphisms included insertions and deletions 1-31 bp in length, none of which were associated with the resistance phenotype based on an association analysis. We also found evidence that some populations have multiple mechanisms of resistance. Our results highlight the astounding genetic diversity in L. multiflorum populations, and the potential for repeated evolution of herbicide resistance across the landscape that challenges weed management approaches and jeopardizes sustainable weed control practices. We provide an in-depth discussion of the evolutionary and practical implications of our results.

Comparison of microbial community assemblages in the rhizosphere of three Amaranthus spp

Weed management remains a major challenge in cropping systems worldwide, with rising interest in ecological based approaches that can be integrated with herbicide use. Soil microbial communities may play important, yet undiscovered, roles in weed success. Little is known about the rhizosphere communities associated with weeds like Amaranthus, commonly known as pigweeds, and considered some of the most problematic weeds in agricultural systems. In a greenhouse experiment that allowed controlled plant growth conditions and a high number of individual plant specimens to analyze for statistical robustness (n = 8 per species), we show that specific bacterial assemblages form in the rhizospheres of A. retroflexus L. (redroot pigweed), A. palmeri S. Watson (Palmer amaranth), and A. tuberculatus (Moq.) J. D. Sauer (waterhemp). Using a relatively rapid and easy approach of T-RFLP community profiling of the 16S rRNA genes, distinct assemblages corresponded to plant species (PERMANOVA F = 14.776, p = 0.001), and further within each species, similar communities (F = 11.449, p = 0.001) were associated with three rhizosphere soil fractions taken in increasing distances away from the root tissue. These results provide the first solid basis for distinct plant-microbe relationships within three closely related Amaranthus species, warranting closer examination of the identities and function of the microorganisms that appear to be selectively recruited from the extant soil community. More intensive efforts to obtain the microbial taxonomic identities via sequencing are underway that can lead to further detailed studies to elucidate important functional plant-microbe interactions that may associate with weed success. Such data provides underlying key information that may ultimately exploit weed-microbe interactions in development of new integrated weed control tactics.

Weak population differentiation and high diversity in Salsola tragus in the inland Pacific Northwest, USA

BACKGROUND

Salsola tragus is a widespread and problematic weed of semi-arid wheat production globally, and in the inland Pacific Northwest region of the USA. The species exhibits high levels of phenotypic diversity across its range and, at least in California USA, previous work has described cryptic diversity comprising a multi-species complex. Such cryptic diversity could suggest the potential for a differential response to management inputs between groups, and have important implications for the spread of herbicide resistance or other adaptive traits within populations. We used a genotyping-by-sequencing approach to characterize the population structure of S. tragus in the inland Pacific Northwest.

RESULTS

Our results indicated that the population in this region is comprised of a single, tetraploid species (S. tragus sensu latu) with weak population structure on a regional scale. Isolation-by-distance appears to be the primary pattern of structure, but an independent set of weakly differentiated clusters of unknown origin were also apparent, along with a mixed mating system and high levels of largely unstructured genetic diversity.

CONCLUSIONS

Despite considerable phenotypic variability within S. tragus in the region, agronomic weed managers can likely consider it as a single entity across the region, rather than a collection of cryptic subgroups with possible differential responses to management inputs or agroecosystem conditions. A lack of strong barriers to migration and gene flow mean that adaptive traits, such as herbicide resistance, can be expected to spread rapidly through populations across the region. © 2022 Society of Chemical Industry.

Reliable and rapid identification of glyphosate-resistance in the invasive weed Amaranthus palmeri in China

BACKGROUND

Glyphosate-resistant invasive plants, including Amaranthus palmeri S. Watson, have greatly challenged management of new invasions. Elucidating their glyphosate resistance levels rapidly and accurately will better inform management strategies. Quantitative real-time PCR (qPCR) has been used to identify glyphosate resistance in A. palmeri by detecting gene copy numbers of 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS), an enzyme inhibited by glyphosate. However, qPCR can only indirectly determine copy numbers because it requires a calibrator sample; it also lacks standardization, thus limiting its usefulness. Droplet digital PCR (ddPCR) is a new method to detect copy number directly and precisely. We evaluated ddPCR as a tool to determine glyphosate-resistance level while using qPCR and glyphosate dose response (GDR) assays as reference technologies to compare performance and efficiency between methods.

RESULTS

We identified seven susceptible and seven resistant populations of A. palmeri using the GDR assay. Resistant levels detected by qPCR and ddPCR were generally consistent with the GDR results. Although detected values obtained by qPCR and ddPCR were highly correlated (R2  = 0.94), ddPCR results had a lower proportion of non-ideal values (36%) with better accuracy (100%) and specificity (100%) than those of qPCR results.

CONCLUSIONS

Our findings demonstrate that ddPCR offers improved accuracy and specificity in detecting EPSPS gene copy numbers and is a robust and rapid method for glyphosate-resistance identification in A. palmeri. Our research is the first to measure glyphosate resistance in A. palmeri by ddPCR assay and will shed light on future applications of ddPCR in identifying herbicide resistance in other invasive weeds. © 2022 Society of Chemical Industry.

Repeated origins, widespread gene flow, and allelic interactions of target-site herbicide resistance mutations

Causal mutations and their frequency in agricultural fields are well-characterized for herbicide resistance. However, we still lack understanding of their evolutionary history: the extent of parallelism in the origins of target-site resistance (TSR), how long these mutations persist, how quickly they spread, and allelic interactions that mediate their selective advantage. We addressed these questions with genomic data from 19 agricultural populations of common waterhemp (Amaranthus tuberculatus), which we show to have undergone a massive expansion over the past century, with a contemporary effective population size estimate of 8 x 10⁷. We found variation at seven characterized TSR loci, two of which had multiple amino acid substitutions, and three of which were common. These three common resistance variants show extreme parallelism in their mutational origins, with gene flow having shaped their distribution across the landscape. Allele age estimates supported a strong role of adaptation from de novo mutations, with a median age of 30 suggesting that most resistance alleles arose soon after the onset of herbicide use. However, resistant lineages varied in both their age and evidence for selection over two different timescales, implying considerable heterogeneity in the forces that govern their persistence. Two such forces are intra- and inter-locus allelic interactions; we report a signal of extended haplotype competition between two common TSR alleles, and extreme linkage with genome-wide alleles with known functions in resistance adaptation. Together, this work reveals a remarkable example of spatial parallel evolution in a metapopulation, with important implications for the management of herbicide resistance.

The genetic architecture and population genomic signatures of glyphosate resistance in Amaranthus tuberculatus

Much of what we know about the genetic basis of herbicide resistance has come from detailed investigations of monogenic adaptation at known target-sites, despite the increasingly recognized importance of polygenic resistance. Little work has been done to characterize the broader genomic basis of herbicide resistance, including the number and distribution of genes involved, their effect sizes, allele frequencies and signatures of selection. In this work, we implemented genome-wide association (GWA) and population genomic approaches to examine the genetic architecture of glyphosate (Round-up) resistance in the problematic agricultural weed Amaranthus tuberculatus. A GWA was able to correctly identify the known target-gene but statistically controlling for two causal target-site mechanisms revealed an additional 250 genes across all 16 chromosomes associated with non-target-site resistance (NTSR). The encoded proteins had functions that have been linked to NTSR, the most significant of which is response to chemicals, but also showed pleiotropic roles in reproduction and growth. Compared to an empirical null that accounts for complex population structure, the architecture of NTSR was enriched for large effect sizes and low allele frequencies, suggesting the role of pleiotropic constraints on its evolution. The enrichment of rare alleles also suggested that the genetic architecture of NTSR may be population-specific and heterogeneous across the range. Despite their rarity, we found signals of recent positive selection on NTSR-alleles by both window- and haplotype-based statistics, and an enrichment of amino acid changing variants. In our samples, genome-wide single nucleotide polymorphisms explain a comparable amount of the total variation in glyphosate resistance to monogenic mechanisms, even in a collection of individuals where 80% of resistant individuals have large-effect TSR mutations, indicating an underappreciated polygenic contribution to the evolution of herbicide resistance in weed populations.

Defence by duplication: The relation between phenotypic glyphosate resistance and EPSPS gene copy number variation in Amaranthus palmeri

Gene copy number variation (CNV) has been increasingly associated with organismal responses to environmental stress, but we know little about the quantitative relation between CNV and phenotypic variation. In this study we quantify the relation between variation in EPSPS (5-enolpyruvylshikimate-3-phosphate synthase) copy number using digital drop PCR and variation in phenotypic glyphosate resistance in 22 populations of Amaranthus palmeri (Palmer Amaranth), a range-expanding agricultural weed. Overall, we detected a significant positive relation between population mean copy number and resistance. The majority of populations exhibited high glyphosate resistance yet maintained low-resistance individuals, resulting in bimodality in many populations. We also investigated threshold models for the relation between copy number and resistance, and found evidence for a threshold of ~15 EPSPS copies: there was a steep increase in resistance below the threshold, followed by a much shallower increase. Across 924 individuals, as copy number increased the range of variation in resistance decreased, yielding an increasing frequency of high phenotypic resistance individuals. Among populations we detected a decline in variation (s.d.) as mean phenotypic resistance increased from moderate to high, consistent with the prediction that as phenotypic resistance increases in populations, stabilizing selection decreases variation in the trait. Our study demonstrates that populations of A. palmeri can harbour wide variation in EPSPS copy number and phenotypic glyphosate resistance, reflecting the history of, and template for future, resistance evolution.

Investigating the origins and evolution of a glyphosate-resistant weed invasion in South America

The global invasion, and subsequent spread and evolution of weeds provides unique opportunities to address fundamental questions in evolutionary and invasion ecology. Amaranthus palmeri is a widespread glyphosate-resistant (GR) weed in the USA. Since 2015, GR populations of A. palmeri have been confirmed in South America, raising questions about introduction pathways and the importance of pre- vs. post-invasion evolution of GR traits. We used RAD-sequencing genotyping to characterize genetic structure of populations from Brazil, Argentina, Uruguay and the USA. We also quantified gene copy number of the glyphosate target, 5-enolpyruvyl-3-shikimate phosphate synthase (EPSPS), and the presence of an extrachromosomal circular DNA (eccDNA) replicon known to confer glyphosate resistance in USA populations. Populations in Brazil, Argentina and Uruguay were only weakly differentiated (pairwise FST  ≤0.043) in comparison to USA populations (mean pairwise FST  =0.161, range =0.068-0.258), suggesting a single major invasion event. However, elevated EPSPS copy number and the EPSPS replicon were identified in all populations from Brazil and Uruguay, but only in a single Argentinean population. These observations are consistent with independent in situ evolution of glyphosate resistance in Argentina, followed by some limited recent migration of the eccDNA-based mechanism from Brazil to Argentina. Taken together, our results are consistent with an initial introduction of A. palmeri into South America sometime before the 1980s, and local evolution of GR in Argentina, followed by a secondary invasion of GR A. palmeri with the unique eccDNA-based mechanism from the USA into Brazil and Uruguay during the 2010s.

Underutilization Versus Nutritional-Nutraceutical Potential of the Amaranthus Food Plant: A Mini-Review

Amaranthus is a C4 plant tolerant to drought, and plant diseases and a suitable option for climate change. This plant could form part of every region’s cultural heritage and can be transferred to the next generation. Moreover, Amaranthus is a multipurpose plant that has been identified as a traditional edible vegetable endowed with nutritional value, besides its fodder, medicinal, nutraceutical, industrial, and ornamental potentials. In recent decade Amaranthus has received increased research interest. Despite its endowment, there is a dearth of awareness of its numerous potential benefits hence, it is being underutilized. Suitable cultivation systems, innovative processing, and value-adding techniques to promote its utilization are scarce. However, a food-based approach has been suggested as a sustainable measure that tackles food-related problem, especially in harsh weather. Thus, in this review, a literature search for updated progress and potential uses of Amaranthus from online databases of peer-reviewed articles and books was conducted. In addition, the nomenclature, nutritional, and nutraceutical value, was reviewed. The species of focus highlighted in the review include, A. blitum, A. caudatus, A. cruentus, A. dubius, A. hypochondriacus, A. spinosus, A. thunbergii, A. tricolor, and A. viridis.

Recent Discovery of Amaranthus palmeri S. Watson in Italy: Characterization of ALS-Resistant Populations and Sensitivity to Alternative Herbicides

Amaranthus palmeri S. Watson (Amaranthaceae Juss.) is a dioecious noxious weed, native to the Americas, which infests summer crops. It causes high crop losses, and rapidly evolves resistance to herbicides. In Europe, A. palmeri was recorded mostly as a casual alien, but in 2018 it was reported infesting a soybean field in Italy, and the next year two more populations were found in the same area. Experiments were conducted on these three populations to evaluate the resistance to ALS-inhibiting herbicides, to determine the main resistance mechanisms involved and assess the efficacy of alternative herbicides with different sites of action than ALS. The three populations were confirmed cross-resistant to ALS-inhibiting herbicides (thifensulfuron-methyl and imazamox). Gene sequencing identified a Trp to Leu substitution at position 574 of ALS gene in resistant plants, proving that the main resistance mechanism for the three populations is target-site related. The presence of other resistance mechanisms cannot be excluded. Metobromuron, metribuzin and glyphosate are still effective on these populations.

High gene flow maintains genetic diversity following selection for high EPSPS copy number in the weed kochia (Amaranthaceae)

Kochia, a major weed species, has evolved resistance to four herbicide modes of action. Herbicide resistance appears to spread quickly, which could result in diminished standing genetic variation, reducing the ability of populations to adapt further. Here we used double digest restriction enzyme associated sequencing to determine the level of gene flow among kochia populations and whether selection for glyphosate resistance reduces genetic variation. Canadian Prairie populations show little to no genetic differentiation (FST = 0.01) and no correlation between genetic and geographic distance (r2 = - 0.02 p = 0.56), indicating high gene flow and no population structure. There is some evidence that kochia populations are genetically depauperate compared to other weed species, but genetic diversity did not differ between glyphosate susceptible and resistant populations or individuals. The inbreeding coefficients suggest there are 23% fewer heterozygotes in these populations than expected with random mating, and no variation was found within the chloroplast. These results indicate that any alleles for herbicide resistance can be expected to spread quickly through kochia populations, but there is little evidence this spread will reduce the species' genetic variation or limit the species' ability to respond to further selection pressure.

Resistance to Fomesafen, Imazamox and Glyphosate in Euphorbia heterophylla from Brazil

Euphorbia heterophylla is a species of weed that was previously controlled by fomesafen, imazamox and glyphosate, but continued use of these herbicides has selected resistant populations from the Rio Grande do Sul (Brazil). One resistant (R) strain and one susceptible (S) strain to fomesafen, imazamox and glyphosate were compared, the latter by recurrent selection. Dose-response tests showed multiple resistance to these herbicides. The required imazamox concentration to inhibit ALS by 50% was approximately 16 times greater in the R population than in the S population. Based on the EPSPS activity results, the R population was 10 fold less sensitive to glyphosate than the S counterpart. In addition, basal EPSPS activity from R plants was 3.3 fold higher than the level detected on S plants. The Proto IX assays showed high resistance to fomesafen in the R population that accumulated less Proto IX than the S population. Malathion assays showed the participation of CytP450 in fomesafen resistance, but a molecular mechanism could also be involved. To our knowledge, this is the first characterisation of multiple resistance to these three groups of herbicides in E. heterophylla in the world.

Homogeneity among glyphosate-resistant Amaranthus palmeri in geographically distant locations

Since the initial report of glyphosate-resistant (GR) Amaranthus palmeri S. Watson in 2006, resistant populations have been reported in 28 states. The mechanism of resistance is amplification of a 399-kb extrachromosomal circular DNA, called the EPSPS replicon, and is unique to glyphosate-resistant plants. The replicon contains a single copy of the 10-kb 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) gene which causes the concomitant increased expression of EPSP synthase, the target enzyme of glyphosate. It is not known whether the resistance by this amplification mechanism evolved once and then spread across the country or evolved independently in several locations. To compare genomic representation and variation across the EPSPS replicon, whole genome shotgun sequencing (WGS) and mapping of sequences from both GR and susceptible (GS) biotypes to the replicon consensus sequence was performed. Sampling of GR biotypes from AZ, KS, GA, MD and DE and GS biotypes from AZ, KS and GA revealed complete contiguity and deep representation with sequences from GR plants, but lack of homogeneity and contiguity with breaks in coverage were observed with sequences from GS biotypes. The high sequence conservation among GR biotypes with very few polymorphisms which were widely distributed across the USA further supports the hypothesis that glyphosate resistance most likely originated from a single population. We show that the replicon from different populations was unique to GR plants and had similar levels of amplification.

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.

Molecular mechanisms of adaptive evolution revealed by global selection for glyphosate resistance

The human-directed, global selection for glyphosate resistance in weeds has revealed a fascinating diversity of evolved resistance mechanisms, including herbicide sequestration in the vacuole, a rapid cell death response, nucleotide polymorphisms in the herbicide target (5-enolpyruvylshikimate-3-phosphate synthase, EPSPS) and increased gene copy number of EPSPS. For this latter mechanism, two distinct molecular genetic mechanisms have been observed, a tandem duplication mechanism and a large extrachromosomal circular DNA (eccDNA) that is tethered to the chromosomes and passed to gametes at meiosis. These divergent mechanisms have a range of consequences for the spread, fitness, and inheritance of resistance traits, and, particularly in the case of the eccDNA, demonstrate how evolved herbicide resistance can generate new insights into plant adaptation to contemporary environmental stress.

The power and potential of genomics in weed biology and management

There have been previous calls for, and efforts focused on, realizing the power and potential of weed genomics for better understanding of weeds. Sustained advances in genome sequencing and assembly technologies now make it possible for individual research groups to generate reference genomes for multiple weed species at reasonable costs. Here, we present the outcomes from several meetings, discussions, and workshops focused on establishing an International Weed Genomics Consortium (IWGC) for a coordinated international effort in weed genomics. We review the 'state of the art' in genomics and weed genomics, including technologies, applications, and on-going weed genome projects. We also report the outcomes from a workshop and a global survey of the weed science community to identify priority species, key biological questions, and weed management applications that can be addressed through greater availability of, and access to, genomic resources. Major focus areas include the evolution of herbicide resistance and weedy traits, the development of molecular diagnostics, and the identification of novel targets and approaches for weed management. There is increasing interest in, and need for, weed genomics, and the establishment of the IWGC will provide the necessary global platform for communication and coordination of weed genomics research. © 2018 Society of Chemical Industry.