Subterranean termite phylogeography reveals multiple postglacial colonization events in southwestern Europe

A long‐standing goal of evolutionary biology is to understand how paleoclimatic and geological events shape the geographical distribution and genetic structure within and among species. Using a diverse set of markers (cuticular hydrocarbons, mitochondrial and nuclear gene sequences, microsatellite loci), we studied Reticulitermes grassei and R. banyulensis, two closely related termite species in southwestern Europe. We sought to clarify the current genetic structure of populations that formed following postglacial dispersal from refugia in southern Spain and characterize the gene flow between the two lineages over the last several million years. Each marker type separately provided a fragmented picture of the evolutionary history at different timescales. Chemical analyses of cuticular hydrocarbons and phylogenetic analyses of mitochondrial and nuclear genes showed clear separation between the species, suggesting they diverged following vicariance events in the Late Miocene. However, the presence of intermediate chemical profiles and mtDNA introgression in some Spanish colonies suggests ongoing gene flow. The current genetic structure of Iberian populations is consistent with alternating isolation and dispersal events during Quaternary glacial periods. Analyses of population genetic structure revealed postglacial colonization routes from southern Spain to France, where populations underwent strong genetic bottlenecks after traversing the Pyrenees resulting in parapatric speciation.     

Postglacial recolonization shaped the genetic diversity of the winter moth (Operophtera brumata) in Europe

Changes in climate conditions, particularly during the Quaternary climatic oscillations, have long been recognized to be important for shaping patterns of species diversity. For species residing in the western Palearctic, two commonly observed genetic patterns resulting from these cycles are as follows: (1) that the numbers and distributions of genetic lineages correspond with the use of geographically distinct glacial refugia and (2) that southern populations are generally more diverse than northern populations (the “southern richness, northern purity” paradigm). To determine whether these patterns hold true for the widespread pest species the winter moth (Operophtera brumata), we genotyped 699 individual winter moths collected from 15 Eurasian countries with 24 polymorphic microsatellite loci. We find strong evidence for the presence of two major genetic clusters that diverged ~18 to ~22 ka, with evidence that secondary contact (i.e., hybridization) resumed ~ 5 ka along a well‐established hybrid zone in Central Europe. This pattern supports the hypothesis that contemporary populations descend from populations that resided in distinct glacial refugia. However, unlike many previous studies of postglacial recolonization, we found no evidence for the “southern richness, northern purity” paradigm. We also find evidence for ongoing gene flow between populations in adjacent Eurasian countries, suggesting that long‐distance dispersal plays an important part in shaping winter moth genetic diversity. In addition, we find that this gene flow is predominantly in a west‐to‐east direction, suggesting that recently debated reports of cyclical outbreaks of winter moth spreading from east to west across Europe are not the result of dispersal.     

Biogeography meets conservation: the genetic structure of the endangered lycaenid butterfly Lycaena helle (Denis & Schiffermüller, 1775)

Cold‐adapted species are thought to have had their largest distribution ranges in central Europe during the glacial periods. Postglacial warming caused severe range shifts of such taxa into higher latitudes and altitudes. We selected the boreomontane butterfly Lycaena helle (Denis & Schiffermüller, 1775) as an example to demonstrate the genetic effects of range changes, and to document the recent status of highly fragmented remnant populations. We analysed five polymorphic microsatellite loci in 1059 individuals sampled at 50 different localities scattered over the European distribution area of the species. Genetic differentiation was strong among the mountain ranges of western Europe, but we did not detect similarly distinct genetic groups following a geographical pattern in the more eastern areas. The Fennoscandian populations form a separate genetic group, and provide evidence for a colonization from southern Finland via northern Scandinavia to south‐central Sweden. Species distribution modelling suggests a large extension of the spatial distribution during the last glacial maximum, but highlights strong retractions to a few mountain areas under current conditions. These findings, combined with our genetic data, suggest a more or less continuous distribution of L. helle throughout central Europe at the end of the last ice age. As a consequence of postglacial warming, the species retreated northwards to Fennoscandia and escaped increasing temperatures through altitudinal shifts. Therefore, the species is today restricted to population remnants located at the mountain tops of western Europe, genetically isolated from each other, and evolved into genetically unique entities. Rising temperatures and advancing habitat destruction threaten this wealth of biodiversity. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 101 , 155–168.     

