Comparative Assessment of Genetic Variability in the Populations of Endemic and Endangered Yellow Catfish, Horabagrus brachysoma (Teleostei: Horabagridae), Based on Allozyme, RAPD, and Microsatellite Markers

The comparative assessment of genetic diversity using allozymes, random amplified polymorphic DNA (RAPD), and microsatellite markers was conducted in endemic and endangered yellow catfish (Horabagrus brachysoma) sampled from three locations in Western Ghats river systems of India. Among the three markers, microsatellites show more polymorphism, having 100% polymorphic loci, whereas allozymes show the least (56%). In RAPD, 60.5% of fragments were polymorphic. Observed heterozygosity and F(ST) values were very high in microsatellites, compared with the other markers. Microsatellite and RAPD markers reported a higher degree of genetic differentiation than allozymes among the populations depicted by pairwise F(ST)/G(ST), AMOVA, Nei's genetic distance, and UPGMA dendrogram. The three classes of markers demonstrated striking genetic differentiation between pairs of H. brachysoma populations. The data emphasize the need for fishery management, conservation, and rehabilitation of this species.     

Waste not,want not: Microsatellites remain an economical and informative technology for conservation genetics

Comparisons of microsatellites and single‐nucleotide polymorphisms (SNPs) have found that SNPs outperform microsatellites in population genetic analyses, questioning the continued utility of microsatellites in population and landscape genetics. Yet, highly polymorphic markers may be of value in species that have reduced genetic variation. This study repeated previous analyses that used microsatellites with SNPs developed from ddRAD sequencing in the black‐capped vireo source‐sink system. SNPs provided greater resolution of genetic diversity, population differentiation, and migrant detection but could not reconstruct parentage relationships due to insufficient heterozygosities. The biological inferences made by both sets of markers were similar: asymmetrical gene flow from source sites to the remaining sink sites. With the landscape genetic analyses, we found different results between the two molecular markers, but associations of the top environmental features (riparian, open habitat, agriculture, and human development) with dispersal estimates were shared between marker types. Despite the higher precision of SNPs, we find that microsatellites effectively uncover population processes and patterns and are superior for parentage analyses in this species with reduced genetic diversity. This study illustrates the continued applicability and relevance of microsatellites in population genetic research.     

If FST does not measure neutral genetic differentiation,then comparing it with QST 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) .     

Postglacial recolonization patterns and genetic relationships among whitefish (Coregonus sp.) populations in Denmark, inferred from mitochondrial DNA and microsatellite markers

The genetic relationships among morphologically and geographically divergent populations of whitefish (genus: Coregonus ) from Denmark and the Baltic Sea region were studied by analysis of microsatellites and polymerase chain reaction–restriction fragment length polymorphism (PCR–RFLP) analysis of mitochondrial DNA (mtDNA) segments. The endangered North Sea houting (classified as C. oxyrhynchus ) differs morphologically and physiologically from other Danish whitefish ( C. lavaretus ). However, limited divergence of North Sea houting was observed both at the level of mtDNA and microsatellites. The implications of these results for the conservation status of North Sea houting are discussed in the light of current definitions of evolutionary significant units. Both mtDNA and microsatellite data indicated that postglacial recolonization by C. lavaretus in Denmark was less likely to have taken place from the Baltic Sea. Instead, the data suggested a recent common origin of all Danish whitefish populations, including North Sea houting, probably by recolonization via the postglacial Elbe River system. Estimates of genetic differentiation among populations based on mtDNA and microsatellites were qualitatively different. In addition, for both classes of markers analyses of genetic differentiation yielded different results, depending on whether molecular distances between alleles or haplotypes were included.     

A comparison of biallelic markers and microsatellites for the estimation of population and conservation genetic parameters in Atlantic salmon (Salmo salar)

