Glacial vicariance and historical biogeography of rock ptarmigan (Lagopus mutus) in the Bering region

In this paper, we address alternative hypotheses for the evolution of subspecies of rock ptarmigan (Lagopus mutus) endemic to the Aleutian Archipelago. To do this we examined patterns of genetic differentiation among populations of rock ptarmigan in the Aleutian Islands and parts of both Alaska and Siberia. Variation in mitochondrial control region sequences of 105 rock ptarmigan from 10 subspecies within the Bering region revealed three major phylogenetic lineages, two of which are endemic to the Aleutian Islands. Accordingly, haplotype and nucleotide diversities of rock ptarmigan within the archipelago are much higher than within mainland Alaska or Siberia. For Aleutian rock ptarmigan, analyses of molecular variance indicated significant genetic structuring and low estimates of gene flow among populations, despite small interisland distances within the archipelago. However, isolation by distance did not describe the pattern of gene flow or differentiation at this scale. Our estimates of divergence times of lineages suggest that Aleutian rock ptarmigan became isolated prior to the most recent Pleistocene glaciation event (late Wisconsin Stade) and that current patterns of genetic variation reflect the postglacial redistribution of divergent lineages and subsequent limited gene flow. In addition, genetic divergence among lineages was concordant with the distribution of plumage types among subspecies. The patterns of genetic variation described here for rock ptarmigan provide evidence for the role of glacial vicariance in contributing to genetic diversity within this and other Bering region species.     

Phylogeography and the origins of range disjunctions in a north temperate fish,the pygmy whitefish (Prosopium coulterii), inferred from mitochondrial and nuclear DNA sequence analysis

Aim To investigate the degree of phylogeographical divergence within pygmy whitefish (Prosopium coulterii) and to test hypotheses concerning the origin of disjunct populations within North America. Location North America from western Alaska to Lake Superior. Methods Mitochondrial (ATPase subunit VI) and nuclear (ITS‐1, ITS‐2) DNA sequence variation was assessed across the species’ North American range to test for the existence of distinct phylogeographical groupings of pygmy whitefish associated with known glacial refugia. Coalescent simulations of the mitochondrial DNA (mtDNA) data were used to test hypotheses of population structure. Results This species is composed of two monophyletic mitochondrial clades across its North American range. The two mtDNA clades differed by an average 3.3% nucleotide sequence divergence. These clades were also distinguished by ITS‐2, but the relationships among lineages were not resolved by the ITS‐1 analysis. Coalescent analyses rejected the null hypothesis of the current disjunct distributions being a result of fragmentation of a single widespread ancestral lineage across a variety of effective population sizes and divergence times. Main conclusions The current range disjunctions of pygmy whitefish in North America probably resulted from isolation, genetic divergence, and selective dispersal from at least two major Pleistocene glacial refugia: Beringia and Cascadia. More recent isolation and dispersal from an upper Mississippi refugium is suggested by relationships within one of the clades and by distributional evidence from co‐distributed species. The Beringian and Cascadian refugia have played major roles in the zoogeography of Nearctic temperate aquatics, but the roles of smaller refugia appear more variable among other species.     

Late Pleistocene divergence between eastern and western populations of wood ducks (Aix sponsa) inferred by the 'isolation with migration' coalescent method

During the Late Pleistocene, glaciers sundered many species into multiple glacial refugia where populations diverged in allopatry. Although deeply divergent mitochondrial DNA (mtDNA) lineages often reflect the number of refugia occupied, it is unlikely that populations that split during the recent Wisconsin glaciations will have reached reciprocal monophyly. We examined mtDNA control region sequences from eastern and western populations of wood ducks (Aix sponsa) to determine whether their current, disjunct distribution is consistent with the occupancy of two glacial refugia. We used the 'isolation with migration' coalescent method (im) to simultaneously estimate effective population sizes, maternal gene flow, and time since divergence. We found 24 unique haplotypes, none of which were shared between the eastern and western populations, but we did not find diagnostic monophyletic lineages suggestive of long-term isolation in multiple glacial refugia. However, a high Phi ST (0.31) indicates that eastern and western populations are well differentiated in mtDNA, and results from im suggest that these populations have been diverging, without extensive gene flow, for 10,000 to 124,000 years. Results from im further suggest that these populations most likely split about 34,000 years ago, and this time of divergence is consistent with the occupancy of multiple glacial refugia during the Late Wisconsin glaciation. Eastern wood ducks are characterized by high genetic diversity, a large effective population size, and a recent population expansion, while western wood ducks have much less genetic diversity, a smaller population size, and have not undergone a recent population expansion.     

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.     

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.     

