Variation for neutral markers is correlated with variation for quantitative traits in the plant pathogenic fungus Mycosphaerella graminicola

We compared genetic variation and population differentiation at RFLP marker loci with seven quantitative characters including fungicide resistance, temperature sensitivity, pycnidial size, pycnidial density, colony size, percentage of leaves covered by pycnidia (PLACP) and percentage of leaves covered by lesions (PLACL) in Mycosphaerella graminicola populations sampled from four regions. Wide variation in population differentiation was found across the quantitative traits assayed. Fungicide resistance, temperature sensitivity, and PLACP displayed a significantly higher Q(ST) than G(ST), consistent with selection for local adaptation, while pycnidial size, pycnidial density and colony size displayed a lower or significantly lower Q(ST) than G(ST), consistent with constraining selection. There was not a statistical difference between Q(ST) and G(ST) in PLACL. We also found a positive and significant correlation between genetic variation in molecular marker loci and quantitative traits at the multitrait scale, suggesting that estimates of overall genetic variation for quantitative traits in M. graminicola could be derived from analysis of the molecular genetic markers.     

Variation at two flowering time genes within and among populations of Arabidopsis thaliana: comparison with markers and traits

Flowering Locus C (FLC) and Frigida are two interacting genes controlling flowering time variation in Arabidopsis thaliana. Variation at these genes was surveyed in 12 A. thaliana populations sampled in France. These populations were also screened for variation at molecular markers [12 microsatellites and 19 cleaved amplified polymorphic sequence (CAPS) markers] and at seven quantitative traits measured with and without vernalization. Seven populations were highly polymorphic at markers (H(S) = 0.57 at microsatellites, 0.24 at CAPS) and showed heritable variation for bolting time and some other traits. Five populations were genetically fixed or nearly fixed. Q(ST) for bolting time without vernalization was significantly higher than F(ST), suggesting local divergent selection. One of the two haplotype groups at FLC (FLC(A)) was very predominant (frequency of 99%). The first exon of Frigida showed elevated nonsynonymous variation, and nine loss-of-function mutations were found throughout the gene. The association between loss-of-function and earlier bolting was confirmed. Overall, 18 Frigida haplotypes were detected. The pattern of variation at Frigida was largely similar to that found at markers and traits, with the same populations being fixed or highly diverse. Metapopulation dynamics is thus probably the main factor shaping genetic variation in A. thaliana. However, F(ST) for functional (FRI) vs. nonfunctional (FRI(Delta)) haplotypes was significantly higher than F(ST) at markers. This suggested that loss-of-function at Frigida is under local selection for flowering time.     

Distinguishing adaptive from nonadaptive genetic differentiation: comparison of Q(ST) and F(ST) at two spatial scales

Genetic differentiation in 20 hierarchically sampled populations of wild barley was analyzed with quantitative traits, allozymes and Random Amplified Polymorphic DNAs (RAPDs), and compared for three marker types at two hierarchical levels. Regional subdivision for both molecular markers was much lower than for quantitative traits. For both allozymes and RAPDs, most loci exhibited minor or no regional differentiation, and the relatively high overall estimates of the latter were due to several loci with exceptionally high regional differentiation. The allozyme- and RAPD-specific patterns of differentiation were concordant in general with one another, but not with quantitative trait differentiation. Divergent selection on quantitative traits inferred from very high regional Q(ST) was in full agreement with our previous results obtained from a test of local adaptation and multilevel selection analysis. In contrast, most variation in allozyme and RAPD variation was neutral, although several allozyme loci and RAPD markers were exceptional in their levels of regional differentiation. However, it is not possible to answer the question whether these exceptional loci are directly involved in the response to selection pressure or merely linked to the selected loci. The fact that Q(ST) and F(ST) did not differ at the population scale, that is, within regions, but differed at the regional scale, for which local adaptation has been previously shown, implies that comparison of the level of subdivision in quantitative traits, as compared with molecular markers, is indicative of adaptive population differentiation only when sampling is carried out at the appropriate scale.     

Does habitat fragmentation reduce fitness and adaptability? A case study of the common frog (Rana temporaria)

