Northern glacial refugia for the pygmy shrew Sorex minutus in Europe revealed by phylogeographic analyses and species distribution modelling

The southern European peninsulas (Iberian, Italian and Balkan) are traditionally recognized as glacial refugia from where many species colonized central and northern Europe after the Last Glacial Maximum (LGM). However, evidence that some species had more northerly refugia is accumulating from phylogeographic, palaeontological and palynological studies, and more recently from species distribution modelling (SDM), but further studies are needed to test the idea of northern refugia in Europe. Here, we take a rarely implemented multidisciplinary approach to assess if the pygmy shrew Sorex minutus, a widespread Eurasian mammal species, had northern refugia during the LGM, and if these influenced its postglacial geographic distribution. First, we evaluated the phylogeographic and population expansion patterns using mtDNA sequence data from 123 pygmy shrews. Then, we used SDM to predict present and past (LGM) potential distributions using two different training data sets, two different algorithms (Maxent and GARP) and climate reconstructions for the LGM with two different general circulation models. An LGM distribution in the southern peninsulas was predicted by the SDM approaches, in line with the occurrence of lineages of S. minutus in these areas. The phylogeographic analyses also indicated a widespread and strictly northern‐central European lineage, not derived from southern peninsulas, and with a postglacial population expansion signature. This was consistent with the SDM predictions of suitable LGM conditions for S. minutus occurring across central and eastern Europe, from unglaciated parts of the British Isles to much of the eastern European Plain. Hence, S. minutus likely persisted in parts of central and eastern Europe during the LGM, from where it colonized other northern areas during the late‐glacial and postglacial periods. Our results provide new insights into the glacial and postglacial colonization history of the European mammal fauna, notably supporting glacial refugia further north than traditionally recognized.     

Glacial refugia of mammals in Europe: evidence from fossil records

• 1 Glacial refugia were core areas for the survival of temperate species during unfavourable environmental conditions and were the sources of postglacial recolonizations. Unfortunately, the locations of glacial refugia of animals and plants are usually described by models, without reference to facts about real geographical ranges at that time.
• 2 Careful consideration of the faunal assemblages of archaeological sites from the Younger Palaeolithic, which are precisely dated to the Last Glacial Maximum (LGM), gives indications about the distribution of species during the LGM (23 000–16 000 bp ) and provides evidence for the locations of glacial refugia for mammalian species in Europe.
• 3 In Europe, 47 LGM sites, dating from 23 000 to 16 000 bp and containing typical temperate mammal species, have been described. The geographical range of these archaeological sites clearly shows a distribution which differs from the hypothesized traditional refuge areas of the temperate fauna. A considerable number of sites situated in the Dordogne in south‐western France and the Carpathian region contain records of red deer Cervus elaphus, roe deer Capreolus capreolus, wild boar Sus scrofa and red fox Vulpes vulpes.
• 4 The faunal composition of the majority of the evaluated Palaeolithic sites in the southern European peninsulas (with the exception of Greece), as well as France and the Carpathian region, indicates the co‐occurrence of these temperate species with cold‐adapted faunal elements such as mammoth Mammuthus primigenius and/or reindeer Rangifer tarandus.
• 5 The survival of species in Central European refugia would have significant consequences for phylogeography and would be revealed by the dominant distribution of haplotypes, originating from this region. A Carpathian refuge could also be the reason for the very early records of small mammals or mustelids from the Late‐Glacial or Interstadials before the LGM in regions like southern Germany.

     

Potential source and sink locations for climate‐driven species range shifts in Europe since the Last Glacial Maximum

Aim To identify potential source and sink locations for climate‐driven species range shifts in Europe since the Last Glacial Maximum (LGM). Location Europe. Methods We developed a new approach combining past‐climate simulations with the concept of analogous climate space. Our index gives a continuous measure of the potential of a location to have acted as a source or a sink for species that have shifted their ranges since the LGM. High glacial source potential is indicated by LGM climatic conditions that are widespread now; high post‐glacial sink potential is indicated by current climatic conditions that were widespread at the LGM. The degree of isolation of source and sink areas was calculated as the median distance to areas with analogous climate conditions. Results We identified areas of high glacial source potential in the previously recognized refugial areas in the southern European peninsulas, but also in large areas in central‐western Europe. The most climatically isolated source areas were located in northern Spain, in north‐western Europe and in eastern Turkey. From here species would have had to cover substantial distances to find current climate conditions analogous to LGM conditions of these areas. Areas with high post‐glacial sink potential were mainly located in Fennoscandia and in central and south‐eastern Europe. Some of the most isolated sink areas were located in the Spanish highlands and around the Baltic Sea. Main conclusions Our species‐independent approach successfully identified previously recognized glacial refugial areas with high source potential for species range shifts in southern Europe and in addition highlighted other potential source areas in central Europe. This study offers new insights into how the distribution of past and current climatic conditions may have influenced past species range shifts and current large‐scale biodiversity patterns.     

