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.     

Phylogeography of the marine isopod Stenosoma nadejda (Rezig, 1989) in North African Atlantic and western Mediterranean coasts reveals complex differentiation patterns and a new species

The transition zone between the Mediterranean and Atlantic basins has been extensively addressed in phylogeographical studies of marine species. However, biases exist towards the analysis of highly dispersive species, and there is a higher sampling effort in European coasts compared to North Africa. This may be hindering a detailed understanding of the historical and contemporary processes that shaped patterns of population genetic structure in the region. In the present study, we investigated the phylogeographical and phylogenetic patterns of mitochondrial cytochrome c oxidase subunit I sequences from a species with direct development and low dispersal abilities, Stenosoma nadejda (Rezig, 1989). The study area included 13 localities along the Atlantic and Mediterranean North African coasts, as well as the Alboran Sea. A new Stenosoma species, from the coasts of Algeria and Alboran Island, was discovered. For S. nadejda, phylogeographical analyses revealed three distinct clades: one in the Iberian Atlantic plus the Alboran Sea, one in the western Mediterranean, and another in the Atlantic coast of Africa. Haplotypes from the Alboran Island were more related to those from the western Mediterranean coast (east of the Almeria–Oran Front). Given the strong differentiation, it is probable that this species survived in multiple glacial refugia during the Pleistocenic glaciations. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 104 , 419–431.     

Phylogeography and historical demography of the Lusitanian snail <Emphasis Type="Italic">Elona quimperiana</Emphasis> reveal survival in unexpected separate glacial refugia

Background  

Present day distributions of Palearctic taxa in northern latitudes mainly result from populations having survived in local patches during the Late Pleistocene and/or from recolonizing populations from southern temperate refugia. If well-studied Mediterranean and eastern European refugia are widely accepted, some recent biogeographical assumptions still remain unclear, such as the occurrence of multiple glacial refugia in Iberia and cryptic refugia in northern Europe during the last glaciations. The Lusitanian snail Elona quimperiana has a remarkably disjunct distribution, limited to northwestern France (Brittany), northwestern Spain and the Basque Country. By describing the phylogeographical structure of this species across its entire range, the present study attempts to identify refugia and subsequent recolonization routes.     

Glacial history of the North Atlantic marine snail, Littorina saxatilis, inferred from distribution of mitochondrial DNA lineages

The North Atlantic intertidal gastropod, Littorina saxatilis (Olivi, 1792), exhibits extreme morphological variation between and within geographic regions and has become a model for studies of local adaptation; yet a comprehensive analysis of the species' phylogeography is lacking. Here, we examine phylogeographic patterns of the species' populations in the North Atlantic and one remote Mediterranean population using sequence variation in a fragment of the mitochondrial cytochrome b gene (607 bp). We found that, as opposed to many other rocky intertidal species, L. saxatilis has likely had a long and continuous history in the Northwest Atlantic, including survival during the last glacial maximum (LGM), possibly in two refugia. In the Northeast Atlantic, several areas likely harboured refugial populations that recolonized different parts of this region after glacial retreat, resulting in strong population structure. However, the outlying monomorphic Venetian population is likely a recent anthropogenic introduction from northern Europe and not a remnant of an earlier wider distribution in the Mediterranean Sea. Overall, our detailed phylogeography of L. saxatilis adds an important piece to the understanding of Pleistocene history in North Atlantic marine biota as well as being the first study to describe the species' evolutionary history in its natural range. The latter contribution is noteworthy because the snail has recently become an important model species for understanding evolutionary processes of speciation; thus our work provides integral information for such endeavours.     

Identifying refugia from climate change using coupled ecological and genetic data in a transitional Mediterranean‐temperate tree species

Populations occurring in areas of overlap between the current and future distribution of a species are particularly important because they can represent “refugia from climate change”. We coupled ecological and range‐wide genetic variation data to detect such areas and to evaluate the impacts of habitat suitability changes on the genetic diversity of the transitional Mediterranean‐temperate tree Fraxinus angustifolia. We sampled and genotyped 38 natural populations comprising 1006 individuals from across Europe. We found the highest genetic diversity in western and northern Mediterranean populations, as well as a significant west to east decline in genetic diversity. Areas of potential refugia that correspond to approximately 70% of the suitable habitat may support the persistence of more than 90% of the total number of alleles in the future. Moreover, based on correlations between Bayesian genetic assignment and climate, climate change may favour the westward spread of the Black Sea gene pool in the long term. Overall, our results suggest that the northerly core areas of the current distribution contain the most important part of the genetic variation for this species and may serve as in situ macrorefugia from ongoing climate change. However, rear‐edge populations of the southern Mediterranean may be exposed to a potential loss of unique genetic diversity owing to habitat suitability changes unless populations can persist in microrefugia that have facilitated such persistence in the past.     

