Glacial refugia and recolonization pathways in the brown seaweed Fucus serratus

The last glacial maximum (20,000-18,000 years ago) dramatically affected extant distributions of virtually all northern European biota. Locations of refugia and postglacial recolonization pathways were examined in Fucus serratus (Heterokontophyta; Fucaceae) using a highly variable intergenic spacer developed from the complete mitochondrial genome of Fucus vesiculosus. Over 1,500 samples from the entire range of F. serratus were analysed using fluorescent single strand conformation polymorphism. A total of 28 mtDNA haplotypes was identified and sequenced. Three refugia were recognized based on high haplotype diversities and the presence of endemic haplotypes: southwest Ireland, the northern Brittany-Hurd Deep area of the English Channel, and the northwest Iberian Peninsula. The Irish refugium was the source for a recolonization sweep involving a single haplotype via northern Scotland and throughout Scandinavia, whereas recolonization from the Brittany-Hurd Deep refugium was more limited, probably because of unsuitable soft-bottom habitat in the Bay of Biscay and along the Belgian and Dutch coasts. The Iberian populations reflect a remnant refugium at the present-day southern boundary of the species range. A generalized skyline plot suggested exponential population expansion beginning in the mid-Pleistocene with maximal growth during the Eems interglacial 128,000-67,000 years ago, implying that the last glacial maximum mainly shaped population distributions rather than demography.     

Post-ice age recolonization and differentiation of Fucus serratus L. (Phaeophyceae; Fucaceae) populations in Northern Europe

The seaweed Fucus serratus is hypothesized to have evolved in the North Atlantic and present populations are thought to reflect recolonization from a southern refugium since the last glacial maximum 18 000-20 000 years bp. We examined genetic structure across several spatial scales by analysing seven microsatellite loci in populations collected from 21 localities throughout the species' range. Spatial auto-correlation analysis of seven microsatellite loci revealed no evidence for spatial clustering of alleles on a scale of 100 m despite limited gamete dispersal in F. serratus of approximately 2 m from parental individuals. Pairwise theta analysis suggested that the minimal panmictic unit for F. serratus was between 0.5 and 2 km. Isolation by distance was significant along some contiguous coastlines. Population differentiation was strong within the Skagerrak-Kattegat-Baltic Seas (SKB) (global theta= 0.17) despite a short history of approximately 7500 years. A neighbour-joining tree based on Reynold's distances computed from the microsatellite data revealed a central assemblage of populations on the Brittany Peninsula surrounded by four well-supported clusters consisting of the SKB, the North Sea (Ireland, Helgoland), and two populations from the northern Spanish coast. Samples from Iceland and Nova Scotia were most closely aligned with northwest Sweden and Brittany, respectively. When sample sizes were standardized (N = 41), allelic diversity was twofold higher for Brittany populations than for populations to the north and threefold higher than southern populations. The Brittany region may be a refugium or a recolonized area, whereas the Spanish populations most likely reflect present-day edge populations that have undergone repeated bottlenecks as a consequence of thermally induced cycles of recolonization and extinction.     

Genetic entities and mating system in hermaphroditic Fucus spiralis and its close dioecious relative F. vesiculosus (Fucaceae, Phaeophyceae)

