Molecular and quantitative trait variation within and among small fragmented populations of the endangered plant species Psilopeganum sinense

Background and Aims

Natural selection and genetic drift are important evolutionary forces in determining genetic and phenotypic differentiation in plant populations. The extent to which these two distinct evolutionary forces affect locally adaptive quantitative traits has been well studied in common plant and animal species. However, we know less about how quantitative traits respond to selection pressures and drift in endangered species that have small population sizes and fragmented distributions. To address this question, this study assessed the relative strengths of selection and genetic drift in shaping population differentiation of phenotypic traits in Psilopeganum sinense, a naturally rare and recently endangered plant species.

Methods

Population differentiation at five quantitative traits (Q_ST) obtained from a common garden experiment was compared with differentiation at putatively neutral microsatellite markers (F_ST) in seven populations of P. sinense. Q_ST estimates were derived using a Bayesian hierarchical variance component method.

Key Results

Trait-specific Q_ST values were equal to or lower than F_ST. Neutral genetic diversity was not correlated with quantitative genetic variation within the populations of P. sinense.

Conclusions

Despite the prevalent empirical evidence for Q_ST > F_ST, the results instead suggest a definitive role of stabilizing selection and drift leading to phenotypic differentiation among small populations. Three traits exhibited a significantly lower Q_ST relative to F_ST, suggesting that populations of P. sinense might have experienced stabilizing selection for the same optimal phenotypes despite large geographical distances between populations and habitat fragmentation. For the other two traits, Q_ST estimates were of the same magnitude as F_ST, indicating that divergence in these traits could have been achieved by genetic drift alone. The lack of correlation between molecular marker and quantitative genetic variation suggests that sophisticated considerations are required for the inference of conservation measures of P. sinense from neutral genetic markers.     

Adaptive differentiation in floral traits in the presence of high gene flow in scarlet gilia (Ipomopsis aggregata)

Plant–pollinator interactions are thought to be major drivers of floral trait diversity. However, the relative importance of divergent pollinator‐mediated selection vs. neutral processes in floral character evolution has rarely been explored. We tested for adaptive floral trait evolution by comparing differentiation at neutral genetic loci to differentiation at quantitative floral traits in a putative Ipomopsis aggregata hybrid zone. Typical I. aggregata subsp. candida displays slender white tubular flowers that are typical of flowers pollinated by hawkmoths, and subsp. collina displays robust red tubular flowers typical of flowers pollinated by hummingbirds; yet, hybrid flower morphs are abundant across the East Slope of the Colorado Rockies. We estimated genetic differentiation (F_ST) for nuclear and chloroplast microsatellite loci and used a half‐sib design to calculate quantitative trait divergence (Q_ST) from collection sites across the morphological hybrid zone. We found little evidence for population structure and estimated mean F_ST to be 0.032. Q_ST values for several floral traits including corolla tube length and width, colour, and nectar volume were large and significantly greater than mean F_ST. We performed multivariate comparisons of neutral loci to genetic correlations within and between populations and found a strong signal for divergent selection, suggesting that specific combinations of floral display and reward traits may be the targets of selection. Our results show little support for historical subspecies categories, yet floral traits are more diverged than expected due to drift alone. Non‐neutral divergence for multivariate quantitative traits suggests that selection by pollinators is maintaining a correlation between display and reward traits.     

F(ST) and Q(ST) under neutrality

A commonly used test for natural selection has been to compare population differentiation for neutral molecular loci estimated by F_ST and for the additive genetic component of quantitative traits estimated by Q_ST. Past analytical and empirical studies have led to the conclusion that when averaged over replicate evolutionary histories, Q_ST = F_ST under neutrality. We used analytical and simulation techniques to study the impact of stochastic fluctuation among replicate outcomes of an evolutionary process, or the evolutionary variance, of Q_ST and F_ST for a neutral quantitative trait determined by n unlinked diallelic loci with additive gene action. We studied analytical models of two scenarios. In one, a pair of demes has recently been formed through subdivision of a panmictic population; in the other, a pair of demes has been evolving in allopatry for a long time. A rigorous analysis of these two models showed that in general, it is not necessarily true that mean Q_ST = F_ST (across evolutionary replicates) for a neutral, additive quantitative trait. In addition, we used finite-island model simulations to show there is a strong positive correlation between Q_ST and the difference Q_ST − F_ST because the evolutionary variance of Q_ST is much larger than that of F_ST. If traits with relatively large Q_ST values are preferentially sampled for study, the difference between Q_ST and F_ST will also be large and positive because of this correlation. Many recent studies have used tests of the null hypothesis Q_ST = F_ST to identify diversifying or uniform selection among subpopulations for quantitative traits. Our findings suggest that the distributions of Q_ST and F_ST under the null hypothesis of neutrality will depend on species-specific biology such as the number of subpopulations and the history of subpopulation divergence. In addition, the manner in which researchers select quantitative traits for study may introduce bias into the tests. As a result, researchers must be cautious before concluding that selection is occurring when Q_ST ≠ F_ST.     

