Molecular and quantitative trait variation within and among populations of the intertidal copepod Tigriopus californicus

Abstract While molecular and quantitative trait variation may be theoretically correlated, empirical studies using both approaches frequently reveal discordant patterns, and these discrepancies can contribute to our understanding of evolutionary processes. Here, we assessed genetic variation in six populations of the copepod Tigriopus californicus. Molecular variation was estimated using five polymorphic microsatellite loci, and quantitative variation was measured using 22-life history and morphometric characters. Within populations, no correlation was found between the levels of molecular variation (heterozygosity) and quantitative variation (heritability). Between populations, quantitative subdivision ( Q _ST) was correlated with molecular subdivision when measured as F _ST but not when measured as R _ST. Unlike most taxa studied to date, the overall level of molecular subdivision exceeded the level of quantitative subdivision ( F _ST= 0.80, R_ST = 0.89, Q _ST = 0.30). Factors that could contribute to this pattern include stabilizing or fluctuating selection on quantitative traits or accelerated rates of molecular evolution.     

Population divergence in the amphicarpic species Amphicarpaea edgeworthii Benth. (Fabaceae): microsatellite markers and leaf morphology

Comparative analyses of the genetic differentiation in microsatellite markers ( F _ST) and leaf morphology characters ( Q _ST) of Amphicarpaea edgeworthii Benth. were conducted to gain insight into the roles of random processes and natural selection in the population divergence. Simple sequence repeat analyses on 498 individuals of 19 natural populations demonstrate that a significant genetic differentiation occurs among populations (mean F _ST = 0.578), and A. edgeworthii is a highly self-fertilized species (mean selfing rate s  = 0.989). The distribution pattern of genetic diversity in this species shows that central populations possess high genetic diversity (e.g. population WL with H _E = 0.673 and population JG with H _E = 0.663), whereas peripheral ones have a low H _E as in population JD (0.011). The morphological divergence of leaf shape was estimated by the elliptical Fourier analysis on the data from 11 natural and four common garden populations. Leaf morphology analyses indicate the morphological divergence does not show strong correlation with the genetic differentiation ( R  = 0.260, P  = 0.069). By comparing the 95% confidence interval of Q _ST with that of F _ST, Q _ST values for five out of 12 quantitative traits are significantly higher than the average F _ST value over eight microsatellite loci. The comparison of F _ST and Q _ST suggests that two kinds of traits can be driven by different evolutionary forces, and the population divergence in leaf morphology is shaped by local selections.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 96 , 505–516.     

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.     

Q ST MEETS THE G MATRIX: THE DIMENSIONALITY OF ADAPTIVE DIVERGENCE IN MULTIPLE CORRELATED QUANTITATIVE TRAITS

The Q _ST– F _ST comparison has become an increasingly common method for inferring adaptive quantitative trait divergence among populations. For cases in which there is divergence in multiple traits, most studies have applied the method by performing multiple univariate Q _ST– F _ST comparisons. However, because traits are often genetically correlated, such univariate analyses are likely to paint a simplified picture of adaptive divergence. Here we show how the multivariate analogue of Q _ST, F_STq, which accounts for genetic correlations among traits, can be used to supply a more detailed picture of multitrait divergence. We apply the method to naturally occurring genetic variation for a suite of sexually selected display traits in Drosophila serrata . The analyses suggest the operation of divergent multivariate selection that has influenced multiple independent axes of genetic variance in a sex-specific manner. Finally, we show how a comparison of the components of F_STq, the average within and among population genetic variance–covariance matrices, G_W and G_B, can be used as an additional test of the null expectation of neutral divergence, and allows for an investigation of whether natural populations have diverged along major or minor axes of genetic variance.     

Genetic Variation and Clonal Diversity of the Two Divergent Types of Clonal Populations of Leymus chinensis Tzvel on the Song Liao Steppe in the West of Northeastern China