To date, molecular markers have not settled the question of the specific status of the closely related, but phylogenetically unresolved, brown seaweeds, hermaphroditic Fucus spiralis and dioecious Fucus vesiculosus, nor their propensity for natural hybridization. To test the degree of species integrity and to assess effect of the mating system on the population genetic structure, 288 individuals coming from parapatric (discontinuous) and sympatric (contiguous) spatial configurations at two sites were genotyped with five microsatellite loci. Using a Bayesian admixture analysis, our results show that F. spiralis and F. vesiculosus comprise clearly distinct genetic entities (clusters) generally characterized by cosexual and unisexual individuals, respectively. Genetic diversity within each entity suggests that F. spiralis reproduces primarily through selfing while F. vesiculosus is characterized by an endogamous breeding regime. Nevertheless, aberrant sexual phenotypes were observed in each cluster, no diagnostic alleles were revealed and 10% of study individuals were intermediate between the two genetic entities. This pattern can be explained by recent divergence of two taxa with retention of ancestral polymorphism or asymmetrical, introgressive hybridization. However, given (i) coincident monomorphism at three loci in spiralis clusters and (ii) that significantly more intermediates were observed in sympatric stations than in parapatric stations, we argue that interspecific gene flow has occurred after divergence of the two taxa. Finally, we show that whether recently separated or recently introgressive, the divergent breeding systems probably contribute to species integrity in these two taxa.     

Unexpected High Genetic Diversity at the Extreme Northern Geographic Limit of Taurulus bubalis (Euphrasen, 1786)

The longspined bullhead (Taurulus bubalis, Euphrasen 1786) belongs to the family Cottidae and is a rocky shore species that inhabits the intertidal zones of the Eastern Atlantic since Iceland, southward to Portugal and also the North Sea and Baltic, northward to the Gulf of Finland, with some occurrences in the northern Mediterranean coasts eastward to the Gulf of Genoa. We analysed the phylogeographic patterns of this species using mitochondrial and nuclear markers in populations throughout most of its distributional range in west Europe. We found that T. bubalis has a relatively shallow genealogy with some differentiation between Atlantic and North Sea. Genetic diversity was homogeneous across all populations studied. The possibility of a glacial refugium near the North Sea is discussed. In many, but not all, marine temperate organisms, patterns of diversity are similar across the species range. If this phenomenon proves to be most common in cold adapted species, it may reflect the availability of glacial refugia not far from their present-day northern limits.     

Phylogeographic analysis of the red seaweed Palmaria palmata reveals a Pleistocene marine glacial refugium in the English Channel

Phylogeography has provided a new approach to the analysis of the postglacial history of a wide range of taxa but, to date, little is known about the effect of glacial periods on the marine biota of Europe. We have utilized a combination of nuclear, plastid and mitochondrial genetic markers to study the biogeographic history of the red seaweed Palmaria palmata in the North Atlantic. Analysis of the nuclear rDNA operon (ITS1-5.8S-ITS2), the plastid 16S-trnI-trnA-23S-5S, rbcL-rbcS and rpl12-rps31-rpl9 regions and the mitochondrial cox2-3 spacer has revealed the existence of a previously unidentified marine refugium in the English Channel, along with possible secondary refugia off the southwest coast of Ireland and in northeast North America and/or Iceland. Coalescent and mismatch analyses date the expansion of European populations from approximately 128,000 BP and suggest a continued period of exponential growth since then. Consequently, we postulate that the penultimate (Saale) glacial maximum was the main event in shaping the biogeographic history of European P. palmata populations which persisted throughout the last (Weichselian) glacial maximum (c. 20,000 BP) in the Hurd Deep, an enigmatic trench in the English Channel.     

Mitochondrial DNA variation in Pleistocene and modern Atlantic salmon from the Iberian glacial refugium

Current understanding of the postglacial colonization of Nearctic and Palearctic species relies heavily on inferences drawn from the phylogeographic analysis of contemporary generic variants. Modern postglacial populations are supposed to be representative of their Pleistocene ancestors, and their current distribution is assumed to reflect the different colonization success and dispersal patterns of refugial lineages. Yet, testing of phylogeographic models against ancestral genomes from glacial refugia has rarely been possible. Here we compare ND1 mitochondrial DNA variation in late Pleistocene (16,000-40,000 years before present), historical and contemporary Atlantic salmon (Salmo salar) populations from northern Spain and other regions of western Europe. Our study demonstrates the presence of Atlantic salmon in the Iberian glacial refugium during the last 40,000 years and points to the Iberian Peninsula as the likely source of the most common haplotype within the Atlantic lineage in Europe. However, our findings also suggest that there may have been significant changes in the genetic structure of the Iberian refugial stock since the last ice age, and question whether modern populations in refugial areas are representative of ice age populations. A common haplotype that persisted in the Iberian Peninsula during the Pleistocene last glacial maximum is now extremely rare or absent from European rivers, highlighting the need for caution when making phylogeographic inferences about the origin and distribution of modern genetic types.     

