Evolutionary aspects of population structure for molecular and quantitative traits in the freshwater snail Radix balthica

Detecting the action of selection in natural populations can be achieved using the QST-FST comparison that relies on the estimation of FST with neutral markers, and QST using quantitative traits potentially under selection. QST higher than FST suggests the action of directional selection and thus potential local adaptation. In this article, we apply the QST-FST comparison to four populations of the hermaphroditic freshwater snail Radix balthica located in a floodplain habitat. In contrast to most studies published so far, we did not detect evidence of directional selection for local optima for any of the traits we measured: QST calculated using three different methods was never higher than FST. A strong inbreeding depression was also detected, indicating that outcrossing is probably predominant over selfing in the studied populations. Our results suggest that in this floodplain habitat, local adaptation of R. balthica populations may be hindered by genetic drift, and possibly altered by uneven gene flow linked to flood frequency.     

Disjunct distribution of highly diverged mitochondrial lineage clade and population subdivision in a marine bivalve with pelagic larval dispersal

Mitochondrial DNA sequence data for 295 individuals of the marine bivalve Macoma balthica (L.) were collected from 10 sites across the European distribution, and from Alaska. The data were used to infer population subdivision history and estimate current levels of gene flow. Inferred historical biogeography was expected to be congruent with colonization of the Atlantic Ocean from the Pacific Ocean after the opening of the Bering Strait 3.5 Ma. In addition, the last glacial maximum, about 18000 years ago, was expected to have been responsible for most of the present-day distribution of molecular variation within Europe, because the area must have been recolonized after confinement to France and the south of the British Isles during the last glacial maximum. Current gene flow was hypothesized to be high, because the larvae of M. balthica spend 2-5 weeks drifting in the water column. The geographical distribution of one highly diverged haplotype clade was found to be disjunct and was encountered exclusively in samples from the Baltic Sea and Alaska. A molecular clock calibration for marine bivalve cytochrome-c-oxidase I dates this clade as having split off from the other haplotypes 9.8-39 Ma. Multiple colonizations of the Atlantic Ocean from the Pacific by M. balthica may explain the strong differences found between Baltic Sea and other European populations of this species. The sympatric occurrence of the highly diverged mitochondrial lineages in western parts of the Baltic Sea points to secondary admixture. With the use of coalescent analysis, population divergence times for French vs. other non-Baltic European populations ('Atlantic population assemblage') were estimated at a minimum of about 110000 years ago, well before the last glacial maximum 18000 years ago. Signatures of population divergence of M. balthica that appear to have originated during the Pleistocene have thus survived the last glacial maximum. Some of the populations within the Atlantic assemblage are currently isolated, while others appear to be connected by gene flow. Apparently, populations of this species can remain highly subdivided in spite of the potential for high gene flow, implying that their population and evolutionary dynamics can be independent.     

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

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

Population structure in the barn swallow,Hirundo rustica: a comparison between neutral DNA markers and quantitative traits

Local adaptation to variable environments can generate clinal variation in morphological traits. Alternatively, similar patterns of clinal variation may be generated simply as a result of genetic drift/migration balance. Teasing apart these different processes is a continuing focus in evolutionary ecology. We compare genetic differentiation at molecular loci and quantitative traits to analyse the effect of these different processes in a morphological latitudinal cline of the barn swallow, Hirundo rustica, breeding across Europe. The results obtained show no structuring at neutral microsatellite loci, which contrasts with positive structuring at five quantitative morphometric traits. This supports the hypothesis that the observed morphometric cline in barn swallows is the result of selection acting in a spatially heterogeneous environment. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 99 , 306–314.     

Contrasting patterns of divergence in quantitative traits and neutral DNA markers: analysis of clinal variation