Morphometric variation in the Nearctic collared lemming (Dicrostonyx)

A. H. Macpherson suggested that much of the current geographic diversity in Canadian Arctic mammals resulted from isolation in refugia during the Wisconsin glacial stage. This study evaluates the refugium hypothesis, insofar as it applies to Nearctic Dicrostonyx , by means of a statistical analysis of geographic variation in 13 skull characters. Overall, geographic variation among samples is not significant, although D. hudsonius and D. unalascensis are geographically and morphologically distinct. Some variation in skull shape is correlated with winter temperature. Partitioning tests on other measures of shape variation indicate some discontinuities consistent with the refugial hypothesis. Discrete samples reflect possible refugial populations in northern North America, Eastern Beringia and two southern periglacial refugia, one in eastern North America and at least one in western North America.     

Phylogeography of endemic ermine (Mustela erminea) in southeast Alaska

The North Pacific Coast (NPC) of North America is a region of high mammalian endemism, possibly due to its highly fragmented landscape and complex glacial history. For example, four island and one mainland subspecies of ermine, Mustela erminea, have been described as endemic to southeast Alaska alone. To better understand the role of past climatic change in generating diversity in the region, we examined DNA sequence variation in the mitochondrial cytochrome b gene of 210 ermine from across North America, with an emphasis on Alaska and British Columbia. We found three distinct (1.5-3.6% uncorrected 'p') lineages of ermine, all of which occur in southeast Alaska. One lineage includes a southeast Alaska endemic and specimens from Alaska (outside of southeast) and Eurasia. A second lineage includes two southeast Alaskan endemics and ermine from western Canada and the coterminous United States. The close relationships of these purported endemics to ermine outside of southeast Alaska suggest that they colonized the region from Beringian and southern glacial refugia, respectively, following deglaciation of the NPC. The third lineage appears restricted to the Prince of Wales Island complex in southeast Alaska (two subspecies) and Graham Island (Haida Gwaii), British Columbia. This restricted distribution suggests that these populations may be derived from relicts that persisted in a coastal refugium during the Wisconsin glaciation. Studies of nuclear genes and adaptive morphological evolution are necessary to further explore discrepancies between the geographical pattern of differentiation based on mtDNA and the existing subspecific taxonomy based on morphology.     

Coalescent-based hypothesis testing supports multiple Pleistocene refugia in the Pacific Northwest for the Pacific giant salamander (Dicamptodon tenebrosus)

Phylogeographic patterns of many taxa are explained by Pleistocene glaciation. The temperate rainforests within the Pacific Northwest of North America provide an excellent example of this phenomenon, and competing phylogenetic hypotheses exist regarding the number of Pleistocene refugia influencing genetic variation of endemic organisms. One such endemic is the Pacific giant salamander, Dicamptodon tenebrosus. In this study, we estimate this species' phylogeny and use a coalescent modeling approach to test five hypotheses concerning the number, location and divergence times of purported Pleistocene refugia. Single refugium hypotheses include: a northern refugium in the Columbia River Valley and a southern refugium in the Klamath-Siskiyou Mountains. Dual refugia hypotheses include these same refugia but separated at varying times: last glacial maximum (20,000 years ago), mid-Pleistocene (800,000 years ago) and early Pleistocene (1.7 million years ago). Phylogenetic analyses and inferences from nested clade analysis reveal distinct northern and southern lineages expanding from the Columbia River Valley and the Klamath-Siskiyou Mountains, respectively. Results of coalescent simulations reject both single refugium hypotheses and the hypothesis of dual refugia with a separation date in the late Pleistocene but not hypotheses predicting dual refugia with separation in early or mid-Pleistocene. Estimates of time since divergence between northern and southern lineages also indicate separation since early to mid-Pleistocene. Tests for expanding populations using mismatch distributions and 'g' distributions reveal demographic growth in the northern and southern lineages. The combination of these results provides strong evidence that this species was restricted into, and subsequently expanded from, at least two Pleistocene refugia in the Pacific Northwest.     

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

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

