Ecological genetics of range size variation in Boechera spp. (Brassicaceae)

Many taxonomic groups contain both rare and widespread species, which indicates that range size can evolve quickly. Many studies have compared molecular genetic diversity, plasticity, or phenotypic traits between rare and widespread species; however, a suite of genetic attributes that unites rare species remains elusive. Here, using two rare and two widespread Boechera (Brassicaceae) species, we conduct a simultaneous comparison of quantitative trait diversity, genetic diversity, and population structure among species with highly divergent range sizes. Consistent with previous studies, we do not find strong associations between range size and within‐population genetic diversity. In contrast, we find that both the degree of phenotypic plasticity and quantitative trait structure (Q_ST) were positively correlated with range size. We also found higher F_ST: Q_ST ratios in rare species, indicative of either a greater response to stabilizing selection or a lack of additive genetic variation. While widespread species occupy more ecological and climactic space and have diverged at both traits and markers, rare species display constrained levels of population differentiation and phenotypic plasticity. Combined, our results provide evidence for a specialization–generalization trade‐off across three orders of magnitude of range size variation in the ecological model genus, Boechera.     

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.     

Population size is weakly related to quantitative genetic variation and trait differentiation in a stream fish

How population size influences quantitative genetic variation and differentiation among natural, fragmented populations remains unresolved. Small, isolated populations might occupy poor quality habitats and lose genetic variation more rapidly due to genetic drift than large populations. Genetic drift might furthermore overcome selection as population size decreases. Collectively, this might result in directional changes in additive genetic variation (V_A) and trait differentiation (Q_ST) from small to large population size. Alternatively, small populations might exhibit larger variation in V_A and Q_ST if habitat fragmentation increases variability in habitat types. We explored these alternatives by investigating V_A and Q_ST using nine fragmented populations of brook trout varying 50‐fold in census size N (179–8416) and 10‐fold in effective number of breeders, N_b (18–135). Across 15 traits, no evidence was found for consistent differences in V_A and Q_ST with population size and almost no evidence for increased variability of V_A or Q_ST estimates at small population size. This suggests that (i) small populations of some species may retain adaptive potential according to commonly adopted quantitative genetic measures and (ii) populations of varying sizes experience a variety of environmental conditions in nature, however extremely large studies are likely required before any firm conclusions can be made.     

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.     

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.     

Natural selection and neutral evolutionary processes contribute to genetic divergence in leaf traits across a precipitation gradient in the tropical oak Quercus oleoides

The impacts of drought are expanding worldwide as a consequence of climate change. However, there is still little knowledge of how species respond to long‐term selection in seasonally dry ecosystems. In this study, we used Q_ST‐F_ST comparisons to investigate (i) the role of natural selection on population genetic differentiation for a set of functional traits related to drought resistance in the seasonally dry tropical oak Quercus oleoides and (ii) the influence of water availability at the site of population origin and in experimental treatments on patterns of trait divergence. We conducted a thorough phenotypic characterization of 1912 seedlings from ten populations growing in field and greenhouse common gardens under replicated watering treatments. We also genotyped 218 individuals from the same set of populations using eleven nuclear microsatellites. Q_ST distributions for leaf lamina area, specific leaf area, leaf thickness and stomatal pore index were higher than F_ST distribution. Results were consistent across growth environments. Genetic differentiation among populations for these functional traits was associated with the index of moisture at the origin of the populations. Together, our results suggest that drought is an important selective agent for Q. oleoides and that differences in length and severity of the dry season have driven the evolution of genetic differences in functional traits.     

Low levels of quantitative and molecular genetic differentiation among natural populations of Medicago ciliaris Kroch. (Fabaceae) of different Tunisian eco-geographical origin

