Phylogeography of western Palaearctic reptiles – Spatial and temporal speciation patterns

A phylogeographic analysis of eight species complexes of European reptiles was performed using different molecular methods. While mitochondrial genes (mainly cytochrome b sequences) enabled conclusions about phylogeography and differentiation, additional application of bisexually inherited markers provided information about speciation stages. As species with similar distribution patterns in southern and Central Europe were selected, matching phylogeographic patterns are useful for drawing general conclusions:

(1) The species complexes are in different stages of speciation. In some cases, cryptic species were detected.

(2) Highest genetic diversity occurs in southern Europe, the Near East and the Caucasus, regions corresponding with glacial refuges in the Iberian, Apennine and Balkan Peninsulas as well as in Turkey and the Caucasus. Often, several microrefugia must have existed in close neighbourhood. Additional microrefugia were located in southern France and in the Carpathian Basin.

(3) North Africa and the Middle East did not serve as glacial refuges for Central or northern European lineages and are typically inhabited by independent clades.

(4) Evidence for multiple range retractions and expansions, which were postulated for the times of Pleistocene climatic oscillations, could be found in the Balkans, but in Central Europe their traces have been wiped out by the last glacial. Only the Holocene invasion has left imprints in the genomes from this area.

(5) Central and northern Europe were recolonized from Balkan and Pontic refugia in the Holocene.

(6) Groups from the Iberian and Apennine Peninsulas rarely conquered other regions. This limitation can be attributed to the barrier function of the Pyrenees and the Alps.

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.     

Maternal genomic variability of the wild boar (Sus scrofa) reveals the uniqueness of East‐Caucasian and Central Italian populations

The phylogeography of the European wild boar was mainly determined by postglacial recolonization patterns from Mediterranean refugia after the last ice age. Here we present the first analysis of SNP polymorphism within the complete mtDNA genome of West Russian (n = 8), European (n = 64), and North African (n = 5) wild boar. Our analyses provided evidence of unique lineages in the East‐Caucasian (Dagestan) region and in Central Italy. A phylogenetic analysis revealed that these lineages are basal to the other European mtDNA sequences. We also show close connection between the Western Siberian and Eastern European populations. Also, the North African samples were clustered with the Iberian population. Phylogenetic trees and migration modeling revealed a high proximity of Dagestan sequences to those of Central Italy and suggested possible gene flow between Western Asia and Southern Europe which was not directly related to Northern and Central European lineages. Our results support the presence of old maternal lineages in two Southern glacial refugia (i.e., Caucasus and the Italian peninsula), as a legacy of an ancient wave of colonization of Southern Europe from an Eastern origin.     

Positioning the Red Deer (Cervus elaphus) Hunted by the Tyrolean Iceman into a Mitochondrial DNA Phylogeny

In the last years several phylogeographic studies of both extant and extinct red deer populations have been conducted. Three distinct mitochondrial lineages (western, eastern and North-African/Sardinian) have been identified reflecting different glacial refugia and postglacial recolonisation processes. However, little is known about the genetics of the Alpine populations and no mitochondrial DNA sequences from Alpine archaeological specimens are available. Here we provide the first mitochondrial sequences of an Alpine Copper Age Cervus elaphus. DNA was extracted from hair shafts which were part of the remains of the clothes of the glacier mummy known as the Tyrolean Iceman or Ötzi (5,350–5,100 years before present). A 2,297 base pairs long fragment was sequenced using a mixed sequencing procedure based on PCR amplifications and 454 sequencing of pooled amplification products. We analyzed the phylogenetic relationships of the Alpine Copper Age red deer''s haplotype with haplotypes of modern and ancient European red deer. The phylogenetic analyses showed that the haplotype of the Alpine Copper Age red deer falls within the western European mitochondrial lineage in contrast with the current populations from the Italian Alps belonging to the eastern lineage. We also discussed the phylogenetic relationships of the Alpine Copper Age red deer with the populations from Mesola Wood (northern Italy) and Sardinia.     

Genetic diversity and phylogeography of the European roe deer: the refuge area theory revisited

