Refugia, differentiation and postglacial migration in arctic-alpine Eurasia, exemplified by the mountain avens (Dryas octopetala L.)

Many arctic-alpine organisms have vast present-day ranges across Eurasia, but their history of refugial isolation, differentiation and postglacial expansion is poorly understood. The mountain avens, Dryas octopetala sensu lato, is a long-lived, wind-dispersed, diploid shrub forming one of the most important components of Eurasian tundras and heaths in terms of biomass. We address differentiation and migration history of the species with emphasis on the western and northern Eurasian parts of its distribution area, also including some East Greenlandic and North American populations (partly referred to as the closely related D. integrifolia M. Vahl). We analysed 459 plants from 52 populations for 155 amplified fragment length polymorphisms (AFLP) markers. The Eurasian plants were separated into two main groups, probably reflecting isolation and expansion from two major glacial refugia, situated south and east of the North European ice sheets, respectively. Virtually all of northwestern Europe as well as East Greenland have been colonized by the Southern lineage, whereas northwest Russia, the Tatra Mountains and the arctic archipelago of Svalbard have been colonized by the Eastern lineage. The data indicate a contact zone between the two lineages in northern Scandinavia and possibly in the Tatra Mountains. The two single populations analysed from the Caucasus and Altai Mountains were most closely related to the Eastern lineage but were strongly divergent from the remaining eastern populations, suggesting survival in separate refugia at least during the last glaciation. The North American populations grouped with those from East Greenland, irrespective of their taxonomic affiliation, but this may be caused by independent hybridization with D. integrifolia and therefore not reflect the true relationship between populations from these areas.     

Genetic variation of Eryngium campestre L. (Apiaceae) in Central Europe

In Germany, Eryngium campestre is restricted to dry habitats along the rivers Rhine and Elbe and to a few areas in Central Germany. This distribution pattern is usually regarded as a typical pattern of postglacial immigration. In the present study, we investigated whether these two geographically distinct distribution areas are genetically differentiated and whether conclusions can be drawn regarding colonization history. To analyse the phylogeographic structure of E. campestre in Central Europe, 278 individuals from 29 populations within Germany and from further reference populations within Europe were analysed. We applied amplified fragment length polymorphisms to examine their genetic relatedness. Our analyses revealed three groups: a Mediterranean group additionally including two Rhine populations; a Rhine–Main group which further includes the westernmost population from the central German dry area; and one group which includes all eastern populations. Our results show that the two geographically distinct areas are genetically differentiated. As genetic diversity within the Elbe populations is very low, we conclude that this area, which was strongly affected through the late glacial maximum, was colonized relatively recently. High genetic diversity in the Rhine populations indicates a contact zone where lineages of different origin met. This would imply that today's patterns of genetic variation were caused through glacial range contractions and expansions. The present study is one of the first studies that deal with the postglacial distribution pattern of a dry grassland plant species in Central Europe and the results suggest that a survival of E. campestre at least during the Dryas cold stage might be possible.     

The genetic basis of color‐related local adaptation in a ring‐like colonization around the Mediterranean

Uncovering the genetic basis of phenotypic variation and the population history under which it established is key to understand the trajectories along which local adaptation evolves. Here, we investigated the genetic basis and evolutionary history of a clinal plumage color polymorphism in European barn owls (Tyto alba). Our results suggest that barn owls colonized the Western Palearctic in a ring‐like manner around the Mediterranean and meet in secondary contact in Greece. Rufous coloration appears to be linked to a recently evolved nonsynonymous‐derived variant of the melanocortin 1 receptor (MC1R) gene, which according to quantitative genetic analyses evolved under local adaptation during or following the colonization of Central Europe. Admixture patterns and linkage disequilibrium between the neutral genetic background and color found exclusively within the secondary contact zone suggest limited introgression at secondary contact. These results from a system reminiscent of ring species provide a striking example of how local adaptation can evolve from derived genetic variation.     

