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.     

Relative information content of polymorphic microsatellites and mitochondrial DNA for inferring dispersal and population genetic structure in the olive sea snake, Aipysurus laevis

Polymorphic microsatellites are widely considered more powerful for resolving population structure than mitochondrial DNA (mtDNA) markers, particularly for recently diverged lineages or geographically proximate populations. Weaker population subdivision for biparentally inherited nuclear markers than maternally inherited mtDNA may signal male-biased dispersal but can also be attributed to marker-specific evolutionary characteristics and sampling properties. We discriminated between these competing explanations with a population genetic study on olive sea snakes, Aipysurus laevis. A previous mtDNA study revealed strong regional population structure for A. laevis around northern Australia, where Pleistocene sea-level fluctuations have influenced the genetic signatures of shallow-water marine species. Divergences among phylogroups dated to the Late Pleistocene, suggesting recent range expansions by previously isolated matrilines. Fine-scale population structure within regions was, however, poorly resolved for mtDNA. In order to improve estimates of fine-scale genetic divergence and to compare population structure between nuclear and mtDNA, 354 olive sea snakes (previously sequenced for mtDNA) were genotyped for five microsatellite loci. F statistics and Bayesian multilocus genotype clustering analyses found similar regional population structure as mtDNA and, after standardizing microsatellite F statistics for high heterozygosities, regional divergence estimates were quantitatively congruent between marker classes. Over small spatial scales, however, microsatellites recovered almost no genetic structure and standardized F statistics were orders of magnitude smaller than for mtDNA. Three tests for male-biased dispersal were not significant, suggesting that recent demographic expansions to the typically large population sizes of A. laevis have prevented microsatellites from reaching mutation-drift equilibrium and local populations may still be diverging.     

Molecular phylogeography and hybridization in members of the circumpolar Potentilla sect. Niveae (Rosaceae)

Glacial events and the formation of ice-free areas serving as refugia for plants and animals are important in shaping present patterns of genetic diversity in arctic areas. Beringia, situated in northeastern Russia and Alaska, has been pointed out as a major refugium. This study focuses on the historical biogeography of the circumpolar taxon Potentilla sect. Niveae. The taxonomy of the group is complex, most likely highly influenced by hybridization and apomixis. cpDNA microsatellites together with AFLP fragments were used to map the genetic variability in the section, from Beringia across the Canadian Arctic to Greenland. The data support the hypothesis that Beringia, as well as parts of adjacent arctic Canada, served as refugia during the Wisconsinan glaciation, and there is some evidence for a northern and a southern migration route out of Beringia. The hair type groups within sect. Niveae are more or less genetically distinct, and hybridization, especially with sect. Multifida, takes place. Haplotype diversity as well as frequency is at its maximum close to the Last Glacial Maximum ice cap edge. This pattern can be explained by merging of previously isolated refugia, by repeated extinction/colonization events close to the ice edge, and by hybridization among sympatric taxonomical lineages.     

Beringia and the global dispersal of modern humans

Until recently, the settlement of the Americas seemed largely divorced from the out‐of‐Africa dispersal of anatomically modern humans, which began at least 50,000 years ago. Native Americans were thought to represent a small subset of the Eurasian population that migrated to the Western Hemisphere less than 15,000 years ago. Archeological discoveries since 2000 reveal, however, that Homo sapiens occupied the high‐latitude region between Northeast Asia and northwest North America (that is, Beringia) before 30,000 years ago and the Last Glacial Maximum (LGM). The settlement of Beringia now appears to have been part of modern human dispersal in northern Eurasia. A 2007 model, the Beringian Standstill Hypothesis, which is based on analysis of mitochondrial DNA (mtDNA) in living people, derives Native Americans from a population that occupied Beringia during the LGM. The model suggests a parallel between ancestral Native Americans and modern human populations that retreated to refugia in other parts of the world during the arid LGM. It is supported by evidence of comparatively mild climates and rich biota in south‐central Beringia at this time (30,000‐15,000 years ago). These and other developments suggest that the settlement of the Americas may be integrated with the global dispersal of modern humans.     

Phylogeographic analyses of a circumarctic coastal and a boreal lacustrine mysid crustacean, and evidence of fast postglacial mtDNA rates

