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

Trait expression and signatures of adaptation in response to nitrogen addition in the common wetland plant Juncus effusus

Wetland ecosystems are known to mitigate high nutrient loadings and thus can improve water quality and prevent potential biodiversity loss caused by eutrophication. Plant traits affect wetland processes directly through effects on accumulation or metabolization of substances, and indirectly by affecting microbial transformation processes in the soil. Understanding the causes and consequences of intraspecific variation in plant functional traits and associated ecosystem processes can aid applied ecological approaches such as wetland restoration and construction. Here we investigated molecular variation and phenotypic variation in response to three levels of nitrogen availability for a regional set of populations of the common wetland plant Juncus effusus. We asked whether trait expression reveals signatures of adaptive differentiation by comparing genetic differentiation in quantitative traits and neutral molecular markers (Q_ST—F_ST comparisons) and relating trait variation to soil conditions of the plant’s origin. Molecular analyses showed that samples clustered into three very distinct genetic lineages with strong population differentiation within and among lineages. Differentiation for quantitative traits was substantial but did not exceed neutral expectations when compared across treatments or for each treatment and lineage separately. However, variation in trait expression could be explained by local soil environmental conditions of sample origin, e.g. for aboveground carbon-to-nitrogen (C:N) ratios, suggesting adaptive differentiation to contribute to trait expression even at regional level.

Epigenetic population differentiation in field- and common garden-grown Scabiosa columbaria plants

Populations often differ in phenotype and these differences can be caused by adaptation by natural selection, random neutral processes, and environmental responses. The most straightforward way to divide mechanisms that influence phenotypic variation is heritable variation and environmental‐induced variation (e.g., plasticity). While genetic variation is responsible for most heritable phenotypic variation, part of this is also caused by nongenetic inheritance. Epigenetic processes may be one of the underlying mechanisms of plasticity and nongenetic inheritance and can therefore possibly contribute to heritable differences through drift and selection. Epigenetic variation may be influenced directly by the environment, and part of this variation can be transmitted to next generations. Field screenings combined with common garden experiments will add valuable insights into epigenetic differentiation, epigenetic memory and can help to reveal part of the relative importance of epigenetics in explaining trait variation. We explored both genetic and epigenetic diversity, structure and differentiation in the field and a common garden for five British and five French Scabiosa columbaria populations. Genetic and epigenetic variation was subsequently correlated with trait variation. Populations showed significant epigenetic differentiation between populations and countries in the field, but also when grown in a common garden. By comparing the epigenetic variation between field and common garden‐grown plants, we showed that a considerable part of the epigenetic memory differed from the field‐grown plants and was presumably environmentally induced. The memory component can consist of heritable variation in methylation that is not sensitive to environments and possibly genetically based, or environmentally induced variation that is heritable, or a combination of both. Additionally, random epimutations might be responsible for some differences as well. By comparing epigenetic variation in both the field and common environment, our study provides useful insight into the environmental and genetic components of epigenetic variation.

Phenotypic Variation and Sexual Size Dimorphism in Dichroplus elongatus (Orthoptera: Acrididae)

Patterns of body size evolution are of particular interest because body size can affect virtually all the physiological and life history traits of an organism. Sexual size dimorphism (SSD), a difference in body size between males and females, is a widespread phenomenon in insects. Much of the variation in SSD is genetically based and likely due to differential selection acting on males and females. The importance of environmental variables and evolutionary processes affecting phenotypeic variation in both sexes may be useful to gain insights into insect ecology and evolution. Dichroplus elongatus Giglio-Tos is a South American grasshopper widely distributed throughout Argentina, Uruguay, most of Chile, and southern Brazil. In this study, we analyzed 122 adult females of D. elongatus collected in eight natural populations from central-east Argentina. Females show large body size variation among the analyzed populations and this variation exhibits a strong relationship with fecundity. Our results have shown that larger females were more fecund than smaller ones. We found that ovariole number varied along a latitudinal gradient, with higher ovariole numbers in populations from warmer locations. A considerable female-biased SSD was detected. SSD for three analyzed morphometric traits scaled isometrically. However, SSD for thorax length displayed a considerable variation across the studied area, indicating a larger relative increase in female size than in male size in warmer environmental conditions.

Commentary: When does understanding phenotypic evolution require identification of the underlying genes?

