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.     

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.     

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.     

Classic clover cline clues

Adaptive clines are striking examples of natural selection in action, yet few have been studied in depth. In this issue of Molecular Ecology, Kooyers & Olsen (2012) introduce modern analyses and thinking towards studies of a classical example of the rapid and repeated evolution of latitudinal and altitudinal clines in cyanogenesis in white clover, Trifolium repens L. Recognizing that adaptive clines represent trade-offs in the selective benefits of traits at different ends of a geographical transect, these researchers focus on whether evidence for selection can be found at regional (coarse) and local (fine) scales. After adjusting for population genetic patterns generated by demographic processes, Kooyers and Olsen provide evidence that the cyanogenesis cline is adaptive across a transect from Louisiana to Wisconsin, USA. Within local populations, divergent selection on coupling dominant and recessive alleles that underlie cyanogenesis is predicted to drive populations to gametic phase disequilibrium (LD), a pattern that has been found in several other studies reviewed by Kooyers and Olsen. The absence of LD within any sampled populations in this study leads the authors to suggest that selective patterns within these clines may be more complex than previously proposed, perhaps even following theoretical predictions of a geographic mosaic.     

Broad Concordance in the Spatial Distribution of Adaptive and Neutral Genetic Variation across an Elevational Gradient in Deer Mice

When species are continuously distributed across environmental gradients, the relative strength of selection and gene flow shape spatial patterns of genetic variation, potentially leading to variable levels of differentiation across loci. Determining whether adaptive genetic variation tends to be structured differently than neutral variation along environmental gradients is an open and important question in evolutionary genetics. We performed exome-wide population genomic analysis on deer mice sampled along an elevational gradient of nearly 4,000 m of vertical relief. Using a combination of selection scans, genotype−environment associations, and geographic cline analyses, we found that a large proportion of the exome has experienced a history of altitude-related selection. Elevational clines for nearly 30% of these putatively adaptive loci were shifted significantly up- or downslope of clines for loci that did not bear similar signatures of selection. Many of these selection targets can be plausibly linked to known phenotypic differences between highland and lowland deer mice, although the vast majority of these candidates have not been reported in other studies of highland taxa. Together, these results suggest new hypotheses about the genetic basis of physiological adaptation to high altitude, and the spatial distribution of adaptive genetic variation along environmental gradients.     

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.     

Shifting latitudinal clines in avian body size correlate with global warming in Australian passerines

Intraspecific latitudinal clines in the body size of terrestrial vertebrates, where members of the same species are larger at higher latitudes, are widely interpreted as evidence for natural selection and adaptation to local climate. These clines are predicted to shift in response to climate change. We used museum specimens to measure changes in the body size of eight passerine bird species from south-eastern Australia over approximately the last 100 years. Four species showed significant decreases in body size (1.8–3.6% of wing length) and a shift in latitudinal cline over that period, and a meta-analysis demonstrated a consistent trend across all eight species. Southern high-latitude populations now display the body sizes typical of more northern populations pre-1950, equivalent to a 7° shift in latitude. Using ptilochronology, we found no evidence that these morphological changes were a plastic response to changes in nutrition, a likely non-genetic mechanism for the pattern observed. Our results demonstrate a generalized response by eight avian species to some major environmental change over the last 100 years or so, probably global warming.     

Frequency of cyanogenesis in tropical rainforests of far north Queensland, Australia

