Reconstruction of a beech population bottleneck using archival demographic information and Bayesian analysis of genetic data

Range expansion and contraction has occurred in the history of most species and can seriously impact patterns of genetic diversity. Historical data about range change are rare and generally appropriate for studies at large scales, whereas the individual pollen and seed dispersal events that form the basis of geneflow and colonization generally occur at a local scale. In this study, we investigated range change in Fagus sylvatica on Mont Ventoux, France, using historical data from 1838 to the present and approximate Bayesian computation (ABC) analyses of genetic data. From the historical data, we identified a population minimum in 1845 and located remnant populations at least 200 years old. The ABC analysis selected a demographic scenario with three populations, corresponding to two remnant populations and one area of recent expansion. It also identified expansion from a smaller ancestral population but did not find that this expansion followed a population bottleneck, as suggested by the historical data. Despite a strong support to the selected scenario for our data set, the ABC approach showed a low power to discriminate among scenarios on average and a low ability to accurately estimate effective population sizes and divergence dates, probably due to the temporal scale of the study. This study provides an unusual opportunity to test ABC analysis in a system with a well-documented demographic history and identify discrepancies between the results of historical, classical population genetic and ABC analyses. The results also provide valuable insights into genetic processes at work at a fine spatial and temporal scale in range change and colonization.     

Intraspecific range dynamics and niche evolution in Candidula land snail species

The range dynamics of a species can either be governed by the spatial tracing of the fundamental environmental niche or by adaptation that allows to occupy new niches. Therefore, the investigation of spatial variation in the realized environmental niche is central to the understanding of species range limit dynamics. However, the study of intraspecific niche variation has been neglected in most phylogeographical studies. We studied the spatial distribution of the realized environmental niche in three land snail species of the genus Candidula , integrating phylogeographical methods, morphometrics, and spatial biodiversity informatics . The phylogeographical analyses showed significant range expansions in all species. These expansions were accompanied in Candidula gigaxii by a shift in the realized environmental niche, the species Candidula unifasciata followed its ancestral niche during expansion while the climate changed in the area of origin and Candidula rugosiuscula tracked the ancestral environmental conditions. The significant niche shifts were associated with potentially adaptive changes of shell morphology. We propose our presented approach as a practicable framework to test hypotheses on intraspecific niche evolution. © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 90 , 303–317.     

The demographic history of Atlantic salmon (Salmo salar) across its distribution range reconstructed from approximate Bayesian computations*

Understanding the dual roles of demographic and selective processes in the buildup of population divergence is one of the most challenging tasks in evolutionary biology. Here, we investigated the demographic history of Atlantic salmon across the entire species range using 2035 anadromous individuals from North America and Eurasia. By combining results from admixture graphs, geo‐genetic maps, and an Approximate Bayesian Computation (ABC) framework, we validated previous hypotheses pertaining to secondary contact between European and Northern American populations, but also identified secondary contacts in European populations from different glacial refugia. We further identified the major sources of admixture from the southern range of North America into more northern populations along with a strong signal of secondary gene flow between genetic regional groups. We hypothesize that these patterns reflect the spatial redistribution of ancestral variation across the entire North American range. Results also support a role for linked selection and differential introgression that likely played an underappreciated role in shaping the genomic landscape of species in the Northern hemisphere. We conclude that studies between partially isolated populations should systematically include heterogeneity in selective and introgressive effects among loci to perform more rigorous demographic inferences of the divergence process.     

Rapid climate change and the rate of adaptation: insight from experimental quantitative genetics

CONTENTS: Summary 752 I. Introduction 752 II. Will migration be enough? 753 III. Can adaptation proceed fast enough? 754 IV. Fitness links demographic and evolutionary processes 755 V. Experimental studies: what do they tell us and how can we improve them? 756 VI. Predicting evolutionary change based on genetic variation and natural selection 757 VII. The chronosequence approach 758 VIII. Resurrection of ancestral propagules 759 IX. The mean and variance in fitness, a link between genetics and demography 760 X. Conclusions 762 Acknowledgements 762 References 762 SUMMARY: Evolution proceeds unceasingly in all biological populations. It is clear that climate-driven evolution has molded plants in deep time and within extant populations. However, it is less certain whether adaptive evolution can proceed sufficiently rapidly to maintain the fitness and demographic stability of populations subjected to exceptionally rapid contemporary climate change. Here, we consider this question, drawing on current evidence on the rate of plant range shifts and the potential for an adaptive evolutionary response. We emphasize advances in understanding based on theoretical studies that model interacting evolutionary processes, and we provide an overview of quantitative genetic approaches that can parameterize these models to provide more meaningful predictions of the dynamic interplay between genetics, demography and evolution. We outline further research that can clarify both the adaptive potential of plant populations as climate continues to change and the role played by ongoing adaptation in their persistence.     

