Isolation and introgression in the Intermountain West: contrasting gene genealogies reveal the complex biogeographic history of the American pika (Ochotona princeps)

Aim We studied the history of colonization, diversification and introgression among major phylogroups in the American pika, Ochotona princeps (Lagomorpha), using comparative and statistical phylogeographic methods. Our goal was to understand how Pleistocene climatic fluctuations have shaped the distribution of diversity at mitochondrial DNA (mtDNA) and nuclear DNA (nDNA) loci in this alpine specialist. Location North America’s Intermountain West. Methods We accumulated mtDNA sequence data (c. 560–1700 bp) from 232 pikas representing 64 localities, and sequenced two nuclear introns (mast cell growth factor, c. 550 bp, n = 148; protein kinase C iota, c. 660 bp, n = 139) from a subset of individuals. To determine the distribution of major mtDNA lineages, we conducted a phylogenetic analysis on the mtDNA sequence data, and we calculated divergence times among the lineages using a Bayesian Markov chain Monte Carlo approach. Relationships among nuclear alleles were explored with minimum spanning networks. Finally, we conducted coalescent simulations of alternative models of population history to test for congruence between nDNA and mtDNA responses to Pleistocene glacial cycles. Results We found that: (1) all individuals could be assigned to one of five allopatric mtDNA lineages; (2) lineages are associated with separate mountain provinces; (3) lineages originated from at least two rounds of differentiation; (4) nDNA and mtDNA markers exhibited overall phylogeographic congruence; and (5) introgression among phylogroups has occurred at nuclear loci since their initial isolation. Main conclusions Pika populations associated with different mountain systems have followed separate but not completely independent evolutionary trajectories through multiple glacial cycles. Range expansion associated with climate cooling (i.e. glaciations) promoted genetic admixture among populations within mountain ranges. It also permitted periodic contact and introgression between phylogroups associated with different mountain systems, the record of which is retained at nDNA but not mtDNA loci. Evidence for different histories at nuclear and mtDNA loci (i.e. periodic introgression versus deep isolation, respectively) emphasizes the importance of multilocus perspectives for reconstructing complete population histories.     

Phylogeography and population history of the crab-eating fox (Cerdocyon thous)

The crab-eating fox is a medium-sized Neotropical canid with generalist habits and a broad distribution in South America. We have investigated its genetic diversity, population structure and demographic history across most of its geographic range by analysing 512 base pairs (bp) of the mitochondrial DNA (mtDNA) control region, 615 bp of the mtDNA cytochrome b gene and 1573 total nucleotides from three different nuclear fragments. MtDNA data revealed a strong phylogeographic partition between northeastern Brazil and other portions of the species' distribution, with complete separation between southern and northern components of the Atlantic Forest. We estimated that the two groups diverged from each other c. 400,000-600,000 years ago, and have had contrasting population histories. A recent demographic expansion was inferred for the southern group, while northern populations seem to have had a longer history of large population size. Nuclear sequence data did not support this north-south pattern of subdivision, likely due at least in part to secondary male-mediated historical gene flow, inferred from multilocus coalescent-based analyses. We have compared the inferred phylogeographic patterns to those observed for other Neotropical vertebrates, and report evidence for a major north-south demographic discontinuity that seems to have marked the history of the Atlantic Forest biota.     

Persistence and diversification of the Holarctic shrew,Sorex tundrensis (Family Soricidae), in response to climate change

Environmental processes govern demography, species movements, community turnover and diversification and yet in many respects these dynamics are still poorly understood at high latitudes. We investigate the combined effects of climate change and geography through time for a widespread Holarctic shrew, Sorex tundrensis. We include a comprehensive suite of closely related outgroup taxa and three independent loci to explore phylogeographic structure and historical demography. We then explore the implications of these findings for other members of boreal communities. The tundra shrew and its sister species, the Tien Shan shrew (Sorex asper), exhibit strong geographic population structure across Siberia and into Beringia illustrating local centres of endemism that correspond to Late Pleistocene refugia. Ecological niche predictions for both current and historical distributions indicate a model of persistence through time despite dramatic climate change. Species tree estimation under a coalescent process suggests that isolation between populations has been maintained across timeframes deeper than the periodicity of Pleistocene glacial cycling. That some species such as the tundra shrew have a history of persistence largely independent of changing climate, whereas other boreal species shifted their ranges in response to climate change, highlights the dynamic processes of community assembly at high latitudes.     

