Genome‐wide SNP data reveal improved evidence for Antarctic glacial refugia and dispersal of terrestrial invertebrates

Antarctica is isolated, surrounded by the Southern Ocean and has experienced extreme environmental conditions for millions of years, including during recent Pleistocene glacial maxima. How Antarctic terrestrial species might have survived these glaciations has been a topic of intense interest, yet many questions remain unanswered, particularly for Antarctica's invertebrate fauna. We examine whether genetic data from a widespread group of terrestrial invertebrates, springtails (Collembola, Isotomidae) of the genus Cryptopygus, show evidence for long‐term survival in glacial refugia along the Antarctic Peninsula. We use genome‐wide SNP analyses (via genotyping‐by‐sequencing, GBS) and mitochondrial data to examine population diversity and differentiation across more than 20 sites spanning >950 km on the Peninsula, and from islands both close to the Peninsula and up to ~1,900 km away. Population structure analysis indicates the presence of strong local clusters of diversity, and we infer that patterns represent a complex interplay of isolation in local refugia coupled with occasional successful long‐distance dispersal events. We identified wind and degree days as significant environmental drivers of genetic diversity, with windier and warmer sites hosting higher diversity. Thus, we infer that refugial areas along the Antarctic Peninsula have allowed populations of indigenous springtails to survive in situ throughout glacial periods. Despite the difficulties of dispersal in cold, desiccating conditions, Cryptopygus springtails on the Peninsula appear to have achieved multiple long‐distance colonization events, most likely through wind‐related dispersal events.     

Molecular data and ecological niche modelling reveal the phylogeographic pattern of Cotinus coggygria (Anacardiaceae) in China's warm‐temperate zone

The phylogeography of common and widespread species helps to elucidate the history of local flora and vegetation. In this study, we selected Cotinus coggygria, a species widely distributed in China's warm‐temperate zone. One chloroplast DNA (cpDNA) region and ecological niche modelling were used to examine the phylogeographic pattern of C. coggygria. The cpDNA data revealed two phylogeographic groups (Southern and Northern) corresponding to the geographic regions. Divergence time analyses revealed that divergence of the two groups occurred at approximately 147,000 years before the present (BP), which coincided with the formation of the downstream area of the Yellow River, indicating that the Yellow River was a weak phylogeographic divide for C. coggygria. The molecular data and ecological niche modelling also indicated that C. coggyria did not experience population expansion after glaciations. This study thus supports the fact that Pleistocene glacial cycles only slightly affected C. coggygria, which survived in situ and occupied multiple localised glacial refugia during glaciations. This finding is contrary to the hypothesis of large‐scale range habitat contraction and retreat into a few main refugia.     

Evolutionary history of Apocheima cinerarius (Lepidoptera: Geometridae), a female flightless moth in northern China

The alterations in the phylogeographical structures of insects in response to the uplift of the Qinghai–Tibet Plateau and the Quaternary glaciations in eastern Asia, particularly in northern China, remain largely unknown. In this study, we selected Apocheima cinerarius, a moth with flightless females, using molecular data (complete mitochondrial genomes and nuclear data) and ecological niche modelling (ENM) to investigate the effects of paleoclimatic changes on the evolutionary history of insects in the area of northern China. The phylogenetic tree of complete mitochondrial genomes indicated that there were two lineages, the western and eastern lineages. The nuclear gene analyses also detected unique haplotypes in each lineage. Time of the most recent common ancestor (TMRCA) of the two lineages was approximately in Early–Middle Pleistocene. Bayesian skyline plots revealed that the western lineage underwent a population expansion event after the Last Glacial Maximum, whereas the eastern lineage underwent expansion between the Last Interglacial and the Last Glacial Maximum. Our results suggest that A. cinerarius expanded eastward from western sites until the moth was distributed across the entire region of northern China. Then, A. cinerarius underwent contraction into isolated glacial refugia followed by subsequent expansion driven by Pleistocene climate changes, which established a narrow sympatric area. Our results indicate that the Quaternary environmental fluctuations had profound influences on the diversification and demography of an insect in northern China, and the same species in north‐western China and north‐eastern China have different demographic histories.     

Pleistocene diversification in Morocco and recent demographic expansion in the Mediterranean pond turtle Mauremys leprosa