Diatom distribution and paleoceanographic reconstruction in the Southern Ocean — Present and last glacial maximum

Diatom assemblage and preservational data are used to reconstruct paleoceanographic conditions at the last glacial maximum (18,000 yrs BP) in the Southern Ocean. From these data, the following points can be made about the last glacial maximum in this region. (1) Contraction and slight northern shift of the belt of well preserved diatoms appears to be related to the northward shift of late spring/early summer sea ice cover. (2) Presence of open-ocean, though poorly preserved, diatom assemblages to the south of this belt strongly suggest that, during many summers, large areas of the Southern Ocean were ice-free. (3) The distribution of theNitzschia kerguelensis factor, both in surface sediments and at the last glacial maximum, indicates that the gyre systems, particularly the Weddell Gyre, were intensified during glacial times. (4) Although data are sparse in the higher-latitude South Atlantic, there is an indication that the Weddell Polynya also existed during glacial times, although it was shifted a few degrees northward.     

Microsatellite analysis of the phylogeography, Pleistocene history and secondary contact hypotheses for the killifish, Fundulus heteroclitus

The mummichog, Fundulus heteroclitus, exhibits extensive latitudinal clinal variation in a number of physiological and biochemical traits, coupled with phylogeographical patterns at mitochondrial and nuclear DNA loci that suggest a complicated history of spatially variable selection and secondary intergradation. This species continues to serve as a model for understanding local and regional adaptation to variable environments. Resolving the influences of historical processes on the distribution of genetic variation within and among extant populations of F. heteroclitus is crucial to a better understanding of how populations evolve in the context of contemporary environments. In this study, we analysed geographical patterns of genetic variation at eight microsatellite loci among 15 populations of F. heteroclitus distributed throughout the North American range of the species from Nova Scotia to Georgia. Genetic variation in Northern populations was lower than in Southern populations and was strongly correlated with latitude throughout the species range. The most common Northern alleles at all eight loci exhibited concordant latitudinal clinal patterns, and the existence of an abrupt transition zone in allele frequencies between Northern and Southern populations was similar to that observed for mitochondrial DNA and allozyme loci. A significant pattern of isolation by distance was observed both within and between northern and southern regions. This pattern was unexpected, particularly for northern populations, given the recent colonization history of post-Pleistocene habitats, and was inconsistent with either a recent northward population expansion or a geographically restricted northern Pleistocene refugium. The data provided no evidence for recent population bottlenecks, and estimates of historical effective population sizes suggest that post-Pleistocene populations have been large throughout the species distribution. These results suggest that F. heteroclitus was broadly distributed throughout most of its current range during the last glacial event and that the abrupt transition in allele frequencies that separate Northern and Southern populations may reflect regional disequilibrium conditions associated with the post-Pleistocene colonization history of habitats in that region.     

RECENT VERSUS RELIC: DISCERNING THE GENETIC SIGNATURE OF FUCUS VESICULOSUS (HETEROKONTOPHYTA; PHAEOPHYCEAE) IN THE NORTHWESTERN ATLANTIC1