Clinal variation in quantitative traits is often attributed to the effects of spatially varying selection. However, identical patterns can be produced by the interplay between purely stochastic processes (i.e. drift in combination with spatially restricted gene flow). One means of distinguishing between adaptive and nonadaptive causes of geographical variation is to compare relative levels of between-population divergence in quantitative traits and neutral DNA markers. Such comparisons can be used to test whether levels of trait divergence attributable to additive genetic effects (as measured by QST) exceed null expectations based on the level of divergence at neutral marker loci (as measured by FST). The purpose of this study was to use an approach based on 'QST vs. FST' contrasts to test for evidence of diversifying selection on body size of an Indian fruit bat, Cynopterus sphinx (Chiroptera: Pteropodidae). Specifically, relative levels of between-population divergence in body size and microsatellite DNA markers were compared to assess whether the observed pattern of clinal size variation could be explained by a neutral model of isolation by distance. QST for body size was calculated using unbiased estimators of within- and between-population variance of principal component scores. The association between body size variation and geographical/environmental distance was tested using pairwise and partial matrix correspondence tests (MCTs). Independent variables (representing causal hypotheses) were constructed as between-locality distance matrices. The effects of neutral genetic divergence were assessed by including a matrix of pairwise FST as an independent variable. Partial MCTs revealed highly significant associations between phenotypic divergence (QST) and both geographical and environmental distance, even when the effects of neutral genetic divergence (FST) were partialled out. Results of the tests confirmed that migration-drift equilibrium is not a sufficient explanation for the latitudinal pattern of clinal size variation in C. sphinx. The geographical patterning of pairwise QST is most likely attributable to spatially varying selection and/or the direct influence of latitudinally ordered environmental effects.     

Contrasting geographic patterns of genetic variation for molecular markers vs. phenotypic traits in the energy grass Miscanthus sinensis

Species and hybrids of Miscanthus are a promising energy crop, but their outcrossing mating systems and perennial life cycles are serious challenges for breeding programs. One approach to accelerating the domestication of Miscanthus is to harness the tremendous genetic variation that is present within this genus using phenotypic data from extensive field trials, high‐density genotyping and sequencing technologies, and rapidly developing statistical methods of relating phenotype to genotype. The success of this approach, however, hinges on detailed knowledge about the population genetic structure of the germplasm used in the breeding program. We therefore used data for 120 single‐nucleotide polymorphism and 52 simple sequence repeat markers to depict patterns of putatively neutral population structure among 244 Miscanthus genotypes grown in a field trial near Aberystwyth (UK) and delineate a population of 145 M . sinensis genotypes that will be used for association mapping and genomic selection. Comparative multivariate analyses of molecular marker and phenotypic data for 17 traits related to phenology, morphology/biomass, and cell wall composition revealed significant geographic patterns in this population. A longitudinal cline accounted for a substantial proportion of molecular marker variation (R² = 0.60, P = 3.4 × 10^?15). In contrast, genetic variation for phenotypic traits tended to follow latitudinal and altitudinal gradients, with several traits appearing to have been affected by divergent selection (i.e., Q_ST >> F_ST). These contrasting geographic trends are unusual relative to other plants and provide opportunities for powerful studies of phenotype–genotype associations and the evolutionary history of M. sinensis.     

Inferring local adaptation from QST-FST comparisons: neutral genetic and quantitative trait variation in European populations of great snipe

We applied a phenotypic QST (PST) vs. FST approach to study spatial variation in selection among great snipe (Gallinago media) populations in two regions of northern Europe. Morphological divergence between regions was high despite low differentiation in selectively neutral genetic markers, whereas populations within regions showed very little neutral divergence and trait differentiation. QST > FST was robust against altering assumptions about the additive genetic proportions of variance components. The homogenizing effect of gene flow (or a short time available for neutral divergence) has apparently been effectively counterbalanced by differential natural selection, although one trait showed some evidence of being under uniform stabilizing selection. Neutral markers can hence be misleading for identifying evolutionary significant units, and adopting the PST-FST approach might therefore be valuable when common garden experiments is not an option. We discuss the statistical difficulties of documenting uniform selection as opposed to divergent selection, and the need for estimating measurement error. Instead of only comparing overall QST and FST values, we advocate the use of partial matrix permutation tests to analyse pairwise QST differences among populations, while statistically controlling for neutral differentiation.     

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.     

Evolution in heterogeneous environments and the potential of maintenance of genetic variation in traits of adaptive significance

