Rapid evolution of an adaptive cyanogenesis cline in introduced North American white clover (Trifolium repens L.)

White clover is polymorphic for cyanogenesis (HCN production after tissue damage), and this herbivore defence polymorphism has served as a classic model for studying adaptive variation. The cyanogenic phenotype requires two interacting biochemical components; the presence/absence of each component is controlled by a simple Mendelian gene (Ac/ac and Li/li). Climate-associated cyanogenesis clines occur in both native (Eurasian) and introduced populations worldwide, with cyanogenic plants predominating in warmer locations. Moreover, previous studies have suggested that epistatic selection may act within populations to maintain cyanogenic (AcLi) plants and acyanogenic plants that lack both components (acli plants) at the expense of plants possessing a single component (Acli and acLi plants). Here, we examine the roles of selection, gene flow and demography in the evolution of a latitudinal cyanogenesis cline in introduced North American populations. Using 1145 plants sampled across a 1650 km transect, we determine the distribution of cyanogenesis variation across the central United States and investigate whether clinal variation is adaptive or an artefact of population introduction history. We also test for the evidence of epistatic selection. We detect a clear latitudinal cline, with cyanogenesis frequencies increasing from 11% to 86% across the transect. Population structure analysis using nine microsatellite loci indicates that the cline is adaptive and not a by-product of demographic history. However, we find no evidence for epistatic selection within populations. Our results provide strong evidence for rapid adaptive evolution in these introduced populations, and they further suggest that the mechanisms maintaining adaptive variation may vary among populations of a species.     

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

The recurrent evolution of adaptive clines within a species can be used to elucidate the selective factors and genetic responses that underlie adaptation. White clover is polymorphic for cyanogenesis (HCN release with tissue damage), and climate-associated cyanogenesis clines have evolved throughout the native and introduced species range. This polymorphism arises through two independently segregating Mendelian polymorphisms for the presence/absence of two required components: cyanogenic glucosides and their hydrolyzing enzyme linamarase. Cyanogenesis is commonly thought to function in herbivore defense; however, the individual cyanogenic components may also serve other physiological functions. To test whether cyanogenesis clines have evolved in response to the same selective pressures acting on the same genetic targets, we examined cyanogenesis cline shape and its environmental correlates in three world regions: southern New Zealand, the central United States and the US Pacific Northwest. For some regional comparisons, cline shapes are remarkably similar despite large differences in the spatial scales over which clines occur (40–1600 km). However, we also find evidence for major differences in both the agents and targets of selection among the sampled clines. Variation in cyanogenesis frequency is best predicted using a combination of minimum winter temperature and aridity variables. Together, our results provide evidence that recurrent adaptive clines do not necessarily reflect shared adaptive processes.     

Adaptive cyanogenesis clines evolve recurrently through geographical sorting of existing gene deletions

Identifying the genetic basis of parallel phenotypic evolution provides insight into the process of adaptation and evolutionary constraint. White clover (Trifolium repens) has evolved climate‐associated adaptive clines in cyanogenesis (the ability to produce hydrogen cyanide upon tissue damage) in several world regions where it has been introduced. Gene‐deletion polymorphisms at the CYP79D15 and Li loci underlie the presence/absence of the cyanogenic phenotype. Both loci have undergone multiple independent gene‐deletion events, which are identifiable through molecular signatures in flanking regions. To investigate whether cyanogenesis clines in introduced populations have evolved through the sorting of standing genetic variation or de novo gene deletions, we examined cyanogenesis gene‐flanking regions in three world regions. In comparison with native Eurasian populations, we find no evidence for novel gene deletion events in any introduced region, which suggests that these adaptive clines have evolved through the geographical sorting of pre‐existing genetic variation.     

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.     

Impact of polymorphic transposable elements on linkage disequilibrium along chromosomes