Comparison of genetic differentiation at marker loci and quantitative traits

The comparison of the degree of differentiation in neutral marker loci and genes coding quantitative traits with standardized and equivalent measures of genetic differentiation (F_ST and Q_ST, respectively) can provide insights into two important but seldom explored questions in evolutionary genetics: (i) what is the relative importance of random genetic drift and directional natural selection as causes of population differentiation in quantitative traits, and (ii) does the degree of divergence in neutral marker loci predict the degree of divergence in genes coding quantitative traits? Examination of data from 18 independent studies of plants and animals using both standard statistical and meta‐analytical methods revealed a number of interesting points. First, the degree of differentiation in quantitative traits (Q_ST) typically exceeds that observed in neutral marker genes (F_ST), suggesting a prominent role for natural selection in accounting for patterns of quantitative trait differentiation among contemporary populations. Second, the F_ST – Q_ST difference is more pronounced for allozyme markers and morphological traits, than for other kinds of molecular markers and life‐history traits. Third, very few studies reveal situations were Q_ST < F_ST, suggesting that selection pressures, and hence optimal phenotypes, in different populations of the same species are unlikely to be often similar. Fourth, there is a strong correlation between Q_ST and F_ST indices across the different studies for allozyme (r=0.81), microsatellite (r=0.87) and combined (r=0.75) marker data, suggesting that the degree of genetic differentiation in neutral marker loci is closely predictive of the degree of differentiation in loci coding quantitative traits. However, these interpretations are subject to a number of assumptions about the data and methods used to derive the estimates of population differentiation in the two sets of traits.     

Selection on life‐history traits and genetic population divergence in rotifers

A combination of founder effects and local adaptation – the Monopolization hypothesis – has been proposed to reconcile the strong population differentiation of zooplankton dwelling in ponds and lakes and their high dispersal abilities. The role genetic drift plays in genetic differentiation of zooplankton is well documented, but the impact of natural selection has received less attention. Here, we compare differentiation in neutral genetic markers (F_ST) and in quantitative traits (Q_ST) in six natural populations of the rotifer Brachionus plicatilis to assess the importance of natural selection in explaining genetic differentiation of life‐history traits. Five life‐history traits were measured in four temperature × salinity combinations in common‐garden experiments. Population differentiation for neutral genetic markers – 11 microsatellite loci – was very high (F_ST = 0.482). Differentiation in life‐history traits was higher in traits related to sexual reproduction than in those related to asexual reproduction. Q_ST values for diapausing egg production (a trait related to sexual reproduction) were higher than their corresponding F_ST in some pairs of populations. Our results indicate the importance of divergent natural selection in these populations and suggest local adaptation to the unpredictability of B. plicatilis habitats.     

Divergence at neutral and non-neutral loci in Drosophila buzzatii populations and their hybrids

The impact of intraspecific hybridisation on fitness and morphological traits depends on the history of natural selection and genetic drift, which may have led to differently coadapted gene-complexes in the parental populations. The divergence at neutral and non-neutral loci between populations can be evaluated by estimating F_ST and Q_ST respectively, and hence give an estimate of drift and selection in the populations. Here we investigate (1) whether divergence between populations in quantitative traits (wing size and shape) can be attributed to selection or drift alone, (2) The impact of intraspecific hybridisation on estimators for divergence at neutral (F_ST) and non-neutral loci (Q_ST) in hybrids, (3) If measurement of shape is more informative than size in order to detect divergence in quantitative traits between populations. The aims were addressed by performing two hybridisations between three populations of Drosophila buzzatii, one between populations from Argentina and the Canary Islands (separated for 200 years), and the other between populations from Argentina and Australia (separated for 80 years). We observed the highest divergence at neutral loci between the Argentinean and Canary Island populations, but highest morphological divergence between the Argentinean and Australian populations, indicating that natural selection is acting on the wings. Divergence based on Q_ST measures in the hybrids was sensitive towards increased phenotypic variance (σ²p) within groups and should be used with care when σ²p of populations differ. Our results indicate that measures of shape give a better estimate of divergence at the underlying quantitative traits loci than measures of size.     

The Influence of Genetic Drift and Selection on Quantitative Traits in a Plant Pathogenic Fungus

