Molecular and quantitative genetic divergence among populations of house mice with known evolutionary histories

Evolutionary biologists have long been interested in the processes influencing population differentiation, but separating the effects of neutral and adaptive evolution has been an obstacle for studies of population subdivision. A recently developed method allows tests of whether disruptive (ie, spatially variable) or stabilizing (ie, spatially uniform) selection is influencing phenotypic differentiation among subpopulations. This method, referred to as the F(ST) vs Q(ST) comparison, separates the total additive genetic variance into within- and among-population components and evaluates this level of differentiation against a neutral hypothesis. Thus, levels of neutral molecular (F(ST)) and quantitative genetic (Q(ST)) divergence are compared to evaluate the effects of selection and genetic drift on phenotypic differentiation. Although the utility of such comparisons appears great, its accuracy has not yet been evaluated in populations with known evolutionary histories. In this study, F(ST) vs Q(ST) comparisons were evaluated using laboratory populations of house mice with known evolutionary histories. In this model system, the F(ST) vs Q(ST) comparisons between the selection groups should reveal quantitative trait differentiation consistent with disruptive selection, while the F(ST) vs Q(ST) comparisons among lines within the selection groups should suggest quantitative trait differentiation in agreement with drift. We find that F(ST) vs Q(ST) comparisons generally produce the correct evolutionary inference at each level in the population hierarchy. Additionally, we demonstrate that when strong selection is applied between populations Q(ST) increases relative to Q(ST) among populations diverging by drift. Finally, we show that the statistical properties of Q(ST), a variance component ratio, need further investigation.     

Experimental demonstration of a causal relationship between heterogeneity of selection and genetic differentiation in quantitative traits

Comparisons of estimates of genetic differentiation at molecular markers (F(ST)) and at quantitative traits (Q(ST)) are a means of inferring the level and heterogeneity of selection in natural populations. However, such comparisons are questionable because they require that the influence of drift and selection on Q(ST) be detectable over possible background influences of environmental or nonadditive genetic effects on Q(ST)-values. Here we test this using an experimental evolution approach in metapopulations of Arabidopsis thaliana experiencing different levels of drift and selection heterogeneity. We estimated the intensity and heterogeneity of selection on morphological and phenological traits via selection differentials. We demonstrate that Q(ST)-values increased with increasing selection heterogeneity when genetic drift was limited. The effect of selection on Q(ST) was thus detectable despite significant genotype-by-environment interactions that most probably biased the estimates of genetic differentiation. Although they cannot be used as a direct validation of the conclusions of prior studies, our results strongly support both the relevance of Q(ST) as an estimator of genetic differentiation and the role of local selection in shaping the genetic differentiation of natural populations.     

Adaptive divergence for a fitness-related trait among invasive Ambrosia artemisiifolia populations in France

The impact of natural selection on the adaptive divergence of invasive populations can be assessed by testing the null hypothesis that the extent of quantitative genetic differentiation (Q(ST) ) would be similar to that of neutral molecular differentiation (F(ST) ). Using eight microsatellite loci and a common garden approach, we compared Q(ST) and F(ST) among ten populations of an invasive species Ambrosia artemisiifolia (common ragweed) in France. In a common garden study with varying water and nutrient levels, we measured Q(ST) for five traits (height, total biomass, reproductive allocation, above- to belowground biomass ratio, and days to flowering). Although low F(ST) indicated weak genetic structure and strong gene flow among populations, we found significant diversifying selection (Q(ST) > F(ST) ) for reproductive allocation that may be closely related to fitness. It suggests that abiotic conditions may have exerted selection pressure on A. artemisiifolia populations to differentiate adaptively, such that populations at higher altitude or latitude evolved greater reproductive allocation. As previous studies indicate multiple introductions from various source populations of A. artemisiifolia in North America, our results suggest that the admixture of introduced populations may have increased genetic diversity and additive genetic variance, and in turn, promoted the rapid evolution and adaptation of this invasive 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.     

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.     

Comparing Bayesian estimates of genetic differentiation of molecular markers and quantitative traits: an application to Pinus sylvestris