Genetic surfing,not allopatric divergence,explains spatial sorting of mitochondrial haplotypes in venomous coralsnakes

Strong spatial sorting of genetic variation in contiguous populations is often explained by local adaptation or secondary contact following allopatric divergence. A third explanation, spatial sorting by stochastic effects of range expansion, has been considered less often though theoretical models suggest it should be widespread, if ephemeral. In a study designed to delimit species within a clade of venomous coralsnakes, we identified an unusual pattern within the Texas coral snake (Micrurus tener): strong spatial sorting of divergent mitochondrial (mtDNA) lineages over a portion of its range, but weak sorting of these lineages elsewhere. We tested three alternative hypotheses to explain this pattern—local adaptation, secondary contact following allopatric divergence, and range expansion. Collectively, near panmixia of nuclear DNA, the signal of range expansion associated sampling drift, expansion origins in the Gulf Coast of Mexico, and species distribution modeling suggest that the spatial sorting of divergent mtDNA lineages within M. tener has resulted from genetic surfing of standing mtDNA variation—not local adaptation or allopatric divergence. Our findings highlight the potential for the stochastic effects of recent range expansion to mislead estimations of population divergence made from mtDNA, which may be exacerbated in systems with low vagility, ancestral mtDNA polymorphism, and male‐biased dispersal.     

Population growth of human Y chromosomes: a study of Y chromosome microsatellites

We use variation at a set of eight human Y chromosome microsatellite loci to investigate the demographic history of the Y chromosome. Instead of assuming a population of constant size, as in most of the previous work on the Y chromosome, we consider a model which permits a period of recent population growth. We show that for most of the populations in our sample this model fits the data far better than a model with no growth. We estimate the demographic parameters of this model for each population and also the time to the most recent common ancestor. Since there is some uncertainty about the details of the microsatellite mutation process, we consider several plausible mutation schemes and estimate the variance in mutation size simultaneously with the demographic parameters of interest. Our finding of a recent common ancestor (probably in the last 120,000 years), coupled with a strong signal of demographic expansion in all populations, suggests either a recent human expansion from a small ancestral population, or natural selection acting on the Y chromosome.     

Contrasting phylogeographical patterns for springtails reflect different evolutionary histories between the Antarctic Peninsula and continental Antarctica

Aim  We examined the genetic structure among populations and regions for the springtails Cryptopygus antarcticus antarcticus and Gomphiocephalus hodgsoni (Collembola) to identify potential historical refugia and subsequent colonization routes, and to examine population growth/expansion and relative ages of population divergence.
Location  Antarctic Peninsula for C. a. antarcticus ; Antarctic continent (southern Victoria Land) for G. hodgsoni .
Methods  Samples were collected from 24 and 28 locations across the Antarctic Peninsula and southern Victoria Land regions for C. a. antarcticus and G. hodgsoni , respectively. We used population genetic, demographic and nested clade analyses based on mitochondrial DNA (cytochrome c oxidase subunit I and subunit II).
Results  Both species were found to have population structures compatible with the presence of historical glacial refugia on Pleistocene (2 Ma–present) time-scales, followed by post-glacial expansion generating contemporary geographically isolated populations. However, G. hodgsoni populations were characterized by a fragmented pattern with several 'phylogroups' (likely ancestral haplotypes present in high frequency) retaining strong ancestral linkages among present-day populations. Conversely, C. a. antarcticus had an excess of rare haplotypes with a much reduced volume of ancestral lineages, possibly indicating historical founder/bottleneck events and widespread expansion.
Main conclusions  We infer that these differences reflect distinct evolutionary histories in each locality despite the resident species having similar life-history characteristics. We suggest that this has predominantly been influenced by variation in the success of colonization events as a result of intrinsic historical glaciological differences between the Antarctic Peninsula and continental Antarctic environments.     