Ancient DNA supports southern survival of Richardson's collared lemming (Dicrostonyx richardsoni) during the last glacial maximum

Collared lemmings (genus Dicrostonyx) are circumpolar Arctic arvicoline rodents associated with tundra. However, during the last glacial maximum (LGM), Dicrostonyx lived along the southern ice margin of the Laurentide ice sheet in communities comprising both temperate and boreal species. To better understand these communities and the fate of these southern individuals, we compare mitochondrial cytochrome b sequence data from three LGM‐age Dicrostonyx fossils from south of the Laurentide ice sheet to sequences from modern Dicrostonyx sampled from across their present‐day range. We test whether the Dicrostonyx populations from LGM‐age continental USA became extinct at the Pleistocene–Holocene transition ~11000 years ago or, alternatively, if they belong to an extant species whose habitat preferences can be used to infer the palaeoclimate along the glacial margin. Our results indicate that LGM‐age Dicrostonyx from Iowa and South Dakota belong to Dicrostonyx richardsoni, which currently lives in a temperate tundra environment west of Hudson Bay, Canada. This suggests a palaeoclimate south of the Laurentide ice sheet that contains elements similar to the more temperate shrub tundra characteristic of extant D. richardsoni habitat, rather than the very cold, dry tundra of the Northern Arctic. While more data are required to determine whether or not the LGM southern population is ancestral to extant D. richardsoni, it seems most probable that the species survived the LGM in a southern refugium.     

Apparently contradictory routes in the expansion of two fish species in the Eastern Atlantic

One specimen of the boreal Clupea harengus (Clupeidae) was caught southerly, in the Gulf of Cádiz (southern Spain), whereas a specimen of the subtropical Epinephelus aeneus (Serranidae) was caught northerly, in Galician waters (north-west of Spain). These catches represent a new southernmost record in the eastern Atlantic for C. harengus and the second northernmost for E. aeneus. The co-occurrence of warm and cold affinity fishes in temperate latitudes is unusual, but the existence of anomalies in the temperature of the Atlantic waters could be behind these findings.     

Past and future range shifts and loss of diversity in dwarf willow (Salix herbacea L.) inferred from genetics, fossils and modelling

Aim Climate change may cause loss of genetic diversity. Here we explore how a multidisciplinary approach can be used to infer effects of past climate change on species distribution and genetic diversity and also to predict loss of diversity due to future climate change. We use the arctic‐alpine plant Salix herbacea L. as a model. Location Europe, Greenland and eastern North America. Methods We analysed 399 samples from 41 populations for amplified fragment length polymorphism (AFLP) to identify current patterns of genetic structure and diversity and likely historical dispersal routes. Macrofossil records were compiled to infer past distribution, and species distribution models were used to predict the Last Glacial Maximum (LGM) and future distribution of climatically suitable areas. Results We found strong genetic differentiation between the populations from Europe/East Greenland and those from Canada/West Greenland, indicating a split probably predating the LGM. Much less differentiation was observed among the four genetic groups identified in Europe, and diversity was high in the Scandinavian as well as in southern alpine populations. Continuous distribution in Central Europe during the last glaciation was inferred based on the fossil records and distribution modelling. A 46–57% reduction in suitable areas was predicted in 2080 compared to present. However, mainly southern alpine populations may go extinct, causing a loss of about 5% of the genetic diversity in the species. Main conclusions From a continuous range in Central Europe during the last glaciation, northward colonization probably occurred as a broad front maintaining diversity as the climate warmed. This explains why potential extinction of southern populations by 2080 will cause a comparatively low loss of the genetic diversity in S. herbacea. For other species with different glacial histories, however, the expected climate‐change induced regional extinction may cause a more severe loss of genetic diversity. We conclude that our multidisciplinary approach may be a useful tool for assessing impact of climate change on loss of genetic diversity.     

Cyclic voles and shrews and non-cyclic mice in a marginal grassland within European temperate forest

