Genetic and phenotypic variation across a hybrid zone between ecologically divergent tree squirrels (Tamiasciurus)

A hybrid zone along an environmental gradient should contain a clinal pattern of genetic and phenotypic variation. This occurs because divergent selection in the two parental habitats is typically strong enough to overcome the homogenizing effects of gene flow across the environmental transition. We studied hybridization between two parapatric tree squirrels (Tamiasciurus spp.) across a forest gradient over which the two species vary in coloration, cranial morphology and body size. We sampled 397 individuals at 29 locations across a 600‐km transect to seek genetic evidence for hybridization; upon confirming hybridization, we examined levels of genetic admixture in relation to maintenance of phenotypic divergence despite potentially homogenizing gene flow. Applying population assignment analyses to microsatellite data, we found that Tamiasciurus douglasii and T. hudsonicus form two distinct genetic clusters but also hybridize, mostly within transitional forest habitat. Overall, based on this nuclear analysis, 48% of the specimens were characterized as T. douglasii, 9% as hybrids and 43% as T. hudsonicus. Hybrids appeared to be reproductively viable, as evidenced by the presence of later‐generation hybrid genotypes. Observed clines in ecologically important phenotypic traits—fur coloration and cranial morphology—were sharper than the cline of putatively neutral mtDNA, which suggests that divergent selection may maintain phenotypic distinctiveness. The relatively recent divergence of these two species (probably late Pleistocene), apparent lack of prezygotic isolating mechanisms and geographic coincidence of cline centres for both genetic and phenotypic variation suggest that environmental factors play a large role in maintaining the distinctiveness of these two species across the hybrid zone.     

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.     

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.     

Molecular and morphological analysis of secondary contact zones of Cottus gobio in Fennoscandia: geographical discordance of character transitions

Northern Europe was postglacially colonized from different directions by distinct phylogeographical lineages of the bullhead Cottus gobio L. (Pisces: Scorpaeniformes). These lineages have then come into contact in coastal habitats of the currently brackish Baltic Sea and in the freshwaters north of it. We studied the patterns of intergradation in the contact zones in four morphometric and six molecular characters. In the north, intergradation between the western (W) and eastern (E) bullhead lineages is found both among rivers (west‐to‐east) and along individual rivers (south‐to‐north). The locations of the transition zones probably relate to the timing of the initial contact, subsequent Baltic shoreline displacement (i.e. emergence of the lower river reaches), and dispersal barriers caused by variations of coastal salinity. The transitions (clines) in different characters are, however, not geographically coincident. Mitochondrial DNA clines are generally found upstream and to the east of the other transitions, and GPI‐1 allozyme clines are mostly shifted downstream in the rivers, and west of the other transitions on the broader scale of the Baltic Sea. The location of the mtDNA clines may best reflect the initial contact between lineages, and the displacement of the other clines could result from dispersal being overall asymmetric (predominantly downstream) and sex‐biased (stronger in males). Alternatively, the non‐coincidence might reflect selection against deleterious cytonuclear character combinations. No clear evidence of reproductive incompatibility between the lineages was seen in local population structures; no remaining genetic correlations were observed locally among traits. In another transition area, a coastal transect in southern Finland, clinal patterns similar to those in the northern contact zone were recorded, but the population compositions could not be explained by simple in situ mixing of any of the putatively pure, invading refugial lineages. Probably, the bullhead stocks that initially came into contact in this southern study area already represented mixtures of the invading lineages. © 2004 The Linnean Society of London, Biological Journal of the Linnean Society, 2004, 81 , 535–552.     

Panmixia on a continental scale in a widely distributed colonial waterbird