Genetic evidence on the demography of speciation in allopatric dolphin species

Under a neutral model, the stochastic lineage sorting that leads to gene monophyly proceeds slowly in large populations. Therefore, in many recent species with large population size, the genome will have mixed support for monophyly unless historical bottlenecks have accelerated coalescence. We use genealogical patterns in mitochondrial DNA and in introns of four nuclear loci to test for historical bottlenecks during the speciation and divergence of two temperate Lagenorhynchus dolphin species isolated by tropical Pacific waters (an antitropical distribution). Despite distinct morphologies, foraging behaviors, and mitochondrial DNAs, these dolphin species are polyphyletic at all four nuclear loci. The abundance of shared polymorphisms between these sister taxa is most consistent with the maintenance of large effective population sizes (5.09 x 10(4) to 10.9 x 10(4)) during 0.74-1.05 million years of divergence. A variety of population size histories are possible, however. We used gene tree coalescent probabilities to explore the rejection region for historical bottlenecks of different intensity given best estimates of effective population size under a strict isolation model of divergence. In L. obliquidens the data are incompatible with a colonization propagule of an effective size of 10 or fewer individuals. Although the ability to reject less extreme historical bottlenecks will require data from additional loci, the intermixed genealogical patterns observed between these dolphin sister species are highly probable only under an extended history of large population size. If similar demographic histories are inferred for other marine antitropical taxa, a parsimonious model for the Pleistocene origin of these distributions would not involve rare breaches of a constant dispersal barrier by small colonization propagules. Instead, a history of large population size in L. obliquidens and L. obscurus contributes to growing biological and environmental evidence that the equatorial barrier became permeable during glacial/interglacial cycles, leading to vicariant isolation of antitropical populations.     

Refugial origins of reindeer (Rangifer tarandus L.) inferred from mitochondrial DNA sequences

The glacial-interglacial cycles of the upper Pleistocene have had a major impact on the recent evolutionary history of Arctic species. To assess the effects of these large-scale climatic fluctuations to a large, migratory Arctic mammal, we assessed the phylogeography of reindeer (Rangifer tarandus) as inferred from mitochondrial DNA (mtDNA) sequence variation in the control region. Phylogenetic relationships among haplotypes seem to reflect historical patterns of fragmentation and colonization rather than clear-cut relationships among extant populations and subspecies. Three major haplogroups were detected, presumably representing three separate populations during the last glacial. The most influential one has contributed to the gene pool of all extant subspecies and seems to represent a large and continuous glacial population extending from Beringia and far into Eurasia. A smaller, more localized refugium was most likely isolated in connection with ice expansion in western Eurasia. A third glacial refugium was presumably located south of the ice sheet in North America, possibly comprising several separate refugial populations. Significant demographic population expansion was detected for the two haplogroups representing the western Eurasian and Beringian glacial populations. The former apparently expanded when the ice cap retreated by the end of the last glacial. The large continuous one, in contrast, seems to have expanded by the end of the last interglacial, indicating that the warm interglacial climate accompanied by marine transgression and forest expansion significantly confined population size on the continental mainland. Our data demonstrate that the current subspecies designation does not reflect the mtDNA phylogeography of the species, which in turn may indicate that morphological differences among subspecies have evolved as adaptive responses to postglacial environmental change.     

Simultaneous extraction of nuclear and mitochondrial DNA from human blood