In this paper, we analyze the genetic variability in four Tunisian natural populations of Medicago ciliaris using 19 quantitative traits and six polymorphic microsatellite loci. We investigated the amplification transferability of 30 microsatellites developed in the model legume M. truncatula to M. ciliaris. Results revealed that about 56.66% of analyzed markers are valuable genetic markers for M. ciliaris. The most genetic diversity at quantitative traits and microsatellite loci was found to occur within populations (>80%). Low differentiations among populations at quantitative traits Q _ST  = 0.146 and molecular markers F _ST  = 0.18 were found. The majority of measured traits exhibited no significant difference in the level of Q _ST and F _ST . Furthermore, significant correlations established between these traits and eco-geographical factors suggested that natural selection should be invoked to explain the level of phenotypic divergence among populations rather than drift. There was no significant correlation between population differentiation at quantitative traits and molecular markers. Significant spatial genetic structure consistent with models of isolation by distance was detected within all studied populations. The site-of-origin environmental factors explain about 9.07% of total phenotypic genetic variation among populations. The eco-geographical factors that influence more the variation of measured traits among populations are the soil texture and altitude. Nevertheless, there were no consistent pattern of associations between gene diversity (He) and environmental factors.     

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.     

Comparison of quantitative and molecular variation in agroforestry populations of the shea tree (Vitellaria paradoxa C.F. Gaertn) in Mali

In this study we investigated the within- and between-population genetic variation using microsatellite markers and quantitative traits of the shea tree, Vitellaria paradoxa, an important agroforestry tree species of the Sudano–Sahelian region in Africa. Eleven populations were sampled across Mali and in northern Côte d’Ivoire. Leaf size and form and growth traits were measured in a progeny test at the nursery stage. Eight microsatellites were used to assess neutral genetic variation. Low levels of heterozygosity were recorded (1.6–3.0 alleles/locus; H_E = 0.25–0.42) and the fixation index (F_IS = −0.227–0.186) was not significantly different from zero suggesting that Hardy–Weinberg equilibrium is encountered in all populations sampled. Quantitative traits exhibited a strong genetic variation between populations and between families within populations. The degree of population differentiation of the quantitative traits (Q_ST = 0.055–0.283, Q_STmean = 0.189) strongly exceeds that in eight microsatellite loci (F_ST = −0.011–0.142, F_STmean = 0.047). Global and pairwise F_ST values were very low and not significantly different from zero suggesting agroforestry practices are amplifying gene flow (Nm = 5.07). The population means for quantitative traits and the rainfall variable were not correlated, showing variation was not linked with this climatic cline. It is suggested that this marked differentiation for quantitative traits, independent of environmental clines and despite a high gene flow, is a result of local adaptation and human selection of shea trees. This process has induced high linkage disequilibrium between underlying loci of polygenic characters.     

Size differentiation in Finnish house sparrows follows Bergmann's rule with evidence of local adaptation

Bergmann's rule predicts that individuals are larger in more poleward populations and that this size gradient has an adaptive basis. Hence, phenotypic divergence in size traits between populations (P_ST) is expected to exceed the level of divergence by drift alone (F_ST). We measured 16 skeletal traits, body mass and wing length in 409 male and 296 female house sparrows Passer domesticus sampled in 12 populations throughout Finland, where the species has its northernmost European distributional margin. Morphometric differentiation across populations (P_ST) was compared with differentiation in 13 microsatellites (F_ST). We find that twelve traits phenotypically diverged more than F_ST in both sexes, and an additional two traits diverged in males. The phenotypic divergence exceeded F_ST in several traits to such a degree that findings were robust also to strong between‐population environmental effects. Divergence was particularly strong in dimensions of the bill, making it a strong candidate for the study of adaptive molecular genetic divergence. Divergent traits increased in size in more northern populations. We conclude that house sparrows show evidence of an adaptive latitudinal size gradient consistent with Bergmann's rule on the modest spatial scale of ca. 600 km.     

Genetic variation of early height growth traits at the xeric limits of Austrocedrus chilensis (Cupressaceae)