Cyclic population dynamics of small mammals are not restricted to the boreal and arctic zones of Eurasia and North America, but long-term data series from lower latitudes are still less common. We demonstrated here the presence of periodic oscillations in small mammal populations in eastern Poland using 22-year (1986–2007) trapping data from marginal meadow and river valley grasslands located in the extensive temperate woodland of Białowieża Primeval Forest. The two most common species inhabiting meadows and river valleys, root vole Microtus oeconomus and common shrew Sorex araneus, exhibited synchronous periodic changes, characterised by a 3-year time lag as indicated by an autocorrelation function. Moreover, the cycles of these two species were synchronous within both habitats. Population dynamics of the striped field mouse Apodemus agrarius was not cyclic. However, this species regularly reached maximum density 1 year before the synchronized peak of root voles and common shrews, which may suggest the existence of interspecific competition. Dynamics of all three species was dominated by direct density-dependent process, whereas delayed density dependent feedback was significant only in the root vole and common shrew. Climatic factors acting in winter and spring (affecting mainly survival and initial reproduction rates) were more important than those acting in summer and autumn and affected significantly only the common shrew. High temperatures in winter and spring had positive effects on autumn-to-autumn changes in abundance of this species, whereas deep snow in combination with high rainfall in spring negatively affected population increase rates in common shrew.     

Natural infections of small mammals with blood parasites on the borderland of boreal and temperate forest zones

Blood parasites of small mammals living in Białowieża Forest (eastern Poland) were investigated between 1996 and 2002. The following haemoparasite species were found:Trypanosoma (Herpetosoma) evotomys in bank voleClethrionomys glareolus; T. (H.) microti in root voleMicrotus oeconomus; Babesia microti in root vole;Hepatozoon erhardovae in bank vole andHepatozoon sp. in root vole. Some non-identifiedBartonella species were found in bank vole, root vole, field voleMicrotus agrestis, yellow-necked mouseApodemus flavicollis, common shrewSorex araneus, Eurasian water shrewNeomys fodiens, and Mediterranean water shrewN. anomalus. The prevalence and diversity of blood parasites were lower in shrews than small rodents. Totally, 52.0% of bank voles, 50.0% of root voles, 32.5% of common shrews, and 41.2% of Eurasian water shrews were infected with any of the blood parasites. Mixed infections were seldom observed in bank vole (17.3% of investigated individuals) and root vole (14.7%). No animals were infected with three or four parasites simultaneously. Infection of Białowieża small mammals with haemoparasites seemed to be similar to those described in other temperate forest regions rather than boreal ones. Infection rates of rodent species seem to be higher in their typical habitats: for bank vole it was the highest in mixed forest, whereas for root vole in sedge swamp. The results suggest that Arvicolidae play a greater role than Muridae or Soricidae in maintenance ofBabesia andHepatozoon foci in natural environments of central Europe.     

Quantifying components of risk for European woody species under climate change

Estimates of species extinction risk under climate change are generally based on differences in present and future climatically suitable areas. However, the locations of potentially suitable future environments (affecting establishment success), and the degree of climatic suitability in already occupied and new locations (affecting population viability) may be equally important determinants of risk. A species considered to be at low risk because its future distribution is predicted to be large, may actually be at high risk if these areas are out of reach, given the species' dispersal and migration rates or if all future suitable locations are only marginally suitable and the species is unlikely to build viable populations in competition with other species. Using bioclimatic models of 17 representative European woody species, we expand on current ways of risk assessment and suggest additional measures based on (a) the distance between presently occupied areas and areas predicted to be climatically suitable in the future and (b) the degree of change in climatic suitability in presently occupied and unoccupied locations. Species of boreal and temperate deciduous forests are predicted to face higher risk from loss of climatically suitable area than species from warmer and drier parts of Europe by 2095 using both the moderate B1 and the severe A1FI emission scenario. However, the average distance from currently occupied locations to areas predicted suitable in the future is generally shorter for boreal species than for southern species. Areas currently occupied will become more suitable for boreal and temperate species than for Mediterranean species whereas new suitable areas outside a species' current range are expected to show greater increases in suitability for Mediterranean species than for boreal and temperate species. Such additional risk measures can be easily derived and should give a more comprehensive picture of the risk species are likely to face under climate change.     

Investigating past range dynamics for a weed of cultivation,Silene vulgaris

