Mitochondrial DNA phylogeography of European hedgehogs

European hedgehog populations belonging to Erinaceus europaeus and E . concolor have been investigated by mitochondrial DNA analysis. A 383 bp fragment of the cytochrome b gene has been sequenced and maximum parsimony and neighbour-joining trees of Tamura–Nei genetic distance values have been constructed. Similar topologies have been produced by both methods, showing a deep divergence between E . europaeus and E . concolor and a further subdivision of each species into a western and an eastern clade. A comparison with previously published allozyme data is made, and concordant and discordant patterns are discussed. The influence of Pleistocene glaciations on the observed pattern of divergence is inferred.     

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.     

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.     

Comparative phylogeography and demographic history of the wood lemming (Myopus schisticolor): implications for late Quaternary history of the taiga species in Eurasia

The association between demographic history, genealogy and geographical distribution of mitochondrial DNA cytochrome b haplotypes was studied in the wood lemming (Myopus schisticolor), a species that is closely associated with the boreal forest of the Eurasian taiga zone from Scandinavia to the Pacific coast. Except for a major phylogeographic discontinuity (0.9% nucleotide divergence) in southeastern Siberia, only shallow regional genetic structure was detected across northern Eurasia. Genetic signs of demographic expansions imply that successive range contractions and expansions on different spatial scales represented the primary historical events that shaped geographical patterns of genetic variation. Comparison of phylogeographic structure across a taxonomically diverse array of other species that are ecologically associated with the taiga forest revealed similar patterns and identified two general aspects. First, the major south-north phylogeographic discontinuity observed in five out of six species studied in southeastern Siberia and the Far East implies vicariant separation in two different refugial areas. The limited distribution range of the southeastern lineages provides no evidence of the importance of the putative southeastern refugial area for postglacial colonization of northern Eurasia by boreal forest species. Second, the lack of phylogeographic structure associated with significant reciprocal monophyly and genetic signatures of demographic expansion in all nine boreal forest animal species studied to date across most of northern Eurasia imply contraction of each species to a single refugial area during the late Pleistocene followed by range expansion on a continental scale. Similar phylogeographic patterns observed in this taxonomically diverse set of organisms with different life histories and dispersal potentials reflect the historical dynamics of their shared environment, the taiga forest in northern Eurasia.     

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.     

Holarctic phylogeography of the root vole (Microtus oeconomus): implications for late Quaternary biogeography of high latitudes

A species-wide phylogeographical study of the root vole (Microtus oeconomus) was performed using the whole 1140 base pair mitochondrial (mt) cytochrome b gene. We examined 83 specimens from 52 localities resulting in 65 unique haplotypes. Our results demonstrate that the root vole is divided into four main mtDNA phylogenetic lineages that seem to have largely allopatric distributions. Net divergence estimates (2.0-3.5%) between phylogroups, as well as relatively high nucleotide diversity estimates within phylogroups, indicate that the distinct phylogeographical structure was initiated by historical events that predated the latest glaciation. European root voles are divided into a Northern and a Central mtDNA phylogroup. The mtDNA data in concert with fossil records imply that root voles remained north of the classical refugial areas in southern Europe during the last glacial period. The currently fragmented populations in central Europe belong to a single mtDNA phylogroup. The Central Asian and the North European lineages are separated by the Ural Mountains, a phylogeographical split also found in collared lemmings (Dicrostonyx) and the common vole (M. arvalis). The Beringian lineage occurs from eastern Russia through Alaska to northwestern Canada. This distribution is congruent with the traditional boundaries of the Beringian refugium and with phylogeographical work on other organisms. In conclusion, similarities between the phylogeographical patterns in the root vole and other rodents, such as Arctic and subarctic lemmings, as well as more temperate vole species, indicate that late Quaternary geological and climatic events played a strong role in structuring northern biotic communities.     

Endemic species may have complex histories: within‐refugium phylogeography of an endangered Iberian vole