Studies examining the effects of anthropogenic habitat fragmentation on both neutral and adaptive genetic variability are still scarce. We compared tadpole fitness-related traits (viz. survival probability and body size) among populations of the common frog (Rana temporaria) from fragmented (F) and continuous (C) habitats that differed significantly in population sizes (C > F) and genetic diversity (C > F) in neutral genetic markers. Using data from common garden experiments, we found a significant positive relationship between the mean values of the fitness related traits and the amount of microsatellite variation in a given population. While genetic differentiation in neutral marker loci (F(ST)) tended to be more pronounced in the fragmented than in the continuous habitat, genetic differentiation in quantitative traits (Q(ST)) exceeded that in neutral marker traits in the continuous habitat (i.e. Q(ST) > F(ST)), but not in the fragmented habitat (i.e. Q(ST) approximately F(ST)). These results suggest that the impact of random genetic drift relative to natural selection was higher in the fragmented landscape where populations were small, and had lower genetic diversity and fitness as compared to populations in the more continuous landscape. The findings highlight the potential importance of habitat fragmentation in impairing future adaptive potential of natural populations.     

Maternal genetic effects on adaptive divergence between anadromous and resident brook charr during early life history

The importance of directional selection relative to neutral evolution may be determined by comparing quantitative genetic variation in phenotype (Q(ST)) to variation at neutral molecular markers (F(ST)). Quantitative divergence between salmonid life history types is often considerable, but ontogenetic changes in the significance of major sources of genetic variance during post-hatch development suggest that selective differentiation varies by developmental stage. In this study, we tested the hypothesis that maternal genetic differentiation between anadromous and resident brook charr (Salvelinus fontinalis Mitchill) populations for early quantitative traits (embryonic size/growth, survival, egg number and developmental time) would be greater than neutral genetic differentiation, but that the maternal genetic basis for differentiation would be higher for pre-resorption traits than post-resorption traits. Quantitative genetic divergence between anadromous (seawater migratory) and resident Laval River (Québec) brook charr based on maternal genetic variance was high (Q(ST) > 0.4) for embryonic length, yolk sac volume, embryonic growth rate and time to first response to feeding relative to neutral genetic differentiation [F(ST) = 0.153 (0.071-0.214)], with anadromous females having positive genetic coefficients for all of the above characters. However, Q(ST) was essentially zero for all traits post-resorption of the yolk sac. Our results indicate that the observed divergence between resident and anadromous brook charr has been driven by directional selection, and may therefore be adaptive. Moreover, they provide among the first evidence that the relative importance of selective differentiation may be highly context-specific, and varies by genetic contributions to phenotype by parental sex at specific points in offspring ontogeny. This in turn suggests that interpretations of Q(ST)-F(ST) comparisons may be improved by considering the structure of quantitative genetic architecture by age category and the sex of the parent used in estimation.     

Comparisons between QST and FST—how wrong have we been?

