Searching for the bull's eye: agents and targets of selection vary among geographically disparate cyanogenesis clines in white clover (Trifolium repens L.)

Introduced species shed friends as well as enemies

Many studies seeking to understand the success of biological invasions focus on species' escape from negative interactions, such as damage from herbivores, pathogens, or predators in their introduced range (enemy release). However, much less work has been done to assess the possibility that introduced species might shed mutualists such as pollinators, seed dispersers, and mycorrhizae when they are transported to a new range. We ran a cross-continental field study and found that plants were being visited by 2.6 times more potential pollinators with 1.8 times greater richness in their native range than in their introduced range. Understanding both the positive and negative consequences of introduction to a new range can help us predict, monitor, and manage future invasion events.

Haplotype-Resolved, Chromosome-Level Assembly of White Clover (Trifolium repens L., Fabaceae)

White clover (Trifolium repens L.; Fabaceae) is an important forage and cover crop in agricultural pastures around the world and is increasingly used in evolutionary ecology and genetics to understand the genetic basis of adaptation. Historically, improvements in white clover breeding practices and assessments of genetic variation in nature have been hampered by a lack of high-quality genomic resources for this species, owing in part to its high heterozygosity and allotetraploid hybrid origin. Here, we use PacBio HiFi and chromosome conformation capture (Omni-C) technologies to generate a chromosome-level, haplotype-resolved genome assembly for white clover totaling 998 Mbp (scaffold N50 = 59.3 Mbp) and 1 Gbp (scaffold N50 = 58.6 Mbp) for haplotypes 1 and 2, respectively, with each haplotype arranged into 16 chromosomes (8 per subgenome). We additionally provide a functionally annotated haploid mapping assembly (968 Mbp, scaffold N50 = 59.9 Mbp), which drastically improves on the existing reference assembly in both contiguity and assembly accuracy. We annotated 78,174 protein-coding genes, resulting in protein BUSCO completeness scores of 99.6% and 99.3% against the embryophyta_odb10 and fabales_odb10 lineage datasets, respectively.

Assessing the Benefits and Costs of the Hydrogen Cyanide Antiherbivore Defense in Trifolium repens

Understanding the evolution of plant defenses against herbivores requires identifying the benefits and costs of defense. Here, we tested the hypothesis that the benefits and costs of hydrogen cyanide (HCN) defense against herbivory on white clover (Trifolium repens) are temperature dependent. We first tested how temperature affected HCN production in vitro, and then examined how temperature influenced the efficacy of HCN defense of T. repens against a generalist slug (Deroceras reticulatum) herbivore using no-choice and choice feeding trial assays. To understand how temperature affected the costs of defense, plants were exposed to freezing, and HCN production, photosynthetic activity, and ATP concentration were quantified. HCN production increased linearly from 5 °C to 50 °C, and cyanogenic plants experienced reduced herbivory compared to acyanogenic plants only at warmer temperatures when fed upon by young slugs. Freezing temperatures induced cyanogenesis in T. repens and decreased chlorophyll fluorescence. Cyanogenic plants experienced lower ATP levels than acyanogenic plants due to freezing. Our study provides evidence that the benefits of HCN defense against herbivores are temperature dependent, and freezing may inhibit ATP production in cyanogenic plants, but the physiological performance of all plants recovered quickly following short-term freezing. These results contribute to understanding how varying environments alter the benefits and costs of defense in a model system for the study of plant chemical defenses against herbivores.

Outlier analyses and genome-wide association study identify glgC and ERD6-like 4 as candidate genes for foliar water-soluble carbohydrate accumulation in Trifolium repens