Fucus vesiculosus L. is one of the most widespread macrophytes in the northwestern Atlantic, ranging from North Carolina (USA) to Greenland (DK). We investigated genetic diversity, population differentiation, patterns of isolation by distance, and putative glacial refugial populations across seven locations from North Carolina (USA) to Cape Breton Island, Nova Scotia (Canada), with microsatellite analyses. Distinct northern versus southern (Delaware–North Carolina) populations were revealed by microsatellite data. Five of six microsatellite loci were fixed in populations in North Carolina, suggesting a recent founder event or a bottleneck, and the same homozygous genotype was found in herbarium materials collected on the North Carolina coast from more than 60 years ago. An additional set of individuals from the northern limit in Greenland was included in our analysis of mitochondrial intergenic spacer (mt IGS) haplotypes in the northwestern Atlantic. Remarkably, 184 of 188 F. vesiculosus specimens from North Carolina to Greenland shared the same haplotype. Recent colonization of the North American shore from Europe is hypothesized based upon the ubiquity of this common haplotype, which was earlier reported from Europe.     

Adaptive traits are maintained on steep selective gradients despite gene flow and hybridization in the intertidal zone

Gene flow among hybridizing species with incomplete reproductive barriers blurs species boundaries, while selection under heterogeneous local ecological conditions or along strong gradients may counteract this tendency. Congeneric, externally-fertilizing fucoid brown algae occur as distinct morphotypes along intertidal exposure gradients despite gene flow. Combining analyses of genetic and phenotypic traits, we investigate the potential for physiological resilience to emersion stressors to act as an isolating mechanism in the face of gene flow. Along vertical exposure gradients in the intertidal zone of Northern Portugal and Northwest France, the mid-low shore species Fucus vesiculosus, the upper shore species Fucus spiralis, and an intermediate distinctive morphotype of F. spiralis var. platycarpus were morphologically characterized. Two diagnostic microsatellite loci recovered 3 genetic clusters consistent with prior morphological assignment. Phylogenetic analysis based on single nucleotide polymorphisms in 14 protein coding regions unambiguously resolved 3 clades; sympatric F. vesiculosus, F. spiralis, and the allopatric (in southern Iberia) population of F. spiralis var. platycarpus. In contrast, the sympatric F. spiralis var. platycarpus (from Northern Portugal) was distributed across the 3 clades, strongly suggesting hybridization/introgression with both other entities. Common garden experiments showed that physiological resilience following exposure to desiccation/heat stress differed significantly between the 3 sympatric genetic taxa; consistent with their respective vertical distribution on steep environmental clines in exposure time. Phylogenetic analyses indicate that F. spiralis var. platycarpus is a distinct entity in allopatry, but that extensive gene flow occurs with both higher and lower shore species in sympatry. Experimental results suggest that strong selection on physiological traits across steep intertidal exposure gradients acts to maintain the 3 distinct genetic and morphological taxa within their preferred vertical distribution ranges. On the strength of distributional, genetic, physiological and morphological differences, we propose elevation of F. spiralis var. platycarpus from variety to species level, as F. guiryi.     

Reconstruction of the colonization route from glacial refugium to the northern distribution range of the European butterfly Polyommatus coridon (Lepidoptera: Lycaenidae)

Genetic analyses are useful tools for reconstructing glacial distribution patterns and postglacial expansion corridors. However, little information is available about the latter. We reconstruct the expansion corridors of the butterfly Polyommatus coridon from its glacial refugium to the northern edge of its current distribution by comparing populations from southern Lower Saxony (central Germany) to other existing genetic data sets. The populations from Lower Saxony clearly belonged to a western lineage that expanded postglacially from the Adriatic‐Mediterranean region. They form part of a southern German group passing through the Burgundian Gap. In the southern German group, populations belong to a western subgroup. Therefore, expansion followed the Rhine valley and through Hesse, finally reaching southern Lower Saxony and western Thuringia in central Germany. Thus, we present a complete colonization route from the glacial refugium to the northern distribution range of P. coridon. Such data are useful for understanding the biogeographic structure and migration corridors for other mobile Mediterranean species.     

Genetic diversity, but not hatching success, is jointly affected by postglacial colonization and isolation in the threatened frog, Rana latastei