Recombination and selection drive the extent of linkage disequilibrium (LD) among loci and therefore affect the reshuffling of adaptive genetic variation. However, it is poorly known to what extent the enrichment of transposable elements (TEs) in recombinationally‐inert regions reflects their inefficient removal by purifying selection and whether the presence of polymorphic TEs can modify the local recombination rate. In this study, we investigate how TEs and recombination interact at fine scale along chromosomes and possibly support linked selection in natural populations. Whole‐genome sequencing data of 304 individuals from nearby alpine populations of Arabis alpina were used to show that the density of polymorphic TEs is specifically correlated with local LD along chromosomes. Consistent with TEs modifying recombination, the characterization of 28 such LD blocks of up to 5.5 Mb in length revealed strong evidence of selective sweeps at a few loci through either site frequency spectrum or haplotype structure. A majority of these blocks were enriched in genes related to ecologically relevant functions such as responses to cold, salt stress or photoperiodism. In particular, the S‐locus (i.e., supergene responsible for strict outcrossing) was identified in a LD block with high levels of polymorphic TEs and evidence of selection. Another such LD block was enriched in cold‐responding genes and presented evidence of adaptive loci related to photoperiodism and flowering being increasingly linked by polymorphic TEs. These results are consistent with the hypothesis that TEs modify recombination landscapes and thus interact with selection in driving blocks of linked adaptive loci in natural populations.     

Extensive linkage disequilibrium and parallel adaptive divergence across threespine stickleback genomes

Population genomic studies are beginning to provide a more comprehensive view of dynamic genome-scale processes in evolution. Patterns of genomic architecture, such as genomic islands of increased divergence, may be important for adaptive population differentiation and speciation. We used next-generation sequencing data to examine the patterns of local and long-distance linkage disequilibrium (LD) across oceanic and freshwater populations of threespine stickleback, a useful model for studies of evolution and speciation. We looked for associations between LD and signatures of divergent selection, and assessed the role of recombination rate variation in generating LD patterns. As predicted under the traditional biogeographic model of unidirectional gene flow from ancestral oceanic to derived freshwater stickleback populations, we found extensive local and long-distance LD in fresh water. Surprisingly, oceanic populations showed similar patterns of elevated LD, notably between large genomic regions previously implicated in adaptation to fresh water. These results support an alternative biogeographic model for the stickleback radiation, one of a metapopulation with appreciable bi-directional gene flow combined with strong divergent selection between oceanic and freshwater populations. As predicted by theory, these processes can maintain LD within and among genomic islands of divergence. These findings suggest that the genomic architecture in oceanic stickleback populations may provide a mechanism for the rapid re-assembly and evolution of multi-locus genotypes in newly colonized freshwater habitats, and may help explain genetic mapping of parallel phenotypic variation to similar loci across independent freshwater populations.     

Widespread phenotypic and genetic divergence along altitudinal gradients in animals

Altitudinal gradients offer valuable study systems to investigate how adaptive genetic diversity is distributed within and between natural populations and which factors promote or prevent adaptive differentiation. The environmental clines along altitudinal gradients tend to be steep relative to the dispersal distance of many organisms, providing an opportunity to study the joint effects of divergent natural selection and gene flow. Temperature is one variable showing consistent altitudinal changes, and altitudinal gradients can therefore provide spatial surrogates for some of the changes anticipated under climate change. Here, we investigate the extent and patterns of adaptive divergence in animal populations along altitudinal gradients by surveying the literature for (i) studies on phenotypic variation assessed under common garden or reciprocal transplant designs and (ii) studies looking for signatures of divergent selection at the molecular level. Phenotypic data show that significant between‐population differences are common and taxonomically widespread, involving traits such as mass, wing size, tolerance to thermal extremes and melanization. Several lines of evidence suggest that some of the observed differences are adaptively relevant, but rigorous tests of local adaptation or the link between specific phenotypes and fitness are sorely lacking. Evidence for a role of altitudinal adaptation also exists for a number of candidate genes, most prominently haemoglobin, and for anonymous molecular markers. Novel genomic approaches may provide valuable tools for studying adaptive diversity, also in species that are not amenable to experimentation.     