A method of simultaneous isolation of nuclear DNA (nDNA) and mitochondrial DNA (mtDNA) from human blood has been proposed by improvising Lahiri's method of isolation of nuclear DNA. The approach presented here provides selectively enriched fractions and eliminates the need for two different methods or separate reagent sets for the extraction of nDNA and mtDNA. It employs an initial nuclear/ cytoplasm partitioning, followed by the similar procedural steps for the two fractions separately. It gives good quality and quantity of the nDNA as well as the mtDNA, suitable for processes like PCR amplification and sequencing and may prove to be useful for people studying population genetics and evolution using molecular markers maximizing the available resources, especially in cases where a large database needs to be generated from limited amount of blood sample. From 3 ml of blood, the yields of mtDNA salvaged from the supernatant were sufficient to set approximately 4x10(5) reactions (starting with 250 fg DNA per reactions) of mtDNA loci which otherwise would have been discarded as per original Lahiri's procedure. The quality of mtDNA from the mitochondrial fraction was suitable for all major downstream processes as confirmed by locus specific PCR amplifications and sequencing. Through this procedure, the wastage of nDNA can be avoided when mtDNA loci is studied.     

Comparative population structure of Cynopterus fruit bats in peninsular Malaysia and southern Thailand

The extent to which response to environmental change is mediated by species-specific ecology is an important aspect of the population histories of tropical taxa. During the Pleistocene glacial cycles and associated sea level fluctuations, the Sunda region in Southeast Asia experienced concurrent changes in landmass area and the ratio of forest to open habitat, providing an ideal setting to test the expectation that habitat associations played an important role in determining species' response to the opportunity for geographic expansion. We used mitochondrial control region sequences and six microsatellite loci to compare the phylogeographic structure and demographic histories of four broadly sympatric species of Old World fruit bats in the genus, Cynopterus. Two forest-associated species and two open-habitat generalists were sampled along a latitudinal transect in Singapore, peninsular Malaysia, and southern Thailand. Contrary to expectations based on habitat associations, the geographic scale of population structure was not concordant across ecologically similar species. We found evidence for long and relatively stable demographic history in one forest and one open-habitat species, and inferred non-coincident demographic expansions in the second forest and open-habitat species. Thus, while these results indicate that Pleistocene climate change did not have a single effect on population structure across species, a correlation between habitat association and response to environmental change was supported in only two of four species. We conclude that interactions between multiple factors, including historical and contemporary environmental change, species-specific ecology and interspecific interactions, have shaped the recent evolutionary histories of Cynopterus fruit bats in Southeast Asia.     

Contemporary patterns in a historical context: phylogeographic history of the pipevine swallowtail,Battus philenor (Papilionidae)

Abstract We examined mitochondrial DNA (mtDNA) variation in pipevine swallowtail butterflies ( Battus philenor ) from throughout its extant range to provide a historical, phylogeographical context for ecological studies of the disjunct population in California. We evaluate current hypotheses regarding host plant use, behavior, and mimetic relationships of B. philenor populations and generate alternative hypotheses. Compared to populations throughout the rest of the species' range, California populations are ecologically distinct in that they lack mimics, lay significantly larger clutches of eggs, and exclusively use a unique, endemic larval host plant. Analysis of molecular variance, tests of population differentiation, and nested clade analysis of mtDNA variation indicate that, despite low levels of population genetic structure across the species' range, there is evidence of recent range expansion from presumed Pleistocene refuge(s) in southeastern North America. Colonization of California appears to have been a recent event. This phylogeographic investigation also suggests that the evolution of life-history adaptations to a novel larval host has occurred rapidly in California and the lack of mimics in California may be attributable to the recency of colonization.     

Multiple nuclear gene sequences identify phylogenetic species boundaries in the rapidly radiating clade of Mexican ambystomatid salamanders