The genetic variation and clonal diversity of two divergent types (grey-green and yellow-green) of clonal populations of Leymus chinensis Tzvel at 14 loci were compared. Total gene diversity (HT) and the coefficient of genetic differentiation (GST) were all higher for the yellow-green type (HT = 0.270; GST =0.186) than for the grey-green type (HT = 0.250; GST = 0.157) of L. chinensis. Rare alleles usually occurred as heterozygotes rather than homozygotes and significant deviations from Hardy-Weinberg equilibrium were found only at a few loci. This indicated that these two types of populations were mainly out-crossing. Clonal diversity, evenness of clones, and mean clone size were not significantly different between the two types. We found that differences between the clone size and genetic variation of the yellow-green type of populations occurred with different climate and habitat population groups. However, for the grey-green type of populations, these genetic variations decreased under conditions of different climate and habitat population groups.     

Spatial and temporal scales of adaptive divergence in marine fishes and the implications for conservation

Knowledge of geographic and temporal scales of adaptive genetic variation is crucial to species conservation, yet understanding of these phenomena, particularly in marine systems, is scant. Until recently, the belief has been that because most marine species have highly dispersive or mobile life stages, local adaptation could occur only on broad geographic scales. This view is supported by comparatively low levels of genetic variation among populations as detected by neutral markers. Similarly, the time scale of adaptive divergence has also been assumed to be very long, requiring thousands of generations. Recent studies of a variety of species have challenged these beliefs. First, there is strong evidence of geographically structured local adaptation in physiological and morphological traits. Second, the proportion of quantitative trait variation at the among-population level ( Q _ST) is much higher than it is for neutral markers ( F _ST) and these two metrics of genetic variation are poorly correlated. Third, evidence that selection is a potent evolutionary force capable of sustaining adaptive divergence on contemporary time scales is summarized. The differing spatial and temporal scales of adaptive v. neutral genetic divergence call for a new paradigm in thinking about the relationship between phenogeography (the geography of phenotypic variation) and phylogeography (the geography of lineages) in marine species. The idea that contemporary selective processes can cause fine-scale spatial and temporal divergence underscores the need for a new emphasis on Darwinian fishery science.     

Subtle population structuring within a highly vagile marine invertebrate, the veined squid Loligo forbesi, demonstrated with microsatellite DNA markers

Microsatellite DNA markers were applied for the first time in a population genetic study of a cephalopod and compared with previous estimates of genetic differentiation obtained using allozyme and mitochondrial DNA (mtDNA) markers. Levels of genetic variation detected with microsatellites were much higher than found with previous markers (mean number of alleles per locus=10.6, mean expected heterozygosity ( H _E)=0.79; allozyme H _E=0.08; mtDNA restriction fragment length polymorphism (RFLP) H _E=0.16). In agreement with previous studies, microsatellites demonstrated genetic uniformity across the population occupying the European shelf seas of the North East Atlantic, and extreme genetic differentiation of the Azores population ( R _ST/ F _ST=0.252/0.245; allozyme F _ST=0.536; mtDNA F _ST=0.789). In contrast to other markers, microsatellites detected more subtle, and significant, levels of differentiation between the populations of the North East Atlantic offshore banks (Rockall and Faroes) and the shelf population ( R _ST=0.048 and 0.057). Breakdown of extensive gene flow among these populations is indicated, with hydrographic (water depth) and hydrodynamic (isolating current regimes) factors suggested as possible barriers to migration. The demonstration of genetic subdivision in an abundant, highly mobile marine invertebrate has implications for the interpretation of dispersal and population dynamics, and consequent management, of such a commercially exploited species. Relative levels of differentiation indicated by the three different marker systems, and the use of measures of differentiation (assuming different mutation models), are discussed.     

Geographic patterns of genetic differentiation and plumage colour variation are different in the pied flycatcher (Ficedula hypoleuca)