Glacial refugia protected and promoted biodiversity during the Pleistocene, not only at a broader scale, but also for many endemics that contracted and expanded their ranges within refugial areas. Understanding the evolutionary history of refugial endemics is especially important in the case of endangered species to recognize the origins of their genetic structure and thus produce better informed conservation practices. The Iberian Peninsula is an important European glacial refugium, rich in endemics of conservation concern, including small mammals, such as the Cabrera vole (Microtus cabrerae). This near‐threatened rodent is characterized by an unusual suite of genetic, life history and ecological traits, being restricted to isolated geographic nuclei in fast‐disappearing Mediterranean subhumid herbaceous habitats. To reconstruct the evolutionary history of the Cabrera vole, we studied sequence variation at mitochondrial, autosomal and sex‐linked loci, using invasive and noninvasive samples. Despite low overall mitochondrial and nuclear nucleotide diversities, we observed two main well‐supported mitochondrial lineages, west and east. Phylogeographic modelling in the context of the Cabrera vole's detailed fossil record supports a demographic scenario of isolation of two populations during the Last Glacial Maximum from a single focus in the southern part of the Iberian Peninsula. In addition, our data suggest subsequent divergence within the east, and secondary contact and introgression of the expanding western population, during the late Holocene. This work emphasizes that refugial endemics may have a phylogeographic history as rich as that of more widespread species, and conservation of such endemics includes the preservation of that genetic legacy.     

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

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

The limits of mtDNA phylogeography: complex patterns of population history in a highly structured Iberian lizard are only revealed by the use of nuclear markers

Phylogeographic analyses based on the sole use of the mitochondrial DNA (mtDNA) molecule reveal only a small part of the evolutionary history of a species or a set of related species. In this study, we have combined the application of slow- and fast-evolving nuclear markers (proteins and microsatellites, respectively) together with the analysis of two-gene genealogies to further understand the history of the Iberian endemic Schreiber's green lizard, Lacerta schreiberi, a species for which a well established phylogeographical scenario is available. In sharp contrast with the observation of four divergent and almost allopatric mtDNA clades, our nuclear data revealed how two groups of populations diverged, persisted and began to admix along the mountains of the Iberian Central System. In addition, the combination of mtDNA and nuclear data showed how the core area of the species distribution responded to ice ages, first by relatively old processes of population expansion to the south followed by episodes of contraction that are at the origin of present-day isolates, and more recently by a postglacial expansion to the Iberian Northwest where new habitats were made available after climatic amelioration. Taken together with recently published results for a variety of other organisms, our results suggest that complex processes of fragmentation, expansion and admixture can only be properly addressed through the use of several and complementary types of molecular markers. Finally, we also suggest that southern European refugia are both hotspots and melting pots of genetic diversity.     

Rangewide phylogeography of a terrestrial slug in Europe: evidence for Alpine refugia and rapid colonization after the Pleistocene glaciations

Intraspecific phylogeographical patterns largely depend on the life history traits of a species. Especially species with a high degree of cold tolerance, limited requirements towards habitat preferences, and relatively low active dispersal capacities may have responded in a different way to the Pleistocene climatological fluctuations than the majority of taxa studied so far. To evaluate this possibility, we studied Arion fuscus (Muller, 1774), a common and widespread European terrestrial slug, from 88 locations (N = 964). Sequence variation was assessed for fragments of the mitochondrial 16S rDNA and COI genes by means of single-strand conformation polymorphisms (SSCP) and subsequent DNA sequencing. Additionally, eight allozyme loci were scored in 843 individuals. Phylogenetic analysis revealed the presence of two major evolutionary lineages, one in the Balkan region and another in the Alps and the rest of Europe. The sequence divergence between the two lineages was limited (3.3%), but gene flow between the regions was absent, suggesting that the two regions have been isolated since the late Pliocene or early Pleistocene. Allozyme differentiation among geographical regions and mitochondrial DNA (mtDNA) lineages was low. The geographical patterns observed in our data showed that (i) haplotype and nucleotide diversities are very low in northern Europe, suggesting that single haplotypes rapidly colonized large areas; (ii) recently expanded haplotype clades have restricted distribution ranges, suggesting that current gene flow is low; and (iii) genetic diversity in the Alps is much higher than in other regions and estimated past gene flow from the Eastern Alps to other regions was high, suggesting that this was a refugial zone during the Pleistocene. This full-range phylogeography suggests the existence of an alternative refugial zone, situated north of the refugial areas currently recognized in most other taxa.     