Comparison of the level of differentiation at neutral molecular markers (estimated as F(ST) or G(ST)) with the level of differentiation at quantitative traits (estimated as Q(ST)) has become a standard tool for inferring that there is differential selection between populations. We estimated Q(ST) of timing of bud set from a latitudinal cline of Pinus sylvestris with a Bayesian hierarchical variance component method utilizing the information on the pre-estimated population structure from neutral molecular markers. Unfortunately, the between-family variances differed substantially between populations that resulted in a bimodal posterior of Q(ST) that could not be compared in any sensible way with the unimodal posterior of the microsatellite F(ST). In order to avoid publishing studies with flawed Q(ST) estimates, we recommend that future studies should present heritability estimates for each trait and population. Moreover, to detect variance heterogeneity in frequentist methods (ANOVA and REML), it is of essential importance to check also that the residuals are normally distributed and do not follow any systematically deviating trends.     

The effects of dominance, regular inbreeding and sampling design on Q(ST), an estimator of population differentiation for quantitative traits

To test whether quantitative traits are under directional or homogenizing selection, it is common practice to compare population differentiation estimates at molecular markers (F(ST)) and quantitative traits (Q(ST)). If the trait is neutral and its determinism is additive, then theory predicts that Q(ST) = F(ST), while Q(ST) > F(ST) is predicted under directional selection for different local optima, and Q(ST) < F(ST) is predicted under homogenizing selection. However, nonadditive effects can alter these predictions. Here, we investigate the influence of dominance on the relation between Q(ST) and F(ST) for neutral traits. Using analytical results and computer simulations, we show that dominance generally deflates Q(ST) relative to F(ST). Under inbreeding, the effect of dominance vanishes, and we show that for selfing species, a better estimate of Q(ST) is obtained from selfed families than from half-sib families. We also compare several sampling designs and find that it is always best to sample many populations (>20) with few families (five) rather than few populations with many families. Provided that estimates of Q(ST) are derived from individuals originating from many populations, we conclude that the pattern Q(ST) > F(ST), and hence the inference of directional selection for different local optima, is robust to the effect of nonadditive gene actions.     

Genetic differences among populations in sexual dimorphism: evidence for selection on males in a dioecious plant

Genetic variation among populations in the degree of sexual dimorphism may be a consequence of selection on one or both sexes. We analysed genetic parameters from crosses involving three populations of the dioecious plant Silene latifolia, which exhibits sexual dimorphism in flower size, to determine whether population differentiation was a result of selection on one or both sexes. We took the novel approach of comparing the ratio of population differentiation of a quantitative trait (Q(ST) ) to that of neutral genetic markers (F(ST) ) for males vs. females. We attributed 72.6% of calyx width variation in males to differences among populations vs. only 6.9% in females. The Q(ST) /F(ST) ratio was 4.2 for males vs. 0.4 for females, suggesting that selection on males is responsible for differentiation among populations in calyx width and its degree of sexual dimorphism. This selection may be indirect via genetic correlations with other morphological and physiological traits.     

Strong divergence in trait means but not in plasticity across hatchery and wild populations of sea-run brown trout Salmo trutta

There is ample evidence that organisms adapt to their native environment when gene flow is restricted. However, evolution of plastic responses across discrete environments is less well examined. We studied divergence in means and plasticity across wild and hatchery populations of sea-run brown trout (Salmo trutta) in a common garden experiment with two rearing environments (hatchery and a nearly natural experimental stream). Since natural and hatchery environments differ, this arrangement provides an experiment in contemporary adaptation across the two environments. A Q(ST) - F(ST) approach was used to investigate local adaptation in survival and growth over the first summer. We found evidence for divergent selection in survival in 1 year and in body length in both years and rearing environments. In general, the hatchery populations had higher survival and larger body size in both environments. Q(ST) in body size did not differ between the rearing environments, and constitutive divergence in the means was in all cases stronger than divergence in the plastic responses. These results suggest that in this system, constitutive changes in mean trait values are more important for local adaptation than increased plasticity. In addition, ex situ rearing conditions induce changes in trait means that are adaptive in the hatchery, but potentially harmful in the wild, suggesting that hatchery rearing is likely to be a suboptimal management strategy for trout populations facing selection in the stream environment.     