The maintenance of genetic variation in traits of adaptive significance has been a major dilemma of evolutionary biology. Considering the pattern of increased genetic variation associated with environmental clines and heterogeneous environments, selection in heterogeneous environments has been proposed to facilitate the maintenance of genetic variation. Some models examining whether genetic variation can be maintained, in heterogeneous environments are reviewed. Genetic mechanisms that constrain evolution in quantitative genetic traits indicate that genetic variation can be maintained but when is not clear. Furthermore, no comprehensive models have been developed, likely due to the genetic and environmental complexity of this issue. Therefore, I have suggested two empirical approaches to provide insight for future theoretical and empirical research. Traditional path analysis has been a very powerful approach for understanding phenotypic selection. However, it requires substantial information on the biology of the study system to construct a causal model and alternatives. Exploratory path analysis is a data driven approach that uses the statistical relationships in the data to construct a set of models. For example, it can be used for understanding phenotypic selection in different environments, where there is no prior information to develop path models in the different environments. Data from Brassica rapa grown in different nutrients indicated that selection changed in the different environments. Experimental evolutionary studies will provide direct tests as to when genetic variation is maintained.     

Life history traits variation in heterogeneous environment: The case of a freshwater snail resistance to pond drying

Ecologists and population geneticists have long suspected that the diversity of living organisms was connected to the structure of their environment. In heterogeneous environments, diversifying selection combined to restricted gene flow may indeed lead to locally adapted populations. The freshwater snail, Galba truncatula, is a good model to address this question because it is present in a heterogeneous environment composed of temporary and permanent waters. In order to test the selective importance of those environments, we proposed here to measure survival of lineages from both habitats during drought episodes. To this purpose, we experimentally submitted adults and juveniles individuals from both habitats to drought. We found a difference in desiccation resistance between temporary and permanents waters only for adults. Adults from temporary habitats were found more resistant to drought. This divergence in desiccation resistance seems to explain the unexpected life history traits differences between habitats observed.     

Stable genetic polymorphism in heterogeneous environments: balance between asymmetrical dispersal and selection in the acorn barnacle

Elucidating the processes responsible for maintaining polymorphism at ecologically relevant genes is intimately related to understanding the interplay between selection imposed by habitat heterogeneity and a species' capacity for dispersal in the face of environmental constraints. In this paper, we used a model-based approach to solve equilibria of balanced polymorphism, given values of fitness and larval dispersal among different habitats in the acorn barnacle Semibalanus balanoides from the Gulf of St Lawrence. Our results showed that allele frequencies observed at both MPI* and GPI* loci represented stable equilibria, given empirical estimates of fitness values, and that considerably more larvae dispersed from one region (north) to the other (south) than vice versa. Dispersal conditions were predicted to be similar for the maintenance of polymorphism at both loci. Moreover, the values of asymmetrical dispersal required by the model to reach stable equilibria were compatible with empirical estimates of larval dispersal and oceanic circulation documented in this system. Overall, this study illustrated the usefulness of a modified and computable version of Bulmer's model (1972) in order to test hypotheses of balanced polymorphism resulting from interactions between spatial selection and asymmetrical dispersal.     

Historical introgression and the role of selective vs. neutral processes in structuring nuclear genetic variation (AFLP) in a circumpolar marine fish, the capelin (Mallotus villosus)

The capelin (Mallotus villosus) is a widespread marine fish species for which previous work has identified geographically distinct mtDNA clades, the frontiers of which are well within adult and larval dispersal capabilities. Here, we use AFLPs to test for the presence of nuclear gene flow among clades. In addition, we evaluate genetic structuring within one clade, the Northwest Atlantic (NWA). We found that each of the mtDNA clades corresponds with a unique nuclear DNA genetic cluster. Within the NWA clade, we detected individuals with small but significant amounts of genetic ancestry from other clades, likely due to historical introgression. Further support for historical introgression comes from analyses of variance in locus-specific differentiation, which support introgression between some clades and divergence without gene flow between others. Within the NWA, we identified two genetic clusters that correspond to sites in geographically adjacent areas. However, these clusters differ primarily at 'outlier' loci, and a genetic subdivision (K=2) was not supported by genetic clustering programs using neutral loci. Significant neutral F(ST) differentiation was found only between sites that otherwise differed at outlier loci. Thus, these populations may be in the initial stages of 'isolation by adaptation'. These results suggest strong between-clade reproductive isolation despite opportunities for gene flow and support the hypothesis that selection can contribute to divergence in otherwise 'open' systems.     