European common frog (Rana temporaria) recolonized Switzerland from multiple glacial refugia in northern Italy via trans‐ and circum‐Alpine routes

The high mountain ranges of Western Europe had a profound effect on the biotic recolonization of Europe from glacial refugia. The Alps present a particularly interesting case because they form an absolute barrier to dispersal for most taxa, obstructing recolonization from multiple refugia in northern Italy. Here, we investigate the effect of the European Alps on the phylogeographic history of the European common frog Rana temporaria. Based on partial cytochrome b and COXI sequences from Switzerland, we find two mitochondrial lineages roughly north and south of the Alpine ridge, with contact zones between them in eastern and western Switzerland. The northern haplogroup falls within the previously identified Western European haplogroup, while the southern haplogroup is unique to Switzerland. We find that the lineages diverged ~110 kya, at approximately the onset of the last glacial glaciation; this indicates that they are from different glacial refugia. Phylogenetic analyses suggest that the northern and southern haplogroups colonized Switzerland via trans‐ and circum‐Alpine routes from at least two separate refugia in northern Italy. Our results illustrate how a complex recolonization history of the central European Alps can arise from the semi‐permeable barrier created by high mountains.     

The influence of Pleistocene glacial refugia on tawny owl genetic diversity and phylogeography in western Europe

The glacial refugia hypothesis indicates that during the height of the Pleistocene glaciations the temperate species that are today widespread in western Europe must have survived in small and climatically favourable areas located in the southern peninsulas of Iberia, Italy and Balkans. One such species is the tawny owl, a relatively sedentary, nonmigratory bird presently distributed throughout Europe. It is a tree-nesting species closely associated with deciduous and mixed coniferous woodlands. In this study I used control region mtDNA sequences from 187 individuals distributed among 14 populations to determine whether current genetic patterns in tawny owl populations were consistent with postglacial expansion from peninsular refugia. European, North African and Asian tawny owls were found to represent three distinct lineages, where North Africa is the sister clade to all European owls. Within Europe, I found three well-supported clades that correspond to each of the three allopatric refugia. Expansion patterns indicate that owls from the Balkan refugium repopulated most of northern Europe, while expansion out of Iberia and Italy had only regional effects leading to admixture in France. Estimates of population divergence times between refugia populations are roughly similar, but one order of magnitude smaller between Greece and northern Europe. Based on a wide range of mutation rates and generation times, divergence between refugia appears to date to the Pleistocene.     

The genetic structure of adders (Vipera berus) in Fennoscandia: congruence between different kinds of genetic markers

In order to elucidate the colonization history of Fennoscandian adders (Vipera berus), the phylogeographical patterns of two nuclear sets of DNA markers (random amplified polymorphic DNA and microsatellite) are compared with that previously obtained from mitochondrial DNA. An eastern and a western lineage within Fennoscandian adders is readily distinguishable using both sets of nuclear markers, corroborating the hypothesis that the lineages stem from separate glacial refugia. Moreover, the same contact zones as were derived from mitochondrial data are clearly identifiable. Both sets of nuclear markers detect a high level of admixture across one zone in northern Finland, with introgression reaching far west into Sweden.     

The postglacial recolonization of Northern Europe by Rana arvalis as revealed by microsatellite and mitochondrial DNA analyses

The postglacial history of the moor frog (Rana arvalis) in Northern Europe was investigated with the aid of eight variable microsatellite loci and a 661 bp sequence of the mitochondrial cytochrome b gene. A division between eastern and western mitochondrial lineages was discovered, supporting two recolonization routes to Fennoscandia since the last glacial maximum. This result was corroborated by the microsatellite data, which revealed a contact zone between the two lineages in Northern Sweden. These findings add to the increasing evidence that an intraspecific genetic biodiversity founded on the existence of eastern and western clades is a common element in Fennoscandian fauna and flora.     

Phylogeography and recolonization of the Swiss Alps by the Valais shrew (Sorex antinorii), inferred with autosomal and sex‐specific markers

