European colonization by the spined loach (Cobitis taenia) from Ponto-Caspian refugia based on mitochondrial DNA variation

In the last 20 years, new species, asexual reproduction, polyploidy and hybridization have all been reported within the genus Cobitis. An understanding of the current distribution and baseline phylogeographical history of 'true' nonhybrid Cobitis species is crucial in order to unravel these discoveries. In the present work, we investigated the phylogeography of the spined loach, Cobitis taenia, using 1126 bp of the mitochondrial cytochrome b gene from 174 individuals collected at 47 sites. In total, 51 haplotypes that differed at 49 positions (4.35%) were detected. We deduce that C. taenia survived European glaciations in at least three refugees in the Ponto-Caspian area. Two of these refugees each provided a major lineage that recolonized Europe in separate directions: one westward to England and the other spreading north into Russia before moving west. A third (minor) lineage that contributed little to the recolonization of Europe was also revealed--remaining near its Black Sea refuge. However, more recent history was difficult to resolve with colonization from a more western refugium during the last glacial maximum (LGM) a distinct possibility. Nested clade analysis indicates a pattern of restricted gene flow with isolation by distance at the first two levels and overall. Unlike many other European freshwater fish species, the Danube is not part of the current distribution of C. taenia, nor was it used as either a refuge or a source of colonization of Europe. Low genetic diversity within C. taenia suggests that its colonization of Europe is relatively recent. Demographic analyses revealed a history of recent expansion and isolation by distance.     

Multiple ice‐age refugia in Pacific cod,Gadus macrocephalus

Pleistocene ice‐ages greatly influenced the historical abundances of Pacific cod, Gadus macrocephalus, in the North Pacific and its marginal seas. We surveyed genetic variation at 11 microsatellite loci and mitochondrial (mt) DNA in samples from twelve locations from the Sea of Japan to Washington State. Both microsatellite (mean H = 0.868) and mtDNA haplotype (mean h = 0.958) diversities were large and did not show any geographical trends. Genetic differentiation between samples was significantly correlated with geographical distance between samples for both microsatellites (F_ST = 0.028, r² = 0.33) and mtDNA (F_ST = 0.027, r² = 0.18). Both marker classes showed a strong genetic discontinuity between northwestern and northeastern Pacific populations that likely represents groups previously isolated during glaciations that are now in secondary contact. Significant differences appeared between samples from the Sea of Japan and Okhotsk Sea that may reflect ice‐age isolations in the northwest Pacific. In the northeast Pacific, a microsatellite and mtDNA partition was detected between coastal and Georgia Basin populations. The presence of two major coastal mtDNA lineages on either side of the Pacific Ocean basin implies at least two ice‐age refugia and separate postglacial population expansions facilitated by different glacial histories. Northward expansions into the Gulf of Alaska were possible 14–15 kyr ago, but deglaciation and colonization of the Georgia Basin probably occurred somewhat later. Population expansions were evident in mtDNA mismatch distributions and in Bayesian skyline plots of the three major lineages, but the start of expansions appeared to pre‐date the last glacial maximum.     

Pleistocene extinctions as drivers of biogeographical patterns on the easternmost Canary Islands

Subtropical islands are often viewed as refuges where Quaternary climatic shifts driving global episodes of extinction were buffered. Island biodiversity, however, may have been impacted by climatic fluctuations at local scales, particularly in spatially heterogeneous island systems. In this study, we generated a conceptual framework for predicting the potential impact of Pleistocene extinctions on the biogeographical pattern of the Canarian spermatophyte flora, with a focus on the easternmost Canarian islands (ECI). Then, we performed an exhaustive bibliographic revision (270 studies) to examine whether taxonomic, phylogenetic and phylogeographical data support our predictions. Although molecular information is limited for many lineages, the available data suggest that the majority of extant ECI plant taxa may be the result of relatively recent (<1 Ma) dispersal from surrounding insular and mainland areas. Different lines of evidence are compatible with the idea of a Pleistocene period of frequent lineage extirpation on ECI. Extinction may thus have provided new ecological opportunities for recent (re)colonization, with some cases of recent establishment mediated by facilitation. Considering background extinction on ECI, we describe five general patterns of colonization for Canarian plant lineages. In addition to factors related to island ontogeny and long‐distance dispersal, we suggest that Pleistocene extinctions may have significantly contributed to extant biogeographical patterns in the Canarian archipelago, such as the biased distribution ranges of island plants and the low endemic richness on ECI. This new scenario provides testable hypotheses for future studies dealing with the phylogeography, taxonomy and conservation of terrestrial biodiversity on the Canarian islands, and possibly, on other near‐shore islands.     