The extant taxa of central and northern Europe are commonly believed to derive from Pleistocene ancestors, who moved to the north from three separate glacial refugia: the Iberian and Italian peninsulae, as well as the southern Balkans. The issue of postglacial dispersal patterns was addressed through the investigation of population structure and phylogeography of the European roe deer, Capreolus capreolus . The genetic diversity in 376 individuals representing 14 allegedly native populations across their European range was assessed, using ten autosomal microsatellite loci and restriction fragment length polymorphisms of the mitochondrial D-loop and NADH dehydrogenase 1 gene segments. Our results suggest the existence of three major genetic lineages of roe deer in Europe. One comprises populations in the south-western limit of the species' distribution (i.e. Iberia), where an internal substructure splits a northern from a southern sublineage. A second lineage includes populations of southern and eastern Europe, as well as a separate sublineage sampled in central-southern Italy, where the existence of the subspecies Capreolus c. italicus was supported. In central-northern Europe, a third lineage is present, which appeared genetically rather homogeneous, although admixed, and equally divergent from both the eastern and western lineages. Current patterns of intraspecific genetic variation suggest that postglacial recolonization routes of this cervid to northern Europe could be due to range expansion from one or more refugia in central-eastern Europe, rather than proceeding from the Mediterranean areas.  © 2006 The Linnean Society of London, Biological Journal of the Linnean Society , 2006, 88 , 85–100.     

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.     

Mitochondrial Haplogroup U5b3: A Distant Echo of the Epipaleolithic in Italy and the Legacy of the Early Sardinians

There are extensive data indicating that some glacial refuge zones of southern Europe (Franco-Cantabria, Balkans, and Ukraine) were major genetic sources for the human recolonization of the continent at the beginning of the Holocene. Intriguingly, there is no genetic evidence that the refuge area located in the Italian Peninsula contributed to this process. Here we show, through phylogeographic analyses of mitochondrial DNA (mtDNA) variation performed at the highest level of molecular resolution (52 entire mitochondrial genomes), that the most likely homeland for U5b3—a haplogroup present at a very low frequency across Europe—was the Italian Peninsula. In contrast to mtDNA haplogroups that expanded from other refugia, the Holocene expansion of haplogroup U5b3 toward the North was restricted by the Alps and occurred only along the Mediterranean coasts, mainly toward nearby Provence (southern France). From there, ∼7,000–9,000 years ago, a subclade of this haplogroup moved to Sardinia, possibly as a result of the obsidian trade that linked the two regions, leaving a distinctive signature in the modern people of the island. This scenario strikingly matches the age, distribution, and postulated geographic source of a Sardinian Y chromosome haplogroup (I2a2-M26), a paradigmatic case in the European context of a founder event marking both female and male lineages.     

Molecular analysis of the Pleistocene history of Saxifraga oppositifolia in the Alps

A recent circumpolar survey of chloroplast DNA (cpDNA) haplotypes identified Pleistocene glacial refugia for the Arctic-Alpine Saxifraga oppositifolia in the Arctic and, potentially, at more southern latitudes. However, evidence for glacial refugia within the ice sheet covering northern Europe during the last glacial period was not detected either with cpDNA or in another study of S. oppositifolia that surveyed random amplified polymorphic DNA (RAPD) variation. If any genotypes survived in such refugia, they must have been swamped by massive postglacial immigration of periglacial genotypes. The present study tested whether it is possible to reconstruct the Pleistocene history of S. oppositifolia in the European Alps using molecular methods. Restriction fragment length polymorphism (RFLP) analysis of cpDNA of S. oppositifolia, partly sampled from potential nunatak areas, detected two common European haplotypes throughout the Alps, while three populations harboured two additional, rare haplotypes. RAPD analysis confirmed the results of former studies on S. oppositifolia; high within, but low among population genetic variation and no particular geographical patterning. Some Alpine populations were not perfectly nested in this common gene pool and contained private RAPD markers, high molecular variance or rare cpDNA haplotypes, indicating that the species could possibly have survived on ice-free mountain tops (nunataks) in some parts of the Alps during the last glaciation. However, the overall lack of a geographical genetic pattern suggests that there was massive immigration of cpDNA and RAPD genotypes by seed and pollen flow during postglacial times. Thus, the glacial history of S. oppositifolia in the Alps appears to resemble closely that suggested previously for the species in northern Europe.     

Molecular evidence for glacial refugia of mountain plants in the European Alps