The genetic structure of Drosophila ananassae populations from Asia, Australia and Samoa

Information about genetic structure and historical demography of natural populations is central to understanding how natural selection changes genomes. Drosophila ananassae is a widespread species occurring in geographically isolated or partially isolated populations and provides a unique opportunity to investigate population structure and molecular variation. We assayed microsatellite repeat-length variation among 13 populations of D. ananassae to assess the level of structure among the populations and to make inferences about their ancestry and historic biogeography. High levels of genetic structure are apparent among all populations, particularly in Australasia and the South Pacific, and patterns are consistent with the hypothesis that the ancestral populations are from Southeast Asia. Analysis of population structure and use of F-statistics and Bayesian analysis suggest that the range expansion of the species into the Pacific is complex, with multiple colonization events evident in some populations represented by lineages that show no evidence of recent admixture. The demographic patterns show isolation by distance among populations and population expansion within all populations. A morphologically distinct sister species, D. pallidosa, collected in Malololelei, Samoa, appears to be more closely related to some of the D. ananassae populations than many of the D. ananassae populations are to one another. The patterns of genotypic diversity suggest that many of the individuals that we sampled may be morphologically indistinguishable nascent species.     

Spatial and temporal patterns of neutral and adaptive genetic variation in the endangered African wild dog (Lycaon pictus)

Deciphering patterns of genetic variation within a species is essential for understanding population structure, local adaptation and differences in diversity between populations. Whilst neutrally evolving genetic markers can be used to elucidate demographic processes and genetic structure, they are not subject to selection and therefore are not informative about patterns of adaptive variation. As such, assessments of pertinent adaptive loci, such as the immunity genes of the major histocompatibility complex (MHC), are increasingly being incorporated into genetic studies. In this study, we combined neutral (microsatellite, mtDNA) and adaptive (MHC class II DLA‐DRB1 locus) markers to elucidate the factors influencing patterns of genetic variation in the African wild dog (Lycaon pictus); an endangered canid that has suffered extensive declines in distribution and abundance. Our genetic analyses found all extant wild dog populations to be relatively small (N_e < 30). Furthermore, through coalescent modelling, we detected a genetic signature of a recent and substantial demographic decline, which correlates with human expansion, but contrasts with findings in some other African mammals. We found strong structuring of wild dog populations, indicating the negative influence of extensive habitat fragmentation and loss of gene flow between habitat patches. Across populations, we found that the spatial and temporal structure of microsatellite diversity and MHC diversity were correlated and strongly influenced by demographic stability and population size, indicating the effects of genetic drift in these small populations. Despite this correlation, we detected signatures of selection at the MHC, implying that selection has not been completely overwhelmed by genetic drift.     

Region‐wide and ecotype‐specific differences in demographic histories of threespine stickleback populations,estimated from whole genome sequences

We analysed 81 whole genome sequences of threespine sticklebacks from Pacific North America, Greenland and Northern Europe, representing 16 populations. Principal component analysis of nuclear SNPs grouped populations according to geographical location, with Pacific populations being more divergent from each other relative to European and Greenlandic populations. Analysis of mitogenome sequences showed Northern European populations to represent a single phylogeographical lineage, whereas Greenlandic and particularly Pacific populations showed admixture between lineages. We estimated demographic history using a genomewide coalescence with recombination approach. The Pacific populations showed gradual population expansion starting >100 Kya, possibly reflecting persistence in cryptic refuges near the present distributional range, although we do not rule out possible influence of ancient admixture. Sharp population declines ca. 14–15 Kya were suggested to reflect founding of freshwater populations by marine ancestors. In Greenland and Northern Europe, demographic expansion started ca. 20–25 Kya coinciding with the end of the Last Glacial Maximum. In both regions, marine and freshwater populations started to show different demographic trajectories ca. 8–9 Kya, suggesting that this was the time of recolonization. In Northern Europe, this estimate was surprisingly late, but found support in subfossil evidence for presence of several freshwater fish species but not sticklebacks 12 Kya. The results demonstrate distinctly different demographic histories across geographical regions with potential consequences for adaptive processes. They also provide empirical support for previous assumptions about freshwater populations being founded independently from large, coherent marine populations, a key element in the Transporter Hypothesis invoked to explain the widespread occurrence of parallel evolution across freshwater stickleback populations.     

