Phylogeography of Pulsatilla vernalis (L.) Mill. (Ranunculaceae): chloroplast DNA reveals two evolutionary lineages across central Europe and Scandinavia

Aim The aim of this study was to test hypotheses regarding some of the main phylogeographical patterns proposed for European plants, in particular the locations of glacial refugia, the post‐glacial colonization routes, and genetic affinities between southern (alpine) and northern (boreal) populations. Location The mountains of Europe (Alps, Balkans, Carpathians, Central Massif, Pyrenees, Scandinavian chain, Sudetes), and central European/southern Scandinavian lowlands. Methods As our model system we used Pulsatilla vernalis, a widely distributed European herbaceous plant occurring both in the high‐mountain environments of the Alps and other European ranges and in lowlands north of these ranges up to Scandinavia. Based on a distribution‐wide sampling of 61 populations, we estimated chloroplast DNA (cpDNA) variation along six regions using polymerase chain reaction–restriction fragment‐length polymorphisms (PCR–RFLPs) (trnH–trnK, trnK–trnK, trnC–trnD, psbC–trnS, psaA–trnS, trnL–trnF) and further sequencing of trnL–trnF and trnH–psbA. In addition, 11 samples of other European species of Pulsatilla were sequenced to survey the genus‐scale cpDNA variation. Results Eleven PCR–RFLP polymorphisms were detected in P. vernalis, revealing seven haplotypes. They formed two distinct genetic groups. Three haplotypes representing both groups dominated and were widely distributed across Europe, whereas the others were restricted to localized regions (central Alps, Tatras/Sudetes mountains) or single populations. Sequencing analysis confirmed the reliability of PCR–RFLPs and homology of haplotypes across their distribution. The chloroplast DNA variation across the section Pulsatilla was low, but P. vernalis did not share haplotypes with other species. Main conclusions The genetic distinctiveness of P. vernalis populations from the south‐western Alps with respect to other Alpine populations, as well as the affinities between the former populations and those from the eastern Pyrenees, is demonstrated, thus providing support for the conclusions of previous studies. Glacial refugia in the Dolomites are also suggested. Isolation is inferred for the high‐mountain populations from the Tatras and Sudetes; this is in contrast to the case for the Balkans, which harboured the common haplotype. Specific microsatellite variation indicates the occurrence of periglacial lowland refugia north of the Alps, acting as a source for the post‐glacial colonization of Scandinavia. The presence of different fixed haplotypes in eastern and western Scandinavia, however, suggests independent post‐glacial colonization of these two areas, with possible founder effects.     

Phylogeographic analyses suggest that a deciduous species (Ostryopsis davidiana Decne., Betulaceae) survived in northern China during the Last Glacial Maximum

Aim Palaeontologial data suggest that all temperate forest species in northern China migrated southwards during the Last Glacial Maximum (LGM) and recolonized post‐glacially within the Holocene. We tested this assumption using phylogeographical studies of a temperate deciduous shrub species (Ostryopsis davidiana Decne., Betulaceae), which has a wide distribution in northern China. Location Northern China. Methods We sequenced two chloroplast DNA (cpDNA) fragments (trnL–trnF and psbA–trnH, together about 1300 bp in length) of 294 plants from 21 populations across the total distribution range of this species. We used maximum parsimony and haplotype network methods to construct phylogenetic relationships among haplotypes. Results The analysis of cpDNA variation identified eight haplotypes. A single haplotype was fixed in all populations except for one population that was polymorphic, having two haplotypes. The population subdivisions were extremely high (G_ST = 0.972 and N_ST = 0.974), suggesting very low gene flow between populations. Haplotypes clustered into two tentative clades, both of which occur in the southern region of the species’ range but only one of which occurs in the northern region. Across the sampled populations, the haplotype distributions were differentiated geographically. Main conclusions Our analyses suggest that multiple refugia were maintained across the range of O. davidiana in both northern (north of the Qing Mountains) and southern (south of the Qing Mountains) regions during the LGM rather than that the species survived only in the south and subsequently colonized northwards. The extremely low within‐population diversity of this species suggests strong bottleneck or founder effects within each fragmented region during the Quaternary climatic oscillations. These findings provide important clues for understanding range shifts and changes in within‐ and/or between‐population genetic diversity of temperate forests in response to past climatic oscillations in northern China.     