Phylogeographic structures of two weakly dispersing Mysis sibling species, one with a circumarctic coastal, the other with a boreal lacustrine-Baltic distribution, were studied from mitochondrial cytochrome c oxidase subunit I (COI) gene sequences. Mysis segerstralei showed high overall diversity and little phylogeographic structure across the Arctic, indicating late-glacial dispersal among coastal and lake populations from Alaska, Siberia and the north of Europe. A strongly divergent refugial lineage was however identified in Beringia. The boreal 'glacial relict'Mysis salemaai in turn displayed clear structuring among postglacially isolated Scandinavian lake populations. The inferred pattern of intralake mitochondrial DNA (mtDNA) monophyly in Scandinavia suggested relatively small population sizes and a remarkably fast postglacial mtDNA divergence rate (0.27% per 10 000 years). Nevertheless, the broader phylogeographic pattern did not support distinct eastern and western glacial refugia in Northern Europe, unlike in some other aquatic taxa. In all, the two species comprised three equidistant mitochondrial lineages (approximately 2% divergence), corresponding to M. salemaai, to the bulk of M. segerstralei, and to the Beringian M. segerstralei lineage. The lack of reciprocal monophyly of the two species in respect to their mitochondrial genealogy could indicate postspeciation mitochondrial introgression, also exemplified by an evidently more recent capture of M. segerstralei mitochondria in a Karelian population of M. salemaai. Overall, the data suggest that the continental boreal M. salemaai has a relatively recent ancestry in arctic coastal waters, whereas two other boreal 'glacial relict'Mysis sibling species in Europe (Mysis relicta) and North America (Mysis diluviana) have colonized inland waters much earlier (approximately 8% COI divergence).     

Effects of Pleistocene glaciations on population structure of North American chestnut-backed chickadees

The postglacial recolonization of northern North America was heavily influenced by the Pleistocene glaciation. In the Pacific Northwest, there are two disjunct regions of mesic temperate forest, one coastal and the other interior. The chestnut-backed chickadee is one of the species associated with this distinctive ecosystem. Using seven microsatellite markers we found evidence of population structure among nine populations of chestnut-backed chickadees. High levels of allelic variation were found in each of the populations. Northern British Columbia and central Alaska populations contained a large number of private alleles compared to other populations, including those from unglaciated regions. The disjunct population in the interior was genetically distinct from the coastal population. Genetic and historical records indicate that the interior population originated from postglacial inland dispersal. Population structuring was found within the continuous coastal population, among which the peripheral populations, specifically those on the Queen Charlotte Islands and the central Alaska mainland, were genetically distinct. The pattern of population structure among contemporary chickadee populations is consistent with a pioneer model of recolonization. The persistence of genetic structure in western North American chestnut-backed chickadees may be aided by their sedentary behaviour, linear distribution, and dependence on cedar-hemlock forests.     

DISCORDANT DISTRIBUTION OF POPULATIONS AND GENETIC VARIATION IN A SEA STAR WITH HIGH DISPERSAL POTENTIAL

Patiria miniata, a broadcast‐spawning sea star species with high dispersal potential, has a geographic range in the intertidal zone of the northeast Pacific Ocean from Alaska to California that is characterized by a large range gap in Washington and Oregon. We analyzed spatial genetic variation across the P. miniata range using multilocus sequence data (mtDNA, nuclear introns) and multilocus genotype data (microsatellites). We found a strong phylogeographic break at Queen Charlotte Sound in British Columbia that was not in the location predicted by the geographical distribution of the populations. However, this population genetic discontinuity does correspond to previously described phylogeographic breaks in other species. Northern populations from Alaska and Haida Gwaii were strongly differentiated from all southern populations from Vancouver Island and California. Populations from Vancouver Island and California were undifferentiated with evidence of high gene flow or very recent separation across the range disjunction between them. The surprising and discordant spatial distribution of populations and alleles suggests that historical vicariance (possibly caused by glaciations) and contemporary dispersal barriers (possibly caused by oceanographic conditions) both shape population genetic structure in this species.     

Phylogeographical lineages of Arctic grayling (Thymallus arcticus) in North America: divergence, origins and affinities with Eurasian Thymallus

