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.     

Putative glacial refugia of Cedrus atlantica deduced from Quaternary pollen records and modern genetic diversity

Aim To investigate the impact of past environmental changes on Cedrus atlantica and its current genetic diversity, and to predict its future distribution. Location Morocco, Algeria and Tunisia. Methods Eleven fossil pollen records from these three countries were used to locate putative glacial refugia and to reconstruct past climate changes. A mechanistic vegetation distribution model was used to simulate the distribution of C. atlantica in the year 2100. In addition, a genetic survey was carried out on modern Moroccan C. atlantica. Results Pollen records indicate that Cedrus was present during the last glacial period, probably in scattered refugia, in Tunisia, Algeria and Morocco. In the Tunisian and Algerian sites, cedar expanded during the late glacial and the early Holocene, then disappeared after c. 8000 yr bp . Reconstructed mean annual precipitation and January temperature show that the last glacial period in Morocco was cooler by 10–15°C and drier by c. 300–400 mm year^?1 than the climate today. Modern chloroplast microsatellites of 15 C. atlantica populations in Morocco confirm the presence of multiple refugia and indicate that cedar recolonized the Moroccan mountains fairly recently. Model simulation indicates that by the year 2100 the potential distribution of C. atlantica will be much restricted with a potential survival area located in the High Atlas. Main conclusions Environmental changes in northern Africa since the last glacial period have had an impact on the geographical distribution of C. atlantica and on its modern genetic diversity. It is possible that by the end of this century C. atlantica may be unable to survive in its present‐day locations. To preserve the species, we suggest that seedlings from modern C. atlantica populations located in the High Atlas mountains, where a high genetic diversity is found, be transplanted into the western High Atlas.     

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

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

Phylogeography of the North American red fox: vicariance in Pleistocene forest refugia

Fossil, archaeological, and morphometric data suggest that indigenous red foxes in North America were derived from vicariance in two disjunct refugia during the last glaciation: one in Beringia and one in the contiguous USA. To test this hypothesis, we conducted a phylogeographical analysis of the North American red fox within its presettlement range. We sequenced portions of the mitochondrial cytochrome b (354 bp) gene and D-loop (342 bp) from 220 historical red fox specimens. Phylogenetic analysis of the cytochrome b gene produced two clades that diverged c . 400 000 years before present ( bp ): a Holarctic and a Nearctic clade. D-loop analyses of the Nearctic clade indicated three distinct subclades (≥ 99% Bayesian posterior probability); two that were more recently derived (rho estimate c . 20 000 bp ) and were restricted to the southwestern mountains and the eastern portion of North America, and one that was older (rho estimate c . 45 000 bp ) and more widespread in North America. Populations that migrated north from the southern refugium following deglaciation were derived from the colonization of North America during or prior to the Illinoian glaciation (300 000–130 000 bp ), whereas populations that migrated south from the northern refugium represent a more recent colonization event during the Wisconsin glaciation (100 000–10 000 bp ). Our findings indicate that Nearctic clade red foxes are phylogenetically distinct from their Holarctic counterparts, and reflect long-term isolation in two disjunct forest refugia during the Pleistocene. The montane lineage, which includes endangered populations, may be ecologically and evolutionarily distinct.     

Imprints of glacial refugia in the modern genetic diversity of Pinus sylvestris

Aim To understand the impact of glacial refugia and migration pathways on the modern genetic diversity of Pinus sylvestris. Location The study was carried out throughout Europe. Methods An extended set of data of pollen and macrofossil remains was used to locate the glacial refugia and reconstruct the migrating routes of P. sylvestris throughout Europe. A vegetation model was used to simulate the extent of the potential refugia during the last glacial period. At the same time a genetic survey was carried out on this species. Results The simulated distribution of P. sylvestris during the last glacial period is coherent with the observed fossil data, which showed a patchy distribution of the refugia between c. 40° N and 50° N. Several migrational fronts were detected within the Iberian and the Italian peninsulas, and outside the Hungarian plain and around the Alps. The modern mitochondrial DNA depicted three different haplotypes for P. sylvestris. Two distinct haplotypes were restricted to northern Spain and Italy, and the third haplotype dominated most of the present‐day remaining distribution range of P. sylvestris in Europe. Main conclusions During the last glacial period P. sylvestris was constrained under severe climatic conditions to survive in scattered and restricted refugial areas. Combining palaeoenvironmental data, vegetation modelling and the genetic data, we have shown that the long‐term isolation in the glacial refugia and the migrational process during the Holocene have played a major role in shaping the modern genetic diversity of P. sylvestris in Europe.     

The roles of rivers and Pleistocene refugia in shaping genetic diversity in Praomys misonnei in tropical Africa

