Microsatellite analysis reveals interpopulation differentiation and gene flow in the endangered tree Changiostyrax dolichocarpa (Styracaceae) with fragmented distribution in central China

Polymorphic simple sequence repeat (SSR) markers were used to investigate the impact of habitat fragmentation on the population structure and gene flow of Changiostyrax dolichocarpa, a critically endangered tree in central China. Intrapopulation genetic diversity, population structure and gene flow in the five extant populations of this species were analysed by eight SSR markers. Intrapopulation genetic diversity results suggest that C. dolichocarpa remnants maintained a relatively high degree of genetic diversity despite severe fragmentation. Low genetic differentiation among populations was found based on Wright's F(ST) and amova analysis. Both the F(ST)-based estimate and private allele method revealed high historical gene flow among the remnant populations. Recent immigrants, detected by assignment tests, tend to decrease from the grandparent generation to the current generation. The potentially highly restricted current gene flow among fragments may render the fragmented populations of C. dolichocarpa at a higher risk of local extinction several generations after fragmentation. Both in situ and ex situ conservation management for the remnant populations of C. dolichocarpa are therefore urgently needed to rescue remaining genetic diversity.     

Testing limits to adaptation along altitudinal gradients in rainforest Drosophila

Given that evolution can generate rapid and dramatic shifts in the ecological tolerance of a species, what prevents populations adapting to expand into new habitat at the edge of their distributions? Recent population genetic models have focused on the relative costs and benefits of migration between populations. On the one hand, migration may limit adaptive divergence by preventing local populations from matching their local selective optima. On the other hand, migration may also contribute to the genetic variance necessary to allow populations to track these changing optima. Empirical evidence for these contrasting effects of gene flow in natural situations are lacking, largely because it remains difficult to acquire. Here, we develop a way to explore theoretical models by estimating genetic divergence in traits that confer stress resistance along similar ecological gradients in rainforest Drosophila. This approach allows testing for the coupling of clinal divergence with local density, and the effects of genetic variance and the rate of change of the optimum on the response to selection. In support of a swamping effect of migration on phenotypic divergence, our data show no evidence for a cline in stress-related traits where the altitudinal gradient is steep, but significant clinal divergence where it is shallow. However, where clinal divergence is detected, sites showing trait means closer to the presumed local optimum have more genetic variation than sites with trait means distant from their local optimum. This pattern suggests that gene flow also aids a sustained response to selection.     

Disentangling the roles of history and local selection in shaping clinal variation of allele frequencies and gene expression in Norway spruce (Picea abies)

Understanding the genetic basis of local adaptation is challenging due to the subtle balance among conflicting evolutionary forces that are involved in its establishment and maintenance. One system with which to tease apart these difficulties is clines in adaptive characters. Here we analyzed genetic and phenotypic variation in bud set, a highly heritable and adaptive trait, among 18 populations of Norway spruce (Picea abies), arrayed along a latitudinal gradient ranging from 47°N to 68°N. We confirmed that variation in bud set is strongly clinal, using a subset of five populations. Genotypes for 137 single-nucleotide polymorphisms (SNPs) chosen from 18 candidate genes putatively affecting bud set and 308 control SNPs chosen from 264 random genes were analyzed for patterns of genetic structure and correlation to environment. Population genetic structure was low (F(ST) = 0.05), but latitudinal patterns were apparent among Scandinavian populations. Hence, part of the observed clinal variation should be attributable to population demography. Conditional on patterns of genetic structure, there was enrichment of SNPs within candidate genes for correlations with latitude. Twenty-nine SNPs were also outliers with respect to F(ST). The enrichment for clinal variation at SNPs within candidate genes (i.e., SNPs in PaGI, PaPhyP, PaPhyN, PaPRR7, and PaFTL2) indicated that local selection in the 18 populations, and/or selection in the ancestral populations from which they were recently derived, shaped the observed cline. Validation of these genes using expression studies also revealed that PaFTL2 expression is significantly associated with latitude, thereby confirming the central role played by this gene in the control of phenology in plants.     

LOCAL ADAPTATION MAINTAINS CLINAL VARIATION IN MELANIN‐BASED COLORATION OF EUROPEAN BARN OWLS (TYTO ALBA)