Many mountain ranges have been strongly glaciated during the Quaternary ice ages, and the locations of glacial refugia of mountain plants have been debated for a long time. A series of detailed molecular studies, investigating intraspecific genetic variation of mountain plants in the European Alps, now allows for a first synopsis. A comparison of the phylogeographic patterns with geological and palaeoenvironmental data demonstrates that glacial refugia were located along the southwestern, southern, eastern and northern border of the Alps. Additional glacial refugia were present in central Alpine areas, where high-elevation plants survived the last glaciation on ice-free mountain tops. The observed intraspecific phylogeographies suggest general patterns of glacial survival, which conform to well-known centres of Alpine species diversity and endemism. This implies that evolutionary or biogeographic processes induced by climatic fluctuations act on gene and species diversity in a similar way.     

Signature of mid‐Pleistocene lineages in the European silver fir (Abies alba Mill.) at its geographic distribution margin

In a conservation and sustainable management perspective, we identify the ecological, climatic, and demographic factors responsible for the genetic diversity patterns of the European silver fir (Abies alba Mill.) at its southwestern range margin (Pyrenees Mountains, France, Europe). We sampled 45 populations throughout the French Pyrenees and eight neighboring reference populations in the Massif Central, Alps, and Corsica. We genotyped 1,620 individuals at three chloroplast and ten nuclear microsatellite loci. We analyzed within‐ and among‐population genetic diversity using phylogeographic reconstructions, tests of isolation‐by‐distance, Bayesian population structure inference, modeling of demographic scenarios, and regression analyses of genetic variables with current and past environmental variables. Genetic diversity decreased from east to west suggesting isolation‐by‐distance from the Alps to the Pyrenees and from the Eastern to the Western Pyrenees. We identified two Pyrenean lineages that diverged from a third Alpine–Corsica–Massif Central lineage 0.8 to 1.1 M years ago and subsequently formed a secondary contact zone in the Central Pyrenees. Population sizes underwent contrasted changes, with a contraction in the west and an expansion in the east. Glacial climate affected the genetic composition of the populations, with the western genetic cluster only observed in locations corresponding to the coldest past climate and highest elevations. The eastern cluster was observed over a larger range of temperatures and elevations. All demographic events shaping the current spatial structure of genetic diversity took place during the Mid‐Pleistocene Transition, long before the onset of the Holocene. The Western Pyrenees lineage may require additional conservation efforts, whereas the eastern lineage is well protected in in situ gene conservation units. Due to past climate oscillations and the likely emergence of independent refugia, east–west oriented mountain ranges may be important reservoir of genetic diversity in a context of past and ongoing climate change in Europe.     

Molecular phylogeography of the viperine snake Natrix maura (Serpentes: Colubridae): Evidence for strong intraspecific differentiation

The molecular phylogeography of the viperine snake, Natrix maura (Linnaeus, 1758), was investigated using complete sequences of the mitochondrial cytochrome b gene and genomic ISSR-PCR fingerprinting. In a total of 120 samples, 44 unique cytochrome b haplotypes were found which defined three major genetic lineages associated with samples from Morocco, Tunisia and Europe, respectively. The same lineages were supported by nuclear data. A possible fourth lineage exists in southern Spain. Genetic distances of cytochrome b sequences between the three main lineages were in the range of 3.9–5.6%, suggesting independent evolution since the early Pliocene. Distinction of the three major lineages at the subspecies or species level is discussed to account taxonomically for the high intraspecific variation in the viperine snake. A more detailed analysis of the European samples based on genetic diversity indices and a network reconstruction suggests a complex Pleistocene history for the viperine snake in Europe. Clear differentiation was found between populations south and north of the central Iberian mountain ranges, suggesting Pleistocene glacial refugia both in the southern and northern Iberian peninsula. In the south, genetic diversity was associated with the main river valleys, whereas northern haplotypes were more broadly distributed, indicating gene flow or postglacial range expansions. Unexpectedly high levels of genetic variation in southeastern France and northwestern Italy would be compatible with the hypothesis of a glacial refugium north of the Pyrenees or in Italy. However, due to the dependence of N. maura on warm climates, the assumption of a northern refugium seems unwarranted. We believe that further sampling in northern Spain is likely to reveal genetically diverse populations which could have served as sources for postglacial recolonization of France and Italy.     