Hierarchical analysis of genetic structure in the habitat‐specialist Eastern Sand Darter (Ammocrypta pellucida)

Quantifying spatial genetic structure can reveal the relative influences of contemporary and historic factors underlying localized and regional patterns of genetic diversity and gene flow – important considerations for the development of effective conservation efforts. Using 10 polymorphic microsatellite loci, we characterize genetic variation among populations across the range of the Eastern Sand Darter (Ammocrypta pellucida), a small riverine percid that is highly dependent on sandy substrate microhabitats. We tested for fine scale, regional, and historic patterns of genetic structure. As expected, significant differentiation was detected among rivers within drainages and among drainages. At finer scales, an unexpected lack of within‐river genetic structure among fragmented sandy microhabitats suggests that stratified dispersal resulting from unstable sand bar habitat degradation (natural and anthropogenic) may preclude substantial genetic differentiation within rivers. Among‐drainage genetic structure indicates that postglacial (14 kya) drainage connectivity continues to influence contemporary genetic structure among Eastern Sand Darter populations in southern Ontario. These results provide an unexpected contrast to other benthic riverine fish in the Great Lakes drainage and suggest that habitat‐specific fishes, such as the Eastern Sand Darter, can evolve dispersal strategies that overcome fragmented and temporally unstable habitats.     

A complex scenario of glacial survival in Mediterranean and continental refugia of a temperate continental vole species (Microtus arvalis) in Europe

The role of glacial refugia in shaping contemporary species distribution is a long-standing question in phylogeography and evolutionary ecology. Recent studies are questioning previous paradigms on glacial refugia and postglacial recolonization pathways in Europe, and more flexible phylogeographic scenarios have been proposed. We used the widespread common vole Microtus arvalis as a model to investigate the origin, locations of glacial refugia, and dispersal pathways, in the group of “Continental” species in Europe. We used a Bayesian spatiotemporal diffusion analysis (relaxed random walk model) of cytochrome b sequences across the species range, including newly collected individuals from 10 Iberian localities and published sequences from 68 localities across 22 European countries. Our data suggest that the species originated in Central Europe, and we revealed the location of multiple refugia (in both southern peninsulas and continental regions) for this continental model species. Our results confirm the monophyly of Iberian voles and the pre-LGM divergence between Iberian and European voles. We found evidence of restricted postglacial dispersal from refugia in Mediterranean peninsulas. We inferred a complex evolutionary and demographic history of M. arvalis in Europe over the last 50,000 years that does not adequately fit previous glacial refugial scenarios. The phylogeography of M. arvalis provides a paradigm of ice-age survival of a temperate continental species in western and eastern Mediterranean peninsulas (sources of endemism) and multiple continental regions (sources of postglacial spread). Our findings also provide support for a major role of large European river systems in shaping geographic boundaries of M. arvalis in Europe.     

Investigating past range dynamics for a weed of cultivation,Silene vulgaris

Since the last glacial maximum (LGM), many plant and animal taxa have expanded their ranges by migration from glacial refugia. Weeds of cultivation may have followed this trend or spread globally following the expansion of agriculture or ruderal habitats associated with human‐mediated disturbance. We tested whether the range expansion of the weed Silene vulgaris across Europe fit the classical model of postglacial expansion from southern refugia, or followed known routes of the expansion of human agricultural practices. We used species distribution modeling to predict spatial patterns of postglacial expansion and contrasted these with the patterns of human agricultural expansion. A population genetic analysis using microsatellite loci was then used to test which scenario was better supported by spatial patterns of genetic diversity and structure. Genetic diversity was highest in southern Europe and declined with increasing latitude. Locations of ancestral demes from genetic cluster analysis were consistent with areas of predicted refugia. Species distribution models showed the most suitable habitat in the LGM on the southern coasts of Europe. These results support the typical postglacial northward colonization from southern refugia while refuting the east‐to‐west agricultural spread as the main mode of expansion for S. vulgaris. We know that S. vulgaris has recently colonized many regions (including North America and other continents) through human‐mediated dispersal, but there is no evidence for a direct link between the Neolithic expansion of agriculture and current patterns of genetic diversity of S. vulgaris in Europe. Therefore, the history of range expansion of S. vulgaris likely began with postglacial expansion after the LGM, followed by more recent global dispersal by humans.     

Phylogeography and paleodistribution models of a widespread birch (Betula platyphylla Suk.) across East Asia: Multiple refugia,multidirectional expansion,and heterogeneous genetic pattern