The number and location of Arctic glacial refugia utilized by taxa during the Pleistocene are continuing uncertainties in Holarctic phylogeography. Arctic grayling (Thymallus arcticus) are widely distributed in freshwaters from the eastern side of Hudson Bay (Canada) west to central Asia. We studied mitochondrial DNA (mtDNA) and microsatellite DNA variation in North American T. arcticus to test for genetic signatures of survival in, and postglacial dispersal from, multiple glacial refugia, and to assess their evolutionary affinities with Eurasian Thymallus. In samples from 32 localities, we resolved 12 mtDNA haplotypes belonging to three assemblages that differed from each other in sequence by between 0.75 and 2.13%: a 'South Beringia' lineage found from western Alaska to northern British Columbia, Canada; a 'North Beringia' lineage found on the north slope of Alaska, the lower Mackenzie River, and to eastern Saskatchewan; and a 'Nahanni' lineage confined to the Nahanni River area of the upper Mackenzie River drainage. Sequence analysis of a portion of the control region indicated monophyly of all North American T. arcticus and their probable origin from eastern Siberian T. arcticus at least 3 Mya. Arctic grayling sampled from 25 localities displayed low allelic diversity and expected heterozygosity (H(E)) across five microsatellite loci (means of 2.1 alleles and 0.27 H(E), respectively) and there were declines in these measures of genetic diversity with distance eastward from the lower Yukon River Valley. Assemblages defined by mtDNA divergences were less apparent at microsatellite loci, but again the Nahanni lineage was the most distinctive. Analysis of molecular variance indicated that between 24% (microsatellite DNA) and 81% (mtDNA) of the variance was attributable to differences among South Beringia, North Beringia and Nahanni lineages. Our data suggest that extant North American Arctic grayling are more diverse phylogeographically than previously suspected and that they consist of at least three major lineages that originated in distinct Pleistocene glacial refugia. T. arcticus probably originated and dispersed from Eurasia to North America in the late to mid-Pliocene, but our data also suggest more recent (mid-late Pleistocene) interactions between lineages across Beringia.     

Reconstruction of caribou evolutionary history in Western North America and its implications for conservation

The role of Beringia as a refugium and route for trans-continental exchange of fauna during glacial cycles of the past 2million years are well documented; less apparent is its contribution as a significant reservoir of genetic diversity. Using mitochondrial DNA sequences and 14 microsatellite loci, we investigate the phylogeographic history of caribou (Rangifer tarandus) in western North America. Patterns of genetic diversity reveal two distinct groups of caribou. Caribou classified as a Northern group, of Beringian origin, exhibited greater number and variability in mtDNA haplotypes compared to a Southern group originating from refugia south of glacial ice. Results indicate that subspecies R. t. granti of Alaska and R. t. groenlandicus of northern Canada do not constitute distinguishable units at mtDNA or microsatellites, belying their current status as separate subspecies. Additionally, the Northern Mountain ecotype of woodland caribou (presently R. t. caribou) has closer kinship to caribou classified as granti or groenlandicus. Comparisons of mtDNA and microsatellite data suggest that behavioural and ecological specialization is a more recently derived life history characteristic. Notably, microsatellite differentiation among Southern herds is significantly greater, most likely as a result of human-induced landscape fragmentation and genetic drift due to smaller population sizes. These results not only provide important insight into the evolutionary history of northern species such as caribou, but also are important indicators for managers evaluating conservation measures for this threatened 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.     

The post‐Pleistocene population genetic structure of a western North American passerine: the chestnut‐backed chickadee Poecile rufescens

The population genetic structure of many high‐latitude species in North America was affected by the last glaciation, and current structure reflects isolation in refugia and colonisation patterns. Large ice‐free areas, both south of the ice sheets and in the north‐west, supported numerous flora and fauna throughout this period. Fossil evidence suggests additional western glacial refugia existed both on Haida Gwaii (the Queen Charlotte Islands) and in northern Idaho. The chestnut‐backed chickadee Poecile rufescens is a songbird found along the western edge of Canada and the United States, with a linear distribution along the coast, and an isolated interior population. Mitochondrial DNA sequence data (control region and ATPase 6–8) from 10 populations (n = 122) were used to test for population genetic structure. The data supported a general north/south separation. Haida Gwaii was found to be genetically distinct from the rest of the populations, and the two northern British Columbia populations separated from all but Alaska. The interior population showed evidence of both historical isolation and secondary colonisation by birds from coastal populations. Neutrality tests suggested a past population expansion in all populations from previously glaciated areas, and a stable population in areas believed to be unglaciated. This pattern supports the use of multiple glacial refugia by the chestnut‐backed chickadee. We could not reject the use of Haida Gwaii or the interior (i.e. Clearwater Basin) as glacial refugia.     

Molecular markers provide insights into contemporary and historic gene flow for a non‐migratory species