Many highly mobile species, such as migratory birds, can move and disperse over long distances, yet exhibit high levels of population genetic structuring. Although movement capabilities may enable dispersal, gene flow may be restricted by behavioural constraints such as philopatry. In the present study, we examined patterns of genetic differentiation across the range of a highly mobile, colonial waterbird. American white pelicans (Pelecanus erythrorhynchos) breed across continental North America and are currently experiencing a range expansion, especially on the eastern range limit. To assess patterns of genetic structuring, we sampled 333 individuals from 19 colonies across their North American range. The use of ten variable microsatellite markers revealed high levels of allelic richness with no population differentiation. Both Bayesian and frequentist approaches to examining genetic structuring revealed a single panmictic population. We found no evidence of genetic structuring across the Continental Divide or between migratory and non‐migratory colonies. The lack of any genetic structure across the range indicates that, unlike other waterbirds with similar life‐history characteristics, extensive gene flow and presumably low philopatry appear to preclude genetic differentiation. The lack of population genetic structure in American white pelicans provides an example of range‐wide panmixia, a rare phenomenon in any terrestrial species. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 102 , 583–592.     

THE EFFECT OF A PARTIAL BARRIER ON THE MOVEMENT OF A HYBRID ZONE

Computer simulations of clines (Brues, 1972; Endler, 1977) as well as theoretical arguments (Nagylaki, 1975), have shown that steps in gene frequencies will be pulled to partial barriers (areas of reduced gene flow) if they form within approximately a cline width of the partial barrier. The behavior of a hybrid zone between two chromosomal taxa (“Moreton” and “Torresian”) of the acridine grasshopper Caledia captiva in southeast Queensland has been analyzed and found to comform qualitatively with a model of altered gene flow patterns. Clines in four enzyme systems were analyzed for 1983 and 1986 along a transect across the hybrid zone. The clines have shifted towards an area of regenerating forest, while homozygote frequencies have increased at this point. This forest barrier has broken the continuity of the spatial distribution of C. captiva, forming population islands in part of the hybrid zone, and thus reducing the amount of gene flow. The distance between the barrier and the original cline is approximately of the order of a cline width, so that they would be expected to interact. Historical information suggests that the secondary contact between the “Moreton” and “Torresian” taxa occurred very recently (1844–1940), due to the intensive land-clearing activities during the European settlement.     

A clinally varying promoter polymorphism associated with adaptive variation in wing size in Drosophila

Body size often shows adaptive clines in many ectotherms across altitude and latitude, but little is known about the genetic basis of these adaptive clines. Here we identify a polymorphism in the Dca (Drosophila cold acclimation) gene in Drosophila melanogaster that influences wing size, affects wing:thorax allometry and also controls a substantial proportion of the clinal wing‐size variation. A polymorphism in the promoter region of Dca had two common alleles showing strong reciprocal clinal variation in frequency with latitude along the east coast of Australia. The Dca‐237 allele increased towards the tropics where wing size is smaller. A within‐population association study highlighted that an increase in the frequency of this allele decreased wing size but did not influence thorax size. A manipulated increase in the level of expression of Dca achieved through UAS‐GAL4 was associated with a decrease in wing size but had no effect on thorax size. This was consistent with higher Dca expression levels in family lines with higher frequency of the Dca‐237 allele. Genetic variation in the promoter region of the Dca gene appears to influence adaptive size variation in the eastern Australian cline of Drosophila melanogaster and accounts for more than 10% of the genetic variation in size within and between populations.     

Complex patterns of local adaptation in heat tolerance in Drosophila simulans from eastern Australia

Latitudinal clines are considered a powerful means of investigating evolutionary responses to climatic selection in nature. However, most clinal studies of climatic adaptation in Drosophila have involved species that contain cosmopolitan inversion polymorphisms that show clinal patterns themselves, making it difficult to determine whether the traits or inversions are under selection. Further, although climatic selection is unlikely to act on only one life stage in metamorphic organisms, a few studies have examined clinal patterns across life stages. Finally, clinal patterns of heat tolerance may also depend on the assay used. To unravel these potentially confounding effects on clinal patterns of thermal tolerance, we examined adult and larval heat tolerance traits in populations of Drosophila simulans from eastern Australia using static and dynamic (ramping 0.06 °C min^?1) assays. We also used microsatellites markers to clarify whether demographic factors or selection are responsible for population differentiation along clines. Significant cubic clinal patterns were observed for adult static basal, hardened and dynamic heat knockdown time and static basal heat survival in larvae. In contrast, static, hardened larval heat survival increased linearly with latitude whereas no clinal association was found for larval ramping survival. Significant associations between adult and larval traits and climatic variables, and low population differentiation at microsatellite loci, suggest a role for climatic selection, rather than demographic processes, in generating these clinal patterns. Our results suggest that adaptation to thermal stress may be species and life‐stage specific, complicating our efforts to understand the evolutionary responses to selection for increasing thermotolerance.     