Widespread tree species cover large geographical areas and play important roles in various vegetation types. Understanding how these species responded to historical climatic changes is important for understanding community assembly mechanisms with evolutionary and conservation implications. However, the location of refugial areas and postglacial history of widespread trees in East Asia remain poorly known. We combined microsatellite data (63 populations, 1756 individuals) and ecological niche modeling to examine the range‐wide population diversity, genetic structure, and historical demography of a pioneer tree species, Asian white birch (Betula platyphylla Suk.) across East Asia. We found a north‐to‐south trend of declining genetic diversity and five clusters, corresponding to geographical regions. Different clusters were inferred to have diverged through Pleistocene climatic oscillations and have different expansion routes, leading to genetic admixture in some populations. Ecological niche models indicated that the distribution of B. platyphylla during the last glacial maximum still had a large latitude span with slight shifts toward southeast, and northern populations had more variable distribution ranges than those in the south during later climatic oscillations. Our results reflect the relatively stable distribution through the last glacial–interglacial cycles and recent multidirectional expansion of B. platyphylla, providing new hypotheses for the response pattern of widespread tree species to climate change. The gradual genetic pattern from northeast to southwest and alternative distribution dynamics possibly resulted from environmental differences caused by latitude and topographic heterogeneity.     

Contrasting population genetic structure among freshwater‐resident and anadromous lampreys: the role of demographic history,differential dispersal and anthropogenic barriers to movement

The tendency of many species to abandon migration remains a poorly understood aspect of evolutionary biology that may play an important role in promoting species radiation by both allopatric and sympatric mechanisms. Anadromy inherently offers an opportunity for the colonization of freshwater environments, and the shift from an anadromous to a wholly freshwater life history has occurred in many families of fishes. Freshwater‐resident forms have arisen repeatedly among lampreys (within the Petromyzontidae and Mordaciidae), and there has been much debate as to whether anadromous lampreys, and their derived freshwater‐resident analogues, constitute distinct species or are divergent ecotypes of polymorphic species. Samples of 543 European river lamprey Lampetra fluviatilis (mostly from anadromous populations) and freshwater European brook lamprey Lampetra planeri from across 18 sites, primarily in the British Isles, were investigated for 13 polymorphic microsatellite DNA loci, and 108 samples from six of these sites were sequenced for 829 bp of mitochondrial DNA (mtDNA). We found contrasting patterns of population structure for mtDNA and microsatellite DNA markers, such that low diversity and little structure were seen for all populations for mtDNA (consistent with a recent founder expansion event), while fine‐scale structuring was evident for nuclear markers. Strong differentiation for microsatellite DNA loci was seen among freshwater‐resident L. planeri populations and between L. fluviatilis and L. planeri in most cases, but little structure was evident among anadromous L. fluviatilis populations. We conclude that postglacial colonization founded multiple freshwater‐resident populations with strong habitat fidelity and limited dispersal tendencies that became highly differentiated, a pattern that was likely intensified by anthropogenic barriers.     

Drifting fronds and drifting alleles: range dynamics,local dispersal and habitat isolation shape the population structure of the estuarine seaweed Fucus ceranoides

Aim The seaweed Fucus ceranoides is restricted to spatially discrete estuarine habitats and lacks planktonic dispersal phases; it is therefore expected to exhibit strong population differentiation. Its cold‐temperate affinities and mtDNA variation imply that the northern part of the species’ range, where F. ceranoides is now ubiquitous, was recently colonized after the onset of the last deglaciation, potentially resulting in areas with greater genetic homogeneity. Here we examine the population structure of F. ceranoides to test these predictions, emphasizing the contrasting genetic signatures of limited dispersal in refugial versus recently colonized regions. Location North‐eastern Atlantic estuaries from Portugal to Norway. Methods A total of 504 individuals from 21 estuarine sites spanning the entire range of F. ceranoides were sampled and genotyped for nine microsatellite loci. Population structure was inferred from several genotypic and allele‐frequency analyses. Geographical patterns of genetic diversity were used to reconstruct the historical biogeography of the species. Results Genetic diversity and regional population differentiation showed a consistent decline with increasing latitude. Southernmost populations harboured most of the endemic variation, whereas the northern populations (> 55° N) were almost fixed for the same alleles across loci. In southern and central regions of its distribution, F. ceranoides showed striking population subdivision, with many of the sampled estuaries corresponding to coherent genetic units that were easily discriminated from one another with standard clustering methods. Main conclusions The geographical pattern of genetic diversity supports the long‐term refugial status of Iberia and a post‐glacial range expansion of F. ceranoides into previously glaciated latitudes. Despite the species’ capacity to colonize newly available habitats, the genetic structure of F. ceranoides outside the recently (re)colonized range reveals that gene flow between populations is extremely low. This study provides a remarkable example of how infrequent and spatially limited dispersal can have contrasting effects at the scales of meta‐population (connectivity) versus range dynamics (habitat tracking), and of how dispersal restrictions can result in either genetic divergence or homogeneity depending on the maturity and demographic conditions of the populations.     