Since the last glacial maximum (LGM), many plant and animal taxa have expanded their ranges by migration from glacial refugia. Weeds of cultivation may have followed this trend or spread globally following the expansion of agriculture or ruderal habitats associated with human‐mediated disturbance. We tested whether the range expansion of the weed Silene vulgaris across Europe fit the classical model of postglacial expansion from southern refugia, or followed known routes of the expansion of human agricultural practices. We used species distribution modeling to predict spatial patterns of postglacial expansion and contrasted these with the patterns of human agricultural expansion. A population genetic analysis using microsatellite loci was then used to test which scenario was better supported by spatial patterns of genetic diversity and structure. Genetic diversity was highest in southern Europe and declined with increasing latitude. Locations of ancestral demes from genetic cluster analysis were consistent with areas of predicted refugia. Species distribution models showed the most suitable habitat in the LGM on the southern coasts of Europe. These results support the typical postglacial northward colonization from southern refugia while refuting the east‐to‐west agricultural spread as the main mode of expansion for S. vulgaris. We know that S. vulgaris has recently colonized many regions (including North America and other continents) through human‐mediated dispersal, but there is no evidence for a direct link between the Neolithic expansion of agriculture and current patterns of genetic diversity of S. vulgaris in Europe. Therefore, the history of range expansion of S. vulgaris likely began with postglacial expansion after the LGM, followed by more recent global dispersal by humans.     

Species distribution modeling and molecular markers suggest longitudinal range shifts and cryptic northern refugia of the typical calcareous grassland species Hippocrepis comosa (horseshoe vetch)

Calcareous grasslands belong to the most diverse, endangered habitats in Europe, but there is still insufficient information about the origin of the plant species related to these grasslands. In order to illuminate this question, we chose for our study the representative grassland species Hippocrepis comosa (Horseshoe vetch). Based on species distribution modeling and molecular markers, we identified the glacial refugia and the postglacial migration routes of the species to Central Europe. We clearly demonstrate that H. comosa followed a latitudinal and due to its oceanity also a longitudinal gradient during the last glacial maximum (LGM), restricting the species to southern refugia situated on the Peninsulas of Iberia, the Balkans, and Italy during the last glaciation. However, we also found evidence for cryptic northern refugia in the UK, the Alps, and Central Germany. Both species distribution modeling and molecular markers underline that refugia of temperate, oceanic species such as H. comosa must not be exclusively located in southern but also in western of parts of Europe. The analysis showed a distinct separation of the southern refugia into a western cluster embracing Iberia and an eastern group including the Balkans and Italy, which determined the postglacial recolonization of Central Europe. At the end of the LGM, H. comosa seems to have expanded from the Iberian refugium, to Central and Northern Europe, including the UK, Belgium, and Germany.     

Multivariate morphometric analysis of the Potamogeton compressus group (Potamogetonaceae)

The Potamogeton compressus group is a complex of three to five closely related species with a circumpolar distribution in the Northern Hemisphere. Multivariate morphometric analyses (principal component analysis, cluster analysis, canonical and classificatory discriminant analyses) were used to elucidate the patterns of variation within this group and to test the morphological differentiation of the species recognized in the current literature. From the entire distribution range, 156 specimens of the group were included in the numerical methods. Results from morphological comparison are discussed in relation to molecular data, reproductive behaviour and geographical distribution. Morphometric analyses provided evidence that this complex can be clearly divided into three groups, one of which was subdivided mainly on the basis of allopatric occurrence and genetic differentiation. These groups correspond to four species accepted here: P. acutifolius (temperate regions of Europe), P. compressus (boreal and temperate regions of Europe and Asia), P. manchuriensis (northeastern China and Russian Far East) and P. zosteriformis (boreal and temperate regions of North America). Two species, P. acutifolius and P. compressus, are partly sympatric, but clearly differentiated morphologically and genetically, and effectively isolated reproductively. Endemic P. manchuriensis is characterized by a unique combination of characters and an occurrence in a limited geographical area. Allopatric P. zosteriformis is weakly differentiated morphologically from P. compressus, but differs markedly in molecular markers correlated with geographical differentiation. It may represent a cryptic species. In contrast, a recently suggested concept of southern Siberian P. henningii was not supported by our analyses. Plants so named are considered here as slender phenotypes of the widespread and variable P. compressus. © 2012 The Linnean Society of London, Botanical Journal of the Linnean Society, 2012, 170 , 112–130.     

Phylogeographical structure of the boreal‐montane orchid Malaxis monophyllos as a result of multi‐directional gene flow