Comparative population genomics of latitudinal variation in Drosophila simulans and Drosophila melanogaster

Examples of clinal variation in phenotypes and genotypes across latitudinal transects have served as important models for understanding how spatially varying selection and demographic forces shape variation within species. Here, we examine the selective and demographic contributions to latitudinal variation through the largest comparative genomic study to date of Drosophila simulans and Drosophila melanogaster, with genomic sequence data from 382 individual fruit flies, collected across a spatial transect of 19 degrees latitude and at multiple time points over 2 years. Consistent with phenotypic studies, we find less clinal variation in D. simulans than D. melanogaster, particularly for the autosomes. Moreover, we find that clinally varying loci in D. simulans are less stable over multiple years than comparable clines in D. melanogaster. D. simulans shows a significantly weaker pattern of isolation by distance than D. melanogaster and we find evidence for a stronger contribution of migration to D. simulans population genetic structure. While population bottlenecks and migration can plausibly explain the differences in stability of clinal variation between the two species, we also observe a significant enrichment of shared clinal genes, suggesting that the selective forces associated with climate are acting on the same genes and phenotypes in D. simulans and D. melanogaster.     

Patterns of divergence in the olive sunbird Cyanomitra olivacea (Aves: Nectariniidae) across the African rainforest–savanna ecotone

In the debate over modes of vertebrate diversification in tropical rainforests, two competing hypotheses of speciation predominate: those that emphasize the role of geographical isolation during glacial periods and those that stress the role of ecology and diversifying selection across ecotones or environmental gradients. To investigate the relative roles of selection versus isolation in refugia, we contrasted genetic and morphologic divergence of the olive sunbird (Cyanomitra olivacea) at 18 sites (approximately 200 individuals) across the forest–savanna ecotone of Central Africa in a region considered to have harboured three hypothesized refugia during glacial periods. Habitats were characterized using bioclimatic and satellite remote‐sensing data. We found relatively high levels of gene flow between ecotone and forest populations and between refugia. Consistent with a pattern of divergence‐with‐gene‐flow, we found morphological characters to be significantly divergent across the gradient [forest versus ecotone (mean ± SD): wing length 60.47 ± 1.81 mm versus 62.18 ± 1.35 mm; tarsus length 15.51 ± 0.82 mm versus 16.00 ± 0.57 mm; upper mandible length 21.77 ± 1.09 mm versus 23.19 ± 0.98 mm, respectively]. Within‐habitat comparisons across forest and ecotone sites showed no significant differences in morphology. The results show that divergence in morphological traits is tied to environmental variables across the gradient and is occurring despite gene flow. The pattern of divergence‐with‐gene‐flow found is similar to that described for other rainforest species across the gradient. These results suggest that neither refugia, nor isolation‐by‐distance have played a major role in divergence in the olive sunbird, although ecological differences along the forest and savanna ecotone may impose significant selection pressures on the phenotype and potentially be important in diversification. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 103 , 821–835.     

Digest: A clinal polymorphism and life‐history adaptations in Drosophila*

Does genetic variation in the insulin/insulin‐like growth factor signaling pathway (IIS) underlie latitudinal life‐history clines in North American Drosophila melanogaster? Durmaz et al. investigated how a clinally varying polymorphism in the IIS gene foxo affects fitness‐related traits by isolating the effects of alternative low and high latitude alleles. The phenotypic effects of the polymorphism—for example, on body size—matched those normally observed across the cline, suggesting that variation in IIS is important for clinal life‐history adaptation.     