Population genomic evidence for plant glacial survival in Scandinavia

Quaternary glaciations have played a major role in shaping the genetic diversity and distribution of plant species. Strong palaeoecological and genetic evidence supports a postglacial recolonization of most plant species to northern Europe from southern, eastern and even western glacial refugia. Although highly controversial, the existence of small in situ glacial refugia in northern Europe has recently gained molecular support. We used genomic analyses to examine the phylogeography of a species that is critical in this debate. Carex scirpoidea Michx subsp. scirpoidea is a dioecious, amphi‐Atlantic arctic–alpine sedge that is widely distributed in North America, but absent from most of Eurasia, apart from three extremely disjunct populations in Norway, all well within the limits of the Weichselian ice sheet. Range‐wide population sampling and variation at 5,307 single nucleotide polymorphisms show that the three Norwegian populations comprise unique evolutionary lineages divergent from Greenland with high between‐population divergence. The Norwegian populations have low within‐population genetic diversity consistent with having experienced genetic bottlenecks in glacial refugia, and host private alleles that probably accumulated in long‐term isolated populations. Demographic analyses support a single, pre‐Weichselian colonization into Norway from East Greenland, and subsequent divergence of the three populations in separate refugia. Other refugial areas are identified in North‐east Greenland, Minnesota/Michigan, Colorado and Alaska. Admixed populations in British Columbia and West Greenland indicate postglacial contact. Taken together, evidence from this study strongly indicates in situ glacial survival in Scandinavia.     

Genetic consequences of glacial survival and postglacial colonization in Norway spruce: combined analysis of mitochondrial DNA and fossil pollen

Norway spruce (Picea abies [L.] Karst.) is a broadly distributed European conifer tree whose history has been intensively studied by means of fossil records to infer the location of full‐glacial refugia and the main routes of postglacial colonization. Here we use recently compiled fossil pollen data as a template to examine how past demographic events have influenced the species’ modern genetic diversity. Variation was assessed in the mitochondrial nad1 gene containing two minisatellite regions. Among the 369 populations (4876 trees) assayed, 28 mitochondrial variants were identified. The patterns of population subdivision superimposed on interpolated fossil pollen distributions indicate that survival in separate refugia and postglacial colonization has led to significant structuring of genetic variation in the southern range of the species. The populations in the northern range, on the other hand, showed a shallow genetic structure consistent with the fossil pollen data, suggesting that the vast northern range was colonized from a single refugium. Although the genetic diversity decreased away from the putative refugia, there were large differences between different colonization routes. In the Alps, the diversity decreased over short distances, probably as a result of population bottlenecks caused by the presence of competing tree species. In northern Europe, the diversity was maintained across large areas, corroborating fossil pollen data in suggesting that colonization took place at high population densities. The genetic diversity increased north of the Carpathians, probably as a result of admixture of expanding populations from two separate refugia.     

The role of disjunction and postglacial population expansion on phylogeographical history and genetic diversity of the circumboreal plant Chamaedaphne calyculata