The pied flycatcher is one of the most phenotypically variable bird species in Europe. The geographic variation in phenotypes has often been attributed to spatial variation in selection regimes that is associated with the presence or absence of the congeneric collared flycatcher. Spatial variation in phenotypes could however also be generated by spatially restricted gene flow and genetic drift. We examined the genetic population structure of pied flycatchers across the breeding range and applied the phenotypic Q _ST ( P _ST)– F _ST approach to detect indirect signals of divergent selection on dorsal plumage colouration in pied flycatcher males. Allelic frequencies at neutral markers were found to significantly differ among populations breeding in central and southern Europe whereas northerly breeding pied flycatchers were found to be one apparently panmictic group of individuals. Pairwise differences between phenotypic ( P _ST) and neutral genetic distances ( F _ST) were positively correlated after removing the most differentiated Spanish and Swiss populations from the analysis, suggesting that genetic drift may have contributed to the observed phenotypic differentiation in some parts of the pied flycatcher breeding range. Differentiation in dorsal plumage colouration however greatly exceeded that observed at neutral genetic markers, which indicates that the observed pattern of phenotypic differentiation is unlikely to be solely maintained by restricted gene flow and genetic drift.     

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

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

THE QUANTITATIVE AND MOLECULAR GENETIC ARCHITECTURE OF A SUBDIVIDED SPECIES

In an effort to elucidate the evolutionary mechanisms that determine the genetic architecture of a species, we have analyzed 17 populations of the microcrustacean Daphnia pulex for levels of genetic variation at the level of life-history characters and molecular markers in the nuclear and mitochondrial genomes. This species is highly subdivided, with approximately 30% of the variation for nuclear molecular markers and 50% of the variation for mitochondrial markers being distributed among populations. The average level of genetic subdivision for quantitative traits is essentially the same as that for nuclear markers, which superficially suggests that the life-history characters are diverging at the neutral rate. However, the existence of a strong correlation between the levels of population subdivision and broadsense heritabilities of individual traits argues against this interpretation, suggesting instead that the among-population divergence of some quantitative traits (most notably body size) is being driven by local adaptation to different environments. The fact that the mean phenotypes of the individual populations are also strongly correlated with local levels of homozygosity indicates that variation in local inbreeding plays a role in population differentiation. Rather than being a passive consequence of local founder effects, levels of homozygosity may be selected for directly for their effects on the phenotype (adaptive inbreeding depression). There is no relationship between the levels of variation within populations for molecular markers and quantitative characters, and this is explained by the fact that the average standing genetic variation for life-history characters in this species is equivalent to only 33 generations of variation generated by mutation.     

Comparison of genetic diversity at microsatellite loci and quantitative traits in hatchery populations of Japanese flounder Paralichthys olivaceus

Relationships of genetic diversity at microsatellite loci and quantitative traits were examined in hatchery-produced populations of Japanese flounder using a relatively straightforward experiment. Five hatchery populations produced by wild-caught and domesticated broodstocks were used to examine the effects of different levels (one to three generations) of domestication on the genetic characteristics of hatchery populations. Allelic richness at seven microsatellite loci in all hatchery populations was lower than that in a wild population. Genetic variation measured by allelic richness and heterozygosity tended to decrease with an increase in generations of domestication. In addition, the degree of genetic differentiation from a wild population increased with an increase in generations of domestication. Significant differences in three morphometric traits (dorsal and anal fin ray counts and vertebral counts) and three physiological traits (high temperature, salinity and formalin tolerance) were observed among the hatchery populations. The degree of phenotypic difference among populations was larger in morphometric traits than in physiological traits. The divergence pattern of some quantitative traits was similar to that observed at microsatellite loci, suggesting that domestication causes the decrease of genetic variation and the increase of genetic differentiation for some quantitative traits concomitantly with those for microsatellite loci. Significant positive correlation was observed between F _ST and the degree of phenotypic difference in the three morphometric traits and formalin tolerance, indicating that genetic variation at microsatellite loci predicts the degree of phenotypic divergence in some quantitative traits.     

Allozyme Diversity in Populations of Cymbidium goeringii (Orchidaceae)