The comparison between quantitative genetic divergence (Q(ST) ) and neutral genetic divergence (F(ST) ) among populations has become the standard test for historical signatures of selection on quantitative traits. However, when the mutation rate of neutral markers is relatively high in comparison with gene flow, estimates of F(ST) will decrease, resulting in upwardly biased comparisons of Q(ST) vs. F(ST) . Reviewing empirical studies, the difference between Q(ST) and F(ST) is positively related to marker heterozygosity. After refuting alternative explanations for this pattern, we conclude that marker mutation rate indeed has had a biasing effect on published Q(ST) -F(ST) comparisons. Hence, it is no longer clear that populations have commonly diverged in response to divergent selection. We present and discuss potential solutions to this bias. Comparing Q(ST) with recent indices of neutral divergence that statistically correct for marker heterozygosity (Hedrick's G'st and Jost's D) is not advised, because these indices are not theoretically equivalent to Q(ST) . One valid solution is to estimate F(ST) from neutral markers with mutation rates comparable to those of the loci underlying quantitative traits (e.g. SNPs). Q(ST) can also be compared to Φ(ST) (Phi(ST) ) of amova, as long as the genetic distance among allelic variants used to estimate Φ(ST) reflects evolutionary history: in that case, neutral divergence is independent of mutation rate. In contrast to their common usage in comparisons of Q(ST) and F(ST) , microsatellites typically have high mutation rates and do not evolve according to a simple evolutionary model, so are best avoided in Q(ST) -F(ST) comparisons.     

Q(St) meets the G matrix: the dimensionality of adaptive divergence in multiple correlated quantitative traits.

The Q(ST)-F(ST) comparison has become an increasingly common method for inferring adaptive quantitative trait divergence among populations. For cases in which there is divergence in multiple traits, most studies have applied the method by performing multiple univariate Q(ST)-F(ST) comparisons. However, because traits are often genetically correlated, such univariate analyses are likely to paint a simplified picture of adaptive divergence. Here we show how the multivariate analogue of Q(ST), F(STq), which accounts for genetic correlations among traits, can be used to supply a more detailed picture of multitrait divergence. We apply the method to naturally occurring genetic variation for a suite of sexually selected display traits in Drosophila serrata. The analyses suggest the operation of divergent multivariate selection that has influenced multiple independent axes of genetic variance in a sex-specific manner. Finally, we show how a comparison of the components of F(STq), the average within and among population genetic variance-covariance matrices, G(W) and G(B), can be used as an additional test of the null expectation of neutral divergence, and allows for an investigation of whether natural populations have diverged along major or minor axes of genetic variance.     

Latitudinal divergence of common frog (Rana temporaria) life history traits by natural selection: evidence from a comparison of molecular and quantitative genetic data

The relative roles of natural selection and direct environmental induction, as well as of natural selection and genetic drift, in creating clinal latitudinal variation in quantitative traits have seldom been assessed in vertebrates. To address these issues, we compared molecular and quantitative genetic differentiation between six common frog (Rana temporaria) populations along an approximately 1600 km long latitudinal gradient across Scandinavia. The degree of population differentiation (QST approximately 0.81) in three heritable quantitative traits (age and size at metamorphosis, growth rate) exceeded that in eight (neutral) microsatellite loci (FST = 0.24). Isolation by distance was clear for both neutral markers and quantitative traits, but considerably stronger for one of the three quantitative traits than for neutral markers. QST estimates obtained using animals subjected to different rearing conditions (temperature and food treatments) revealed some environmental dependency in patterns of population divergence in quantitative traits, but in general, these effects were weak in comparison to overall patterns. Pairwise comparisons of FST and QST estimates across populations and treatments revealed that the degree of quantitative trait differentiation was not generally predictable from knowledge of that in molecular markers. In fact, both positive and negative correlations were observed depending on conditions where the quantitative genetic variability had been measured. All in all, the results suggest a very high degree of genetic subdivision both in neutral marker genes and genes coding quantitative traits across a relatively recently (< 9000 years) colonized environmental gradient. In particular, they give evidence for natural selection being the primary agent behind the observed latitudinal differentiation in quantitative traits.     

Genetic structure and mating system of Euterpe edulis Mart. Populations: a comparative analysis using microsatellite and allozyme markers

A comparative study between microsatellite and allozyme markers was conducted on the genetic structure and mating system in natural populations of Euterpe edulis Mart. Three cohorts, including seedlings, saplings, and adults, were examined in 4 populations using 10 allozyme loci and 10 microsatellite loci. As expected, microsatellite markers had a much higher degree of polymorphism than allozymes, but estimates of multilocus outcrossing rate ( = 1.00), as well as estimates of genetic structure (F(IS), G(ST)), were similar for the 2 sets of markers. Estimates of R(ST), for microsatellites, were higher than those of G(ST), but results of both statistics revealed a close agreement for the genetic structure of the species. This study provides support for the important conclusion that allozymes are still useful and reliable markers to estimate population genetic parameters. Effects of sample size on estimates from hypervariable loci are also discussed in this paper.     

Comparative studies of quantitative trait and neutral marker divergence: a meta-analysis

Comparative studies of quantitative genetic and neutral marker differentiation have provided means for assessing the relative roles of natural selection and random genetic drift in explaining among-population divergence. This information can be useful for our fundamental understanding of population differentiation, as well as for identifying management units in conservation biology. Here, we provide comprehensive review and meta-analysis of the empirical studies that have compared quantitative genetic (Q(ST)) and neutral marker (F(ST)) differentiation among natural populations. Our analyses confirm the conclusion from previous reviews - based on ca. 100% more data - that the Q(ST) values are on average higher than F(ST) values [mean difference 0.12 (SD 0.27)] suggesting a predominant role for natural selection as a cause of differentiation in quantitative traits. However, although the influence of trait (life history, morphological and behavioural) and marker type (e.g. microsatellites and allozymes) on the variance of the difference between Q(ST) and F(ST) is small, there is much heterogeneity in the data attributable to variation between specific studies and traits. The latter is understandable as there is no reason to expect that natural selection would be acting in similar fashion on all populations and traits (except for fitness itself). We also found evidence to suggest that Q(ST) and F(ST) values across studies are positively correlated, but the significance of this finding remains unclear. We discuss these results in the context of utility of the Q(ST)-F(ST) comparisons as a tool for inferring natural selection, as well as associated methodological and interpretational problems involved with individual and meta-analytic studies.     

Population Structure in Daphnia Obtusa: Quantitative Genetic and Allozymic Variation

Quantitative genetic analyses for body size and for life history characters within and among populations of Daphnia obtusa reveal substantial genetic variance at both hierarchical levels for all traits measured. Simultaneous allozymic analysis on the same population samples indicate a moderate degree of differentiation: G(ST) = 0.28. No associations between electrophoretic genotype and phenotypic characters were found, providing support for the null hypothesis that the allozymic variants are effectively neutral. Therefore, G(ST) can be used as the null hypothesis that neutral phenotypic evolution within populations led to the observed differentiation for the quantitative traits, which I call Q(ST). The results of this study provide evidence that natural selection has promoted diversification for body size among populations, and has impeded diversification for relative fitness. Analyses of population differentiation for clutch size, age at reproduction, and growth rate indicate that neutral phenotypic evolution cannot be excluded as the cause.     

Divergent selection as revealed by P(ST) and QTL-based F(ST) in three-spined stickleback (Gasterosteus aculeatus) populations along a coastal-inland gradient

Three measures of divergence, estimated at nine putatively neutral microsatellite markers, 14 quantitative traits, and seven quantitative trait loci (QTL) were compared in eight populations of the three-spined stickleback (Gasterosteus aculeatus L.) living in the Scheldt river basin (Belgium). Lowland estuarine and polder populations were polymorphic for the number of lateral plates, whereas upland freshwater populations were low-plated. The number of short gill rakers and the length of dorsal and pelvic spines gradually declined along a coastal-inland gradient. Plate number, short gill rakers and spine length showed moderate to strong signals of divergent selection between lowland and upland populations in comparison between P(ST) (a phenotypic alternative for Q(ST)) and neutral F(ST). However, such comparisons rely on the unrealistic assumption that phenotypic variance equals additive genetic variance, and that nonadditive genetic effects and environmental effects can be minimized. In order to verify this assumption and to confirm the phenotypic signals of divergence, we tested for divergent selection at the underlying QTL. For plate number, strong genetic evidence for divergent selection between lowland and upland populations was obtained based on an intron marker of the Eda gene, of which the genotype was highly congruent with plate morph. Genetic evidence for divergent selection on short gill rakers was limited to some population pairs where F(ST) at only one of two QTL was detected as an outlier, although F(ST) at both loci correlated significantly with P(ST). No genetic confirmation was obtained for divergent selection on dorsal spine length, as no outlier F(ST)s were detected at dorsal spine QTL, and no significant correlations with P(ST) were observed.     

Molecular and quantitative genetic differentiation across Europe in yellow dung flies

Relating geographic variation in quantitative traits to underlying population structure is crucial for understanding processes driving population differentiation, isolation and ultimately speciation. Our study represents a comprehensive population genetic survey of the yellow dung fly Scathophaga stercoraria, an important model organism for evolutionary and ecological studies, over a broad geographic scale across Europe (10 populations from the Swiss Alps to Iceland). We simultaneously assessed differentiation in five quantitative traits (body size, development time, growth rate, proportion of diapausing individuals and duration of diapause), to compare differentiation in neutral marker loci (F(ST)) to that of quantitative traits (Q(ST)). Despite long distances and uninhabitable areas between sampled populations, population structuring was very low but significant (F(ST) = 0.007, 13 microsatellite markers; F(ST) = 0.012, three allozyme markers; F(ST) = 0.007, markers combined). However, only two populations (Iceland and Sweden) showed significant allelic differentiation to all other populations. We estimated high levels of gene flow [effective number of migrants (Nm) = 6.2], there was no isolation by distance, and no indication of past genetic bottlenecks (i.e. founder events) and associated loss of genetic diversity in any northern or island population. In contrast to the low population structure, quantitative traits were strongly genetically differentiated among populations, following latitudinal clines, suggesting that selection is responsible for life history differentiation in yellow dung flies across Europe.     

If FST does not measure neutral genetic differentiation,then comparing it with QST is misleading. Or is it?

The comparison between neutral genetic differentiation (F(ST) ) and quantitative genetic differentiation (Q(ST) ) is commonly used to test for signatures of selection in population divergence. However, there is an ongoing discussion about what F(ST) actually measures, even resulting in some alternative metrics to express neutral genetic differentiation. If there is a problem with F(ST) , this could have repercussions for its comparison with Q(ST) as well. We show that as the mutation rate of the neutral marker increases, F(ST) decreases: a higher within-population heterozygosity (He) yields a lower F(ST) value. However, the same is true for Q(ST) : a higher mutation rate for the underlying QTL also results in a lower Q(ST) estimate. The effect of mutation rate is equivalent in Q(ST) and F(ST) . Hence, the comparison between Q(ST) and F(ST) remains valid, if one uses neutral markers whose mutation rates are not too high compared to those of quantitative traits. Usage of highly variable neutral markers such as hypervariable microsatellites can lead to serious biases and the incorrect inference that divergent selection has acted on populations. Much of the discussion on F(ST) seems to stem from the misunderstanding that it measures the differentiation of populations, whereas it actually measures the fixation of alleles. In their capacity as measures of population differentiation, Hedrick's G'(ST) and Jost's D reach their maximum value of 1 when populations do not share alleles even when there remains variation within populations, which invalidates them for comparisons with Q(ST) .     

Evolutionary inference from QST

Q(ST) is a commonly used metric of the degree of genetic differentiation among populations displayed by quantitative traits. Typically, Q(ST) is compared to F(ST) measured on putatively neutral loci; if Q(ST)=F(ST), this is taken as evidence of spatially heterogeneous and diversifying selection. This paper reviews the uses, assumptions and statistics of Q(ST) and F(ST) comparisons. Unfortunately, Q(ST)/F(ST) comparisons are statistically challenging. For a single trait, Q(ST) must be compared not to the mean F(ST) but to the distribution of F(ST) values. The sources of biases and sampling error for Q(ST) are reviewed, and a new method for comparing Q(ST) and F(ST) is suggested. Simulation results suggest that the distribution of neutral F(ST) and Q(ST) values are little affected by various deviations from the island model. Consequently, the distributions of Q(ST) and F(ST) are well approximated by the Lewontin-Krakauer prediction, even with realistic deviations from the island-model assumptions.     

Comparative Assessment of Genetic Variability in the Populations of Endemic and Endangered Yellow Catfish, Horabagrus brachysoma (Teleostei: Horabagridae), Based on Allozyme, RAPD, and Microsatellite Markers

The comparative assessment of genetic diversity using allozymes, random amplified polymorphic DNA (RAPD), and microsatellite markers was conducted in endemic and endangered yellow catfish (Horabagrus brachysoma) sampled from three locations in Western Ghats river systems of India. Among the three markers, microsatellites show more polymorphism, having 100% polymorphic loci, whereas allozymes show the least (56%). In RAPD, 60.5% of fragments were polymorphic. Observed heterozygosity and F(ST) values were very high in microsatellites, compared with the other markers. Microsatellite and RAPD markers reported a higher degree of genetic differentiation than allozymes among the populations depicted by pairwise F(ST)/G(ST), AMOVA, Nei's genetic distance, and UPGMA dendrogram. The three classes of markers demonstrated striking genetic differentiation between pairs of H. brachysoma populations. The data emphasize the need for fishery management, conservation, and rehabilitation of this species.     

Genetic isolation of fragmented populations is exacerbated by drift and selection

Reduced genetic variation at marker loci in small populations has been well documented, whereas the relationship between quantitative genetic variation and population size has attracted little empirical investigation. Here we demonstrate that both neutral and quantitative genetic variation are reduced in small populations of a fragmented plant metapopulation, and that both drift and selective change are enhanced in small populations. Measures of neutral genetic differentiation (F(ST)) and quantitative genetic differentiation (Q(ST)) in two traits were higher among small demes, and Q(ST) between small populations exceeded that expected from drift alone. This suggests that fragmented populations experience both enhanced genetic drift and divergent selection on phenotypic traits, and that drift affects variation in both neutral markers and quantitative traits. These results highlight the need to integrate natural selection into conservation genetic theory, and suggests that small populations may represent reservoirs of genetic variation adaptive within a wide range of environments.     

Adaptive genetic differentiation in a predominantly self-pollinating species analyzed by transplanting into natural environment, crossbreeding and Q(ST)-F(ST) test

? Both genetic drift and natural selection result in genetic/phenotypic differentiation over space. I analyzed the role of local adaptation in the genetic differentiation of populations of the annual grass Hordeum spontaneum sampled along an aridity gradient. ? The study included the introduction of plants having desert vs nondesert origin into natural (desert) environment, analysis of population differentiation in allozymes and random amplified polymorphic DNA (RAPD) markers vs phenotypic traits (Q(ST) -F(ST) comparison), and planting interpopulation hybrids under simulated desert conditions in a glasshouse. ? The results of the home advantage test, Q(ST) -F(ST) comparison and crossbreeding were consistent with local adaptation; that is, that differentiation of the desert plants from plants of nondesert origin in phenotypic traits was adaptive, giving them home advantage. Each method used provided additional, otherwise unavailable, information, meaning that they should be viewed as complementary rather than alternative approaches. ? Gene flow from adjacent populations (i.e. populations experiencing the desert environment) via seeds (but not pollen) had a positive effect on fitness by enhancing natural selection and counteracting drift. At the same time, the effect of genes from the species distributional core (nondesert plants) by either seed or pollen had a negative fitness effect despite its enriching effect on neutral diversity. The pattern of outbreeding depression observed in interpopulation hybrids (F(1) ) and their segregating progeny (F(2) ) was inconsistent with underdominance, but indicated the presence of additive, dominance and epistatic effects.     

Population genetic structure in a Mediterranean pine (Pinus pinaster Ait.): a comparison of allozyme markers and quantitative traits

F-statistics were employed to analyse quantitative and allozyme variation among 19 native populations of maritime pine (Pinus pinaster Ait.). Fourteen polymorphic allozyme loci were used to provide an empirical basis for constructing a null hypothesis to test natural selection as a determinant of quantitative evolution in stem form, total height growth and survival at 30 years old. Hidden biases, that may result in a difference between quantitative (Q(ST)) and allozyme (F(ST)) differentiation which are not because of the action of natural selection, were avoided by comparing pairs of populations using linear models. All quantitative traits showed higher differentiation than allozymes. The highest divergence was found in stem form, whereas divergences in total height and survival were significantly lower. Differential adaptation to regional and local patterns of precipitation, temperature and soil type seem to be the best explanation of the different structure found in quantitative traits and allozyme loci. Possible bias in the estimation of Q(ST) due to the level of quantitative within-population diversity and the role of adaptation of maritime pine after the last glaciation to highly diverse ecological conditions are discussed with special reference to the actual geographical structure of gene diversity in the species' native range.     

Decoupled mitochondrial and chloroplast DNA population structure reveals Holocene collapse and population isolation in a threatened Mexican-endemic conifer

Chihuahua spruce (Picea chihuahuana Martínez) is a montane subtropical conifer endemic to the Sierra Madre Occidental in northwestern México. Range-wide variation was investigated using maternally inherited mitochondrial (mtDNA) and paternally inherited chloroplast (cpDNA) DNA markers. Among the 16 mtDNA regions analysed, only two mitotypes were detected, while the study of six cpDNA microsatellite markers revealed eight different chlorotypes. The average cpDNA diversity (H = 0.415) was low but much higher than that for mtDNA (H = 0). The distribution of mitotypes revealed two clear nonoverlapping areas (G(ST) = N(ST) = 1), one including northern populations and the second one including the southern and central stands, suggesting that these two regions may represent different ancestral populations. The cpDNA markers showed lower population differentiation (G(ST) = 0.362; R(ST) = 0.230), implying that the two ancestral populations continued to exchange pollen after their initial geographic separation. A lack of a phylogeographic structure was revealed by different spatial analyses of cpDNA (G(ST) > R(ST); and samova), and reduced cpDNA gene flow was noted among populations (Nm = 0.873). Some stands deviated significantly from the mutation-drift equilibrium, suggesting recent bottlenecks. Altogether, these various trends are consistent with the hypothesis of a population collapse during the Holocene warming and suggest that most of the modern P. chihuahuana populations are now effectively isolated with their genetic diversity essentially modelled by genetic drift. The conservation efforts should focus on most southern populations and on the northern and central stands exhibiting high levels of genetic diversity. Additional mtDNA sequence analysis confirmed that P. martinezii (Patterson) is not conspecific with P. chihuahuana, and thus deserves separate conservation efforts.