The study of the genetic variation of early height growth traits in seedlings helps to predict the possible outcomes of tree populations in the face of climate change. Second‐year height growth of 10 geographically marginal populations of Patagonian cypress (Austrocedrus chilensis (D. Don) Pic. Ser. et Bizzarri) (Cupressaceae) was characterized under greenhouse conditions. Variation among and within an average of 15 open‐pollinated families (comprising 21 seedlings per family) for each population was analysed for six size and timing traits obtained from fitted Boltzmann growth curves. The among‐family and among‐population variances were 4.03% and 2.74% of the total phenotypic variation, while the residual variance was 84.57% on average. Genetic differentiation among populations was low, except for the maximum growth rate (Q_ST = 0.35) and for growth initiation (Q_ST = 1). For most traits, genetic variation and heritability were variable across populations, except for growth initiation, which showed in general null intra‐population levels of genetic variance. Although no direct associations were found between the additive genetic variation and latitude or altitude, the north range of the distribution was more variable for the pool of the analysed traits. Although most extreme‐marginal populations of A. chilensis would be very limited in their ability to evolve if climate in north‐west Patagonia turns drier and warmer, their long‐term persistence could largely rely on a phenotypic diversification strategy.     

Are QST–FST comparisons for natural populations meaningful?

Comparisons between putatively neutral genetic differentiation amongst populations, F_ST, and quantitative genetic variation, Q_ST, are increasingly being used to test for natural selection. However, we find that approximately half of the comparisons that use only data from wild populations confound phenotypic and genetic variation. We urge the use of a clear distinction between narrow‐sense Q_ST, which can be meaningfully compared with F_ST, and phenotypic divergence measured between populations, P_ST, which is inadequate for comparisons in the wild. We also point out that an unbiased estimate of Q_ST can be found using the so‐called ‘animal model’ of quantitative genetics.     

Local selection modifies phenotypic divergence among Rana temporaria populations in the presence of gene flow

In ectotherms, variation in life history traits among populations is common and suggests local adaptation. However, geographic variation itself is not a proof for local adaptation, as genetic drift and gene flow may also shape patterns of quantitative variation. We studied local and regional variation in means and phenotypic plasticity of larval life history traits in the common frog Rana temporaria using six populations from central Sweden, breeding in either open‐canopy or partially closed‐canopy ponds. To separate local adaptation from genetic drift, we compared differentiation in quantitative genetic traits (Q_ST) obtained from a common garden experiment with differentiation in presumably neutral microsatellite markers (F_ST). We found that R. temporaria populations differ in means and plasticities of life history traits in different temperatures at local, and in F_ST at regional scale. Comparisons of differentiation in quantitative traits and in molecular markers suggested that natural selection was responsible for the divergence in growth and development rates as well as in temperature‐induced plasticity, indicating local adaptation. However, at low temperature, the role of genetic drift could not be separated from selection. Phenotypes were correlated with forest canopy closure, but not with geographical or genetic distance. These results indicate that local adaptation can evolve in the presence of ongoing gene flow among the populations, and that natural selection is strong in this system.     

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.     

Parasite‐mediated selection in a natural metapopulation of Daphnia magna

Parasite‐mediated selection varying across time and space in metapopulations is expected to result in host local adaptation and the maintenance of genetic diversity in disease‐related traits. However, nonadaptive processes like migration and extinction‐(re)colonization dynamics might interfere with adaptive evolution. Understanding how adaptive and nonadaptive processes interact to shape genetic variability in life‐history and disease‐related traits can provide important insights into their evolution in subdivided populations. Here we investigate signatures of spatially fluctuating, parasite‐mediated selection in a natural metapopulation of Daphnia magna. Host genotypes from infected and uninfected populations were genotyped at microsatellite markers, and phenotyped for life‐history and disease traits in common garden experiments. Combining phenotypic and genotypic data a Q_ST–F_ST‐like analysis was conducted to test for signatures of parasite mediated selection. We observed high variation within and among populations for phenotypic traits, but neither an indication of host local adaptation nor a cost of resistance. Infected populations have a higher gene diversity (Hs) than uninfected populations and Hs is strongly positively correlated with fitness. These results suggest a strong parasite effect on reducing population level inbreeding. We discuss how stochastic processes related to frequent extinction‐(re)colonization dynamics as well as host and parasite migration impede the evolution of resistance in the infected populations. We suggest that the genetic and phenotypic patterns of variation are a product of dynamic changes in the host gene pool caused by the interaction of colonization bottlenecks, inbreeding, immigration, hybrid vigor, rare host genotype advantage and parasitism. Our study highlights the effect of the parasite in ameliorating the negative fitness consequences caused by the high drift load in this metapopulation.     