In the last decade a number of studies has illustrated quite different phylogeographical patterns amongst plants with a northern present‐day geographical distribution, spanning the entire circumboreal region and/or circumarctic region and southern mountains. These works, employing several marker systems, have brought to light the complex evolutionary histories of this group. Here I focus on one circumboreal plant species, Chamaedaphne calyculata (leatherleaf), to unravel its phylogeographical history and patterns of genetic diversity across its geographical range. A survey of 29 populations with combined analyses of chloroplast DNA (cpDNA), internal transcribed spacer (ITS) and AFLP markers revealed structuring into two groups: Eurasian/north‐western North American, and north‐eastern North American. The present geographical distribution of C. calyculata has resulted from colonization from two putative refugial areas: east Beringia and south‐eastern North America. The variation of chloroplast DNA (cpDNA) and ITS sequences strongly indicated that the evolutionary histories of the Eurasian/north‐western North American and the north‐eastern North American populations were independent of each other because of a geographical disjunction in the distribution area and ice‐sheet history between north‐eastern and north‐western North America. Mismatch analysis using ITS confirmed that the present‐day population structure is the result of rapid expansion, probably since the last glacial maximum. The AFLP data revealed low genetic diversity of C. calyculata (P = 19.5%, H = 0.085) over the whole geographical range, and there was no evidence of loss of genetic diversity within populations in the continuous range, either at the margins or in formerly glaciated and nonglaciated regions. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 105 , 761–775.     

Contrasting patterns of nuclear microsatellite genetic structure of Fraxinus mandshurica var. japonica between northern and southern populations in Japan

Aim The aim of this study is to detect extant patterns of population genetic structure of Fraxinus mandshurica var. japonica in Japan, and to provide insights into the post‐glacial history of this species during the Holocene. Location Hokkaido and Honshu islands, Japan (including the Oshima and Shimokita peninsulas). Methods We examined nine polymorphic nuclear microsatellite loci to assess genetic variation within and among 15 populations across almost the entire range of the species in Japan. Extant patterns of geographical structure were analysed using Bayesian clustering, Monmonier’s algorithm, analysis of molecular variance, Mantel tests and principal coordinates analysis. Recent bottlenecks within populations and regional genetic variation were also assessed. Results Northern populations (Hokkaido Island and the Shimokita Peninsula) formed a single homogeneous deme, maintaining the highest level of allelic diversity on the Oshima Peninsula. By contrast, southern populations (Honshu Island) demonstrated strong substructure on both coasts. Specifically, populations on the Pacific side of Honshu exhibited significant bottlenecks and erosion of allelic diversity but preserved distinct subclusters diverging from widespread subclusters on the Japan Sea side of this island. Main conclusions Genetic evidence and life history traits suggest that F. mandshurica occupied cryptic northern refugia on the Oshima Peninsula during the Last Glacial Maximum, which is reflected in the species’ extant northern distribution. Strong geographical structure in southern populations, in agreement with fossil pollen records, suggests geographical isolation by mountain ranges running north–south along Honshu. Given that this tree species is cold‐adapted and found in riparian habitats, populations on the Pacific side of Honshu probably contracted into higher‐elevation swamps during warm post‐glacial periods, leading to a reduction of effective population sizes and rare allelic richness.     

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.     

Multiple refugia and barriers explain the phylogeography of the Valais shrew,Sorex antinorii (Mammalia: Soricomorpha)

The aim of the present study was to investigate the genetic structure of the Valais shrew (Sorex antinorii) by a combined phylogeographical and landscape genetic approach, and thereby to infer the locations of glacial refugia and establish the influence of geographical barriers. We sequenced part of the mitochondrial cytochrome b (cyt b) gene of 179 individuals of S. antinorii sampled across the entire species' range. Six specimens attributed to S. arunchi were included in the analysis. The phylogeographical pattern was assessed by Bayesian molecular phylogenetic reconstruction, population genetic analyses, and a species distribution modelling (SDM)‐based hindcasting approach. We also used landscape genetics (including isolation‐by‐resistance) to infer the determinants of current intra‐specific genetic structure. The phylogeographical analysis revealed shallow divergence among haplotypes and no clear substructure within S. antinorii. The starlike structure of the median‐joining network is consistent with population expansion from a single refugium, probably located in the Apennines. Long branches observed on the same network also suggest that another refugium may have existed in the north‐eastern part of Italy. This result is consistent with SDM, which also suggests several habitable areas for S. antinorii in the Italian peninsula during the LGM. Therefore S. antinorii appears to have occupied disconnected glacial refugia in the Italian peninsula, supporting previous data for other species showing multiple refugia within southern refugial areas. By coupling genetic analyses and SDM, we were able to infer how past climatic suitability contributed to genetic divergence of populations. The genetic differentiation shown in the present study does not support the specific status of S. arunchi. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 105 , 864–880.     