Abstract: Using 14 allozyme loci, we investigated levels of genetic diversity within populations, and degree of genetic divergence among 24 populations of Cymbidium goeringii (Orchidaceae) in Korea and Japan. Cymbidium goeringii maintains high levels of genetic diversity both at population (mean expected heterozygosity, H _e = 0.238) and species levels (0.260). Means of H _e found in 24 populations were not significantly different from each other. About 90 % of the total variation in the species is common to all populations (mean G _ST = 0.108). No unique allele was found in any population. The indirect estimate of gene flow based on the mean G _ST was high ( Nm = 2.06). Nei's genetic identities for pairs of populations had high values (mean = 0.974 [SD = 0.013]). The Mantel-Z test showed a significant correlation between genetic distance and geographic distance. However, the mean G _ST value between 17 populations in Korea and seven Japanese populations was relatively low (0.029), even though the land connection between the southern Korean peninsula and southern Japanese archipelagos has not existed since the middle Pleistocene. Large numbers of small seeds of C. goeringii might travel long distances by wind from populations to populations both in Korea and Japan, increasing genetic diversity within populations and maintaining low genetic differentiation among populations.     

History, chance and adaptation during biological invasion: separating stochastic phenotypic evolution from response to selection

Introduced species often exhibit changes in genetic variation, population structure, selection regime and phenotypic traits as they colonize and expand into new ranges. For these reasons, species invasions are increasingly recognized as promising systems for studying adaptive evolution over contemporary time scales. However, changes in phenotypic traits during invasion occur under non-equilibrium demographic conditions and may reflect the influences of prior evolutionary history and chance events, as well as selection. We briefly review the evidence for phenotypic evolution and the role of selection during invasion. While there is ample evidence for evolutionary change, it is less clear if selection is the primary mechanism. We then discuss the likelihood that stochastic events shift phenotypic distributions during invasion, and argue that hypotheses of adaptation should be tested against appropriate null models. We suggest two experimental frameworks for separating stochastic evolution from adaptation: statistically accounting for phenotypic variation among putative invasion sources identified by using phylogenetic or assignment methods and by comparing estimates of differentiation within and among ranges for both traits and neutral markers ( Q _ST vs. F _ST). Designs that incorporate a null expectation can reveal the role of history and chance in the evolutionary process, and provide greater insights into evolution during species invasions.     

Ancestry and divergence of subtropical montane forest isolates: molecular biogeography of the genus Abies (Pinaceae) in southern México and Guatemala

The genus Abies has a complex history in southern México and Guatemala. In this region, four closely related species, Abies flinckii , A. guatemalensis , A. hickelii , and A. religiosa , are distributed in fragmented and isolated montane populations. Range-wide genetic variation was investigated across species using cytoplasmic DNA markers with contrasted inheritance. Variation at two maternally inherited mitochondrial DNA markers was low. All species shared two of the nine mitotypes detected, while the remaining seven mitochondrial DNA types were restricted to a few isolated stands. Mitochondrial genetic differentiation across taxa was high ( G _ST = 0.933), it was not related to the taxonomic identity ( amova ; P  > 0.05) of the populations, and it was not phylogeographically structured ( G _ST ≈  N _ST). In contrast, variation at three paternally inherited chloroplast DNA microsatellites was high. Chloroplast genetic differentiation was lower ( G _ST = 0.402; R _ST = 0.547) than for mitochondrial DNA, but it was significantly related to taxonomy ( amova ; P  < 0.001), and exhibited a significant phylogeographical structure ( G _ST <  R _ST). Different analyses of population structure indicated that A. flinckii was the most divergent taxon, while the remaining three species formed a relatively homogeneous group. However, a small number of the populations of these three taxa, all located at the limits of their respective ranges or in the Transverse Volcanic Belt, diverged from this main cluster. These trends suggest that the Mesoamerican Abies share a recent common ancestor and that their divergence and speciation is mainly driven by genetic drift and isolation during the warm interglacial periods.     

Microsatellite analysis of genetic population structure in an endangered salmonid: the coastal cutthroat trout (Oncorhynchus clarki clarki)