Genetic drift and selection are ubiquitous evolutionary forces acting to shape genetic variation in populations. While their relative importance has been well studied in plants and animals, less is known about their relative importance in fungal pathogens. Because agro-ecosystems are more homogeneous environments than natural ecosystems, stabilizing selection may play a stronger role than genetic drift or diversifying selection in shaping genetic variation among populations of fungal pathogens in agro-ecosystems. We tested this hypothesis by conducting a Q _ST/F _ST analysis using agricultural populations of the barley pathogen Rhynchosporium commune. Population divergence for eight quantitative traits (Q _ST) was compared with divergence at eight neutral microsatellite loci (F _ST) for 126 pathogen strains originating from nine globally distributed field populations to infer the effects of genetic drift and types of selection acting on each trait. Our analyses indicated that five of the eight traits had Q _ST values significantly lower than F _ST, consistent with stabilizing selection, whereas one trait, growth under heat stress (22°C), showed evidence of diversifying selection and local adaptation (Q _ST>F _ST). Estimates of heritability were high for all traits (means ranging between 0.55–0.84), and average heritability across traits was negatively correlated with microsatellite gene diversity. Some trait pairs were genetically correlated and there was significant evidence for a trade-off between spore size and spore number, and between melanization and growth under benign temperature. Our findings indicate that many ecologically and agriculturally important traits are under stabilizing selection in R. commune and that high within-population genetic variation is maintained for these traits.     

Genetic and Morphological Divergence in Three Strains of Brook Trout Salvelinus fontinalis Commonly Stocked in Lake Superior

Fitness related traits often show spatial variation across populations of widely distributed species. Comparisons of genetic variation among populations in putatively neutral DNA markers and in phenotypic traits susceptible to selection (Q_ST F_ST analysis) can be used to determine to what degree differentiation among populations can be attributed to selection or genetic drift. Traditionally, Q_ST F_ST analyses require a large number of populations to achieve sufficient statistical power; however, new methods have been developed that allow Q_ST F_ST comparisons to be conducted on as few as two populations if their pedigrees are informative. This study compared genetic and morphological divergence in three strains of brook trout Salvelinus fontinalis that were historically or currently used for stocking in the Lake Superior Basin. Herein we examined if morphological divergence among populations showed temporal variation, and if divergence could be attributed to selection or was indistinguishable from genetic drift. Multivariate Q_ST F_ST analysis showed evidence for divergent selection between populations. Univariate analyses suggests that the pattern observed in the multivariate analyses was largely driven by divergent selection for length and weight, and moreover by divergence between the Assinica strain and each of the Iron River and Siskiwit strains rather than divergent selection between each population pair. While it could not be determined if divergence was due to natural selection or inadvertent artificial selection in hatcheries, selected differences were consistent with patterns of domestication commonly found in salmonids.     

Testing for Spatially Divergent Selection: Comparing QST to FST