Increasing water-soluble carbohydrate (WSC) content in white clover is important for improving nutritional quality and reducing environmental impacts from pastoral agriculture. Elucidation of genes responsible for foliar WSC variation would enhance genetic improvement by enabling molecular breeding approaches. The aim of the present study was to identify single nucleotide polymorphisms (SNPs) associated with variation in foliar WSC in white clover. A set of 935 white clover individuals, randomly sampled from five breeding pools selectively bred for divergent (low or high) WSC content, were assessed with 14,743 genotyping-by-sequencing SNPs, using three outlier detection methods: PCAdapt, BayeScan and KGD-F_ST. These analyses identified 33 SNPs as discriminating between high and low WSC populations and putatively under selection. One SNP was located in the intron of ERD6-like 4, a gene coding for a sugar transporter located on the vacuole membrane. A genome-wide association study using a subset of 605 white clover individuals and 5,757 SNPs, identified a further 12 SNPs, one of which was associated with a starch biosynthesis gene, glucose-1-phosphate adenylyltransferase, glgC. Our results provide insight into genomic regions underlying WSC accumulation in white clover, identify candidate genomic regions for further functional validation studies, and reveal valuable information for marker-assisted or genomic selection in white clover.

Evolution in response to climate in the native and introduced ranges of a globally distributed plant

The extent to which species can adapt to spatiotemporal climatic variation in their native and introduced ranges remains unresolved. To address this, we examined how clines in cyanogenesis (hydrogen cyanide [HCN] production-an antiherbivore defense associated with decreased tolerance to freezing) have shifted in response to climatic variation in space and time over a 60-year period in both the native and introduced ranges of Trifolium repens. HCN production is a polymorphic trait controlled by variation at two Mendelian loci (Ac and Li). Using phenotypic assays, we estimated within-population frequencies of HCN production and dominant alleles at both loci (i.e., Ac and Li) from 10,575 plants sampled from 131 populations on five continents, and then compared these frequencies to those from historical data collected in the 1950s. There were no clear relationships between changes in the frequency of HCN production, Ac, or Li and changes in temperature between contemporary and historical samples. We did detect evidence of continued evolution to temperature gradients in the introduced range, whereby the slope of contemporary clines for HCN and Ac in relation to winter temperature became steeper than historical clines and more similar to native clines. These results suggest that cyanogenesis clines show no clear changes through time in response to global warming, but introduced populations continue to adapt to their contemporary environments.

Global urban environmental change drives adaptation in white clover

Urbanization transforms environments in ways that alter biological evolution. We examined whether urban environmental change drives parallel evolution by sampling 110,019 white clover plants from 6169 populations in 160 cities globally. Plants were assayed for a Mendelian antiherbivore defense that also affects tolerance to abiotic stressors. Urban-rural gradients were associated with the evolution of clines in defense in 47% of cities throughout the world. Variation in the strength of clines was explained by environmental changes in drought stress and vegetation cover that varied among cities. Sequencing 2074 genomes from 26 cities revealed that the evolution of urban-rural clines was best explained by adaptive evolution, but the degree of parallel adaptation varied among cities. Our results demonstrate that urbanization leads to adaptation at a global scale.

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.

Genetic trade-offs underlie divergent life history strategies for local adaptation in white clover

Local adaptation is common in plants, yet characterization of its underlying genetic basis is rare in herbaceous perennials. Moreover, while many plant species exhibit intraspecific chemical defence polymorphisms, their importance for local adaptation remains poorly understood. We examined the genetic architecture of local adaptation in a perennial, obligately-outcrossing herbaceous legume, white clover (Trifolium repens). This widespread species displays a well-studied chemical defence polymorphism for cyanogenesis (HCN release following tissue damage) and has evolved climate-associated cyanogenesis clines throughout its range. Two biparental F₂  mapping populations, derived from three parents collected in environments spanning the U.S. latitudinal species range (Duluth, MN, St. Louis, MO and Gainesville, FL), were grown in triplicate for two years in reciprocal common garden experiments in the parental environments (6,012 total plants). Vegetative growth and reproductive fitness traits displayed trade-offs across reciprocal environments, indicating local adaptation. Genetic mapping of fitness traits revealed a genetic architecture characterized by allelic trade-offs between environments, with 100% and 80% of fitness QTL in the two mapping populations showing significant QTL×E interactions, consistent with antagonistic pleiotropy. Across the genome there were three hotspots of QTL colocalization. Unexpectedly, we found little evidence that the cyanogenesis polymorphism contributes to local adaptation. Instead, divergent life history strategies in reciprocal environments were major fitness determinants: selection favoured early investment in flowering at the cost of multiyear survival in the southernmost site versus delayed flowering and multiyear persistence in the northern environments. Our findings demonstrate that multilocus genetic trade-offs contribute to contrasting life history characteristics that allow for local adaptation in this outcrossing herbaceous perennial.

