Subantarctic flowering plants: pre‐glacial survivors or post‐glacial immigrants?

Aim The aim here was to assess whether the present‐day assemblage of subantarctic flowering plants is the result of a rapid post‐Last Glacial Maximum (LGM) colonization or whether subantarctic flowering plants survived on the islands in glacial refugia throughout the LGM. Location The circumpolar subantarctic region, comprising six remote islands and island groups between latitudes 46° and 55° S, including South Georgia in the South Atlantic Ocean, the Prince Edward Islands, Îles Crozet, Îles Kerguelen, the Heard Island group in the South Indian Ocean and Macquarie Island in the South Pacific Ocean. Methods Floristic affinities between the subantarctic islands were assessed by cluster analysis applied to an up‐to‐date dataset of the phanerogamic flora in order to test for the existence of provincialism within the subantarctic. A review of the primary literature on the palaeobotany, geology and glacial history of the subantarctic islands was carried out and supplemented with additional palaeobotanical data and new field observations from South Georgia, Île de la Possession (Îles Crozet) and Îles Kerguelen. Results First, a strong regionalism was observed, with different floras characterizing the islands in each of the ocean basins, and endemic species being present in the South Indian Ocean and South Pacific Ocean provinces. Second, the majority of the plant species were present at the onset of accumulation of post‐glacial organic sediment and there is no evidence for the natural arrival of new immigrants during the subsequent period. Third, a review of geomorphological data suggested that the ice cover was incomplete during the LGM on the majority of the islands, and ice‐free biological refugia were probably present even on the most glaciated islands. Main conclusions Several independent lines of evidence favour the survival of a native subantarctic phanerogamic flora in local refugia during the LGM rather than a post‐LGM colonization from more distant temperate landmasses in the Southern Hemisphere.     

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.     

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

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

Palaeoclimatic models help to understand current distribution of Caucasian forest species

Spatial and temporal constraints on dispersal explain the absence of species from areas with potentially suitable conditions. Previous studies have shown that post‐glacial recolonization has shaped the current ranges of many species, yet it is not completely clear to what extent interspecific differences in range size depend on different dispersal rates. The inferred boundaries of glacial refugia are difficult to validate, and may bias spatial distribution models (SDMs) that consider post‐glacial dispersal constraints. We predicted the current distribution of 12 Caucasian forest plants and animals, factoring in the effective geographical distance from inferred glacial refugia as an additional predictor. To infer glacial refugia, we tested the transferability of the current SDMs based on the distribution of climatic variables, and projected the most transferable ones onto two climate scenarios simulated for the Last Glacial Maximum (LGM). We then calculated least‐cost distances from the inferred refugia, using elevation as a friction surface, and recalculated the current SDMs incorporating the distances as an additional variable. We compared the predictive powers of the initial with the final SDMs. The palaeoclimatic simulation that best matched the distribution of species was assumed to represent the closest fit to the true palaeoclimate. SDMs incorporating refugial distance performed significantly better for all but one studied species, and the Model for Interdisciplinary Research on Climate (MIROC) climatic simulation provided a more convincing pattern of the LGM climate than the Community Climate System Model (CCSM) simulation. Our results suggest that the projection of suitable habitat models onto past climatic conditions may yield realistic boundaries of glacial refugia, and that the current distribution of forest species in the study region is strongly associated with locations of former refugia. We inferred six major forest refugia throughout western Asia: (1) Colchis; (2) western Anatolia; (3) western Taurus; (4) the upper reaches of the Tigris River; (5) the Levant; and (6) the southern Caspian basin. The boundaries of the modelled refugia were substantially broader than the refugia boundaries inferred solely from pollen records. Thus, our method could be used to: (1) improve models of current species distributions by considering the dispersal histories of the species; and (2) validate alternative reconstructions of palaeoclimate with current distribution data. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 105 , 231–248.     

History and evolution of alpine plants endemic to the Qinghai-Tibetan Plateau: Aconitum gymnandrum (Ranunculaceae)