ALTITUDINAL CLINAL VARIATION IN WING SIZE AND SHAPE IN AFRICAN DROSOPHILA MELANOGASTER: ONE CLINE OR MANY?

Geographical patterns of morphological variation have been useful in addressing hypotheses about environmental adaptation. In particular, latitudinal clines in phenotypes have been studied in a number of Drosophila species. Some environmental conditions along latitudinal clines—for example, temperature—also vary along altitudinal clines, but these have been studied infrequently and it remains unclear whether these environmental factors are similar enough for convergence or parallel evolution. Most clinal studies in Drosophila have dealt exclusively with univariate phenotypes, allowing for the detection of clinal relationships, but not for estimating the directions of covariation between them. We measured variation in wing shape and size in D. melanogaster derived from populations at varying altitudes and latitudes across sub‐Saharan Africa. Geometric morphometrics allows us to compare shape changes associated with latitude and altitude, and manipulating rearing temperature allows us to quantify the extent to which thermal plasticity recapitulates clinal effects. Comparing effect vectors demonstrates that altitude, latitude, and temperature are only partly associated, and that the altitudinal shape effect may differ between Eastern and Western Africa. Our results suggest that selection responsible for these phenotypic clines may be more complex than just thermal adaptation.     

Nonrandom larval dispersal can steepen marine clines

Sharp and stable clinal variation is enigmatic when found in species with high gene flow. Classical population genetic models treat gene flow as a random homogenizing force countering local adaptation across habitat discontinuities. Under this view, dispersal over large spatial scales will lower the effectiveness of adaptation by natural selection at finer spatial scales. Thus, random gene flow will create a shallow phenotypic cline across an ecotone in response to a steep selection gradient. In sedentary marine species that disperse primarily as larvae, nonrandom dispersal patterns are expected due to coastal hydrodynamics. Surprisingly sharp phenotypic and genotypic clines have been documented in marine species with high gene flow. We are interested in the extent to which nonrandom dispersal could accentuate such clines. We model a linear species range in which populations have stable and uniform densities along a selection gradient; in contrast to random dispersal, convergent advection of larvae can amplify phenotypic differentiation if coupled with a semipermeable dispersal barrier in the convergence zone. The migration load caused by directional dispersal pushes the phenotypic mean away from the local trait optimum in downstream populations, that is, near the convergence zone. A dispersal barrier is possible as a result of colliding currents if the water and larvae are mostly displaced offshore, away from suitable settlement habitat. Disjunctions in a quantitative trait were enlarged in the convergence zone by faster current flows or a more complete dispersal barrier. With advection of larvae per generation one-third as far as the average dispersal distance by diffusion, convergence on a dispersal barrier with 40% permeability generated a trait disjunction across the convergence zone of two phenotypic standard deviations. Without directional dispersal, similar clines also developed across a habitat gap, where population density was low, or across dispersal barriers with less than 1% permeability. These findings suggest that the types of hydrographic phenomena often associated with marine transition zones can strongly affect the balance between gene flow and selection and generate surprisingly steep clines given the large-scale gene flow expected from larvae.     

Parallel effects of the inversion In(3R)Payne on body size across the North American and Australian clines in Drosophila melanogaster

Chromosomal inversions are thought to play a major role in climatic adaptation. In D. melanogaster, the cosmopolitan inversion In(3R)Payne exhibits latitudinal clines on multiple continents. As many fitness traits show similar clines, it is tempting to hypothesize that In(3R)P underlies observed clinal patterns for some of these traits. In support of this idea, previous work in Australian populations has demonstrated that In(3R)P affects body size but not development time or cold resistance. However, similar data from other clines of this inversion are largely lacking; finding parallel effects of In(3R)P across multiple clines would considerably strengthen the case for clinal selection. Here, we have analysed the phenotypic effects of In(3R)P in populations originating from the endpoints of the latitudinal cline along the North American east coast. We measured development time, egg‐to‐adult survival, several size‐related traits (femur and tibia length, wing area and shape), chill coma recovery, oxidative stress resistance and triglyceride content in homokaryon lines carrying In(3R)P or the standard arrangement. Our central finding is that the effects of In(3R)P along the North American cline match those observed in Australia: standard arrangement lines were larger than inverted lines, but the inversion did not influence development time or cold resistance. Similarly, In(3R)P did not affect egg‐to‐adult survival, oxidative stress resistance and lipid content. In(3R)P thus seems to specifically affect size traits in populations from both continents. This parallelism strongly suggests an adaptive pattern, whereby the inversion has captured alleles associated with growth regulation and clinal selection acts on size across both continents.     