Biallelic markers such as single nucleotide polymorphisms (SNPs) and insertion/deletion polymorphisms have become increasingly popular markers for various population genetics applications. However, the effort required to develop biallelic markers in nonmodel organisms is still substantial. In this study, we compared the estimation of various population genetic parameters (genetic divergence and structuring, isolation-by-distance, genetic diversity) using a limited number of biallelic markers (in total 7 loci) to those estimated with 14 microsatellite loci in 21 Atlantic salmon (Salmo salar) populations from northern Europe. Pairwise FST values were significantly correlated between biallelic loci and microsatellite datasets, as was overall heterozygosity when both anadromous and nonanadromous populations were analyzed together. However, when the anadromous and nonanadromous samples were analyzed separately, only genetic divergence correlations remained significant. Biallelic markers alone were not sufficient for reliable neighbor-joining clustering of populations but gave highly similar isolation-by-distance signals when compared with microsatellites. Finally, although several population prioritization measures for conservation exhibited significant correlation between different marker types, the specific populations highlighted as being most valuable for conservation purposes varied depending on the marker type and conservation criteria applied. This study demonstrates that a relatively small set of biallelic markers can be sufficient for obtaining concordant results in most of the analyses compared with microsatellites, although estimates of genetic distance are generally more concordant than estimates of genetic diversity. This suggests that a relatively small number of biallelic markers can provide useful information for various population genetic applications. However, we emphasize that the use of much higher number of loci is preferable, especially when the genetic differences between populations are subtle or individual multilocus genotype-based analyses are to be performed.     

Species distinction in Irish populations of Quercus petraea and Q. robur: morphological versus molecular analyses

BACKGROUND AND AIMS: Populations of oak (Quercus petraea and Q. robur) were investigated using morphological and molecular (AFLP) analyses to assess species distinction. The study aimed to describe species distinction in Irish oak populations and to situate this in a European context. METHODS: Populations were sampled from across the range of the island of Ireland. Leaf morphological characters were analysed through clustering and ordination methods. Putative neutral molecular markers (AFLPs) were used to analyse the molecular variation. Cluster and ordination analyses were also performed on the AFLP markers in addition to calculations of genetic diversity and F-statisitcs. KEY RESULTS: A notable divergence was uncovered between the morphological and molecular analyses. The morphological analysis clearly differentiated individuals into their respective species, whereas the molecular analysis did not. Twenty species-specific AFLP markers were observed from 123 plants in 24 populations but none of these was species-diagnostic. Principal Coordinate Analysis of the AFLP data revealed a clustering, across the first two axes, of individuals according to population rather than according to species. High F(ST) values calculated from AFLP markers also indicated population differentiation (F(ST) = 0.271). Species differentiation accounted for only 13 % of the variation in diversity compared with population differentiation, which accounted for 27 %. CONCLUSIONS: The results show that neutral molecular variation is partitioned more strongly between populations than between species. Although this could indicate that the populations of Q. petraea and Q. robur studied may not be distinct species at a molecular level, it is proposed that the difficulty in distinguishing the species in Irish oak populations using AFLP markers is due to population differentiation masking species differences. This could result from non-random mating in small, fragmented woodland populations. Hybridization and introgression between the species could also have a significant role.     

Genetic diversity and differentiation in Eryngium alpinum L. (Apiaceae): comparison of AFLP and microsatellite markers