BACKGROUND AND AIMS: Plant cyanogenesis is the release of toxic cyanide from endogenous cyanide-containing compounds, typically cyanogenic glycosides. Despite a large body of phytochemical, taxonomic and ecological work on cyanogenic species, little is known of their frequency in natural plant communities. This study aimed to investigate the frequency of cyanogenesis in Australian tropical rainforests. Secondary aims were to quantify the cyanogenic glycoside content of tissues, to investigate intra-plant and intra-population variation in cyanogenic glycoside concentration and to appraise the potential chemotaxonomic significance of any findings in relation to the distribution of cyanogenesis in related taxa. METHODS: All species in six 200 m(2) plots at each of five sites across lowland, upland and highland tropical rainforest were screened for cyanogenesis using Feigl-Anger indicator papers. The concentrations of cyanogenic glycosides were accurately determined for all cyanogenic individuals. KEY RESULTS: Over 400 species from 87 plant families were screened. Overall, 18 species (4.5 %) were cyanogenic, accounting for 7.3 % of total stem basal area. Cyanogenesis has not previously been reported for 17 of the 18 species, 13 of which are endemic to Australia. Several species belong to plant families or orders in which cyanogenesis has been little reported, if at all (e.g. Elaeocarpaceae, Myrsinaceae, Araliaceae and Lamiaceae). A number of species contained concentrations of cyanogenic glycosides among the highest ever reported for mature leaves-up to 5.2 mg CN g(-1) d. wt, for example, in leaves of Elaeocarpus sericopetalus. There was significant variation in cyanogenic glycoside concentration within individuals; young leaves and reproductive tissues typically had higher cyanogen content. In addition, there was substantial variation in cyanogenic glycoside content within populations of single species. CONCLUSIONS: This study expands the limited knowledge of the frequency of cyanogenesis in natural plant communities, includes novel reports of cyanogenesis among a range of taxa and characterizes patterns in intra-plant and intra-population variation of cyanogensis.     

Cyanogenesis of Wild Lima Bean (Phaseolus lunatus L.) Is an Efficient Direct Defence in Nature

In natural systems plants face a plethora of antagonists and thus have evolved multiple defence strategies. Lima bean (Phaseolus lunatus L.) is a model plant for studies of inducible indirect anti-herbivore defences including the production of volatile organic compounds (VOCs) and extrafloral nectar (EFN). In contrast, studies on direct chemical defence mechanisms as crucial components of lima beans'' defence syndrome under natural conditions are nonexistent. In this study, we focus on the cyanogenic potential (HCNp; concentration of cyanogenic glycosides) as a crucial parameter determining lima beans'' cyanogenesis, i.e. the release of toxic hydrogen cyanide from preformed precursors. Quantitative variability of cyanogenesis in a natural population of wild lima bean in Mexico was significantly correlated with missing leaf area. Since existing correlations do not by necessity mean causal associations, the function of cyanogenesis as efficient plant defence was subsequently analysed in feeding trials. We used natural chrysomelid herbivores and clonal lima beans with known cyanogenic features produced from field-grown mother plants. We show that in addition to extensively investigated indirect defences, cyanogenesis has to be considered as an important direct defensive trait affecting lima beans'' overall defence in nature. Our results indicate the general importance of analysing ‘multiple defence syndromes’ rather than single defence mechanisms in future functional analyses of plant defences.     

Secondary contact zones of closely‐related Erebia butterflies overlap with narrow phenotypic and parasitic clines

Zones of secondary contact between closely related taxa are a common legacy of the Quaternary ice ages. Despite their abundance, the factors that keep species apart and prevent hybridization are often unknown. Here, we study a very narrow contact zone between three closely related butterfly species of the Erebia tyndarus species complex. Using genomic data, we first determined whether gene flow occurs and then assessed whether it might be hampered by differences in chromosome number between some species. We found interspecific gene flow between sibling species that differ in karyotype by one chromosome. Conversely, only F1 hybrids occurred between two species that have the same karyotype, forming a steep genomic cline. In a second step, we fitted clines to phenotypic, ecological and parasitic data to identify the factors associated with the genetic cline. We found clines for phenotypic data and the prevalence of the endosymbiont parasite Wolbachia to overlap with the genetic cline, suggesting that they might be drivers for separating the two species. Overall, our results highlight that some gene flow is possible between closely related species despite different chromosome numbers, but that other barriers restrict such gene flow.     

Different cell size and cell number contribution in two newly established and one ancient body size cline of Drosophila subobscura