Molecular Markers Reveal Limited Population Genetic Structure in a North American Corvid,Clark’s Nutcracker (Nucifraga columbiana)

The genetic impact of barriers and Pleistocene glaciations on high latitude resident species has not been widely investigated. The Clark’s nutcracker is an endemic North American corvid closely associated with Pinus-dominated forests. The nutcracker’s encompasses known barriers to dispersal for other species, and glaciated and unglaciated areas. Clark’s nutcrackers also irruptively disperse long distances in search of pine seed crops, creating the potential for gene flow among populations. Using the highly variable mitochondrial DNA control region, seven microsatellite loci, and species distribution modeling, we examined the effects of glaciations and dispersal barriers on population genetic patterns and population structure of nutcrackers. We sequenced 900 bp of mitochondrial control region for 169 individuals from 15 populations and analysed seven polymorphic microsatellite loci for 13 populations across the Clark’s nutcracker range. We used species distribution modeling and a range of phylogeographic analyses to examine evolutionary history. Clark’s nutcracker populations are not highly differentiated throughout their range, suggesting high levels of gene flow among populations, though we did find some evidence of isolation by distance and peripheral isolation. Our analyses suggested expansion from a single refugium after the last glacial maximum, but patterns of genetic diversity and paleodistribution modeling of suitable habitat were inconclusive as to the location of this refugium. Potential barriers to dispersal (e.g. mountain ranges) do not appear to restrict gene flow in Clark’s nutcracker, and postglacial expansion likely occurred quickly from a single refugium located south of the ice sheets.     

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.     

Isolation,characterization and cross‐salmonid amplification of 31 microsatellite loci in the lake whitefish (Coregonus clupeaformis,Mitchill)

Coregonine fish represent the most successful evolutionary lineage of salmonids with Coregonus as the most speciose salmonid genus inhabiting numerous postglacial lakes across the northern hemisphere. We isolated and characterized 31 polymorphic microsatellite loci in Coregonus clupeaformis with an average number of 5.3 alleles per locus (range three to eight) and an overall expected heterozygosity of 0.74 ± 0.11. Two loci revealed significant linkage associations through analyses of mapping families. Six additional salmonid taxa assessed for cross‐species amplification revealed between 18 and 26 positive amplifications and between two and 12 polymorphic loci per species.     

Isolation and characterization of microsatellite markers in the tetraploid birch,Betula pubescens ssp. tortuosa

The mountain birch, Betula pubescens ssp. tortuosa is a tetraploid tree species which forms the tree limit in northern Fennoscandia. We identified nine polymorphic microsatellite loci in order to characterize the genetic structure of populations at different elevations close to the treeline. The microsatellite loci were highly polymorphic, with 14–42 alleles per locus and an average expected heterozygosity of 0.73 ± 0.25 under random chromosome segregation and 0.68 ± 0.23 under random chromatid segregation.     

Influence of northern limit range on genetic diversity and structure in a widespread North American tree,sugar maple (Acer saccharum Marshall)