Using one male‐inherited, one female‐inherited and eight biparentally inherited markers, we investigate the population genetic structure of the Valais shrew (Sorex antinorii) in the Swiss Alps. Bayesian analysis on autosomal microsatellites suggests a clear genetic differentiation between two groups of populations. This geographically based structure is consistent with two separate postglacial recolonization routes of the species into Switzerland from Italian refugia after the last Pleistocene glaciations. Sex‐specific markers also confirm genetic structuring among western and eastern areas, since very few haplotypes for either Y chromosome or mtDNA genome are shared between the two regions. Overall, these results suggest that two already well‐differentiated genetic lineages colonized the Swiss Alps and came into secondary contact in the Rhône Valley. Low level of admixture between the two lineages is likely explained by the mountainous landscape structure of lateral valleys orthogonal to the main Rhône valley.     

Phylogeography of the white-tailed eagle, a generalist with large dispersal capacity

Aim Late Pleistocene glacial changes had a major impact on many boreal and temperate taxa, and this impact can still be detected in the present‐day phylogeographic structure of these taxa. However, only minor effects are expected in species with generalist habitat requirements and high dispersal capability. One such species is the white‐tailed eagle, Haliaeetus albicilla, and we therefore tested for the expected weak population structure at a continental level in this species. This also allowed us to describe phylogeographic patterns, and to deduce Ice Age refugia and patterns of postglacial recolonization of Eurasia. Location Breeding populations from the easternmost Nearctic (Greenland) and across the Palaearctic (Iceland, continental Europe, central and eastern Asia, and Japan). Methods Sequencing of a 500 base‐pair fragment of the mitochondrial DNA control region in 237 samples from throughout the distribution range. Results Our analysis revealed pronounced phylogeographic structure. Overall, low genetic variability was observed across the entire range. Haplotypes clustered in two distinct haplogroups with a predominantly eastern or western distribution, and extensive overlap in Europe. These two major lineages diverged during the late Pleistocene. The eastern haplogroup showed a pattern of rapid population expansion and colonization of Eurasia around the end of the Pleistocene. The western haplogroup had lower diversity and was absent from the populations in eastern Asia. These results suggest survival during the last glaciation in two refugia, probably located in central and western Eurasia, followed by postglacial population expansion and admixture. Relatively high genetic diversity was observed in northern regions that were ice‐covered during the last glacial maximum. This, and phylogenetic relationships between haplotypes encountered in the north, indicates substantial population expansion at high latitudes. Areas of glacial meltwater runoff and proglacial lakes could have provided suitable habitats for such population growth. Main conclusions This study shows that glacial climate fluctuations had a substantial impact on white‐tailed eagles, both in terms of distribution and demography. These results suggest that even species with large dispersal capabilities and relatively broad habitat requirements were strongly affected by the Pleistocene climatic shifts.     

Phylogeography and postglacial dispersal of arctic charr Salvelinus alpinus in North America

Restriction fragment length polymorphism analysis of mitochondrial DNA (mtDNA) was used to reconstruct postglacial dispersal routes of arctic charr Salvelinus alpinus in North America. Twelve of 35 restriction enzymes detected polymorphisms among representative populations, revealing two distinct lineages with an estimated nucleotide divergence of 1.32%. Subsequent screening of 869 fish from 54 populations with four diagnostic restriction enzymes showed that these lineages have largely allopatric distributions, suggesting their dispersal from separate northern and eastern glacial refugia. In addition, geographical and genetic structure among eastern populations suggested the existence of a second eastern refuge. Among the three lineages, the most divergent (Arctic) lineage occurred from Alaska east to northern Labrador. Quebec, New Brunswick, and New England were colonized by a second (Laurentian) lineage, and Labrador by a third group. Contact between refugial groups was only detected in two Labrador populations. The Arctic lineage was highly differentiated from eastern North American and European haplotypes, and probably diverged during the early Pleistocene. By contrast, the Laurentian and Labrador groups were similar to Old World charr, suggesting a shared ancestry during the mid-Pleistocene. In addition, the close relationship between Labrador and Laurentian charr indicates their probable divergence during the Wisconsinan glaciation.     