Ecological parameters vary in space, and the resulting heterogeneity of selective forces can drive adaptive population divergence. Clinal variation represents a classical model to study the interplay of gene flow and selection in the dynamics of this local adaptation process. Although geographic variation in phenotypic traits in discrete populations could be remainders of past adaptation, maintenance of adaptive clinal variation requires recurrent selection. Clinal variation in genetically determined traits is generally attributed to adaptation of different genotypes to local conditions along an environmental gradient, although it can as well arise from neutral processes. Here, we investigated whether selection accounts for the strong clinal variation observed in a highly heritable pheomelanin‐based color trait in the European barn owl by comparing spatial differentiation of color and of neutral genes among populations. Barn owl's coloration varies continuously from white in southwestern Europe to reddish‐brown in northeastern Europe. A very low differentiation at neutral genetic markers suggests that substantial gene flow occurs among populations. The persistence of pronounced color differentiation despite this strong gene flow is consistent with the hypothesis that selection is the primary force maintaining color variation among European populations. Therefore, the color cline is most likely the result of local adaptation.     

Landscape-level spatial genetic structure in Quercus acutissima (Fagaceae)

Quercus acutissima (Fagaceae), a deciduous broad-leaved tree, is an important forest element in hillsides of South Korea. We used allozyme loci, Wright's F statistics, and multilocus spatial autocorrelation statistics to examine the distribution of genetic diversity within and among three local populations and the spatial genetic structure at a landscape scale (15 ha, 250 × 600 m) on Oenaro Island, South Korea. Levels of genetic diversity in Q. acutissima populations were comparable to mean values for other oak species. A moderate but significant deficit of heterozygotes (mean F(IS) = 0.069) was detected within local populations and low but significant differentiation was observed among populations (F(ST) = 0.010). Spatial autocorrelation analyses revealed little evidence of significant genetic structure at spatial scales of 100-120 m. The failure to detect genetic structure within populations may be due to intraspecific competition or random mortality among saplings, resulting in extensive thinning within maternal half-sib groups. Alternatively, low genetic differentiation at the landscape scale indicates substantial gene flow among local populations. Although wind-borne pollen may be the primary source of gene flow in Q. acutissima, these results suggest that acorn movement by animals may be more extensive than previously anticipated. Comparison of these genetic data for Oenaro Island with a disturbed isolated inland population suggests that population-to-population differences in internal genetic structure may be influenced by local variation in regeneration environment (e.g., disturbance).     

Natural selection drives clinal life history patterns in the perennial sunflower species, Helianthus maximiliani

In plants, ecologically important life history traits often display clinal patterns of population divergence. Such patterns can provide strong evidence for spatially varying selection across environmental gradients but also may result from nonselective processes, such as genetic drift, population bottlenecks and spatially restricted gene flow. Comparison of population differentiation in quantitative traits (measured as Q(ST) ) with neutral molecular markers (measured as F(ST) ) provides a useful tool for understanding the relative importance of adaptive and nonadaptive processes in the formation and maintenance of clinal variation. Here, we demonstrate the existence of geographic variation in key life history traits in the diploid perennial sunflower species Helianthus maximiliani across a broad latitudinal transect in North America. Strong population differentiation was found for days to flowering, growth rate and multiple size-related traits. Differentiation in these traits greatly exceeds neutral predictions, as determined both by partial Mantel tests and by comparisons of global Q(ST) values with theoretical F(ST) distributions. These findings indicate that clinal variation in these life history traits likely results from local adaptation driven by spatially heterogeneous environments.     

Effect of gene flow on spatial genetic structure in the riparian canopy tree Cercidiphyllum japonicum revealed by microsatellite analysis

Few studies have analyzed pollen and seed movements at local scale, and genetic differentiation among populations covering the geographic distribution range of a species. We carried out such a study on Cercidiphyllum japonicum; a dioecious broad-leaved tree of cool-temperate riparian forest in Japan. We made direct measurement of pollen and seed movements in a site, genetic structure at the local scale, and genetic differentiation between populations covering the Japanese Archipelago. Parentage analysis of seedlings within a 20-ha study site indicated that at least 28.8% of seedlings were fertilized by pollen from trees outside the study site. The average pollination distance within the study site was 129 m, with a maximum of 666 m. The genotypes of 30% of seedlings were incompatible with those of the nearest female tree, and the maximum seed dispersal distance within the study site was over 300 m. Thus, long-distance gene dispersal is common in this species. The correlation between genetic relatedness and spatial distance among adult trees within the population was not significant, indicating an absence of fine-scale genetic structure perhaps caused by high levels of pollen flow and overlapping seed shadows. Six populations sampled throughout the distribution of C. japonicum in Japan showed significant isolation-by-distance but low levels of genetic differentiation (F(ST) = 0.043), also indicating long-distance gene flow in C. japonicum. Long-distance gene flow had a strong influence on the genetic structure at different spatial scales, and contributes to the maintenance of genetic diversity in C. japonicum.     