Adaptive evolution is fundamentally a genetic process. Over the past three decades, characterizing the genes underlying adaptive phenotypic change has revealed many important aspects of evolutionary change. At the same time, natural selection is often fundamentally an ecological process that can often be studied without identifying the genes underlying the variation on which it acts. This duality has given rise to disagreement about whether, and under what circumstances, it is necessary to identify specific genes associated with phenotypic change. This issue is of practical concern, especially for researchers who study nonmodel organisms, because of the often enormous cost and labor required to "go for the genes." We here consider a number of situations and questions commonly addressed by researchers. Our conclusion is that although gene identification can be crucial for answering some questions, there are others for which definitive answers can be obtained without finding underlying genes. It should thus not be assumed that considerations of "empirical completeness" dictate that gene identification is always desirable.

Divergent selection along climatic gradients in a rare central European endemic species, Saxifraga sponhemica

BACKGROUND AND AIMS

The effects of habitat fragmentation on quantitative genetic variation in plant populations are still poorly known. Saxifraga sponhemica is a rare endemic of Central Europe with a disjunct distribution, and a stable and specialized habitat of treeless screes and cliffs. This study therefore used S. sponhemica as a model species to compare quantitative and molecular variation in order to explore (1) the relative importance of drift and selection in shaping the distribution of quantitative genetic variation along climatic gradients; (2) the relationship between plant fitness, quantitative genetic variation, molecular genetic variation and population size; and (3) the relationship between the differentiation of a trait among populations and its evolvability.

METHODS

Genetic variation within and among 22 populations from the whole distribution area of S. sponhemica was studied using RAPD (random amplified polymorphic DNA) markers, and climatic variables were obtained for each site. Seeds were collected from each population and germinated, and seedlings were transplanted into a common garden for determination of variation in plant traits.

KEY RESULTS

In contrast to previous results from rare plant species, strong evidence was found for divergent selection. Most population trait means of S. sponhemica were significantly related to climate gradients, indicating adaptation. Quantitative genetic differentiation increased with geographical distance, even when neutral molecular divergence was controlled for, and QST exceeded FST for some traits. The evolvability of traits was negatively correlated with the degree of differentiation among populations (QST), i.e. traits under strong selection showed little genetic variation within populations. The evolutionary potential of a population was not related to its size, the performance of the population or its neutral genetic diversity. However, performance in the common garden was lower for plants from populations with reduced molecular genetic variation, suggesting inbreeding depression due to genetic erosion.

CONCLUSIONS

The findings suggest that studies of molecular and quantitative genetic variation may provide complementary insights important for the conservation of rare species. The strong differentiation of quantitative traits among populations shows that selection can be an important force for structuring variation in evolutionarily important traits even for rare endemic species restricted to very specific habitats.

Molecular and quantitative genetic differentiation in Sitobion avenae populations from both sides of the Qinling Mountains

Quantitative trait differences are often assumed to be correlated with molecular variation, but the relationship is not certain, and empirical evidence is still scarce. To address this issue, we sampled six populations of the cereal aphid Sitobion avenae from areas north and south of the Qinling Mountains, and characterized their molecular variation at seven microsatellite loci and quantitative variation at nine life-history traits. Our results demonstrated that southern populations had slightly longer developmental times of nymphs but much higher lifetime fecundity, compared to northern populations. Of the nine tested quantitative characters, eight differed significantly among populations within regions, as well as between northern and southern regions. Genetic differentiation in neutral markers was likely to have been caused by founder events and drift. Increased subdivision for quantitative characters was found in northern populations, but reduced in southern populations. This phenomenon was not found for molecular characters, suggesting the decoupling between molecular and quantitative variation. The pattern of relationships between FST and QST indicated divergent selection and suggested that local adaptation play a role in the differentiation of life-history traits in tested S. avenae populations, particularly in those traits closely related to reproduction. The main role of natural selection over genetic drift was also supported by strong structural differences in G-matrices among S. avenae populations. However, cluster analyses did not result in two groups corresponding to northern and southern regions. Genetic differentiation between northern and southern populations in neutral markers was low, indicating considerable gene flow between them. The relationship between molecular and quantitative variation, as well as its implications for differentiation and evolution of S. avenae populations, was discussed.

Phenotypic divergence during the invasion of Phyla canescens in Australia and France: evidence for selection-driven evolution

