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.     

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.     

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

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

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.     

Molecular and quantitative genetic differentiation across Europe in yellow dung flies

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

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.     

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

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

Population Structure in Daphnia Obtusa: Quantitative Genetic and Allozymic Variation

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

Adaptive population differentiation in phenology across a latitudinal gradient in European aspen (Populus tremula, L.): a comparison of neutral markers, candidate genes and phenotypic traits

A correct timing of growth cessation and dormancy induction represents a critical ecological and evolutionary trade-off between survival and growth in most forest trees (Rehfeldt et al. 1999; Horvath et al. 2003; Howe et al. 2003). We have studied the deciduous tree European Aspen (Populus tremula) across a latitudinal gradient and compared genetic differentiation in phenology traits with molecular markers. Trees from 12 different areas covering 10 latitudinal degrees were cloned and planted in two common gardens. Several phenology traits showed strong genetic differentiation and clinal variation across the latitudinal gradient, with Q(ST) values generally exceeding 0.5. This is in stark contrast to genetic differentiation at several classes of genetic markers (18 neutral SSRs, 7 SSRs located close to phenology candidate genes and 50 SNPs from five phenology candidate genes) that all showed F(ST) values around 0.015. We thus find strong evidence for adaptive divergence in phenology traits across the latitudinal gradient. However, the strong population structure seen at the quantitative traits is not reflected in underlying candidate genes. This result fit theoretical expectations that suggest that genetic differentiation at candidate loci is better described by F(ST) at neutral loci rather than by Q(ST) at the quantitative traits themselves.     

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

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

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.     

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.     

Genome scanning for interspecific differentiation between two closely related oak species [Quercus robur L. and Q. petraea (Matt.) Liebl.

Interspecific differentiation values (G(ST)) between two closely related oak species (Quercus petraea and Q. robur) were compiled across different studies with the aim to explore the distribution of differentiation at the genome level. The study was based on a total set of 389 markers (isozymes, AFLPs, SCARs, microsatellites, and SNPs) for which allelic frequencies were estimated in pairs of populations sampled throughout the sympatric distribution of the two species. The overall distribution of G(ST) values followed an L-shaped curve with most markers exhibiting low species differentiation (G(ST) < 0.01) and only a few loci reaching >10% levels. Twelve percent of the loci exhibited significant G(ST) deviations to neutral expectations, suggesting that selection contributed to species divergence. Coding regions expressed higher differentiation than noncoding regions. Among the 389 markers, 158 could be mapped on the 12 linkage groups of the existing Q. robur genetic map. Outlier loci with large G(ST) values were distributed over 9 linkage groups. One cluster of three outlier loci was found within 0.51 cM; but significant autocorrelation of G(ST) was observed at distances <2 cM. The size and distribution of genomic regions involved in species divergence are discussed in reference to hitchhiking effects and disruptive selection.     

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.     

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.     

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.     

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.     

Partitioning adaptive differentiation across a patchy landscape: shade avoidance traits in impatiens capensis

Adaptation to different habitat types across a patchy landscape may either arise independently in each patch or occur due to repeated colonization of each patch by the same specialized genotype. We tested whether open- and closed-canopy forms of Impatiens capensis, an herbaceous annual plant of eastern North America, have evolved repeatedly by comparing hierarchical measures of F(ST) estimated from AFLPs to morphological differentiation measured by Q(ST) for five pairs of populations found in open and closed habitats in five New England regions. Morphological differentiation between habitats (Q(HT)) in elongation traits was greater than marker divergence (F(HT)), suggesting adaptive differentiation. Genotypes from open- and closed-canopy habitats differed in shade avoidance traits in several population pairs, whereas patterns of AFLP differentiation suggest this differentiation does not have a single origin. These results suggest that open- and closed-canopy habitats present different selective pressures, but that the outcome of diversifying selection may differ depending on specific closed- and open-canopy habitats and on starting genetic variation. Hierarchical partitioning of F(ST) and Q(ST) makes it possible to distinguish global stabilizing selection on traits across a landscape from diversifying selection between habitat types within regions.     

Adaptive differentiation of quantitative traits in the globally distributed weed, wild radish (Raphanus raphanistrum)

Weedy species with wide geographical distributions may face strong selection to adapt to new environments, which can lead to adaptive genetic differentiation among populations. However, genetic drift, particularly due to founder effects, will also commonly result in differentiation in colonizing species. To test whether selection has contributed to trait divergence, we compared differentiation at eight microsatellite loci (measured as F(ST)) to differentiation of quantitative floral and phenological traits (measured as Q(ST)) of wild radish (Raphanus raphanistrum) across populations from three continents. We sampled eight populations: seven naturalized populations and one from its native range. By comparing estimates of Q(ST) and F(ST), we found that petal size was the only floral trait that may have diverged more than expected due to drift alone, but inflorescence height, flowering time, and rosette formation have greatly diverged between the native and nonnative populations. Our results suggest the loss of a rosette and the evolution of early flowering time may have been the key adaptations enabling wild radish to become a major agricultural weed. Floral adaptation to different pollinators does not seem to have been as necessary for the success of wild radish in new environments.     

Thermal adaptation in the fungal pathogen Mycosphaerella graminicola

Genetic differentiation in thermal adaptation can result from a trade-off between the performance of organisms across different temperatures or from the accumulation of deleterious mutations. In this experiment, we assayed thermal sensitivity of 138 genetically distinct Mycosphaerella graminicola isolates sampled from five host populations in four locations under two temperature regimes (22 and 15 °C) and found significant differences in growth rate and response to temperature among populations. On average, genetic differentiation accounted for more than 50% of phenotypic variation in thermal adaptation while plasticity contributed less than a quarter of phenotypic variation. Populations originating from warm places performed better under the high-temperature regime and had a larger positive response to increasing temperature. Pairwise population differentiation (Q(ST) ) in temperature sensitivity, measured by taking the ratio of growth rates at 22 to 15 °C, was positively and significantly correlated to the pairwise difference in annual mean temperature at the collection sites. Because overall Q(ST) in temperature sensitivity was significantly higher than overall G(ST) in neutral restriction fragment length polymorphism loci, we believe that the primary mechanism underlying this thermal adaptation is antagonistic pleiotropy. Our results indicate that temperature sensitivity is a better indicator of thermal adaptation than growth rate at individual temperatures.