Truffles in the sky: The impact of stochastic and deterministic drivers on Rhizopogon communities of the Madrean sky island archipelago

Temperate forests across the globe are migrating as a result of global warming, but little is known about how changes in climate and host geographic distributions will affect the obligate symbionts of these forests. Analysis of past events can provide insight into how these symbionts have responded to previous climatic changes and inform predictions for contemporary and future climate change events. The Madrean Sky Islands Archipelago (MSIA) comprises mountain “islands” whose Pine-Oak forests appear in stark contrast to the surrounding “sea” of Sonoran Desert vegetation. The archipelago formed due to a post-Pleistocene warming climate that resulted in the expansion of the Sonoran Desert and migration of forests to fragmented, higher elevation areas. Rhizopogon (Boletales) consists of obligate ectomycorrhizal (EcM) symbionts that form truffle sporocarps and associate exclusively with Pinaceae. As such, the MSIA-Rhizopogon system represents a natural experiment of how fungal symbionts responded to climatic change and host migration. Rhizopogon was sampled from nine islands at two sites (one Pinus site and one Pseudotsuga site) per island, and diversity was characterized using the ITS rRNA gene determined from both sporocarps and bioassay-based EcM root tips derived from soil samples collected at each site. We described the biodiversity of Rhizopogon within and among sky islands of the MSIA, and tested whether symbiont species richness and community structure were determined by host association, island identity, geographic distance, or some interaction among these factors. Twenty-five OTUs at 99% similarity in the genus Rhizopogon were identified across nine sky islands with a range of 5–15 OTUs per island. While differential host association with Pinus and Pseudotsuga was a significant driver of community composition, our results supported an even stronger island effect. Furthermore, Rhizopogon communities associated with Pinus forest sites were characterized by random phylogenetic structures across sky islands and are not structured by geographic distance. Our results supported a strong isolation effect that involved historical habitat fragmentation of sky islands in response to past climate changes, and that both host association and stochastic processes, e.g., ecological drift, played a role in shaping Rhizopogon communities of the MSIA.     

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.     

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 late Quaternary phylogeography for Metrosideros (Myrtaceae) in New Zealand inferred from chloroplast DNA haplotypes

Because of its remoteness and mid-latitude position, New Zealand lacks access to the tropical climates that might have ensured straightforward survival for frost-sensitive species during glacial times. Nevertheless, the New Zealand lowland flora retains a substantial complement of plants sourced in the tropics. While there have been extinction pulses for elements of the frost-sensitive flora under glacial/stadial regimes, the surviving remnants have been able to recolonize large areas of habitat during successive warm climate periods. Refugia for such species in stadial New Zealand are likely to have been localized and ecologically suboptimal. To examine these relationships we have applied chloroplast DNA sequence data to the investigation of phylogeographical pattern for five endemic species of Metrosideros subg. Metrosideros , a wide-ranging group of mostly frost-sensitive woody plants in New Zealand. The results of this research verify the location of two generally mooted stadial refugia for the country and provide support for the existence of a third. A simple pattern of chloroplast haplotype diversity was recorded in extra-refugial areas, compared with a greater complexity in the vicinity of the identified refugia. This pattern was independently repeated in both main islands. The proposed refugia correspond to contemporary localities of high average winter temperatures. The sharing of chloroplast haplotypes between the different species of Metrosideros examined suggests that there has been a history of repeated hybridization and introgression for these plants, possibly initiated by periods of refugial confinement.  © 2004 The Linnean Society of London, Biological Journal of the Linnean Society , 2004, 83 , 399–412.     

Caucasus Mountains divide postulated postglacial colonization routes in the white‐breasted hedgehog,Erinaceus concolor