Hairy woodpeckers Picoides villosus are a common, year round resident with distinct plumage and morphological variation across North America. We genotyped 335 individuals at six variable microsatellite loci and analyzed 322 mtDNA control region sequences in order to examine the role of contemporary and historical barriers to gene flow. In addition we combined genetic analyses with ecological niche modelling to test if hairy woodpeckers were isolated in northern refugia (Alaska, Newfoundland and the Queen Charlotte Islands) during the last glacial maximum. Genetic analyses revealed that gene flow among North American hairy woodpecker populations is restricted, but not to the extent predicted for a sedentary species. Populations clustered into two main genetic groups, east and west of the Great Plains in the south and the Rocky Mountains in the north. Contact zones between the two main genetic groups exist in central British Columbia and Washington, but are narrow. Within each group we found additional population structure with genetic breaks between subgroups in the geographic west corresponding to breaks in forested habitat and physical barriers like open expanses of water. Population genetic patterns for hairy woodpeckers have resulted from isolation in multiple southern refugia with the current distribution of genetic groups resulting from post‐glacial expansion and subsequent reduction in gene flow. While populations in Alaska, Newfoundland and the Queen Charlotte Islands are genetically distinct from other populations, we found no evidence of these areas acting as refugia throughout the Pleistocene. Atlantic Canada populations contained unique haplotypes raising the possibility of a separate colonization from the rest of eastern Canada. The endemic subspecies on the island of Newfoundland is not genetically distinct from their closest mainland population unlike the Queen Charlotte Island subspecies.     

Less Pollen-Mediated Gene Flow for More Signatures of Glacial Lineages: Congruent Evidence from Balsam Fir cpDNA and mtDNA for Multiple Refugia in Eastern and Central North America

The phylogeographic structure and postglacial history of balsam fir (Abies balsamea), a transcontinental North American boreal conifer, was inferred using mitochondrial DNA (mtDNA) and chloroplast DNA (cpDNA) markers. Genetic structure among 107 populations (mtDNA data) and 75 populations (cpDNA data) was analyzed using Bayesian and genetic distance approaches. Population differentiation was high for mtDNA (dispersed by seeds only), but also for cpDNA (dispersed by seeds and pollen), indicating that pollen gene flow is more restricted in balsam fir than in other boreal conifers. Low cpDNA gene flow in balsam fir may relate to low pollen production due to the inherent biology of the species and populations being decimated by recurrent spruce budworm epidemics, and/or to low dispersal of pollen grains due to their peculiar structural properties. Accordingly, a phylogeographic structure was detected using both mtDNA and cpDNA markers and population structure analyses supported the existence of at least five genetically distinct glacial lineages in central and eastern North America. Four of these would originate from glacial refugia located south of the Laurentide ice sheet, while the last one would have persisted in the northern Labrador region. As expected due to reduced pollen-mediated gene flow, congruence between the geographic distribution of mtDNA and cpDNA lineages was higher than in other North American conifers. However, concordance was not complete, reflecting that restricted but nonetheless detectable cpDNA gene flow among glacial lineages occurred during the Holocene. As a result, new cpDNA and mtDNA genome combinations indicative of cytoplasmic genome capture were observed.     

Phylogeographic and population genetic analyses reveal Pleistocene isolation followed by high gene flow in a wide ranging,but endangered,freshwater mussel

Freshwater organisms of North America have had their contemporary genetic structure shaped by vicariant events, especially Pleistocene glaciations. Life history traits promoting dispersal and gene flow continue to shape population genetic structure. Cumberlandia monodonta, a widespread but imperiled (IUCN listed as endangered) freshwater mussel, was examined to determine genetic diversity and population genetic structure throughout its range. Mitochondrial DNA sequences and microsatellite loci were used to measure genetic diversity and simulate demographic events during the Pleistocene using approximate Bayesian computation (ABC) to test explicit hypotheses explaining the evolutionary history of current populations. A phylogeny and molecular clock suggested past isolation created two mtDNA lineages during the Pleistocene that are now widespread. Two distinct groups were also detected with microsatellites. ABC simulations indicated the presence of two glacial refugia and post-glacial admixture of them followed by simultaneous dispersal throughout the current range of the species. The Ouachita population is distinct from others and has the lowest genetic diversity, indicating that this is a peripheral population of the species. Gene flow within this species has maintained high levels of genetic diversity in most populations; however, all populations have experienced fragmentation. Extirpation from the center of its range likely has isolated remaining populations due to the geographic distances among them.     

