F(ST) and Q(ST) under neutrality

Ethanol resistance in Drosophila melanogaster has increased in parallel cold-adapted populations and shows a variable genetic architecture within and between populations

Understanding the genetic properties of adaptive trait evolution is a fundamental crux of biological inquiry that links molecular processes to biological diversity. Important uncertainties persist regarding the genetic predictability of adaptive trait change, the role of standing variation, and whether adaptation tends to result in the fixation of favored variants. Here, we use the recurrent evolution of enhanced ethanol resistance in Drosophila melanogaster during this species' worldwide expansion as a promising system to add to our understanding of the genetics of adaptation. We find that elevated ethanol resistance has evolved at least three times in different cooler regions of the species' modern range-not only at high latitude but also in two African high-altitude regions. Applying a bulk segregant mapping framework, we find that the genetic architecture of ethanol resistance evolution differs substantially not only between our three resistant populations, but also between two crosses involving the same European population. We then apply population genetic scans for local adaptation within our quantitative trait locus regions, and we find potential contributions of genes with annotated roles in spindle localization, membrane composition, sterol and alcohol metabolism, and other processes. We also apply simulation-based analyses that confirm the variable genetic basis of ethanol resistance and hint at a moderately polygenic architecture. However, these simulations indicate that larger-scale studies will be needed to more clearly quantify the genetic architecture of adaptive evolution and to firmly connect trait evolution to specific causative loci.

Evidence of local adaptation despite strong drift in a Neotropical patchily distributed bromeliad

Both genetic drift and divergent selection are predicted to be drivers of population differentiation across patchy habitats, but the extent to which these forces act on natural populations to shape traits is strongly affected by species' ecological features. In this study, we infer the genomic structure of Pitcairnia lanuginosa, a widespread herbaceous perennial plant with a patchy distribution. We sampled populations in the Brazilian Cerrado and the Central Andean Yungas and discovered and genotyped SNP markers using double-digest restriction-site associated DNA sequencing. In addition, we analyzed ecophysiological traits obtained from a common garden experiment and compared patterns of phenotypic and genetic divergence (P_ST-F_ST comparisons) in a subset of populations from the Cerrado. Our results from molecular analyses pointed to extremely low genetic diversity and a remarkable population differentiation, supporting a major role of genetic drift. Approximately 0.3% of genotyped SNPs were flagged as differentiation outliers by at least two distinct methods, and Bayesian generalized linear mixed models revealed a signature of isolation by environment in addition to isolation by distance for high-differentiation outlier SNPs among the Cerrado populations. P_ST-F_ST comparisons suggested divergent selection on two ecophysiological traits linked to drought tolerance. We showed that these traits vary among populations, although without any particular macro-spatial pattern, suggesting local adaptation to differences in micro-habitats. Our study shows that selection might be a relevant force, particularly for traits involved in drought stress, even for populations experiencing strong drift, which improves our knowledge on eco-evolutionary processes acting on non-continuously distributed species.

Bi-directional Selection in Upland Rice Leads to Its Adaptive Differentiation from Lowland Rice in Drought Resistance and Productivity

Drought resistance is required in rice breeding to address the challenge of frequent droughts. However, the evolutionary mechanism of rice drought resistance is not fully understood. We investigated the genetic differentiation between upland and lowland rice domesticated in agro-ecosystems with contrasting water-soil conditions using genome-wide SNPs. We estimated morphological differences among upland and lowland rice in drought resistance and productivity through common garden experiments. Upland rice had better drought resistance but poorer productivity. The negative correlations between traits of drought resistance and productivity are attributed to the underlying genetic trade-offs through tight linkages (e.g., DCA1 and OsCesA7) or pleiotropic effects (e.g., LAX1). The genetic trade-offs are common and greatly shape the evolution of drought resistance in upland rice. In genomic regions associated with both productivity and drought resistance, signs of balancing selection were detected in upland rice, while signs of directional selection were detected in lowland rice, potentially contributing to their adaptive differentiation. Signs of balancing selection in upland rice resulted from bi-directional selection during its domestication in drought-prone upland agro-ecosystems. Using genome-wide association analysis, we identified several valuable quantitative trait loci associated with drought resistance, for which highly differentiated genes should be considered candidates. Bi-directional selection breaking tight linkages by accumulating recombination events would be applicable in breeding water-saving and drought-resistance rice.