How Quaternary climatic oscillations affected range distributions and intraspecific divergence of alpine plants on the Qinghai‐Tibetan Plateau (QTP) remains largely unknown. Here, we report a survey of chloroplast DNA (cpDNA) and nuclear ribosomal internal transcribed spacer (ITS) DNA variation aimed at exploring the phylogeographical history of the QTP alpine endemic Aconitum gymnandrum. We sequenced three cpDNA fragments (rpl20–rps12 intergenic spacer, the trnV intron and psbA‐trnH spacer) and also the nuclear (ITS) region in 245 individuals from 23 populations sampled throughout the species’ range. Two distinct lineages, with eastern and western geographical distributions respectively, were identified from a phylogenetic analysis of ITS sequence variation. Based on a fast substitution rate, these were estimated to have diverged from each other in the early Pleistocene approximately 1.45 Ma. The analysis of cpDNA variation identified nine chlorotypes that clustered into two major clades that were broadly congruent in geographical distribution with the two ITS lineages. The east–west split of cpDNA divergence was supported by an amova which partitioned approximately half of the total variance between these two groups of populations. Analysis of the spatial distribution of chlorotypes showed that each clade was subdivided into two groups of populations such that a total of four population groups existed in the species. It is suggested that these different groups derive from four independent glacial refugia that existed during the Last Glacial Maximum (LGM), and that three of these refugia were located at high altitude on the QTP platform itself at that time. Coalescent simulation of chlorotype genealogies supported both an early Pleistocene origin of the two main cpDNA clades and also the ‘four‐refugia’ hypothesis during the LGM. Two previous phylogeographical studies of QTP alpine plants indicated that such plants retreated to refugia at the eastern/south‐eastern plateau edge during the LGM and/or previous glacial maxima. However, the results for A. gymnandrum suggest that at least some of these cold‐tolerant species may have also survived centrally on the QTP platform throughout the Quaternary.     

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

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

Rethinking long‐term vegetation dynamics: multiple glacial refugia and local expansion of a species complex

Evidence is accumulating that some arcto‐boreal plant taxa persisted through the last glacial maximum (LGM) in Alaska and adjacent Canada. However, the spatial patterns of glacial persistence and associated postglacial colonization remain largely unknown. In this study, we investigated the LGM refugia of an alder (Alnus) species complex (n = 3 taxa) and assess the spatiotemporal dynamics of Alnus in this vast region. Specifically, we conducted high‐throughput DNA sequencing (ddRADseq) on Alnus foliar samples collected from a dense population network to investigate patterns of genetic structure and infer the presence of glacial lineages. Species distribution modeling (SDM) was used to investigate the probability and possible locations of glacial persistence. These analyses were integrated and then compared with fossil pollen data to identify the locations of refugial populations and spatial patterns of postglacial colonization. Our genetic analyses revealed two glacial lineages with separate geographic origins for each Alnus taxon, suggesting that the genus persisted in multiple LGM refugia. Non‐overlapping hindcast distributions based on SDMs further support the presence of multiple, spatially distinct refugia. These ddRADseq and SDM results, in conjunction with reassessment of fossil pollen records, suggest that Alnus expanded from several population nuclei that existed during the LGM and coalesced during the Holocene to form its present range. These results challenge the unidirectional model for postglacial vegetation expansion, implying that climate buffering associated with landscape heterogeneity and adaptation to millennial‐scale environmental variability played important roles in driving late‐Quaternary population dynamics.     

Resilient forest faunal communities in South Africa: a legacy of palaeoclimatic change and extinction filtering?