Delimiting the boundaries of species involved in radiations is critical to understanding the tempo and mode of lineage formation. Single locus gene trees may or may not reflect the underlying pattern of population divergence and lineage formation, yet they constitute the vast majority of the empirical data in species radiations. In this study we make use of an expressed sequence tag (EST) database to perform nuclear (nDNA) and mitochondrial (mtDNA) genealogical tests of species boundaries in Ambystoma ordinarium, a member of an adaptive radiation of metamorphic and paedomorphic salamanders (the Ambystoma tigrinum complex) that have diversified across terrestrial and aquatic environments. Gene tree comparisons demonstrate extensive nonmonophyly in the mtDNA genealogy of A. ordinarium, while seven of eight independent nuclear loci resolve the species as monophyletic or nearly so, and diagnose it as a well-resolved genealogical species. A differential introgression hypothesis is supported by the observation that western A. ordinarium localities contain mtDNA haplotypes that are identical or minimally diverged from haplotypes sampled from a nearby paedomorphic species, Ambystoma dumerilii, while most nDNA trees place these species in distant phylogenetic positions. These results provide a strong example of how historical introgression can lead to radical differences between gene trees and species histories, even among currently allopatric species with divergent life history adaptations and morphologies. They also demonstrate how EST-based nuclear resources can be used to more fully resolve the phylogenetic history of species radiations.     

Genetic structuring and recent demographic history of red pandas (Ailurus fulgens) inferred from microsatellite and mitochondrial DNA

Clarification of the genetic structure and population history of a species can shed light on the impacts of landscapes, historical climate change and contemporary human activities and thus enables evidence‐based conservation decisions for endangered organisms. The red panda (Ailurus fulgens) is an endangered species distributing at the edge of the Qinghai‐Tibetan Plateau and is currently subject to habitat loss, fragmentation and population decline, thus representing a good model to test the influences of the above‐mentioned factors on a plateau edge species. We combined nine microsatellite loci and 551 bp of mitochondrial control region (mtDNA CR) to explore the genetic structure and demographic history of this species. A total of 123 individuals were sampled from 23 locations across five populations. High levels of genetic variation were identified for both mtDNA and microsatellites. Phylogeographic analyses indicated little geographic structure, suggesting historically wide gene flow. However, microsatellite‐based Bayesian clustering clearly identified three groups (Qionglai‐Liangshan, Xiaoxiangling and Gaoligong‐Tibet). A significant isolation‐by‐distance pattern was detected only after removing Xiaoxiangling. For mtDNA data, there was no statistical support for a historical population expansion or contraction for the whole sample or any population except Xiaoxiangling where a signal of contraction was detected. However, Bayesian simulations of population history using microsatellite data did pinpoint population declines for Qionglai, Xiaoxiangling and Gaoligong, demonstrating significant influences of human activity on demography. The unique history of the Xiaoxiangling population plays a critical role in shaping the genetic structure of this species, and large‐scale habitat loss and fragmentation is hampering gene flow among populations. The implications of our findings for the biogeography of the Qinghai‐Tibetan Plateau, subspecies classification and conservation of red pandas are discussed.     

A reevaluation of the Native American mtDNA genome diversity and its bearing on the models of early colonization of Beringia

The Americas were the last continents to be populated by humans, and their colonization represents a very interesting chapter in our species' evolution in which important issues are still contentious or largely unknown. One difficult topic concerns the details of the early peopling of Beringia, such as for how long it was colonized before people moved into the Americas and the demography of this occupation. A recent work using mitochondrial genome (mtDNA) data presented evidence for a so called "three-stage model" consisting of a very early expansion into Beringia followed by approximately 20,000 years of population stability before the final entry into the Americas. However, these results are in disagreement with other recent studies using similar data and methods. Here, we reanalyze their data to check the robustness of this model and test the ability of Native American mtDNA to discriminate details of the early colonization of Beringia. We apply the Bayesian Skyline Plot approach to recover the past demographic dynamic underpinning these events using different mtDNA data sets. Our results refute the specific details of the "three-stage model", since the early stage of expansion into Beringia followed by a long period of stasis could not be reproduced in any mtDNA data set cleaned from non-Native American haplotypes. Nevertheless, they are consistent with a moderate population bottleneck in Beringia associated with the Last Glacial Maximum followed by a strong population growth around 18,000 years ago as suggested by other recent studies. We suggest that this bottleneck erased the signals of ancient demographic history from recent Native American mtDNA pool, and conclude that the proposed early expansion and occupation of Beringia is an artifact caused by the misincorporation of non-Native American haplotypes.     