Role of selection and gene flow in population differentiation at the edge vs. interior of the species range differing in climatic conditions

Evaluating the relative importance of neutral and adaptive processes as determinants of population differentiation across environments is a central theme of evolutionary biology. We applied the QST-FST comparison flanked by a direct test for local adaptation to infer the role of climate-driven selection and gene flow in population differentiation of an annual grass Avena sterilis in two distinct parts of the species range, edge and interior, which represent two globally different climates, desert and Mediterranean. In a multiyear reciprocal transplant experiment, the plants of desert and Mediterranean origin demonstrated home advantage, and population differentiation in several phenotypic traits related to reproduction exceeded neutral predictions, as determined by comparisons of QST values with theoretical FST distributions. Thus, variation in these traits likely resulted from local adaptation to desert and Mediterranean environments. The two separate common garden experiments conducted with different experimental design revealed that two population comparisons, in contrast to multi-population comparisons, are likely to detect population differences in virtually every trait, but many of these differences reflect effects of local rather than regional environment. We detected a general reduction in neutral (SSR) genetic variation but not in adaptive quantitative trait variation in peripheral desert as compared with Mediterranean core populations. On the other hand, the molecular data indicated intensive gene flow from the Mediterranean core towards desert periphery. Although species range position in our study (edge vs. interior) was confounded with climate (desert vs. Mediterranean), the results suggest that the gene flow from the species core does not have negative consequences for either performance of the peripheral plants or their adaptive potential.

Decanalization of wing development accompanied the evolution of large wings in high-altitude Drosophila

In higher organisms, the phenotypic impacts of potentially harmful or beneficial mutations are often modulated by complex developmental networks. Stabilizing selection may favor the evolution of developmental canalization--that is, robustness despite perturbation--to insulate development against environmental and genetic variability. In contrast, directional selection acts to alter the developmental process, possibly undermining the molecular mechanisms that buffer a trait's development, but this scenario has not been shown in nature. Here, we examined the developmental consequences of size increase in highland Ethiopian Drosophila melanogaster. Ethiopian inbred strains exhibited much higher frequencies of wing abnormalities than lowland populations, consistent with an elevated susceptibility to the genetic perturbation of inbreeding. We then used mutagenesis to test whether Ethiopian wing development is, indeed, decanalized. Ethiopian strains were far more susceptible to this genetic disruption of development, yielding 26 times more novel wing abnormalities than lowland strains in F2 males. Wing size and developmental perturbability cosegregated in the offspring of between-population crosses, suggesting that genes conferring size differences had undermined developmental buffering mechanisms. Our findings represent the first observation, to our knowledge, of morphological evolution associated with decanalization in the same tissue, underscoring the sensitivity of development to adaptive change.

Sympatric divergence and clinal variation in multiple coloration traits of Ficedula flycatchers

Geographic variation in phenotypes plays a key role in fundamental evolutionary processes such as local adaptation, population differentiation and speciation, but the selective forces behind it are rarely known. We found support for the hypothesis that geographic variation in plumage traits of the pied flycatcher Ficedula hypoleuca is explained by character displacement with the collared flycatcher Ficedula albicollis in the contact zone. The plumage traits of the pied flycatcher differed strongly from the more conspicuous collared flycatcher in a sympatric area but increased in conspicuousness with increasing distance to there. Phenotypic differentiation (PST ) was higher than that in neutral genetic markers (FST ), and the effect of geographic distance remained when statistically controlling for neutral genetic differentiation. This suggests that a cline created by character displacement and gene flow explains phenotypic variation across the distribution of this species. The different plumage traits of the pied flycatcher are strongly to moderately correlated, indicating that they evolve non-independently from each other. The flycatchers provide an example of plumage patterns diverging in two species that differ in several aspects of appearance. The divergence in sympatry and convergence in allopatry in these birds provide a possibility to study the evolutionary mechanisms behind the highly divergent avian plumage patterns.

Spatially varying selection shapes life history clines among populations of Drosophila melanogaster from sub-Saharan Africa