Quaternary climatic oscillations and geographic barriers have strongly influenced the distribution and diversification of thermophilic species occurring in the Mediterranean Basin. The Western Mediterranean pond turtle, Mauremys leprosa, is widely distributed throughout the Iberian Peninsula, southern France and most of the Maghreb region, with two subspecies currently recognized. In this work, we used 566 samples, including 259 new individuals, across the species range, and sequenced two mitochondrial markers (cytochrome b gene and control region; 163 samples in a concatenated mtDNA dataset) and one nuclear intron (R35; 23 samples representing all identified sublineages) to study the evolutionary history of M. leprosa. We combined phylogenetic methods and phylogeographic continuous diffusion models with spatial analysis. Our results (1) show a high level of genetic structure in Morocco originated during the Pleistocene; (2) reveal two independent population expansion waves from Morocco to Tunisia and to southern Europe, which later expanded throughout the Iberian Peninsula, and (3) identify several secondary contact zones in Morocco. Our study also sheds new light on the role of geographical features (Moroccan mountains ranges and the Strait of Gibraltar) and Pleistocene climatic oscillations in shaping genetic diversity and structure of M. leprosa, and underlines the importance of the Maghreb as a differentiation centre harbouring distinct glacial refugia.     

Mitochondrial phylogeography,contact zones and taxonomy of grass snakes (Natrix natrix,N. megalocephala)

Grass snakes (Natrix natrix) represent one of the most widely distributed snake species of the Palaearctic region, ranging from the North African Maghreb region and the Iberian Peninsula through most of Europe and western Asia eastward to the region of Lake Baikal in Central Asia. Within N. natrix, up to 14 distinct subspecies are regarded as valid. In addition, some authors recognize big‐headed grass snakes from western Transcaucasia as a distinct species, N. megalocephala. Based on phylogenetic analyses of a 1984‐bp‐long alignment of mtDNA sequences (ND4+tRNAs, cyt b) of 410 grass snakes, a nearly range‐wide phylogeography is presented for both species. Within N. natrix, 16 terminal mitochondrial clades were identified, most of which conflict with morphologically defined subspecies. These 16 clades correspond to three more inclusive clades from (i) the Iberian Peninsula plus North Africa, (ii) East Europe and Asia and (iii) West Europe including Corso‐Sardinia, the Apennine Peninsula and Sicily. Hypotheses regarding glacial refugia and postglacial range expansions are presented. Refugia were most likely located in each of the southern European peninsulas, Corso‐Sardinia, North Africa, Anatolia and the neighbouring Near and Middle East, where the greatest extant genetic diversity occurs. Multiple distinct microrefugia are inferred for continental Italy plus Sicily, the Balkan Peninsula, Anatolia and the Near and Middle East. Holocene range expansions led to the colonization of more northerly regions and the formation of secondary contact zones. Western Europe was invaded from a refuge within southern France, while Central Europe was reached by two distinct range expansions from the Balkan Peninsula. In Central Europe, there are two contact zones of three distinct mitochondrial clades, and one of these contact zones was theretofore completely unknown. Another contact zone is hypothesized for Eastern Europe, which was colonized, like north‐western Asia, from the Caucasus region. Further contact zones were identified for southern Italy, the Balkans and Transcaucasia. In agreement with previous studies using morphological characters and allozymes, there is no evidence for the distinctiveness of N. megalocephala. Therefore, N. megalocephala is synonymized with N. natrix.     

Refugia,colonization and diversification of an arid‐adapted bird: coincident patterns between genetic data and ecological niche modelling

Phylogeographical studies are common in boreal and temperate species from the Palaearctic, but scarce in arid‐adapted species. We used nuclear and mitochondrial markers to investigate phylogeography and to estimate chronology of colonization events of the trumpeter finch Bucanetes githagineus, an arid‐adapted bird. We used 271 samples from 16 populations, most of which were fresh samples but including some museum specimens. Microsatellite data showed no clear grouping according to the sampling locations. Microsatellite and mitochondrial data showed the clearest differentiation between Maghreb and Canary Islands and between Maghreb and Western Sahara. Mitochondrial data suggest differentiation between different Maghreb populations and among Maghreb and Near East populations, between Iberian Peninsula and Canary Islands, as well as between Western Sahara and Maghreb. Our coalescence analyses indicate that the trumpeter finch colonized North Africa during the humid Marine Isotope Stage 5 (MIS5) period of the Sahara region 125 000 years ago. We constructed an ecological niche model (ENM) to estimate the geographical distribution of climatically suitable habitats for the trumpeter finch. We tested whether changes in the species range in relation to glacial–interglacial cycles could be responsible for observed patterns of genetic diversity and structure. Modelling results matched with those from genetic data as the species' potential range increases in interglacial scenarios (in the present climatic scenario and during MIS5) and decreases in glacial climates (during the last glacial maximum, LGM, 21 000 years ago). Our results suggest that the trumpeter finch responded to Pleistocene climatic changes by expanding and contracting its range.     