Molecular phylogeography and the evolution and conservation of Amazonian mammals

The phylogeographic structure of 15 genera of Amazonian marsupials and rodents is summarized based on comparative sequence of the mitochondrial cytochrome b gene. The data are limited in geographical coverage, with samples widely scattered throughout Amazonia from the base of the Andes in Peru to the Guianan coast and eastern Brazil. We use this approach to define species boundaries, based minimally on the principle of reciprocal monophyly, in conjunction with morphological or other genetic discontinuities. The taxa so defined are older than previously appreciated, with many lineages dating from 1 to more than 3 Myr, and thus apparently predating the early Pleistocene. We relate patterns of concordant geographical shifts with underlying tectonic history and to current positions of major rivers. Finally, we provide comments on the utility of these data and patterns to conservation, articulating a need to incorporate phylogeographic information as part of the rationale in establishing conservation priorities at the organismal and geographical area levels.     

Evolutionary history of the European whitefish Coregonus lavaretus (L.) species complex as inferred from mtDNA phylogeography and gill-raker numbers

We compared mitochondrial DNA and gill-raker number variation in populations of the European whitefish Coregonus lavaretus (L.) species complex to illuminate their evolutionary history, and discuss mechanisms behind diversification. Using single-strand conformation polymorphism (SSCP) and sequencing 528 bp of combined parts of the cytochrome oxidase b (cyt b) and NADH dehydrogenase subunit 3 (ND3) mithochondrial DNA (mtDNA) regions, we documented phylogeographic relationships among populations and phylogeny of mtDNA haplotypes. Demographic events behind geographical distribution of haplotypes were inferred using nested clade analysis (NCA) and mismatch distribution. Concordance between operational taxonomical groups, based on gill-raker numbers, and mtDNA patterns was tested. Three major mtDNA clades were resolved in Europe: a North European clade from northwest Russia to Denmark, a Siberian clade from the Arctic Sea to southwest Norway, and a South European clade from Denmark to the European Alps, reflecting occupation in different glacial refugia. Demographic events inferred from NCA were isolation by distance, range expansion, and fragmentation. Mismatch analysis suggested that clades which colonized Fennoscandia and the Alps expanded in population size 24 500-5800 years before present, with minute female effective population sizes, implying small founder populations during colonization. Gill-raker counts did not commensurate with hierarchical mtDNA clades, and poorly with haplotypes, suggesting recent origin of gill-raker variation. Whitefish designations based on gill-raker numbers were not associated with ancient clades. Lack of congruence in morphology and evolutionary lineages implies that the taxonomy of this species complex should be reconsidered.     

Comparative phylogeography of five avian species: implications for Pleistocene evolutionary history in the Qinghai‐Tibetan plateau