Clines in life history traits, presumably driven by spatially varying selection, are widespread. Major latitudinal clines have been observed, for example, in Drosophila melanogaster, an ancestrally tropical insect from Africa that has colonized temperate habitats on multiple continents. Yet, how geographic factors other than latitude, such as altitude or longitude, affect life history in this species remains poorly understood. Moreover, most previous work has been performed on derived European, American and Australian populations, but whether life history also varies predictably with geography in the ancestral Afro-tropical range has not been investigated systematically. Here, we have examined life history variation among populations of D. melanogaster from sub-Saharan Africa. Viability and reproductive diapause did not vary with geography, but body size increased with altitude, latitude and longitude. Early fecundity covaried positively with altitude and latitude, whereas lifespan showed the opposite trend. Examination of genetic variance-covariance matrices revealed geographic differentiation also in trade-off structure, and QST -FST analysis showed that life history differentiation among populations is likely shaped by selection. Together, our results suggest that geographic and/or climatic factors drive adaptive phenotypic differentiation among ancestral African populations and confirm the widely held notion that latitude and altitude represent parallel gradients.

So far away, yet so close: strong genetic structure in Homonota uruguayensis (Squamata, Phyllodactylidae), a species with restricted geographic distribution in the Brazilian and Uruguayan Pampas

The Pampas is a biologically rich South American biome, but is poorly represented in phylogeographic studies. While the Pleistocene glacial cycles may have affected the evolutionary history of species distributed in forested biomes, little is known about their effects on the habitats that remained stable through glacial cycles. The South American Pampas have been covered by grasslands during both glacial and interglacial periods and therefore represent an interesting system to test whether the genetic structure in such environments is less pronounced. In this study, we sampled Pampean populations of Homonota uruguayensis from Southern Brazil and Uruguay to assess the tempo and mode of population divergence, using both morphological measurements and molecular markers. Our results indicate that, in spite of its narrow geographic distribution, populations of H. uruguayensis show high levels of genetic structure. We found four major well-supported mtDNA clades with strong geographic associations. Estimates of their divergence times fell between 3.16 and 1.82 million years before the present. Populations from the central portion of the species distribution, on the border between Uruguay and Brazil, have high genetic diversity and may have undergone a population expansion approximately 250,000 years before the present. The high degree of genetic structure is reflected in the analyses of morphological characters, and most individuals could be correctly assigned to their parental population based on morphology alone. Finally, we discuss the biogeographic and conservation implications of these findings.

Genetic differentiation revealed by selective loci of drought-responding EST-SSRs between upland and lowland rice in China

Upland and lowland rice (Oryza sativa L.) represent two of the most important rice ecotypes adapted to ago-ecosystems with contrasting soil-water conditions. Upland rice, domesticated in the water-limited environment, contains valuable drought-resistant characters that can be used in water-saving breeding. Knowledge about the divergence between upland and lowland rice will provide valuable cues for the evolution of drought-resistance in rice. Genetic differentiation between upland and lowland rice was explored by 47 Simple Sequence Repeats (SSRs) located in drought responding expressed sequence tags (ESTs) among 377 rice landraces. The morphological traits of drought-resistance were evaluated in the field experiments. Different outlier loci were detected in the japonica and indica subspecies, respectively. Considerable genetic differentiation between upland and lowland rice on these outlier loci was estimated in japonica (Fst = 0.258) and indica (Fst = 0.127). Furthermore, populations of the upland and lowland ecotypes were clustered separately on these outlier loci. A significant correlation between genetic distance matrices and the dissimilarity matrices of drought-resistant traits was determined, indicating a certain relationship between the upland-lowland rice differentiation and the drought-resistance. Divergent selections occur between upland and lowland rice on the drought-resistance as the Qsts of some drought-resistant traits are significantly higher than the neutral Fst. In addition, the upland- and lowland-preferable alleles responded differently among ecotypes or allelic types under osmotic stress. This shows the evolutionary signature of drought resistance at the gene expression level. The findings of this study can strengthen our understanding of the evolution of drought-resistance in rice with significant implications in the improvement of rice drought-resistance.

Variation and signatures of selection on the human face

There has been much debate about why humans throughout the world differ in facial form. Previous studies of human skull morphology found levels of among-population differentiation that were comparable to those of neutral genetic markers, suggesting that genetic drift (neutral processes) played an important role in influencing facial differentiation. However, variation in soft-tissue morphology has not been studied in detail. In this study, we analyzed high-resolution 3D images of soft-tissue facial form in four Eurasian populations: Han Chinese, Tibetans, Uyghur and Europeans. A novel method was used to establish a high-density alignment across all of the faces, allowing facial diversity to be examined at an unprecedented resolution. These data exhibit signatures of population structure and history. However, among-population differentiation was higher for soft-tissue facial form than for genome-wide genetic loci, and high-resolution analyses reveal that the nose, brow area and cheekbones exhibit particularly strong signals of differentiation (Qst estimates: 0.3-0.8) between Europeans and Han Chinese. Our results suggest that local adaptation and/or sexual selection have been important in shaping human soft-tissue facial morphology.

