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.     

Phylogeographic analysis reveals a deep lineage split within North Atlantic Littorina saxatilis

Phylogeographic studies provide critical insight into the evolutionary histories of model organisms; yet, to date, range-wide data are lacking for the rough periwinkle Littorina saxatilis, a classic example of marine sympatric speciation. Here, we use mitochondrial DNA (mtDNA) sequence data to demonstrate that L. saxatilis is not monophyletic for this marker, but is composed of two distinct mtDNA lineages (I and II) that are shared with sister species Littorina arcana and Littorina compressa. Bayesian coalescent dating and phylogeographic patterns indicate that both L. saxatilis lineages originated in the eastern North Atlantic, around the British Isles, at approximately 0.64 Ma. Both lineages are now distributed broadly across the eastern, central and western North Atlantic, and show strong phylogeographic structure among regions. The Iberian Peninsula is genetically distinct, suggesting prolonged isolation from northeastern North Atlantic populations. Western North Atlantic populations of L. saxatilis lineages I and II predate the last glacial maximum and have been isolated from eastern North Atlantic populations since that time. This identification of two distinct, broadly distributed mtDNA lineages further complicates observed patterns of repeated incipient ecological speciation in L. saxatilis, because the sympatric origins of distinct ecotype pairs on eastern North Atlantic shores may be confounded by admixture of divergent lineages.     

Northern and Southern expansions of Atlantic brown trout (Salmo trutta) populations during the Pleistocene

The phylogeography of Atlantic brown trout ( Salmo trutta ) was analysed using mitochondrial DNA control region complete sequences of 774 individuals from 57 locations. Additionally, the available haplotype information from 100 published populations was incorporated in the analysis. Combined information from nested clade analysis, haplotype trees, mismatch distributions, and coalescent simulations was used to characterize population groups in the Atlantic basin. A major clade involved haplotypes assigned to the Atlantic (AT) lineage, but another major clade should be considered as a distinct endemic lineage restricted to the Iberian Peninsula. The phylogeography of the Atlantic populations showed the mixed distribution of several Atlantic clades in glaciated areas of Northern Europe, whereas diverged haplotypes dominated the coastal Iberian rivers. Populations inhabiting the Atlantic rivers of southern France apparently contributed to postglacial colonization of northern basins, but also comprised the source of southern expansions during the Pleistocene.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 97 , 904–917.     

A new cytochrome b phylogroup of the common vole (Microtus arvalis) endemic to the Balkans and its implications for the evolutionary history of the species

The phylogeographic architecture of the common vole, Microtus arvalis, has been well‐studied using mitochondrial DNA and used to test hypotheses relating to glacial refugia. The distribution of the five described cytochrome b (cyt b) lineages in Europe west of Russia has been interpreted as a consequence of postglacial expansion from both southern and central European refugia. A recently proposed competing model suggests that the ‘cradle’ of the M. arvalis lineages is in western central Europe from where they dispersed in different directions after the Last Glacial Maximum. In the present study, we report a new cyt b lineage of the common vole from the Balkans that is not closely related to any other lineage and whose presence might help resolve these issues of glacial refugia. The Balkan phylogroup occurs along the southern distributional border of M. arvalis in central and eastern Bosnia and Herzegovina, Montenegro, and eastern Serbia. Further north and west in Slovenia, Bosnia and Herzegovina, and Serbia, common voles belong to the previously‐described Eastern lineage, whereas both lineages are sympatric in one site in Bosnia and Herzegovina. The Balkan phylogroup most reasonably occupied a glacial refugium already known for various Balkan endemic species, in contrast to the recently proposed model. South‐east Europe is an absolutely crucial area for understanding the postglacial colonization history of small mammals in Europe and the present study adds to the very few previous detailed phylogeographic studies of this region. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 788–796.     

Pleistocene divergence of Dinaric Drusus endemics (Trichoptera, Limnephilidae) in multiple microrefugia within the Balkan Peninsula

