A dynamic history of admixture from Mediterranean and Carpathian glacial refugia drives genomic diversity in the bank vole

Understanding the historical contributions of differing glacial refugia is key to evaluating the roles of microevolutionary forces, such as isolation, introgression, and selection in shaping genomic diversity in present‐day populations. In Europe, where both Mediterranean and extra‐Mediterranean (e.g., Carpathian) refugia of the bank vole (Clethrionomys glareolus) have been identified, mtDNA indicates that extra‐Mediterranean refugia were the main source of colonization across the species range, while Mediterranean peninsulas harbor isolated, endemic lineages. Here, we critically evaluate this hypothesis using previously generated genomic data (>6,000 SNPs) for over 800 voles, focusing on genomic contributions to bank voles in central Europe, a key geographic area in considering range‐wide colonization. The results provide clear evidence that both extra‐Mediterranean (Carpathian) and Mediterranean (Spanish, Calabrian, and Balkan) refugia contributed to the ancestry and genomic diversity of bank vole populations across Europe. Few strong barriers to dispersal and frequent admixture events in central Europe have led to a prominent mid‐latitude peak in genomic diversity. Although the genomic contribution of the centrally located Carpathian refugium predominates, populations in different parts of Europe have admixed origins from Mediterranean (28%–47%) and the Carpathian (53%–72%) sources. We suggest that the admixture from Mediterranean refugia may have provisioned adaptive southern alleles to more northern populations, facilitating the end‐glacial spread of the admixed populations and contributing to increased bank vole diversity in central Europe. This study adds critical details to the complex end‐glacial colonization history of this well‐studied organism and underscores the importance of genomic data in phylogeographic interpretation.     

High genomic diversity in the bank vole at the northern apex of a range expansion: The role of multiple colonizations and end‐glacial refugia

The history of repeated northern glacial cycling and southern climatic stability has long dominated explanations for how genetic diversity is distributed within temperate species in Eurasia and North America. However, growing evidence indicates the importance of cryptic refugia for northern colonization dynamics. An important geographic region to assess this is Fennoscandia, where recolonization at the end of the last glaciation was restricted to specific routes and temporal windows. We used genomic data to analyse genetic diversity and colonization history of the bank vole (Myodes glareolus) throughout Europe (>800 samples) with Fennoscandia as the northern apex. We inferred that bank voles colonized Fennoscandia multiple times by two different routes; with three separate colonizations via a southern land‐bridge route deriving from a “Carpathian” glacial refugium and one via a north‐eastern route from an “Eastern” glacial refugium near the Ural Mountains. Clustering of genome‐wide SNPs revealed high diversity in Fennoscandia, with eight genomic clusters: three of Carpathian origin and five Eastern. Time estimates revealed that the first of the Carpathian colonizations occurred before the Younger Dryas (YD), meaning that the first colonists survived the YD in Fennoscandia. Results also indicated that introgression between bank and northern red‐backed voles (Myodes rutilus) took place in Fennoscandia just after end‐glacial colonization. Therefore, multiple colonizations from the same and different cryptic refugia, temporal and spatial separations and interspecific introgression have shaped bank vole genetic variability in Fennoscandia. Together, these processes drive high genetic diversity at the apex of the northern expansion in this emerging model species.     

Niche differentiation in a postglacial colonizer,the bank vole Clethrionomys glareolus