Genotype matrix mapping: searching for quantitative trait loci interactions in genetic variation in complex traits.

In order to reveal quantitative trait loci (QTL) interactions and the relationship between various interactions in complex traits, we have developed a new QTL mapping approach, named genotype matrix mapping (GMM), which searches for QTL interactions in genetic variation. The central approach in GMM is the following. (1) Each tested marker is given a virtual matrix, named a genotype matrix (GM), containing intersecting lines and rows equal to the total allele number for that marker in the population analyzed. (2) QTL interactions are then estimated and compared through virtual networks among the GMs. To evaluate the contribution of marker combinations to a quantitative phenotype, the GMM method divides the samples into two non-overlapping subclasses, S(0) and S(1); the former contains the samples that have a specific genotype pattern to be evaluated, and the latter contains samples that do not. Based on this division, the F-measure is calculated as an index of significance. With the GMM method, we extracted significant marker combinations consisting of one to three interacting markers. The results indicated there were multiple QTL interactions affecting the phenotype (flowering date). GMM will be a valuable approach to identify QTL interactions in genetic variation of a complex trait within a variety of organisms.     

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.     

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.     

Contrasting patterns of floral and molecular variation across a cline in Mimulus aurantiacus

Abstract Steep clines in ecologically important traits may be caused by divergent natural selection. However, processes that do not necessarily invoke ongoing selection, such as secondary contact or restricted gene flow, can also cause patterns of phenotypic differentiation over short spatial scales. Distinguishing among all possible scenarios is difficult, but an attainable goal is to establish whether scenarios that imply selection need to be invoked. We compared the extent of morphological and genetic differentiation between geographically structured red and yellow floral races of Mimulus aurantiacus (bush monkeyflower; Phrymacea). Flower color was assessed in a common garden as well as in the field to determine whether variation was genetic and to quantify the extent of geographical differentiation. Population genetic differentiation at marker loci was measured for both chloroplast and nuclear genomes, and the degree of population structure within and among the floral races was evaluated. Flower color shows both a strong genetic basis and a sharp geographic transition, with pure red-flowered populations in western San Diego County and pure yellow-flowered populations to the east. In the zone of contact, both pure and intermediate phenotypes occur. Patterns of genetic differentiation at marker loci are far less pronounced, as little of the variation is partitioned according to the differences in flower color. Phenotypic differentiation (Q_ST) between populations with different flower colors is much greater than neutral genetic differentiation (F_ST). When comparisons are made between populations of the same flower color, the opposite trend is evident. Limited neutral genetic structure between the floral races, combined with sharp differentiation in flower color, is consistent with the hypothesis that current or recent natural selection maintains the cline in flower color.     

On the neutrality of molecular genetic markers: pedigree analysis of genetic variation in fragmented populations

Many studies employ molecular markers to infer ecological and evolutionary processes, assuming that variation found at genetic loci offers a reliable representation of stochastic events in natural populations. Increasingly, evidence emerges that molecular markers might not always be selectively neutral. However, only a few studies have analysed how deviations from neutrality could affect estimates of genetic variation, using populations with known genealogy. We monitored changes in allozyme variation over eight generations in captive metapopulations of the butterfly Bicyclus anynana. Population demography was recorded by individually marking 35 000 butterflies and constructing pedigrees. We designed a computer program that simulated the inheritance of founder allozyme alleles in butterfly pedigrees. We thus tested whether the observed transmission of allozyme alleles could be explained by random genetic drift alone, or whether there was evidence for positive or negative selection. This analysis showed that in the smallest metapopulations the loss of allozyme variation exceeded the neutral rate. Possibly, linkage disequilibria between deleterious mutations and marker alleles resulted in background selection and a faster erosion of allozyme variation. In larger metapopulations, one locus (MDH) showed a significant heterozygote excess and smaller than expected loss in heterozygosity, observations consistent with (associative) overdominance. This study demonstrates that the neutrality of molecular markers cannot always be assumed, particularly in small populations with a high mutation load.