Unravelling genetics at the top: mountain islands or isolated belts?

Background and Aims

In mountain plant populations, local adaptation has been described as one of the main responses to climate warming, allowing plants to persist under stressful conditions. This is especially the case for marginal populations at their lowest elevation, as they are highly vulnerable. Adequate levels of genetic diversity are required for selection to take place, while high levels of altitudinal gene flow are seen as a major limiting factor potentially precluding local adaptation processes. Thus, a compromise between genetic diversity and gene flow seems necessary to guarantee persistence under oncoming conditions. It is therefore critical to determine if gene flow occurs preferentially between mountains at similar altitudinal belts, promoting local adaptation at the lowest populations, or conversely along altitude within each mountain.

Methods

Microsatellite markers were used to unravel genetic diversity and population structure, inbreeding and gene flow of populations at two nearby altitudinal gradients of Silene ciliata, a Mediterranean high-mountain cushion plant.

Key Results

Genetic diversity and inbreeding coefficients were similar in all populations. Substantial gene flow was found both along altitudinal gradients and horizontally within each elevation belt, although greater values were obtained along altitudinal gradients. Gene flow may be responsible for the homogeneous levels of genetic diversity found among populations. Bayesian cluster analyses also suggested that shifts along altitudinal gradients are the most plausible scenario.

Conclusions

Past population shifts associated with glaciations and interglacial periods in temperate mountains may partially explain current distributions of genetic diversity and population structure. In spite of the predominance of gene flow along the altitudinal gradients, local genetic differentiation of one of the lower populations together with the detection of one outlier locus might support the existence of different selection forces at low altitudes.     

Population differentiation and gene flow within a metapopulation of a threatened tree, Magnolia stellata (Magnoliaceae)

We examined genetic differentiation among eight local populations of a metapopulation of Magnolia stellata using 10 nuclear and three chloroplast microsatellite (nSSR and cpSSR) markers and evaluated the influence of historical gene flow on population differentiation. The coefficient of genetic differentiation among populations for nSSR (F(ST) = 0.053) was less than half that for cpSSR (0.137). An isolation-by-distance pattern was detected for nSSRs, but not cpSSRs. These results suggest that pollen flow, as well as seed dispersal, has significantly reduced genetic differentiation among populations. We also examined patterns of contemporary pollen flow by paternity analysis of seeds from nine seed parents in one of the populations using the nSSR markers and found it to be greatly restricted by the distance between parents. Although most pollen flow occurred within the population, pollen flow from outside the population accounted for 2.5% of the total. When historical and contemporary pollen flows among populations were compared, the levels of pollen flow seem to have declined recently. We conclude that to conserve M. stellata, it is important to preserve the whole population by maintaining its metapopulation structure and the gene flow among its populations.     

Nuclear and chloroplast microsatellites reveal extreme population differentiation and limited gene flow in the Aegean endemic Brassica cretica (Brassicaceae)

Nuclear and chloroplast microsatellite markers were used to study population structure and gene flow among seven Cretan populations of the Aegean endemic plant species Brassica cretica (Brassicaceae). Both nuclear and chloroplast markers revealed exceptionally high levels of population differentiation (overall F(ST)=0.628 and 1.000, respectively) and relatively little within-population diversity (overall H(S)=0.211 and 0.000, respectively). Maximum-likelihood estimates of directional migration rates were low among all pairs of populations (average Nm=0.286). There was no evidence that differences in flower colour between populations had any influence on historical levels of gene flow. In addition, a haplotype network showed that all five chloroplast haplotypes found in the sample were closely related. Together, these results suggest that current patterns of diversification in B. cretica are mainly a result of genetic drift during the last half million years. The main conclusions from the present study are consistent with the prevailing hypothesis that plant diversification in the Aegean region is driven by random rather than adaptive differentiation among isolated populations.     

Reproductive biology and population genetic structure of Kandelia candel (Rhizophoraceae), a viviparous mangrove species