Rapid divergence and postglacial colonization in western North American Steller's jays (Cyanocitta stelleri)

Post-Pleistocene avian colonization of deglaciated North America occurred from multiple refugia, including a coastal refugium in the northwest. The location of a Pacific Coastal refugium is controversial; however, multiple lines of evidence suggest that it was located near the Queen Charlotte Islands (also known as Haida Gwaii). The Queen Charlotte Islands contain a disproportionately large number of endemic plants and animals including the Steller's jay Cyanocitta stelleri carlottae. Using five highly variable microsatellite markers, we studied population structure among eight populations of Steller's jay (N = 150) from geographical areas representing three subspecies in western North America: C. s. carlottae, C. s. stelleri and C. s. annectens. Microsatellite analyses revealed genetic differentiation between each of the three subspecies, although more extensive sampling of additional C. s. annectens populations is needed to clarify the level of subspecies differentiation. High levels of population structure were found among C. s. stelleri populations with significant differences in all but two pairwise comparisons. A significant isolation by distance pattern was observed amongst populations in the Pacific Northwest and Alaska. In the C. s. carlottae population, there was evidence of reduced genetic variation, higher number of private alleles than northern C. s. stelleri populations and higher levels of divergence between Queen Charlotte Island and other populations. We were unable to reject the hypothesis that the Queen Charlotte Islands served as a refugium during the Pleistocene. Steller's jay may have colonized the Queen Charlotte Islands near the end of the last glaciation or persisted throughout the Pleistocene, and this subspecies may thus represent a glacial relic. The larger number of private alleles, despite reduced genetic variation, morphological distinctiveness and high divergence from other populations suggests that the Queen Charlotte Island colonization pre-dates that of the mainland. Furthermore, our results show rapid divergence in Steller's jay populations on the mainland following the retreat of the ice sheets.     

Glacial vicariance in the Pacific Northwest: evidence from a lodgepole pine mitochondrial DNA minisatellite for multiple genetically distinct and widely separated refugia

The Canadian side of the Pacific Northwest was almost entirely covered by ice during the last glacial maximum, which has induced vicariance and genetic population structure for several plant and animal taxa. Lodgepole pine ( Pinus contorta Dougl. ex. Loud.) has a wide latitudinal and longitudinal distribution in the Pacific Northwest. Our main objective was to identify relictual signatures of glacial vicariance in the population structure of the species and search for evidence of distinct glacial refugia in the Pacific Northwest. A maternally inherited mitochondrial DNA minisatellite-like marker was used to decipher haplotype diversity in 91 populations of lodgepole pine located across the natural range. Overall population differentiation was sizeable ( G _ST = 0.365 and R _ST = 0.568). Four relatively homogeneous groups of populations, possibly representative of as many genetically distinct glacial populations, were identified for the two main subspecies, ssp. latifolia and ssp. contorta . For ssp. contorta , one glacial lineage is suggested to have been located at high latitudes and possibly off the coast of mainland British Columbia (BC), while the other is considered to have been located south of the ice sheet along the Pacific coast. For ssp. latifolia , two genetically distinct glacial populations probably occurred south of the ice sheet: in the area bounded by the Cascades and Rocky Mountains ranges, and on the eastern side of the Rockies. A possible fifth refugium located in the Yukon may have also been present for ssp. latifolia . Zones of contact between these ancestral lineages were also apparent in interior and northern BC. These results indicate the role of the Queen Charlotte Islands and the Alexander Archipelago as a refugial zone for some Pacific Northwest species and the vicariant role played by the Cascades and the American Rocky Mountains during glaciation.     

Mountainous genus Anterastes (Orthoptera,Tettigoniidae): autochthonous survival across several glacial ages via vertical range shifts