Geometric morphometric analysis of Colombian Anopheles albimanus (Diptera: Culicidae) reveals significant effect of environmental factors on wing traits and presence of a metapopulation

Anopheles albimanus is a major malaria mosquito vector in Colombia. In the present study, wing variability (size and shape) in An. albimanus populations from Colombian Maracaibo and Chocó bio-geographical eco-regions and the relationship of these phenotypic traits with environmental factors were evaluated. Microsatellite and morphometric data facilitated a comparison of the genetic and phenetic structure of this species. Wing size was influenced by elevation and relative humidity, whereas wing shape was affected by these two variables and also by rainfall, latitude, temperature and eco-region. Significant differences in mean shape between populations and eco-regions were detected, but they were smaller than those at the intra-population level. Correct assignment based on wing shape was low at the population level (<58%) and only slightly higher (>70%) at the eco-regional level, supporting the low population structure inferred from microsatellite data. Wing size was similar among populations with no significant differences between eco-regions. Population relationships in the genetic tree did not agree with those from the morphometric data; however, both datasets consistently reinforced a panmictic population of An. albimanus. Overall, site-specific population differentiation is not strongly supported by wing traits or genotypic data. We hypothesize that the metapopulation structure of An. albimanus throughout these Colombian eco-regions is favoring plasticity in wing traits, a relevant characteristic of species living under variable environmental conditions and colonizing new habitats.

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 (QST) obtained from a common garden experiment was compared with differentiation at putatively neutral microsatellite markers (FST) in seven populations of P. sinense. QST estimates were derived using a Bayesian hierarchical variance component method.

KEY RESULTS

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

CONCLUSIONS

Despite the prevalent empirical evidence for QST > FST, 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 QST relative to FST, 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, QST estimates were of the same magnitude as FST, 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.

Among- and within-population variation in flowering time of Iberian Arabidopsis thaliana estimated in field and glasshouse conditions

The study of the evolutionary and population genetics of quantitative traits requires the assessment of within- and among-population patterns of variation. We carried out experiments including eight Iberian Arabidopsis thaliana populations (10 individuals per population) in glasshouse and field conditions. We quantified among- and within-population variation for flowering time and for several field life-history traits. Individuals were genotyped with microsatellites, single nucleotide polymorphisms and four well-known flowering genes (FRI, FLC, CRY2 and PHYC). Phenotypic and genotypic data were used to conduct Q(ST)-F(ST) comparisons. Life-history traits varied significantly among- and within-populations. Flowering time also showed substantial within- and among-population variation as well as significant genotype × environment interactions among the various conditions. Individuals bearing FRI truncations exhibited reduced recruitment in field conditions and differential flowering time behavior across experimental conditions, suggesting that FRI contributes to the observed significant genotype × environment interactions. Flowering time estimated in field conditions was the only trait showing significantly higher quantitative genetic differentiation than neutral genetic differentiation values. Overall, our results show that these A. thaliana populations are genetically more differentiated for flowering time than for neutral markers, suggesting that flowering time is likely to be under divergent selection.

Q(ST) < F(ST) As a signature of canalization

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

Genetic basis of adaptation in Arabidopsis thaliana: local adaptation at the seed dormancy QTL DOG1

Local adaptation provides an opportunity to study the genetic basis of adaptation and investigate the allelic architecture of adaptive genes. We study delay of germination 1 (DOG1), a gene controlling natural variation in seed dormancy in Arabidopsis thaliana and investigate evolution of dormancy in 41 populations distributed in four regions separated by natural barriers. Using F(ST) and Q(ST) comparisons, we compare variation at DOG1 with neutral markers and quantitative variation in seed dormancy. Patterns of genetic differentiation among populations suggest that the gene DOG1 contributes to local adaptation. Although Q(ST) for seed dormancy is not different from F(ST) for neutral markers, a correlation with variation in summer precipitation supports that seed dormancy is adaptive. We characterize dormancy variation in several F(2) -populations and show that a series of functionally distinct alleles segregate at the DOG1 locus. Theoretical models have shown that the number and effect of alleles segregatin at quantitative trait loci (QTL) have important consequences for adaptation. Our results provide support to models postulating a large number of alleles at quantitative trait loci involved in adaptation.