Large scale geographic clines of parasite damage to Populus tremula L.

According to the geographic mosaic theory of coevolution (GMTC), clines of traits reflecting local co‐adaptation (including resistance genes) should be common between a host and its parasite and should persist across time. To test the GMTC‐assumption of persistent clinal patterns we compared the natural prevalence of two parasites on aspen Populus tremula trees: mining moths of the genus Phyllocnistis and leaf rust Melampsora spp. Damage data were collated from the Swedish National Forest Damage Inventory (2004–2006). In addition, occurrence of the parasites was scored in field conditions in two common gardens in the north and south of Sweden over five growing seasons (2004–2008), then related to biomass (stem height and diameter) and to concentrations of eleven leaf phenolics. Phyllocnistis mainly occurred in the northern garden, a distribution range which was confirmed by the countrywide inventory, although Phyllocnistis was more abundant on southern clones, providing evidence for possible local maladaptation. Melampsora occurred all over the country and in both gardens, but built up more quickly on northern clones, which suggests a centre of local clone maladaptation in the north. Stem growth also followed a clinal pattern as did the concentration of three phenolic compounds: benzoic acid, catechin and cinnamic acid. However, only benzoic acid was related to parasite presence: negatively to Phyllocnistis and positively to Melampsora and it could thus be a potential trait under selection. In conclusion, clines of Phyllocnistis were stronger and more persistent compared to Melampsora, which showed contrasting clines of varying strength. Our data thus support the assumption of the GMTC model that clines exist in the border between hot and cold spots and that they may be less persistent for parasites with an elevated gene flow, and/or for parasites which cover relatively larger hot spots surrounded by fewer cold spots.     

Latitudinal clines in body size,but not in thermal tolerance or heat‐shock cognate 70 (HSC70), in the highly‐dispersing intertidal gastropod Littorina keenae (Gastropoda: Littorinidae)

Natural populations of widely‐distributed animals often exhibit clinal variation in phenotypic traits or in allele frequencies of a particular gene over their geographical range. A planktotrophic intertidal snail, Littorina keenae is broadly distributed along the north‐eastern Pacific coast through a large latitudinal range (24°50′N–43°18′N). We tested for latitudinal clines in two complex phenotypic traits – thermal tolerance and body size – and one single locus trait – heat shock cognate 70 (HSC70) – in L. keenae along almost its entire geographical range. We found only weak evidence for a latitudinal cline in the thermal tolerance and no evidence for a cline in allele frequencies at HSC70. However, as predicted by Bergmann's rule, we detected a strong latitudinal cline that accounted for 60% of the variance in body size (R² = 0.598; P < 0.001). In contrast, body size did not significantly affect thermal tolerance. HSC70 showed no genetic differentiation among the populations, supporting our previous mitochondrial gene‐based estimate of high gene flow during this snail's free‐swimming larval stage. Given that L. keenae experiences panmixia along its species range, the observed size cline may be partially or entirely caused by a phenotypically plastic response to local thermal environments rather than by genetic divergence in body size among populations in response to locally optimizing natural selection. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 494–505.     

Genomics of clinal variation in Drosophila: disentangling the interactions of selection and demography