The pollination biology, mating system, and population genetic structure of Kandelia candel were investigated. Field observations on its pollination and reproductive biology suggested that this species is pollinator dependent for fruit set, and bee activities can lead to substantial geitonogamous selfing. Quantitative analysis of the mating system parameters was performed using progeny arrays assayed for allozyme markers. Multilocus outcrossing rates (t(m)) were estimated to be 0.697 ± 0.091 and 0.797 ± 0.062 in two populations. In comparison to other plant species with mixed-mating system, the level of allozyme variation was very low in the 13 populations sampled along the coastlines of Hong Kong. At the species level, the proportion of polymorphic loci was 20%, number of alleles per locus was 1.2, and heterozygosity was 0.0362. The total gene diversity was primarily distributed within populations (H(S )= 0.0339), and the coefficient of genetic differentiation among populations was low (G(ST )= 0.064). This pattern of population genetic structure suggests that gene flow, primarily in the form of water-dispersed seedlings in viviparous mangrove species, is not as limited as previously thought. However, microgeographic pattern in allele frequency at the marker loci could still be detected between the western and eastern coastal populations.     

Comparative organization of chloroplast, mitochondrial and nuclear diversity in plant populations

Plants offer excellent models to investigate how gene flow shapes the organization of genetic diversity. Their three genomes can have different modes of transmission and will hence experience varying levels of gene flow. We have compiled studies of genetic structure based on chloroplast DNA (cpDNA), mitochondrial DNA (mtDNA) and nuclear markers in seed plants. Based on a data set of 183 species belonging to 103 genera and 52 families, we show that the precision of estimates of genetic differentiation (G(ST)) used to infer gene flow is mostly constrained by the sampling of populations. Mode of inheritance appears to have a major effect on G(ST). Maternally inherited genomes experience considerably more subdivision (median value of 0.67) than paternally or biparentally inherited genomes (approximately 0.10). G(ST) at cpDNA and mtDNA markers covary narrowly when both genomes are maternally inherited, whereas G(ST) at paternally and biparentally inherited markers also covary positively but more loosely and G(ST) at maternally inherited markers are largely independent of values based on nuclear markers. A model-based gross estimate suggests that, at the rangewide scale, historical levels of pollen flow are generally at least an order of magnitude larger than levels of seed flow (median of the pollen-to-seed migration ratio: 17) and that pollen and seed gene flow vary independently across species. Finally, we show that measures of subdivision that take into account the degree of similarity between haplotypes (N(ST) or R(ST)) make better use of the information inherent in haplotype data than standard measures based on allele frequencies only.     

Contrasting patterns of pollen and seed flow influence the spatial genetic structure of sweet vernal grass (Anthoxanthum odoratum) populations

The spatial genetic structure of plant populations is determined by a combination of gene flow, genetic drift, and natural selection. Gene flow in most plants can result from either seed or pollen dispersal, but detailed investigations of pollen and seed flow among populations that have diverged following local adaptation are lacking. In this study, we compared pollen and seed flow among 10 populations of sweet vernal grass (Anthoxanthum odoratum) on the Park Grass Experiment. Overall, estimates of genetic differentiation that were based on chloroplast DNA (cpDNA) and, which therefore resulted primarily from seed flow, were lower (average F(ST) = 0.058) than previously published estimates that were based on nuclear DNA (average F(ST) = 0.095). Unlike nuclear DNA, cpDNA showed no pattern of isolation by adaptation; cpDNA differentiation was, however, inversely correlated with the number of additions (nutrients and lime) that each plot had received. We suggest that natural selection is restricting pollen flow among plots, whereas nutrient additions are increasing seed flow and genetic diversity by facilitating the successful germination and growth of immigrant seeds. This study highlights the importance of considering all potential gene flow mechanisms when investigating determinants of spatial genetic structure, and cautions against the widespread assumption that pollen flow is more important than seed flow for population connectivity in wind-pollinated species.     

Population structure of an endangered species living in contrasted habitats: Parnassia palustris (Saxifragaceae)

In endangered species, it is critical to analyse the level at which populations interact (i.e. dispersal) as well as the levels of inbreeding and local adaptation to set up conservation policies. These parameters were investigated in the endangered species Parnassia palustris living in contrasted habitats. We analysed population structure in 14 populations of northern France for isozymes, cpDNA markers and phenotypic traits related to fitness. Within population genetic diversity and inbreeding coefficients were not correlated to population size. Populations seem not to have undergone severe recent bottleneck. Conversely to pollen migration, seed migration seems limited at a regional scale, which could prevent colonization of new sites even if suitable habitats appear. Finally, the habitat type affects neither within-population genetic diversity nor genetic and phenotypic differentiation among populations. Thus, even if unnoticed local adaptation to habitats exists, it does not influence gene flow between populations.     