Phylogeography of the burnet moth Zygaena transalpina complex: molecular and morphometric differentiation suggests glacial refugia in Southern France,Western France and micro‐refugia within the Alps

Patterns of common recolonization routes from glacial refugia to Central Europe during the Pleistocene are generalized to paradigms of postglacial recolonization in Europe. Recent studies indicate, however, that the actual phylogeographic history of many species might be more complex and cannot be simplified to generalized patterns. Burnet moths of the Zygaena transalpina complex represent a group of closely related taxa, which are considered as a typical example for these generalized patterns. At present, three groups are recognized that are assumed to have spread from three classical refugia in Western Europe, Italy and the Balkans to Central Europe. Here, we re‐investigate their phylogeography using a combined molecular and morphometric approach. Phylogenetic and nested clade phylogeographic analyses of 476 samples from 55 localities taken from Southern and Central Europe reveal that the Zygaena transalpina complex consists of three distinct haplotype clusters, which geographically roughly correspond to possible refugia in Western Europe, Italy and the Balkans. A synthesis of the data with a geometric morphometry dataset of 425 specimens from 46 localities corroborates this molecular result but differs in several aspects. Important new aspects are multiple refugia of the western ‘hippocrepidis’ branch and micro‐habitats within the Alps of the central ‘transalpina’ branch. Further, our results display a more complex phylogeographic pattern for this species complex, which is not tractable with a rigid, generalized pattern.     

Phylogeography and historical demography of the Lusitanian snail <Emphasis Type="Italic">Elona quimperiana</Emphasis> reveal survival in unexpected separate glacial refugia

Background  

Present day distributions of Palearctic taxa in northern latitudes mainly result from populations having survived in local patches during the Late Pleistocene and/or from recolonizing populations from southern temperate refugia. If well-studied Mediterranean and eastern European refugia are widely accepted, some recent biogeographical assumptions still remain unclear, such as the occurrence of multiple glacial refugia in Iberia and cryptic refugia in northern Europe during the last glaciations. The Lusitanian snail Elona quimperiana has a remarkably disjunct distribution, limited to northwestern France (Brittany), northwestern Spain and the Basque Country. By describing the phylogeographical structure of this species across its entire range, the present study attempts to identify refugia and subsequent recolonization routes.     

Multiple Glacial Refugia of the Low-Dispersal Ground Beetle Carabus irregularis: Molecular Data Support Predictions of Species Distribution Models

Classical glacial refugia such as the southern European peninsulas were important for species survival during glacial periods and acted as sources of post-glacial colonisation processes. Only recently, some studies have provided evidence for glacial refugia north of the southern European peninsulas. In the present study, we combined species distribution models (SDMs) with phylogeographic analyses (using mitochondrial DNA = mtDNA) to investigate if the cold-adapted, stenotopic and flightless ground beetle species, Carabus irregularis, survived the Last Glacial Maximum (LGM) in classical and/or other refugia. SDMs (for both a western European and for a Carpathian subgroup) were calculated with MAXENT on the basis of 645 species records to predict current and past distribution patterns. Two mtDNA loci (CO1 and ND5, concatenated sequence length: 1785 bp) were analyzed from 91 C. irregularis specimens to reconstruct the phylogeography of Central and eastern European populations and to estimate divergence times of the given lineages. Strong intra-specific genetic differentiation (inter-clade Φ_ST values ranged from 0.92 to 0.99) implied long-term isolation of major clades and subsclades. The high divergence between the nominate subspecies and the Carpathian subspecies C. i. montandoni points to two independent species rather than subspecies (K-2P distance 0.042 ± 0.004; supposed divergence of the maternal lineages dated back 1.6 to 2.5 million years BP) differing not only morphologically but also genetically and ecologically from each other. The SDMs also inferred classical as well as other refugia for C. irregularis, especially north of the Alps, in southeastern Europe and in the Carpathians. The coincidences between the results of both methods confirm the assumption of multiple glacial refugia for the studied species and the usefulness of combining methodological approaches for the understanding of the history of low-dispersal insect species.     

Glacial in situ survival in the Western Alps and polytopic autopolyploidy in Biscutella laevigata L. (Brassicaceae)