Q_ST is a standardized measure of the genetic differentiation of a quantitative trait among populations. The distribution of Q_ST''s for neutral traits can be predicted from the F_ST for neutral marker loci. To test for the neutral differentiation of a quantitative trait among populations, it is necessary to ask whether the Q_ST of that trait is in the tail of the probability distribution of neutral traits. This neutral distribution can be estimated using the Lewontin–Krakauer distribution and the F_ST from a relatively small number of marker loci. We develop a simulation method to test whether the Q_ST of a given trait is consistent with the null hypothesis of selective neutrality over space. The method is most powerful with small mean F_ST, strong selection, and a large number (>10) of measured populations. The power and type I error rate of the new method are far superior to the traditional method of comparing Q_ST and F_ST.IN 1993, Spitze (1993) and Prout and Barker (1993) introduced Q_ST, a quantitative genetic analog of Wright''s F_ST. Just as F_ST gives a standardized measure of the genetic differentiation among populations for a genetic locus, Q_ST measures the amount of genetic variance among populations relative to the total genetic variance. In the years since, Q_ST has been frequently used to test for the effects of spatially divergent (or less commonly, spatially uniform) selection (see reviews in Lynch et al. 1999; Merilä and Crnokrak 2001; McKay and Latta 2002; Howe et al. 2003; Leinonen et al. 2008; Whitlock 2008). In principle, the average Q_ST of a neutral additive quantitative trait is expected to be equal to the mean value of F_ST for neutral genetic loci. F_ST can be readily measured on commonly available genetic markers, and Q_ST can be measured as well with an appropriate breeding design in a common-garden setting. As a result, Q_ST promises to be an index of the effect of selection on the quantitative trait. If Q_ST is higher than F_ST, then this is taken as evidence of spatially divergent selection on the trait. If Q_ST is much smaller than F_ST, then this has been taken as evidence of spatially uniform stabilizing selection, which makes the trait diverge less than expected by chance.The comparison with F_ST is essential to rule out genetic drift as an alternative mechanism for phenotypic divergence among populations. Because finite populations may diverge genetically in the absence of selection, divergence must be greater than expected by drift alone if we are to conclusively demonstrate that divergent selection has played a role in genetic differentiation among populations. Therefore it has become common practice to use F_ST of putatively neutral markers as a control for the effects of genetic drift and to compare observed Q_ST values for traits to these neutral F_ST values.These comparisons follow two separate methods, to address related but distinct questions. First, many studies of quantitative genetic differentiation measure the Q_ST of many traits and the F_ST of many loci, followed by a comparison of the mean Q_ST to the mean F_ST. Such a comparison may judge whether the conditions are suitable in that species for local adaptation, that is, whether selective differences between populations are large enough relative to gene flow to allow adaptive differentiation (Whitlock 2008). We do not consider this sort of comparison in this article.The other type of comparison asks whether the Q_ST of a single trait is greater than expected by drift, as measured by F_ST. This type of comparison is most common, but it is statistically difficult. Unfortunately, as emphasized in a recent review by Whitlock (2008), there is great variation in the expected F_ST among neutral loci and among the Q_ST of different neutral traits (see Figure 1). The majority of this variation results from evolutionary differences between loci and not sampling error in the observations. Rogers and Harpending (1983) imply that the distribution of Q_ST of a single neutral trait should be approximately equivalent to that for F_ST of a single neutral locus, and this has been confirmed by simulation for traits determined by additive loci compared to biallelic marker loci (Whitlock 2008). The two distributions are similar, but there is great heterogeneity among traits or loci. As a result, to show that selection is acting on a trait, it is necessary to show that the value of Q_ST has a low probability of being observed given the distribution of neutral Q_ST.Open in a separate windowFigure 1.—The distribution of F_ST for neutral loci and Q_ST for neutral quantitative traits. The histograms show the results of simulations of a set of 10 local populations each of 100 individuals, connected by 5% migration following island model assumptions. The solid line shows the distribution predicted by the Lewontin–Krakauer distribution. The distribution of Q_ST for neutral traits is very similar to the distribution of F_ST for single neutral loci, as can be seen by their mutual good fit to the Lewontin–Krakauer distribution (Figure modified from Whitlock 2008).Comparing Q_ST to the distribution inferred from F_ST is difficult for two reasons. First, typical data sets rarely include enough loci to directly infer the distribution of F_ST without extra inferential steps. In our approach, we use the distribution of Q_ST predicted from the mean F_ST and the χ² distribution by Lewontin and Krakauer (1973) to bridge this gap. Whitlock (2008) has shown that this distribution is appropriate for nearly all realistic situations for traits determined by additive genetic effects. Second, Q_ST for a trait is rarely measured with high precision, so the position of a given estimated Q_ST value in the distribution cannot be known without error.To test the null hypothesis that the spatial distribution of a particular trait is not affected by selection, we wish to compare the observed of that trait (marked with a hat to indicate it is an estimate) to the distribution of Q_ST expected for neutral traits. Unfortunately, calculating the distribution of Q_ST for neutral traits is not straightforward, because the estimate of Q_ST for a particular trait is variable for several reasons. The estimate of Q_ST is subject to measurement error, caused by the finite samples of families and individuals in the quantitative genetic experiment. These cause error in the estimate of the additive genetic variance within populations (V_A,within) and the genetic variance among populations (V_G,among), which translate into error of the estimate of Q_ST. In addition, there is another source of variation in Q_ST among neutral traits, caused by the idiosyncrasies of the evolutionary process in each local population in the study. The true value of Q_ST for the set of populations being studied can vary tremendously around its expectation, even for neutral traits, because by chance a finite set of populations may drift in a similar direction (Whitlock 2008). As a result, measurements of Q_ST can vary because of both statistical and evolutionary variation.Fortunately, these two sources of variation are fairly well understood individually. The sampling error for the estimates of the variance components can be estimated from standard approaches, and this variation can be well approximated using information from the mean squares of the analysis of the breeding experiment (O''Hara and Merilä 2005). The variation in neutral Q_ST that results from heterogeneity of evolutionary history can be approximated by the Lewontin–Krakauer distribution (Lewontin and Krakauer 1973), if information is available on the mean Q_ST of neutral traits (Whitlock 2008). This approximation does not depend on the demographic details of the populations in question (Whitlock 2008), but the effects of deviations from assumptions of additive gene effect have not yet been tested. The mean of the distribution of values of Q_ST for neutral traits is usually not known, but fortunately the mean of the distribution of F_ST of neutral loci is expected to be approximately equal to the mean Q_ST of neutral traits (Spitze 1993), and this does not depend on demographic details (Whitlock 1999). Therefore the mean F_ST measured from a series of genetic markers thought to be selectively neutral can be combined with the Lewontin–Krakauer distribution to predict the distribution of true neutral Q_ST across the range of possible evolutionary trajectories.Given that the mean value of of neutral traits is expected to equal the mean F_ST of neutral markers under certain assumptions (discussed later), we will use as a test statistic and compare the observed quantity to the zero value proposed by the null hypothesis. We will use a traditional hypothesis testing approach, which means that we need to specify the sampling distribution of under the assumption of neutrality. Traditionally, the sampling distribution of is inferred from the data on the trait itself, for example, using bootstrapping to infer the sampling distribution. This is appropriate when calculating a confidence interval for Q_ST but is a biased measure of the sampling variance of neutral Q_ST. The variance of the sampling distribution of varies with its expected value; larger values of true Q_ST have more variable sampling distributions than traits with smaller true Q_ST. This association between Q_ST and its sampling error is quite strong, as shown in Figure 2. As a result, if the sampling properties of neutral are inferred from a trait with high Q_ST, the estimate of the variance of the null distribution will be too high, and the hypothesis test comparing to F_ST will be conservative. On the other hand, if a low Q_ST is used to estimate the variance of the null distribution, the estimated error will be too small, and the test will reject true null hypotheses too often.Open in a separate windowFigure 2.—The width of the estimated sampling distribution of varies with mean Q_ST. The solid line shows the sampling distribution of Q_ST when the true mean Q_ST value is 0.05. The dotted line shows the sampling distribution that would be estimated for Q_ST from a trait that by chance was at the first percentile of this distribution, and the dashed line shows the sampling distribution that would be inferred from a value taken at the 99th percentile. If the Q_ST of a trait differs from the expectation by chance, then the width of the sampling distribution will also be estimated with substantial error. In particular, the error variance of is overestimated with Q_ST estimates that are too high and underestimated for small Q_ST values.We address this problem by using F_ST from putatively neutral maker loci in combination with estimates of the additive genetic variance within populations to predict the sampling variance that would be expected for the Q_ST of a neutral trait. We show that the power and type I error rate of this test are greatly superior to traditional methods.     