Population structure attributable to reproductive time: isolation by time and adaptation by time

Many populations are composed of a mixture of individuals that reproduce at different times, and these times are often heritable. Under these conditions, gene flow should be limited between early and late reproducers, even within populations having a unimodal temporal distribution of reproductive activity. This temporal restriction on gene flow might be called "isolation by time" (IBT) to acknowledge its analogy with isolation by distance (IBD). IBD and IBT are not exactly equivalent, however, owing to differences between dispersal in space and dispersal in time. We review empirical studies of natural populations that provide evidence for IBT based on heritabilities of reproductive time and on molecular genetic differences associated with reproductive time. When IBT is present, variation in selection through the reproductive season may lead to adaptive temporal variation in phenotypic traits [adaptation by time (ABT)]. We introduce a novel theoretical model that shows how ABT increases as (i) selection on the trait increases; (ii) environmental influences on reproductive time decrease; (iii) the heritability of reproductive time increases; and (iv) the temporal distribution of reproductive activity becomes increasingly uniform. We then review empirical studies of natural populations that provide evidence for ABT by documenting adaptive temporal clines in phenotypic traits. The best evidence for IBT and ABT currently comes from salmonid fishes and flowering plants, but we expect that future work will show these processes are more widespread.     

Latitudinal variation for two enzyme loci and an inversion polymorphism in Drosophila melanogaster from Central and South America

Abstract.— Many organisms show latitudinal variation for various genetically determined traits. Such clines may involve neutral variation and originate from historical events or their maintenance may be explained by selection. For Drosophila melanogaster , latitudinal variation for allozymes, inversions, and quantitative traits has been found on several continents. We sampled D. melanogaster populations in Panama and along a transect of 40 latitudinal degrees on the west coast of South America. Negative correlations with latitude were found for Adh^s and _αGpdh^F allele frequencies and for the frequency of the cosmopolitan inversion In(2L)t in Adh^sαpdh^F chromosomes. A positive correlation existed between wing length and latitude. Significant correlations were found between these traits and climatic variables like temperature and rainfall. The observed clines show considerable resemblance to those found on other continents. Gametic disequilibrium between Adh^s and αGpdh^F occurred predominantly at higher latitudes and was caused by the presence of In(2L)t . The reasons for the clinal distributions are discussed and it is argued that selection is the most likely explanation. However, the exact nature of the selective force and the interactions of allozymes with each other and with In(2L)t are complex and not fully understood. In tropical regions In(2L)t -containing genotypes have higher fitness than ST/ST and Adh and αGpdh hitchhike with the inversion, but there is also evidence for balancing selection at the Adh locus.     

The adaptive significance of chromosomal inversion polymorphisms in Drosophila melanogaster

Chromosomal inversions, structural mutations that reverse a segment of a chromosome, cause suppression of recombination in the heterozygous state. Several studies have shown that inversion polymorphisms can form clines or fluctuate predictably in frequency over seasonal time spans. These observations prompted the hypothesis that chromosomal rearrangements might be subject to spatially and/or temporally varying selection. Here, we review what has been learned about the adaptive significance of inversion polymorphisms in the vinegar fly Drosophila melanogaster, the species in which they were first discovered by Sturtevant in 1917. A large body of work provides compelling evidence that several inversions in this system are adaptive; however, the precise selective mechanisms that maintain them polymorphic in natural populations remain poorly understood. Recent advances in population genomics, modelling and functional genetics promise to greatly improve our understanding of this long‐standing and fundamental problem in the near future.     

Disentangling the roles of history and local selection in shaping clinal variation of allele frequencies and gene expression in Norway spruce (Picea abies)