For many European species, the mountains of the Alps and the Pyrenees have acted as significant barriers to northwards colonization from southern glacial refugia. To the east, the Caucasus Mountains would seem to have been a similar barrier to the white‐breasted hedgehog (Erinaceus concolor). A deep divergence among hedgehog mitochondrial sequences to the north and south of the Caucasus Mountains suggests two colonization routes, originating from separate refugial regions and divided by this mountain barrier. From a Balkan refugium, hedgehogs have colonized northwards into Russia and to the northern foothills of the Caucasus Mountains. The origins of hedgehogs colonizing the southern parts of the Caucasus are not entirely clear, although fossil and climatic data suggest a glacial refugium on the southern shores of the Black Sea. Divergence within the southern group indicates a long‐standing fragmentation within such a refugium or the presence of further cryptic refugia in Turkey and the Near East. The Caucasus barrier would seem to have been an important factor in structuring the late Pleistocene distribution of species.     

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

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

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.     

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

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

Land masses and oceanic currents drive population structure of Heritiera littoralis,a widespread mangrove in the Indo‐West Pacific

Phylogeographic forces driving evolution of sea‐dispersed plants are often influenced by regional and species characteristics, although not yet deciphered at a large spatial scale for many taxa like the mangrove species Heritiera littoralis. This study aimed to assess geographic distribution of genetic variation of this widespread mangrove in the Indo‐West Pacific region and identify the phylogeographic factors influencing its present‐day distribution. Analysis of five chloroplast DNA fragments’ sequences from 37 populations revealed low genetic diversity at the population level and strong genetic structure of H. littoralis in this region. The estimated divergence times between the major genetic lineages indicated that glacial level changes during the Pleistocene epoch induced strong genetic differentiation across the Indian and Pacific Oceans. In comparison to the strong genetic break imposed by the Sunda Shelf toward splitting the lineages of the Indian and Pacific Oceans, the genetic differentiation between Indo‐Malesia and Australasia was not so prominent. Long‐distance dispersal ability of H. littoralis propagules helped the species to attain transoceanic distribution not only across South East Asia and Australia, but also across the Indian Ocean to East Africa. However, oceanic circulation pattern in the South China Sea was found to act as a barrier creating further intraoceanic genetic differentiation. Overall, phylogeographic analysis in this study revealed that glacial vicariance had profound influence on population differentiation in H. littoralis and caused low genetic diversity except for the refugia populations near the equator which might have persisted through glacial maxima. With increasing loss of suitable habitats due to anthropogenic activities, these findings therefore emphasize the urgent need for conservation actions for all populations throughout the distribution range of H. littoralis.     

Native range genetic variation in Arabidopsis thaliana is strongly geographically structured and reflects Pleistocene glacial dynamics

Despite Arabidopsis thaliana 's pre-eminence as a model organism, major questions remain regarding the geographic structure of its genetic variation due to the geographically incomplete sample set available for previous studies. Many of these questions are addressed here with an analysis of genome-wide variation at 10 loci in 475 individuals from 167 globally distributed populations, including many from critical but previously un-sampled regions. Rooted haplotype networks at three loci suggest that A. thaliana arose in the Caucasus region. Identification of large-scale metapopulations indicates clear east–west genetic structure, both within proposed Pleistocene refugia and post-Pleistocene colonized regions. The refugia themselves are genetically differentiated from one another and display elevated levels of within-population genetic diversity relative to recolonized areas. The timing of an inferred demographic expansion coincides with the Eemian interglacial (approximately 120 000 years ago). Taken together, these patterns are strongly suggestive of Pleistocene range dynamics. Spatial autocorrelation analyses indicate that isolation by distance is pervasive at all hierarchical levels, but that it is reduced in portions of Europe.     

Mitochondrial haplotype diversity among Portuguese brown trout Salmo trutta L. populations: relevance to the post-pleistocene recolonization of northern Europe

Mitochondrial haplotype diversity in seven Portuguese populations of brown trout, Salmo trutta L., was investigated by sequencing the 5' end of the mitochondrial DNA (mtDNA) control region. Five new haplotypes were described for this species, each two to three mutational steps distant from the common north Atlantic haplotype. Significant population subdivision of mtDNA haplotypes was also apparent. Based on these results, as well as on published data describing the distribution of both mtDNA haplotypes and allozyme alleles throughout Europe, the postglacial recolonization of northern Europe was re-evaluated. It is argued that the available data do not support the contribution of two major glacial refugia (southwest Atlantic and Ponto-Caspian Basin) to this postglacial recolonization, as proposed in a recently published model. The unique genetic architecture of Portuguese brown trout within the Atlantic-basin clade of this species represents a highly valuable genetic resource that should be protected from introgression with nonendemic strains of hatchery fish.     

The influence of Pleistocene glacial refugia on tawny owl genetic diversity and phylogeography in western Europe