Rapid adaptive divergence between ecotypes of an aquatic isopod inferred from FST–QST analysis

Divergent natural selection is often thought to be the principal factor driving phenotypic differentiation between populations. We studied two ecotypes of the aquatic isopod Asellus aquaticus which have diverged in parallel in several Swedish lakes. In these lakes, isopods from reed belts along the shores colonized new stonewort stands in the centre of the lakes and rapid phenotypic changes in size and pigmentation followed after colonization. We investigated if selection was likely to be responsible for these observed phenotypic changes using indirect inferences of selection (F_ST–Q_ST analysis). Average Q_ST for seven quantitative traits were higher than the average F_ST between ecotypes for putatively neutral markers (AFLPs). This suggests that divergent natural selection has played an important role during this rapid diversification. In contrast, the average Q_ST between the different reed ecotype populations was not significantly different from the mean F_ST. Genetic drift could therefore not be excluded as an explanation for the minor differences between allopatric populations inhabiting the same source habitat. We complemented this traditional F_ST–Q_ST approach by comparing the F_ST distributions across all loci (n = 67–71) with the Q_ST for each of the seven traits. This analysis revealed that pigmentation traits had diverged to a greater extent and at higher evolutionary rates than size‐related morphological traits. In conclusion, this extended and detailed type of F_ST–Q_ST analysis provides a powerful method to infer adaptive phenotypic divergence between populations. However, indirect inferences about the operation of divergent selection should be analyzed on a per‐trait basis and complemented with detailed ecological information.     

Molecular and Quantitative Genetic Differentiation in Sitobion avenae Populations from Both Sides of the Qinling Mountains

Quantitative trait differences are often assumed to be correlated with molecular variation, but the relationship is not certain, and empirical evidence is still scarce. To address this issue, we sampled six populations of the cereal aphid Sitobion avenae from areas north and south of the Qinling Mountains, and characterized their molecular variation at seven microsatellite loci and quantitative variation at nine life-history traits. Our results demonstrated that southern populations had slightly longer developmental times of nymphs but much higher lifetime fecundity, compared to northern populations. Of the nine tested quantitative characters, eight differed significantly among populations within regions, as well as between northern and southern regions. Genetic differentiation in neutral markers was likely to have been caused by founder events and drift. Increased subdivision for quantitative characters was found in northern populations, but reduced in southern populations. This phenomenon was not found for molecular characters, suggesting the decoupling between molecular and quantitative variation. The pattern of relationships between F_ST and Q_ST indicated divergent selection and suggested that local adaptation play a role in the differentiation of life-history traits in tested S. avenae populations, particularly in those traits closely related to reproduction. The main role of natural selection over genetic drift was also supported by strong structural differences in G-matrices among S. avenae populations. However, cluster analyses did not result in two groups corresponding to northern and southern regions. Genetic differentiation between northern and southern populations in neutral markers was low, indicating considerable gene flow between them. The relationship between molecular and quantitative variation, as well as its implications for differentiation and evolution of S. avenae populations, was discussed.     