Pleistocene climate fluctuations have shaped the patterns of genetic diversity observed in extant species. In contrast to Europe and North America where the effects of recent glacial cycles on genetic diversity have been well studied, the genetic legacy of the Pleistocene for the Qinghai‐Tibetan (Tibetan) plateau, a region where glaciation was not synchronous with the North Hemisphere ice sheet maxima, remains poorly understood. Here, we compared the phylogeographical patterns of five avian species on the Qinghai‐Tibetan plateau by three mitochondrial DNA fragments: the Tibetan snow finch (Montifringilla adamsi), the Blanford’s snow finch (Pyrgilauda blanfordi), the horned lark (Eremophila alpestris), the twite (Carduelis flavirostris) and the black redstart (Phoenicurus ochruros). Our results revealed the three species mostly distributed on the platform region of the plateau that experienced population expansion following the retreat of the extensive glaciation period (0.5–0.175 Ma). These results are at odds with the results from avian species of Europe and North America, where population expansions occurred after Last Glacial Maximum (LGM, 0.023–0.018 Ma). A single refugium was identified in a restricted semi‐continuous area around the eastern margin of the plateau, instead of multiple independent refugia for European and North American species. For the other two species distributed on the edges of the plateau (the twite and black redstart), populations were maintained at stable levels. Edge areas are located on the eastern margin, which might have had little or no ice cover during the glaciation period. Thus, milder climate may have mitigated demographic stresses for edge species relative to the extremes experienced by platform counterparts, the present‐day ranges of which were heavily ice covered during the glaciation period. Finally, various behavioural and ecological characteristics, including dispersal capacities, habitat preference and altitude specificity along with evolutionary history might have helped to shape different phylogeographical structures appearing in these five species.     

Phylogeography of the Russian flying squirrel (Pteromys volans): implication of refugia theory in arboreal small mammal of Eurasia

A phylogeographical study of the Russian (Siberian) flying squirrel (Pteromys volans) was carried out using the complete mitochondrial (mt) cytochrome b gene sequences with special reference to the refugia theory for the arboreal traits of this species. We examined 31 specimens from 24 localities, resulting in 28 haplotypes. One breeding specimen with a unique haplotype from Hokkaido, Japan was included in the phylogenetic analysis. There were three mtDNA lineages: Hokkaido, Far Eastern, and northern Eurasia. Divergence data among lineages demonstrated that the Hokkaido group separated from the other groups during the Holsteinian interglacial. The phylogeographical pattern of P. volans is different from that previously reported for terrestrial rodents associated with treeless habitats. Unlike grasslands, forests decreased during glaciation and moved southward because of the cold and arid environmental conditions. The glacial refugia of P. volans would have been associated with forest dynamics in the Pleistocene.     

Postglacial expansion of the southern red-backed vole (Clethrionomys gapperi) in North America

Dynamic climatic oscillations of the Pleistocene dramatically changed the distributions of high latitude species. Molecular investigations of a variety of organisms show that processes of postglacial colonization of boreal regions were more complex than initially thought. Phylogeographical and coalescent analyses were conducted on partial sequences of the cytochrome b gene (600 bp) from 64 individuals of Clethrionomys gapperi from North Carolina, Pennsylvania, Minnesota, Idaho, Washington, British Columbia, Northwest Territories, and Alaska to test hypotheses relating to Pleistocene refugia and postglacial colonization routes. Three divergent clades (east, west, central) were identified with highest net divergence (dA = 5.2%) between the eastern and western clades. Populations from the recently deglaciated higher latitudes of Canada and Alaska are closely related to lower latitude populations of the central clade (dA = 1.2%) suggesting recent expansion from this midwestern region. No representatives from the east or west clade were found at latitudes higher than 50 degrees N, indicating that postglacial colonization occurred through a midcontinental route. The high latitude population from the Northwest Territories exhibited demographic patterns and genetic diversity consistent with a stable noncolonizing population. This population is found near the Mackenzie range, where the two continental ice sheets were believed to have coalesced. Molecular variation observed in this population may be the result of leading edge population diversifying in the continental corridor or may reflect the signal of a high latitude refugial population.     

Comparative phylogeography of coastal limpets across a marine disjunction in New Zealand

Cook Strait, which separates the North and South Island of New Zealand, has been a transient, but re-occurring feature of the New Zealand land mass throughout the Pleistocene, maintaining its current width and depth for the past 5000 years. Historic land fragmentation coupled with the complex hydrography of the Greater Cook Strait region has created both biogeographic and phylogeographic disjunctions between the North and South Island in several marine species. Here we use mitochondrial cytochrome b DNA sequences of three endemic intertidal limpets, Cellana ornata, Cellana radians and Cellana flava to assess intraspecific phylogeographic patterns across Cook Strait and to look for interspecific concordance of ecological and evolutionary processes among closely related taxa. We sequenced 328-359 bp in 85-321 individuals from 8-31 populations spanning the biogeographic range of the three species. Intraspecific phylogeographic analyses show moderate to strong genetic discontinuity among North and South Island populations due to allopatric fragmentation. This pattern was broadly concordant across the three species and the observed divergence among this group of intertidal limpets (0.3-2.0%) is similar to that of previously studied subtidal organisms. For each species, divergence time calculations suggest contemporary North and South Island lineages diverged from their respective most recent common ancestor approximately 200 000 to 300 000 years before present (bp), significantly earlier than previous estimates in other coastal marine taxa that arose from a miscalculation of divergence time.     