Post‐glacial colonization routes coincide with a life‐history breakpoint along a latitudinal gradient

Although adaptive divergence along environmental gradients has repeatedly been demonstrated, the role of post‐glacial colonization routes in determining phenotypic variation along gradients has received little attention. Here, we used a hierarchical Q_ST–F_ST approach to separate the roles of adaptive and neutral processes in shaping phenotypic variation in moor frog (Rana arvalis) larval life histories along a 1,700 km latitudinal gradient across northern Europe. This species has colonized Scandinavia via two routes with a contact zone in northern Sweden. By using neutral SNP and common garden phenotypic data from 13 populations at two temperatures, we showed that most of the variation along the gradient occurred between the two colonizing lineages. We found little phenotypic divergence within the lineages; however, all phenotypic traits were strongly diverged between the southern and northern colonization routes, with higher growth and development rates and larger body size in the north. The Q_ST estimates between the colonization routes were four times higher than F_ST, indicating a prominent role for natural selection. Q_ST within the colonization routes did not generally differ from F_ST, but we found temperature‐dependent adaptive divergence close to the contact zone. These results indicate that lineage‐specific variation can account for much of the adaptive divergence along a latitudinal gradient.     

The effects of selection,drift and genetic variation on life‐history trait divergence among insular populations of natterjack toad,Bufo calamita

Although loss of genetic variation is frequently assumed to be associated with loss of adaptive potential, only few studies have examined adaptation in populations with little genetic variation. On the Swedish west coast, the northern fringe populations of the natterjack toad Bufo calamita inhabit an atypical habitat consisting of offshore rock islands. There are strong among‐population differences in the amount of neutral genetic variation, making this system suitable for studies on mechanisms of trait divergence along a gradient of within‐population genetic variation. In this study, we examined the mechanisms of population divergence using Q_ST–F_ST comparisons and correlations between quantitative and neutral genetic variation. Our results suggest drift or weak stabilizing selection across the six populations included in this study, as indicated by low Q_ST–F_ST values, lack of significant population × temperature interactions and lack of significant differences among the islands in breeding pond size. The six populations included in this study differed in both neutral and quantitative genetic variation. Also, the correlations between neutral and quantitative genetic variation tended to be positive, however, the relatively small number of populations prevents any strong conclusions based on these correlations. Contrary to the majority of Q_ST–F_ST comparisons, our results suggest drift or weak stabilizing selection across the examined populations. Furthermore, the low heritability of fitness‐related traits may limit evolutionary responses in some of the populations.     

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.     

Separating Effects of Gene Flow and Natural Selection along an Environmental Gradient

Genetic differentiation along environmental clines is often observed as a result of interplay between gene flow and natural selection. In order to understand the relative roles of these processes in shaping this differentiation, we designed a study in which we used two approaches that have not previously been combined, the Q _ST–F _ST comparison and crossbreeding. We examined (1) interpopulation phenotypic and genetic (AFLP) variation, and (2) performance of interpopulation hybrids in a common annual Senecio glaucus. Fitness of interpopulation hybrids (F1 and F2) was tested under simulated population natural conditions in terms of aridity and analyzed for a relationship with (1) spatial distance and (2) environmental differences (amount of annual rainfall). While phenotypic variation corresponded to the clinal changes in aridity along population locations, viz. narrower and longer leaves, longer leaf outgrowths and advanced flowering in more arid environments, the F _ST < 0.1 calculated from AFLP data suggested intensive interpopulation gene flow, with little if any contribution of genetic drift. Performance of hybrids in simulated natural environments revealed heterosis in F1, but a hybrid breakdown in F2 generation. These effects were related to both the spatial distance between hybrid parents and their population rainfall differences. The detected clinal phenotypic variation and outbreeding depression in F2 strongly support presence of aridity-induced clinal natural selection, which is matched by the observed Q _ST ≫ F _ST. From this we conclude that Q _ST–F _ST comparison can detect effect of diversifying selection when patterns of phenotypic variation across sampled locations can be reliably predicted from environmental variation.