The phylogeographical clade trade: tracing the impact of human-mediated dispersal on the colonization of northern Europe by the oak gallwasp Andricus kollari

Human dispersal of organisms is an important process modifying natural patterns of biodiversity. Such dispersal generates new patterns of genetic diversity that overlie natural phylogeographical signatures, allowing discrimination between alternative dispersal mechanisms. Here we use allele frequency and DNA sequence data to distinguish between alternative scenarios (unassisted range expansion and long range introduction) for the colonization of northern Europe by an oak-feeding gallwasp, Andricus kollari. Native to Mediterranean latitudes from Portugal to Iran, this species became established in northern Europe following human introduction of a host plant, the Turkey oak Quercus cerris. Colonization of northern Europe is possible through three alternative routes: (i) unassisted range expansion from natural populations in the Iberian Peninsula; (ii) unassisted range expansion from natural populations in Italy and Hungary; or (iii) descent from populations imported to the UK as trade goods from the eastern Mediterranean in the 1830s. We show that while populations in France were colonized from sources in Italy and Hungary, populations in the UK and neighbouring parts of coastal northern Europe encompass allozyme and sequence variation absent from the known native range. Further, these populations show demographic signatures expected for large stable populations, rather than signatures of rapid population growth from small numbers of founders. The extent and spatial distribution of genetic diversity in the UK suggests that these A. kollari populations are derived from introductions of large numbers of individuals from each of two genetically divergent centres of diversity in the eastern Mediterranean. The strong spatial patterning in genetic diversity observed between different regions of northern Europe, and between sites in the UK, is compatible with leptokurtic models of population establishment.     

Population Genomics of Sub-Saharan Drosophila melanogaster: African Diversity and Non-African Admixture

Drosophila melanogaster has played a pivotal role in the development of modern population genetics. However, many basic questions regarding the demographic and adaptive history of this species remain unresolved. We report the genome sequencing of 139 wild-derived strains of D. melanogaster, representing 22 population samples from the sub-Saharan ancestral range of this species, along with one European population. Most genomes were sequenced above 25X depth from haploid embryos. Results indicated a pervasive influence of non-African admixture in many African populations, motivating the development and application of a novel admixture detection method. Admixture proportions varied among populations, with greater admixture in urban locations. Admixture levels also varied across the genome, with localized peaks and valleys suggestive of a non-neutral introgression process. Genomes from the same location differed starkly in ancestry, suggesting that isolation mechanisms may exist within African populations. After removing putatively admixed genomic segments, the greatest genetic diversity was observed in southern Africa (e.g. Zambia), while diversity in other populations was largely consistent with a geographic expansion from this potentially ancestral region. The European population showed different levels of diversity reduction on each chromosome arm, and some African populations displayed chromosome arm-specific diversity reductions. Inversions in the European sample were associated with strong elevations in diversity across chromosome arms. Genomic scans were conducted to identify loci that may represent targets of positive selection within an African population, between African populations, and between European and African populations. A disproportionate number of candidate selective sweep regions were located near genes with varied roles in gene regulation. Outliers for Europe-Africa F_ST were found to be enriched in genomic regions of locally elevated cosmopolitan admixture, possibly reflecting a role for some of these loci in driving the introgression of non-African alleles into African populations.     

Contrasted polymorphism patterns in a large sample of populations from the evolutionary genetics model Drosophila simulans

African populations of Drosophila simulans are thought to be ancestral in this model species and are increasingly used for testing general hypotheses in evolutionary genetics. It is often assumed that African populations are more likely to be at a neutral mutation drift equilibrium than other populations. Here we examine population structuring and the demographic profile in nine populations of D. simulans. We surveyed sequence variation in four X-linked genes (runt, sevenless, Sex-lethal, and vermilion) that have been used in a parallel study in the closely related species D. melanogaster. We found that an eastern group of populations from continental Africa and Indian Ocean islands (Kenya, Tanzania, Madagascar, and Mayotte Island) is widespread, shows little differentiation, and has probably undergone demographic expansion. The other two African populations surveyed (Cameroon and Zimbabwe) show no evidence of population expansion and are markedly differentiated from each other as well as from the populations from the eastern group. Two other populations, Europe and Antilles, are probably recent invaders to these areas. The Antilles population is probably derived from Europe through a substantial bottleneck. The history of these populations should be taken into account when drawing general conclusions from variation patterns.     