The genetic population structure of coastal cutthroat trout ( Oncorhynchus clarki clarki ) in Washington state was investigated by analysis of variation in allele frequencies at six highly polymorphic microsatellite loci for 13 anadromous populations, along with one outgroup population from the Yellowstone subspecies ( O. clarki bouvieri) (mean heterozygosity = 67%; average number of alleles per locus = 24). Tests for genetic differentiation revealed highly significant differences in genotypic frequencies for pairwise comparisons between all populations within geographical regions and overall population subdivision was substantial ( F _ST = 0.121, R _ST = 0.093), with 44.6% and 55.4% of the among-population diversity being attributable to differences between streams ( F _SR = 0.054) and between regions ( F _RT = 0.067), respectively. Analysis of genetic distances and geographical distances did not support a simple model of isolation by distance for these populations. With one exception, neighbour-joining dendrograms from the Cavalli-Sforza and Edwards' chord distances and maximum likelihood algorithms clustered populations by physiogeographic region, although overall bootstrap support was relatively low (53%). Our results suggest that coastal cutthroat trout populations are ultimately structured genetically at the level of individual streams. It appears that the dynamic balance between gene flow and genetic drift in the subspecies favours a high degree of genetic differentiation and population subdivision with the simultaneous maintenance of high heterozygosity levels within local populations. Results are discussed in terms of coastal cutthroat trout ecology along with implications for the designation of evolutionarily significant units pursuant to the US Endangered Species Act of 1973 and analogous conservation units.     

Variation for neutral markers is correlated with variation for quantitative traits in the plant pathogenic fungus Mycosphaerella graminicola

We compared genetic variation and population differentiation at RFLP marker loci with seven quantitative characters including fungicide resistance, temperature sensitivity, pycnidial size, pycnidial density, colony size, percentage of leaves covered by pycnidia (PLACP) and percentage of leaves covered by lesions (PLACL) in Mycosphaerella graminicola populations sampled from four regions. Wide variation in population differentiation was found across the quantitative traits assayed. Fungicide resistance, temperature sensitivity, and PLACP displayed a significantly higher Q(ST) than G(ST), consistent with selection for local adaptation, while pycnidial size, pycnidial density and colony size displayed a lower or significantly lower Q(ST) than G(ST), consistent with constraining selection. There was not a statistical difference between Q(ST) and G(ST) in PLACL. We also found a positive and significant correlation between genetic variation in molecular marker loci and quantitative traits at the multitrait scale, suggesting that estimates of overall genetic variation for quantitative traits in M. graminicola could be derived from analysis of the molecular genetic markers.     

Population genetics of a cyclostome species pair, river lamprey (Lampetra fluviatilis L.) and brook lamprey (Lampetra planeri Bloch)

Horizontal agarose gel electrophoresis of 24 allozyme loci in four species of Central European lampreys (321 Lampetra planeri , 83 L. fluviatilis , 11 Eudontomyzon mariae and nine Petromyzon marinus ) was used to study the 'paired species' L. fluviatilis and L. planeri . The genetic differentiation of the anadromous river lamprey ( L. fluviatilis ) from the stationary brook lamprey ( L. planeri ) was within the range of ingroup differentiation of the latter, but L. fluviatilis exhibited much greater population cohesion over a more extended geographic range: G _ST = 0.0537 versus G _ST = 0.3398, N _em = 4.402 versus N _em = 0.4856, mean genetic among-stock distances D = 0.0047 versus D = 0.0257. L. planeri populations coexisting geographically with L. fluviatilis in the Rhine and Elbe river systems were genetically more cohesive than L. planeri stocks from the Danubian basin where L. fluviatilis is absent. Danubian L. planeri populations exhibit a lower degree of heterozygosity than brook lampreys from the Rhine river system, but comprise deeper genetic lineages ( G _ST = 0.4629 versus G _ST = 0.2434), despite being sampled from a much more restricted area. Isolation-by-distance is observed for L. planeri from the Danubian but not from the Atlantic drainage basins. Transspecific gene flow between L. planeri from Atlantic drainage basins and the long-distant migrating L. fluviatilis is inferred, raising doubt on the validity of two separate biospecies. E. mariae and P. marinus are clearly differentiated from Lampetra spp. at several allozyme loci.     