Both postglacial colonization and habitat fragmentation can reduce the genetic diversity of populations, which in turn can affect fitness. However, since these processes occur at different spatial and temporal scales, the consequences of either process may differ. To disentangle the relative role of isolation and postglacial colonization in determining genetic diversity and fitness, we studied microsatellite diversity of 295 individuals from 10 populations and measured the hatch rate of 218 clutches from eight populations of a threatened frog, R. latastei. The populations that were affected by fragmentation to a greater extent suffered higher embryo mortality and reduced hatch rate, while no effects of distance from glacial refugium on hatch rate were detected. Altogether, distance from glacial refugium and isolation explained > 90% of variation in genetic diversity. We found that the genetic diversity was lowest in populations both isolated and far from the glacial refugium, and that distance from refugium seems to have the primary role in determining genetic diversity. The relationship between genetic diversity and hatch rate was not significant. However, the proportion of genetic diversity lost through recent isolation had a significant, negative effect on fitness. It is possible that selection at least partially purged the negative effects of the ancestral loss of genetic diversity.     

Unique mitochondrial DNA lineages in Irish stickleback populations: cryptic refugium or rapid recolonization?

Repeated recolonization of freshwater environments following Pleistocene glaciations has played a major role in the evolution and adaptation of anadromous taxa. Located at the western fringe of Europe, Ireland and Britain were likely recolonized rapidly by anadromous fishes from the North Atlantic following the last glacial maximum (LGM). While the presence of unique mitochondrial haplotypes in Ireland suggests that a cryptic northern refugium may have played a role in recolonization, no explicit test of this hypothesis has been conducted. The three‐spined stickleback is native and ubiquitous to aquatic ecosystems throughout Ireland, making it an excellent model species with which to examine the biogeographical history of anadromous fishes in the region. We used mitochondrial and microsatellite markers to examine the presence of divergent evolutionary lineages and to assess broad‐scale patterns of geographical clustering among postglacially isolated populations. Our results confirm that Ireland is a region of secondary contact for divergent mitochondrial lineages and that endemic haplotypes occur in populations in Central and Southern Ireland. To test whether a putative Irish lineage arose from a cryptic Irish refugium, we used approximate Bayesian computation (ABC). However, we found no support for this hypothesis. Instead, the Irish lineage likely diverged from the European lineage as a result of postglacial isolation of freshwater populations by rising sea levels. These findings emphasize the need to rigorously test biogeographical hypothesis and contribute further evidence that postglacial processes may have shaped genetic diversity in temperate fauna.     

Signatures of vicariance, postglacial dispersal and spawning philopatry: population genetics of the walleye Sander vitreus

Population genetic relationships reveal the signatures of current processes such as reproductive behaviour and migration, as well as historic events including vicariance and climate change. We analyse population structure of native walleye Sander vitreus across North America, encompassing 10 nuclear DNA microsatellite loci, 26 spawning sites and 921 samples from watersheds across the Great Lakes, Lake Winnipeg, upper Mississippi River, Ohio River and Mobile Bay of the Gulf Coast. Geographical patterning is assessed using phylogenetic trees, pairwise F _ST analogues, hierarchical partitioning, Mantel regression, Bayesian assignment and Monmonier geographical networks. Results reveal congruent divergences among population groups, corresponding to historic isolation in glacial refugia, dispersal patterns and basin divisions. Broad-scale relationships show genetic isolation with geographical distance, but reproductive groups within basins do not – with some having pronounced differences. Greatest divergence distinguishes outlying Gulf Coastal and northwest populations, the latter tracing to dispersal from the Missourian refugium to former glacial Lake Agassiz, and basin isolation ∼7000 ya. Genetic barriers in the Great Lakes separate groups in Lakes Superior, Huron's Georgian Bay, Erie and Ontario, reflecting contributions from Mississippian and Atlantic refugia, and changes in connectivity patterns. Walleye genetic patterns thus reflect vicariance among watersheds and glacial refugia, followed by re-colonization pathways and changing drainage connections that established modern-day northern populations, whose separations are maintained through spawning site fidelity. Conservation management practices should preserve genetic identity and unique characters among these divergent walleye populations.     