Of glaciers and refugia: a decade of study sheds new light on the phylogeography of northwestern North America

Glacial cycles have played a dominant role in shaping the genetic structure and distribution of biota in northwestern North America. The two major ice age refugia of Beringia and the Pacific Northwest were connected by major mountain chains and bordered by the Pacific Ocean. As a result, numerous refugial options were available for the regions taxa during glacial advances. We reviewed the importance of glaciations and refugia in shaping northwestern North America’s phylogeographic history. We also tested whether ecological variables were associated with refugial history. The recurrent phylogeographic patterns that emerged were the following: (i) additional complexity, i.e. refugia within refugia, in both Beringia and the Pacific Northwest; and (ii) strong evidence for cryptic refugia in the Alexander Archipelago and Haida Gwaii, the Canadian Arctic and within the ice‐sheets. Species with contemporary ranges that covered multiple refugia, or those with high dispersal ability, were significantly more likely to have resided in multiple refugia. Most of the shared phylogeographic patterns can be attributed to multiple refugial locales during the last glacial maximum or major physiographic barriers like rivers and glaciers. However, some of the observed patterns are much older and appear connected to the orogeny of the Cascade‐Sierra chain or allopatric differentiation during historic glacial advances. The emergent patterns from this review suggest we should refine the classic Beringian‐southern refugial paradigm for northwestern North American biota and highlight the ecological and evolutionary consequences of colonization from multiple refugia.     

Pleistocene glaciations and contemporary genetic diversity in a Beringian fish,the broad whitefish,Coregonus nasus (Pallas): inferences from microsatellite DNA variation

The contemporary distribution of genetic variation within and among high latitude populations cannot be fully understood without taking into consideration how species responded to the impacts of Pleistocene glaciations. Broad whitefish, Coregonus nasus, a species endemic to northwest North America and the Arctic coast of Russia, was undoubtedly impacted by such events because its geographic distribution suggests that it survived solely within the Beringian refuge from where it dispersed post‐glacially to achieve its current range. We used microsatellite DNA to investigate the role of glaciations in promoting intraspecific genetic variation in broad whitefish (N = 14 localities, 664 fish) throughout their North American range and in one Russian sample. Broad whitefish exhibited relatively high intrapopulation variation (average of 11.7 alleles per locus, average H_E = 0.61) and moderate levels of interpopulation divergence (overall F_ST = 0.10). The main regions assayed in our study (Russia, Alaska, Mackenzie River and Travaillant Lake systems) were genetically differentiated from each other and there were declines in genetic diversity with distance from putative refugia. Additionally, Mackenzie River system populations showed less developed and more variable patterns of isolation‐by‐distance than populations occupying former Alaskan portions of Beringia. Finally, our data suggest that broad whitefish dispersed from Beringia using coastal environments and opportunistically via headwater stream connections that once existed between Yukon and Mackenzie River drainages. Our results illustrate the importance of history (e.g. glaciation) and contemporary dispersal ecology in shaping the current genetic population structure of Arctic faunas.     

Multilocus genetic analyses and spatial modeling reveal complex population structure and history in a widespread resident North American passerine (Perisoreus canadensis)

An increasing body of studies of widely distributed, high latitude species shows a variety of refugial locations and population genetic patterns. We examined the effects of glaciations and dispersal barriers on the population genetic patterns of a widely distributed, high latitude, resident corvid, the gray jay (Perisoreus canadensis), using the highly variable mitochondrial DNA (mtDNA) control region and microsatellite markers combined with species distribution modeling. We sequenced 914 bp of mtDNA control region for 375 individuals from 37 populations and screened seven loci for 402 individuals from 27 populations across the gray jay range. We used species distribution modeling and a range of phylogeographic analyses (haplotype diversity, Φ_ST, SAMOVA, F_ST, Bayesian clustering analyses) to examine evolutionary history and population genetic structure. MtDNA and microsatellite markers revealed significant genetic differentiation among populations with high concordance between markers. Paleodistribution models supported at least five potential areas of suitable gray jay habitat during the last glacial maximum and revealed distributions similar to the gray jay's contemporary during the last interglacial. Colonization from and prolonged isolation in multiple refugia is evident. Historical climatic fluctuations, the presence of multiple dispersal barriers, and highly restricted gene flow appear to be responsible for strong genetic diversification and differentiation in gray jays.     