Adaptive genetic differentiation in a predominantly self-pollinating species analyzed by transplanting into natural environment, crossbreeding and Q(ST)-F(ST) test

• Both genetic drift and natural selection result in genetic/phenotypic differentiation over space. I analyzed the role of local adaptation in the genetic differentiation of populations of the annual grass Hordeum spontaneum sampled along an aridity gradient. • The study included the introduction of plants having desert vs nondesert origin into natural (desert) environment, analysis of population differentiation in allozymes and random amplified polymorphic DNA (RAPD) markers vs phenotypic traits (Q(ST) -F(ST) comparison), and planting interpopulation hybrids under simulated desert conditions in a glasshouse. • The results of the home advantage test, Q(ST) -F(ST) comparison and crossbreeding were consistent with local adaptation; that is, that differentiation of the desert plants from plants of nondesert origin in phenotypic traits was adaptive, giving them home advantage. Each method used provided additional, otherwise unavailable, information, meaning that they should be viewed as complementary rather than alternative approaches. • Gene flow from adjacent populations (i.e. populations experiencing the desert environment) via seeds (but not pollen) had a positive effect on fitness by enhancing natural selection and counteracting drift. At the same time, the effect of genes from the species distributional core (nondesert plants) by either seed or pollen had a negative fitness effect despite its enriching effect on neutral diversity. The pattern of outbreeding depression observed in interpopulation hybrids (F(1) ) and their segregating progeny (F(2) ) was inconsistent with underdominance, but indicated the presence of additive, dominance and epistatic effects.

Does selection or genetic drift explain geographic differentiation of morphological characters in house sparrows Passer domesticus?

Understanding the relative influence of genetic drift and selection is fundamental in evolutionary biology. The theory of neutrality predicts that the genetic differentiation of a quantitative trait (QST) equals the genetic differentiation at neutral molecular markers (FST) if the quantitative trait has not been under selection. Thus, the relative magnitude of observed QST and expected QST under neutral expectations suggests the importance of selection and genetic drift for any observed phenotypic divergence. Because QST is based on additive genetic variance, estimating QST based on phenotypic measurements is problematic due to unknown environmental effects. To account for this, we used a model where the environmental component was allowed to vary when estimating QST. The model was used on data from 14 house sparrow (Passer domesticus) populations in Norway. In accordance with the significant phenotypic inter-population differences our analyses suggested that directional selection may have favoured different optimal phenotypes for some morphological traits across populations. In particular, different body mass and male ornamental phenotypes seemed to have been favoured. The conclusions are, however, dependent on assumptions regarding the proportion of the observed inter-population variation that is due to additive genetic differences, showing the importance of collecting such information in natural populations. By estimating QST, allowing the additive genetic proportion of phenotypic inter-population variation to vary, and by making use of recent statistical methods to compare observed QST with neutral expectations, we can use data that are relatively easy to collect to identify adaptive variation in natural populations.

Whither Pst? The approximation of Qst by Pst in evolutionary and conservation biology

Local adaptation through natural selection can be inferred in case the additive genetic divergence in a quantitative trait across populations (Q(st)) exceeds the neutral expectation based on differentiation of neutral alleles across these populations (e.g. F(st)). As such, measuring Q(st) in relation to neutral differentiation presents a first-line investigation applicable in evolutionary biology (selection on functional genes) and conservation biology (identification of locally adapted coding genes). However, many species, especially those in need of conservation actions, are not amenable for the kind of breeding design required to estimate either narrow- or broad-sense Q(st). In such cases, Q(st) has been approximated by the phenotypic divergence in a trait across populations (P(st)). I here argue that the critical aspect for how well P(st) approximates Q(st) depends on the extent that additive genetic effects determine variation between populations relative to within populations. I review how the sensitivity of conclusions regarding local adaptation based on P(st) have been evaluated in the literature and find that many studies make a anticonservative null assumption in estimating P(st) and/or use a nonconservative approach to explore sensitivity of their conclusions. Data from two studies that have provided a second, independent assessment of selection in their system suggest that P(st)-F(st) comparisons should be interpreted very conservatively. I conclude with recommendations for improving the robustness of the inferences drawn from comparing P(st) with neutral differentiation.

If F(ST) does not measure neutral genetic differentiation, then comparing it with Q(ST) is misleading. Or is it?