Latitudinal genetic clines in body size occur in many ectotherms including Drosophila species. In the wing of D. melanogaster, these clines are generally based on latitudinal variation in cell number. In contrast, differences in wing area that evolve by thermal selection in the laboratory are in general based on cell size. To investigate possible reasons for the different cellular bases of these two types of evolutionary response, we compared the newly established North and South American wing size clines of Drosophila subobscura. The new clines are based on latitudinal variation in cell area in North America and cell number in South America. The ancestral European cline is also based on latitudinal variation in cell number. The difference in the cellular basis of wing size variation in the American clines, which are roughly the same age, together with the similar cellular basis of the new South American cline and the ancient European one, suggest that the antiquity of a cline does not explain its cellular basis. Furthermore, the results indicate that wing size as a whole, rather than its cellular basis, is under selection. The different cellular bases of different size clines are most likely explained either entirely by chance or by different patterns of genetic variance--or its expression--in founding populations.     

QTL mapping reveals a striking coincidence in the positions of genomic regions associated with adaptive variation in body size in parallel clines of Drosophila melanogaster on different continents

Latitudinal genetic clines in body size are common in many ectotherm species and are attributed to climatic adaptation. Here, we use Quantitative Trait Loci (QTL) mapping to identify genomic regions associated with adaptive variation in body size in natural populations of Drosophila melanogaster from extreme ends of a cline in South America. Our results show that there is a significant association between the positions of QTL with strong effects on wing area in South America and those previously reported in a QTL mapping study of Australian cline end populations (P < 0.05). In both continents, the right arm of the third chromosome is associated with QTL with the strongest effect on wing area. We also show that QTL peaks for wing area and thorax length are associated with the same genomic regions, indicating that the clinal variation in the body size traits may have a similar genetic basis. The consistency of the results found for the South American and Australian cline end populations indicate that the genetic basis of the two clines may be similar and future efforts to identify the genes producing the response to selection should be focused on the genomic regions highlighted by the present work.     

Cline coupling and uncoupling in a stickleback hybrid zone

Strong ecological selection on a genetic locus can maintain allele frequency differences between populations in different environments, even in the face of hybridization. When alleles at divergent loci come into tight linkage disequilibrium, selection acts on them as a unit and can significantly reduce gene flow. For populations interbreeding across a hybrid zone, linkage disequilibria between loci can force clines to share the same slopes and centers. However, strong ecological selection on a locus can also pull its cline away from the others, reducing linkage disequilibrium and weakening the barrier to gene flow. We looked for this “cline uncoupling” effect in a hybrid zone between stream resident and anadromous sticklebacks at two genes known to be under divergent natural selection (Eda and ATP1a1) and five morphological traits that repeatedly evolve in freshwater stickleback. These clines were all steep and located together at the top of the estuary, such that we found no evidence for cline uncoupling. However, we did not observe the stepped shape normally associated with steep concordant clines. It thus remains possible that these clines cluster together because their individual selection regimes are identical, but this would be very surprising given their diverse roles in osmoregulation, body armor, and swimming performance.     

Genetic Screening Identifies Cyanogenesis-Deficient Mutants of Lotus japonicus and Reveals Enzymatic Specificity in Hydroxynitrile Glucoside Metabolism

Cyanogenesis, the release of hydrogen cyanide from damaged plant tissues, involves the enzymatic degradation of amino acid–derived cyanogenic glucosides (α-hydroxynitrile glucosides) by specific β-glucosidases. Release of cyanide functions as a defense mechanism against generalist herbivores. We developed a high-throughput screening method and used it to identify cyanogenesis deficient (cyd) mutants in the model legume Lotus japonicus. Mutants in both biosynthesis and catabolism of cyanogenic glucosides were isolated and classified following metabolic profiling of cyanogenic glucoside content. L. japonicus produces two cyanogenic glucosides: linamarin (derived from Val) and lotaustralin (derived from Ile). Their biosynthesis may involve the same set of enzymes for both amino acid precursors. However, in one class of mutants, accumulation of lotaustralin and linamarin was uncoupled. Catabolic mutants could be placed in two complementation groups, one of which, cyd2, encoded the β-glucosidase BGD2. Despite the identification of nine independent cyd2 alleles, no mutants involving the gene encoding a closely related β-glucosidase, BGD4, were identified. This indicated that BGD4 plays no role in cyanogenesis in L. japonicus in vivo. Biochemical analysis confirmed that BGD4 cannot hydrolyze linamarin or lotaustralin and in L. japonicus is specific for breakdown of related hydroxynitrile glucosides, such as rhodiocyanoside A. By contrast, BGD2 can hydrolyze both cyanogenic glucosides and rhodiocyanosides. Our genetic analysis demonstrated specificity in the catabolic pathways for hydroxynitrile glucosides and implied specificity in their biosynthetic pathways as well. In addition, it has provided important tools for elucidating and potentially modifying cyanogenesis pathways in plants.     