Understanding the genetic basis of local adaptation is challenging due to the subtle balance among conflicting evolutionary forces that are involved in its establishment and maintenance. One system with which to tease apart these difficulties is clines in adaptive characters. Here we analyzed genetic and phenotypic variation in bud set, a highly heritable and adaptive trait, among 18 populations of Norway spruce (Picea abies), arrayed along a latitudinal gradient ranging from 47°N to 68°N. We confirmed that variation in bud set is strongly clinal, using a subset of five populations. Genotypes for 137 single-nucleotide polymorphisms (SNPs) chosen from 18 candidate genes putatively affecting bud set and 308 control SNPs chosen from 264 random genes were analyzed for patterns of genetic structure and correlation to environment. Population genetic structure was low (F(ST) = 0.05), but latitudinal patterns were apparent among Scandinavian populations. Hence, part of the observed clinal variation should be attributable to population demography. Conditional on patterns of genetic structure, there was enrichment of SNPs within candidate genes for correlations with latitude. Twenty-nine SNPs were also outliers with respect to F(ST). The enrichment for clinal variation at SNPs within candidate genes (i.e., SNPs in PaGI, PaPhyP, PaPhyN, PaPRR7, and PaFTL2) indicated that local selection in the 18 populations, and/or selection in the ancestral populations from which they were recently derived, shaped the observed cline. Validation of these genes using expression studies also revealed that PaFTL2 expression is significantly associated with latitude, thereby confirming the central role played by this gene in the control of phenology in plants.     

HYBRID ZONES WITH DOBZHANSKY-TYPE EPISTATIC SELECTION

Dobzhansky's model of epistatic selection assumes that viable genotypes form “clusters” in genotype space so that populations can evolve from one state to a reproductively isolated state following a “ridge” of well-fit genotypes without crossing any deep adaptive valleys. Recently, the importance of Dobzhansky-type models in evolutionary studies has been reemphasized by Gavrilets (1997a) and Gavrilets and Gravner (1997) who argue that the existence of “ridges” of well-fit genotypes connecting reproductively isolated genotypes is actually a general property of multidimensional adaptive landscapes. Using rigorous techniques and numerical simulations, I analyze clines in the frequencies of selected and neutral alleles maintained by a balance of migration and Dobzhansky-type epistatic selection acting on two diallelic loci. I show that Dobzhansky-type epistatic selection can build up a very strong barrier to neutral gene flow. I describe properties of clines that are indicative of Dobzhansky-type selection.     

Adaptive potential of northernmost tree populations to climate change, with emphasis on Scots pine (Pinus sylvestris L.)

The adaptive potential of the northernmost Pinus sylvestris L. (and other northern tree) populations is considered by examining first the current patterns of quantitative genetic adaptive traits, which show high population differentiation and clines. We then consider the postglacial history of the populations using both paleobiological and genetic data. The current patterns of diversity at nuclear genes suggest that the traces of admixture are mostly visible in mitochondrial DNA variation patterns. There is little evidence of increased diversity due to admixture between an eastern and western colonization lineage, but no signal of reduced diversity (due to sequential bottlenecks) either. Quantitative trait variation in the north is not associated with the colonizing lineages. The current clines arose rapidly and may be based on standing genetic variation. The initial phenotypic response of Scots pine in the north is predicted to be increased survival and growth. The genetic responses are examined based on quantitative genetic predictions of sustained selection response and compared with earlier simulation results that have aimed at more ecological realism. The phenotypic responses of increased growth and survival reduce the opportunity for selection and delay the evolutionary responses. The lengthening of the thermal growing period also causes selection on the critical photoperiod in the different populations. Future studies should aim at including multiple ecological and genetic factors in evaluating potential responses.     

Developmental stability and adaptive radiation in the Spalax ehrenbergi superspecies in the Near-East