Lack of concordance between genetic diversity estimates at the molecular and quantitative-trait levels

In many applications of population genetics, particularly in the field of conservation biology, estimates of molecular diversity are used as surrogate indicators of less easily acquired measures of genetic variation for quantitative traits. The general validity of this approach to inferring levels of quantitative genetic variation within populations is called into question by the demonstration that estimates of molecular and quantitative-genetic variation are essentially uncorrelated in natural populations of Daphnia, one of the few organisms for which multiple estimates of both quantities are available. On the other hand, molecular measures of population subdivision seem to give conservatively low estimates of the degree of genetic subdivision at the level of quantitative traits. This suggests that although molecular markers provide little information on the level of genetic variation for quantitative traits within populations, they may be valid indicators of population subdivision for such characters.     

Scanning the genome for gene single nucleotide polymorphisms involved in adaptive population differentiation in white spruce

Conifers are characterized by a large genome size and a rapid decay of linkage disequilibrium, most often within gene limits. Genome scans based on noncoding markers are less likely to detect molecular adaptation linked to genes in these species. In this study, we assessed the effectiveness of a genome-wide single nucleotide polymorphism (SNP) scan focused on expressed genes in detecting local adaptation in a conifer species. Samples were collected from six natural populations of white spruce ( Picea glauca ) moderately differentiated for several quantitative characters. A total of 534 SNPs representing 345 expressed genes were analysed. Genes potentially under natural selection were identified by estimating the differentiation in SNP frequencies among populations ( F _ST) and identifying outliers, and by estimating local differentiation using a Bayesian approach. Both average expected heterozygosity and population differentiation estimates ( H _E = 0.270 and F _ST = 0.006) were comparable to those obtained with other genetic markers. Of all genes, 5.5% were identified as outliers with F _ST at the 95% confidence level, while 14% were identified as candidates for local adaptation with the Bayesian method. There was some overlap between the two gene sets. More than half of the candidate genes for local adaptation were specific to the warmest population, about 20% to the most arid population, and 15% to the coldest and most humid higher altitude population. These adaptive trends were consistent with the genes' putative functions and the divergence in quantitative traits noted among the populations. The results suggest that an approach separating the locus and population effects is useful to identify genes potentially under selection. These candidates are worth exploring in more details at the physiological and ecological levels.     

Does natural selection promote population divergence? A comparative analysis of population structure using amplified fragment length polymorphism markers and quantitative traits

Divergent natural selection is considered an important force in plant evolution leading to phenotypic differentiation between populations exploiting different environments. Extending an earlier greenhouse study of population differentiation in the selfing annual plant Senecio vulgaris, we estimated the degree of population divergence in several quantitative traits related to growth and life history and compared these estimates with those based on presumably neutral molecular markers (amplified fragment length polymorphisms; AFLPs). This approach allowed us to disentangle the effects of divergent selection from that of other evolutionary forces (e.g. genetic drift). Five populations were examined from each of two habitat types (ruderal and agricultural habitats). We found a high proportion of total genetic variance to be among populations, both for AFLP markers (phiST = 0.49) and for quantitative traits (range of QST: 0.26-0.77). There was a strong correlation between molecular and quantitative genetic differentiation between pairs of populations (Mantel's r = 0.59). However, estimates of population differentiation in several quantitative traits exceeded the neutral expectation (estimated from AFLP data), suggesting that divergent selection contributed to phenotypic differentiation, especially between populations from ruderal and agricultural habitats. Estimates of within-population variation in AFLP markers and quantitative genetic were poorly correlated, indicating that molecular marker data may be of limited value to predict the evolutionary potential of populations of S. vulgaris.