Rapid increase in dispersal during range expansion in the invasive ladybird Harmonia axyridis

The evolutionary trajectories associated with demographic, genetic and spatial disequilibrium have become an issue of growing interest in population biology. Invasive species provide unique opportunities to explore the impact of recent range expansion on life‐history traits, making it possible to test for a spatial arrangement of dispersal abilities along the expanding range, in particular. We carried out controlled experiments in laboratory conditions to test the hypothesis of an increase in dispersal capacity with range expansion in Harmonia axyridis, a ladybird that has been invading Europe since 2001. We found a marked increase in the flight speed of the insects from the core to the front of the invasion range in two independent sampling transects. By contrast, we found that two other traits associated with dispersal (endurance and motivation to fly off) did not follow the same spatial gradient. Our results provide a striking illustration of the way in which predictable directional genetic changes may occur rapidly for some traits associated with dispersal during biological invasions. We discuss the consequences of our results for invasion dynamics and the evolutionary outcomes of spatially expanding populations.     

Adaptive responses to cool climate promotes persistence of a non-native lizard

Successful establishment and range expansion of non-native species often require rapid accommodation of novel environments. Here, we use common-garden experiments to demonstrate parallel adaptive evolutionary response to a cool climate in populations of wall lizards (Podarcis muralis) introduced from southern Europe into England. Low soil temperatures in the introduced range delay hatching, which generates directional selection for a shorter incubation period. Non-native lizards from two separate lineages have responded to this selection by retaining their embryos for longer before oviposition—hence reducing the time needed to complete embryogenesis in the nest—and by an increased developmental rate at low temperatures. This divergence mirrors local adaptation across latitudes and altitudes within widely distributed species and suggests that evolutionary responses to climate can be very rapid. When extrapolated to soil temperatures encountered in nests within the introduced range, embryo retention and faster developmental rate result in one to several weeks earlier emergence compared with the ancestral state. We show that this difference translates into substantial survival benefits for offspring. This should promote short- and long-term persistence of non-native populations, and ultimately enable expansion into areas that would be unattainable with incubation duration representative of the native range.     

The role of variation and plasticity in parental care during the adaptive radiation of three‐spine sticklebacks

Phenotypic plasticity might influence evolutionary processes such as adaptive radiations. Plasticity in parental care might be especially effective in facilitating adaptive radiations if it allows populations to persist in novel environments. Here, we test the hypothesis that behavioral plasticity by parents in response to predation risk facilitated the adaptive radiation of three‐spine sticklebacks. We compared the behavior of fathers across multiple ancestral (marine) and derived (freshwater) stickleback populations that differ in time since establishment. We measured behavioral plasticity in fathers in response to a predator found only in freshwater environments, simulating conditions marine males experience when colonizing freshwater. The antipredator behavior of males from newly established freshwater populations was intermediate between marine populations and well‐established freshwater populations. In contrast to our predictions, on average, there was greater behavioral plasticity in derived freshwater populations than in ancestral marine populations. However, we found greater individual variation in behavioral reaction norms in marine populations compared to well‐established freshwater populations, with newly established freshwater populations intermediate. This suggests that standing variation in behavioral reaction norms within ancestral populations might provide different evolutionary trajectories, and illustrates how plasticity can contribute to adaptive radiations.     

Spatial pattern of chloroplast DNA variation of Cyclobalanopsis glauca in Taiwan and East Asia