Comparison of quantitative and molecular genetic variation of native vs. invasive populations of purple loosestrife (Lythrum salicaria L., Lythraceae)

Study of adaptive evolutionary changes in populations of invasive species can be advanced through the joint application of quantitative and population genetic methods. Using purple loosestrife as a model system, we investigated the relative roles of natural selection, genetic drift and gene flow in the invasive process by contrasting phenotypical and neutral genetic differentiation among native European and invasive North American populations ( Q _ST −  F _ST analysis). Our results indicate that invasive and native populations harbour comparable levels of amplified fragment length polymorphism variation, a pattern consistent with multiple independent introductions from a diverse European gene pool. However, it was observed that the genetic variation reduced during subsequent invasion, perhaps by founder effects and genetic drift. Comparison of genetically based quantitative trait differentiation ( Q _ST) with its expectation under neutrality ( F _ST) revealed no evidence of disruptive selection ( Q _ST >  F _ST) or stabilizing selection ( Q _ST <  F _ST). One exception was found for only one trait (the number of stems) showing significant sign of stabilizing selection across all populations. This suggests that there are difficulties in distinguishing the effects of nonadaptive population processes and natural selection. Multiple introductions of purple loosestrife may have created a genetic mixture from diverse source populations and increased population genetic diversity, but its link to the adaptive differentiation of invasive North American populations needs further research.     

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.     

The effect of neutral nonadditive gene action on the quantitative index of population divergence

For neutral additive genes, the quantitative index of population divergence (Q(ST)) is equivalent to Wright's fixation index (F(ST)). Thus, divergent or convergent selection is usually invoked, respectively, as a cause of the observed increase (Q(ST) > F(ST)) or decrease (Q(ST) < F(ST)) of Q(ST) from its neutral expectation (Q(ST) = F(ST)). However, neutral nonadditive gene action can mimic the additive expectations under selection. We have studied theoretically the effect of consecutive population bottlenecks on the difference F(ST) - Q(ST) for two neutral biallelic epistatic loci, covering all types of marginal gene action. With simple dominance, Q(ST) < F(ST) for only low to moderate frequencies of the recessive alleles; otherwise, Q(ST) > F(ST). Additional epistasis extends the condition Q(ST) < F(ST) to a broader range of frequencies. Irrespective of the type of nonadditive action, Q(ST) < F(ST) generally implies an increase of both the within-line additive variance after bottlenecks over its ancestral value (V(A)) and the between-line variance over its additive expectation (2F(ST)V(A)). Thus, both the redistribution of the genetic variance after bottlenecks and the F(ST) - Q(ST) value are governed largely by the marginal properties of single loci. The results indicate that the use of the F(ST) - Q(ST) criterion to investigate the relative importance of drift and selection in population differentiation should be restricted to pure additive traits.     

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.     

Patterns of genetic variability within and among populations of wild radish,Raphanus raphanistrum (Brassicaceae)

The genetic structure of populations is an important determinant of the evolutionary potential of a species. Colonizing plants tend to be characterized by low within- and high among-population variability. Genetic differentiation of both floral traits and isozymes was studied in six populations of wild radish (Raphanus raphanistrum). Evidence for differentiation in both sets of traits was found, but patterns of differentiation of floral traits did not coincide with isozyme differentiation. Contrary to most colonizing species, wild radish showed high within- and only moderate among-population variability at isozyme loci. In addition, levels of differentiation did not correspond to geographic distance between the populations. These results are likely due at least in part to the self-incompatibility system of this species, long-distance movement of large numbers of wild radish seeds by humans, and introgression from cultivated radish (R. sativus).     

Divergent selection on flowering phenology but not on floral morphology between two closely related orchids

