De novo genome assembly of white clover (Trifolium repens L.) reveals the role of copy number variation in rapid environmental adaptation.

BACKGROUND: White clover (Trifolium repens) is a globally important perennial forage legume. This species also serves as an eco-evolutionary model system for studying within-species chemical defense variation; it features a well-studied polymorphism for cyanogenesis (HCN release following tissue damage), with higher frequencies of cyanogenic plants favored in warmer locations worldwide. Using a newly generated haplotype-resolved genome and two other long-read assemblies, we tested the hypothesis that copy number variants (CNVs) at cyanogenesis genes play a role in the ability of white clover to rapidly adapt to local environments. We also examined questions on subgenome evolution in this recently evolved allotetraploid species and on chromosomal rearrangements in the broader IRLC legume clade. RESULTS: Integration of PacBio HiFi, Omni-C, Illumina, and linkage map data yielded a completely de novo genome assembly for white clover (created without a priori sequence assignment to subgenomes). We find that white clover has undergone extensive transposon diversification since its origin but otherwise shows highly conserved genome organization and composition with its diploid progenitors. Unlike some other clover species, its chromosomal structure is conserved with other IRLC legumes. We further find extensive evidence of CNVs at the major cyanogenesis loci; these contribute to quantitative variation in the cyanogenic phenotype and to local adaptation across wild North American populations. CONCLUSIONS: This work provides a case study documenting the role of CNVs in local adaptation in a plant species, and it highlights the value of pan-genome data for identifying contributions of structural variants to adaptation in nature.

Variable expression of cyanide detoxification and tolerance genes in cyanogenic and acyanogenic white clover (Trifolium repens).

PREMISE: ß-Cyanoalanine synthase (ß-CAS) and alternative oxidase (AOX) play important roles in the ability of plants to detoxify and tolerate hydrogen cyanide (HCN). These functions are critical for all plants because HCN is produced at low levels during basic metabolic processes, and especially for cyanogenic species, which release high levels of HCN following tissue damage. However, expression of ß-CAS and Aox genes has not been examined in cyanogenic species, nor compared between cyanogenic and acyanogenic genotypes within a species. METHODS: We used a natural polymorphism for cyanogenesis in white clover to examine ß-CAS and Aox gene expression in relation to cyanogenesis-associated HCN exposure. We identified all ß-CAS and Aox gene copies present in the genome, including members of the Aox1, Aox2a, and Aox2d subfamilies previously reported in legumes. Expression levels were compared between cyanogenic and acyanogenic genotypes and between damaged and undamaged leaf tissue. RESULTS: ß-CAS and Aox2a expression was differentially elevated in cyanogenic genotypes, and tissue damage was not required to induce this increased expression. Aox2d, in contrast, appeared to be upregulated as a generalized wounding response. CONCLUSIONS: These findings suggest a heightened constitutive role for HCN detoxification (via elevated ß-CAS expression) and HCN-toxicity mitigation (via elevated Aox2a expression) in plants that are capable of cyanogenesis. As such, freezing-induced cyanide autotoxicity is unlikely to be the primary selective factor in the evolution of climate-associated cyanogenesis clines.

Dual-species origin of an adaptive chemical defense polymorphism.

Allopolyploid speciation and chemical defense diversification are two of the most characteristic features of plant evolution; although the former has likely shaped the latter, this has rarely been documented. Here we document allopolyploidy-mediated chemical defense evolution in the origin of cyanogenesis (HCN release upon tissue damage) in white clover (Trifolium repens). We combined linkage mapping of the loci that control cyanogenesis (Ac, controlling production of cyanogenic glucosides; and Li, controlling production of their hydrolyzing enzyme linamarase) with genome sequence comparisons between white clover, a recently evolved allotetraploid, and its diploid progenitors (Trifolium pallescens, Trifolium occidentale). The Ac locus (a three-gene cluster comprising the cyanogenic glucoside pathway) is derived from T. occidentale; it maps to linkage group 2O (occidentale subgenome) and is orthologous to a highly similar cluster in the T. occidentale reference genome. By contrast, Li maps to linkage group 4P (pallescens subgenome), indicating an origin in the other progenitor species. These results indicate that cyanogenesis evolved in white clover as a product of the interspecific hybridization that created the species. This allopolyploidization-derived chemical defense, together with subsequent selection on intraspecific cyanogenesis variation, appears to have contributed to white clover's ecological success as a globally distributed weed species.

Micro- and macroevolutionary adaptation through repeated loss of a complete metabolic pathway.

There is growing evidence for the convergent evolution of physically linked gene clusters encoding chemical defense pathways. Metabolic clusters are proposed to evolve because they ensure co-inheritance of all required genes where the defense is favored, and prevent inheritance of toxic partial pathways where it is not. This hypothesis rests on the assumption that clusters evolve in species where selection favors intraspecific polymorphism for the defense; however, they have not been examined in polymorphic species. We examined metabolic cluster evolution in relation to an adaptive polymorphism for cyanogenic glucoside (CNglc) production in clover. Using 163 accessions, we performed CNglc assays, BAC sequencing, Southern hybridizations and molecular evolutionary analyses. We find that the CNglc pathway forms a 138-kb cluster in white clover, and that the adaptive polymorphism occurs through presence/absence of the complete cluster. Component genes are orthologous to those in the distantly related legume Lotus japonicus. These findings provide empirical support for the co-inheritance hypothesis, and they indicate that adaptive CNglc variation in white clover evolves through recurrent deletions of the entire pathway. They further indicate that the shared ancestor of many important legume crops was likely cyanogenic and that this defense was lost repeatedly over the last 50 Myr.

Long-term balancing selection at the Phosphorus Starvation Tolerance 1 (PSTOL1) locus in wild, domesticated and weedy rice (Oryza).