Plantago spp. as Models for Studying the Ecology and Evolution of Species Interactions across Environmental Gradients

AbstractA central challenge in ecology and evolutionary biology is to understand how variation in abiotic and biotic factors combine to shape the distribution, abundance, and diversity of focal species. Environmental gradients, whether natural (e.g., latitude, elevation, ocean proximity) or anthropogenic (e.g., land-use intensity, urbanization), provide compelling settings for addressing this challenge. However, not all organisms are amenable to the observational and experimental approaches required for untangling the factors that structure species along gradients. Here we highlight herbaceous plants in the genus Plantago as models for studying the ecology and evolution of species interactions along abiotic gradients. Plantago lanceolata and P. major are native to Europe and Asia but distributed globally, and they are established models for studying population ecology and interactions with herbivores, pathogens, and soil microbes. Studying restricted range congeners in comparison with those cosmopolitan species can provide insight into abiotic and biotic determinants of range size and population structure. We highlight one such species, P. rugelii, which is endemic to eastern North America. We give an overview of the literature on these focal Plantago species and explain why they are logical candidates for studies of species interactions across environmental gradients. Finally, we emphasize collaborative and community science approaches that can facilitate such research and note the amenability of Plantago for authentic research projects in science education.

Predicting the strength of urban-rural clines in a Mendelian polymorphism along a latitudinal gradient

Cities are emerging as models for addressing the fundamental question of whether populations evolve in parallel to similar environments. Here, we examine the environmental factors that drive the evolution of parallel urban‐rural clines in a Mendelian trait—the cyanogenic antiherbivore defense of white clover (Trifolium repens). Previous work suggested urban‐rural gradients in frost and snow depth could drive the evolution of reduced hydrogen cyanide (HCN) frequencies in urban populations. Here, we sampled over 700 urban and rural clover populations across 16 cities along a latitudinal transect in eastern North America. In each population, we quantified changes in the frequency of genotypes that produce HCN, and in a subset of the cities we estimated the frequency of the alleles at the two genes (CYP79D15 and Li) that epistatically interact to produce HCN. We then tested the hypothesis that cold climatic conditions are necessary for the evolution of cyanogenesis clines by comparing the strength of clines among cities located along a latitudinal gradient of winter temperature and frost exposure. Overall, half of the cities exhibited urban‐rural clines in the frequency of HCN, whereby urban populations evolved lower HCN frequencies. Clines did not evolve in cities with the lowest temperatures and greatest snowfall, supporting the hypothesis that snow buffers plants against winter frost and constrains the formation of clines. By contrast, the strongest clines occurred in the warmest cities where snow and frost are rare, suggesting that alternative selective agents are maintaining clines in warmer cities. Some clines were driven by evolution at only CYP79D15, consistent with stronger and more consistent selection on this locus than on Li. Together, our results demonstrate that urban environments often select for similar phenotypes, but different selective agents and targets underlie the evolutionary response in different cities.

Modern spandrels: the roles of genetic drift, gene flow and natural selection in the evolution of parallel clines

Urban environments offer the opportunity to study the role of adaptive and non-adaptive evolutionary processes on an unprecedented scale. While the presence of parallel clines in heritable phenotypic traits is often considered strong evidence for the role of natural selection, non-adaptive evolutionary processes can also generate clines, and this may be more likely when traits have a non-additive genetic basis due to epistasis. In this paper, we use spatially explicit simulations modelled according to the cyanogenesis (hydrogen cyanide, HCN) polymorphism in white clover (Trifolium repens) to examine the formation of phenotypic clines along urbanization gradients under varying levels of drift, gene flow and selection. HCN results from an epistatic interaction between two Mendelian-inherited loci. Our results demonstrate that the genetic architecture of this trait makes natural populations susceptible to decreases in HCN frequencies via drift. Gradients in the strength of drift across a landscape resulted in phenotypic clines with lower frequencies of HCN in strongly drifting populations, giving the misleading appearance of deterministic adaptive changes in the phenotype. Studies of heritable phenotypic change in urban populations should generate null models of phenotypic evolution based on the genetic architecture underlying focal traits prior to invoking selection's role in generating adaptive differentiation.