Closely related species often differ in traits that influence reproductive success, suggesting that divergent selection on such traits contribute to the maintenance of species boundaries. Gymnadenia conopsea ss. and Gymnadenia densiflora are two closely related, perennial orchid species that differ in (a) floral traits important for pollination, including flowering phenology, floral display, and spur length, and (b) dominant pollinators. If plant–pollinator interactions contribute to the maintenance of trait differences between these two taxa, we expect current divergent selection on flowering phenology and floral morphology between the two species. We quantified phenotypic selection via female fitness in one year on flowering start, three floral display traits (plant height, number of flowers, and corolla size) and spur length, in six populations of G. conopsea s.s. and in four populations of G. densiflora. There was indication of divergent selection on flowering start in the expected direction, with selection for earlier flowering in two populations of the early‐flowering G. conopsea s.s. and for later flowering in one population of the late‐flowering G. densiflora. No divergent selection on floral morphology was detected, and there was no significant stabilizing selection on any trait in the two species. The results suggest ongoing adaptive differentiation of flowering phenology, strengthening this premating reproductive barrier between the two species. Synthesis: This study is among the first to test whether divergent selection on floral traits contribute to the maintenance of species differences between closely related plants. Phenological isolation confers a substantial potential for reproductive isolation, and divergent selection on flowering time can thus greatly influence reproductive isolation and adaptive differentiation.     

High levels of genetic differentiation of Oryza officinalis Wall. ex Watt. from China

In order to determine the population genetic structure of wild rice (Oryza officinalis Wall. ex Watt.), an endangered tropical and subtropical species, allozyme diversity encoded by 24 loci was analyzed electrophoretically in 145 individuals of eight natural populations from Hainan, Guangxi, and Yunnan provinces, China. A fairly high genetic differentiation (F(ST) = 0.882 and mean I = 0.786) was found among the studied populations. Our results suggest that restricted gene flow may play a significant role in shaping such a population genetic structure. In addition, high genetic differentiation among populations within a geographically limited region may stem from a reduced population size and consequent genetic drift.     

Population differentiation in G matrix structure due to natural selection in Rana temporaria

The additive genetic variance-covariance matrix (G) is a concept central to discussions about evolutionary change over time in a suite of traits. However, at the moment we do not know how fast G itself changes as a consequence of selection or how sensitive it is to environmental influences. We investigated possible evolutionary divergence and environmental influences on G using data from a factorial common-garden experiment where common frog (Rana temporaria) tadpoles from two divergent populations were exposed to three different environmental treatments. G-matrices were estimated using an animal model approach applied to data from a NCII breeding design. Matrix comparisons using both Flury and multivariate analysis of variance methods revealed significant differences in G matrices both between populations and between treatments within populations, the former being generally larger than the latter. Comparison of levels of population differentiation in trait means using Q(ST) indices with that observed in microsatellite markers (F(ST)) revealed that the former values generally exceeded the neutral expectation set by F(ST). Hence, the results suggest that intraspecific divergence in G matrix structure has occurred mainly due to natural selection.     

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.     

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.     

闽江水系绒螯蟹的表型性状差异与分子遗传差异的比较与分析

表型性状差异(differentiation in phenotypictraits,PST)和分子遗传差异(differentiation at neutral molecular markers,FST)是近期进化生物学的研究热点之一。闽江水系是我国中华绒螯蟹与合浦绒螯蟹的主要混杂地域,是研究绒螯蟹遗传与进化的理想地之一。为探讨闽江水系绒螯蟹的PST和FST,以2009和2010年度闽江水系的133个绒螯蟹样本为材料,进行了14个表型数量性状差异分析和6个微卫星标记的遗传差异分析。结果发现:除3个表型性状不存在显著差异外,其他表型性状在不同年份间均存在极显著差异(P<0.01);2009年绒螯蟹的平均期望杂合度极显著高于2010年绒螯蟹(P=0.008),而平均等位基因丰富度、观测杂合度和近交系数均不存在显著差异(P=0.136～0.675);年份间的平均FST为0.1429;通过对PST与FST的比较发现,除第二步足掌节长度(F2)性状外,其他表型数量性状的PST值均高于FST值,表明这些性状均受到了较明显的选择压力。该文研究结果为绒螯蟹的分子进化研究积累了资料,也为其他水产生物的PST和FST比较研究提供了参考。