Strong population genetic structure and contrasting demographic histories for the small-spotted catshark (Scyliorhinus canicula) in the Mediterranean Sea

Coastal and demersal chondrichthyans, such as the small-spotted catshark, are expected to exhibit genetic differentiation in areas of complex geomorphology like the Mediterranean Basin because of their limited dispersal ability. To test this hypothesis, we used a fragment of the mitochondrial cytochrome c oxidase subunit I gene and 12 nuclear microsatellite loci in order to investigate the genetic structure and historical demography of this species, and to identify potential barriers to gene flow. Samples were collected from the Balearic Islands, the Algerian Basin, the Ionian Sea, the Corinthian Gulf and various locations across the Aegean Sea. Additional sequences from the Atlantic and the Levantine Basin retrieved from GenBank were included in the mitochondrial DNA analysis. Both mitochondrial and nuclear microsatellite DNA data revealed a strong genetic subdivision, mainly between the western and eastern Mediterranean, whereas the Levantine Basin shared haplotypes with both areas. The geographic isolation of the Mediterranean basins seems to enforce the population genetic differentiation of the species, with the deep sea acting as a strong barrier to its dispersal. Contrasting historical demographic patterns were also observed in different parts of the species'' distribution, most notably a population growth trend in the western Mediterranean/Atlantic area and a slight decreasing one in the Aegean Sea. The different effects of the Pleistocene glacial periods on the habitat availability may explain the contrasting demographic patterns observed. The current findings suggest that the small-spotted catshark exhibits several genetic stocks in the Mediterranean, although further study is needed.     

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.     

Is Niagara Falls a barrier to gene flow in riverine fishes? A test using genome‐wide SNP data from seven native species

Since the early Holocene, fish population genetics in the Laurentian Great Lakes have been shaped by the dual influences of habitat structure and post‐glacial dispersal. Riverscape genetics theory predicts that longitudinal habitat corridors and unidirectional downstream water‐flow drive the downstream accumulation of genetic diversity, whereas post‐glacial dispersal theory predicts that fish genetic diversity should decrease with increasing distance from glacial refugia. This study examines populations of seven native fish species codistributed above and below the 58 m high Niagara Falls – a hypothesized barrier to gene flow in aquatic species. A better understanding of Niagara Falls’ role as a barrier to gene flow and dispersal is needed to identify drivers of Great Lakes genetic diversity and guide strategies to limit exotic species invasions. We used genome‐wide SNPs and coalescent models to test whether populations are: (a) genetically distinct, consistent with the Niagara Falls barrier hypothesis; (b) more genetically diverse upstream, consistent with post‐glacial expansion theory, or downstream, consistent with the riverscape habitat theory; and (c) have migrated either upstream or downstream past Niagara Falls. We found that genetic diversity is consistently greater below Niagara Falls and the falls are an effective barrier to migration, but two species have probably dispersed upstream past the falls after glacial retreat yet before opening of the Welland Canal. Models restricting migration to after opening of the Welland Canal were generally rejected. These results help explain how river habitat features affect aquatic species’ genetic diversity and highlight the need to better understand post‐glacial dispersal pathways.     

A dynamic history of admixture from Mediterranean and Carpathian glacial refugia drives genomic diversity in the bank vole