Genetically diverse but with surprisingly little geographical structure: the complex history of the widespread herb Carex nigra (Cyperaceae)

Aim Based on extensive range‐wide sampling, we address the phylogeographical history of one of the most widespread and taxonomically complex sedges in Europe, Carex nigra s. lat. We compare the genetic structure of the recently colonized northern areas (front edge) and the long‐standing southern areas (rear edge), and assess the potential genetic basis of suggested taxonomic divisions at the rank of species and below. Location Amphi‐Atlantic, central and northern Europe, circum‐Mediterranean mountain ranges, central Siberia, Himalayas. Methods A total of 469 individuals sampled from 83 populations, covering most of the species’ range, were analysed with amplified fragment length polymorphism (AFLP) and chloroplast DNA (cpDNA) markers. Bayesian clustering, principal coordinates analysis, and estimates of diversity and differentiation were used for the analysis of AFLP data. CpDNA data were analysed with statistical parsimony networks and maximum parsimony and Bayesian inference of phylogenetic trees. Results Overall genetic diversity was high, but differentiation among populations was limited. Major glacial refugia were inferred in the Mediterranean Basin and in western Russia; in addition, there may have been minor refugia in the North Atlantic region. In the southern part of the range, we found high levels, but geographically quite poorly structured genetic diversity, whereas the levels of genetic diversity varied among different areas in the north. North American populations were genetically very similar to the European populations. Main conclusions The data are consistent with extensive gene flow, which has obscured the recent history of the taxon. The limited differentiation in the south probably results from the mixing of lineages expanding from several local refugia. Northward post‐glacial colonization resulted in a leading‐edge pattern of low diversity in the Netherlands, Belgium, Scotland and Iceland, whereas the observed high diversity levels in Fennoscandia suggest broad‐fronted colonization from the south as well as from the east. The patterns found in the American populations are consistent with post‐glacial colonization, possibly even with anthropogenic introduction from Europe. Our data also suggest that the tussock‐forming populations of C. nigra, often referred to as a distinct species (Carex juncella), represent an ecotype that has originated repeatedly from different populations with creeping rhizomes.     

Intervarietal and intravarietal genetic structure in Douglas‐fir: nuclear SSRs bring novel insights into past population demographic processes,phylogeography, and intervarietal hybridization

Douglas‐fir (Pseudotsuga menziesii) is one of numerous wide‐range forest tree species represented by subspecies/varieties, which hybridize in contact zones. This study examined the genetic structure of this North American conifer and its two hybridizing varieties, coastal and Rocky Mountain, at intervarietal and intravarietal level. The genetic structure was subsequently associated with the Pleistocene refugial history, postglacial migration and intervarietal hybridization/introgression. Thirty‐eight populations from the USA and Canada were genotyped for 13 nuclear SSRs and analyzed with simulations and traditional population genetic structuring methods. Eight genetic clusters were identified. The coastal clusters embodied five refugial populations originating from five distinct refugia. Four coastal refugial populations, three from California and one from western Canada, diverged during the Pleistocene (56.9–40.1 ka). The three Rocky Mountain clusters reflected distinct refugial populations of three glacial refugia. For Canada, ice covered during the Last Glacial Maximum, we present the following three findings. (1) One refugial population of each variety was revealed in the north of the distribution range. Additional research including paleodata is required to support and determine whether both northern populations originated from cryptic refugia situated south or north of the ice‐covered area. (2) An interplay between intravarietal gene flow of different refugial populations and intervarietal gene flow by hybridization and introgression was identified. (3) The Canadian hybrid zone displayed predominantly introgressants of the Rocky Mountain into the coastal variety. This study provides new insights into the complex Quaternary dynamics of this conifer essential for understanding its evolution (outside and inside the native range), adaptation to future climates and for forest management.     