Getting ahead of the curve: cities as surrogates for global change

Urbanization represents an unintentional global experiment that can provide insights into how species will respond and interact under future global change scenarios. Cities produce many conditions that are predicted to occur widely in the future, such as warmer temperatures, higher carbon dioxide (CO2) concentrations and exacerbated droughts. In using cities as surrogates for global change, it is challenging to disentangle climate variables-such as temperature-from co-occurring or confounding urban variables-such as impervious surface-and then to understand the interactive effects of multiple climate variables on both individual species and species interactions. However, such interactions are also difficult to replicate experimentally, and thus the challenges of cities are also their unique advantage. Here, we review insights gained from cities, with a focus on plants and arthropods, and how urban findings agree or disagree with experimental predictions and historical data. We discuss the types of hypotheses that can be best tested in cities, caveats to urban research and how to further validate cities as surrogates for global change. Lastly, we summarize how to achieve the goal of using urban species responses to predict broader regional- and ecosystem-level patterns in the future.

Freeze-induced cyanide toxicity does not maintain the cyanogenesis polymorphism in white clover (Trifolium repens)

PREMISE OF THE STUDY

The maintenance of adaptive polymorphisms within species requires fitness trade-offs reflecting selection for each morph. Cyanogenesis, the ability to produce hydrogen cyanide (HCN) after tissue damage, occurs in >3000 plant species and exists as a discrete polymorphism in white clover. This polymorphism is spatially distributed in recurrent clines, with higher frequencies of cyanogenic plants in warmer climates. The HCN autotoxicity hypothesis proposes that cyanogenic plants are selected against where frosts are common, as freezing liberates HCN and could impair cellular respiration.

METHODS

We tested the HCN autotoxicity hypothesis using a freezing chamber to examine survival, tissue damage, and physiological recovery as assessed via chlorophyll fluorescence following mild and severe freezing treatments. We utilized 65 genotypes from a single polymorphic population to eliminate effects of population structure.

KEY RESULTS

Cyanogenic plants did not differ from acyanogenic plants in survival, tissue damage, or recovery following freezing. However, plants producing either of the two required cyanogenic precursors had lower survival and tissue damage after freezing than plants lacking both precursors.

CONCLUSIONS

These results suggest that freezing-induced HCN toxicity is unlikely to be responsible for the maintenance of the cyanogenesis polymorphism in white clover. However, energetic trade-offs associated with costs of producing the cyanogenic precursors may confer a fitness benefit to acyanogenic plants under stressful climatic conditions. The lack of evidence for HCN toxicity suggests that cyanogenic clover uses physiological mechanisms mediated by β-cyanoalanine synthase and alternative oxidase to maintain cellular function in the presence of HCN.

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.

Continent-Wide Climatic Variation Drives Local Adaptation in North American White Clover

Climate-associated clines in adaptive polymorphisms are commonly cited as evidence of local adaptation within species. However, the contribution of the clinally varying trait to overall fitness is often unknown. To address this question, we examined survival, vegetative growth, and reproductive output in a central US common garden experiment using 161 genotypes of white clover (Trifolium repens L.) originating from 15 locations across North America. White clover is polymorphic for cyanogenesis (hydrogen cyanide release upon tissue damage), a chemical defense against generalist herbivores, and climate-associated cyanogenesis clines have repeatedly evolved across the species range. Over a 12-month experiment, we observed striking correlations between the population of origin and plant performance in the common garden, with climatic distance from the common garden site predicting fitness more accurately than geographic distance. Assessments of herbivore leaf damage over the 2015 growing season indicated marginally lower herbivory on cyanogenic plants; however, this effect did not result in increased fitness in the common garden location. Linear mixed modeling suggested that while cyanogenesis variation had little predictive value for vegetative growth, it is as important as climatic variation for predicting reproductive output in the central United States. Together, our findings suggest that knowledge of climate similarity, as well as knowledge of locally favored adaptive traits, will help to inform transplantation strategies for restoration ecology and other conservation efforts in the face of climate change.