Incomplete lineage sorting or secondary admixture: disentangling historical divergence from recent gene flow in the Vinous‐throated parrotbill (Paradoxornis webbianus)

Although Pleistocene glaciations had a major impact on the population genetic patterns of many species in North America and Europe, it remains unclear how these climatic fluctuations contributed to species diversification in East Asia. One reason for this is the difficulty of distinguishing genetic admixture following secondary contact from incomplete lineage sorting, both of which can generate similar patterns of genetic variation. Using a combination of multilocus analyses and coalescent simulation, we explore how these two processes occurred in the Pleistocene evolutionary history of a widespread East Asian bird, the Vinous‐throated parrotbill, Paradoxornis webbianus. Maximum likelihood (ML) tree identified two major mitochondrial lineages, which are geographically separated in most parts of its range, but are sympatric at a few sampling sites. NJ tree and Structure analysis of microsatellite data set revealed an extensive level of admixture and little population structure, suggesting recent admixture between two formerly separated groups. Networks from nuclear DNA data sets, however, did not indicate any geographically isolated groups but rather a panmictic population, thus support incomplete lineage sorting. By using coalescent simulation approaches, we show that both processes did occur, although at different temporal scales. During the Pleistocene glaciations, probably around 0.1–0.5 Ma (the Marine Isotope Stage 6, MIS6), P. webbianus contracted into two separate refugia, and subsequently accumulated genetic divergence. During the interglacial MIS5, the species expanded into previously glaciated areas allowing the once separated groups to come into contact and become admixed. Taken together, our results indicate the current genetic variation within P. webbianus is a combination pattern of widespread distribution in pre‐Pleistocene, then contraction and fragmentation into separated refugia during glacial advance, followed by recently postglacial expansion and admixture.     

Phylogeographic analysis reveals northerly refugia for the riverine amphibian Triturus dobrogicus (Caudata: Salamandridae)

We investigated the recent evolutionary history of the Danube crested newt, Triturus dobrogicus through reconstructions of: (1) the number and position of refugia at the last glacial maximum, (2) the role of major central European rivers in pattern of post‐glacial dispersal, and (3) the present‐day distribution pattern. We analysed sequences of mitochondrial DNA (ND2, 1065 bp) and six microsatellite loci in 363 T. dobrogicus individuals from 58 populations covering the range of the species. Our analyses suggested that T. dobrogicus survived the last glacial maximum in two separate refugia positioned in northwestern Pannonia and in Southern Pannonia from where its range expanded along the Danube and Tisza Rivers. Our findings also confirmed that rivers played an important role in shaping the evolutionary history of amphibian species in Central Europe. We compared the T. dobrogicus range with another lowland amphibian, the fire‐bellied toad Bombina bombina, using species distribution modelling. In line with these models, the isolated mountains inside Pannonia are occupied not by B. bombina, but by its mountain‐dwelling sister‐species B. variegata. However, in contrast to the model, crested newts (including T. dobrogicus) are absent from these mountains. We attribute this biogeographical discrepancy to the positioning of the species’ refugia at the last glacial maximum.     

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.     

Distributional and demographic consequences of Pleistocene climate fluctuations for a marine demersal fish in the north-eastern Atlantic