Species‐level environmental niche modeling has been crucial in efforts to understand how species respond to climate variation and change. However, species often exhibit local adaptation and intraspecific niche differences that may be important to consider in predicting responses to climate. Here, we explore whether phylogeographic lineages of the bank vole originating from different glacial refugia (Carpathian, Western, Eastern, and Southern) show niche differentiation, which would suggest a role for local adaptation in biogeography of this widespread Eurasian small mammal. We first model the environmental requirements for the bank vole using species‐wide occurrences (210 filtered records) and then model each lineage separately to examine niche overlap and test for niche differentiation in geographic and environmental space. We then use the models to estimate past [Last Glacial Maximum (LGM) and mid‐Holocene] habitat suitability to compare with previously hypothesized glacial refugia for this species. Environmental niches are statistically significantly different from each other for all pairs of lineages in geographic and environmental space, and these differences cannot be explained by habitat availability within their respective ranges. Together with the inability of most of the lineages to correctly predict the distributions of other lineages, these results support intraspecific ecological differentiation in the bank vole. Model projections of habitat suitability during the LGM support glacial survival of the bank vole in the Mediterranean region and in central and western Europe. Niche differences between lineages and the resulting spatial segregation of habitat suitability suggest ecological differentiation has played a role in determining the present phylogeographic patterns in the bank vole. Our study illustrates that models pooling lineages within a species may obscure the potential for different responses to climate change among populations.     

Beyond the Mediterranean peninsulas: evidence of central European glacial refugia for a temperate forest mammal species, the bank vole (Clethrionomys glareolus)

This study details the phylogeographic pattern of the bank vole, Clethrionomys glareolus, a European rodent species strongly associated with forest habitat. We used sequences of 1011 base pairs of the mitochondrial DNA cytochrome b gene from 207 bank voles collected in 62 localities spread throughout its distribution area. Our results reveal the presence of three Mediterranean (Spanish, Italian and Balkan) and three continental (western, eastern and 'Ural') phylogroups. The endemic Mediterranean phylogroups did not contribute to the post-glacial recolonization of much of the Palaearctic range of species. Instead, the major part of this region was apparently recolonized by bank voles that survived in glacial refugia in central Europe. Moreover, our phylogeographic analyses also reveal differentiated populations of bank voles in the Ural mountains and elsewhere, which carry the mitochondrial DNA of another related vole species, the ruddy vole (Clethrionomys rutilus). In conclusion, this study demonstrates a complex phylogeographic history for a forest species in Europe which is sufficiently adaptable that, facing climate change, survives in relict southern and northern habitats. The high level of genetic diversity characterizing vole populations from parts of central Europe also highlights the importance of such regions as a source of intraspecific genetic biodiversity.     

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.     

Post‐glacial colonization of eastern Europe from the Carpathian refugium: evidence from mitochondrial DNA of the common vole Microtus arvalis

There is now considerable evidence for the survival of temperate species within glacial refugia that were situated at relatively high latitudes, notably the Carpathian Basin and Dordogne region in Europe. However, the prevalence of fossil remains in such locations is rarely matched by molecular evidence for their contribution to subsequent geographical and demographic expansion of the species in question. One obstacle to this has been insufficient analysis of modern samples from the relevant areas, in particular the parts of eastern Europe that surround the Carpathian refugium. In the present study, we examine the patterns of variation in mitochondrial DNA of the common vole (Microtus arvalis), obtained from existing museum specimens and from newly‐collected samples obtained in this area. We show that common voles from one of six extant mitochondrial DNA lineages have colonized most of the species' range in eastern Europe. We contend that the post‐glacial dispersal of this lineage most likely originated from the Carpathian refugium, adding support to the argument that such northern refugia made an important contribution to existing genetic diversity in Europe. © 2015 The Linnean Society of London, Biological Journal of the Linnean Society, 2015, ●● , ●●–●●.     

Evolutionary history of the bank vole Myodes glareolus: a morphometric perspective

The bank vole experienced a complex history during the Quaternary. Repeated isolation in glacial refugia led to the differentiation of several lineages in less than 300 000 years. We investigated if such a recent differentiation led to a significant divergence of phenotypic characters between European lineages, which might provide insight into processes of intraspecific differentiation. The size and shape of the first and third upper molars, and first lower molar, of bank voles genetically attributed to different lineages were quantified using an outline analysis of their occlusal surface. The three teeth present similar trends of decreasing size towards high latitudes. This trend, the inverse of Bergmann's rule, is interpreted as the result of a balance between metabolic efficiency and food availability, favouring small body size in cold regions. Molar shape appeared to differ between lineages despite genetic evidence of suture zones. A mosaic pattern of evolution between the different teeth was evidenced. The analysis of such phenotypic features appears as a valuable complement to genetic analyses, providing a complementary insight into evolutionary processes, such as selective pressures, that have driven the differentiation of the lineages. It may further allow the integration of the paleontological dimension of the bank vole phylogeographic history. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 681–694.     