BACKGROUND: The ability to grow in phosphorus-depleted soils is an important trait for rice cultivation in many world regions, especially in the tropics. The Phosphorus Starvation Tolerance 1 (PSTOL1) gene has been identified as underlying the ability of some cultivated rice varieties to grow under low-phosphorus conditions; however, the gene is absent from other varieties. We assessed PSTOL1 presence/absence in a geographically diverse sample of wild, domesticated and weedy rice and sequenced the gene in samples where it is present. RESULTS: We find that the presence/absence polymorphism spans cultivated, weedy and wild Asian rice groups. For the subset of samples that carry PSTOL1, haplotype sequences suggest long-term selective maintenance of functional alleles, but with repeated evolution of loss-of-function alleles through premature stops and frameshift mutations. The loss-of-function alleles have evolved convergently in multiple rice species and cultivated rice varieties. Greenhouse assessments of plant growth under low- and high-phosphorus conditions did not reveal significant associations with PSTOL1 genotype variation; however, the striking signature of balancing selection at this locus suggests that further phenotypic characterizations of PSTOL1 allelic variants is warranted and may be useful for crop improvement. CONCLUSIONS: These findings suggest balancing selection for both functional and non-functional PSTOL1 alleles that predates and transcends Asian rice domestication, a pattern that may reflect fitness tradeoffs associated with geographical variation in soil phosphorus content.

Recurrent gene deletions and the evolution of adaptive cyanogenesis polymorphisms in white clover (Trifolium repens L.).

Understanding the molecular evolution of genes that underlie intraspecific polymorphisms can provide insights into the process of adaptive evolution. For adaptive polymorphisms characterized by gene presence/absence (P/A) variation, underlying loci commonly show signatures of long-term balancing selection, with gene-presence and gene-absence alleles maintained as two divergent lineages. We examined the molecular evolution of two unlinked P/A polymorphisms that underlie a well-documented adaptive polymorphism for cyanogenesis (hydrogen cyanide release with tissue damage) in white clover. Both cyanogenic and acyanogenic plants occur in this species, and the ecological forces that maintain this chemical defence polymorphism have been studied for several decades. Using a sample of 65 plants, we investigated the molecular evolution of sequences flanking the two underlying cyanogenesis genes: Ac/ac (controlling the presence/absence of cyanogenic glucosides) and Li/li (controlling the presence/absence of their hydrolysing enzyme, linamarase). A combination of genome walking, PCR assays, DNA sequence analysis and Southern blotting was used to test whether these adaptive P/A polymorphisms show evidence of long-term balancing selection, or whether gene-absence alleles have evolved repeatedly through independent deletion events. For both loci, we detect no signatures of balancing selection in the closest flanking genomic sequences. Instead, we find evidence for variation in the size of the deletions characterizing gene-absence alleles. These observations strongly suggest that both of these polymorphisms have been evolving through recurrent gene deletions over time. We discuss the genetic mechanisms that could account for this surprising pattern and the implications of these findings for mechanisms of rapid adaptive evolution in white clover.

Evidence on the molecular basis of the Ac/ac adaptive cyanogenesis polymorphism in white clover (Trifolium repens L).

White clover is polymorphic for cyanogenesis, with both cyanogenic and acyanogenic plants occurring in nature. This chemical defense polymorphism is one of the longest-studied and best-documented examples of an adaptive polymorphism in plants. It is controlled by two independently segregating genes: Ac/ac controls the presence/absence of cyanogenic glucosides; and Li/li controls the presence/absence of their hydrolyzing enzyme, linamarase. Whereas Li is well characterized at the molecular level, Ac has remained unidentified. Here we report evidence that Ac corresponds to a gene encoding a cytochrome P450 of the CYP79D protein subfamily (CYP79D15), and we describe the apparent molecular basis of the Ac/ac polymorphism. CYP79D orthologs catalyze the first step in cyanogenic glucoside biosynthesis in other cyanogenic plant species. In white clover, Southern hybridizations indicate that CYP79D15 occurs as a single-copy gene in cyanogenic plants but is absent from the genomes of ac plants. Gene-expression analyses by RT-PCR corroborate this finding. This apparent molecular basis of the Ac/ac polymorphism parallels our previous findings for the Li/li polymorphism, which also arises through the presence/absence of a single-copy gene. The nature of these polymorphisms may reflect white clover's evolutionary origin as an allotetraploid derived from cyanogenic and acyanogenic diploid progenitors.

Adaptive gains through repeated gene loss: parallel evolution of cyanogenesis polymorphisms in the genus Trifolium (Fabaceae).

Variation in cyanogenesis (hydrogen cyanide release following tissue damage) was first noted in populations of white clover more than a century ago, and subsequent decades of research have established this system as a classic example of an adaptive chemical defence polymorphism. Here, we document polymorphisms for cyanogenic components in several relatives of white clover, and we determine the molecular basis of this trans-specific adaptive variation. One hundred and thirty-nine plants, representing 13 of the 14 species within Trifolium section Trifoliastrum, plus additional species across the genus, were assayed for cyanogenic components (cyanogenic glucosides and their hydrolysing enzyme, linamarase) and for the presence of underlying cyanogenesis genes (CYP79D15 and Li, respectively). One or both cyanogenic components were detected in seven species, all within section Trifoliastrum; polymorphisms for the presence/absence (PA) of components were detected in six species. In a pattern that parallels our previous findings for white clover, all observed biochemical polymorphisms correspond to gene PA polymorphisms at CYP79D15 and Li. Relationships of DNA sequence haplotypes at the cyanogenesis loci and flanking genomic regions suggest independent evolution of gene deletions within species. This study thus provides evidence for the parallel evolution of adaptive biochemical polymorphisms through recurrent gene deletions in multiple species.