Complex fine‐scale phylogeographical patterns in a putative refugial region for Fagus sylvatica (Fagaceae)

Broad‐scale plastid (chloroplast) DNA studies of beech (Fagus sylvatica) populations suggest the existence of glacial refugia and introgression zones in south‐eastern Europe. We choose a possible refugium of beech in northern Greece, Mt. Paggeo, which hosts a private plastid haplotype for beech, to conduct a fine‐scale genetic study. We attempt to confirm or reject the hypothesis of the existence of a small‐scale refugium and to gain an understanding of the ecological and topographical factors affecting the spatial distribution of plastid haplotypes in the area. Our results reveal a high haplotype diversity on Mt. Paggeo, but the overall distribution of haplotypes shows no significant correlation with the ecological characteristics of the beech forests. However, the private haplotype is found at high frequencies in beech forests located in or near ravines, having a high spatial overlap with a relict vegetation type occurring in ecological conditions found mainly in ravines. This result emphasizes the importance of topography in the existence of glacial refugia in the wider area. Furthermore, haplotypes originating from two more widespread beech lineages in Greece are found on Mt. Paggeo, indicating a possible mixing of populations originating from a local refugium with populations from remote refugia that possibly migrated into the area after the last glaciation. © 2014 The Linnean Society of London, Botanical Journal of the Linnean Society, 2014, 174 , 516–528.     

Identification of glacial refugia in south-eastern North America by phylogeographical analyses of a forest understorey plant, Trillium cuneatum

Aim We examine several hypotheses emerging from biogeographical and fossil records regarding glacial refugia of a southern thermophilic plant species. Specifically, we investigated the glacial history and post‐glacial colonization of a forest understorey species, Trillium cuneatum. We focused on the following questions: (1) Did T. cuneatum survive the Last Glacial Maximum (LGM) in multiple refugia, and (if so) where were they located, and is the modern genetic structure congruent with the fossil record‐based reconstruction of refugia for mesic deciduous forests? (2) What are the post‐glacial colonization patterns in the present geographical range? Location South‐eastern North America. Methods We sampled 45 populations of T. cuneatum throughout its current range. We conducted phylogeographical analyses based on maternally inherited chloroplast DNA (cpDNA haplotypes) and used TCS software to reconstruct intraspecific phylogeny. Results We detected six cpDNA haplotypes, geographically highly structured into non‐overlapping areas. With one exception, none of the populations had mixed haplotype composition. TCS analysis resulted in two intraspecific cpDNA lineages, with one clade subdivided further by shallower diversification. Main conclusions Our investigation revealed that T. cuneatum survived the LGM in multiple refugia, belonging to two (western, eastern) genealogical lineages geographically structured across south‐eastern North America. The western clade is confined to the south‐western corner of T. cuneatum’s modern range along the Lower Mississippi Valley, where fossil records document a major refugium of mesic deciduous forest. For the eastern clade, modern patterns of cpDNA haplotype distribution suggest cryptic vicariance, in the form of forest contractions and subsequent expansions associated with Pleistocene glacial cycles, rather than simple southern survival and subsequent northward colonization. The north–south partitioning of cpDNA haplotypes was unexpected, suggesting that populations of this rather southern thermophilic species may have survived in more northern locations than initially expected based on LGM climate reconstruction, and that the Appalachian Mountains functioned as a barrier to the dispersal of propagules originating in more southern refugia. Furthermore, our results reveal south‐west to north‐east directionality in historical migration through the Valley and Ridge region of north‐west Georgia.     

Multiple glacial refugia and postglacial colonization routes inferred for a woodland geophyte,Cyclamen purpurascens: patterns concordant with the Pleistocene history of broadleaved and coniferous tree species