Bank Voles in Southern Eurasia: Vicariance and Adaptation

Phylogeographic lineages are interpreted as the product of repeated isolation in glacial refugia, leading to vicariant differentiation. Being restricted to a given geographic area could also promote adaptive divergence in response to local conditions. The role of phylogeny and climate in the evolution of the bank vole (Myodes glareolus) was investigated here, focusing on molar tooth shape, a morphological feature related to the exploitation of food resources. A balanced role of phylogeny and climate was demonstrated. Response to environmental factors led to morphological convergence of bank voles from different lineages living in similar environments, and to within-lineage divergence in extreme environments. An important interaction of climate and phylogeny was found, suggesting that each lineage is living in a particular environment. This lineage-specific adaptation to a range of environmental conditions may have conditioned the potential of post-glacial recolonization of each lineage. Morphological covariation with environmental conditions further highlights the potential of adaptation of this species.     

Colonization history of north European field voles (Microtus agrestis) revealed by mitochondrial DNA

The genetic structure of field vole (Microtus agrestis) populations from northern Europe was examined by restriction fragment length polymorphisms of mitochondrial DNA (mtDNA) in 150 individuals from 67 localities. A total of 83 haplotypes was observed, most of which were rare and highly localized geographically. Overall nucleotide diversity was high (134%), but showed a tendency to decrease with higher latitude. Two major mtDNA lineages differing by 2% in nucleotide sequence were identified. A southern mtDNA lineage was observed in field voles from Britain, Denmark and southern and central Sweden, whereas voles from Finland and northern Sweden belonged to a northern lineage. The strict phylogeographic pattern suggests that the present population generic structure in field voles reflects glacial history: the two groups are derived from different glacial refugia, and recolonized Fennoscandia from two directions. A 150–200-km-wide secondary contact zone between the two mtDNA groups was found in northern Sweden. Distinct phylogeographic substructuring was observed within both major mtDNA groups.     

The effects of drift and selection on latitudinal genetic variation in Scandinavian common toads (Bufo bufo) following postglacial recolonisation

Clinal variation is paramount for understanding the factors shaping genetic diversity in space and time. During the last glacial maximum, northern Europe was covered by glacial ice that rendered the region uninhabitable for most taxa. Different evolutionary processes during and after the recolonisation of this area from different glacial refugia have affected the genetic landscape of the present day European flora and fauna. In this study, we focus on the common toad (Bufo bufo) in Sweden and present evidence suggesting that these processes have resulted in two separate lineages of common toad, which colonised Sweden from two directions. Using ddRAD sequencing data for demographic modelling, structure analyses, and analysis of molecular variance (AMOVA), we provide evidence of a contact zone located between Uppland and Västerbotten in central Sweden. Genetic diversity was significantly higher in southern Sweden compared to the north, in accordance with a pattern of decreased genetic diversity with increasing distance from glacial refugia. Candidate genes under putative selection are identified through outlier detection and gene–environment association methods. We provide evidence of divergent selection related to stress response and developmental processes in these candidate genes. The colonisation of Sweden by two separate lineages may have implications for how future conservation efforts should be directed by identifying management units and putative local adaptations.Subject terms: Evolutionary genetics, Genomics, Genetic variation     

Contrasting Patterns of Genetic Structuring in Natural Populations of Arabidopsis lyrata Subsp. petraea across Different Regions in Northern Europe