The Balkan Peninsula is one of three major European refugial areas. It has high biodiversity and endemism, but data on the age and origin of its fauna, especially endemics, are limited. Mitochondrial sequence data (COI and 16S genes) were used to study the population structure and phylogeography of the caddisfly Drusus croaticus and the phylogeny and divergence of seven other Drusus species, mostly range-restricted endemics of the Dinaric region of the Balkan Peninsula. The divergence of D. croaticus populations in Croatia and allopatric Drusus species in Bosnia dated to the Pleistocene, showing the importance of this time period for the origin and diversification of Balkan endemic taxa. The divergence of more distantly related species dated to the Late Miocene/Early Pliocene. Population genetic and phylogeographic analysis of 115 individuals from 11 populations of D. croaticus revealed a high level of genetic differentiation and absence of gene flow between populations separated by more than 10 km. The existence of allopatrically fragmented lineages in D. croaticus and the endemic Bosnian species is most likely the result of long-term isolation in multiple microrefugia, probably due to the specific habitat requirements and life-history traits of Drusinae coupled with the topographic complexity and historical changes in geomorphology of the region. Overall, these findings shed light on the processes generating the high genetic complexity of this refugial region that parallels the 'refugia within refugia' pattern widely reported from the Iberian refugium.     

Phylogeographic structure and past history of the circum-Mediterranean species Tomicus destruens Woll. (Coleoptera: Scolytinae)

Phylogeographic studies are often focused on temperate European species with relict footholds in the Mediterranean region. Past climatic oscillations usually induced range contractions and expansions from refugial areas located in southern Europe, and spatial distribution of genetic diversity show that northward expansions were usually pioneer-like. Actually, few studies have focused on circum-Mediterranean species, which probably were not influenced in the same way by climatic oscillations. We present the phylogeography of the bark beetle Tomicus destruens, which is restricted to the whole Mediterranean basin and the Atlantic coasts of North Africa and Portugal. We systematically sequenced 617 bp of the mitochondrial genes COI and COII for 42 populations (N = 219). Analysis revealed 53 haplotypes geographically structured in two clades, namely eastern and western clades, that diverged during the Pleistocene. A contact zone was identified along the Adriatic coast of Italy. Interestingly, we found contrasting levels of genetic structure within each clade. The eastern group was characterized by a significant phylogeographic pattern and low levels of gene flow, whereas the western group barely showed a spatial structure in haplotype distribution. Moreover, the main pine hosts were different between groups, with the Aleppo-brutia complex in the east and the maritime pine in the west. Potential roles of host species, climatic parameters and geographical barriers are discussed and the phylogeographic patterns are compared to classical models of postglacial recolonization in Europe.     

Variation in mitochondrial DNA reveals multiple distant glacial refugia in black spruce (Picea mariana), a transcontinental North American conifer

Range-wide genetic variation of black spruce (Picea mariana) was studied using polymerase chain reaction-random fragment length polymorphism markers of the mitochondrial genome. Four polymorphic mitochondrial DNA (mtDNA) loci were surveyed and two or three alleles were detected at each locus, resulting in 10 multilocus mtDNA types or mitotypes. A significant subdivision of population genetic diversity was detected (GST = 0.671; NST = 0.726), suggesting low levels of gene flow among populations. The distribution of mitotypes was not random (NST > GST; P < 0.05) and revealed four partially overlapping zones, presumably representative of different glacial populations. Comparison of the genetic structure derived from mtDNA markers and the colonization paths previously deduced from the fossil and pollen records allow us to infer at least three southern and one northeastern glacial populations for black spruce. The patterns revealed in this study suggest that black spruce shares its biogeographical history with other forest-associated North American species.     

Refugia within refugia – patterns in endemism and genetic divergence are linked to Late Quaternary climate stability in the Iberian Peninsula

In Europe, southern peninsulas served as major refugia during Pleistocene cold periods. However, growing evidence has revealed complex patterns of glacial survival within these southern regions, with multiple glacial refugia within each larger refugial area. We investigated the extent to which patterns of endemism and phylogeographic are concordant across animal species in the Iberian Peninsula, one of the most important unglaciated areas in Europe during the Pleistocene, can be explained in terms of climatic stability. We found that historical climatic stability (notably climate velocity measures integrating macroclimatic shifts with local spatial topoclimate gradients) was often among the most important predictors of endemic species richness for different taxonomic groups using models that also incorporated measures of modern climate. Furthermore, for some taxonomic groups, climatic stability was also correlated with patterns of spatial concordance in interpopulation genetic divergence across multiple taxa, and private haplotypes were more frequently found in relatively stable areas. Overall, our results suggest that both endemism patterns and cross‐taxa concordant phylogeographic patterns across the Iberian Peninsula to some extent are linked to spatial variation in Late Quaternary climate stability, in agreement with the proposed ‘refugia‐within‐refugia’ scenario. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 13–28.     

Phylogeography of red deer (Cervus elaphus) in Europe