The phylogeographical history of mid‐altitude woodland herbs that depend on moist and shaded forest habitats is poorly understood. Here, we analysed the genetic structure of Cyclamen purpurascens, a mountainous calcicolous perennial, to test hypotheses regarding its glacial survival in single or multiple refugia and postglacial colonization routes, and to explore how they are congruent with the histories inferred for temperate trees and other mountainous herbs. We gathered AFLP data and chloroplast DNA sequences (trnD‐trnT region) from 68 populations spanning the entire distribution range (the Jura Mountains, Alps, western Carpathians, Dinarides). Both genetic markers revealed two main phylogeographical groups (phylogroups) in C. purpurascens. Additionally, AFLP data detected a more detailed structure of five phylogroups: two widespread, showing east?west geographical separation, and three local ones, restricted to somewhat disjunct, marginal regions of the species range. We suggest that C. purpurascens survived the last glaciation in two main regions, the foothills of the Southern Limestone Alps and the Karst area of the north‐western Dinarides, and possibly also in microrefugia in the Western Carpathians. The glacial persistence and colonization routes of this woodland herb are highly concordant with those inferred for several temperate trees, especially the European beech. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 105 , 741–760.     

Genetic and morphological variation in a Mediterranean glacial refugium: evidence from Italian pygmy shrews,Sorex minutus (Mammalia: Soricomorpha)

At the Last Glacial Maximum (LGM), the southern European peninsulas were important refugia for temperate species. Current genetic subdivision of species within these peninsulas may reflect past population subdivision at the LGM, as in ‘refugia within refugia’, and/or at other time periods. In the present study, we assess whether pygmy shrew populations from different regions within Italy are genetically and morphologically distinct. One maternally and two paternally inherited molecular markers (cytochrome b and Y‐chromosome introns, respectively) were analysed using several phylogenetic methods. A geometric morphometric analysis was performed on mandibles to evaluate size and shape variability between populations. Mandible shape was also explored with a functional approach that considered the mandible as a first‐order lever affecting bite force. We found genetically and morphologically distinct European, Italian, and southern Italian groups. Mandible size increased with decreasing latitude and southern Italian pygmy shrews exhibited mandibles with the strongest bite force. It is not clear whether or not the southern Italian and Italian groups of pygmy shrews occupied different refugia within the Italian peninsula at the LGM. It is likely, however, that geographic isolation earlier than the LGM on islands at the site of present‐day Calabria was important in generating the distinctive southern Italian group of pygmy shrews, and also the genetic groups in other small vertebrates that we review here. Calabria is an important hotspot for genetic diversity, and is worthy of conservation attention. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 774–787.     

Phylogeography of Douglas‐fir based on mitochondrial and chloroplast DNA sequences: testing hypotheses from the fossil record

The integration of fossil and molecular data can provide a synthetic understanding of the ecological and evolutionary history of an organism. We analysed range‐wide maternally inherited mitochondrial DNA and paternally inherited chloroplast DNA sequence data with coalescent simulations and traditional population genetic methods to test hypotheses of population divergence generated from the fossil record of Douglas‐fir (Pseudotsuga menziesii), an ecologically and economically important western North American conifer. Specifically, we tested (i) the hypothesis that the Pliocene orogeny of the Cascades and Sierra Nevada caused the divergence of coastal and Rocky Mountain Douglas‐fir varieties; and (ii) the hypothesis that multiple glacial refugia existed on the coast and in the Rocky Mountains. We found that Douglas‐fir varieties diverged about 2.11 Ma (4.37 Ma–755 ka), which could be consistent with a Pliocene divergence. Rocky Mountain Douglas‐fir probably resided in three or more glacial refugia. More variable molecular markers would be required to detect the two coastal refugia suggested in the fossil record. Comparison of mitochondrial DNA and chloroplast DNA variation revealed that gene flow via pollen linked populations isolated from seed exchange. Postglacial colonization of Canada from coastal and Rocky Mountain refugia near the ice margin at the Last Glacial Maximum produced a wide hybrid zone among varieties that formed almost exclusively by pollen exchange and chloroplast DNA introgression, not seed exchange. Postglacial migration rates were 50–165 m/year, insufficient to track projected 21st century warming in some regions. Although fossil and genetic data largely agree, each provides unique insights.     

Multiple glacial refugia for cool‐temperate deciduous trees in northern East Asia: the Mongolian oak as a case study