Level and partitioning of genetic diversity is expected to vary between contrasting habitats, reflecting differences in strength of ecological and evolutionary processes. Therefore, it is necessary to consider processes acting on different time scales when trying to explain diversity patterns in different parts of species'' distributions. To explore how historical and contemporary factors jointly may influence patterns of genetic diversity and population differentiation, we compared genetic composition in the perennial herb Arabidopsis lyrata ssp. petraea from the northernmost parts of its distribution range on Iceland to that previously documented in Scandinavia. Leaf tissue and soil were sampled from ten Icelandic populations of A. lyrata. Seedlings were grown from soil samples, and tissue from above-ground and seed bank individuals were genotyped with 21 microsatellite markers. Seed bank density in Icelandic populations was low but not significantly different from that observed in Norwegian populations. While within-population genetic diversity was relatively high on Iceland (H _E = 0.35), among-population differentiation was low (F _ST = 0.10) compared to Norwegian and Swedish populations. Population differentiation was positively associated with geographical distance in both Iceland and Scandinavia, but the strength of this relationship varied between regions. Although topography and a larger distribution range may explain the higher differentiation between mountainous Norwegian relative to lowland populations in Sweden, these factors cannot explain the lower differentiation in Icelandic compared to Swedish populations. We propose that low genetic differentiation among Icelandic populations is not caused by differences in connectivity, but is rather due to large historical effective population sizes. Thus, rather than contemporary processes, historical factors such as survival of Icelandic lineages in northern refugia during the last glacial period may have contributed to the observed pattern.     

Phylogeographic structure in long‐tailed voles (Rodentia: Arvicolinae) belies the complex Pleistocene history of isolation,divergence, and recolonization of Northwest North America's fauna

Quaternary climate fluctuations restructured biodiversity across North American high latitudes through repeated episodes of range contraction, population isolation and divergence, and subsequent expansion. Identifying how species responded to changing environmental conditions not only allows us to explore the mode and tempo of evolution in northern taxa, but also provides a basis for forecasting future biotic response across the highly variable topography of western North America. Using a multilocus approach under a Bayesian coalescent framework, we investigated the phylogeography of a wide‐ranging mammal, the long‐tailed vole, Microtus longicaudus. We focused on populations along the North Pacific Coast to refine our understanding of diversification by exploring the potentially compounding roles of multiple glacial refugia and more recent fragmentation of an extensive coastal archipelago. Through a combination of genetic data and species distribution models (SDMs), we found that historical climate variability influenced contemporary genetic structure, with multiple isolated locations of persistence (refugia) producing multiple divergent lineages (Beringian or northern, southeast Alaska or coastal, and southern or continental) during glacial advances. These vole lineages all occur along the North Pacific Coast where the confluence of numerous independent lineages in other species has produced overlapping zones of secondary contact, collectively a suture zone. Finally, we detected high levels of neoendemism due to complex island geography that developed in the last 10,000 years with the rising sea levels of the Holocene.     

Population genetic structure of the bank vole Myodes glareolus within its glacial refugium in peninsular Italy

It is now well established that Southern European peninsulas have been major glacial refugia for temperate species during Pleistocene climatic oscillations. However, substantial environmental changes occurred also within these peninsulas throughout the Pleistocene, raising questions about the role and interplay of various microevolutionary processes in shaping patterns of intraspecific diversity within these areas. Here, we investigate the patterns of genetic variation in the bank vole Myodes glareolus within the Italian peninsula. By using a panel of 13 microsatellite loci, we found more intraspecific variation than expected based on previous assessments. Indeed, both Bayesian and ordination‐based clustering analyses of variation recovered five main geographic/genetic clusters along the peninsula, with three clusters geographically restricted to the southern portion of the peninsula. This clustering is supported by previous evidences of some morphological distinctiveness among these populations. This pattern can be explained by a refugia‐within‐refugia scenario, with the occurrence of multiple sub‐refugia for the bank vole within the Italian peninsula, likely promoted by the major palaeo‐environmental changes which affected forested habitats within this area during the Pleistocene. Moreover, our results support a scenario whereby the high levels of intraspecific diversity observed within major Pleistocene refugia are better explained by dynamic microevolutionary processes occurred within these areas, rather than by long‐term demographic stability of refugial population. Finally, the narrow and isolated distribution of some of the identified lineages suggests the need for future assessments of their conservation and taxonomic status.     