Genetic structure in the European endemic seabird,Phalacrocorax aristotelis,shaped by a complex interaction of historical and contemporary,physical and nonphysical drivers

Geographically separated populations tend to be less connected by gene flow, as a result of physical or nonphysical barriers preventing dispersal, and this can lead to genetic structure. In this context, highly mobile organisms such as seabirds are interesting because the small effect of physical barriers means nonphysical ones may be relatively more important. Here, we use microsatellite and mitochondrial data to explore the genetic structure and phylogeography of Atlantic and Mediterranean populations of a European endemic seabird, the European shag, Phalacrocorax aristotelis, and identify the primary drivers of their diversification. Analyses of mitochondrial markers revealed three phylogenetic lineages grouping the North Atlantic, Spanish/Corsican and eastern Mediterranean populations, apparently arising from fragmentation during the Pleistocene followed by range expansion. These traces of historical fragmentation were also evident in the genetic structure estimated by microsatellite markers, despite significant contemporary gene flow among adjacent populations. Stronger genetic structure, probably promoted by landscape, philopatry and local adaptation, was found among distant populations and those separated by physical and ecological barriers. This study highlights the enduring effect of Pleistocene climatic changes on shag populations, especially within the Mediterranean Basin, and suggests a role for cryptic northern refugia, as well as known southern refugia, on the genetic structure of European seabirds. Finally, it outlines how contemporary ecological barriers and behavioural traits may maintain population divergence, despite long‐distance dispersal triggered by extreme environmental conditions (e.g. population crashes).     

Phylogeography of northern Dolly Varden Salvelinus malma (Salmoniformes: Salmonidae) from Asia and North America: An analysis based on the mitochondrial DNA genealogy

Spawning in habitats affected by Pleistocene glacial advances over most of its natural range, northern Dolly Varden Salvelinus malma malma typifies Arctic fauna distributed in northeastern Asia and northwestern North America. We reconstructed a genealogy of mtDNA haplotypes from 27 Alaskan and Asian populations to study the influence of historical events on the phylogeography and contemporary population genetic structure. Analysis of molecular variance partitioned most of the mtDNA variability to the intrapopulation component (72.5%) with much reduced differences between populations (21.1%) and regions (6.4%). Similar patterns of variation apparent from hierarchical diversity and nested clade phylogeographical analysis (NCPA) of mtDNA haplotypes identify weak spatial differentiation and low levels of divergence. These findings suggest (1) that demographic history has been influenced by historical range expansions and recent isolation by distance, (2) that present populations from Asia and North America were colonized from one main Beringian Refugium, and (3) that this taxon’s ancestral population probably experienced a bottleneck in the Beringian Refugium during the late Pleistocene (Wisconsin) glacial period. The genealogical and NCPA analyses, and mismatch distribution of S. m. malma mtDNA haplotypes do not confirm the assumptions about presence of the two refugia on the territories of the Beringian Land, in which allopatric S. m. malma ancestral populations evolved, and independent origin of the Sea of Okhotsk populations.     

Population structure of two rabies hosts relative to the known distribution of rabies virus variants in Alaska

For pathogens that infect multiple species, the distinction between reservoir hosts and spillover hosts is often difficult. In Alaska, three variants of the arctic rabies virus exist with distinct spatial distributions. We tested the hypothesis that rabies virus variant distribution corresponds to the population structure of the primary rabies hosts in Alaska, arctic foxes (Vulpes lagopus) and red foxes (Vulpes vulpes) to possibly distinguish reservoir and spillover hosts. We used mitochondrial DNA (mtDNA) sequence and nine microsatellites to assess population structure in those two species. mtDNA structure did not correspond to rabies virus variant structure in either species. Microsatellite analyses gave varying results. Bayesian clustering found two groups of arctic foxes in the coastal tundra region, but for red foxes it identified tundra and boreal types. Spatial Bayesian clustering and spatial principal components analysis identified 3 and 4 groups of arctic foxes, respectively, closely matching the distribution of rabies virus variants in the state. Red foxes, conversely, showed eight clusters comprising two regions (boreal and tundra) with much admixture. These results run contrary to previous beliefs that arctic fox show no fine‐scale spatial population structure. While we cannot rule out that the red fox is part of the maintenance host community for rabies in Alaska, the distribution of virus variants appears to be driven primarily by the arctic fox. Therefore, we show that host population genetics can be utilized to distinguish between maintenance and spillover hosts when used in conjunction with other approaches.