Although the high‐latitude range margins in Europe and North America are intensively studied, attention is gradually turned towards the taxa/populations inhabiting glacial refugia. Here, we evaluate the genealogical history of the cold‐adapted Anatolio‐Balkan genus Anterastes especially to test the possible effects of intrarefugial vertical range shifts during climatic oscillations of the Quaternary. Using concatenated data from sequences of COI+16S and ITS1–5.8S–ITS2, intrageneric relationships and the time of speciation events were estimated. Thirteen different demographic analyses were performed using a data set produced from sequences of 16S. Different phylogenetic analyses recovered similar lineages with high resolution. The molecular chronogram estimated speciation events in a period ranging from 5.60 to 1.22 Myr. Demographic analyses applied to 13 populations and five lineages suggested constant population size. Genetic diversity is significantly reduced in a few populations, while not in others. Fixation indices suggested extremely diverged populations. In the light of these data, the following main conclusions were raised: (i) although glacial refugia are the biodiversity hotspots, species level radiation of the cold‐adapted lineages is mainly prior to the Mid‐Pleistocene transition; (ii) heterogeneous topography provides refugial habitats and allows populations to survive through vertical range shifts during climatic fluctuations; (iii) prolonged isolation of refugial populations do not always result in reduced intrapopulation diversity, but in high level of genetic differentiation; (iv) the cold‐adapted lineages with low dispersal ability might have not colonised the area out of Anatolian refugium during interglacial periods; and (v) populations of invertebrates may have restricted ranges, but this does not mean that they have small effective population size.     

Postglacial colonisation patterns and the role of isolation and expansion in driving diversification in a passerine bird

Pleistocene glacial cycles play a major role in diversification and speciation, although the relative importance of isolation and expansion in driving diversification remains debated. We analysed mitochondrial DNA sequence data from 15 great reed warbler (Acrocephalus arundinaceus) populations distributed over the vast Eurasian breeding range of the species, and revealed unexpected postglacial expansion patterns from two glacial refugia. There were 58 different haplotypes forming two major clades, A and B. Clade A dominated in Western Europe with declining frequencies towards Eastern Europe and the Middle East, but showed a surprising increase in frequency in Western and Central Asia. Clade B dominated in the Middle East, with declining frequencies towards north in Central and Eastern Europe and was absent from Western Europe and Central Asia. A parsimonious explanation for these patterns is independent postglacial expansions from two isolated refugia, and mismatch distribution analyses confirmed this suggestion. Gene flow analyses showed that clade A colonised both Europe and Asia from a refugium in Europe, and that clade B expanded much later and colonised parts of Europe from a refugium in the Middle East. Great reed warblers in the eastern parts of the range have slightly paler plumage than western birds (sometimes treated as separate subspecies; A. a. zarudnyi and A. a. arundinaceus, respectively) and our results suggest that the plumage diversification took place during the easterly expansion of clade A. This supports the postglacial expansion hypothesis proposing that postglacial expansions drive diversification in comparatively short time periods. However, there is no indication of any (strong) reproductive isolation between clades and our data show that the refugia populations became separated during the last glaciation. This is in line with the Pleistocene speciation hypothesis invoking that much longer periods of time in isolation are needed for speciation to occur.     

Allopatric divergence and secondary contacts in <Emphasis Type="Italic">Euphorbia spinosa</Emphasis> L: Influence of climatic changes on the split of the species

Euphorbia spinosa, a perennial xerophilous shrub naturally distributed across the Italian peninsula, was selected for examination of the role of the Ligurian Alps and Apennines in glacial survival. The Italian Peninsula is considered to be one of the principal glacial refugia in Europe, but few plant population genetic and phylogeography studies have been undertaken in this region. The combined analysis of chloroplast and nuclear loci (ITS, cpSSR and ISSR) enabled us to detect extensive DNA variation and proved to be a very powerful tool for the reconstruction of the phylogeography. Molecular data support the hypothesis of a long-term separation of the Northwestern (Maritime Alps, Sardinia, Corsica, Northern Apennines) and Southeastern (Southern Apennines and Balkan area) lineages in glacial refugia. The existence of allopatrically fragmented lineages is most probably the result of isolation in different glacial refugia, possibly due to the Last Glacial Maximum cooling and the topographic complexity of the Italian peninsula. The most plausible hypothesis assumes the formation of two migration paths during more recent periods: the first one starting with southward migration and the second one moving northwards. The Central Apennines should be considered the confluence of migration routes radiating from separate refugia according to this hypothesis.     

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

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

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.     

Phylogeography of Potamon ibericum (Brachyura: Potamidae) identifies Quaternary glacial refugia within the Caucasus biodiversity hot spot