A preliminary study of genetic variation in populations of Monstera adansonii var. klotzschiana (Araceae) from North-East Brazil, estimated with AFLP molecular markers

BACKGROUND AND AIMS: This study sought genetic evidence of long-term isolation in populations of Monstera adansonii var. klotzschiana (Araceae), a herbaceous, probably outbreeding, humid forest hemi-epiphyte, in the brejo forests of Ceará (north-east Brazil), and clarification of their relationships with populations in Amazonia and the Atlantic forest of Brazil. METHODS: Within-population genetic diversity and between-population dissimilarity were estimated using AFLP molecular markers in 75 individuals from eight populations located in Ceará, the Brazilian Atlantic Forest and Amazonia. KEY RESULTS: The populations showed a clinal pattern of weak genetic differentiation over a large geographical region (F(ST) = 0.1896). A strong correlation between genetic and geographical distance (Mantel test: r = 0.6903, P = 0.002) suggests a historical pattern of isolation by distance. Genetic structure analysis revealed at least two distinct gene pools in the data. The two isolated Ceará populations are significantly different from each other (pairwise Phi(PT) = 0.137, P = 0.003) and as diverse (Nei's gene diversity, average H(e) = 0.1832, 0.1706) as those in the Atlantic and Amazon forest regions. The population in southern Brazil is less diverse (Nei's gene diversity, average H(e) = 0.127) than the rest. The Ceará populations are related to those of the Atlantic forest rather than those from Amazonia (AMOVA, among-groups variation = 11.95 %, P = 0.037). CONCLUSIONS: The gene pools detected within an overall pattern of clinal variation suggest distinct episodes of gene flow, possibly correlated with past humid forest expansions. The Ceará populations show no evidence of erosion of genetic diversity, although this was expected because of their isolation. Their genetic differentiation and relatively high diversity reinforce the importance of conserving the endangered brejo forests.     

Contrasting patterns of genetic variation in the two sympatric geckos Gekko tawaensis and G. japonicus (Reptilia: Squamata) from western Japan,as revealed by allozyme analyses

Allozyme variation in two congeneric sympatric geckos, Gekko tawaensis and G. japonicus, from western Japan was examined. These species show similar densities and spatial arrangements of populations in this region, and their genetic structures are thus expected to have been formed under the influences of comparable geohistorical, environmental, and demographic factors. Results of the analyses, however, revealed strikingly different genetic patterns in the two species. Populations of G. tawaensis invariably showed a remarkably lowered heterozygosity (0-0.017) compared to G. japonicus (0.089-0.124). On the other hand, the genetic heterogeneity among populations is much greater in G. tawaensis (F(ST)=0.726) than in G. japonicus (F(ST)=0.101). The Mantel test failed to detect any significant correlations between log (estimated migration rate) and log (geographic distance) in either species, or between matrices of interpopulation pairwise F(ST) for the two species. These results suggest that, in each species, formation of the current genetic structure in western Japan has been chiefly influenced by stochastic factors, rather than the geohistorical architecture of this region. The high F(ST) and low heterozygosity in G. tawaensis suggest the effects of severe local fragmentation. On the other hand, the relatively low F(ST) and high heterozygosity in G. japonicus imply extensive gene flow among populations. Absence of significant correlations between the estimated migration rate and geographic distance in G. japonicus may suggest that such gene flow is promoted by human-mediated transport of this primarily house-dwelling lizard.     

Landscape Genetic Structure of a Streamside Tree Species Euptelea pleiospermum (Eupteleaceae): Contrasting Roles of River Valley and Mountain Ridge

We used landscape genetics and statistical models to test how landscape features influence connectivity or create barriers to dispersal for a mountain riparian tree species, Euptelea pleiospermum. Young leaves from 1078 individuals belonging to 36 populations at elevations of 900–2000 m along upper reaches of four rivers were genotyped using eight nuclear microsatellite markers. We found no evidence for the unidirectional dispersal hypothesis in E. pleiospermum within each river. The linear dispersal pattern along each river valley is mostly consistent with the “classical metapopulaton” model. Mountain ridges separating rivers were genetic barriers for this wind-pollinated tree species with anemochorous seeds, whereas river valleys provided important corridors for dispersal. Gene flow among populations along elevational gradients within each river prevails over gene flow among populations at similar elevations but from different rivers. This pattern of gene flow is likely to promote elevational range shifts of plant populations and to hinder local adaptation along elevational gradients. This study provides a paradigm to determine which of the two strategies (migration or adaptation) will be adopted by mountain riparian plants under climate warming.     