Mitochondrial DNA under siege in avian phylogeography

Mitochondrial DNA (mtDNA) has been the workhorse of research in phylogeography for almost two decades. However, concerns with basing evolutionary interpretations on mtDNA results alone have been voiced since the inception of such studies. Recently, some authors have suggested that the potential problems with mtDNA are so great that inferences about population structure and species limits are unwarranted unless corroborated by other evidence, usually in the form of nuclear gene data. Here we review the relative merits of mitochondrial and nuclear phylogeographical studies, using birds as an exemplar class of organisms. A review of population demographic and genetic theory indicates that mitochondrial and nuclear phylogeographical results ought to concur for both geographically unstructured populations and for populations that have long histories of isolation. However, a relatively common occurrence will be shallow, but geographically structured mtDNA trees--without nuclear gene corroboration--for populations with relatively shorter periods of isolation. This is expected because of the longer coalescence times of nuclear genes (approximately four times that of mtDNA); such cases do not contradict the mtDNA inference of recent isolation and evolutionary divergence. Rather, the nuclear markers are more lagging indicators of changes in population structure. A review of the recent literature on birds reveals the existence of relatively few cases in which nuclear markers contradict mitochondrial markers in a fashion not consistent with coalescent theory. Preliminary information from nuclear genes suggests that mtDNA patterns will prove to be robust indicators of patterns of population history and species limits. At equilibrium, mitochondrial loci are generally a more sensitive indicator of population structure than are nuclear loci, and mitochondrial estimates of F(ST)-like statistics are generally expected to exceed nuclear ones. Hence, invoking behavioural or ecological explanations of such differences is not parsimonious. Nuclear genes will prove important for quantitative estimates of the depths of haplotype trees, rates of population growth and values of gene flow.     

PERCHED AT THE MITO‐NUCLEAR CROSSROADS: DIVERGENT MITOCHONDRIAL LINEAGES CORRELATE WITH ENVIRONMENT IN THE FACE OF ONGOING NUCLEAR GENE FLOW IN AN AUSTRALIAN BIRD

Relationships among multilocus genetic variation, geography, and environment can reveal how evolutionary processes affect genomes. We examined the evolution of an Australian bird, the eastern yellow robin Eopsaltria australis, using mitochondrial (mtDNA) and nuclear (nDNA) genetic markers, and bioclimatic variables. In southeastern Australia, two divergent mtDNA lineages occur east and west of the Great Dividing Range, perpendicular to latitudinal nDNA structure. We evaluated alternative scenarios to explain this striking discordance in landscape genetic patterning. Stochastic mtDNA lineage sorting can be rejected because the mtDNA lineages are essentially distinct geographically for > 1500 km. Vicariance is unlikely: the Great Dividing Range is neither a current barrier nor was it at the Last Glacial Maximum according to species distribution modeling; nuclear gene flow inferred from coalescent analysis affirms this. Female philopatry contradicts known female‐biased dispersal. Contrasting mtDNA and nDNA demographies indicate their evolutionary histories are decoupled. Distance‐based redundancy analysis, in which environmental temperatures explain mtDNA variance above that explained by geographic position and isolation‐by‐distance, favors a nonneutral explanation for mitochondrial phylogeographic patterning. Thus, observed mito‐nuclear discordance accords with environmental selection on a female‐linked trait, such as mtDNA, mtDNA–nDNA interactions or genes on W‐chromosome, driving mitochondrial divergence in the presence of nuclear gene flow.     