Understanding the historical contributions of differing glacial refugia is key to evaluating the roles of microevolutionary forces, such as isolation, introgression, and selection in shaping genomic diversity in present‐day populations. In Europe, where both Mediterranean and extra‐Mediterranean (e.g., Carpathian) refugia of the bank vole (Clethrionomys glareolus) have been identified, mtDNA indicates that extra‐Mediterranean refugia were the main source of colonization across the species range, while Mediterranean peninsulas harbor isolated, endemic lineages. Here, we critically evaluate this hypothesis using previously generated genomic data (>6,000 SNPs) for over 800 voles, focusing on genomic contributions to bank voles in central Europe, a key geographic area in considering range‐wide colonization. The results provide clear evidence that both extra‐Mediterranean (Carpathian) and Mediterranean (Spanish, Calabrian, and Balkan) refugia contributed to the ancestry and genomic diversity of bank vole populations across Europe. Few strong barriers to dispersal and frequent admixture events in central Europe have led to a prominent mid‐latitude peak in genomic diversity. Although the genomic contribution of the centrally located Carpathian refugium predominates, populations in different parts of Europe have admixed origins from Mediterranean (28%–47%) and the Carpathian (53%–72%) sources. We suggest that the admixture from Mediterranean refugia may have provisioned adaptive southern alleles to more northern populations, facilitating the end‐glacial spread of the admixed populations and contributing to increased bank vole diversity in central Europe. This study adds critical details to the complex end‐glacial colonization history of this well‐studied organism and underscores the importance of genomic data in phylogeographic interpretation.     

Unravelling the evolutionary history and future prospects of endemic species restricted to former glacial refugia

The contemporary distribution and genetic composition of biodiversity bear a signature of species’ evolutionary histories and the effects of past climatic oscillations. For many European species, the Mediterranean peninsulas of Iberia, Italy and the Balkans acted as glacial refugia and the source of range recolonization, and as a result, they contain disproportionately high levels of diversity. As these areas are particularly threatened by future climate change, it is important to understand how past climatic changes affected their biodiversity. We use an integrated approach, combining markers with different evolutionary rates and combining phylogenetic analysis with approximate Bayesian computation and species distribution modelling across temporal scales. We relate phylogeographic processes to patterns of genetic variation in Myotis escalerai, a bat species endemic to the Iberian Peninsula. We found a distinct population structure at the mitochondrial level with a strong geographic signature, indicating lineage divergence into separate glacial refugia within the Iberian refugium. However, microsatellite markers suggest higher levels of gene flow resulting in more limited structure at recent time frames. The evolutionary history of M. escalerai was shaped by the effects of climatic oscillations and changes in forest cover and composition, while its future is threatened by climatically induced range contractions and the role of ecological barriers due to competition interactions in restricting its distribution. This study warns that Mediterranean peninsulas, which provided refuge for European biodiversity during past glaciation events, may become a trap for limited dispersal and ecologically limited endemic species under future climate change, resulting in loss of entire lineages.     

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.     

Ecological patterns strongly impact the biogeography of western Palaearctic longhorn beetles (Coleoptera: Cerambycoidea)

We aim to unravel the biogeographic structuring of western Palaearctic longhorn beetles with focus on the location of different refugia, barriers to dispersal and postglacial range expansions with their particular filters. The interaction of different ecological features with these structures is analysed. The western Palaearctic was divided into 95 geographic entities. We produced presence-only matrices for all 955 Cerambycoidea species autochthonous to this area and derived species richness distributions and extracted faunal regions and faunal elements by cluster analyses and principal component analyses. Similar analyses were performed for sub-families and ecological groups. Longhorn beetles show a strong biogeographic structuring in the western Palaearctic. Species numbers strongly decrease to the north and west. Less mobile species and root feeders mostly contribute to the fauna of the Mediterranean region, whilst mobile species are more widespread. Feeders on broad-leaved trees dominate in western Europe, whilst feeders on coniferous trees are most important in northern Europe. Our results support multiple refugia in the Mediterranean region and underline the importance of Provence, Crimea and Crete as such refugia. Crete even might be an area of old endemism. The Atlanto- and the Ponto-Mediterranean regions are more strongly structured than assumed in classical biogeography. Mediterranean assemblages are mostly composed of non-flying species, root feeders and species with small distributions not found outside their glacial refugia. Tree feeders left their glacial retreats with their host plants. These range dynamics result in biogeographic structures with several dispersal barriers and filters composed of mountains, sea straits and climatic conditions.     

Phylogeographic insights into <Emphasis Type="Italic">Artemisia crithmifolia</Emphasis> (Asteraceae) reveal several areas of the Iberian Atlantic coast as refugia for genetic diversity