Phenotypic divergence of the common toad (Bufo bufo) along an altitudinal gradient: evidence for local adaptation

Variation in the environment can induce different patterns of genetic and phenotypic differentiation among populations. Both neutral processes and selection can influence phenotypic differentiation. Altitudinal phenotypic variation is of particular interest in disentangling the interplay between neutral processes and selection in the dynamics of local adaptation processes but remains little explored. We conducted a common garden experiment to study the phenotypic divergence in larval life-history traits among nine populations of the common toad (Bufo bufo) along an altitudinal gradient in France. We further used correlation among population pairwise estimates of quantitative trait (Q_ST) and neutral genetic divergence (F_ST from neutral microsatellite markers), as well as altitudinal difference, to estimate the relative role of divergent selection and neutral genetic processes in phenotypic divergence. We provided evidence for a neutral genetic differentiation resulting from both isolation by distance and difference in altitude. We found evidence for phenotypic divergence along the altitudinal gradient (faster development, lower growth rate and smaller metamorphic size). The correlation between pairwise Q_STs–F_STs and altitude differences suggested that this phenotypic differentiation was most likely driven by altitude-mediated selection rather than by neutral genetic processes. Moreover, we found different divergence patterns for larval traits, suggesting that different selective agents may act on these traits and/or selection on one trait may constrain the evolution on another through genetic correlation. Our study highlighted the need to design more integrative studies on the common toad to unravel the underlying processes of phenotypic divergence and its selective agents in the context of environmental clines.     

Generalization of the QST framework in hierarchically structured populations: Impacts of inbreeding and dominance

Q_ST is a differentiation parameter based on the decomposition of the genetic variance of a trait. In the case of additive inheritance and absence of selection, it is analogous to the genic differentiation measured on individual loci, F_ST. Thus, Q_ST?F_ST comparison is used to infer selection: selective divergence when Q_ST > F_ST, or convergence when Q_ST < F_ST. The definition of Q‐statistics was extended to two‐level hierarchical population structures with Hardy–Weinberg equilibrium. Here, we generalize the Q‐statistics framework to any hierarchical population structure. First, we developed the analytical definition of hierarchical Q‐statistics for populations not at Hardy–Weinberg equilibrium. We show that the Q‐statistics values obtained with the Hardy–Weinberg definition are lower than their corresponding F‐statistics when F_IS > 0 (higher when F_IS < 0). Then, we used an island model simulation approach to investigate the impact of inbreeding and dominance on the Q_ST?F_ST framework in a hierarchical population structure. We show that, while differentiation at the lower hierarchical level (Q_SR) is a monotonic function of migration, differentiation at the upper level (Q_RT) is not. In the case of additive inheritance, we show that inbreeding inflates the variance of Q_RT, which can increase the frequency of Q_RT > F_RT cases. We also show that dominance drastically reduces Q‐statistics below F‐statistics for any level of the hierarchy. Therefore, high values of Q‐statistics are good indicators of selection, but low values are not in the case of dominance.     

Natural Selection and Neutral Evolution Jointly Drive Population Divergence between Alpine and Lowland Ecotypes of the Allopolyploid Plant Anemone multifida (Ranunculaceae)

Population differentiation can be driven in large part by natural selection, but selectively neutral evolution can play a prominent role in shaping patters of population divergence. The decomposition of the evolutionary history of populations into the relative effects of natural selection and selectively neutral evolution enables an understanding of the causes of population divergence and adaptation. In this study, we examined heterogeneous genomic divergence between alpine and lowland ecotypes of the allopolyploid plant, Anemone multifida. Using peak height and dominant AFLP data, we quantified population differentiation at non-outlier (neutral) and outlier loci to determine the potential contribution of natural selection and selectively neutral evolution to population divergence. We found 13 candidate loci, corresponding to 2.7% of loci, with signatures of divergent natural selection between alpine and lowland populations and between alpine populations (Fst  = 0.074–0.445 at outlier loci), but neutral population differentiation was also evident between alpine populations (F_ST  = 0.041–0.095 at neutral loci). By examining population structure at both neutral and outlier loci, we determined that the combined effects of selection and neutral evolution are associated with the divergence of alpine populations, which may be linked to extreme abiotic conditions and isolation between alpine sites. The presence of outlier levels of genetic variation in structured populations underscores the importance of separately analyzing neutral and outlier loci to infer the relative role of divergent natural selection and neutral evolution in population divergence.     