Aim The Pleistocene glaciations were the most significant historical event during the evolutionary life span of most extant species. However, little is known about the consequences of these climate changes for the distribution and demography of marine animals of the north‐eastern Atlantic. The present study focuses on the phylogeographic and demographic patterns of the sand goby, Pomatoschistus minutus (Teleostei: Gobiidae), a small marine demersal fish. Location North‐eastern Atlantic, Mediterranean, Irish, North and Baltic seas. Methods Analysis was carried out by sequencing the mtDNA cytochrome b gene of sand gobies from 12 localities throughout the species’ range, and using this information in combination with published data of allozyme markers and mtDNA control region sequences. Several phylogenetic methods and a network analysis were used to explore the phylogeographic pattern. The historical demography of P. minutus was studied through a mismatch analysis and a Bayesian skyline plot. Results Reciprocal monophyly was found between a Mediterranean Sea (MS) clade and an Atlantic Ocean (AO) clade, both with a Middle Pleistocene origin. The AO Clade contains two evolutionary significant units (ESUs): the Iberian Peninsula (IB) Group and the North Atlantic (NA) Group. These two groups diverged during Middle Pleistocene glacial cycles. For the NA Group there is evidence for geographic sorting of the ancestral haplotypes with recent radiations in the Baltic Sea, Irish Sea, North Sea and Bay of Biscay. The demographic histories of the Mediterranean Clade and the two Atlantic ESUs were influenced mainly by expansions dated as occurring during the Middle Pleistocene glaciations and post‐Eem, respectively. Main conclusions The pre‐LGM (Last Glacial Maximum) subdivision signals were not erased for P. minutus during the LGM. Middle Pleistocene glaciations yielded isolated and differently evolving sets of populations. In contrast to the case for most other taxa, only the northern Atlantic group contributed to the post‐glacial recolonization. The historical demography of Mediterranean sand gobies was influenced mainly by Middle Pleistocene glaciations, in contrast to that of the Atlantic populations, which was shaped by Late Pleistocene expansions.     

Molecular markers provide evidence for a broad‐fronted recolonisation of the widespread calcareous grassland species Sanguisorba minor from southern and cryptic northern refugia

• Calcareous grasslands belong to the most species‐rich and endangered habitats in Europe. However, little is known about the origin of the species typically occurring in these grasslands. In this study we analysed the glacial and post‐glacial history of Sanguisorba minor, a typical plant species frequently occurring in calcareous grasslands.
• The study comprised 38 populations throughout the whole distribution range of the species across Europe. We used molecular markers (AFLP) and applied Bayesian cluster analysis as well as spatial principal components analysis (sPCA) to identify glacial refugia and post‐glacial migration routes to Central Europe.
• Our study revealed significant differences in the level of genetic variation and the occurrence of rare fragments within populations of S. minor and a distinct separation of eastern and western lineages. The analyses uncovered traditional southern but also cryptic northern refugia and point towards a broad fronted post‐glacial recolonisation.
• Based on these results we postulate that incomplete lineage sorting may have contributed to the detected pattern of genetic variation and that S. minor recolonised Central Europe post‐glacially from Iberia and northern glacial refugia in France, Belgium or Germany. Our results highlight the importance of refugial areas for the conservation of intraspecific variation in calcareous grassland species.

     

Interglacial genetic diversification of Moussonia deppeana (Gesneriaceae), a hummingbird‐pollinated,cloud forest shrub in northern Mesoamerica

Recent empirical work on cloud forest‐adapted species supports the role of both old divergences across major geographical areas and more recent divergences attributed to Pleistocene climate changes. The shrub Moussonia deppeana is distributed in northern Mesoamerica, with geographically disjunct populations. Based on sampling throughout the species range and employing plastid and nuclear markers, we (i) test whether the fragmented distribution is correlated with main evolutionary lineages, (ii) reconstruct its phylogeographical history to infer the history of cloud forest in northern Mesoamerica and (iii) evaluate a set of refugia/vicariance scenarios for the region and demographic patterns of the populations whose ranges expanded and tracked cloud forest conditions during the Last Glacial Maximum. We found a deep evolutionary split in M. deppeana about 6–3 Ma, which could be consistent with a Pliocene divergence. Comparison of variation in plastid and nuclear markers revealed several lineages mostly congruent with their isolated geographical distribution and restricted gene flow among groups. Results of species distribution modelling and coalescent simulations fit a model of multiple refugia diverging during interglacial cycles. The demographic history of M. deppeana is not consistent with an expanding–contracting cloud forest archipelago model during the Last Glacial Maximum. Instead, our data suggest that populations persisted across the geographical range throughout the glacial cycles, and experienced isolation and divergence during interglacial periods.     

Population structure within an alpine archipelago: strong signature of past climate change in the New Zealand rock wren (Xenicus gilviventris)