A comparison of nuclear and mitochondrial cline shapes in a hybrid zone in the Sceloporus grammicus complex (Squamata; Phrynosomatidae)

The F5 and FM2 chromosome races of the Sceloporus grammicus complex form a hybrid zone in the Mexican state of Hidalgo. Previous studies of this zone have assessed genetic structure by averaging estimates of shape and width across three diagnostic chromosome markers. This approach is likely to mask subtle differences in cline shape among loci (e.g. selected vs. neutral), and obscure any displacement of cline centres (if present). Here we use maximum likelihood methods to construct the best fitting individual clines for three chromosomal markers, and also add two new markers; the mitochondrial DNA (mtDNA) locus, and the nuclear ribosomal DNA (rDNA) repeat. For each locus, hybrid zone models were fitted by cline shape and width, and the position and number of segments describing the centre of the zone. Pairwise comparisons between all clines revealed concordance between chromosomes 2 and 6, but significant discordance in cline structure among all other paired combinations. The concordance of chromosomes 2 and 6 suggests that these clines are maintained by genome-wide forces. The discordance of the chromosome 1 cline suggests an influence of asymmetric introgression, while the mtDNA cline is probably influenced by selection and drift. The rDNA locus reveals a pattern best explained by either extreme asymmetric introgression or gene conversion. The structure of zone indicates that genome-wide processes and locus specific selective forces as well as drift, are operating to different degrees on different loci. The locus-by-locus approach used here permits a finer discrimination among possible mechanisms responsible for the maintenance of the individual clines.     

Learned Vocal Variation Is Associated with Abrupt Cryptic Genetic Change in a Parrot Species Complex

Contact zones between subspecies or closely related species offer valuable insights into speciation processes. A typical feature of such zones is the presence of clinal variation in multiple traits. The nature of these traits and the concordance among clines are expected to influence whether and how quickly speciation will proceed. Learned signals, such as vocalizations in species having vocal learning (e.g. humans, many birds, bats and cetaceans), can exhibit rapid change and may accelerate reproductive isolation between populations. Therefore, particularly strong concordance among clines in learned signals and population genetic structure may be expected, even among continuous populations in the early stages of speciation. However, empirical evidence for this pattern is often limited because differences in vocalisations between populations are driven by habitat differences or have evolved in allopatry. We tested for this pattern in a unique system where we may be able to separate effects of habitat and evolutionary history. We studied geographic variation in the vocalizations of the crimson rosella (Platycercus elegans) parrot species complex. Parrots are well known for their life-long vocal learning and cognitive abilities. We analysed contact calls across a ca 1300 km transect encompassing populations that differed in neutral genetic markers and plumage colour. We found steep clinal changes in two acoustic variables (fundamental frequency and peak frequency position). The positions of the two clines in vocal traits were concordant with a steep cline in microsatellite-based genetic variation, but were discordant with the steep clines in mtDNA, plumage and habitat. Our study provides new evidence that vocal variation, in a species with vocal learning, can coincide with areas of restricted gene flow across geographically continuous populations. Our results suggest that traits that evolve culturally can be strongly associated with reduced gene flow between populations, and therefore may promote speciation, even in the absence of other barriers.     

Neutral and adaptive drivers of genomic change in introduced brook trout (Salvelinus fontinalis) populations revealed by pooled sequencing