Aim This study aims to elucidate the phylogeography of the murid rodent Praomys misonnei and to document whether or not rain forest refugia and rivers structure patterns of diversity within this species. Location Tropical Africa, from Ghana to Kenya. Methods Patterns of genetic structure and signatures of population history (cytochrome b gene) were assessed in a survey of 229 individuals from 54 localities. Using maximum likelihood, Bayesian, network and genetic structure analyses, we inferred intra‐specific relationships and tested hypotheses for historical patterns of gene flow within P. misonnei. Results Our phylogenetic analyses reveal a strong phylogeographical structure. We identified four major geographical clades within P. misonnei: one clade in Ghana and Benin, a Nigerian clade, a West Central African clade and a Central and East African clade. Several subclades were identified within these four major clades. A signal of population expansion was detected in most clades or subclades. Coalescence within all of the major clades of P. misonnei occurred during the Middle Pleistocene and/or the beginning of Late Pleistocene. Main conclusions Our results suggest a role for both Pleistocene refugia and rivers in structuring genetic diversity in P. misonnei. This forest‐dwelling rodent may have been isolated in a number of forest fragments during arid periods and expanded its range during wetter periods. Potential forest refugia may have been localized in Benin–Ghana, south‐western Cameroon, southern Gabon, northern Gabon and eastern Democratic Republic of Congo–western Uganda. The Niger and/or the Cross Rivers, the Oubangui‐Congo, Sanaga, Ogooue and/or Ivindo Rivers probably stopped the re‐expansion of the species from relict areas.     

Staying out in the cold: glacial refugia and mitochondrial DNA phylogeography in ancient European brown bears

Models for the development of species distribution in Europe typically invoke restriction in three temperate Mediterranean refugia during glaciations, from where recolonization of central and northern Europe occurred. The brown bear, Ursus arctos, is one of the taxa from which this model is derived. Sequence data generated from brown bear fossils show a complex phylogeographical history for western European populations. Long-term isolation in separate refugia is not required to explain our data when considering the palaeontological distribution of brown bears. We propose continuous gene flow across southern Europe, from which brown bear populations expanded after the last glaciation.     

Chloroplast DNA supports a hypothesis of glacial refugia over postglacial recolonization in disjunct populations of black pine (Pinus nigra) in western Europe

European black pine (Pinus nigra Arn.) is a widely distributed Mediterranean conifer. To test the hypothesis that fragmented populations in western Europe survived in situ during the last glacial rather than having been re-colonized in the postglacial period, genetic variation was assessed using a suite of 10 chloroplast DNA microsatellites. Among 311 individuals analysed, 235 haplotypes were detected revealing high levels of chloroplast haplotype diversity in most populations. Bayesian analysis using a model of linked loci, with no prior assumption of population structure, assigned individuals to 10 clusters that corresponded well with the six predefined sampling regions, while an analysis carried out at the population level and assuming unlinked loci, recovered the original six sampling regions. This regional structure was supported by a biogeographical analysis that detected five barriers, with the two most significant separating Alps from Corsica and southern Italy, and southern Spain from the Pyrenees. No signals of demographic expansion were detected, and comparisons of R(ST) with pR(ST) suggested that a stepwise mutational model was important in regional differentiation, but not in population-within-region differentiation. These tests support long-term persistence of the species within the six regions. The temporal depth estimate, assuming a high mutation rate in coalescent modelling, placed the deepest split between the Alps and the other regions at about 150 000 years ago, and the most recent split of Pyrenees from southern France at about 30 000 years ago. Taken together, the data suggest that chloroplast DNA is structured in black pine and disjunct populations in western Europe are likely to have been present during the Last Glacial Maximum.     

Comparative phylogeography of codistributed aquatic insects in western North America: insights into dispersal and regional patterns of genetic structure

相似文献   

Glacial refugia and modern genetic diversity of 22 western North American tree species

North American tree species, subspecies and genetic varieties have primarily evolved in a landscape of extensive continental ice and restricted temperate climate environments. Here, we reconstruct the refugial history of western North American trees since the last glacial maximum using species distribution models, validated against 3571 palaeoecological records. We investigate how modern subspecies structure and genetic diversity corresponds to modelled glacial refugia, based on a meta-analysis of allelic richness and expected heterozygosity for 473 populations of 22 tree species. We find that species with strong genetic differentiation into subspecies had widespread and large glacial refugia, whereas species with restricted refugia show no differentiation among populations and little genetic diversity, despite being common over a wide range of environments today. In addition, a strong relationship between allelic richness and the size of modelled glacial refugia (r² = 0.55) suggest that population bottlenecks during glacial periods had a pronounced effect on the presence of rare alleles.     

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

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