Divergent selection along climatic gradients in a rare central European endemic species,Saxifraga sponhemica

Background and Aims The effects of habitat fragmentation on quantitative genetic variation in plant populations are still poorly known. Saxifraga sponhemica is a rare endemic of Central Europe with a disjunct distribution, and a stable and specialized habitat of treeless screes and cliffs. This study therefore used S. sponhemica as a model species to compare quantitative and molecular variation in order to explore (1) the relative importance of drift and selection in shaping the distribution of quantitative genetic variation along climatic gradients; (2) the relationship between plant fitness, quantitative genetic variation, molecular genetic variation and population size; and (3) the relationship between the differentiation of a trait among populations and its evolvability.Methods Genetic variation within and among 22 populations from the whole distribution area of S. sponhemica was studied using RAPD (random amplified polymorphic DNA) markers, and climatic variables were obtained for each site. Seeds were collected from each population and germinated, and seedlings were transplanted into a common garden for determination of variation in plant traits.Key Results In contrast to previous results from rare plant species, strong evidence was found for divergent selection. Most population trait means of S. sponhemica were significantly related to climate gradients, indicating adaptation. Quantitative genetic differentiation increased with geographical distance, even when neutral molecular divergence was controlled for, and Q_ST exceeded F_ST for some traits. The evolvability of traits was negatively correlated with the degree of differentiation among populations (Q_ST), i.e. traits under strong selection showed little genetic variation within populations. The evolutionary potential of a population was not related to its size, the performance of the population or its neutral genetic diversity. However, performance in the common garden was lower for plants from populations with reduced molecular genetic variation, suggesting inbreeding depression due to genetic erosion.Conclusions The findings suggest that studies of molecular and quantitative genetic variation may provide complementary insights important for the conservation of rare species. The strong differentiation of quantitative traits among populations shows that selection can be an important force for structuring variation in evolutionarily important traits even for rare endemic species restricted to very specific habitats.     

Contrasting Levels of Variation in Neutral and Quantitative Genetic Loci on Island Populations of Moor Frogs (<Emphasis Type="Italic">Rana arvalis</Emphasis>)

Reduced levels of genetic variability and a prominent differentiation in both neutral marker genes and phenotypic traits are typical for many island populations as compared to their mainland conspecifics. However, whether genetic diversity in neutral marker genes reflects genetic variability in quantitative traits, and thus, their evolutionary potential, remains typically unclear. Moreover, the phenotypic differentiation on islands could be attributable to phenotypic plasticity, selection or drift; something which seldom has been tested. Using eight polymorphic microsatellite loci and quantitative genetic breeding experiments we conducted a detailed comparison on genetic variability and differentiation between Nordic islands (viz. Gotland, Öland and Læsø) and neighbouring mainland populations of moor frogs (Rana arvalis). As expected, the neutral variation was generally lower in island than in mainland populations. But as opposed to this, higher levels of additive genetic variation (V _A) in body size and tibia length were found on the island of Gotland as compared to the mainland population. When comparing the differentiation seen in neutral marker genes (F _ST) with the differentiation in genes coding quantitative traits (Q _ST) two different evolutionary scenarios were found: while selection might explain a smaller size of moor frogs on Gotland, the differentiation seen in tibia length could be explained by genetic drift. These results highlight the limited utility of microsatellite loci alone in inferring the causes behind an observed phenotypic differentiation, or in predicting the amount of genetic variation in ecologically important quantitative traits.     

QST < FST As a signature of canalization

A key aim of evolutionary biology – inferring the action of natural selection on wild species – can be achieved by comparing neutral genetic differentiation between populations (F_ST) with quantitative genetic variation (Q_ST). Each of the three possible outcomes of comparisons of Q_ST and F_ST (Q_ST > F_ST, Q_ST = F_ST, Q_ST < F_ST) is associated with an inference (diversifying selection, genetic drift, uniform selection, respectively). However, published empirical and theoretical studies have focused on the Q_ST > F_ST outcome. We believe that this reflects the absence of a straightforward biological interpretation of the Q_ST < F_ST pattern. We here report recent evidence of this neglected evolutionary pattern, provide guidelines to its interpretation as either a canalization phenomenon or a consequence of uniform selection and discuss the significant importance this issue will have for the area of evolutionary biology.     