The study of fragmentation processes should be tackled using different approaches, in order to obtain solid and robust evidence that could help in identifying potential barriers and threats for species. In this study, we have evaluated the spatial fragmentation patterns in Artemisia crithmifolia (Asteraceae) along its current distribution along the Atlantic coastlines of Europe, from Portugal to the Netherlands and the UK. Niche modelling analyses considering current and past climatic conditions, combined with plastid markers, have been used to evaluate the disconnected distribution pattern of the species, genetic flow between populations and habitat suitability for present and past populations. Plastid markers showed very low variability, while suggesting that the westernmost populations from the Atlantic coast of the Iberian Peninsula are partially genetically isolated, due to the existence of a potential barrier for gene flow on the northern Spanish coast. However, niche models identify this area as a potential habitat for the species, even during the last glacial maximum (20,000 years ago). The combination of both techniques allows identification of potential refugia for the species, highlighting the most likely recolonisation routes and distribution patterns which resulted in the overall low levels of genetic diversity. Anthropogenic activities (urban sprawl, beach expansions, etc.) are most likely behind local extinctions, thus preventing the establishment and expansion of new populations.     

Northern areas as refugia for temperate species under current climate warming: Atlantic salmon (Salmo salar L.) as a model in northern Europe

In this work, patterns of geographical genetic diversity in Atlantic salmon Salmo salar were studied across the whole Atlantic Arc; whether these patterns (and thus genetic population structure) were affected by water temperatures was also evaluated. Salmo salar populations were characterized using microsatellite loci and then analysed with reference to ocean surface temperature data from across the region. Analysis showed the presence of a latitudinal cline of genetic variability (higher in northern areas) and water temperatures (sea surface temperatures) determining genetic population structure (the latter in combination with genetic drift in southern populations). Under the current global change scenario, northern areas of Europe would constitute refugia for diversity in the future. This is effectively the inverse of what appears to have happened in glacial refugia during the last glacial maximum. From this perspective, the still abundant and large northern populations S. salar should be considered as precious as the small almost relict southern ones and given appropriate protection. Careful management of the species, coordinated across countries and latitudes, is needed in order to avoid its extinction in Europe.     

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

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

Rapid colonization of new habitats in the Wadden Sea by the ovoviviparousLittorina saxatilis (Olivi)

The intertidal periwinkleLittorina saxatilis completely lacks larval dispersal and adult vagility is low. Although this suggests a low dispersal rate,L. saxatilis is frequently found in recently established habitats “exotic” to the Wadden Sea. Populations occur on man-made structures like dikes, breakwater and groynes, some of which are not older than several years. Furthermore,L. saxatilis is found on marsh grassSpartina anglica, introduced to the Wadden Sea in the 1920s and 1930s, as well as on mats of green macroalgae, which have become an abundant feature on the tidal flats since the late 1970s. Seagrass beds are likely to be the original habitat ofL. saxatilis in the Wadden Sea. Since seagrass populations have dramatically declined over the last decades, colonization of new habitat types enabledL. saxatilis to maintain its Wadden Sea populations despite a changing environment. Colonizers can reach new habitats by means of passive transport, especially by rafting on macrophytes and by aerial dispersal attached to birds. In thew Wadden Sea, the ovoviviparously reproducingL. saxatilis has demonstrated its ability to successfully found new populations with only a few individuals. No reduction of genetic variablility (founder effect) was observed in recently established populations.     

Divergence and molecular variation in common whelk Buccinum undatum (Gastropoda: Buccinidae) in Iceland: a trans‐Atlantic comparison

The dispersal and history of species affects their genetic population structure at both small and large geographical scales. The common whelk, Buccinum undatum, is a widespread subtidal gastropod in the North Atlantic that has no planktonic larvae and has thus limited dispersal capacity. The snail, which has been harvested by humans for centuries, is highly variable in morphology. To evaluate the population structure in the rich fishing grounds in western Iceland and its divergence from samples across the Atlantic, genetic patterns based on sequence variation in two mitochondrial (mt)DNA genes (COI and 16S) and five microsatellites were studied and compared with variation in populations from both sides of the Atlantic. Significant differences in allele and haplotype frequencies were found among samples separated by short distances along the coast of Iceland. Partition of the variation showed larger variance among samples obtained from distant regions than from neighbouring sites and genetic distances were correlated with geographical distance among populations in Europe. Phylogeographic patterns in mtDNA reveal different monophyletic lineages on both sides of the Atlantic, which predate the onset of the Ice Age and which may constitute cryptic species. Similar micro‐ and macrogeographical patterns were observed for the mtDNA and microsatellite markers, despite high frequencies of null alleles. Bayesian skyline reconstructions of the demographic history and mismatch distributions suggest that, although sizes of some populations were unaffected by Ice Age glaciations, others show signs of expansion after the Last Glacial Maximum. These phylogeographical patterns are consistent with patterns expected for low dispersal species that have survived in allopatric glacial refugial populations on both sides of the Atlantic and in deep‐sea refugia within each continent. The observed genetic structure has implications for conservation and sustainable management of the harvested populations. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 145–159.     