Aim To investigate the phylogeographical patterns of red deer (Cervus elaphus) in Europe, and to disentangle the influence of ancient (e.g. Pleistocene ice ages) from more recent processes (e.g. human translocations). Location Europe. Methods In this study we provide by far the most extensive analysis of genetic structure in European red deer, based on analyses of variation at two mitochondrial markers (cyt b and D‐loop) in a large number of individuals from 39 locations. Relationships of mitochondrial DNA haplotypes were determined using minimum spanning networks and phylogenetic analyses. Population structure was examined by analyses of molecular variance. Historical processes shaping the present patterns were inferred from nested clade analysis and nucleotide diversity statistics. Results Within Europe, we detected three deeply divergent mitochondrial DNA lineages. The three lineages displayed a phylogeographical pattern dividing individuals into western European, eastern European and Mediterranean (Sardinia, Spain and Africa) groups, suggesting contraction into three separate refugia during the last glaciation. Few haplotypes were shared among these three groups, a finding also confirmed by F_ST values. Calculations of divergence times suggest that the groups probably split during the Pleistocene. Main conclusions The observed pattern is interpreted to result from isolation in different refugia during the last glaciation. The western and eastern European lineages could be linked to an Iberian and Balkan refugium, respectively. The third lineage might originate from a Sardinian or African refugium. We link local phylogeographical patterns observed in Europe to the post‐glacial recolonization process, shaped by the geographical localization of refugia and barriers to gene flow. Regardless of the importance of red deer as a game species and the tradition of translocating red deer in Europe, we detected few individuals that did not match the trichotomous pattern, suggesting that translocations have occurred mainly at smaller spatial scales.     

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.     

Phylogeography of endemic ermine (Mustela erminea) in southeast Alaska

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

Contrasting phylogeographic patterns in Sphagnum fimbriatum and Sphagnum squarrosum (Bryophyta, Sphagnopsida) in Europe

The chloroplast phylogeography of two peat mosses (Sphagnum fimbriatum and Sphagnum squarrosum) with similar distributions but different life history characteristics was investigated in Europe. Our main aim was to test whether similar distributions reflect similar phylogeographic and phylodemographic processes. Accessions covering the European distributions of the species were collected and approx. 2000 bp of the chloroplast genome of each species was sequenced. Maximum parsimony, statistical parsimony and phylodemographic analyses were used to address the question of whether these species with similar distributions show evidence of similar phylogeographic and phylodemographic processes. The chloroplast haplotypes of the currently spreading species S. fimbriatum showed strong geographic structure, whereas those of S. squarrosum, which has stable historical population sizes, showed only very weak geographic affinity and were widely distributed. We hypothesize that S. fimbriatum survived the last glaciations along the Atlantic coast of Europe, whereas S. squarrosum had numerous, scattered refugia in Europe. The dominance of one haplotype of S. fimbriatum across almost all of Europe suggests rapid colonization after the last glacial maximum. We hypothesize that high colonizing ability is an inherent characteristic of the species and its recent expansion in Europe is a response to climate change.     

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.     

A mitochondrial DNA minisatellite reveals the postglacial history of jack pine (Pinus banksiana), a broad-range North American conifer

Jack pine (Pinus banksiana Lamb.) is a broadly distributed North American conifer and its current range was covered by the Laurentian ice sheet during the last glacial maximum. To infer about the history and postglacial colonization of this boreal species, range-wide genetic variation was assessed using a new and highly variable minisatellite-like marker of the mitochondrial genome. Among the 543 trees analysed, 14 distinct haplotypes were detected, which corresponded to different repeat numbers of the 32-nucleotide minisatellite-like motif. Several haplotypes were rare with limited distribution, suggesting recent mutation events during the Holocene. At the population level, an average of 2.6 haplotypes and a mean haplotype diversity (H) of 0.328 were estimated. Population subdivision of genetic diversity was quite high with G(ST) and R(ST) values of 0.569 and 0.472, respectively. Spatial analyses identified three relatively homogeneous groups of populations presumably representative of genetically distinct glacial populations, one west and one east of the Appalachian Mountains in the United States and a third one presumably on the unglaciated northeastern coastal area in Canada. These results indicate the significant role of the northern part of the US Appalachian Mountains as a factor of vicariance during the ice age. A fourth distinct group of populations was observed in central Québec where the continental glacier retreated last. It included populations harbouring haplotypes present into the three previous groups, and it had higher level of haplotype diversity per population (H = 0.548) and lower population differentiation (G(ST) = 0.265), which indicates a zone of suture or secondary contact between the migration fronts of the three glacial populations. Introgression from Pinus contorta Dougl. var. latifolia Engelm. was apparent in one western population from Alberta. Altogether, these results indicate that the mitochondrial DNA variation of jack pine is geographically highly structured and it correlates well with large-scale patterns emerging from recent phylogeographical studies of other tree boreal species in North America.     