Rapid adaptive evolution has been advocated as a mechanism that promotes invasion. Demonstrating adaptive evolution in invasive species requires rigorous analysis of phenotypic shifts driven by selection. Here, we document selection-driven evolution of Phyla canescens, an Argentine weed, in two invaded regions (Australia and France). Invasive populations possessed similar or higher diversity than native populations, and displayed mixed lineages from different sources, suggesting that genetic bottlenecks in both countries might have been alleviated by multiple introductions. Compared to native populations, Australian populations displayed more investment in sexual reproduction, whereas French populations possessed enhanced vegetative reproduction and growth. We partitioned evolutionary forces (selection vs. stochastic events) using two independent methods. Results of both analyses suggest that the pattern of molecular and phenotypic variability among regions was consistent with selection-driven evolution, rather than stochastic events. Our findings indicate that selection has shaped the evolution of P. canescens in two different invaded regions.

Elucidating the role of genetic drift and natural selection in cork oak differentiation regarding drought tolerance

Drought is the main selection agent in Mediterranean ecosystems and it has been suggested as an important evolutionary force responsible for population diversification in these types of environments. However, population divergence in quantitative traits can be driven by either natural selection, genetic drift or both. To investigate the roles of these forces on among-population divergence in ecophysiological traits related to drought tolerance (carbon isotope discrimination, specific leaf area, leaf size and leaf nitrogen content), we compared molecular and quantitative genetic differentiation in a common garden experiment including thirteen cork oak (Quercus suber L.) populations across a gradient of rainfall and temperature. Population differentiation for height, specific leaf area, leaf size and nitrogen leaf content measured during a dry year far exceeded the molecular differentiation measured by six nuclear microsatellites. Populations from dry-cool sites showed the lowest nitrogen leaf content and the smallest and thickest leaves contrasting with those from humid-warm sites. These results suggest (i) these traits are subjected to divergence selection and (ii) the genetic differences among populations are partly due to climate adaptation. By contrast, the low among-population divergence found in basal diameter, annual growth and carbon isotopic discrimination (a surrogate for water use efficiency) suggests low or no divergence selection for these traits. Among-population differentiation for neutral markers was not a good predictor for differentiation regarding the quantitative traits studied here, except for leaf size. The correlation observed between the genetic differentiation for leaf size and that for molecular markers was exclusively due to the association between leaf size and the microsatellite QpZAG46, which suggests a possible linkage between QpZAG46 and genes encoding for leaf size.

Speciation with gene flow in the large white-headed gulls: does selection counterbalance introgression?

We investigated the role of selection in generating and maintaining species distinctness in spite of ongoing gene flow, using two zones of secondary contact between large gull species in Europe (Larus argentatus and Larus cachinnans) and North America (Larus glaucescens and Larus occidentalis). We used the pattern of neutral genetic differentiation at nine microsatellite loci (F(ST)) as an indicator of expected changes under neutral processes and compared it with phenotypic differentiation (P(ST)) for a large number of traits (size, plumage melanism and coloration of bare parts). Even assuming very low heritability, interspecific divergence between L. glaucescens and L. occidentalis in plumage melanism and orbital ring colour clearly exceeded neutral differentiation. Similarly, melanism of the central primaries was highly divergent between L. argentatus and L. cachinnans. Such divergence is unlikely to have arisen randomly and is therefore attributed to spatially varying selection. Variation in plumage melanism in both transects agrees with Gloger's rule, which suggests that latitude (and associated sun and humidity gradients) could be the selective pressure shaping differentiation in plumage melanism. We suggest that strong species differentiation in orbital ring colour results from sexual selection. We conclude that these large gull species, along with other recently diverged species that hybridize after coming into secondary contact, may differ only in restricted regions of the genome that are undergoing strong disruptive selection because of their phenotypic effects.

F(ST) and Q(ST) under neutrality

A commonly used test for natural selection has been to compare population differentiation for neutral molecular loci estimated by F(ST) and for the additive genetic component of quantitative traits estimated by Q(ST). Past analytical and empirical studies have led to the conclusion that when averaged over replicate evolutionary histories, Q(ST) = F(ST) under neutrality. We used analytical and simulation techniques to study the impact of stochastic fluctuation among replicate outcomes of an evolutionary process, or the evolutionary variance, of Q(ST) and F(ST) for a neutral quantitative trait determined by n unlinked diallelic loci with additive gene action. We studied analytical models of two scenarios. In one, a pair of demes has recently been formed through subdivision of a panmictic population; in the other, a pair of demes has been evolving in allopatry for a long time. A rigorous analysis of these two models showed that in general, it is not necessarily true that mean Q(ST) = F(ST) (across evolutionary replicates) for a neutral, additive quantitative trait. In addition, we used finite-island model simulations to show there is a strong positive correlation between Q(ST) and the difference Q(ST) - F(ST) because the evolutionary variance of Q(ST) is much larger than that of F(ST). If traits with relatively large Q(ST) values are preferentially sampled for study, the difference between Q(ST) and F(ST) will also be large and positive because of this correlation. Many recent studies have used tests of the null hypothesis Q(ST) = F(ST) to identify diversifying or uniform selection among subpopulations for quantitative traits. Our findings suggest that the distributions of Q(ST) and F(ST) under the null hypothesis of neutrality will depend on species-specific biology such as the number of subpopulations and the history of subpopulation divergence. In addition, the manner in which researchers select quantitative traits for study may introduce bias into the tests. As a result, researchers must be cautious before concluding that selection is occurring when Q(ST) not equal F(ST).