Pleistocene glacial cycles drove lineage diversification and fusion in the Yosemite toad (Anaxyrus canorus)

Species endemic to alpine environments can evolve via steep ecological selection gradients between lowland and upland environments. Additionally, many alpine environments have faced repeated glacial episodes over the past two million years, fracturing these endemics into isolated populations. In this “glacial pulse” model of alpine diversification, cycles of allopatry and ecologically divergent glacial refugia play a role in generating biodiversity, including novel admixed (“fused”) lineages. We tested for patterns of glacial pulse lineage diversification in the Yosemite toad (Anaxyrus [Bufo] canorus), an alpine endemic tied to glacially influenced meadow environments. Using double‐digest RADseq on populations densely sampled from a portion of the species range, we identified nine distinct lineages with divergence times ranging from 18 to 724 thousand years ago (ka), coinciding with multiple Sierra Nevada glacial events. Three lineages have admixed origins, and demographic models suggest these fused lineages have persisted throughout past glacial cycles. Directionality indices supported the hypothesis that some lineages recolonized Yosemite from east of the ice sheet, whereas other lineages remained in western refugia. Finally, refugial niche reconstructions suggest that low‐ and high‐elevation lineages have convergently adapted to similar climatic niches. Our results suggest glacial cycles and refugia may be important crucibles of adaptive diversity across deep evolutionary time.     

A complex scenario of glacial survival in Mediterranean and continental refugia of a temperate continental vole species (Microtus arvalis) in Europe

The role of glacial refugia in shaping contemporary species distribution is a long-standing question in phylogeography and evolutionary ecology. Recent studies are questioning previous paradigms on glacial refugia and postglacial recolonization pathways in Europe, and more flexible phylogeographic scenarios have been proposed. We used the widespread common vole Microtus arvalis as a model to investigate the origin, locations of glacial refugia, and dispersal pathways, in the group of “Continental” species in Europe. We used a Bayesian spatiotemporal diffusion analysis (relaxed random walk model) of cytochrome b sequences across the species range, including newly collected individuals from 10 Iberian localities and published sequences from 68 localities across 22 European countries. Our data suggest that the species originated in Central Europe, and we revealed the location of multiple refugia (in both southern peninsulas and continental regions) for this continental model species. Our results confirm the monophyly of Iberian voles and the pre-LGM divergence between Iberian and European voles. We found evidence of restricted postglacial dispersal from refugia in Mediterranean peninsulas. We inferred a complex evolutionary and demographic history of M. arvalis in Europe over the last 50,000 years that does not adequately fit previous glacial refugial scenarios. The phylogeography of M. arvalis provides a paradigm of ice-age survival of a temperate continental species in western and eastern Mediterranean peninsulas (sources of endemism) and multiple continental regions (sources of postglacial spread). Our findings also provide support for a major role of large European river systems in shaping geographic boundaries of M. arvalis in Europe.     

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.     

Multiple genetic divergences and population expansions of a Mediterranean sandfly,Phlebotomus ariasi,in Europe during the Pleistocene glacial cycles

Phlebotomus ariasi is one of the two sandflies transmitting the causative agent of zoonotic leishmaniasis, Leishmania infantum, in France and Iberia, and provides a rare case study of the postglacial re-colonization of France by a Mediterranean species. Four DNA sequences were analysed—mitochondrial cytochrome b (cyt b), nuclear elongation factor-1α (EF-1α) and two anonymous nuclear loci—for 14–15 French populations and single populations from northeast Spain, northwest Spain, Portugal and Morocco. The presence of cryptic sibling species was not revealed by phylogenetic analyses and testing for reproductive isolation between sympatric populations defined by the two most divergent cyt b haplogroups. No locus was shown to be under positive directional or balancing selection and, therefore, molecular variation was explained demographically. Each nuclear locus showed shallow isolation by distance from Portugal to the French Pyrenees, but for both cyt b and EF-1α there was then a step change to the upland Massif Central, where leading-edge populations showed low diversity at all loci. Multiple genetic divergences and population expansions were detected by analyses of cyt b and dated to the Pleistocene. Endemicity of one cyt b sub-lineage suggested the presence of a refuge north of the Pyrenees during the last glacial period. Monopolization of the Massif Central by genetically differentiated populations of P. ariasi might possibly hinder the northwards spread of leishmaniasis.     

Multiple Pleistocene refugia and post‐glacial colonization in the European chub (Squalius cephalus) revealed by combined use of nuclear and mitochondrial markers