Aim To examine the influence of climatic extinction filtering during the last glacial maximum (LGM; c. 18,000 yr bp ) and of the subsequent recolonization of forest faunas on contemporary assemblage composition in southern African forests. Location South Africa, Mozambique, Swaziland, Zimbabwe. Methods Data comprised presence/absence by quarter‐degree grid cell for forest‐dependent and forest‐associated birds, non‐volant mammals and frogs. Twenty‐one forest subregions were assigned to one of three previously identified forest types: Afrotemperate, scarp, and Indian Ocean coastal belt. Differences among forest types were examined through patterns and gradients of species richness and endemism, assemblage similarity, species turnover, and coefficients of species dispersal direction. The influence of contemporary environment on assemblage composition was investigated using partial canonical correspondence analysis. Several alternative biogeographical hypotheses for the recolonization of forest faunas were tested. Results Afrotemperate faunas are relatively species‐poor, have low species turnover, and are unsaturated and infiltrated by generalist species. In northern and central regions, communities are supplemented by recolonization from scarp forest refugia, and among frogs by autochthanous speciation in localized refugia. Scarp faunas are relatively species‐rich, contain many forest‐dependent species, have high species turnover, and overlap with coastal and Afrotemperate faunas. Coastal forests are relatively species‐rich with high species turnover. Main conclusions Afrotemperate communities were affected most by climatic extinction filtering events. Scarp forests were Afrotemperate refugia during the LGM and are a contemporary overlap zone between Afrotemperate and coastal forest. Coastal faunas derive from post‐LGM colonization along the eastern seaboard from tropical East African refugia. The greatest diversity is achieved in scarp and coastal forest faunas in northern KwaZulu–Natal province. This historical centre of diversity has influenced the faunal diversity of nearly all other forests in South Africa. The response of vertebrate taxa to large‐scale, historical processes is dependent on their relative mobility: forest birds best illustrate patterns resulting from post‐glacial faunal dispersal, while among mammals and frogs the legacy of climatic extinction filtering remains stronger.     

Chloroplast evidence for geographic stasis of the Australian bird‐dispersed shrub Tasmannia lanceolata (Winteraceae)

Few chloroplast‐based genetic studies have been undertaken for plants of mesic temperate forests in the southern hemisphere and fossil‐based models have provided evidence of vegetation history only at the broadest scales in this region. This study investigates the chloroplast DNA phylogeography of Tasmannia lanceolata (Winteraceae), a fleshy‐fruited, bird‐dispersed shrub that is widespread in the mountains of southeastern Australia and Tasmania. Thirty haplotypes were identified after sequencing 3206 bp of chloroplast DNA in each of 244 individuals collected across the species’ range. These haplotypes showed unexpectedly strong phylogeographic structuring, including a phylogeographic break within a continuous part of the species’ range, with the distribution of four major clades mostly not overlapping, and geographic structuring of haplotypes within these clades. This strong geographic patterning of chloroplast DNA provided evidence for the survival of T. lanceolata in multiple putative wet forest refugia as well as evidence for additional wet forest species refugia in southeastern Australia. In western Tasmania lower haplotype diversity below the LGM tree line compared to above the LGM tree line suggests that glacial refugia at high altitudes may have been important for T. lanceolata. The level of geographic structuring in T. lanceolata is similar to gravity dispersed southern hemisphere plants such as Nothofagus and Eucalyptus. Behavioural traits of the birds transporting seed may have had a strong bearing on the limited transport of T. lanceolata seed, although factors limiting establishment, possibly including selection, may also have been important.     

Mitochondrial DNA of pre‐last glacial maximum red deer from NW Spain suggests a more complex phylogeographical history for the species

The major climatic oscillations that characterized the Quaternary had a great influence on the evolution and distribution of several species. During cold periods, the distribution of temperate‐adapted species became fragmented with many surviving in southern refugia (Iberian, Italian, and Balkan Peninsulas). Red deer was one of the species that contracted its original range to southern refugia. Currently, two main lineages have been described for the species: western and eastern. We have analyzed fossils pre‐dating the last glacial maximum (LGM) from Liñares cave (NW Spain) that belongs to the peripheral range of the western clade, and fossils from the Danish Holocene belonging to the central part of the same clade. Phylogenetic analyses place our samples in the western clade. However, some specimens from Liñares represent an early split in the tree along with other pre‐LGM western samples from previous studies. Despite low bootstrap values in the Bayesian phylogenies, haplotype networks connect these foreign haplotypes to the eastern clade. We suggest a mixed phylogeographical model to explain this pattern with range expansions from the east during the expansion phase after the cold periods in marine isotope stage 3. We find slight isolation by distance in post‐LGM populations that could be a consequence of the recolonization from southern refugia after the LGM.     

Multilocus phylogeography and population structure of common eiders breeding in North America and Scandinavia