Understanding the drivers of successful species invasions is important for conserving native biodiversity and for mitigating the economic impacts of introduced species. However, whole‐genome resolution investigations of the underlying contributions of neutral and adaptive genetic variation in successful introductions are rare. Increased propagule pressure should result in greater neutral genetic variation, while environmental differences should elicit selective pressures on introduced populations, leading to adaptive differentiation. We investigated neutral and adaptive variation among nine introduced brook trout (Salvelinus fontinalis) populations using whole‐genome pooled sequencing. The populations inhabit isolated alpine lakes in western Canada and descend from a common source, with an average of ~19 (range of 7–41) generations since introduction. We found some evidence of bottlenecks without recovery, no strong evidence of purifying selection, and little support that varying propagule pressure or differences in local environments shaped observed neutral genetic variation differences. Putative adaptive loci analysis revealed nonconvergent patterns of adaptive differentiation among lakes with minimal putatively adaptive loci (0.001%–0.15%) that did not correspond with tested environmental variables. Our results suggest that (i) introduction success is not always strongly influenced by genetic load; (ii) observed differentiation among introduced populations can be idiosyncratic, population‐specific, or stochastic; and (iii) conservatively, in some introduced species, colonization barriers may be overcome by support through one aspect of propagule pressure or benign environmental conditions.     

The influence of gene flow on latitudinal clines of photoperiodic adult diapause in the Drosophila auraria species-complex

Temperate species belonging to the Drosophila auraria species complex, D. auraria Peng, D. biauraria Bock & Wheeler, D. triauraria Bock & Wheeler and D. subauraria Kimura, enter reproductive diapause to pass the winter in response to short daylengths. These species from Japan showed latitudinal clines in critical daylength which is longer in populations from higher latitudes. The slopes of these clines coincided well with that of the cline which is approximately predicted from climatic data, suggesting that these clines result from adaptation of the species to the latitudinal gradient of climatic conditions. Between the mainlands and the surrounding islands of Japan, the slopes of clines did not differ significantly, but the deviation from the regression line was usually smaller in mainland populations. It is assumed that gene flow reduces the genetic variation among mainland populations and results in the development of smooth clines. In the plain of east China, D. triauraria did not show clinal variation in critical daylength, although the development of the cline is expected from climatic data. Extensive gene flow among Chinese populations is considered to prevent the development of a cline.     

Parapatric diversification after post-glacial range expansion in the gall fly Urophora cardui (Tephritidae)

Aim Primary and secondary genetic clines in post-glacial colonized regions have different implications for biogeographic distributions and the origin of species. Primary clines arise in situ after colonization as adaptive responses to environmental gradients, while secondary clines are caused by contact between vicariant lineages. Here we analyse primary versus secondary origin of a genetic cline in the tephritid fly Urophora cardui in Jutland, Denmark, in a post-glacial landscape. Location Western Palaearctic. Methods Phylogeographic and demographic analyses of U. cardui based on mitochondrial DNA (mtDNA) genealogies, hierarchical genetic variance tests based on allozymes and distribution analysis of a rare allele from the Jutland cline. Results There was no phylogeographic divergence between the Jutland population of U. cardui north of the cline and neighbouring western European regional populations, which all shared the common western European mtDNA haplotype H1. At nuclear loci, by contrast, the North Jutland population was diverged above the mean level of divergence among regional populations and had no loss of genetic variation. A rare allozyme allele that was frequent in the cline area (up to 45%) and was missing north of the cline also occurred at low frequency (0–14%) elsewhere in the sampling range. Shallow phylogeographic divergence was observed between Russian and western European populations and between English and continental populations. Main conclusions The genetic variation patterns support primary cline evolution and parapatric divergence in Jutland following a demographic expansion of a western European ancestral source population of U. cardui, and suggest cryptic refugia and/or selection in other European population assemblages. The patterns of intra-specific regional divergence are discussed with respect to the interpretation of cryptic refugia in Europe after the most recent ice age.