Aim To analyse patterns of nuclear and mitochondrial genetic variation in the European chub, Squalius cephalus (Linnaeus, 1758), in order to understand the evolutionary history of this species and to test biogeographical hypotheses for the existence of co‐distributed European freshwater fish species. Location Rivers in Europe (Finland, Poland, Czech Republic, France, Bulgaria, Spain, Italy). Methods We genotyped 12 polymorphic microsatellite markers derived from 310 individuals collected from across the distribution of S. cephalus in Europe (including a total of 15 populations) and sequenced mitochondrial DNA (mtDNA) from a subset of 75 individuals. Sequences of mtDNA cytochrome b were analysed using both phylogenetic (median‐joining networks) and population genetic methods (tests for demographic history, mismatch distributions, Bayesian coalescent analysis). Geographical structure in microsatellite loci was examined using a distance method (F_ST), factorial correspondence analysis (FCA) and a Bayesian clustering method (structure ). Results The mtDNA network showed a clear split into four different haplogroup lineages: Western (separated into Atlantic and Danubian sublineages), Eastern, Aegean (occurring in two distinct sublineages in the Balkans and in Spain) and Adriatic. Our results indicate recent population expansion in the Eastern and Western Atlantic lineages and the admixture of two previously separate sublineages (Atlantic and Danubian) in the Western lineage. Bayesian structure analysis as well as FCA results roughly corresponded to the mtDNA‐based structure, separating the sampled individuals into almost non‐overlapping groups. Main conclusions Our results support hypotheses suggesting origins of extant lineages of freshwater fishes in multiple refugia and the subsequent post‐glacial colonization of Europe via different routes. We confirmed the previously proposed two‐step expansion scenario from the Danube refuge, the existence of a secondary (Atlantic) refuge during the last glaciation (probably in the Rhone River) and population expansion of this lineage. Conspicuous divergences among Mediterranean populations reflect their different origin, as well as their low contribution to the recent genetic pool of chub in central Europe.     

Circum Black Sea phylogeography of Barbus freshwater fishes: divergence in the Pontic glacial refugium

We used DNA sequence variation at the mitochondrial cytochrome b gene (1141 bp) to assess the phylogeography of Barbus fishes in the Black Sea region. Our aim was to test whether the recent ( approximately 22 000-7500 years ago) freshwater phase of the Black Sea was a conduit for gene flow among freshwater fishes that today are found in streams entering the saltwater Black Sea. Deep phylogeographical breaks suggestive of allopatric divergence were observed between four regional groups of populations. Coalescent simulation used to distinguish between this and an alternative scenario that the phylogeographical structure was due to random lineage sorting showed that the contemporary populations were unlikely (P < 0.001) to have been founded by a single ancestral population. Divergences between the lineages (0.86-2.54%) were dated to the Middle to Late Pleistocene using distances and a molecular clock corrected for superimposed substitutions. Taken together, this evidence suggests that multiple refugial populations survived over several later glaciations in the vicinity of the Black Sea. This Pontic refugium served as the primary source for the postglacial expansion throughout Europe as far as the Atlantic basin. However, only one of the phylogeographical lineages contributed to this dispersion, whereas the others remain restricted to the Black Sea region and followed independent evolutionary trajectories.     

On glacial refugia,genetic diversity,and microevolutionary processes: deep phylogeographical structure in the endemic newt Lissotriton italicus1

We investigated the phylogeographical patterns of Lissotriton italcus, a newt endemic to the Italian peninsula, aiming to determine why hotspots of intraspecific diversity so ‘hot’. We found two main mitochindrial DNA lineages (net sequence divergence of 6.8% at two fragments of total length of 1897 bp): one restricted to part of the Calabrian peninsula (i.e the southernmost portion of the species range) and the other widespread throughout the rest of the species range. Both lineages, which had a parapatric distribution, showed evidence of further subdivisions, with an overall number of eight terminal haplogroups, most of whose times to the most recent common ancestors were estimated at the Late Pleistocene. Analysis of molecular variance suggested that partitioning populations according to the geographical distribution of these haplogroups can explain 97% of the observed genetic variation. These results suggest that L. italicus underwent repeated cycles of allopatric fragmentation throughout the Pleistocene, and that it likely survived the Late Pleistocene paleoenvironmental changes within eight separate refugia. Thus, the current hotspot of intraspecific diversity of L. italicus (within the Calabrian peninsula) has not been moulded by long‐term stability of large populations but rather by multiple events of allopatric fragmentation and divergence. When compared with the patterns recently identified in other species, these results suggest that the occurrence of phases of allopatric divergence (eventually followed by secondary admixture) could be a common, albeit probably underrated feature in the history of formation of hotspots of intraspecific diversity. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 105 , 42–55.