Assaying polymorphism at DNA level for genetic diversity diagnostics of the safflower (Carthamus tinctorius L.) world germplasm resources

Carthamus tinctorius (2n = 2x = 24), commonly known as safflower, is widely cultivated in agricultural production systems of Asia, Europe, Australia, and the Americas as a source of high quality vegetable and industrial oil. Twenty-two RAPD primers, 18 SSR primers, and 10 AFLP primer combinations were used to assess: (1) the genetic diversity of 85 accessions (originating from 24 countries) representing global germplasm variability of safflower and (2) the interrelationships among safflower 'centers of similarity' or 'regional gene pools' proposed earlier. The RAPD and SSR primers and AFLP primer combinations revealed 57.6, 68.0, and 71.2% polymorphism, respectively, among 111, 72, and 330 genetic loci amplified from the accessions. The sum of effective number of alleles (66.44), resolving power (59.16), and marker index (51.3) explicitly revealed the relative superiority of AFLP as a marker system in uncovering variation in safflower. Overall, AFLP markers could recognize 'centers of similarity' or 'regional gene pools'. Analysis of molecular variance and Shannon's information index provided corroborating evidences for the present and previous studies that concluded fragmentation of safflower gene pool into many gene pools. Divergent directional selection is likely to have played an important role in shaping the diversity. From the practical applications standpoint, the diversity of Iran-Afghanistan gene pool is very high, equivalent to the total diversity of the species. The Far East gene pool is the least diverse. The present comprehensive input, first of its own kind in safflower, will assist marker based improvement programmes in the crop.

Adaptive responses of independent traits to the same environmental gradient in the intertidal snail Bembicium vittatum

The snail Bembicium vittatum occupies a wide range of intertidal habitats in the Houtman Abrolhos Islands, Western Australia. Allozyme variation reflects patterns of connectivity, which are independent of local habitat. In contrast, heritable differences in shell shape among 83 shore sites vary with habitat, indicating local adaptation. Here we examine dimorphisms of colour and spotting of the shell in the same populations, as a test of consistency and complexity of patterns of local adaptation. Within populations, the frequency of spotted shells is higher in dark shells. Despite this association, spatial variations of colour and spotting are only weakly correlated. As predicted for traits associated with local adaptation, subdivision is greater for colour, spotting and shape than for allozymes. Colour and shape are associated with local habitat, such that populations on vertical shores have higher frequencies of dark and relatively flatter shells than those on gently sloping shores. These associations are repeatable between three separate groups of islands. Spotting shows a weaker, but significant association with the same gradient. Although shape does not differ between colour morphs within populations, the proportion of dark shells is strongly associated with shape. Thus, the independent shell traits are apparently adapted to a common, biologically significant gradient, even though the adaptive mechanisms probably differ for colour and shape. The parallel variations of independent traits highlight both the complexity of local adaptation and the potential to reveal evolutionarily significant environmental contrasts by examining adaptively relevant traits.

Genetic variation and population structure in black-grass (Alopecurus myosuroides Huds.), a successful, herbicide-resistant, annual grass weed of winter cereal fields