We investigated the phylogeographical structure of the boreal‐montane orchid Malaxis monophyllos in its Eurasian geographical range. We analysed four sequences of plastid DNA (trnL, trnL–trnF, rps16 and accD‐psaI), resulting in 19 haplotypes and revealing a high level of intraspecific diversity (H_D = 0.702 and π = 0.196 × 10⁻²), but showing a lack of phylogeographical structure. This pattern might be caused by multiple phenomena and processes, e.g. broad‐fronted recolonization with accompanying multi‐directional gene flow between populations and expansion from at least two refugial areas. Despite the lack of phylogeographical structure, three centres of haplotype diversity were indicated in the European part of the range of M. monophyllos. According to these data, alpine and lowland glacial refugia located between the ice sheets in the European Alps and the Scandinavian glaciers seem most likely to be in Europe. Moreover, models of climatically suitable areas during the Last Glacial Maximum (LGM) confirmed the Alps as a possible refuge, and indicated an opportunity for the persistence of M. monophyllos populations in Beringia and parts of Siberia. Using two models [Model for Interdisciplinary Research on Climate (MIROC) and Community Climate System Model (CCSM)], we predicted a significant reduction in climatically suitable areas for M. monophyllos in the future (2080). Our study also demonstrated that the biological features of M. monophyllos, including breeding system and dispersal mode, seem to be crucial in understanding its phylogeographical pattern. Our results also highlighted the importance of anthropogenic habitats as reservoirs of genetic diversity and alternative habitats for this species in the context of declining natural populations. © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 178 , 138–154.     

Phylogeography of European moose (Alces alces) based on contemporary mtDNA data and archaeological records

Phylogeography can help to determine LGM refugia and postglacial migration routes. However, the locations of LGM refugial areas in eastern Europe are not clear. Moose (Alces alces) is presently a common species in central and north-eastern Europe, but there are no studies showing its phylogenetic pattern and genetic diversity across its whole continuous range. Moose never became extinct in the eastern part of its range, and the eastern mtDNA lineage has the largest effective population size. The present study shows the phylogeographic pattern and genetic diversity of European moose and compares the results of mtDNA analyses with the archaeological record of the species to identify its LGM refugia and postglacial migration routes. I combined the mtDNA control region sequences obtained in all studies of moose in Europe and western Asia. The genetic data were then compared with the archaeological records of the species dated to the LGM. I found that the European moose lineage inhabits Europe and western Asia. It is composed of two clades: the eastern and the central-western, consisting of a total of six discrete haplogroups. The most complex, the eastern clade, has the largest range. Some of the haplogroups have narrow or scattered distributions and two are common in almost the whole range. Genetic diversity hotspots were detected in contact zones of different mtDNA haplogroups rather than in the LGM refugial areas of moose. Archaeological records dated to the LGM were found in several localities in central, southern and eastern Europe as well as in western Asia. The range of the moose during the LGM was much larger than previously thought. The eastern clade survived the LGM in western Siberia, the Ural Mountains and Russian plain. LGM refugia of moose were also located in the Caucasus, Carpathians, Balkans and northern Italy.     

Late‐glacial recolonization and phylogeography of European red deer (Cervus elaphus L.)

The Pleistocene was an epoch of extreme climatic and environmental changes. How individual species responded to the repeated cycles of warm and cold stages is a major topic of debate. For the European fauna and flora, an expansion–contraction model has been suggested, whereby temperate species were restricted to southern refugia during glacial times and expanded northwards during interglacials, including the present interglacial (Holocene). Here, we test this model on the red deer (Cervus elaphus) a large and highly mobile herbivore, using both modern and ancient mitochondrial DNA from the entire European range of the species over the last c. 40 000 years. Our results indicate that this species was sensitive to the effects of climate change. Prior to the Last Glacial Maximum (LGM) haplogroups restricted today to South‐East Europe and Western Asia reached as far west as the UK. During the LGM, red deer was mainly restricted to southern refugia, in Iberia, the Balkans and possibly in Italy and South‐Western Asia. At the end of the LGM, red deer expanded from the Iberian refugium, to Central and Northern Europe, including the UK, Belgium, Scandinavia, Germany, Poland and Belarus. Ancient DNA data cannot rule out refugial survival of red deer in North‐West Europe through the LGM. Had such deer survived, though, they were replaced by deer migrating from Iberia at the end of the glacial. The Balkans served as a separate LGM refugium and were probably connected to Western Asia with genetic exchange between the two areas.     