Chromosomal species of the mole rat, Spalax ehrenbergi, in Israel have been shown to display distinct adaptive strategies to increasing aridity. This adaptive radiation appeared to be associated with an increase in allozymic heterozygosity. In the present study, the developmental stability (DS) estimated by fluctuating asymmetry (FA) of dental traits was used to assess the suitability of habitat and the efficiency of adaptation to local environmental conditions among populations and chromosomal species. Although FA levels were highly heterogeneous among populations, they were not found to differ between species. DS of populations appeared, however, to be impaired at higher altitudes and in indurate soils. Since these environmental features were largely covariant, the effect of each one could not be precisely determined. Interestingly, while aridity is considered as the major selective force acting on populations southwards, DS was not altered under arid conditions, suggesting that mole rat populations were adapted to their local conditions of aridity. However, the cline of aridity is matched to several environmental and genetic clines among which are the increasing heterozygosity and recombination rate among species southwards. In studies of natural populations, the potential complementary effects of environmental and genetics on DS have to be considered and hamper the interpretation of habitat suitability expressed by DS in terms of adaptive strategies.     

Müllerian mimicry of a quantitative trait despite contrasting levels of genomic divergence and selection

Hybrid zones, where distinct populations meet and interbreed, give insight into how differences between populations are maintained despite gene flow. Studying clines in genetic loci and adaptive traits across hybrid zones is a powerful method for understanding how selection drives differentiation within a single species, but can also be used to compare parallel divergence in different species responding to a common selective pressure. Here, we study parallel divergence of wing colouration in the butterflies Heliconius erato and H. melpomene, which are distantly related Müllerian mimics which show parallel geographic variation in both discrete variation in pigmentation, and quantitative variation in structural colour. Using geographic cline analysis, we show that clines in these traits are positioned in roughly the same geographic region for both species, which is consistent with direct selection for mimicry. However, the width of the clines varies markedly between species. This difference is explained in part by variation in the strength of selection acting on colour traits within each species, but may also be influenced by differences in the dispersal rate and total strength of selection against hybrids between the species. Genotyping‐by‐sequencing also revealed weaker population structure in H. melpomene, suggesting the hybrid zones may have evolved differently in each species, which may also contribute to the patterns of phenotypic divergence in this system. Overall, we conclude that multiple factors are needed to explain patterns of clinal variation within and between these species, although mimicry has probably played a central role.     

PLANT RESPONSES TO CLIMATE IN THE CAPE FLORISTIC REGION OF SOUTH AFRICA: EVIDENCE FOR ADAPTIVE DIFFERENTIATION IN THE PROTEACEAE

Local adaptation along environmental gradients may drive plant species radiation within the Cape Floristic Region (CFR), yet few studies examine the role of ecologically based divergent selection within CFR clades. In this study, we ask whether populations within the monophyletic white protea clade (Protea section Exsertae, Proteaceae) differ in key functional traits along environmental gradients and whether differences are consistent with local adaptation. Using seven taxa, we measured trait–environment associations and selection gradients across 35 populations of wild adults and their offspring grown in two common gardens. Focal traits were leaf size and shape, specific leaf area (SLA), stomatal density, growth, and photosynthetic rate. Analyses on wild and common garden plants revealed heritable trait differences that were associated with gradients in rainfall seasonality, drought stress, cold stress, and less frequently, soil fertility. Divergent selection between gardens generally matched trait–environment correlations and literature‐based predictions, yet variation in selection regimes among wild populations generally did not. Thus, selection via seedling survival may promote gradient‐wide differences in SLA and leaf area more than does selection via adult fecundity. By focusing on the traits, life stages, and environmental clines that drive divergent selection, our study uniquely demonstrates adaptive differentiation among plant populations in the CFR.     

Historical connectivity,contemporary isolation and local adaptation in a widespread but discontinuously distributed species endemic to Taiwan,Rhododendron oldhamii (Ericaceae)