Molecular phylogeography of the Amazonian Chagas disease vectors Rhodnius prolixus and R. robustus

The phylogeographical structure of the closely related species Rhodnius prolixus and R. robustus is presented based on a 663-base pair (bp) fragment of the mitochondrial cytochrome b gene. Twenty haplotypes were recovered from 84 samples examined, representing 26 populations from seven Latin American countries. The resulting phylogenetic tree is composed of five major reciprocally monophyletic clades, one representing R. prolixus and four representing R. robustus. While R. prolixus is a very homogeneous assemblage, R. robustus has deeper clades and is paraphyletic, with the clade comprising R. robustus from Venezuela (Orinoco region) more closely related to the R. prolixus clade than to the other R. robustus populations from the Amazon region. The R. robustus paraphyly was supported further by the analysis of a nuclear gene (D2 region of the 28S RNA) for a subset of specimens. The data support the view that R. robustus represents a species complex. Levels of sequence divergence between clades within each region are compatible with a Pleistocene origin. Nucleotide diversity (pi) for all R. prolixus populations was extremely low (0.0008), suggesting that this species went through a recent bottleneck, and was subsequently dispersed by man.     

Phylogeography of Pinus tabulaeformis Carr. (Pinaceae), a dominant species of coniferous forest in northern China

How coniferous trees in northern China changed their distribution ranges in response to Quaternary climatic oscillations remains largely unknown. Here we report a study of the phylogeography of Pinus tabulaeformis, an endemic and dominant species of coniferous forest in northern China. We examined sequence variation of maternally inherited, seed‐dispersed mitochondrial DNA (mtDNA) (nad5 intron 1 and nad4/3–4) and paternally inherited, pollen‐ and seed‐dispersed chloroplast DNA (cpDNA) (rpl16 and trnS‐trnG) within and among 30 natural populations across the entire range of the species. Six mitotypes and five chlorotypes were recovered among 291 trees surveyed. Population divergence was high for mtDNA variation (G_ST = 0.738, N_ST = 0.771) indicating low levels of seed‐based gene flow and significant phylogeographical structure (N_ST > G_ST, P < 0.05). The spatial distribution of mitotypes suggests that five distinct population groups exist in the species: one in the west comprising seven populations, a second with a north–central distribution comprising 15 populations, a third with a southern and easterly distribution comprising five populations, a fourth comprising one central and one western population, and a fifth comprising a single population located in the north‐central part of the species’ range. Each group apart from the fourth group is characterized by a distinct mitotype, with other mitotypes, if present, occurring at low frequency. It is suggested, therefore, that most members of each group apart from Group 4 are derived from ancestors that occupied different isolated refugia in a previous period of range fragmentation of the species, possibly at the time of the Last Glacial Maximum. Possible locations for these refugia are suggested. A comparison of mitotype diversity between northern and southern subgroups within the north‐central group of populations (Group 2) showed much greater uniformity in the northern part of the range both within and between populations. This could indicate a northward migration of the species from a southern refugium in this region during the postglacial period, although alternative explanations cannot be ruled out. Two chlorotypes were distributed across the geographical range of the species, resulting in lower levels of among‐population chlorotype variation. The geographical pattern of variation for all five chlorotypes provided some indication of the species surviving past glaciations in more than one refugium, although differentiation was much less marked, presumably due to the greater dispersal of cpDNA via pollen.