Distribution and parameters of genetic polymorphism in northern red-backed vole (<Emphasis Type="Italic">Clethrionomys rutilus</Emphasis>) and bank vole (<Emphasis Type="Italic">Clethrionomys glareolus</Emphasis>) in West Siberia

This article presents data on the genetic variability of the northern red-backed vole and the bank vole that live sympatrically in West Siberia. The two species of voles have comparable, relatively high indices of genetic variability of inter simple sequences repeats DNA. The proportion of polymorphic DNA markers is 95–98%, and the Nei’s genetic diversity index is 0.33–0.35. A total of 47–58% of allozyme loci in the voles are polymorphic, and the average heterozygosity per locus is 0.058 in the northern red-backed vole and 0.054 in the bank vole. Interpopulation differentiation is less pronounced in the red-backed vole (F _ST 0.293) compared to the bank vole (F _ST 0.475). Individuals of the hybrid line of the bank vole with the mitochondrial haplotype of the red-backed vole have been found by PCR typing of cytochrome b gene fragment of mtDNA. The distribution boundary of the hybrid line of bank voles goes farther to the northeast than was shown in earlier works. The proportion of hybrid specimens range from 2 to 34%. The indices of genetic variability in the hybrid line of the bank vole are lower than those of the parental species.     

Genome‐wide set of SNPs reveals evidence for two glacial refugia and admixture from postglacial recolonization in an alpine ungulate

Past glaciation events have played a major role in shaping the genetic diversity and distribution of wild sheep in North America. The advancement of glaciers can isolate populations in ice‐free refugia, where they can survive until the recession of ice sheets. The major Beringian refugium is thought to have held thinhorn sheep (Ovis dalli) populations during times of glacial advance. While isolation in the major refugium can account for much of the genetic and morphological diversity seen in extant thinhorn sheep populations, mounting evidence suggests the persistence of populations in smaller minor refugia. We investigated the refugial origins of thinhorn sheep using ~10 000 SNPs obtained via a cross‐species application of the domestic sheep ovine HD BeadChip to genotype 52 thinhorn sheep and five bighorn sheep (O. canadensis) samples. Phylogenetic inference revealed a distinct lineage of thinhorn sheep inhabiting British Columbia, which is consistent with the survival of a group of thinhorn sheep in a minor refugium separate from the Beringian refugium. Isolation in separate glacial refugia probably mediated the evolution of the two thinhorn sheep subspecies, the white Dall's sheep (O. d. dalli), which persisted in Beringia, and the dark Stone's sheep (O. d. stonei), which utilized the minor refugium. We also found the first genetic evidence for admixture between sheep from different glacial refugia in south‐central Yukon as a consequence of post glacial expansion and recolonization. These results show that glaciation events can have a major role in the evolution of species inhabiting previously glaciated habitats and the need to look beyond established refugia when examining the evolutionary history of such species.     

The Celtic fringe of Britain: insights from small mammal phylogeography

Recent genetic studies have challenged the traditional view that the ancestors of British Celtic people spread from central Europe during the Iron Age and have suggested a much earlier origin for them as part of the human recolonization of Britain at the end of the last glaciation. Here we propose that small mammals provide an analogue to help resolve this controversy. Previous studies have shown that common shrews (Sorex araneus) with particular chromosomal characteristics and water voles (Arvicola terrestris) of a specific mitochondrial (mt) DNA lineage have peripheral western/northern distributions with striking similarities to that of Celtic people. We show that mtDNA lineages of three other small mammal species (bank vole Myodes glareolus, field vole Microtus agrestis and pygmy shrew Sorex minutus) also form a ‘Celtic fringe’. We argue that these small mammals most reasonably colonized Britain in a two-phase process following the last glacial maximum (LGM), with climatically driven partial replacement of the first colonists by the second colonists, leaving a peripheral geographical distribution for the first colonists. We suggest that these natural Celtic fringes provide insight into the same phenomenon in humans and support its origin in processes following the end of the LGM.     

