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.     

Decoupling of differentiation between traits and their underlying genes in response to divergent selection

We dissected the relationship between genetic differentiation (Q(ST)) for a trait and its underlying genes (G(STq), differentiation for a quantitative locus) in an evolutionary context, with the aim of identifying the conditions in which these two measurements are decoupled. We used two parameters (θ(B) and θ(W)) scaling the contributions of inter- and intrapopulation allelic covariation between genes controlling the trait of interest. We monitored the changes in θ(B) and θ(W), Q(ST) and G(STq) over successive generations of divergent and stabilizing selection, in simulations for an outcrossing species with extensive gene flow. The dynamics of these parameters are characterized by two phases. Initially, during the earliest generations, differentiation of the trait increases very rapidly and the principal and immediate driver of Q(ST) is θ(B). During subsequent generations, G(STq) increases steadily and makes an equal contribution to Q(ST). These results show that selection first captures beneficial allelic associations at different loci at different populations, and then targets changes in allelic frequencies. The same patterns are observed when environmental change modifies divergent selection, as shown by the very rapid response of θ(B) to the changes of selection regimes. We compare our results with previous experimental findings and consider their relevance to the detection of molecular signatures of natural selection.     

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.     

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.     

Adaptive gradients and isolation-by-distance with postglacial migration in Picea sitchensis

Fossil pollen records suggest rapid migration of tree species in response to Quaternary climate warming. Long-distance dispersal and high gene flow would facilitate rapid migration, but would initially homogenize variation among populations. However, contemporary clinal variation in adaptive traits along environmental gradients shown in many tree species suggests that local adaptation can occur during rapid migration over just a few generations in interglacial periods. In this study, we compared growth performance and pollen genetic structure among populations to investigate how populations of Sitka spruce (Picea sitchensis) have responded to local selection along the historical migration route. The results suggest strong adaptive divergence among populations (average Q(ST)=0.61), corresponding to climatic gradients. The population genetic structure, determined by microsatellite markers (R(ST)=0.09; F(ST)=0.11), was higher than previous estimates from less polymorphic genetic markers. The significant correlation between geographic and pollen haplotype genetic (R(ST)) distances (r=0.73, P<0.01) indicates that the current genetic structure has been shaped by isolation-by-distance, and has developed in relatively few generations. This suggests relatively limited gene flow among populations on a recent timescale. Gene flow from neighboring populations may have provided genetic diversity to founder populations during rapid migration in the early stages of range expansion. Increased genetic diversity subsequently enhanced the efficiency of local selection, limiting gene flow primarily to among similar environments and facilitating the evolution of adaptive clinal variation along environmental gradients.     

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.     

Strong linkage disequilibrium near the selected Yr17 resistance gene in a wheat experimental population

Dynamic management (DM) is a method of genetic resources conservation that aims at maintaining evolutionary process in subdivided populations cultivated in contrasted environments. Such populations are often submitted to strong natural selection as it was the case for experimental wheat populations maintained under DM. Understanding impacts of selection on genetic diversity around selected genes is necessary for the middle-term maintenance of genetic variability in DM populations. Evolution of diversity at six neutral markers located near the yellow rust resistance gene Yr17 has been studied for the parental lines and for generations 1, 5, 10 and 17 in one of the DM populations. Yr17 provided complete resistance to yellow rust in France until 1997 and thus was suspected to be under strong selection. The gene is located on a fragment introgressed in winter wheat from a wild species. The presence of the gene was estimated using a marker closely related to the gene. We showed that the Yr17 gene has been selected between generations 5 and 10. Generally, selection tends to reduce diversity around selected genes, generating linkage disequilibrium (LD) between a gene and adjacent markers. Here, the major effect of the Yr17 gene selection was a reduction of multilocus diversity and the maintenance of strong pre-existing LD in the zone surrounding the gene for a distance of 20 cM. As expected, the presence of the exogenous introgression was responsible for restrictions to recombination which contributed to the maintenance of strong correlations between loci. However, we found a noticeable number of recombinations around the gene indicating a progressive incorporation of the fragment into the wheat genome.     

Joint effect of selection, linkage and partial inbreeding on additive genetic variance in an infinite population

Under the inifinitesimal model of gene effects, selection reduces the additive genetic variance by inducing negative linkage disequilibrium among selected genes. If the selected genes are linked, the decay of linkage disequilibrium is delayed, and the reduction of additive genetic variance is enhanced. Inbreeding in an infinite population also alters the additive genetic variance through the generation of positive association among genes within a locus. In the present study, the joint effect of selection, linkage and partial inbreeding (partial selfing or partial full-sib mating) on the additive genetic variance was modeled. The recurrence relations of the additive genetic variance between successive generations and the prediction equation of the asymptotic additive genetic variance were derived. Numerical computation showed that although partially inbred populations initially maintain larger genetic variances, the accumulated effect of selection overrides the effect of inbreeding. Stochastic simulation was carried out to check the precision of prediction, showing that the obtained equations give a satisfactory prediction during initial generations. However, the predicted values always overestimate the simulated values, especially in later generations. Based on these results, possible extensions and perspectives of the assumed model were discussed.     

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.     

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.     