Evidence of a highly complex phylogeographic structure on a specialist river bird species, the dipper (Cinclus cinclus)

This study details the phylogeographic pattern of the white-throated dipper (Cinclus cinclus), a Palearctic, temperate, passerine bird that is exclusively associated with flowing water. Our results reveal a complex phylogeographic structure with at least five distinct lineages for the Western Palearctic region. As for many species of the Western Palearctic fauna and flora, this genetic structure is probably linked to the isolation of populations in different southern refuges during glacial periods. Furthermore, the isolation of populations in Scandinavia and/or Eastern regions, but also in Morocco and probably in Corsica, was accentuated by ecological and biogeographic barriers during Quaternary interglacial periods. During glacial periods, Italy, Sicily and the Balkano-Carpathian region acted as major refuge zones for the dipper. At the end of the last ice age, Western Europe was repopulated by dippers from an Italian refuge, while Eastern Europe was recolonised by Balkano-Carpathian birds. A large contact zone between these two lineages was evidenced and extends from Luxembourg to Hungary. Finally, our results indicate the need to clarify the taxonomic status of the dipper, especially concerning the European subspecies whose validity appears uncertain.     

Phylogeography and genetic structuring of European nine-spined sticklebacks (Pungitius pungitius)-mitochondrial DNA evidence

As a consequence of colonisation from different glacial refugia, many northern European taxa are split into distinct western and eastern lineages. However, as for the nine-spined stickleback (Pungitius pungitius), the exact location of the contact zone between lineages often remains poorly known. We assessed the genetic differentiation and diversity in the nine-spined stickleback within Europe using 1037 base pairs of cytochrome b sequence for 320 individuals from 57 locations, including pond, lake, river, and coastal habitats. Our main aims were (i) to locate the contact zone between the previously recognized western and eastern lineages, (ii) investigate latitudinal patterns in genetic diversity, (iii) compare genetic diversity among different habitat types, and (iv) date the known split between eastern and western lineages. The data revealed the split between eastern and western to be located across the Danish Straits and roughly following the Norway/Sweden border to the North. Reference sites from Canada form their own clades, and one of the Canadian sites was found to have a haplotype common to the Eastern European lineage, possibly representing an ancestral polymorphism. The split between the two European clades was dated to approximately 1.48 million years ago (Mya), and between Canada and Europe to approximately 1.62 Mya. After controlling for habitat effects, nucleotide (but not haplotype) diversity across populations decreased with increasing latitude. Coastal populations showed significantly higher haplotype diversity (but not nucleotide diversity) than pond populations, but there were no detectable differences in haplotype diversity among different freshwater habitat types (viz. river, lake and pond populations), or between coastal and lake/river populations. Sequences were found to cluster according to their geographic proximity, rather than by habitat type, and all habitat types were found within each major clade, implying that colonisation and adaptation between the coastal and freshwater environments in different regions must have occurred in parallel.     

Cryptic speciation, habitat differentiation, and geography in Hamatocaulis vernicosus (Calliergonaceae, Bryophyta)

A molecular study of the wetland pleurocarpous moss Hamatocaulis vernicosus (Mitt.) Hedenäs, based on the ITS region, the trnL-trnF region, and rpl16 revealed two cryptic species. One of these is widespread in Europe, and was also sampled from several localities in Minnesota (United States). The other cryptic species occurs only south of the boreal zone in Europe, and was in addition sampled from single localities in Peru and northernmost Asiatic Russian Federation. Although it seems likely that the southern species colonized parts of Europe from southern refugia after the last glaciation, its wide distribution outside Europe indicates that the speciation event is unrelated to the European glacial history. No correlations could be found between variation in habitat preferences (pH and electric conductivity) and molecular variation.     