Elucidation of the evolutionary processes that constrain or facilitate adaptive divergence is a central goal in evolutionary biology, especially in non-model organisms. We tested whether changes in dynamics of gene flow (historical vs contemporary) caused population isolation and examined local adaptation in response to environmental selective forces in fragmented Rhododendron oldhamii populations. Variation in 26 expressed sequence tag-simple sequence repeat loci from 18 populations in Taiwan was investigated by examining patterns of genetic diversity, inbreeding, geographic structure, recent bottlenecks, and historical and contemporary gene flow. Selection associated with environmental variables was also examined. Bayesian clustering analysis revealed four regional population groups of north, central, south and southeast with significant genetic differentiation. Historical bottlenecks beginning 9168–13,092 years ago and ending 1584–3504 years ago were revealed by estimates using approximate Bayesian computation for all four regional samples analyzed. Recent migration within and across geographic regions was limited. However, major dispersal sources were found within geographic regions. Altitudinal clines of allelic frequencies of environmentally associated positively selected outliers were found, indicating adaptive divergence. Our results point to a transition from historical population connectivity toward contemporary population isolation and divergence on a regional scale. Spatial and temporal dispersal differences may have resulted in regional population divergence and local adaptation associated with environmental variables, which may have played roles as selective forces at a regional scale.     

Cellular Basis and Developmental Timing in a Size Cline of Drosophila Melanogaster

We examined 20 Drosophila melanogaster populations collected from a 2600-km north-south transect in Australia. In laboratory culture at constant temperature and standard larval density, a genetic cline in thorax length and wing area was found, with both traits increasing with latitude. The cline in wing area was based on clines in both cell size and cell number, but was primarily determined by changes in cell number. Body size and larval development time were not associated among populations. We discuss our results in the context of selection processes operating in natural and experimental populations.     

Genetic differentiation in life-history traits of introduced and native common ragweed (Ambrosia artemisiifolia) populations

Introduced species represent opportunities to observe evolution over contemporary time scales, and as exotics encounter new environments, adaptive responses can occur, potentially contributing to invasion. Here, we compare 22 native North American populations and 12 introduced European populations of common ragweed (Ambrosia artemisiifolia) in five common gardens (control, herbivory, light stress, nutrient stress and drought). We found evidence for improved growth and reproduction of the introduced populations in most environments, particularly in the light stress. However, under drought conditions, the introduced plants experienced more rapid wilting and mortality than their native counterparts, evidence consistent with a life-history trade-off between rapid growth and drought tolerance. Moreover, we found parallel latitudinal clines in flowering time and correlations between fitness components and the local climate of the source populations in both ranges. Together these data provide evidence for adaptation to local environmental conditions in the native and introduced range of common ragweed.     

Spatiotemporal environmental heterogeneity and the maintenance of the tailspot polymorphism in the variable platyfish (Xiphophorus variatus)

Genetic variation is critical for adaptive evolution. Despite its importance, there is still limited evidence in support of some prominent theoretical models explaining the maintenance of genetic polymorphism within populations. We examined 84 populations of Xiphophorus variatus, a livebearing fish with a genetic polymorphism associated with physiological performance, to test: (1) whether niche differentiation explains broad‐scale maintenance of polymorphism, (2) whether polymorphism is maintained among populations by local adaptation and migration, or (3) whether heterogeneity in explicit environmental variables could be linked to levels of polymorphism within populations. We found no evidence of climatic niche differentiation that could generate or maintain broad geographic variation in polymorphism. Subsequently, hierarchical partitioning of genetic richness and partial mantel tests revealed that 76% of the observed genetic richness was partitioned within populations with no effect of geographic distance on polymorphism. These results strongly suggest a lack of migration‐selection balance in the maintenance of polymorphism, and model selection confirmed a significant relationship between environmental heterogeneity and genetic richness within populations. Few studies have demonstrated such effects at this scale, and additional studies in other taxa should examine the generality of gene‐by‐environment interactions across populations to better understand the dynamics and scale of balancing selection.