Geographical variation in genetic structure of an Atlantic Coastal Forest frog reveals regional differences in habitat stability

Climatic oscillations throughout the Pleistocene combined with geological and topographic complexity resulted in extreme habitat heterogeneity along the Atlantic coast of Brazil. Inferring how these historic landscape patterns have structured the current diversity of the region's biota is important both for our understanding of the factors promoting diversification, as well as the conservation of this biodiversity hotspot. Here we evaluate potential historical scenarios of diversification in the Atlantic Coastal Forest of Brazil by investigating the population genetic structure of a frog endemic to the region. Using mitochondrial and nuclear sequences, we generated a Bayesian population-level phylogeny of the Thoropa miliaris species complex. We found deep genetic divergences among five geographically distinct clades. Southern clades were monophyletic and nested within paraphyletic northern clades. Analyses of historical demographic patterns suggest an overall north to south population expansion, likely associated with regional differences in habitat stability during the Pliocene and early Pleistocene. However, genetic structure among southern populations is less pronounced and likely represents more recent vicariant events resulting from Holocenic sea-level oscillations. Our analyses corroborate that the Atlantic Coastal Forest has been a biogeographically dynamic landscape and suggest that the high diversity of its fauna and flora resulted from a combination of climatic and geologic events from the Pliocene to the present.     

Phylogeography of the greater horseshoe bat, Rhinolophus ferrumequinum: contrasting results from mitochondrial and microsatellite data

Phylogeographical studies are typically based on haplotype data, occasionally on nuclear markers such as microsatellites, but rarely combine both. This is unfortunate because the use of markers with contrasting modes of inheritance and rates of evolution might provide a more accurate and comprehensive understanding of a species' history. Here we present a detailed study of the phylogeography of the greater horseshoe bat, Rhinolophus ferrumequinum , using 1098 bp of the mitochondrial ND2 gene from 45 localities from across its Palaearctic range to infer population history. In addition, we re-analysed a large microsatellite data set available for this species and compared the results of both markers to infer population relationships and the historical processes influencing them. We show that mtDNA, the most popular marker in phylogeography studies, yielded a misleading result, and would have led us to conclude erroneously that a single expansion had taken place in Europe. Only by combining the mitochondrial and microsatellite data sets are we able to reconstruct the species' history and show two colonization events in Europe, one before the Last Glacial Maximum (LGM) and one after it. Combining markers also revealed the importance of Asia Minor as an ancient refugium for this species and a source population for the expansion of the greater horseshoe bat into Europe before the LGM.     

Phylogeography of the vairone (Leuciscus souffia,Risso 1826) in Central Europe

The vairone Leuciscus souffia is a cyprinid fish that inhabits river systems in and around the Alps. The complete mitochondrial DNA control region (945 bp) was sequenced in 295 vairone from 22 populations in Central Europe. A total of 51 haplotypes were identified with a maximum pairwise distance between haplotypes of 5.6%. Phylogenetic analyses revealed two major clades in L. souffia, an ‘Italian’ clade, and an ‘Alpine’ clade. Two hybrid zones exist, in the Mediterranean Alps and in the Soca basin. The position of the sister species of L. souffia, L. turskyi, to the ‘Alpine’ and the ‘Italian’ clade could not be resolved unambiguously. However, a linearized tree analysis indicated that L. turskyi represents a third lineage, that originated at the same time as the ‘Alpine’ and the ‘Italian’ clades of L. souffia. In the ‘Alpine’ clade two groups were resolved, a subclade with haplotypes from the Rhône and Var basins and a cluster with haplotypes from the Danube and Rhine systems. Our data suggest a long history of the vairone in Central Europe, predating Pleistocene glacial cycles. Two main refugia during glaciations must have existed, one in Italy and another one most probably in the Danube system. However, age estimates based on molecular clock calibrations suggest the survival of ‘Alpine’ haplotypes in several drainages during the last glaciation cycles. The Rhine system was only recently colonized.