Identifying refugia from climate change

This article highlights how the loose definition of the term ‘refugia’ has led to discrepancies in methods used to assess the vulnerability of species to the current trend of rising global temperatures. The term ‘refugia’ is commonly used without distinguishing between macrorefugia and microrefugia, ex situ refugia and in situ refugia, glacial and interglacial refugia or refugia based on habitat stability and refugia based on climatic stability. It is not always clear which definition is being used, and this makes it difficult to assess the appropriateness of the methods employed. For example, it is crucial to develop accurate fine‐scale climate grids when identifying microrefugia, but coarse‐scale macroclimate might be adequate for determining macrorefugia. Similarly, identifying in situ refugia might be more appropriate for species with poor dispersal ability but this may overestimate the extinction risk for good dispersers. More care needs to be taken to properly define the context when referring to refugia from climate change so that the validity of methods and the conservation significance of refugia can be assessed.     

Environmental confirmation of multiple ice age refugia for Pacific cod, Gadus macrocephalus

The concept of species surviving through quaternary climatic extremes by retreating to glacial refugia, and then evolving genetically during re-population movements of the following interglacial, has been in the literature for over 40 years. Recently, advances in genetic analysis have enabled this concept to be validated and theories regarding population expansions and contractions to be built. For the major Northern Hemisphere species of cod, Gadus morhua (Atlantic cod) and Gadus macrocephalus (Pacific cod), genetic analysis has suggested retreat to separate refugia on both sides of their respective ocean basins during the last glacial period. Ecosystem niche modelling has previously confirmed that environmental conditions during the last glacial were compatible with the existence of these separate refugia for Atlantic cod. Here it is shown that such modelling also confirms a reduced core glacial distribution for G. macrocephalus, but probable refugia on either side of the Pacific. Existing mitochondrial DNA analyses suggest two separate glacial populations in the northwest Pacific, which modelling confirms, with predicted separate marine refugia in the land-locked Sea of Japan basin and the Sea of Okhotsk. Existing mitochondrial DNA for the northeast Pacific populations is less conclusive regarding whether there were one or two separate refugia off this coast, and their location. Using environmental niche models this study shows the glacial NE Pacific environment could support two marine refugia, one centred in the Aleutians/Gulf of Alaska and the other off British Columbia. The intervening Cordilleran Ice Sheet, and the glacially sub-aerial and ice-free Queen Charlotte Islands shelf, is hypothesised to have constrained exchange between glacial stocks either side of the Islands. The postulated southern marine refugium is off-shore from an established terrestrial refugium, suggesting greater dependence of species and ecosystems during environmental change. An earth system approach to evolutionary change could enhance understanding of past and future ecosystems.     

Across the southern Andes on fin: glacial refugia,drainage reversals and a secondary contact zone revealed by the phylogeographical signal of Galaxias platei in Patagonia

We employed DNA sequence variation at two mitochondrial (control region, COI) regions from 212 individuals of Galaxias platei (Pisces, Galaxiidae) collected throughout Patagonia (25 lakes/rivers) to examine how Andean orogeny and the climatic cycles throughout the Quaternary affected the genetic diversity and phylogeography of this species. Phylogenetic analyses revealed four deep genealogical lineages which likely represent the initial division of G. platei into eastern and western lineages by Andean uplift, followed by further subdivision of each lineage into separate glacial refugia by repeated Pleistocene glacial cycles. West of the Andes, refugia were likely restricted to the northern region of Patagonia with small relicts in the south, whereas eastern refugia appear to have been much larger and widespread, consisting of separate northern and southern regions that collectively spanned most of Argentinean Patagonia. The retreat of glacial ice following the last glacial maximum allowed re‐colonization of central Chile from nonlocal refugia from the north and east, representing a region of secondary contact between all four glacial lineages. Northwestern glacial relicts likely followed pro‐glacial lakes into central Chilean Patagonia, whereas catastrophic changes in drainage direction (Atlantic → Pacific) for several eastern palaeolakes were the likely avenues for invasions from the east. These mechanisms, combined with evidence for recent, rapid and widespread population growth could explain the extensive contemporary distribution of G. platei throughout Patagonia.