Genetic consequences of Pleistocene glaciations for the tundra vole (Microtus oeconomus) in Beringia

Repeated glacial events during the Pleistocene fragmented and displaced populations throughout the northern continents. Different models of the effects of these climate-driven events predict distinct phylogeographic and population genetic outcomes for high-latitude faunas. The role of glaciations in (i) promoting intraspecific genetic differentiation and (ii) influencing genetic diversity was tested within a phylogeographic framework using the rodent Microtus oeconomus. The spatial focus for the study was Beringia, which spans eastern Siberia and northwestern North America, and was a continental crossroads and potential high arctic refugium during glaciations. Variation in mitochondrial DNA (cytochrome b and control region; 214 individuals) and nuclear DNA (ALDH1 intron; 63 individuals) was investigated across the Beringian region. Close genetic relationships among populations on either side of the Bering Strait are consistent with a history of periodic land connections between North America and Asia. A genetic discontinuity observed in western Beringia between members of a Central Asian clade and a Beringian clade is geographically congruent with glacial advances and with phylogeographic discontinuities identified in other organisms. Divergent island populations in southern Alaska were probably initially isolated by glacial vicariance, but subsequent differentiation has resulted from insularity. Tests of the genetic effects of postglacial colonization were largely consistent with expansion accompanied by founder effect bottlenecking, which yields reduced diversity in populations from recently deglaciated areas. Evidence that populations in the Beringian clade share a history of expansion from a low-diversity ancestral population suggests that Beringia was colonized by a small founder population from central Asia, which subsequently expanded in isolation.     

Convergent adaptation to a marginal habitat by homoploid hybrids and polyploid ecads in the seaweed genus Fucus

Hybridization and polyploidy are two major sources of genetic variability that can lead to adaptation in new habitats. Most species of the brown algal genus Fucus are found along wave-swept rocky shores of the Northern Hemisphere, but some species have adapted to brackish and salt marsh habitats. Using five microsatellite loci and mtDNA RFLP, we characterize two populations of morphologically similar, muscoides-like Fucus inhabiting salt marshes in Iceland and Ireland. The Icelandic genotypes were consistent with Fucus vesiculosus x Fucus spiralis F1 hybrids with asymmetrical hybridization, whereas the Irish ones consisted primarily of polyploid F. vesiculosus.     

Phylogeography of the West Indian manatee (Trichechus manatus): how many populations and how many taxa?

To resolve the population genetic structure and phylogeography of the West Indian manatee ( Trichechus manatus ), mitochondrial (mt) DNA control region sequences were compared among eight locations across the western Atlantic region. Fifteen haplotypes were identified among 86 individuals from Florida, Puerto Rico, the Dominican Republic, Mexico, Colombia, Venezuela, Guyana and Brazil. Despite the manatee's ability to move thousands of kilometres along continental margins, strong population separations between most locations were demonstrated with significant haplotype frequency shifts. These findings are consistent with tagging studies which indicate that stretches of open water and unsuitable coastal habitats constitute substantial barriers to gene flow and colonization. Low levels of genetic diversity within Florida and Brazilian samples might be explained by recent colonization into high latitudes or bottleneck effects. Three distinctive mtDNA lineages were observed in an intraspecific phylogeny of T. manatus , corresponding approximately to: (i) Florida and the West Indies; (ii) the Gulf of Mexico to the Caribbean rivers of South America; and (iii) the northeast Atlantic coast of South America. These lineages, which are not concordant with previous subspecies designations, are separated by sequence divergence estimates of d = 0.04–0.07, approximately the same level of divergence observed between T. manatus and the Amazonian manatee ( T. inunguis , n = 16). Three individuals from Guyana, identified as T. manatus , had mtDNA haplotypes which are affiliated with the endemic Amazon form T. inunguis . The three primary T. manatus lineages and the T. inunguis lineage may represent relatively deep phylogeographic partitions which have been bridged recently due to changes in habitat availability (after the Wisconsin glacial period, 10 000 BP ), natural colonization, and human-mediated transplantation.     