Evolutionary history of the scorpionfly Dicerapanorpa magna (Mecoptera,Panorpidae)

Climate changes can have fundamental impacts on the distributional patterns of montane species, and range shifts frequently lead to allopatric divergence followed by the establishment of secondary contact zones. Many European and North American organisms have retreated to southern refugia during glacial periods and colonized northward during postglacial periods, but little is known about the evolutionary response of cold‐adapted insects to Pleistocene climate changes in eastern Asia. The scorpionfly Dicerapanorpa magna (Chou), with cold temperate habitat preference and weak dispersal ability, provides a good model system to explore how climate changes have influenced the distribution and divergence of cold‐adapted insects in eastern Asia. This study reconstructed the demographic dynamics and evolutionary history of D. magna with phylogeographic approaches, and predicted the species’ suitable areas under the Last Glacial Maximum (LGM) and current scenarios with the ecological niche modelling analysis. The mitochondrial cytochrome c oxidase subunit I resolved three phylogenetic lineages in D. magna dating back to Pleistocene, corresponding well with the geographically isolated Qinling, Bashan and Minshan Mountains. The ecological niche modelling recovered the suitable habitats for D. magna were the Qinling and Bashan Mountains under LGM and current conditions. The three lineages of D. magna might be in a process of incipient speciation, and likely derived their current distribution from separate glacial origins, followed by vicariance and divergence.     

Postglacial colonization of Europe by the barbastelle bat: agreement between molecular data and past predictive modelling

The barbastelle (Barbastella barbastellus) is a rare forest bat with a wide distribution in Europe. Here, we combine results from the analysis of two mtDNA fragments with species distribution modelling to determine glacial refugia and postglacial colonization routes. We also investigated whether niche conservatism occurs in this species. Glacial refugia were identified in the three southern European peninsulas: Iberia, Italy and the Balkans. These latter two refugia played a major role in the postglacial colonization process, with their populations expanding to England and central Europe, respectively. Palaeo‐distribution models predicted that suitable climatic conditions existed in the inferred refugia during the last glacial maximum (LGM). Nevertheless, the overlap between the current and the LGM distributions was almost inexistent in Italy and in the Balkans, meaning that B. barbastellus populations were forced to shift range between glacial and interglacial periods, a process that probably caused some local extinctions. In contrast, Iberian populations showed a ‘refugia within refugium’ pattern, with two unconnected areas containing stable populations (populations that subsisted during both glacial and interglacial phases). Moreover, the match between LGM models and the refugial areas determined by molecular analysis supported the hypothesis of niche conservatism in B. barbastellus. We argue that geographic patterns of genetic structuring, altogether with the modelling results, indicate the existence of four management units for conservation: Morocco, Iberia, Italy and UK, and Balkans and central Europe. In addition, all countries sampled possessed unique gene pools, thus stressing the need for the conservation of local populations.     

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.     

Spatial patterns of genetic diversity across European subspecies of the mountain hare, Lepus timidus L

Fossil evidence shows that populations of species that currently inhabit arctic and boreal regions were not isolated in refugia during glacial periods, but instead maintained populations across large areas of central Europe. These species commonly display little reduction in genetic diversity in northern areas of their range, in contrast to many temperate species. The mountain hare currently inhabits both temperate and arctic-boreal regions. We used nuclear microsatellite and mtDNA sequence data to examine population structure and alternate phylogeographic hypotheses for the mountain hare, that is, temperate type (lower genetic diversity in northern areas) and arctic-boreal type (high northern genetic diversity). Both data sets revealed concordant patterns. Highest allelic richness, expected heterozygosity and mtDNA haplotype diversity were identified in the most northerly subspecies, indicating that this species more closely maps to phylogeographic patterns observed in arctic-boreal rather than temperate species. With regard to population structure, the Alpine and Fennoscandian subspecies were most genetically similar (F(ST) approximately 0.1). These subspecies also clustered together on the mtDNA tree and were assigned with highest likelihood to a common Bayesian cluster. This is consistent with fossil evidence for intermediate populations in the central European plain, persisting well into the postglacial period. In contrast, the geographically close Scottish and Irish populations occupied separate Bayesian clusters, distinct clades on the mtDNA maximum likelihood tree and were genetically divergent from each other (F(ST) > 0.4) indicating the influence of genetic drift, long isolation (possibly dating from the late glacial era) and/or separate postglacial colonisation routes.