Genetic drift and uniform selection shape evolution of most traits in <Emphasis Type="Italic">Eugenia dysenterica</Emphasis> DC. (Myrtaceae)

Knowing how microevolutionary processes, such as genetic drift and natural selection, shape variation in adaptive traits is strategic for conservation measures. One way to estimate local adaptation is to compare divergences in quantitative traits (Q_ST) and neutral loci (F_ST). Therefore, we have assessed the pattern of phenotypic and molecular genetic divergence among natural subpopulations of the fruit tree Eugenia dysenterica DC. A provenance and progeny test was performed to assess the quantitative traits of the subpopulations collected in a wide distribution area of the species in the Brazilian Cerrado. The sampled environments are in a biodiversity hotspot with heterogeneous soil and climate conditions. By associating quantitative trait variation in initial seedling development with neutral microsatellite marker variation, we tested the local adaptation of the traits by the Q_ST–F_ST contrast. Genetic drift was prevalent in the phenotypic differentiation among the subpopulations, although the traits seedling emergence time and root green mass, which are relevant for adaptation to the Cerrado climate, showed signs of uniform selection. Our results suggest that E. dysenterica has a spatial genetic structure divided into two large groups, separated by a line that divides the Cerrado biome in a southwestern to northeastern direction. This structure must be taken into account for managing E. dysenterica genetic resources both for conservation and breeding purposes.     

Heritability and quantitative genetic divergence of serotiny,a fire-persistence plant trait

Background and Aims

Although it is well known that fire acts as a selective pressure shaping plant phenotypes, there are no quantitative estimates of the heritability of any trait related to plant persistence under recurrent fires, such as serotiny. In this study, the heritability of serotiny in Pinus halepensis is calculated, and an evaluation is made as to whether fire has left a selection signature on the level of serotiny among populations by comparing the genetic divergence of serotiny with the expected divergence of neutral molecular markers (Q_ST–F_ST comparison).

Methods

A common garden of P. halepensis was used, located in inland Spain and composed of 145 open-pollinated families from 29 provenances covering the entire natural range of P. halepensis in the Iberian Peninsula and Balearic Islands. Narrow-sense heritability (h²) and quantitative genetic differentiation among populations for serotiny (Q_ST) were estimated by means of an ‘animal model’ fitted by Bayesian inference. In order to determine whether genetic differentiation for serotiny is the result of differential natural selection, Q_ST estimates for serotiny were compared with F_ST estimates obtained from allozyme data. Finally, a test was made of whether levels of serotiny in the different provenances were related to different fire regimes, using summer rainfall as a proxy for fire regime in each provenance.

Key Results

Serotiny showed a significant narrow-sense heritability (h²) of 0·20 (credible interval 0·09–0·40). Quantitative genetic differentiation among provenances for serotiny (Q_ST = 0·44) was significantly higher than expected under a neutral process (F_ST = 0·12), suggesting adaptive differentiation. A significant negative relationship was found between the serotiny level of trees in the common garden and summer rainfall of their provenance sites.

Conclusions

Serotiny is a heritable trait in P. halepensis, and selection acts on it, giving rise to contrasting serotiny levels among populations depending on the fire regime, and supporting the role of fire in generating genetic divergence for adaptive traits.     

Evidence for adaptive phenotypic differentiation in Baltic Sea sticklebacks

The evidence for adaptive phenotypic differentiation in mobile marine species remains scarce, partly due to the difficulty of obtaining quantitative genetic data to demonstrate the genetic basis of the observed phenotypic differentiation. Using a combination of phenotypic and molecular genetic approaches, we elucidated the relative roles of natural selection and genetic drift in explaining lateral plate number differentiation in threespine sticklebacks (Gasterosteus aculeatus) across the entire Baltic Sea basin (approximately 392 000 km²). We found that phenotypic differentiation (P_ST = 0.213) in plate number exceeded that in neutral markers (F_ST = 0.008), suggesting an adaptive basis for the observed differentiation. Because a close correspondence was found between plate phenotype and genotype at a quantitative trait loci (QTL; STN381) tightly linked to the gene (Ectodysplasin) underlying plate variation, the evidence for adaptive differentiation was confirmed by comparison of F_ST at the QTL (F_STQ = 0.089) with F_ST at neutral marker loci. Hence, the results provide a comprehensive demonstration of adaptive phenotypic differentiation in a high‐gene‐flow marine environment with direct, rather than inferred, verification for the genetic basis of this differentiation. In general, the results illustrate the utility of P_ST–F_ST–F_STQ comparisons in uncovering footprints of natural selection and evolution and add to the growing evidence for adaptive genetic differentiation in high‐gene‐flow marine environments, including that of the relatively young Baltic Sea.