Post-glacial expansion into the Paleozoic Plateau: evidence of an Ozarkian refugium for the Ozark minnow Notropis nubilus (Teleostei: Cypriniformes)

Genetic variation was examined within the Ozark minnow Notropis nubilus using complete mtDNA cytochrome b gene sequences from 160 individuals representing 30 localities to test hypotheses on the origin of the distribution. Phylogenetic analyses revealed three strongly supported clades of haplotypes consistent with geographic distributions: a clade from the Western Ozarks, a clade from the Southern Ozarks and a clade from the Northern Ozarks and upper Mississippi River basin. The estimated mean ages of these clades indicated that they diverged during pre-Illinoian glacial cycles extending from the late Pliocene into the early Pleistocene. Results of demographic analyses based on coalescent approaches supported the hypothesis that the Paleozoic Plateau was not a refugium for N. nubilus during periodic glacial advances. There is evidence of a genetic signature of northern expansion into the Paleozoic Plateau from a Southern Ozarkian refugium. Populations expanded out of drainages in the Northern Ozarks into the Paleozoic Plateau during the late Pleistocene. Subsequently, the two regions were isolated due to the recent extirpation of intervening populations caused by the loss of suitable habitat.     

Multilocus phylogeography and population structure of common eiders breeding in North America and Scandinavia

Aim Glacial refugia during the Pleistocene had major impacts on the levels and spatial apportionment of genetic diversity of species in northern latitude ecosystems. We characterized patterns of population subdivision, and tested hypotheses associated with locations of potential Pleistocene refugia and the relative contribution of these refugia to the post‐glacial colonization of North America and Scandinavia by common eiders (Somateria mollissima). Specifically, we evaluated localities hypothesized as ice‐free areas or glacial refugia for other Arctic vertebrates, including Beringia, the High Arctic Canadian Archipelago, Newfoundland Bank, Spitsbergen Bank and north‐west Norway. Location Alaska, Canada, Norway and Sweden. Methods Molecular data from 12 microsatellite loci, the mitochondrial DNA (mtDNA) control region, and two nuclear introns were collected and analysed for 15 populations of common eiders (n = 716) breeding throughout North America and Scandinavia. Population genetic structure, historical population fluctuations and gene flow were inferred using F‐statistics, analyses of molecular variance, and multilocus coalescent analyses. Results Significant inter‐population variation in allelic and haplotypic frequencies were observed (nuclear DNA F_ST = 0.004–0.290; mtDNA Φ_ST = 0.051–0.927). Whereas spatial differentiation in nuclear genes was concordant with subspecific designations, geographic proximity was more predictive of inter‐population variance in mitochondrial DNA haplotype frequency. Inferences of historical population demography were consistent with restriction of common eiders to four geographic areas during the Last Glacial Maximum: Belcher Islands, Newfoundland Bank, northern Alaska and Svalbard. Three of these areas coincide with previously identified glacial refugia: Newfoundland Bank, Beringia and Spitsbergen Bank. Gene‐flow and clustering analyses indicated that the Beringian refugium contributed little to common eider post‐glacial colonization of North America, whereas Canadian, Scandinavian and southern Alaskan post‐glacial colonization is likely to have occurred in a stepwise fashion from the same glacial refugium. Main conclusions Concordance of proposed glacial refugia used by common eiders and other Arctic species indicates that Arctic and subarctic refugia were important reservoirs of genetic diversity during the Pleistocene. Furthermore, suture zones identified at MacKenzie River, western Alaska/Aleutians and Scandinavia coincide with those identified for other Arctic vertebrates, suggesting that these regions were strong geographic barriers limiting dispersal from Pleistocene refugia.