In East Asia, temperate forests are predicted to have retracted southward to c. 30° N during the last glacial maximum (LGM) based on fossil pollen data, whereas phylogeographic studies have often suggested glacial in situ survival of cool‐temperate deciduous trees in their modern northern ranges. Here we report a study of the genetic diversity and structure of 29 natural Mongolian oak (Quercus mongolica) populations using 19 nuclear simple sequence repeat (nSSR) loci and four chloroplast DNA fragments. Bayesian clustering analysis with nSSRs revealed five groups, which were inferred by approximate Bayesian computation (ABC) to have diverged in multiple refugia through multiple glacial–interglacial cycles. Analysis of chloroplast DNA variation revealed four lineages that were largely but incompletely geographically disjunct. Ecological niche modelling (ENMs) indicated a southward range shift of the oak's distribution at the LGM, although high suitability scores were also evident in the Changbai Mts. (Northeast China), the Korean Peninsula, areas surrounding the Bohai Sea, and along the coast of the Russian Far East. In addition, endemic chloroplast DNA haplotypes and nuclear lineages occurred in high‐latitude northern areas where the ENM predicted no suitable habitat. The combined evidence from nuclear and chloroplast DNA, and the results of the ENM clearly demonstrate that multiple northern refugia, including cryptic ones, were maintained across the current distributional range of the Mongolian oak during the LGM or earlier glacial periods. Though spatially limited, postglacial expansions from these refugia have led to a pattern of decreased genetic diversity with increasing latitude.     

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

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

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.     

Rhytidium rugosum (Bryophyta) colonized Scandinavia from at least two glacial refugial source populations

The Scandinavian post‐glacial history of the moss Rhytidium rugosum is traced on the basis of information from the nuclear markers ITS and gpd for 229 Scandinavian and 81 other specimens. Some haplotypes, groups or lineages identified in a NeighborNet split network are predominantly northern Scandinavian, whereas others are southern. With the distributions of individual haplotypes and the timing of the deglaciation in different parts of Scandinavia, this implies colonization from the south and from the north or north‐east. High haplotype and nucleotide diversity and the occurrence of certain private haplotypes in the north suggest that the species may have survived the Last Glacial Maximum in local refugia. Slightly higher numbers of private haplotypes in Scandinavia than in central or north‐eastern Europe also favour an explanation with at least some local glacial survival. Low diversity in the southernmost contiguous region of the Scandinavian mountain range is probably a result of recent land uplift and late colonization. The Scandinavian lowland regional populations probably represent remains of an earlier widespread population that became increasingly restricted to small and isolated areas when the vegetation closed during the post‐glacial period. Some of the lowland populations require extensive management to survive. © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 179 , 635–657.     

Molecular markers provide evidence for a broad‐fronted recolonisation of the widespread calcareous grassland species Sanguisorba minor from southern and cryptic northern refugia

• Calcareous grasslands belong to the most species‐rich and endangered habitats in Europe. However, little is known about the origin of the species typically occurring in these grasslands. In this study we analysed the glacial and post‐glacial history of Sanguisorba minor, a typical plant species frequently occurring in calcareous grasslands.
• The study comprised 38 populations throughout the whole distribution range of the species across Europe. We used molecular markers (AFLP) and applied Bayesian cluster analysis as well as spatial principal components analysis (sPCA) to identify glacial refugia and post‐glacial migration routes to Central Europe.
• Our study revealed significant differences in the level of genetic variation and the occurrence of rare fragments within populations of S. minor and a distinct separation of eastern and western lineages. The analyses uncovered traditional southern but also cryptic northern refugia and point towards a broad fronted post‐glacial recolonisation.
• Based on these results we postulate that incomplete lineage sorting may have contributed to the detected pattern of genetic variation and that S. minor recolonised Central Europe post‐glacially from Iberia and northern glacial refugia in France, Belgium or Germany. Our results highlight the importance of refugial areas for the conservation of intraspecific variation in calcareous grassland species.