Due to climate change, the ranges of many North American tree species are expected to shift northward. Sugar maple (Acer saccharum Marshall) reaches its northern continuous distributional limit in northeastern North America at the transition between boreal mixed‐wood and temperate deciduous forests. We hypothesized that marginal fragmented northern populations from the boreal mixed wood would have a distinct pattern of genetic structure and diversity. We analyzed variation at 18 microsatellite loci from 23 populations distributed along three latitudinal transects (west, central, and east) that encompass the continuous–discontinuous species range. Each transect was divided into two zones, continuous (temperate deciduous) and discontinuous (boreal mixed wood), based on sugar maple stand abundance. Respective positive and negative relationships were found between the distance of each population to the northern limit (D_north), and allelic richness (A_R) and population differentiation (F_ST). These relations were tested for each transect separately; the pattern (discontinuous–continuous) remained significant only for the western transect. structure analysis revealed the presence of four clusters. The most northern populations of each transect were assigned to a distinct group. Asymmetrical gene flow occurred from the southern into the four northernmost populations. Southern populations in Québec may have originated from two different postglacial migration routes. No evidence was found to validate the hypothesis that northern populations were remnants of a larger population that had migrated further north of the species range after the retreat of the ice sheet. The northernmost sugar maple populations possibly originated from long‐distance dispersal.     

Glacial refugia and postglacial expansion of the alpine–prealpine plant species Polygala chamaebuxus

The shrubby milkwort (Polygala chamaebuxus L.) is widely distributed in the Alps, but occurs also in the lower mountain ranges of Central Europe such as the Franconian Jura or the Bohemian uplands. Populations in these regions may either originate from glacial survival or from postglacial recolonization. In this study, we analyzed 30 populations of P. chamaebuxus from the whole distribution range using AFLP (Amplified Fragment Length Polymorphism) analysis to identify glacial refugia and to illuminate the origin of P. chamaebuxus in the lower mountain ranges of Central Europe. Genetic variation and the number of rare fragments within populations were highest in populations from the central part of the distribution range, especially in the Southern Alps (from the Tessin Alps and the Prealps of Lugano to the Triglav Massiv) and in the middle part of the northern Alps. These regions may have served, in accordance with previous studies, as long‐term refugia for the glacial survival of the species. The geographic pattern of genetic variation, as revealed by analysis of molecular variance, Bayesian cluster analysis and a PopGraph genetic network was, however, only weak. Instead of postglacial recolonization from only few long‐term refugia, which would have resulted in deeper genetic splits within the data set, broad waves of postglacial expansion from several short‐term isolated populations in the center to the actual periphery of the distribution range seem to be the scenario explaining the observed pattern of genetic variation most likely. The populations from the lower mountain ranges in Central Europe were more closely related to the populations from the southwestern and northern than from the nearby eastern Alps. Although glacial survival in the Bohemian uplands cannot fully be excluded, P. chamaebuxus seems to have immigrated postglacially from the southwestern or central‐northern parts of the Alps into these regions during the expansion of the pine forests in the early Holocene.     

Drastic post‐LGM expansion and lack of historical genetic structure of a subtropical fig‐pollinating wasp (Ceratosolen sp. 1) of Ficus septica in Taiwan

The climatic oscillations of the last glacial period have had profound influences on the demography and levels of genetic diversity of extant species. Molecular evidence of glacial effects on temperate species has been well documented, whereas little is known regarding that on subtropical species. Here we present analyses based on partial sequences of the mitochondrial cytochrome c oxidase subunit I (COI) gene (1052 bp) and genotypes at 15 microsatellite loci to investigate the historical demography, population structure and ongoing gene flow of an undescribed fig‐pollinating wasp (Ceratosolen sp. 1) of Ficus septica in subtropical Taiwan. Reconstructed historical demography based on the coalescent tree of COI sequences suggests that C. sp. 1 has undergone a drastic population expansion which was tightly coupled with climatic changes since the last glacial maximum (LGM). The magnitude of the population size change was approximately 500‐fold, indicating that the population of this wasp and its host was likely highly compressed during the last glacial period. The lack of significant population differentiation (F_ST < 0.02, for all pairwise F_ST values) may be due to rapid postglacial expansion facilitated by long‐distance dispersal, although a low frequency of first‐generation migrants was detected. Our results clearly demonstrate how recent climatic changes since the LGM and dispersal ability have jointly shaped the genetic composition of a subtropical fig‐pollinating wasp.     

Microsatellite analysis of the natterjack toad (Bufo calamita) in Denmark: populations are islands in a fragmented landscape