Past climatic changes and especially the ice ages have had a great impact on both the distribution and the genetic composition of plant populations, but whether they promoted speciation is still controversial. The autopolyploid complex Biscutella laevigata is a classical example of polyploidy linked to glaciations and is an interesting model to explore migration and speciation driven by climate changes in a complex alpine landscape. Diploid taxa survived the last glacial maximum in several never-glaciated areas and autotetraploids are clearly dominant in the central parts of the Alps; however, previous range-wide studies failed to identify their diploid ancestor(s). This study highlights the phylogeographical relationships of maternal lineages in the Western Alps and investigates the polyploidy process using plastid DNA sequences (trnS-trnG and trnK-intron) combined with plastid DNA length polymorphism markers, which were transferable among Brassicaceae species. Twenty-one distinct plastid DNA haplotypes were distinguished in 67 populations densely sampled in the Western Alps and main lineages were identified by a median-joining network. The external Alps harboured high levels of genetic diversity, while the Central Alps contained only a subset of haplotypes due to postglacial recolonization. Several haplotypes were restricted to local peripheral refugia and evidence of in situ survival in central nunataks was detected by the presence of highly differentiated haplotypes swamped by frequent ones. As hierarchical genetic structure pointed to an independent evolution of the species in different biogeographical districts, and since tetraploids displayed haplotypes belonging to different lineages restricted to either the northern or the southern parts of the Alpine chain, polytopic autopolyploidy was also apparent in the Western Alps.     

Glacial refugia and postglacial expansion of the alpine–prealpine plant species Polygala chamaebuxus

The shrubby milkwort (Polygala chamaebuxus L.) is widely distributed in the Alps, but occurs also in the lower mountain ranges of Central Europe such as the Franconian Jura or the Bohemian uplands. Populations in these regions may either originate from glacial survival or from postglacial recolonization. In this study, we analyzed 30 populations of P. chamaebuxus from the whole distribution range using AFLP (Amplified Fragment Length Polymorphism) analysis to identify glacial refugia and to illuminate the origin of P. chamaebuxus in the lower mountain ranges of Central Europe. Genetic variation and the number of rare fragments within populations were highest in populations from the central part of the distribution range, especially in the Southern Alps (from the Tessin Alps and the Prealps of Lugano to the Triglav Massiv) and in the middle part of the northern Alps. These regions may have served, in accordance with previous studies, as long‐term refugia for the glacial survival of the species. The geographic pattern of genetic variation, as revealed by analysis of molecular variance, Bayesian cluster analysis and a PopGraph genetic network was, however, only weak. Instead of postglacial recolonization from only few long‐term refugia, which would have resulted in deeper genetic splits within the data set, broad waves of postglacial expansion from several short‐term isolated populations in the center to the actual periphery of the distribution range seem to be the scenario explaining the observed pattern of genetic variation most likely. The populations from the lower mountain ranges in Central Europe were more closely related to the populations from the southwestern and northern than from the nearby eastern Alps. Although glacial survival in the Bohemian uplands cannot fully be excluded, P. chamaebuxus seems to have immigrated postglacially from the southwestern or central‐northern parts of the Alps into these regions during the expansion of the pine forests in the early Holocene.     

Tracing the evolutionary history of the mole,Talpa europaea,through mitochondrial DNA phylogeography and species distribution modelling