Naturally subdivided populations such as those occupying high‐altitude habitat patches of the ‘alpine archipelago’ can provide significant insight into past biogeographical change and serve as useful models for predicting future responses to anthropogenic climate change. Among New Zealand's alpine taxa, phylogenetic studies support two major radiations: the first correlating with geological forces (Pliocene uplift) and the second with climatic processes (Pleistocene glaciations). The rock wren (Xenicus gilviventris) is a threatened alpine passerine belonging to the endemic New Zealand wren family (Acanthisittidae). Rock wren constitute a widespread, naturally fragmented population, occurring in patches of suitable habitat over c. 900 m in altitude throughout the length of the South Island, New Zealand. We investigated the relative role of historical geological versus climatic processes in shaping the genetic structure of rock wren (N = 134) throughout their range. Using microsatellites combined with nuclear and mtDNA sequence data, we identify a deep north–south divergence in rock wren (3.7 ± 0.5% at cytochrome b) consistent with the glacial refugia hypothesis whereby populations were restricted in isolated refugia during the Pleistocene c. 2 Ma. This is the first study of an alpine vertebrate to test and provide strong evidence for the glacial refugia hypothesis as an explanation for the low endemicity central zone known as the biotic ‘gap’ in the South Island of New Zealand.     

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.     

Patterns of chloroplast DNA variation in <Emphasis Type="Italic">Cycas debaoensis</Emphasis> (Cycadaceae): conservation implications

Climate changes during glacial periods have had significant effects on the current distribution of plant species. Palaeontologial data suggest that modern cycads originated in southwest China. Cycas debaoensis (Cycadaceae) is an endangered species restricted to a small area of southwest China. This species has been classified into two types: sand and karst, according to the soil matrix they grow on. To determine the locations of its glacial refugia and its genetic structure, we examined chloroplast sequence variation of the atpB-rbcL and psbA-trnH intergenic spacers. Four chloroplast DNA haplotypes were obtained from 120 individuals collected from 11 populations covering the entire extant distribution range of the species. Significant population subdivision was detected (G _ST = 0.684 and F _ST = 0.74160), suggesting low levels of gene flow between regions and populations. There was marked haplotype differentiation between populations in the sand and karst regions, with only one haplotype being present in both. The molecular phylogenetic data, together with the geographic distribution of the haplotypes, suggest that C. debaoensis experienced range contraction during glacial periods, and that the current populations are still confined to the areas of the original refugia. These results implied that isolated refugia might have maintained in both sand and karst regions during the last glacial maximum and even earlier glaciations. The low within-population diversity of C. debaoensis suggested that there were strong bottleneck events or founder effects within each separate region during the Quaternary climatic oscillations. These findings are important for the conservation of this endangered species.     

Tracking the ice: Subterranean harvestmen distribution matches ancient glacier margins

Biogeographic studies often underline the role of glacial dynamism during Pleistocene (1.806–0.011 Mya) in shaping the distribution of subterranean species. Accordingly, it is presumed that present‐day distribution of most specialized cold‐adapted (cryophilic) cave‐dwelling species should bear the signatures of past climatic events. To test this idea, we modelled the distribution of specialized cold‐adapted subterranean alpine harvestmen (Arachnida: Opiliones: Ischyropsalididae: Ischyropsalis). We found that the distance from the glacier margins during Last Glacial Maximum (LGM; about 22,000 years ago) was the most important predictor of their present‐day distribution. In particular, the peak in the probability of occurrence of alpine subterranean Ischyropsalis was found to be in close proximity to the LGM glacier, with a sharp drop at a distance of 30 km from the ice margin. In light of the role played by past climatic events in determining the species current range, we briefly discuss their biogeographic history and the role played by glacial refugia dynamics in determining the current distribution of these species. We argue that low dispersal harvestmen such as our model species can be used as biological indicators for tracking past glaciations and other similar biogeographic events.     

Chromosome polymorphisms track trans‐Atlantic divergence and secondary contact in Atlantic salmon