The comparison between neutral genetic differentiation (F(ST) ) and quantitative genetic differentiation (Q(ST) ) is commonly used to test for signatures of selection in population divergence. However, there is an ongoing discussion about what F(ST) actually measures, even resulting in some alternative metrics to express neutral genetic differentiation. If there is a problem with F(ST) , this could have repercussions for its comparison with Q(ST) as well. We show that as the mutation rate of the neutral marker increases, F(ST) decreases: a higher within-population heterozygosity (He) yields a lower F(ST) value. However, the same is true for Q(ST) : a higher mutation rate for the underlying QTL also results in a lower Q(ST) estimate. The effect of mutation rate is equivalent in Q(ST) and F(ST) . Hence, the comparison between Q(ST) and F(ST) remains valid, if one uses neutral markers whose mutation rates are not too high compared to those of quantitative traits. Usage of highly variable neutral markers such as hypervariable microsatellites can lead to serious biases and the incorrect inference that divergent selection has acted on populations. Much of the discussion on F(ST) seems to stem from the misunderstanding that it measures the differentiation of populations, whereas it actually measures the fixation of alleles. In their capacity as measures of population differentiation, Hedrick's G'(ST) and Jost's D reach their maximum value of 1 when populations do not share alleles even when there remains variation within populations, which invalidates them for comparisons with Q(ST) .

Adaptation to local ultraviolet radiation conditions among neighbouring Daphnia populations

Understanding the historical processes that generated current patterns of phenotypic diversity in nature is particularly challenging in subdivided populations. Populations often exhibit heritable genetic differences that correlate with environmental variables, but the non-independence among neighbouring populations complicates statistical inference of adaptation. To understand the relative influence of adaptive and non-adaptive processes in generating phenotypes requires joint evaluation of genetic and phenotypic divergence in an integrated and statistically appropriate analysis. We investigated phenotypic divergence, population-genetic structure and potential fitness trade-offs in populations of Daphnia melanica inhabiting neighbouring subalpine ponds of widely differing transparency to ultraviolet radiation (UVR). Using a combination of experimental, population-genetic and statistical techniques, we separated the effects of shared population ancestry and environmental variables in predicting phenotypic divergence among populations. We found that native water transparency significantly predicted divergence in phenotypes among populations even after accounting for significant population structure. This result demonstrates that environmental factors such as UVR can at least partially account for phenotypic divergence. However, a lack of evidence for a hypothesized trade-off between UVR tolerance and growth rates in the absence of UVR prevents us from ruling out the possibility that non-adaptive processes are partially responsible for phenotypic differentiation in this system.

Influence of mutation rate on estimators of genetic differentiation--lessons from Arabidopsis thaliana

BACKGROUND

The analysis of molecular variation within and between populations is crucial to establish strategies for conservation as well as to detect the footprint of spatially heterogeneous selection. The traditional estimator of genetic differentiation (F(ST)) has been shown to be misleading if genetic diversity is high. Alternative estimators of F(ST) have been proposed, but their robustness to variation in mutation rate is not clearly established. We first investigated the effect of mutation and migration rate using computer simulations and examined their joint influence on Q(ST), a measure of genetic differentiation for quantitative traits. We further used experimental data in natural populations of Arabidopsis thaliana to characterize the effect of mutation rate on various estimates of population differentiation. Since natural species exhibit various degrees of self-fertilisation, we also investigated the effect of mating system on the different estimators.

RESULTS

If mutation rate is high and migration rate low, classical measures of genetic differentiation are misleading. Only Phi(ST), an estimator that takes the mutational distances between alleles into account, is independent of mutation rate, for all migration rates. However, the performance of Phi(ST) depends on the underlying mutation model and departures from this model cause its performance to degrade. We further show that Q(ST) has the same bias. We provide evidence that, in A. thaliana, microsatellite variation correlates with mutation rate. We thereby demonstrate that our results on estimators of genetic differentiation have important implications, even for species that are well established models in population genetics and molecular biology.

CONCLUSIONS

We find that alternative measures of differentiation like F'(ST) and D are not suitable for estimating effective migration rate and should not be used in studies of local adaptation. Genetic differentiation should instead be measured using an estimator that takes mutation rate into account, such as Phi(ST). Furthermore, in systems where migration between populations is low, such as A. thaliana, Q(ST) < F(ST) cannot be taken as evidence for homogenising selection as has been traditionally thought.