Limited hybridization along a large contact zone between two genetic lineages of the butterfly Erebia medusa (Satyrinae, Lepidoptera) in Central Europe

Genetic lineages evolving during glacial isolation frequently come into contact as the result of postglacial range expansions. Hybridization often occurs along these contact zones. In Europe, the high mountain systems of the Alps and Pyrenees are well known for their hybrid belts. This article studies the contact zone of the Woodland Ringlet Erebia medusa in a Hercynian mountain area in the Czech-German border region not exceeding 1500 m a.s.l. Hybrid populations between an eastern and a western genetic lineage were detected by amova -based tests, principal component analysis and neighbour joining analysis. Over most of the range of the Czech-German Border Mts, the ridges separate the western and the eastern genetic lineage from each other. However, two important hybrid areas were detected: (1) the watershed of the Ohře river in the north-west of that area, a major valley system passing through these mountains and (2) the high plateaux of the Šumava Mts in the south-east, an extended area of high elevation. The location of hybrid populations in geographical vicinity to non-hybrid populations and the generally low F _IS (2.1%) make reduced fitness of hybrid individuals little likely. The hybrid populations have intermediate genetic diversity between the genetically poor western and the genetically rich eastern lineage populations.     

Phylogeography of field voles (Microtus agrestis) in Eurasia inferred from mitochondrial DNA sequences

In a distribution-wide phylogeographic survey of the field vole (Microtus agrestis), 75 specimens from 56 localities across Eurasia were examined for DNA sequence variation along the whole 1140 base pair (bp) mitochondrial (mt) cytochrome b gene. The species is subdivided into three main mtDNA phylogeographic groups - western, eastern and southern - with largely allopatric distributions. The western phylogeographical group is found in west and central Europe and spread most probably from a glacial refugium in the Carpathians. The eastern group covers a large range from Lithuania to central Asia, and probably originated from a southeast European source (e.g. the southern Urals or the Caucasus). The southern group occupies an area from Portugal to Hungary, with division into two distinct mtDNA sublineages that presumably derive from separate glacial refugia in the Iberian Peninsula. Molecular clock estimates suggest that the western and eastern field vole populations separated during the last glaciation, whereas the southern population dates back 0.5-0.9 Myr. High levels of mtDNA variation indicate relatively large population sizes and subdivisions within phylogeographic groups during the last glaciation. We report a possible new suture zone in east Europe.     

Evolutionary history of the Dactylorhiza maculata polyploid complex (Orchidaceae)

Taxonomic complexity may be associated with migration history and polyploidy. We used plastid and nuclear DNA markers to investigate the evolutionary history of the systematically challenging Dactylorhiza maculata polyploid complex. A total of 1833 individuals from 298 populations from throughout Europe were analysed. We found that gene flow was limited between the two major taxa, diploid ssp. fuchsii (including ssp. saccifera) and tetraploid ssp. maculata. A minimum of three autotetraploid lineages were discerned: (1) southern/western ssp. maculata; (2) northern/eastern ssp. maculata; and (3) Central European ssp. fuchsii. The two ssp. maculata lineages, which probably pre‐date the last glaciation, form a contact zone with high genetic diversity in central Scandinavia. Intermediate plastid haplotypes in the contact zone hint at recombination. Central Europe may have been a source area for the postglacial migration for the southern/western lineage of ssp. maculata, as well as for ssp. fuchsii. The northern/eastern lineage of ssp. maculata may have survived the LGM in central Russia west of the Urals. The tetraploid lineage of ssp. fuchsii is indistinguishable from diploid ssp. fuchsii, and is probably of postglacial origin. The Mediterranean region and the Caucasus have not contributed to the northward migration of either ssp. fuchsii or ssp. maculata. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 101 , 503–525.     

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.     

Phylogeography of endemic ermine (Mustela erminea) in southeast Alaska