Black-grass (Alopecurus myosuroides) is an allogamous grass weed common in cereal fields of northern Europe, which developed resistance to a widely used family of herbicides, the ACCase-inhibiting herbicides. Resistance is caused by mutations at the ACCase gene and other, metabolism-based, mechanisms. We investigated the genetic structure of 36 populations of black-grass collected in one region of France (Côte d'Or), using 116 amplified fragment length polymorphism (AFLP) loci and sequence data at the ACCase gene. The samples were characterized for their level of herbicide resistance and genotyped for seven known ACCase mutations conferring resistance. All samples contained herbicide-resistant plants, and 19 contained ACCase mutations. The genetic diversity at AFLP loci was high (H(T) = 0.246), while differentiation among samples was low (F(ST) = 0.023) and no isolation by distance was detected. Genetic diversity within samples did not vary with the frequency of herbicide resistance. A Bayesian algorithm was used to infer population structure. The two genetic clusters inferred were not associated with any geographical structure or with herbicide resistance. A high haplotype diversity (H(d) = 0.873) and low differentiation (G(ST) = 0.056) were observed at ACCase. However, haplotype diversity within samples decreased with the frequency of ACCase-based resistance. We suggest that the genetic structure of black-grass is affected by its recent expansion as a weed. Our data demonstrate that the strong selection imposed by herbicides did not modify the genome-wide genetic structure of an allogamous weed that probably has large effective population sizes. Our study gives keys to a better understanding of the evolution of successful, noxious weeds in modern agriculture.

High quantitative and no molecular differentiation of a freshwater snail (Galba truncatula) between temporary and permanent water habitats

We investigate the variation in quantitative and molecular traits in the freshwater snail Galba truncatula, from permanent and temporary water habitats. Using a common garden experiment, we measured 20 quantitative traits and molecular variation using seven microsatellites in 17 populations belonging to these two habitats. We estimated trait means in each habitat. We also estimated the distributions of overall genetic quantitative variation (QST), and of molecular variation (FST), within and between habitats. Overall, we observed a lack of association between molecular and quantitative variance. Among habitats, we found QST>FST, an indication of selection for different optima. Individuals from temporary water habitat matured older, at a larger size and were less fecund than individuals from permanent water habitat. We discuss these findings in the light of several theories for life-history traits evolution.

Changes in mating system and social structure of the ant Petalomyrmex phylax are associated with range expansion in Cameroon

Past climate shifts have led to major oscillations in species distributions. Hence historical contingencies and selective processes occurring during such phases may be determinants for understanding the forces that have shaped extant phenotypes. In the plant-ant Petalomyrmex phylax (Formicinae), we observed spatial variation in number of queens in mature colonies, from several queens (high polygyny) in the median part of its distribution to a moderate number of queens (weak polygyny) or even only a single queen (monogyny) in the southwesternmost populations. This variation did not correlate with indicators of variation in current nest site availability and colony turnover, the supposedly determinant selective forces acting on gyny in ants. We show here that the variation in social structure correlates with a historical process corresponding to a progressive colonization of coastal southern Cameroon by the ant. Using microsatellite markers, we observed a clear pattern of isolation by distance except for the southernmost populations. Measures of genetic variability that do not take into account allele size were at equilibrium in all except the southernmost populations, suggesting recent foundation of the latter. Measures of genetic diversity taking into account allele size showed a clinal north-south decrease in variance of allele size. We propose that southern populations have yet to regain allele size variance after bottlenecks associated with the foundation of new populations, and that this variance is regained over time. Hence variation in social structure mirrors an old but still active southward colonization process or metapopulation dynamics, possibly in association with an expansion of the rain forest habitat during the late Holocene. A low number of queens in ant colonies is typically associated with strong dispersal capacity. We therefore suggest that the initial founders of new populations belong to the monogynous to weakly polygynous phenotype, and that queen number progressively increases in older populations.

Population decline despite high genetic diversity in the new allopolyploid species Senecio cambrensis (Asteraceae)

Senecio cambrensis (Welsh groundsel) is a new allohexaploid species, which originated in Wales, UK, in the early part of the 20th century following hybridization between the native tetraploid groundsel (Senecio vulgaris) and the introduced diploid Oxford ragwort (Senecio squalidus). A survey of the number of populations and flowering individuals per population of S. cambrensis in Wales was conducted at peak flowering time in June 2002, 2003 and 2004. The results show a dramatic decrease in both population number and population size of the species since the 1980s when the last population census was conducted. A survey of amplified fragment length polymorphism (AFLP) variation showed that this decline has occurred despite the fact that S. cambrensis contains a high level of genetic diversity with each individual screened possessing a unique multilocus phenotype. The level of variance within the species was similar to that found in one parent (S. vulgaris) and slightly greater than that among samples of the other parent (S. squalidus). Only a small proportion (5%) of AFLP diversity was partitioned among populations indicating a lack of population structure and possibly high levels of gene flow via seed dispersal in what is predominantly a selfing species. Senecio cambrensis showed closer similarity in AFLP phenotype to S. vulgaris than to S. squalidus. Possible causes of this and also the high level of AFLP diversity found in S. cambrensis are discussed. It is suggested that intergenomic recombination following occasional multivalent formation during meiosis in S. cambrensis is likely to be an important cause of both phenomena, although other causes are not ruled out.

Lack of genetic structure among ecologically adapted populations of an Australian rainforestDrosophilaspecies as indicated by microsatellite markers and mitochondrial DNA sequences

Although fragmented rainforest environments represent hotspots for invertebrate biodiversity, few genetic studies have been conducted on rainforest invertebrates. Thus, it is not known if invertebrate species in rainforests are highly genetically fragmented, with the potential for populations to show divergent selection responses, or if there are low levels of gene flow sufficient to maintain genetic homogeneity among fragmented populations. Here we use microsatellite markers and DNA sequences from the mitochondrial ND5 locus to investigate genetic differences among Drosophila birchii populations from tropical rainforests in Queensland, Australia. As found in a previous study, mitochondrial DNA diversity was low with no evidence for population differentiation among rainforest fragments. The pattern of mitochondrial haplotype variation was consistent with D. birchii having undergone substantial past population growth. Levels of nuclear genetic variation were high in all populations while F(ST) values were very low, even for flies from geographically isolated areas of rainforest. No significant differentiation was observed between populations on either side of the Burdekin Gap (a long-term dry corridor), although there was evidence for higher gene diversity in low-latitude populations. Spatial autocorrelation coefficients were low and did not differ significantly from random, except for one locus which revealed a clinal-like pattern. Comparisons of microsatellite differentiation contrasted with previously established clinal patterns in quantitative traits in D. birchii, and indicate that the patterns in quantitative traits are likely to be due to selection. These results suggest moderate gene flow in D. birchii over large distances. Limited population structure in this species appears to be due to recent range expansions or cycles of local extinctions followed by recolonizations/expansions. Nevertheless, patterns of local adaptation have developed in D. birchii that may result in populations showing different selection responses when faced with environmental change.

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.

The effects of dominance, regular inbreeding and sampling design on Q(ST), an estimator of population differentiation for quantitative traits

To test whether quantitative traits are under directional or homogenizing selection, it is common practice to compare population differentiation estimates at molecular markers (F(ST)) and quantitative traits (Q(ST)). If the trait is neutral and its determinism is additive, then theory predicts that Q(ST) = F(ST), while Q(ST) > F(ST) is predicted under directional selection for different local optima, and Q(ST) < F(ST) is predicted under homogenizing selection. However, nonadditive effects can alter these predictions. Here, we investigate the influence of dominance on the relation between Q(ST) and F(ST) for neutral traits. Using analytical results and computer simulations, we show that dominance generally deflates Q(ST) relative to F(ST). Under inbreeding, the effect of dominance vanishes, and we show that for selfing species, a better estimate of Q(ST) is obtained from selfed families than from half-sib families. We also compare several sampling designs and find that it is always best to sample many populations (>20) with few families (five) rather than few populations with many families. Provided that estimates of Q(ST) are derived from individuals originating from many populations, we conclude that the pattern Q(ST) > F(ST), and hence the inference of directional selection for different local optima, is robust to the effect of nonadditive gene actions.

Distinguishing adaptive from nonadaptive genetic differentiation: comparison of QST and FST 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.

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 QUANTITATIVE AND NEUTRAL GENETIC VARIATION IN LOCALLY ADAPTED POPULATIONS OF THE NATTERJACK TOAD, BUFO CALAMITA

The relative importance of natural selection and genetic drift in determining patterns of phenotypic diversity observed in nature is still unclear. The natterjack toad (Bufo calamita) is one of a few amphibian species capable of breeding in saline ponds, even though water salinity represents a considerable stress for them. Results from two common-garden experiments showed a pattern of geographic variation in embryonic salinity tolerance among populations from either fresh or brackish environments, consistent with the hypothesis of local adaptation. Full-sib analysis showed increased variation in survival among sibships within population for all populations as osmotic stress was increased (broad-sense heritability increased as salinity raised). Nevertheless, toads native to the brackish water environment had the highest overall survival under brackish conditions. Levels of population genetic differentiation for salinity tolerance were higher than those of neutral genetic differentiation, the latter obtained through the analysis of eight microsatellite loci. Microsatellite markers also revealed little population differentiation, lack of an isolation-by-distance pattern, and moderate gene flow connecting the populations. Therefore, environmental stress tolerance appears to have evolved in absence of geographic isolation, and consequently we reject the null hypothesis of neutral differentiation.

Intra-specific variability of ozone sensitivity in Centaurea jacea L., a potential bioindicator for elevated ozone concentrations

Brown knapweed (Centaurea jacea L.) has been suggested as a potential bioindicator for tropospheric ozone (O3), but little is known about the intra-specific variation in O3 sensitivity in this wild species. The aim of this study was to quantify the differences in O3 sensitivity among and within five populations, and to relate the differences to morphological, phenological, and genetic characteristics. These parameters were periodically recorded in two consecutive experiments on a total of 357 plants from five different European countries (Norway, Hungary, Switzerland, Italy, Slovenia). They were grown from seed in natural soil under ambient conditions at a site with seasonally elevated O3 concentrations (Cadenazzo, southern Switzerland). The populations differed significantly both in frequency and extent of O3 injury, as well as in phenological development. The observed degree of O3 injury was highest in the Slovenian and the Swiss populations, while only few Hungarian and Norwegian plants showed slight symptoms of injury. Plants were generally most sensitive to O3 when reaching the reproductive stage, and insensitive at the rosette stage. Amplified fragment length polymorphism analysis (AFLP) demonstrated genetic distinctiveness of the five C. jacea populations. All individuals of four of the five populations were correctly assigned to the respective populations based on principal component analysis. Cluster analysis quite accurately reflected the geographic origin of each population. Overall, the analysis revealed a high degree of intra-specific variability in O3 sensitivity in C. jacea, and underlined the important influence of the climate-dependent population-specific plant development on O3 sensitivity. These observations may constrain the development of a standardized biomonitoring system.

POPULATION DIFFERENTIATION IN G MATRIX STRUCTURE DUE TO NATURAL SELECTION IN RANA TEMPORARIA

The additive genetic variance-covariance matrix (G) is a concept central to discussions about evolutionary change over time in a suite of traits. However, at the moment we do not know how fast G itself changes as a consequence of selection or how sensitive it is to environmental influences. We investigated possible evolutionary divergence and environmental influences on G using data from a factorial common-garden experiment where common frog (Rana temporaria) tadpoles from two divergent populations were exposed to three different environmental treatments. G-matrices were estimated using an animal model approach applied to data from a NCII breeding design. Matrix comparisons using both Flury and multivariate analysis of variance methods revealed significant differences in G matrices both between populations and between treatments within populations, the former being generally larger than the latter. Comparison of levels of population differentiation in trait means using Q(ST) indices with that observed in microsatellite markers (F(ST)) revealed that the former values generally exceeded the neutral expectation set by F(ST). Hence, the results suggest that intraspecific divergence in G matrix structure has occurred mainly due to natural selection.

Correlation between molecular markers and adaptively significant genetic variation in Bromus tectorum (Poaceae), an inbreedingannual grass

Single sequence repeat (SSR) and amplified fragment length polymorphic (AFLP) molecular marker genotypes in cheatgrass (Bromus tectorum) were compared to published data on phenotypic variation in seed dormancy, vernalization requirement, and resistance to the pathogen Ustilago bullata. Several features of cheatgrass facilitated this study: it is a recent invader in the western United States, has considerable phenotypic polymorphism, and is an obligate self-pollinator. Forty self-pollinating lines from four populations common to the three phenotypic data sets were analyzed for molecular genetic variation using seven SSR loci and 31 AFLP loci. We examined correlations between distance matrices using the Mantel test for each pair of studies. The two molecular data sets were significantly correlated (r = 0.636). The AFLP markers often distinguished among several lines with identical SSR genotypes. The AFLP data were also significantly correlated with the phenotypic data (r values from 0.4640 to 0.5658), but the SSR data were much more highly correlated (r values from 0.677 to 0.844). The difference between molecular marker systems was especially notable when an outlier population from Potosi Pass, Nevada, was excluded from the analysis. These results suggest that SSR markers may be good surrogates for phenotypic traits in population genetic studies of strongly inbreeding species such as cheatgrass.

Low genetic differentiation among seasonal cohorts in Senecio vulgaris as revealed by amplified fragment length polymorphism analysis

Common groundsel, Senecio vulgaris (Asteraceae), is a highly selfing semelparous ephemeral weed that belongs to the few plant species in central Europe capable of growing, flowering and fruiting all year round. In temperate climates, flowering S. vulgaris cohorts were found to appear up to three times per year. Using amplified fragment length polymorphism (AFLP) molecular markers we examined temporal genetic differentiation among spring, summer and autumn cohorts at each of seven sites located in two regions in Switzerland. Strong genetic differentiation among cohorts may indicate the existence of seasonal races of S. vulgaris, reproductively isolated by nonoverlapping flowering phenologies. Analysis of molecular variance (amova) revealed that < 2.5% of the AFLP variation resided among cohorts within sites, whereas there was significant genetic differentiation among plants from different sites (15.6%) and among individuals within cohorts (81.9%). Significant genetic differentiation was also observed between the two regions. Isolation-by-distance was found on a regional scale, but not on a local scale. Gene flow was estimated to be approximately 15-fold higher among cohorts within sites than among sites. We further found, on average, similar levels of genetic diversity within the three seasonal cohorts. The results of this study demonstrate that season of growth represents a weak barrier for genetic exchange among S. vulgaris populations and does not affect molecular variance. Therefore, there is no evidence for the existence of seasonally specialized races of S. vulgaris. We discuss some implications of the results for the biological control of S. vulgaris using a native rust fungus.

MOLECULAR AND QUANTITATIVE TRAIT VARIATION WITHIN AND AMONG POPULATIONS OF THE INTERTIDAL COPEPOD TIGRIOPUS CALIFORNICUS

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.

Latitudinal divergence of common frog (Rana temporaria) life history traits by natural selection: evidence from a comparison of molecular and quantitative genetic data

The relative roles of natural selection and direct environmental induction, as well as of natural selection and genetic drift, in creating clinal latitudinal variation in quantitative traits have seldom been assessed in vertebrates. To address these issues, we compared molecular and quantitative genetic differentiation between six common frog (Rana temporaria) populations along an approximately 1600 km long latitudinal gradient across Scandinavia. The degree of population differentiation (QST approximately 0.81) in three heritable quantitative traits (age and size at metamorphosis, growth rate) exceeded that in eight (neutral) microsatellite loci (FST = 0.24). Isolation by distance was clear for both neutral markers and quantitative traits, but considerably stronger for one of the three quantitative traits than for neutral markers. QST estimates obtained using animals subjected to different rearing conditions (temperature and food treatments) revealed some environmental dependency in patterns of population divergence in quantitative traits, but in general, these effects were weak in comparison to overall patterns. Pairwise comparisons of FST and QST estimates across populations and treatments revealed that the degree of quantitative trait differentiation was not generally predictable from knowledge of that in molecular markers. In fact, both positive and negative correlations were observed depending on conditions where the quantitative genetic variability had been measured. All in all, the results suggest a very high degree of genetic subdivision both in neutral marker genes and genes coding quantitative traits across a relatively recently (< 9000 years) colonized environmental gradient. In particular, they give evidence for natural selection being the primary agent behind the observed latitudinal differentiation in quantitative traits.

A format for databasing and comparison of AFLP fingerprint profiles

BACKGROUND

Amplified fragment length polymorphism (AFLP) is a PCR-based technique that involves restriction of genomic DNA followed by ligation of adaptors to the fragments generated and selective PCR amplification of a subset of these fragments. The amplified fragments are separated on a sequencing gel and visualized by autoradiography or fluorescent sequencing equipment. AFLP allows high-resolution genotyping but the lack of a format for databasing and comparison of AFLP fingerprint profiles limits its wider applications in profiling large numbers of biological samples.

RESULTS

A scheme is described to represent a DNA fingerprint profile with a nucleotide sequence-like format in which the information line contains the minimal necessary details to interpret an AFLP DNA fingerprint profile. They include technique used, information on restriction enzymes, primer combination, biological source for DNA materials, fragment sizing and annotation. The bodylines contain information on size and relative intensity of DNA fragments by a string of defined alphabets or symbols. Algorithms for normalizing raw data, binning of fragments and comparing AFLP DNA fingerprint profiles are described. Firstly, the peak heights are normalized against their average and then represented by five symbols according to their relative intensities. Secondly, a binning algorithm based loosely on common springs and rubber bands is applied, which positions sequence fragments into their best possible integer approximations. A BLAST-like reward-penalty concept is used to compare AFLP fingerprint profiles by matching peaks using two metrics: score and percentage of similarity. A software package was developed based on our scheme and proposed algorithms. Example of use this software is given in evaluating novelty of a new tropical orchid cultivar by comparing its AFLP fingerprint profile against those of related commercial cultivars in a database.

CONCLUSIONS

AFLP DNA fingerprint profiles can be databased and compared effectively with software developed based on our scheme and algorithms. It will facilitate wider use of this DNA fingerprinting technique in areas such as forensic study, intellectual property protection for biological materials and biodiversity management. Moreover, the same concepts can be applied to databasing and comparing DNA fingerprint profiles obtained with other DNA fingerprint techniques.