Aim Glacial refugia during the Pleistocene had major impacts on the levels and spatial apportionment of genetic diversity of species in northern latitude ecosystems. We characterized patterns of population subdivision, and tested hypotheses associated with locations of potential Pleistocene refugia and the relative contribution of these refugia to the post‐glacial colonization of North America and Scandinavia by common eiders (Somateria mollissima). Specifically, we evaluated localities hypothesized as ice‐free areas or glacial refugia for other Arctic vertebrates, including Beringia, the High Arctic Canadian Archipelago, Newfoundland Bank, Spitsbergen Bank and north‐west Norway. Location Alaska, Canada, Norway and Sweden. Methods Molecular data from 12 microsatellite loci, the mitochondrial DNA (mtDNA) control region, and two nuclear introns were collected and analysed for 15 populations of common eiders (n = 716) breeding throughout North America and Scandinavia. Population genetic structure, historical population fluctuations and gene flow were inferred using F‐statistics, analyses of molecular variance, and multilocus coalescent analyses. Results Significant inter‐population variation in allelic and haplotypic frequencies were observed (nuclear DNA F_ST = 0.004–0.290; mtDNA Φ_ST = 0.051–0.927). Whereas spatial differentiation in nuclear genes was concordant with subspecific designations, geographic proximity was more predictive of inter‐population variance in mitochondrial DNA haplotype frequency. Inferences of historical population demography were consistent with restriction of common eiders to four geographic areas during the Last Glacial Maximum: Belcher Islands, Newfoundland Bank, northern Alaska and Svalbard. Three of these areas coincide with previously identified glacial refugia: Newfoundland Bank, Beringia and Spitsbergen Bank. Gene‐flow and clustering analyses indicated that the Beringian refugium contributed little to common eider post‐glacial colonization of North America, whereas Canadian, Scandinavian and southern Alaskan post‐glacial colonization is likely to have occurred in a stepwise fashion from the same glacial refugium. Main conclusions Concordance of proposed glacial refugia used by common eiders and other Arctic species indicates that Arctic and subarctic refugia were important reservoirs of genetic diversity during the Pleistocene. Furthermore, suture zones identified at MacKenzie River, western Alaska/Aleutians and Scandinavia coincide with those identified for other Arctic vertebrates, suggesting that these regions were strong geographic barriers limiting dispersal from Pleistocene refugia.     

Strong paleoclimatic legacies in current plant functional diversity patterns across Europe

Numerous studies indicate that environmental changes during the late Quaternary have elicited long‐term disequilibria between species diversity and environment. Despite its importance for ecosystem functioning, the importance of historical environmental conditions as determinants of FD (functional diversity) remains largely unstudied. We quantified the geographic distributions of plant FD (richness and dispersion) across Europe using distribution and functional trait information for 2702 plant species. We then compared the importance of historical and contemporary factors to determine the relevance of past conditions as predictors of current plant FD in Europe. For this, we compared the strength of the relationships between FD with temperature and precipitation stability since the LGM (Last Glacial Maximum), accessibility to LGM refugia, and contemporary environmental conditions (climate, productivity, soil, topography, and land use). Functional richness and dispersion exhibited geographic patterns with strong associations to the environmental history of the region. The effect size of accessibility to LGM refugia and climate stability since the LGM was comparable to that of the contemporary predictors. Both functional richness and dispersion increased with temperature stability since the LGM and accessibility to LGM refugia. Functional richness' geographic pattern was primarily associated with accessibility to LGM refugia growing degree‐days, land use heterogeneity, diversity of soil types, and absolute minimum winter temperature. Functional dispersion's geographic pattern was primarily associated with accessibility to LGM refugia growing degree‐days and absolute minimum winter temperature. The high explained variance and model support of historical predictors are consistent with the idea that long‐term variability in environmental conditions supplements contemporary factors in shaping FD patterns at continental scales. Given the importance of FD for ecosystem functioning, future climate change may elicit not just short‐term shifts in ecosystem functioning, but also long‐term functional disequilibria.     

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.     

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

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

Likelihood‐based inference of population history from low‐coverage de novo genome assemblies