This study examined the spatial pattern of chloroplast DNA (cpDNA) variation in Cyclobalanopsis glauca (Thunb. ex Murray) Oerst. (Fagaceae) in 140 trees from Taiwan (25 populations), Japan (three), Ryukyus (two), Hong Kong (one) and Mainland China (one). By sequencing three cpDNA intergenic spacer fragments using universal primers (trnT-trnL, trnV-trnM, including the trnV intron, and petG-trnP), we found a total of 1,980 bp and 15 polymorphic sites. Among them, 12 sites were caused by point mutation, and three resulted from insertion. This gives rise to a total of 13 cpDNA haplotypes. The level of differentiation among the populations studied is relatively high (GST = 0.612). Two ancestral haplotypes (A and B) are distributed widely in East Asia. Interestingly, all the derived cpDNA variations are found only in Taiwan but not in other areas. The Central Mountain Ridge (CMR) of Taiwan creates an unsurpassed barrier to the east-west gene flow of C. glauca. Among the populations on the west of CMR, only three separated populations, Yangmingshan, Wushe and Chinshuiying, have high haplotype diversity, each consisting of sister haplotypes all mutated from the same ancestral haplotype. Thus, they have probably originated from de novo mutation after the last glaciation. This inference agrees with the observation that no spatial autocorrelation existed on the west side. Two unrelated dominant lineages on the east of the CMR (haplotypes D and F) showed significant spatial genetic structure. Estimate of NST - GST was -0.090 and differed significantly from zero. Thus at the local scale, the phylogeographical component of the genetic structure is significant on the east of the CMR. Accompanied by published palynological records of the last glaciation, this study suggests the possibility that these two types were colonized northward from the southeastern part of Taiwan. 'Star-like' genealogy is characterized, with all the haplotypes coalescing rapidly and as a general outcome of population expansion (Page & Holmes 1998). A neutrality test also suggested a demographic expansion recovered from a bottleneck. We therefore inferred that the southeastern part of Taiwan might be a potential refugium for C. glauca.     

Causes and conditions of formation of the aromorphic organization

Analysis of data concerning the origin of terrestrial vertebrates allows to conclude that the organism undergone the aromorphosis (=arogenesis) not only improves the general degree of organization, but also becomes able to survive in much more complex and diversified environment than its ancestors and, consequently, expands the ancestor's adaptive zone. It was shown that aromorphosis occurs when the ancestral taxon specializes in a narrow adaptive zone having vacant neighboring zone which might be occupied. The aromorphic adaptations are formed very slowly in the course of so-called coherent evolution by accumulation of certain multiply adaptations and by their coordination with each other. In the end, the set of these adaptations overrides the oscillation range of environmental conditions in the ancestral adaptive zone and can be considered as a pre-adaptation to adaptive zone expansion. Some features which are characteristic for the modem 'aromorphic' groups such as Tetrapoda, Avia and Mammalia, appeared independently and asynchronously in many related taxa. However, appearance of a few detached features does not help to widen an ancestral adaptive zone. Taxon becomes really aromorphic only due to the simultaneous formation of all adaptations allowing the expansion of ancestral adaptive zone, i.e. successful exploration and occupation of neighborning adaptive zone.     

Genome-wide genetic variation coupled with demographic and ecological niche modeling of the dusky-footed woodrat (Neotoma fuscipes) reveal patterns of deep divergence and widespread Holocene expansion across northern California

Understanding how species have responded to past climate change may help refine projections of how species and biotic communities will respond to future change. Here, we integrate estimates of genome-wide genetic variation with demographic and niche modeling to investigate the historical biogeography of an important ecological engineer: the dusky-footed woodrat, Neotoma fuscipes. We use RADseq to generate a genome-wide dataset for 71 individuals from across the geographic distribution of the species in California. We estimate population structure using several model-based methods and infer the demographic history of regional populations using a site frequency spectrum-based approach. Additionally, we use ecological niche modeling to infer current and past (Last Glacial Maximum) environmental suitability across the species’ distribution. Finally, we estimate the directionality and possible spatial origins of regional population expansions. Our analyses indicate this species is subdivided into three regionally distinct populations, with the deepest divergence occurring ~1.7 million years ago across the modern-day San Francisco-Bay Delta region; a common biogeographic barrier for the flora and fauna of California. Our models of environmental suitability through time coincide with our estimates of population expansion, with relative long-term stability in the southern portion of the range, and more recent expansion into the northern end of the range. Our study illustrates how the integration of genome-wide data with spatial and demographic modeling can reveal the timing and spatial extent of historic events that determine patterns of biotic diversity and may help predict biotic response to future change.Subject terms: Evolutionary genetics, Molecular evolution     

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.