Drift and selection influence geographic variation at immune loci of prairie-chickens

Previous studies of immunity in wild populations have focused primarily on genes of the major histocompatibility complex (MHC); however, studies of model species have identified additional immune-related genes that also affect fitness. In this study, we sequenced five non-MHC immune genes in six greater prairie-chicken (Tympanuchus cupido) populations that have experienced varying degrees of genetic drift as a consequence of population bottlenecks and fragmentation. We compared patterns of geographic variation at the immune genes with six neutral microsatellite markers to investigate the relative effects of selection and genetic drift. Global F(ST) outlier tests identified positive selection on just one of five immune genes (IAP-1) in one population. In contrast, at other immune genes, standardized G'(ST) values were lower than those at microsatellites for a majority of pairwise population comparisons, consistent with balancing selection or with species-wide positive or purifying selection resulting in similar haplotype frequencies across populations. The effects of genetic drift were also evident as summary statistics (e.g., Tajima's D) did not differ from neutrality for the majority of cases, and immune gene diversity (number of haplotypes per gene) was correlated positively with population size. In summary, we found that both genetic drift and selection shaped variation at the five immune genes, and the strength and type of selection varied among genes. Our results caution that neutral forces, such as drift, can make it difficult to detect current selection on genes.     

长期人工选择对群体遗传结构的影响——一项计算机模拟研究 Ⅱ.长期选择对群体遗传方差的作用

本文给出了显性与超显性模型下加性方差的分剖公式,为研究选择作用下基因间关系的变化提供了有力的方法。并模拟研究了群体大小、连锁强度与遗传力水平对遗传方差变化的影响。小群体中遗传方差在世代间波动很大;大群体中则稳定下降、波动较小。选择作用下平衡加性方差下降很快,特别是高遗传力性状。紧密连锁在小群体中一方面降低选择反应,一方面维持了更多的加性方差,从而使得预测长期选择反应甚为困难。     

Contrasting evolutionary forces driving population structure at expressed sequence tag polymorphisms, allozymes and quantitative traits in white spruce

Patterns of variation in quantitative characters and genetic markers were compared among six regional populations of white spruce [Picea glauca (Moench) Voss]. Although some phenotypic characters were correlated with latitude (r = 0.791), longitude (r = -0.796) and precipitation during the growing season (r = 0.789), variability at genetic markers was not correlated with geographical or bioclimatic variables, and followed neutral expectations. Estimates of genetic diversity and population differentiation for 14 allozymes (translated regions of coding genes) were essentially indistinguishable from those observed for 11 expressed sequence tag polymorphisms (ESTPs) from untranslated regions of coding genes. Variation among populations for quantitative traits such as eighth year height (Q(ST) = 0.082), thirteenth year height (Q(ST) = 0.069), total wood density (Q(ST) = 0.102) and date of budset (Q(ST) = 0.246), was greater than for allozymes (G(ST) = 0.014) and ESTPs (G(ST) = 0.019). These trends suggest a strong adaptive response in quantitative traits, contrasting to allozymes and ESTPs where no selective response could be detected and where populations appeared to be essentially in a migration-drift equilibrium.     

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.     

Genetic variability at neutral markers,quantitative trait land trait in a subdivided population under selection

Genetic variability in a subdivided population under stabilizing and diversifying selection was investigated at three levels: neutral markers, QTL coding for a trait, and the trait itself. A quantitative model with additive effects was used to link genotypes to phenotypes. No physical linkage was introduced. Using an analytical approach, we compared the diversity within deme (H(S)) and the differentiation (F(ST)) at the QTL with the genetic variance within deme (V(W)) and the differentiation (Q(ST)) for the trait. The difference between F(ST) and Q(ST) was shown to depend on the relative amounts of covariance between QTL within and between demes. Simulations were used to study the effect of selection intensity, variance of optima among demes, and migration rate for an allogamous and predominantly selfing species. Contrasting dynamics of the genetic variability at markers, QTL, and trait were observed as a function of the level of gene flow and diversifying selection. The highest discrepancy among the three levels occurred under highly diversifying selection and high gene flow. Furthermore, diversifying selection might cause substantial heterogeneity among QTL, only a few of them showing allelic differentiation, while the others behave as neutral markers.