Urbanization drives the evolution of parallel clines in plant populations

Urban ecosystems are an increasingly dominant feature of terrestrial landscapes. While evidence that species can adapt to urban environments is accumulating, the mechanisms through which urbanization imposes natural selection on populations are poorly understood. The identification of adaptive phenotypic changes (i.e. clines) along urbanization gradients would facilitate our understanding of the selective factors driving adaptation in cities. Here, we test for phenotypic clines in urban ecosystems by sampling the frequency of a Mendelian-inherited trait-cyanogenesis-in white clover (Trifolium repens L.) populations along urbanization gradients in four cities. Cyanogenesis protects plants from herbivores, but reduces tolerance to freezing temperatures. We found that the frequency of cyanogenic plants within populations decreased towards the urban centre in three of four cities. A field experiment indicated that spatial variation in herbivory is unlikely to explain these clines. Rather, colder minimum winter ground temperatures in urban areas compared with non-urban areas, caused by reduced snow cover in cities, may select against cyanogenesis. In the city with no cline, high snow cover might protect plants from freezing damage in the city centre. Our study suggests that populations are adapting to urbanization gradients, but regional climatic patterns may ultimately determine whether adaptation occurs.

Replicate altitudinal clines reveal that evolutionary flexibility underlies adaptation to drought stress in annual Mimulus guttatus

Examining how morphology, life history and physiology vary along environmental clines can reveal functional insight into adaptations to climate and thus inform predictions about evolutionary responses to global change. Widespread species occurring over latitudinal and altitudinal gradients in seasonal water availability are excellent systems for investigating multivariate adaptation to drought stress. Under common garden conditions, we characterized variation in 27 traits for 52 annual populations of Mimulus guttatus sampled from 10 altitudinal transects. We also assessed variation in the critical photoperiod for flowering and surveyed neutral genetic markers to control for demography when analyzing clinal patterns. Many drought escape (e.g. flowering time) and drought avoidance (e.g. specific leaf area, succulence) traits exhibited geographic or climatic clines, which often remained significant after accounting for population structure. Critical photoperiod and flowering time in glasshouse conditions followed distinct clinal patterns, indicating different aspects of seasonal phenology confer adaptation to unique agents of selection. Although escape and avoidance traits were negatively correlated range-wide, populations from sites with short growing seasons produced both early flowering and dehydration avoidance phenotypes. Our results highlight how abundant genetic variation in the component traits that build multivariate adaptations to drought stress provides flexibility for intraspecific adaptation to diverse climates.

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.

Aridity shapes cyanogenesis cline evolution in white clover (Trifolium repens L.)

Adaptive differentiation between populations is often proposed to be the product of multiple interacting selective pressures, although empirical support for this is scarce. In white clover, populations show adaptive differentiation in frequencies of cyanogenesis, the ability to produce hydrogen cyanide after tissue damage. This polymorphism arises through independently segregating polymorphisms for the presence/absence of two required cyanogenic components, cyanogenic glucosides and their hydrolysing enzyme. White clover populations worldwide have evolved a series of recurrent, climate-associated clines, with higher frequencies of cyanogenic plants in warmer locations. These clines have traditionally been hypothesized to reflect a fitness trade-off between chemical defence in herbivore-rich areas (warmer climates) and energetic costs of producing cyanogenic components in areas of low herbivore pressure (cooler climates). Recent observational studies suggest that cyanogenic components may also be beneficial in water-stressed environments. We investigated fitness trade-offs associated with temperature-induced water stress in the cyanogenesis system using manipulative experiments in growth chambers and population surveys across a longitudinal precipitation gradient in the central United States. We find that plants producing cyanogenic glucosides have higher relative fitness in treatments simulating a moderate, persistent drought stress. In water-neutral treatments, there are energetic costs to producing cyanogenic components, but only in treatments with nutrient stress. These fitness trade-offs are consistent with cyanogenesis frequencies in natural populations, where we find clinal variation in the proportion of plants producing cyanogenic glucosides along the precipitation gradient. These results suggest that multiple selective pressures interact to maintain this adaptive polymorphism and that modelling adaptation will require knowledge of environment-specific fitness effects.