Genetic diversity and structure of 12 populations of Eryngium alpinum L. were investigated using 63 dominant amplified fragment length polymorphism (AFLP) and seven codominant microsatellite (48 alleles) markers. Within-population diversity estimates obtained with both markers were not correlated, but the microsatellite-based fixation index Fis was correlated with both AFLP diversity indices (number of polymorphic bands and Nei's expected heterozygosity). Only AFLP diversity indices increased with the size of populations, although they did not significantly differ among them (Kruskall-Wallis test). The discrepancy between AFLPs and microsatellites may be explained by a better coverage of the genome with numerous AFLPs, the higher mutation rates of microsatellites or the absence of significant difference among within-population diversity estimates. Genetic differentiation was higher with AFLPs (theta=0.40) than with microsatellites (theta=0.23), probably due to the higher polymorphism of microsatellites. Thus, we considered global qualitative patterns rather than absolute estimates to compare the performance of both types of markers. On a large geographic scale, the Mantel test and multivariate analysis showed that genetic patterns were more congruent with the spatial arrangement of populations when inferred from microsatellites than from AFLPs, suggesting higher homoplasy of AFLP markers. On a small spatial scale, AFLPs managed to discriminate individuals from neighboring populations whereas microsatellites did not (multivariate analysis), and the percentage of individuals correctly assigned to their population of origin was higher with AFLPs than with microsatellites. However, dominant AFLPs cannot be used to study heterozygosity-related topics. Thus, distinct molecular markers should be used depending on the biological question and the geographical scale investigated.     

Allozyme and microsatellite loci provide discordant estimates of population differentiation in the endangered dusky grouper (Epinephelus marginatus) within the Mediterranean Sea

The dusky grouper, Epinephelus marginatus, inhabits coastal reefs in the Mediterranean Sea and Atlantic Ocean. A decline in the abundance of this long-lived protogynous hermaphrodite has led to its listing as an endangered species in the Mediterranean, and heightened management concerns regarding its genetic variability and population substructure. To address these concerns, we analysed genetic variation at seven microsatellite and 28 allozyme loci in dusky groupers sampled from seven areas (for microsatellites) and three areas (for allozymes) in the west-central Mediterranean. Levels of genetic variability were higher for microsatellites than for allozymes (mean H(E) = 0.78 and 0.07, respectively), but similar to those observed in other marine fishes with comparable markers. Both microsatellites and allozymes revealed significant genetic differentiation among all areas analysed with each class of marker, but the magnitude of differentiation revealed by allozymes over three locales (F(ST) = 0.214) was greater than that detected with microsatellites over seven areas, or over the three areas shared with the allozyme analysis (F(ST) = 0.018 and approximately 0, respectively). A large proportion of the allozyme differentiation was due to a single locus (ADA*) possibly influenced by selection, but allozyme differentiation over the three areas was still highly significant (F(ST) = 0.06, P < 0.0001), and the 95% confidence intervals for allozyme and microsatellite F(ST) did not overlap when this locus was excluded. There was no evidence of isolation by distance with either class of markers. Our results lead us to conclude that dusky groupers are not panmictic in the Mediterranean Sea and suggest that they should be managed on a local basis. However, more work is needed to elucidate genetic relationships among populations.     

Variation at two flowering time genes within and among populations of Arabidopsis thaliana: comparison with markers and traits

Flowering Locus C (FLC) and Frigida are two interacting genes controlling flowering time variation in Arabidopsis thaliana. Variation at these genes was surveyed in 12 A. thaliana populations sampled in France. These populations were also screened for variation at molecular markers [12 microsatellites and 19 cleaved amplified polymorphic sequence (CAPS) markers] and at seven quantitative traits measured with and without vernalization. Seven populations were highly polymorphic at markers (H(S) = 0.57 at microsatellites, 0.24 at CAPS) and showed heritable variation for bolting time and some other traits. Five populations were genetically fixed or nearly fixed. Q(ST) for bolting time without vernalization was significantly higher than F(ST), suggesting local divergent selection. One of the two haplotype groups at FLC (FLC(A)) was very predominant (frequency of 99%). The first exon of Frigida showed elevated nonsynonymous variation, and nine loss-of-function mutations were found throughout the gene. The association between loss-of-function and earlier bolting was confirmed. Overall, 18 Frigida haplotypes were detected. The pattern of variation at Frigida was largely similar to that found at markers and traits, with the same populations being fixed or highly diverse. Metapopulation dynamics is thus probably the main factor shaping genetic variation in A. thaliana. However, F(ST) for functional (FRI) vs. nonfunctional (FRI(Delta)) haplotypes was significantly higher than F(ST) at markers. This suggested that loss-of-function at Frigida is under local selection for flowering time.     

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.     

Morphological and molecular diversity among Italian populations of Quercus petraea (Fagaceae)

Quercus petraea (sessile oak) has a scattered distribution in southern and central Italy. The objective of this work was to evaluate the level and distribution of diversity in five Italian populations of Q. petraea by using morphological markers and hypervariable molecular markers such as microsatellites. Forty-eight morphological traits and six nuclear and three plastid loci were scored for each population. Evidence for differentiation in both sets of traits was found, but patterns of differentiation of morphological traits did not coincide with microsatellite differentiation. Morphological variation was correlated with ecological conditions at the site of origin. Analysis of molecular variance revealed significant genetic variation among populations (P < 0.001), both at the nuclear and plastid levels. There was a slight, but significant, correlation between nuclear genetic distance and geographic distance. The relatively high genetic diversity in the populations analysed indicates that the maintenance of their evolutionary potential is possible if population sizes are maintained or increased. Low levels of haplotype diversity found within the small southernmost population (Piano Costantino) indicates that genetic erosion may increase the extinction risk for this population.     

Regional population dynamics define the local genetic structure in Sorbus torminalis

Recent changes in sylvicultural practices in Central Europe have created forests with closed canopies, and tree species preferring open and sunny forests have declined in area and abundance. This led to increased isolation of populations of many rare insect-pollinated, fleshy-fruited species with a naturally scattered distribution. To gain insight into the regional population dynamics of such species, we investigated the consequences of spatial isolation, population size and density on the genetic structure of Sorbus torminalis and simultaneously considered the relationship between fecundity and habitat quality. Genotype data for biparentally (ISSRs) and maternally inherited (cpDNA PCR-RFLPs) molecular markers were generated for 26 Swiss populations of S. torminalis. We applied analyses of molecular variance (amova) to both marker types and separated the relative contributions of pollen and seed dispersal to historical gene flow. amova detected significant differentiation among populations (Phi(ST ISSR) = 0.107; Phi(ST cpDNA) = 0.370) in both marker types. The relative rate of pollen to seed gene flow was low (r = 2.919) and significantly different from equality. Isolation by distance was weak within Eastern and Western Switzerland, although populations were substantially differentiated. Within-population molecular variance was not explained by population size, whereas habitat quality (openness) positively influenced the percentage of fruiting trees and the degree of fruiting per tree, indicating that more open forests enhance sexual reproduction. Our findings of significant genetic differentiation in the absence of clear geographical structuring can be explained by the distinct ecology of S. torminalis and nondirectional colonization events in metapopulation-like dynamics.     

Autosomal, mtDNA, and Y-chromosome diversity in Amerinds: pre- and post-Columbian patterns of gene flow in South America

To evaluate sex-specific differences in gene flow between Native American populations from South America and between those populations and recent immigrants to the New World, we examined the genetic diversity at uni- and biparental genetic markers of five Native American populations from Colombia and in published surveys from native South Americans. The Colombian populations were typed for five polymorphisms in mtDNA, five restriction sites in the beta-globin gene cluster, the DQA1 gene, and nine autosomal microsatellites. Elsewhere, we published results for seven Y-chromosome microsatellites in the same populations. Autosomal polymorphisms showed a mean G(ST) of 6.8%, in agreement with extensive classical marker studies of South American populations. MtDNA and Y-chromosome markers resulted in G(ST) values of 0.18 and 0.165, respectively. When only Y chromosomes of confirmed Amerind origin were used in the calculations (as defined by the presence of allele T at locus DYS199), G(ST) increased to 0.22. G(ST) values calculated from published data for other South American natives were 0.3 and 0.29 for mtDNA and Amerind Y chromosomes, respectively. The concordance of these estimates does not support an important difference in migration rates between the sexes throughout the history of South Amerinds. Admixture analysis of the Colombian populations suggests an asymmetric pattern of mating involving mostly immigrant men and native women.     

Genetic structure and mating system of Euterpe edulis Mart. Populations: a comparative analysis using microsatellite and allozyme markers

A comparative study between microsatellite and allozyme markers was conducted on the genetic structure and mating system in natural populations of Euterpe edulis Mart. Three cohorts, including seedlings, saplings, and adults, were examined in 4 populations using 10 allozyme loci and 10 microsatellite loci. As expected, microsatellite markers had a much higher degree of polymorphism than allozymes, but estimates of multilocus outcrossing rate ( = 1.00), as well as estimates of genetic structure (F(IS), G(ST)), were similar for the 2 sets of markers. Estimates of R(ST), for microsatellites, were higher than those of G(ST), but results of both statistics revealed a close agreement for the genetic structure of the species. This study provides support for the important conclusion that allozymes are still useful and reliable markers to estimate population genetic parameters. Effects of sample size on estimates from hypervariable loci are also discussed in this paper.     

Genetic isolation of fragmented populations is exacerbated by drift and selection

Reduced genetic variation at marker loci in small populations has been well documented, whereas the relationship between quantitative genetic variation and population size has attracted little empirical investigation. Here we demonstrate that both neutral and quantitative genetic variation are reduced in small populations of a fragmented plant metapopulation, and that both drift and selective change are enhanced in small populations. Measures of neutral genetic differentiation (F(ST)) and quantitative genetic differentiation (Q(ST)) in two traits were higher among small demes, and Q(ST) between small populations exceeded that expected from drift alone. This suggests that fragmented populations experience both enhanced genetic drift and divergent selection on phenotypic traits, and that drift affects variation in both neutral markers and quantitative traits. These results highlight the need to integrate natural selection into conservation genetic theory, and suggests that small populations may represent reservoirs of genetic variation adaptive within a wide range of environments.     

G'ST and D do not replace FST

The genetic differentiation among populations is affected by mutation as well as by migration, drift and selection. For loci with high mutation rates, such as microsatellites, the amount of mutation can influence the values of indices of differentiation such as G(ST) and F(ST). For many purposes, this effect is undesirable, and as a result, new indices such as G'(ST) and D have been proposed to measure population differentiation. This paper shows that these new indices are not effective measures of the causes or consequences of population structure. Both G'(ST) and D depend heavily on mutation rate, but both are insensitive to any population genetic process when the mutation rate is high relative to the migration rate. Furthermore, D is specific to the locus being measured, and so little can be inferred about the population demography from D. However, at equilibrium, D may provide an index of whether a particular marker is more strongly affected by mutation than by migration. I argue that F(ST) is a more important summary of the effects of population structure than D and that R(ST) or other measures that explicitly account for the mutation process are much better than G(ST), G'(ST), or D for highly mutable markers. Markers with lower mutation rates will often be easier to interpret.     

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

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

Quantitative and molecular genetic variation in sympatric populations of Medicago laciniata and M. truncatula (Fabaceae): relationships with eco-geographical factors

Medicago laciniata is restricted to south of the Mediterranean basin and it extends in Tunisia from the inferior semi-arid to Saharan stages, whereas M. truncatula is a widespread species in such areas. The genetic variability in four Tunisian sympatric populations of M. laciniata and M. truncatula was analysed using 19 quantitative traits and 20 microsatellites. We investigated the amplification transferability of 52 microsatellites developed in M. truncatula to M. laciniata. Results indicate that about 78.85% of used markers are valuable genetic markers for M. laciniata. M. laciniata displayed significantly lower quantitative differentiation among populations (QST=0.12) than did M. truncatula (QST=0.45). However, high molecular differentiations, with no significant difference, were observed in M. laciniata (FST=0.48) and M. truncatula (FST=0.47). Several quantitative traits exhibited significantly smaller QST than FST for M. laciniata, consistent with constraining selection. For M. truncatula, the majority of traits displayed no statistical difference in the level of QST and FST. Furthermore, these traits are significantly associated with eco-geographical factors, consistent with selection for local adaptation rather than genetic drift. In both species, there was no significant correlation between genetic variation at quantitative traits and molecular markers. The site-of-origin explains about 5.85% and 11.27% of total quantitative genetic variability among populations of M. laciniata and M. truncatula, respectively. Established correlations between quantitative traits and eco-geographical factors were generally more moderate for M. laciniata than for M. truncatula, suggesting that the two species exhibit different genetic bases of local adaptation to varying environmental conditions. Nevertheless, no consistent patterns of associations were found between gene diversity (He) and environmental factors in either species.     

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

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

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.