Evolutionary relationship between marble trout of the northern and the southern Adriatic basin

Marble trout (Salmo marmoratus) populate two geographically separated areas in the northern and southern parts of the Adriatic Sea drainage. Although morphologically similar, each population is distinguished by a different set of unrelated mitochondrial haplotypes, suggesting that they have evolved from different ancestors. Due to a possible discordance between mitochondrial and species phylogeny, we performed phylogenetic analysis based on 22 nuclear loci. The results inferred from Maximum-likelihood and Bayesian Inference analysis revealed that northern and southern populations are closely related, forming a monophyletic group. This observation is concordant with the present marble trout classification, which considers both populations as conspecific. On the other hand, our findings are in marked contrast to those of previous mtDNA-based studies and highlight potential dangers of making phylogenetic inferences from mtDNA alone. Reasons for discordance between mtDNA and nDNA phylogeny are discussed with incomplete lineage sorting proposed as the most parsimonious explanation for mtDNA divergence in marble trout.     

Population genetic consequences of geographic disjunction: a prairie plant isolated on Great Lakes alvars

Species may often exhibit geographic variation in population genetic structure due to contemporary and historical variation in population size and gene flow. Here, we test the predictions that populations on the margins of a species' distribution contain less genetic variation and are more differentiated than populations towards the core of the range by comparing patterns of genetic variation at five microsatellite loci between disjunct and core populations of the perennial, allohexaploid herb Geum triflorum. We sampled nine populations isolated on alvar habitat within the eastern Great Lakes region in North America, habitats that include disjunct populations of several plant species, and compared these to 16 populations sampled from prairie habitat throughout the core of the species' distribution in midwestern Canada and the USA. Alvar populations exhibited much lower within-population diversity and contained only a subset of alleles found in prairie populations. We detected isolation by distance across the species' range and within alvar and prairie regions separately. As predicted, genetic differentiation was higher among alvar populations than among prairie populations, even after controlling for the geographic distance between sampled populations. Low diversity and high differentiation can be accounted for by the greater contemporary spatial isolation of alvar populations. However, the genetic structure of alvar populations may also have been influenced by postglacial range expansion and contraction. Our results are consistent with alvar populations being founded during an expansion of prairie habitat during the warmer, hypsithermal period approximately 5000 bp and subsequently becoming stranded on isolated alvar habitat as the climate grew cooler and wetter.     

Indo-Pacific population structure and evolutionary history of the coconut crab Birgus latro

Mitochondrial DNA variation was used to examine population structure in a widespread, marine-dispersed species, Birgus latro . Crabs were collected from eight locations throughout the species' Indo-Pacific distribution. Purified mtDNA from 160 individuals was cut with five restriction enzymes, revealing high haplotype diversity (0.96) and moderate nucleotide diversity (0.75%). Island populations from the Indian Ocean (Christmas I.) and Pacific Ocean were significantly different ( G _ST= 0.37) and had distinct mtDNA lineages with a net sequence divergence of 1.4%. Pacific island populations had diverged in a manner consistent with isolation by distance, with only the most peripheral populations being significantly different. The results for mtDNA are largely concordant with those from allozymes, although estimates of gene flow between the Indian and Pacific Oceans were much lower when based on mtDNA. The mtDNA phylogeny also permitted a deeper examination of the evolutionary and demographic history of Birgus latro . Long-term separation of populations is evident in the complete phylogenetic subdivision of mtDNA lineages between the Indian and Pacific Ocean populations sampled. The starlike phylogeny of alleles from the Pacific suggests a rapid population expansion in the Pacific during the Pleistocene. Including information about allele phylogeny, as well as distribution and frequency, obscured contemporary population structure, but provided unique insights into the evolutionary history of the species.     

Consistent geographic structure among multiple nuclear sequences and cpDNA polymorphisms of Cardamine nipponica Franch. et Savat. (Brassicaceae)