Short‐read sequencing technologies have in principle made it feasible to draw detailed inferences about the recent history of any organism. In practice, however, this remains challenging due to the difficulty of genome assembly in most organisms and the lack of statistical methods powerful enough to discriminate between recent, nonequilibrium histories. We address both the assembly and inference challenges. We develop a bioinformatic pipeline for generating outgroup‐rooted alignments of orthologous sequence blocks from de novo low‐coverage short‐read data for a small number of genomes, and show how such sequence blocks can be used to fit explicit models of population divergence and admixture in a likelihood framework. To illustrate our approach, we reconstruct the Pleistocene history of an oak‐feeding insect (the oak gallwasp Biorhiza pallida), which, in common with many other taxa, was restricted during Pleistocene ice ages to a longitudinal series of southern refugia spanning the Western Palaearctic. Our analysis of sequence blocks sampled from a single genome from each of three major glacial refugia reveals support for an unexpected history dominated by recent admixture. Despite the fact that 80% of the genome is affected by admixture during the last glacial cycle, we are able to infer the deeper divergence history of these populations. These inferences are robust to variation in block length, mutation model and the sampling location of individual genomes within refugia. This combination of de novo assembly and numerical likelihood calculation provides a powerful framework for estimating recent population history that can be applied to any organism without the need for prior genetic resources.     

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.

     

Phylogeography and niche modelling of the relict plant Amborella trichopoda (Amborellaceae) reveal multiple Pleistocene refugia in New Caledonia

Amborella trichopoda Baill. (Amborellaceae, Amborellales), the sole living member of the sister group to all other extant angiosperms, is endemic to New Caledonia. We addressed the intraspecific phylogeography of Amborella by investigating whether its present population genetic structure could be related to its current and past habitats. We found moderate range‐wide genetic diversity based on nuclear microsatellite data and detected four well‐differentiated, geographically distinct genetic groups using Bayesian clustering analyses. We modelled the ecological niche of Amborella based on the current climatic and environmental conditions. The predictive ability of the model was very good throughout the Central East mainland zone, but Amborella was predicted in the northern part of the island where this plant has not been reported. Furthermore, no significant barrier was detected based on habitat suitability that could explain the genetic differentiation across the area. Conversely, we found that the main genetic clusters could be related to the distribution of the suitable habitat at the last glacial maximum (LGM, c. 21 000 years BP), when Amborella experienced a dramatic 96.5% reduction in suitable area. At least two lineages survived in distinct putative refugia located in the Massif des Lèvres and in the vicinity of Mount Aoupinié. Our findings finally confirmed the importance of LGM rainforest refugia in shaping the current intra‐ and interspecific diversity in New Caledonian plants and revealed the possibility of an as yet unreported refugium. The combination of niche modelling and population genetics thereby offered novel insight into the biogeographical history of an emblematic taxon.     

Genome‐wide set of SNPs reveals evidence for two glacial refugia and admixture from postglacial recolonization in an alpine ungulate

Past glaciation events have played a major role in shaping the genetic diversity and distribution of wild sheep in North America. The advancement of glaciers can isolate populations in ice‐free refugia, where they can survive until the recession of ice sheets. The major Beringian refugium is thought to have held thinhorn sheep (Ovis dalli) populations during times of glacial advance. While isolation in the major refugium can account for much of the genetic and morphological diversity seen in extant thinhorn sheep populations, mounting evidence suggests the persistence of populations in smaller minor refugia. We investigated the refugial origins of thinhorn sheep using ~10 000 SNPs obtained via a cross‐species application of the domestic sheep ovine HD BeadChip to genotype 52 thinhorn sheep and five bighorn sheep (O. canadensis) samples. Phylogenetic inference revealed a distinct lineage of thinhorn sheep inhabiting British Columbia, which is consistent with the survival of a group of thinhorn sheep in a minor refugium separate from the Beringian refugium. Isolation in separate glacial refugia probably mediated the evolution of the two thinhorn sheep subspecies, the white Dall's sheep (O. d. dalli), which persisted in Beringia, and the dark Stone's sheep (O. d. stonei), which utilized the minor refugium. We also found the first genetic evidence for admixture between sheep from different glacial refugia in south‐central Yukon as a consequence of post glacial expansion and recolonization. These results show that glaciation events can have a major role in the evolution of species inhabiting previously glaciated habitats and the need to look beyond established refugia when examining the evolutionary history of such species.     

Glacial refugia of mammals in Europe: evidence from fossil records

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

     

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

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