New insights on postglacial colonization in western Europe: the phylogeography of the Leisler's bat (Nyctalus leisleri)

Despite recent advances in the understanding of the interplay between a dynamic physical environment and phylogeography in Europe, the origins of contemporary Irish biota remain uncertain. Current thinking is that Ireland was colonized post-glacially from southern European refugia, following the end of the last glacial maximum (LGM), some 20 000 years BP. The Leisler''s bat (Nyctalus leisleri), one of the few native Irish mammal species, is widely distributed throughout Europe but, with the exception of Ireland, is generally rare and considered vulnerable. We investigate the origins and phylogeographic relationships of Irish populations in relation to those across Europe, including the closely related species N. azoreum. We use a combination of approaches, including mitochondrial and nuclear DNA markers, in addition to approximate Bayesian computation and palaeo-climatic species distribution modelling. Molecular analyses revealed two distinct and diverse European mitochondrial DNA lineages, which probably diverged in separate glacial refugia. A western lineage, restricted to Ireland, Britain and the Azores, comprises Irish and British N. leisleri and N. azoreum specimens; an eastern lineage is distributed throughout mainland Europe. Palaeo-climatic projections indicate suitable habitats during the LGM, including known glacial refugia, in addition to potential novel cryptic refugia along the western fringe of Europe. These results may be applicable to populations of many species.     

A complex scenario of glacial survival in Mediterranean and continental refugia of a temperate continental vole species (Microtus arvalis) in Europe

The role of glacial refugia in shaping contemporary species distribution is a long-standing question in phylogeography and evolutionary ecology. Recent studies are questioning previous paradigms on glacial refugia and postglacial recolonization pathways in Europe, and more flexible phylogeographic scenarios have been proposed. We used the widespread common vole Microtus arvalis as a model to investigate the origin, locations of glacial refugia, and dispersal pathways, in the group of “Continental” species in Europe. We used a Bayesian spatiotemporal diffusion analysis (relaxed random walk model) of cytochrome b sequences across the species range, including newly collected individuals from 10 Iberian localities and published sequences from 68 localities across 22 European countries. Our data suggest that the species originated in Central Europe, and we revealed the location of multiple refugia (in both southern peninsulas and continental regions) for this continental model species. Our results confirm the monophyly of Iberian voles and the pre-LGM divergence between Iberian and European voles. We found evidence of restricted postglacial dispersal from refugia in Mediterranean peninsulas. We inferred a complex evolutionary and demographic history of M. arvalis in Europe over the last 50,000 years that does not adequately fit previous glacial refugial scenarios. The phylogeography of M. arvalis provides a paradigm of ice-age survival of a temperate continental species in western and eastern Mediterranean peninsulas (sources of endemism) and multiple continental regions (sources of postglacial spread). Our findings also provide support for a major role of large European river systems in shaping geographic boundaries of M. arvalis in Europe.     

Multiple cryptic refugia of forest grass Bromus benekenii in Europe as revealed by ISSR fingerprinting and species distribution modelling

Despite not having been fully recognized, the cryptic northern refugia of temperate forest vegetation in Central and Western Europe are one of the most important in the Holocene history of the vegetation on the subcontinent. We have studied a forest grass Bromus benekenii in 39 populations in Central, Western and Southern Europe with the use of PCR-ISSR fingerprinting. The indices of genetic population diversity, multivariate, and Bayesian analyses, supplemented with species distribution modelling have enabled at least three putative cryptic northern refugial areas to be recognized: in Western Europe—the Central and Rhenish Massifs, in Central Europe—the Bohemia–Moravia region and in the Eastern/Western Carpathians. Central Poland is the regional genetic melting-pot where several migratory routes might have met. Southern Poland had a different postglacial history and was under the influence of an Eastern/Western Carpathian cryptic refugium. More forest species should be checked in a west–east gradient in Europe to corroborate the hypothesis on the Western European glacial refugia.