Pleistocene glaciations drove repeated range contractions and expansions shaping contemporary intraspecific diversity. Atlantic salmon (Salmo salar) in the western and eastern Atlantic diverged >600,000 years before present, with the two lineages isolated in different southern refugia during glacial maxima, driving trans‐Atlantic genomic and karyotypic divergence. Here, we investigate the genomic consequences of glacial isolation and trans‐Atlantic secondary contact using 108,870 single nucleotide polymorphisms genotyped in 80 North American and European populations. Throughout North America, we identified extensive interindividual variation and discrete linkage blocks within and between chromosomes with known trans‐Atlantic differences in rearrangements: Ssa01/Ssa23 translocation and Ssa08/Ssa29 fusion. Spatial genetic analyses suggest independence of rearrangements, with Ssa01/Ssa23 showing high European introgression (>50%) in northern populations indicative of post‐glacial trans‐Atlantic secondary contact, contrasting with low European ancestry genome‐wide (3%). Ssa08/Ssa29 showed greater intrapopulation diversity, suggesting a derived chromosome fusion polymorphism that evolved within North America. Evidence of potential selection on both genomic regions suggests that the adaptive role of rearrangements warrants further investigation in Atlantic salmon. Our study highlights how Pleistocene glaciations can influence large‐scale intraspecific variation in genomic architecture of northern species.     

Population genomic evidence for plant glacial survival in Scandinavia

Quaternary glaciations have played a major role in shaping the genetic diversity and distribution of plant species. Strong palaeoecological and genetic evidence supports a postglacial recolonization of most plant species to northern Europe from southern, eastern and even western glacial refugia. Although highly controversial, the existence of small in situ glacial refugia in northern Europe has recently gained molecular support. We used genomic analyses to examine the phylogeography of a species that is critical in this debate. Carex scirpoidea Michx subsp. scirpoidea is a dioecious, amphi‐Atlantic arctic–alpine sedge that is widely distributed in North America, but absent from most of Eurasia, apart from three extremely disjunct populations in Norway, all well within the limits of the Weichselian ice sheet. Range‐wide population sampling and variation at 5,307 single nucleotide polymorphisms show that the three Norwegian populations comprise unique evolutionary lineages divergent from Greenland with high between‐population divergence. The Norwegian populations have low within‐population genetic diversity consistent with having experienced genetic bottlenecks in glacial refugia, and host private alleles that probably accumulated in long‐term isolated populations. Demographic analyses support a single, pre‐Weichselian colonization into Norway from East Greenland, and subsequent divergence of the three populations in separate refugia. Other refugial areas are identified in North‐east Greenland, Minnesota/Michigan, Colorado and Alaska. Admixed populations in British Columbia and West Greenland indicate postglacial contact. Taken together, evidence from this study strongly indicates in situ glacial survival in Scandinavia.     

A TEST OF THE GLACIAL REFUGIUM HYPOTHESIS USING PATTERNS OF MITOCHONDRIAL AND NUCLEAR DNA SEQUENCE VARIATION IN ROCK PTARMIGAN (LAGOPUS MUTUS)

The glacial refugium hypothesis (GRH) proposes that glaciers promoted differentiation and generation of intraspecific diversity by isolating populations in ice-free refugia. We tested three predictions of this hypothesis for the evolutionary divergence of rock ptarmigan (Lagopus mutus) during the Wisconsin glaciation of the late Pleistocene. To do this, we examined subspecies distributions, population genetic structure, and phylogenetic relationships in 26 populations across North America and the Bering Sea region. First, we analyzed sequence variation in the mitochondrial control region, in a nuclear intron (Gapdh), and in an internal transcribed spacer (ITS1). Control region sequences of 154 rock ptarmigan revealed strong population and phylogeographic structure. Variation in intron sequences of 114 rock ptarmigan also revealed significant population structure compatible with results for the control region. Rock ptarmigan were invariant for ITS1. Second, we show that five known Nearctic refugia and an Icelandic refugium are concordant with the current distribution of morphologically distinct subspecies; five of these six refugia are geographically concordant with the distribution of closely related control region haplotypes. Third, our estimates of the time since phylogenetic lineages diverged predated the last glacial maximum for all but two lineages. In addition, all lines of evidence suggest that two unknown refugia in the Bering Sea region supported rock ptarmigan during the Wisconsin glaciation. Overall, our results are most consistent with the hypothesis that isolated populations of rock ptarmigan diverged in multiple refugia during the Wisconsin and that geographic variation reflects patterns of recolonization of the Nearctic after the ice receded. The GRH may therefore offer the most plausible explanation for similar biogeographic patterns in a variety of Nearctic vertebrates.