The North Pacific Coast (NPC) of North America is a region of high mammalian endemism, possibly due to its highly fragmented landscape and complex glacial history. For example, four island and one mainland subspecies of ermine, Mustela erminea, have been described as endemic to southeast Alaska alone. To better understand the role of past climatic change in generating diversity in the region, we examined DNA sequence variation in the mitochondrial cytochrome b gene of 210 ermine from across North America, with an emphasis on Alaska and British Columbia. We found three distinct (1.5-3.6% uncorrected 'p') lineages of ermine, all of which occur in southeast Alaska. One lineage includes a southeast Alaska endemic and specimens from Alaska (outside of southeast) and Eurasia. A second lineage includes two southeast Alaskan endemics and ermine from western Canada and the coterminous United States. The close relationships of these purported endemics to ermine outside of southeast Alaska suggest that they colonized the region from Beringian and southern glacial refugia, respectively, following deglaciation of the NPC. The third lineage appears restricted to the Prince of Wales Island complex in southeast Alaska (two subspecies) and Graham Island (Haida Gwaii), British Columbia. This restricted distribution suggests that these populations may be derived from relicts that persisted in a coastal refugium during the Wisconsin glaciation. Studies of nuclear genes and adaptive morphological evolution are necessary to further explore discrepancies between the geographical pattern of differentiation based on mtDNA and the existing subspecific taxonomy based on morphology.     

East Asian allopatry and north Eurasian sympatry in Long‐tailed Tit lineages despite similar population dynamics during the late Pleistocene

Eurasia is a large continent characterized by heterogeneous environments. Glacial cycles during the late Pleistocene have had variable impacts on the avifauna across Eurasia. Bird populations from South‐East Asia show stability through the Last Glacial Maximum (LGM), while populations from Europe exhibit evidence of post‐LGM expansion. We investigated the phylogeography of the Long‐tailed Tit (Aegithalos caudatus), which spans the longitudinal breadth of Eurasia to test how climatic history and regional topographical complexity affected populations and diversification within the species complex. Our results show that two lineages from central and southern China (lineages C and D) segregate geographically, while lineages across northern Eurasia (lineage A and B) show substantial sympatry. Bayesian estimates for the timing of diversification suggest that the four lineages diverged during the middle Pleistocene, splitting in parallel and undergoing concurrent demographic histories since divergence. A. caudatus lineages experienced similar and synchronous population size dynamics during glacial cycles before the LGM. We conclude that the difference in geo‐topologic complexity may be an important factor that led to the variation in secondary admixture between northern Eurasian and eastern Asian lineages.     

Speciation in Mediterranean refugia and post‐glacial expansion of Zerynthia polyxena (Lepidoptera,Papilionidae)

Migration of populations to and from glacial refugia is responsible for various cases of speciation and subspeciation in Europe. The pattern of distribution and the degree of diversification between lineages originated by isolation in different glacial refugia usually depends on ecological traits, especially to their dispersal ability. Zerynthia polyxena is a philopatric species, scattered in small populations and rarely colonizing mountain areas. These characteristics probably caused repeated isolation during the Quaternary and may have favoured diversification. Actually two studies, based on both morphological and genetic data, suggest the existence of two highly distinct lineages in Europe having in Northern Italy their contact zone. In this study, I applied geometric morphometrics to male genitalia and demonstrated that (i) two morphotypes exist in Europe approximately facing on the two sides of the Po River; (ii) the two lineages probably survived glaciations in Italy and the Balkan Peninsula, respectively; then the Balkans lineage expanded to Central and Eastern Europe; (iii) no hybrid populations seem to exist in the contact area and, in one locality at least, the two lineages live in sympatry without any evidence of intermediates. These results suggest that (i) two sister species of Zerynthia exist in Europe. Accordingly, Papilio cassandra Geyer, 1828 is reinstated, as Zerynthia cassandra stat. rev., as the species to which the Zerynthia from Italy South of the Po River belong. Male genitalia differences with Zerynthia polyxena are described.