Refugia are critical for the maintenance of biodiversity during the periods of Quaternary climatic oscillations. The long‐term persistence of refugial populations in a large continuous refugium has resulted in a homogenous pattern of genetic structure among populations, while highly structured evolutionary lineages characterize the restriction of refugial populations to smaller subrefugia. These mechanisms have resulted in the identification of hot spots of biodiversity within putative glacial refugia. We studied phylogeography of Potamon ibericum (Brachyura: Potamidae) in the drainages of the western Caucasus biodiversity hot spot (i.e., Colchis and the Caucasus) to infer spatial genetic structure and potential refugia for a freshwater crab in this region. These areas have traditionally considered as a refugium due to the presence of Tertiary relict species. We integrated population genetic data and historical demographic analysis from cytochrome oxidase subunit I sequences and paleoclimatic data from species distribution modeling (SDM). The results revealed the lack of phylogeographic structure and provided evidence for demographic expansion. The SDM presented a rather homogenous and large refugium that extended from northeast Turkey to Colchis during the last glacial period. In contrast to these findings, previous phylogeographic study on P. ibericum of the eastern Caucasus biodiversity hot spot (i.e., Hyrcania) identified multiple independent refugia. By combining these results, we explain the significance of this important western Palearctic hot spot of biological diversity in shaping the geographic distribution of intraspecific genetic diversity in a freshwater taxon.     

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.     

Islands in the ice: colonisation routes for rock ptarmigan to the Svalbard archipelago

The ancestors of rock ptarmigan Lagopus muta originated in the Beringia area well before the disruption of the Beringian land bridge. They spread westwards to Siberia and eastwards to the North American arctic and Greenland. The distribution of rock ptarmigan has been affected by glaciations restricting it to geographically limited refugia. Today the species has a circumpolar distribution in arctic tundra and alpine habitats, with up to 30 subspecies recognised based on morphological characters. We sequenced the mitochondrial control region for 72 individuals and genotyped 69 individuals for 12 microsatellite loci to investigate neutral genetic variation within and among five rock ptarmigan populations, Greenland, Iceland, Scandinavia, Svalbard and Taimyr. Genetic structure among the studied samples was high, overall F_ST estimated from microsatellite loci was 0.18 and only one out of 16 mtDNA haplotypes was found in more than one population. Genetic variation ( h, π, H_e, allelic richness) was slightly lower in the Svalbard population than in other populations, suggesting a low effective population size, possibly due to isolation following colonisation. An unrooted network and a phylogenetic tree showed that the Scandinavian population has diverged from the other populations by at least ten mutational steps, probably due to independent colonisation of Europe and subsequent long-term isolation, and rules out Scandinavia as a source for colonisation of Svalbard. Alignment with partial control region sequences from other studies showed that the haplotype that was central in our network and found on Svalbard and in Taimyr, most likely corresponds to a haplotype found in Siberia, Alaska and the Canadian Arctic, but not in Greenland, Scandinavia and Iceland. This suggests an eastern origin of rock ptarmigan in Svalbard, although this question cannot be settled conclusively.     

Phylogeography of ninespine sticklebacks (Pungitius pungitius) in North America: glacial refugia and the origins of adaptive traits

The current geographical distribution of the ninespine stickleback (Pungitius pungitius) was shaped in large part by the glaciation events of the Pleistocene epoch (2.6 Mya–10 Kya). Previous efforts to elucidate the phylogeographical history of the ninespine stickleback in North America have focused on a limited set of morphological traits, some of which are likely subject to widespread convergent evolution, thereby potentially obscuring relationships among populations. In this study, we used genetic information from both mitochondrial DNA (mtDNA) sequences and nuclear microsatellite markers to determine the phylogenetic relationships among ninespine stickleback populations. We found that ninespine sticklebacks in North America probably dispersed from at least three glacial refugia—the Mississippi, Bering, and Atlantic refugia—not two as previously thought. However, by applying a molecular clock to our mtDNA data, we found that these three groups diverged long before the most recent glacial period. Our new phylogeny serves as a critical framework for examining the evolution of derived traits in this species, including adaptive phenotypes that evolved multiple times in different lineages. In particular, we inferred that loss of the pelvic (hind fin) skeleton probably evolved independently in populations descended from each of the three putative North American refugia.