The European natterjack toad (Bufo calamita) has declined rapidly in recent years, primarily due to loss of habitat, and in Denmark it is estimated that 50% of the isolated populations are lost each decade. To efficiently manage and conserve this species and its genetic diversity, knowledge of the genetic structure is crucial. Based on nine polymorphic microsatellite loci, the genetic diversity, genetic structure and gene flow were investigated at 12 sites representing 5–10% of the natterjack toad localities presently known in Denmark. The expected heterozygosity (H _E) within each locality was generally low (range: 0.18–0.43). Further analyses failed to significantly correlate genetic diversity with population size, degree of isolation and increasing northern latitude, indicating a more complex combination of factors in determining the present genetic profile. Genetic differentiation was high (overall θ = 0.29) and analyses based on a Bayesian clustering method revealed that the dataset constituted 11 genetic clusters, defining nearly all sampling sites as distinct populations. Contemporary gene flow among populations was undetectable in nearly all cases, and the failure to detect a pattern of isolation by distance within major regions supported this apparent lack of a gene flow continuum. Indications of a genetic bottleneck were found in three populations. The analyses suggest that the remaining Bufo calamita populations in Denmark are genetically isolated, and represent independent units in a highly fragmented gene pool. Future conservation management of this species is discussed in light of these results. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Levels of gene flow among populations of a wolf spider in a recently fragmented habitat: current versus historical rates

Gene flow among populations is important for countering the deleterious effects of random genetic drift and inbreeding, as well as spreading beneficial mutations. Wind-driven aerial dispersal is known to occur in numerous plants and invertebrates. Its evolution suggests that historically, suitable habitat patches were dense enough to make such undirected dispersal evolutionarily advantageous. Using microsatellite markers we assessed the population genetic structure of seven populations of a wolf spider (Rabidosa rabida) capable of ballooning. Historically, each spider population received a mean of 1.5 migrants per generation from the other six populations. Over the past several generations the number of migrants reaching a population is only 0.2. This statistically significant reduction in gene flow coincides with high levels of habitat fragmentation and suggests that undirected aerial dispersal is ineffective in this fragmented landscape. Further, individuals within populations showed signficantly elevated levels of homozygosity relative to Hardy–Weinberg expectations, suggesting that cursorial dispersal may be very limited and genetic structure within populations exists. Inbreeding coefficients averaged 0.18 over all seven populations with very little variation among populations (s = 0.02). Fitness was lower in smaller populations relative to larger ones. Altered landscapes pose evolutionary dilemmas for many metapopulations and species that depend on undirected movement for dispersal may be particularly vulnerable to habitat fragmentation.     

Phylogeography of Pinus tabulaeformis Carr. (Pinaceae), a dominant species of coniferous forest in northern China

How coniferous trees in northern China changed their distribution ranges in response to Quaternary climatic oscillations remains largely unknown. Here we report a study of the phylogeography of Pinus tabulaeformis, an endemic and dominant species of coniferous forest in northern China. We examined sequence variation of maternally inherited, seed‐dispersed mitochondrial DNA (mtDNA) (nad5 intron 1 and nad4/3–4) and paternally inherited, pollen‐ and seed‐dispersed chloroplast DNA (cpDNA) (rpl16 and trnS‐trnG) within and among 30 natural populations across the entire range of the species. Six mitotypes and five chlorotypes were recovered among 291 trees surveyed. Population divergence was high for mtDNA variation (G_ST = 0.738, N_ST = 0.771) indicating low levels of seed‐based gene flow and significant phylogeographical structure (N_ST > G_ST, P < 0.05). The spatial distribution of mitotypes suggests that five distinct population groups exist in the species: one in the west comprising seven populations, a second with a north–central distribution comprising 15 populations, a third with a southern and easterly distribution comprising five populations, a fourth comprising one central and one western population, and a fifth comprising a single population located in the north‐central part of the species’ range. Each group apart from the fourth group is characterized by a distinct mitotype, with other mitotypes, if present, occurring at low frequency. It is suggested, therefore, that most members of each group apart from Group 4 are derived from ancestors that occupied different isolated refugia in a previous period of range fragmentation of the species, possibly at the time of the Last Glacial Maximum. Possible locations for these refugia are suggested. A comparison of mitotype diversity between northern and southern subgroups within the north‐central group of populations (Group 2) showed much greater uniformity in the northern part of the range both within and between populations. This could indicate a northward migration of the species from a southern refugium in this region during the postglacial period, although alternative explanations cannot be ruled out. Two chlorotypes were distributed across the geographical range of the species, resulting in lower levels of among‐population chlorotype variation. The geographical pattern of variation for all five chlorotypes provided some indication of the species surviving past glaciations in more than one refugium, although differentiation was much less marked, presumably due to the greater dispersal of cpDNA via pollen.