The glacial refugia hypothesis indicates that during the height of the Pleistocene glaciations the temperate species that are today widespread in western Europe must have survived in small and climatically favourable areas located in the southern peninsulas of Iberia, Italy and Balkans. One such species is the tawny owl, a relatively sedentary, nonmigratory bird presently distributed throughout Europe. It is a tree-nesting species closely associated with deciduous and mixed coniferous woodlands. In this study I used control region mtDNA sequences from 187 individuals distributed among 14 populations to determine whether current genetic patterns in tawny owl populations were consistent with postglacial expansion from peninsular refugia. European, North African and Asian tawny owls were found to represent three distinct lineages, where North Africa is the sister clade to all European owls. Within Europe, I found three well-supported clades that correspond to each of the three allopatric refugia. Expansion patterns indicate that owls from the Balkan refugium repopulated most of northern Europe, while expansion out of Iberia and Italy had only regional effects leading to admixture in France. Estimates of population divergence times between refugia populations are roughly similar, but one order of magnitude smaller between Greece and northern Europe. Based on a wide range of mutation rates and generation times, divergence between refugia appears to date to the Pleistocene.     

Persistent genetic signatures of colonization in Brachionus manjavacas rotifers in the Iberian Peninsula

Recent phylogeographical assessments have consistently shown that continental zooplankton display high levels of population subdivision, despite the high dispersal capacity of their diapausing propagules. As such, there is an apparent paradox between observed cosmopolitanism in the zooplankton that is associated with long-distance dispersal, and strong phylogeographical structures at a regional scale. Such population dynamics, far from migration-drift equilibrium, have been shown in the rotifer species complex Brachionus plicatilis, a group of over a dozen species inhabiting salt lakes and coastal lagoons worldwide. Here we present the mitochondrial DNA phylogeography of one of these species, Brachionus manjavacas, in the Iberian Peninsula, where it often co-occurs with the morphologically similar species B. plicatilis sensu stricto. We obtained sequences from 233 individuals from diapausing eggs and clonal cultures from 16 lakes in the Iberian Peninsula, and a Tunisian lake. Two strongly supported deep mitochondrial DNA clades were found (A and B). Phylogenetic and nested clade analysis showed that clade A has a strong phylogeographical structure, with a strong similarity of phylogeographical patterns between B. manjavacas clade A and B. plicatilis s.s. These include (i) signatures of allopatric fragmentation between central and southern populations, and (ii) range expansions in the Iberian Peninsula, both likely to have occurred during the Pleistocene. We find evidence for a glacial refugium in the Guadiana basin. Clades A and B co-occurred in several of these lakes because of range expansion and secondary contact between both clades. The co-occurrence between B. plicatilis s.s. and B. manjavacas is not recent, and both species might have experienced similar environmental challenges during the Pleistocene. The strong correlation of genetic and geographical distance found suggests that historical events can lead to such correlation, mirroring the effects of 'isolation by distance' in equilibrium populations.     

Phylogeography of the burnet moth Zygaena transalpina complex: molecular and morphometric differentiation suggests glacial refugia in Southern France,Western France and micro‐refugia within the Alps

Patterns of common recolonization routes from glacial refugia to Central Europe during the Pleistocene are generalized to paradigms of postglacial recolonization in Europe. Recent studies indicate, however, that the actual phylogeographic history of many species might be more complex and cannot be simplified to generalized patterns. Burnet moths of the Zygaena transalpina complex represent a group of closely related taxa, which are considered as a typical example for these generalized patterns. At present, three groups are recognized that are assumed to have spread from three classical refugia in Western Europe, Italy and the Balkans to Central Europe. Here, we re‐investigate their phylogeography using a combined molecular and morphometric approach. Phylogenetic and nested clade phylogeographic analyses of 476 samples from 55 localities taken from Southern and Central Europe reveal that the Zygaena transalpina complex consists of three distinct haplotype clusters, which geographically roughly correspond to possible refugia in Western Europe, Italy and the Balkans. A synthesis of the data with a geometric morphometry dataset of 425 specimens from 46 localities corroborates this molecular result but differs in several aspects. Important new aspects are multiple refugia of the western ‘hippocrepidis’ branch and micro‐habitats within the Alps of the central ‘transalpina’ branch. Further, our results display a more complex phylogeographic pattern for this species complex, which is not tractable with a rigid, generalized pattern.