Our understanding of the effect of Pleistocene climatic changes on the biodiversity of European mammals mostly comes from phylogeographical studies of non‐subterranean mammals, whereas the influence of glaciation cycles on subterranean mammals has received little attention. The lack of data raises the question of how and to what extent the current amount and distribution of genetic variation in subterranean mammals is the result of Pleistocene range contractions/expansions. The common mole (Talpa europaea) is a strictly subterranean mammal, widespread across Europe, and represents one of the best candidates for studying the influence of Quaternary climatic oscillation on subterranean mammals. Cytochrome b sequences, as obtained from a sampling covering the majority of the distribution area, were used to evaluate whether Pleistocene climate change influenced the evolution of T. europaea and left a trace in the genetic diversity comparable to that observed in non‐subterranean small mammals. Subsequently, we investigated the occurrence of glacial refugia by comparing the results of phylogeographical analysis with species distribution modelling. We found three differentiated mitochondrial DNA lineages: two restricted to Spain and Italy and a third that was widespread across Europe. Phylogenetic inferences and the molecular clock suggest that the Spanish moles represent a highly divergent and ancient lineage, highlighting for the first time the paraphyly of T. europaea. Furthermore, our analyses suggest that the genetic break between the Italian and the European lineages predates the last glacial phase. Historical demography and spatial principal component analysis further suggest that the Last Glacial Maximum left a signature both in the Italian and in the European lineages. Genetic data combined with species distribution models support the presence of at least three putative glacial refugia in southern Europe (France, Balkan Peninsula and Black Sea) during thelast glacial maximum that likely contributed to post‐glacial recolonization of Europe. By contrast, the Italian lineage remained trapped in the Italian peninsula and, according to the pattern observed in other subterranean mammals, did not contribute to the recolonization of northern latitudes. © 2015 The Linnean Society of London, Biological Journal of the Linnean Society, 2015, 114 , 495–512.     

Glacial refugia,recolonization patterns and diversification forces in Alpine‐endemic Megabunus harvestmen

The Pleistocene climatic fluctuations had a huge impact on all life forms, and various hypotheses regarding the survival of organisms during glacial periods have been postulated. In the European Alps, evidence has been found in support of refugia outside the ice shield (massifs de refuge) acting as sources for postglacial recolonization of inner‐Alpine areas. In contrast, evidence for survival on nunataks, ice‐free areas above the glacier, remains scarce. Here, we combine multivariate genetic analyses with ecological niche models (ENMs) through multiple timescales to elucidate the history of Alpine Megabunus harvestmen throughout the ice ages, a genus that comprises eight high‐altitude endemics. ENMs suggest two types of refugia throughout the last glacial maximum, inner‐Alpine survival on nunataks for four species and peripheral refugia for further four species. In some geographic regions, the patterns of genetic variation are consistent with long‐distance dispersal out of massifs de refuge, repeatedly coupled with geographic parthenogenesis. In other regions, long‐term persistence in nunataks may dominate the patterns of genetic divergence. Overall, our results suggest that glacial cycles contributed to allopatric diversification in Alpine Megabunus, both within and at the margins of the ice shield. These findings exemplify the power of ENM projections coupled with genetic analyses to identify hypotheses about the position and the number of glacial refugia and thus to evaluate the role of Pleistocene glaciations in driving species‐specific responses of recolonization or persistence that may have contributed to observed patterns of biodiversity.     

Phylogeography and Pleistocene refugia of the Little Owl Athene noctua inferred from mtDNA sequence data

Pleistocene glaciations greatly affected the distribution of genetic diversity in animal populations. The Little Owl is widely distributed in temperate regions and could have survived the last glaciations in southern refugia. To describe the phylogeographical structure of European populations, we sequenced the mitochondrial cytochrome c oxidase I (COI) and control region (CR1) in 326 individuals sampled from 22 locations. Phylogenetic analyses of COI identified two deeply divergent clades: a western haplogroup distributed in western and northwestern Europe, and an eastern haplogroup distributed in southeastern Europe. Faster evolving CR1 sequences supported the divergence between these two main clades, and identified three subgroups within the eastern clade: Balkan, southern Italian and Sardinian. Divergence times estimated from COI with fossil calibrations indicate that the western and eastern haplogroups split 2.01–1.71 Mya. Slightly different times for splits were found using the standard 2% rate and 7.3% mtDNA neutral substitution rate. CR1 sequences dated the origin of endemic Sardinian haplotypes at 1.04–0.26 Mya and the split between southern Italian and Balkan haplogroups at 0.72–0.21 Mya, coincident with the onset of two Pleistocene glaciations. Admixture of mtDNA haplotypes was detected in northern Italy and in central Europe. These findings support a model of southern Mediterranean and Balkan refugia, with postglacial expansion and secondary contacts for Little Owl populations. Central and northern Europe was predominantly recolonized by Little Owls from Iberia, whereas expansion out of the Balkans was more limited. Northward expansion of the Italian haplogroup was probably prevented by the Alps, and the Sardinian haplotypes remained confined to the island. Results showed a clear genetic pattern differentiating putative subspecies. Genetic distances between haplogroups were comparable with those recorded between different avian species.     

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.