Molecular phylogeography has inferred the history of differentiation between regions and/or among populations following the Pleistocene climatic oscillations, mostly based on the genetic structure of organelle DNA. However, such genetic structure only reflects the history of a single gene, and studies based on single-copy genes of nuclear DNA (nDNA) are required for phylogeography, although their efficiency remains unclear. To examine the utility of nDNA loci, the genetic structures of three genes from Cardamine nipponica, which is closely related to the model species Arabidopsis thaliana, were elucidated: the nDNA genes DET1, PHYA, PHYE, as well as chloroplast DNA (cpDNA). In 279 individuals collected from throughout the range of the species, strong genetic differentiation between northern and central Japan was found for all loci. This result suggested that populations in central Japan experienced a different history from those in northern Japan during the Pleistocene climatic oscillations. In addition, the evidence of refugia at the edges of the distribution, where the genetic structure was less influenced by colonization following range expansion, was shown for several loci. The specific genetic structure within the southernmost populations of northern Japan suggested that this region was also isolated during range expansion. Hence, the consistent history among loci and a more detailed history from several loci indicated that cpDNA can represent the history of vicariance and demonstrated the efficiency of single-copy nuclear genes in phylogeography.     

Genetic differentiation across a latitudinal gradient in two co-occurring butterfly species: revealing population differences in a context of climate change

Genetic differentiation within a species' range is determined by natural selection, genetic drift, and gene flow. Selection and drift enhance genetic differences if populations are sufficiently isolated, while gene flow precludes differentiation and local adaptation. Over large geographical areas, these processes can create a variety of scenarios, ranging from admixture to a high degree of population differentiation. Genetic differences among populations may signal functional differences within a species' range, potentially leading to population or ecotype-specific responses to global change. We investigated differentiation within the geographical range of two butterfly species along a broad latitudinal gradient. This gradient is the primary axis of climatic variation, and many ecologists expect populations at the poleward edge of this gradient to expand under climate change. Our study species inhabit a shared ecosystem and differ in body size and resource specialization; both also find their poleward range limit on an island. We find evidence for divergence of peripheral populations from the core in both taxa, suggesting the potential for genetic distinctiveness at the leading edge of climate change. We also find differences between the species in the extent of peripheral differentiation with the smaller and more specialized species showing greater population divergence (microsatellites and mtDNA) and reduced gene flow (mtDNA). Finally, gene flow estimates in both species differed strongly between two marker types. These findings suggest caution in assuming that populations are invariant across latitude and thus will respond as a single ecotype to climatic change.     

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.     

Contrasting signatures of population growth for mitochondrial DNA and Y chromosomes among human populations in Africa

A history of Pleistocene population expansion has been inferred from the frequency spectrum of polymorphism in the mitochondrial DNA (mtDNA) of many human populations. Similar patterns are not typically observed for autosomal and X-linked loci. One explanation for this discrepancy is a recent population bottleneck, with different rates of recovery for haploid and autosomal loci as a result of their different effective population sizes. This hypothesis predicts that mitochondrial and Y chromosomal DNA will show a similar skew in the frequency spectrum in populations that have experienced a recent increase in effective population size. We test this hypothesis by resequencing 6.6 kb of noncoding Y chromosomal DNA and 780 basepairs of the mtDNA cytochrome c oxidase subunit III (COIII) gene in 172 males from 5 African populations. Four tests of population expansion are employed for each locus in each population: Fu's Fs statistic, the R(2) statistic, coalescent simulations, and the mismatch distribution. Consistent with previous results, patterns of mtDNA polymorphism better fit a model of constant population size for food-gathering populations and a model of population expansion for food-producing populations. In contrast, none of the tests reveal evidence of Y chromosome growth for either food-gatherers or food-producers. The distinct mtDNA and Y chromosome polymorphism patterns most likely reflect sex-biased demographic processes in the recent history of African populations. We hypothesize that males experienced smaller effective population sizes and/or lower rates of migration during the Bantu expansion, which occurred over the last 5,000 years.