Clines in phenotypes and genotype frequencies across environmental gradients are commonly taken as evidence for spatially varying selection. Classical examples include the latitudinal clines in various species of Drosophila, which often occur in parallel fashion on multiple continents. Today, genomewide analysis of such clinal systems provides a fantastic opportunity for unravelling the genetics of adaptation, yet major challenges remain. A well‐known but often neglected problem is that demographic processes can also generate clinality, independent of or coincident with selection. A closely related issue is how to identify true genic targets of clinal selection. In this issue of Molecular Ecology, three studies illustrate these challenges and how they might be met. Bergland et al. report evidence suggesting that the well‐known parallel latitudinal clines in North American and Australian D. melanogaster are confounded by admixture from Africa and Europe, highlighting the importance of distinguishing demographic from adaptive clines. In a companion study, Machado et al. provide the first genomic comparison of latitudinal differentiation in D. melanogaster and its sister species D. simulans. While D. simulans is less clinal than D. melanogaster, a significant fraction of clinal genes is shared between both species, suggesting the existence of convergent adaptation to clinaly varying selection pressures. Finally, by drawing on several independent sources of evidence, Bo?i?evi? et al. identify a functional network of eight clinal genes that are likely involved in cold adaptation. Together, these studies remind us that clinality does not necessarily imply selection and that separating adaptive signal from demographic noise requires great effort and care.     

DIVERSIFICATION AND GENE FLOW IN NASCENT LINEAGES OF ISLAND AND MAINLAND NORTH AMERICAN TREE SQUIRRELS (TAMIASCIURUS)

Pleistocene climate cycles and glaciations had profound impacts on taxon diversification in the Boreal Forest Biome. Using population genetic analyses with multilocus data, we examined diversification, isolation, and hybridization in two sibling species of tree squirrels (Tamiasciurus douglasii and Tamiasciurus hudsonicus) with special attention to the geographically and genetically enigmatic population of T. hudsonicus on Vancouver Island, Canada. The two species differentiated only about 500,000 years ago, in the Late Pleistocene. The island population is phylogenetically nested within T. hudsonicus according to our nuclear analysis but within T. douglasii according to mitochondrial DNA. This conflict is more likely due to historical hybridization than to incomplete lineage sorting, and it appears that bidirectional gene flow occurred between the island population and both species on the mainland. This interpretation of our genetic analyses is consistent with our bioclimatic modeling, which demonstrates that both species were able to occupy this region throughout the Late Pleistocene. The divergence of the island population 40,000 years ago suggests that tree squirrels persisted in a refugium on Vancouver Island at the last glacial maximum, 20,000 years ago. Our observations demonstrate how Pleistocene climate change and habitat shifts have created incipient divergence in the presence of gene flow.     

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.     

Natural selection in a postglacial range expansion: the case of the colour cline in the European barn owl

Gradients of variation—or clines—have always intrigued biologists. Classically, they have been interpreted as the outcomes of antagonistic interactions between selection and gene flow. Alternatively, clines may also establish neutrally with isolation by distance (IBD) or secondary contact between previously isolated populations. The relative importance of natural selection and these two neutral processes in the establishment of clinal variation can be tested by comparing genetic differentiation at neutral genetic markers and at the studied trait. A third neutral process, surfing of a newly arisen mutation during the colonization of a new habitat, is more difficult to test. Here, we designed a spatially explicit approximate Bayesian computation (ABC) simulation framework to evaluate whether the strong cline in the genetically based reddish coloration observed in the European barn owl (Tyto alba) arose as a by‐product of a range expansion or whether selection has to be invoked to explain this colour cline, for which we have previously ruled out the actions of IBD or secondary contact. Using ABC simulations and genetic data on 390 individuals from 20 locations genotyped at 22 microsatellites loci, we first determined how barn owls colonized Europe after the last glaciation. Using these results in new simulations on the evolution of the colour phenotype, and assuming various genetic architectures for the colour trait, we demonstrate that the observed colour cline cannot be due to the surfing of a neutral mutation. Taking advantage of spatially explicit ABC, which proved to be a powerful method to disentangle the respective roles of selection and drift in range expansions, we conclude that the formation of the colour cline observed in the barn owl must be due to natural selection.