Local adaptation at the range peripheries of Sitka spruce

High‐dispersal rates in heterogeneous environments and historical rapid range expansion can hamper local adaptation; however, we often see clinal variation in high‐dispersal tree species. To understand the mechanisms of the species’ distribution, we investigated local adaptation and adaptive plasticity in a range‐wide context in Sitka spruce, a wind‐pollinated tree species that has recently expanded its range after glaciations. Phenotypic traits were observed using growth chamber experiments that mimicked temperature and photoperiodic regimes from the limits of the species realized niche. Bud phenology exhibited parallel reaction norms among populations; however, putatively adaptive plasticity and strong divergent selection were seen in bud burst and bud set timing respectively. Natural selection appears to have favoured genotypes that maximize growth rate during available frost‐free periods in each environment. We conclude that Sitka spruce has developed local adaptation and adaptive plasticity throughout its range in response to current climatic conditions despite generally high pollen flow and recent range expansion.     

Genetic diversity of north-east Asian cattle based on microsatellite data

In order to assess the genetic variability and population structure of north-east Asian cattle, 13 microsatellite loci were analysed for a total of 200 individuals including Korean, Chinese, Japanese Black and European Holstein cattle. Observed and expected heterozygosity, two estimators (F(ST) and G(ST)) of gene differentiation, and Nei's DA distance were evaluated. Based on expected mean heterozygosity, the lowest genetic diversity was exhibited in Japanese Black cattle (H(E)=0.471), and the highest in Chinese cattle (H(E)=0.744). Korean cattle revealed a relatively high degree of genetic diversity (H(E)=0.728). Average proportion of genetic variation because of interpopulation subdivision among north-east Asian cattle varied between 10.9 and 9.9%, depending on the estimator used. N-J tree based on Nei's DA genetic distance showed that Korean and Chinese cattle are closely related, whereas Japanese Black cattle are clearly distinct from the other two populations, forming a north-east Asian outgroup.     

The influence of altitude and topography on genetic structure in the long-toed salamander (Ambystoma macrodactulym)

A primary goal of molecular ecology is to understand the influence of abiotic factors on the spatial distribution of genetic variation. Features including altitudinal clines, topography and landscape characteristics affect the proportion of suitable habitat, influence dispersal patterns, and ultimately structure genetic differentiation among populations. We studied the effects of altitude and topography on genetic variation of long-toed salamanders (Ambystoma macrodactylum), a geographically widespread amphibian species throughout northwestern North America. We focused on 10 low altitude sites (< 1200 m) and 11 high-altitude sites in northwestern Montana and determined multilocus genotypes for 549 individuals using seven microsatellite loci. We tested four hypotheses: (1) gene flow is limited between high- and low-altitude sites; and, (2) gene flow is limited among high-altitude sites due to harsh habitat and extreme topographical relief between sites; (3) low-altitude sites exhibit higher among-site gene flow due to frequent flooding events and low altitudinal relief; and (4) there is a negative correlation between altitude and genetic variation. Overall F(ST) values were moderate (0.08611; P < 0.001). Pairwise F(ST) estimates between high and low populations and a population graphing method supported the hypothesis that low-altitude and high-altitude sites, taken together, are genetically differentiated from each other. Also as predicted, gene flow is more prominent among low-altitude sites than high-altitude sites; low-altitude sites had a significantly lower F(ST) (0.03995; P < 0.001) than high altitude sites (F(ST) = 0.10271; P < 0.001). Use of Bayesian analysis of population structure (BAPS) resulted in delineation of 10 genetic groups, two among low-altitude populations and eight among high-altitude populations. In addition, within high altitude populations, basin-level genetic structuring was apparent. A nonequilibrium algorithm for detecting current migration rates supported these population distinctions. Finally, we also found a significant negative correlation between genetic diversity and altitude. These results are consistent with the hypothesis that topography and altitudinal gradients shape the spatial distribution of genetic variation in a species with a broad geographical range and diverse life history. Our study sheds light on which key factors limit dispersal and ultimately species' distributions.