Population differentiation of polygenic score predictions under stabilizing selection

Given the many small-effect loci uncovered by genome-wide association studies (GWAS), polygenic scores have become central to genomic medicine, and have found application in diverse settings including evolutionary studies of adaptation. Despite their promise, polygenic scores have been found to suffer from limited portability across human populations. This at first seems in conflict with the observation that most common genetic variation is shared among populations. We investigate one potential cause of this discrepancy: stabilizing selection on complex traits. Counterintuitively, while stabilizing selection constrains phenotypic evolution, it accelerates the loss and fixation of alleles underlying trait variation within populations (GWAS loci). Thus even when populations share an optimum phenotype, stabilizing selection erodes the variance contributed by their shared GWAS loci, such that predictions from GWAS in one population explain less of the phenotypic variation in another. We develop theory to quantify how stabilizing selection is expected to reduce the prediction accuracy of polygenic scores in populations not represented in GWAS samples. In addition, we find that polygenic scores can substantially overstate average genetic differences of phenotypes among populations. We emphasize stabilizing selection around a common optimum as a useful null model to connect patterns of allele frequency and polygenic score differentiation. This article is part of the theme issue 'Celebrating 50 years since Lewontin's apportionment of human diversity'.

Evolution of polygenic traits under global vs local adaptation

Observations about the number, frequency, effect size, and genomic distribution of alleles associated with complex traits must be interpreted in light of evolutionary process. These characteristics, which constitute a trait's genetic architecture, can dramatically affect evolutionary outcomes in applications from agriculture to medicine, and can provide a window into how evolution works. Here, I review theoretical predictions about the evolution of genetic architecture under spatially homogeneous, global adaptation as compared with spatially heterogeneous, local adaptation. Due to the tension between divergent selection and migration, local adaptation can favor "concentrated" genetic architectures that are enriched for alleles of larger effect, clustered in a smaller number of genomic regions, relative to expectations under global adaptation. However, the evolution of such architectures may be limited by many factors, including the genotypic redundancy of the trait, mutation rate, and temporal variability of environment. I review the circumstances in which predictions differ for global vs local adaptation and discuss where progress can be made in testing hypotheses using data from natural populations and lab experiments. As the field of comparative population genomics expands in scope, differences in architecture among traits and species will provide insights into how evolution works, and such differences must be interpreted in light of which kind of selection has been operating.

Hidden paths to endless forms most wonderful: parasite-blind diversification of host quality

Evolutionary diversification can occur in allopatry or sympatry, can be driven by selection or unselected, and can be phenotypically manifested immediately or remain latent until manifested in a newly encountered environment. Diversification of host-parasite interactions is frequently studied in the context of intrinsically selective coevolution, but the potential for host-parasite interaction phenotypes to diversify latently during parasite-blind host evolution is rarely considered. Here, we use a social bacterium experimentally adapted to several environments in the absence of phage to analyse allopatric diversification of host quality-the degree to which a host population supports a viral epidemic. Phage-blind evolution reduced host quality overall, with some bacteria becoming completely resistant to growth suppression by phage. Selective-environment differences generated only mild divergence in host quality. However, selective environments nonetheless played a major role in shaping evolution by determining the degree of stochastic diversification among replicate populations within treatments. Ancestral motility genotype was also found to strongly shape patterns of latent host-quality evolution and diversification. These outcomes show that (i) adaptive landscapes can differ in how they constrain stochastic diversification of a latent phenotype and (ii) major effects of selection on biological diversification can be missed by focusing on trait means. Collectively, our findings suggest that latent-phenotype evolution should inform host-parasite evolution theory and that diversification should be conceived broadly to include latent phenotypes.

Genetic population structure constrains local adaptation in sticklebacks

Repeated and independent adaptation to specific environmental conditions from standing genetic variation is common. However, if genetic variation is limited, the evolution of similar locally adapted traits may be restricted to genetically different and potentially less optimal solutions or prevented from happening altogether. Using a quantitative trait locus (QTL) mapping approach, we identified the genomic regions responsible for the repeated pelvic reduction (PR) in three crosses between nine-spined stickleback populations expressing full and reduced pelvic structures. In one cross, PR mapped to linkage group 7 (LG7) containing the gene Pitx1, known to control pelvic reduction also in the three-spined stickleback. In the two other crosses, PR was polygenic and attributed to 10 novel QTL, of which 90% were unique to specific crosses. When screening the genomes from 27 different populations for deletions in the Pitx1 regulatory element, these were only found in the population in which PR mapped to LG7, even though the morphological data indicated large-effect QTL for PR in several other populations as well. Consistent with the available theory and simulations parameterized on empirical data, we hypothesize that the observed variability in genetic architecture of PR is due to heterogeneity in the spatial distribution of standing genetic variation caused by >2× stronger population structuring among freshwater populations and >10× stronger genetic isolation by distance in the sea in nine-spined sticklebacks as compared to three-spined sticklebacks.

Different Genomic Changes Underlie Adaptive Evolution in Populations of Contrasting History

Experimental evolution is a powerful tool to understand the adaptive potential of populations under environmental change. Here, we study the importance of the historical genetic background in the outcome of evolution at the genome-wide level. Using the natural clinal variation of Drosophila subobscura, we sampled populations from two contrasting latitudes (Adraga, Portugal and Groningen, Netherlands) and introduced them in a new common environment in the laboratory. We characterized the genome-wide temporal changes underlying the evolutionary dynamics of these populations, which had previously shown fast convergence at the phenotypic level, but not at chromosomal inversion frequencies. We found that initially differentiated populations did not converge either at genome-wide level or at candidate SNPs with signs of selection. In contrast, populations from Portugal showed convergence to the control population that derived from the same geographical origin and had been long-established in the laboratory. Candidate SNPs showed a variety of different allele frequency change patterns across generations, indicative of an underlying polygenic basis. We did not detect strong linkage around candidate SNPs, but rather a small but long-ranging effect. In conclusion, we found that history played a major role in genomic variation and evolution, with initially differentiated populations reaching the same adaptive outcome through different genetic routes.

The Genomic Architecture of Flowering Time Varies Across Space and Time in Mimulus guttatus

The degree to which genomic architecture varies across space and time is central to the evolution of genomes in response to natural selection. Bulked-segregant mapping combined with pooled sequencing provides an efficient means to estimate the effect of genetic variants on quantitative traits. We develop a novel likelihood framework to identify segregating variation within multiple populations and generations while accommodating estimation error on a sample- and SNP-specific basis. We use this method to map loci for flowering time within natural populations of Mimulus guttatus, collecting the early- and late-flowering plants from each of three neighboring populations and two consecutive generations. Structural variants, such as inversions, and genes from multiple flowering-time pathways exhibit the strongest associations with flowering time. We find appreciable variation in genetic effects on flowering time across both time and space; the greatest differences evident between populations, where numerous factors (environmental variation, genomic background, and private polymorphisms) likely contribute to heterogeneity. However, the changes across years within populations clearly identify genotype-by-environment interactions as an important influence on flowering time variation.

LONG-TERM EXPERIMENTAL EVOLUTION IN ESCHERICHIA COLI. III. VARIATION AMONG REPLICATE POPULATIONS IN CORRELATED RESPONSES TO NOVEL ENVIRONMENTS

Twelve populations of Escherichia coli were founded from a single clone and propagated for 2000 generations in identical glucose-limited environments. During this time, the mean fitnesses of the evolving populations relative to their common ancestor improved greatly, but their fitnesses relative to one another diverged only slightly. Although the populations showed similar fitness increases, they may have done so by different underlying adaptations, or they may have diverged in other respects by random genetic drift. Therefore, we examined the relative fitnesses of independently derived genotypes in two other sugars, maltose and lactose, to determine whether they were homogeneous or heterogeneous in these environments. The genetic variation among the derived lines in fitness on maltose and lactose was more than 100-times greater than their variation in fitness on glucose. Moreover, the glucose-adapted genotypes, on average, showed significant adaptation to lactose, but not to maltose. That pathways for use of maltose and glucose are virtually identical in E. coli, except for their distinct mechanisms of uptake, suggests that the derived genotypes have adapted primarily by improved glucose transport. From consideration of the number of generations of divergence, the mutation rate in E. coli, and the proportion of its genome required for growth on maltose (but not glucose), we hypothesize that pleiotropy involving the selected alleles, rather than random genetic drift of alleles at other loci, was the major cause of the variation among the derived genotypes in fitness on these other sugars.

THE GENETICS OF CENTRAL AND MARGINAL POPULATIONS OF DROSOPHILA SERRATA. I. GENETIC VARIATION FOR STRESS RESISTANCE AND SPECIES BORDERS

A selection experiment was used to determine if levels of genetic variance in an ecologically important trait, desiccation resistance, were different in central and marginal populations. Four populations of Drosophila serrata were sampled from central and marginal areas of its distribution, along a 3000-km stretch of Australia's east coast. Rainfall patterns along this stretch of coastline change from a tropical cycle in the north to a temperate cycle in the south. Replicate lines from the four populations underwent selection for desiccation resistance for 14 generations. Realized heritabilities calculated after 10 and 14 generations of selection indicated that the four populations differed significantly in the level of genetic variation for desiccation resistance available to selection. Populations from the more southern marginal areas had lower realized heritabilities than more northern central populations. However, a corresponding increase in mean desiccation resistance toward the margin was not found. A mechanism by which D. serrata seemed to have responded to selection was a reduction in the extent that metabolic rate was increased when flies were exposed to low humidity. This response indicates genetic variation for the control of metabolic rate. In contrast, increased desiccation resistance was not associated with lipid or glycogen levels. Increased resistance to desiccation was accompanied by increased starvation resistance, but radiation resistance was not affected. Selection did not affect the degree that replicate lines or populations had diverged.

Evolving the future: toward a science of intentional change

Humans possess great capacity for behavioral and cultural change, but our ability to manage change is still limited. This article has two major objectives: first, to sketch a basic science of intentional change centered on evolution; second, to provide examples of intentional behavioral and cultural change from the applied behavioral sciences, which are largely unknown to the basic sciences community. All species have evolved mechanisms of phenotypic plasticity that enable them to respond adaptively to their environments. Some mechanisms of phenotypic plasticity count as evolutionary processes in their own right. The human capacity for symbolic thought provides an inheritance system having the same kind of combinatorial diversity as does genetic recombination and antibody formation. Taking these propositions seriously allows an integration of major traditions within the basic behavioral sciences, such as behaviorism, social constructivism, social psychology, cognitive psychology, and evolutionary psychology, which are often isolated and even conceptualized as opposed to one another. The applied behavioral sciences include well-validated examples of successfully managing behavioral and cultural change at scales ranging from individuals to small groups to large populations. However, these examples are largely unknown beyond their disciplinary boundaries, for lack of a unifying theoretical framework. Viewed from an evolutionary perspective, they are examples of managing evolved mechanisms of phenotypic plasticity, including open-ended processes of variation and selection. Once the many branches of the basic and applied behavioral sciences become conceptually unified, we are closer to a science of intentional change than one might think.

The population genomics of repeated evolution in the blind cavefish Astyanax mexicanus

Distinct populations of Astyanax mexicanus cavefish offer striking examples of repeatable convergence or parallelism in their independent evolutions from surface to cave phenotypes. However, the extent to which the repeatability of evolution occurred at the genetic level remains poorly understood. To address this, we first characterized the genetic diversity of 518 single-nucleotide polymorphisms (SNPs), obtained through RAD tag sequencing and distributed throughout the genome, in seven cave and three groups of surface populations. The cave populations represented two distinct lineages (old and new). Thirty-one SNPs were significantly differentiated between surface and old cave populations, two SNPs were differentiated between surface and new cave populations, and 44 SNPs were significantly differentiated in both old and new cave populations. In addition, we determined whether these SNPs map to the same locations of previously described quantitative trait loci (QTL) between surface and cave populations. A total of 25 differentiated SNPs co-map with several QTL, such as one containing a fibroblast growth factor gene (Fgf8) involved in eye development and lens size. Further, the identity of many SNPs that co-mapped with QTL was the same in independently derived cave populations. These conclusions were further confirmed by haplotype analyses of SNPs within QTL regions. Our findings indicate that the repeatability of evolution at the genetic level is substantial, suggesting that ancestral standing genetic variation significantly contributed to the population genetic variability used in adaptation to the cave environment.

How to run far: multiple solutions and sex-specific responses to selective breeding for high voluntary activity levels

The response to uniform selection may occur in alternate ways that result in similar performance. We tested for multiple adaptive solutions during artificial selection for high voluntary wheel running in laboratory mice. At generation 43, the four replicate high runner (HR) lines averaged 2.85-fold more revolutions per day as compared with four non-selected control (C) lines, and females ran 1.11-fold more than males, with no sex-by-linetype interaction. Analysis of variance indicated significant differences among C lines but not among HR for revolutions per day. By contrast, average speed varied significantly among HR lines, but not among C, and showed a sex-by-linetype interaction, with the HR/C ratio being 2.02 for males and 2.45 for females. Time spent running varied among both HR and C lines, and showed a sex-by-linetype interaction, with the HR/C ratio being 1.52 for males but only 1.17 for females. Thus, females (speed) and males (speed, but also time) evolved differently, as did the replicate selected lines. Speed and time showed a trade-off among HR but not among C lines. These results demonstrate that uniform selection on a complex trait can cause consistent responses in the trait under direct selection while promoting divergence in the lower-level components of that trait.

Conditions for mutation-order speciation

Two models for speciation via selection have been proposed. In the well-known model of 'ecological speciation', divergent natural selection between environments drives the evolution of reproductive isolation. In a second 'mutation-order' model, different, incompatible mutations (alleles) fix in different populations adapting to the same selective pressure. How to demonstrate mutation-order speciation has been unclear, although it has been argued that it can be ruled out when gene flow occurs because the same, most advantageous allele will fix in all populations. However, quantitative examination of the interaction of factors influencing the likelihood of mutation-order speciation is lacking. We used simulation models to study how gene flow, hybrid incompatibility, selective advantage, timing of origination of new mutations and an initial period of allopatric differentiation affect population divergence via the mutation-order process. We find that at least some population divergence can occur under a reasonably wide range of conditions, even with moderate gene flow. However, strong divergence (e.g. fixation of different alleles in different populations) requires very low gene flow, and is promoted when (i) incompatible mutations have similar fitness advantages, (ii) less fit mutations arise slightly earlier in evolutionary time than more fit alternatives, and (iii) allopatric divergence occurs prior to secondary contact.

Experimental evolution with E. coli in diverse resource environments. I. Fluctuating environments promote divergence of replicate populations

BACKGROUND

Environmental conditions affect the topology of the adaptive landscape and thus the trajectories followed by evolving populations. For example, a heterogeneous environment might lead to a more rugged adaptive landscape, making it more likely that replicate populations would evolve toward distinct adaptive peaks, relative to a uniform environment. To date, the influence of environmental variability on evolutionary dynamics has received relatively little experimental study.

RESULTS

We report findings from an experiment designed to test the effects of environmental variability on the adaptation and divergence of replicate populations of E. coli. A total of 42 populations evolved for 2000 generations in 7 environmental regimes that differed in the number, identity, and presentation of the limiting resource. Regimes were organized in two sets, having the sugars glucose and maltose singly and in combination, or glucose and lactose singly and in combination. Combinations of sugars were presented either simultaneously or as temporally fluctuating resource regimes. This design allowed us to compare the effects of resource identity and presentation on the evolutionary trajectories followed by replicate populations. After 2000 generations, the fitness of all populations had increased relative to the common ancestor, but to different extents. Populations evolved in glucose improved the least, whereas populations evolving in maltose or lactose increased the most in their respective sets. Among-population divergence also differed across regimes, with variation higher in those groups that evolved in fluctuating environments than in those that faced constant resource regimens. This divergence under the fluctuating conditions increased between 1000 and 2000 generations, consistent with replicate populations evolving toward distinct adaptive peaks.

CONCLUSIONS

These results support the hypothesis that environmental heterogeneity can give rise to more rugged adaptive landscapes, which in turn promote evolutionary diversification. These results also demonstrate that this effect depends on the form of environmental heterogeneity, with greater divergence when the pairs of resources fluctuated temporally rather than being presented simultaneously.

Variation at range margins across multiple spatial scales: environmental temperature, population genetics and metabolomic phenotype

Range margins are spatially complex, with environmental, genetic and phenotypic variations occurring across a range of spatial scales. We examine variation in temperature, genes and metabolomic profiles within and between populations of the subalpine perennial plant Arabidopsis lyrata ssp. petraea from across its northwest European range. Our surveys cover a gradient of fragmentation from largely continuous populations in Iceland, through more fragmented Scandinavian populations, to increasingly widely scattered populations at the range margin in Scotland, Wales and Ireland. Temperature regimes vary substantially within some populations, but within-population variation represents a larger fraction of genetic and especially metabolomic variances. Both physical distance and temperature differences between sites are found to be associated with genetic profiles, but not metabolomic profiles, and no relationship was found between genetic and metabolomic population structures in any region. Genetic similarity between plants within populations is the highest in the fragmented populations at the range margin, but differentiation across space is the highest there as well, suggesting that regional patterns of genetic diversity may be scale dependent.

Spatio-temporal variation in the strength and mode of selection acting on major histocompatibility complex diversity in water vole (Arvicola terrestris) metapopulations

Patterns of spatio-temporal genetic variation at a class II major histocompatibility complex (MHC) locus and multiple microsatellite loci were analysed within and between three water vole metapopulations in Scotland, UK. Comparisons of MHC and microsatellite spatial genetic differentiation, based on standardised tests between two demographically asynchronous zones within a metapopulation, suggested that spatial MHC variation was affected by balancing selection, directional selection and random genetic drift, but that the relative effects of these microevolutionary forces vary temporally. At the metapopulation level, between-year differentiation for MHC loci was significantly correlated with that of microsatellites, signifying that neutral factors such as migration and drift were primarily responsible for overall temporal genetic change at the metapopulation scale. Between metapopulations, patterns of genetic differentiation implied that, at large spatial scales, MHC variation was primarily affected by directional selection and drift. Levels of MHC heterozygosity in excess of Hardy-Weinberg expectations were consistent with overdominant balancing selection operating on MHC variation within metapopulations. However, this effect was not constant among all samples, indicating temporal variation in the strength of selection relative to other factors. The results highlight the benefit of contrasting variation at MHC with neutral markers to separate the effects of stochastic and deterministic microevolutionary forces, and add to a growing body of evidence showing that the mode and relative strength of selection acting on MHC diversity varies both spatially and temporally.

Rapid divergence of genetic variance-covariance matrix within a natural population

The matrix of genetic variances and covariances (G matrix) represents the genetic architecture of multiple traits sharing developmental and genetic processes and is central for predicting phenotypic evolution. These predictions require that the G matrix be stable. Yet the timescale and conditions promoting G matrix stability in natural populations remain unclear. We studied stability of the G matrix in a 20-year evolution field experiment, where a population of the cosmopolitan parthenogenetic soil nematode Acrobeloides nanus was subjected to drift and divergent selection (benign and stress environments). Selection regime did not influence the level of absolute genetic constraints: under both regimes, two genetic dimensions for three life-history traits were identified. A substantial response to selection in principal components structure and in general matrix pattern was indicated by three statistical methods. G structure was also influenced by drift, with higher divergence under benign conditions. These results show that the G matrix might evolve rapidly in natural populations. The observed high dynamics of G structure probably represents the general feature of asexual species and limits the predictive power of G in phenotypic evolution analyses.

Evolutionary dynamics of molecular markers during local adaptation: a case study in Drosophila subobscura

BACKGROUND

Natural selection and genetic drift are major forces responsible for temporal genetic changes in populations. Furthermore, these evolutionary forces may interact with each other. Here we study the impact of an ongoing adaptive process at the molecular genetic level by analyzing the temporal genetic changes throughout 40 generations of adaptation to a common laboratory environment. Specifically, genetic variability, population differentiation and demographic structure were compared in two replicated groups of Drosophila subobscura populations recently sampled from different wild sources.

RESULTS

We found evidence for a decline in genetic variability through time, along with an increase in genetic differentiation between all populations studied. The observed decline in genetic variability was higher during the first 14 generations of laboratory adaptation. The two groups of replicated populations showed overall similarity in variability patterns. Our results also revealed changing demographic structure of the populations during laboratory evolution, with lower effective population sizes in the early phase of the adaptive process. One of the ten microsatellites analyzed showed a clearly distinct temporal pattern of allele frequency change, suggesting the occurrence of positive selection affecting the region around that particular locus.

CONCLUSION

Genetic drift was responsible for most of the divergence and loss of variability between and within replicates, with most changes occurring during the first generations of laboratory adaptation. We also found evidence suggesting a selective sweep, despite the low number of molecular markers analyzed. Overall, there was a similarity of evolutionary dynamics at the molecular level in our laboratory populations, despite distinct genetic backgrounds and some differences in phenotypic evolution.

The Hill–Robertson effect: evolutionary consequences of weak selection and linkage in finite populations

The 'Hill-Robertson (HR) effect' describes that linkage between sites under selection will reduce the overall effectiveness of selection in finite populations. Here we discuss the major concepts associated with the HR effect and present results of computer simulations focusing on the linkage effects generated by multiple sites under weak selection. Most models of linkage and selection forecast differences in effectiveness of selection between chromosomes or chromosomal regions involving a number of genes. The abundance and physical clustering of weakly selected mutations across genomes, however, justify the investigation of HR effects at a very local level and we pay particular attention to linkage effects among selected sites of the same gene. Overall, HR effects caused by weakly selected mutations predict differences in effectiveness of selection between genes that differ in exon-intron structures and across genes. Under this scenario, introns might play an advantageous role reducing intragenic HR effects. Finally, we summarize observations that are consistent with local HR effects in Drosophila, discuss potential consequences on population genetic studies and suggest future lines of research.

Do Species Converge during Adaptation? A Case Study inDrosophila

Adaptation to novel environments is a crucial theme in evolutionary biology, particularly because ex situ conservation forces populations to adapt to captivity. Here we analyze the evolution of life-history traits in two closely related species, Drosophila subobscura Collin and Drosophila madeirensis Monclus, during adaptation to the laboratory. Drosophila madeirensis, an endemic species from Madeira, is here shown to have less ability to adapt to the laboratory. Early fecundity was the only trait where this species showed a significant improvement with time. By comparison, D. subobscura improved in most traits, and its early fecundity increased faster than that of D. madeirensis. Our findings suggest that different species, even closely related ones, may adapt at different rates to the same environment.

Genetic structure among and within peripheral and central populations of three endangered floodplain violets

Understanding the partitioning of genetic variance in peripheral and central populations may shed more light on the effects of genetic drift and gene flow on population genetic structure and, thereby, improve attempts to conserve genetic diversity. We analysed genetic structure of peripheral and central populations of three insect-pollinated violets (Viola elatior, Viola pumila, Viola stagnina) to evaluate to what extent these patterns can be explained by gene flow and genetic drift. Amplified fragment length polymorphism was used to analyse 930 individuals of 50 populations. Consistent with theoretical predictions, peripheral populations were smaller and more isolated, differentiation was stronger, and genetic diversity and gene flow lower in peripheral populations of V. pumila and V. stagnina. In V. elatior, probably historic fragmentation effects linked to its specific habitat type were superimposed on the plant geographic (peripheral-central) patterns, resulting in lower relative importance of gene flow in central populations. Genetic variation between regions (3-6%), among (30-37%) and within populations (60-64%) was significant. Peripheral populations lacked markers that were rare and localized in central populations. Loss of widespread markers in peripheral V. stagnina populations indicated genetic erosion. Autocorrelation within populations was statistically significant up to a distance of 10-20 m. Higher average genetic similarity in peripheral populations than in central ones indicated higher local gene flow, probably owing to management practices. Peripheral populations contributed significantly to genetic variation and contained unique markers, which made them valuable for the conservation of genetic diversity.

Geographical variation in selection, from phenotypes to molecules

Molecular technologies now allow researchers to isolate quantitative trait loci (QTLs) and measure patterns of gene sequence variation within chromosomal regions containing important polymorphisms. I develop a simulation model to investigate gene sequence evolution within genomic regions that harbor QTLs. The QTLs influence a trait experiencing geographical variation in selection, which is common in nature and produces obvious differentiation at the phenotypic level. Counter to expectations, the simulations suggest that selection can substantially affect quantitative genetic variation without altering the amount and pattern of molecular variation at sites closely linked to the QTLs. Even with large samples of gene sequences, the likelihood of rejecting neutrality is often low. The exception is situations where strong selection is combined with low migration among demes, conditions that may be common in many plant species. The results have implications for gene sequence surveys and, perhaps more generally, for interpreting the apparently weak connection between levels of molecular and quantitative trait variation within species.

Molecular and quantitative genetic divergence among populations of house mice with known evolutionary histories

Evolutionary biologists have long been interested in the processes influencing population differentiation, but separating the effects of neutral and adaptive evolution has been an obstacle for studies of population subdivision. A recently developed method allows tests of whether disruptive (ie, spatially variable) or stabilizing (ie, spatially uniform) selection is influencing phenotypic differentiation among subpopulations. This method, referred to as the F(ST) vs Q(ST) comparison, separates the total additive genetic variance into within- and among-population components and evaluates this level of differentiation against a neutral hypothesis. Thus, levels of neutral molecular (F(ST)) and quantitative genetic (Q(ST)) divergence are compared to evaluate the effects of selection and genetic drift on phenotypic differentiation. Although the utility of such comparisons appears great, its accuracy has not yet been evaluated in populations with known evolutionary histories. In this study, F(ST) vs Q(ST) comparisons were evaluated using laboratory populations of house mice with known evolutionary histories. In this model system, the F(ST) vs Q(ST) comparisons between the selection groups should reveal quantitative trait differentiation consistent with disruptive selection, while the F(ST) vs Q(ST) comparisons among lines within the selection groups should suggest quantitative trait differentiation in agreement with drift. We find that F(ST) vs Q(ST) comparisons generally produce the correct evolutionary inference at each level in the population hierarchy. Additionally, we demonstrate that when strong selection is applied between populations Q(ST) increases relative to Q(ST) among populations diverging by drift. Finally, we show that the statistical properties of Q(ST), a variance component ratio, need further investigation.

Divergent evolution during an experimental adaptive radiation

How repeatable a process is evolution? Comparative studies of multicellular eukaryotes and experimental studies with unicellular prokaryotes document the repeated evolution of adaptive phenotypes during similar adaptive radiations, suggesting that the outcome of adaptive radiation is broadly reproducible. The goal of this study was to test this hypothesis by using phenotypic traits to infer the genetic basis of adaptation to simple carbon-limited environments in an extensive adaptive radiation. We used a clone of the bacterium Pseudomonas fluorescens to found two sets of experimental lines. The first set of lines was allowed to adapt to one of 23 novel environments for 1100 generations while the second set of lines was allowed to accumulate mutations by drift for 2000 generations. All lines were then assayed in the 95 environments provided by Biolog microplates to determine the phenotypic consequences of selection and drift. Replicate selection lines propagated in a common environment evolved similar adaptive components of their phenotype but showed extensive variation in non-adaptive phenotypic traits. This variation in non-adaptive phenotypic traits primarily resulted from the ascendance of different beneficial mutations in different lines. We argue that these results reconcile experimental and comparative approaches to studying adaptation by demonstrating that the convergent phenotypic evolution that occurs during adaptive radiation may be associated with radically different sets of beneficial mutations.

Contribution of individual random mutations to genotype-by-environment interactions in Escherichia coli

Numerous studies have shown genotype-by-environment (GxE) interactions for traits related to organismal fitness. However, the genetic architecture of the interaction is usually unknown because these studies used genotypes that differ from one another by many unknown mutations. These mutations were also present as standing variation in populations and hence had been subject to prior selection. Based on such studies, it is therefore impossible to say what fraction of new, random mutations contributes to GxE interactions. In this study, we measured the fitness in four environments of 26 genotypes of Escherichia coli, each containing a single random insertion mutation. Fitness was measured relative to their common progenitor, which had evolved on glucose at 37 degrees C for the preceding 10,000 generations. The four assay environments differed in limiting resource and temperature (glucose, 28 degrees C; maltose, 28 degrees C; glucose, 37 degrees C; and maltose, 37 degrees C). A highly significant interaction between mutation and resource was found. In contrast, there was no interaction involving temperature. The resource interaction reflected much higher among mutation variation for fitness in maltose than in glucose. At least 11 mutations (42%) contributed to this GxE interaction through their differential fitness effects across resources. Beneficial mutations are generally thought to be rare but, surprisingly, at least three mutations (12%) significantly improved fitness in maltose, a resource novel to the progenitor. More generally, our findings demonstrate that GxE interactions can be quite common, even for genotypes that differ by only one mutation and in environments differing by only a single factor.

Variability in phenotypic traits in core and peripheral populations of wild barley Hordeum spontaneum Koch

Populations of wild barley, H. spontaneum Koch., were collected in two countries, Israel and Turkmenistan, in environments representing two similar sharp clines of aridity. This allowed us to use the same criteria to define species core and periphery in two regions. Plants from 10 Israeli and 19 Turkmenian populations were grown in a field trial with three water treatments and compared for amount and structure of variation in phenological and morphological traits. Extent of variation was similar in populations at species border (periphery by aridity criterion) or at species border and near it (by habitat), and in populations inhabiting favorable environments away from the border (core). In contrast, two regions (Israel and Turkmenistan) exhibited different amount of variation in phenotypic traits. Israeli populations were more diverse than Turkmenian populations in all except one phenotypic traits and the variation was differently structured among regionally specific core and periphery. While Turkmenian core and periphery defined by either criterion did not differ for all except one trait, corresponding Israeli populations exhibited opposite patterns of variation for trait complexes. Israeli core (Mediterranean) populations were twice as variable than peripheral (desert) populations in parameters of reproductive growth (stem, spike and awn length) and grain filling (spikelet weight) and half as variable in the length of flag and penultimate leaves and onset of reproduction. Possible modes of regionally specific natural selection as a cause of regional/local variation in phenotypic triats are discussed.

Perspective: reverse evolution

For some time, the reversibility of evolution was primarily discussed in terms of comparative patterns. Only recently has this problem been studied using experimental evolution over shorter evolutionary time frames. This has raised questions of definition, experimental procedure, and the hypotheses being tested. Experimental evolution has provided evidence for multiple population genetic mechanisms in reverse evolution, including pleiotropy and mutation accumulation. It has also pointed to genetic factors that might prevent reverse evolution, such as a lack of genetic variability, epistasis, and differential genotype-by-environment interactions. The main focus of this perspective is on laboratory studies and their relevance to the genetics of reverse evolution. We discuss reverse evolution experiments with Drosophila, bacterial, and viral populations. Field studies of the reverse evolution of melanism in the peppered moth are also reviewed.

Nuclear and cytoplasmic contributions to intraspecific divergence in an annual legume

The genetic architecture of trait differentiation was evaluated between two ecologically distinct populations of Chamaecrista fasciculata. Individuals from Maryland and Illinois populations were crossed to create 10 types of seed: Maryland and Illinois parents, reciprocal F1 and F2 hybrids, and backcrosses to Maryland and to Illinois on reciprocal F1 hybrids. Reciprocal crosses created hybrid generation seeds with both Maryland and Illinois cytoplasmic backgrounds. Experimental individuals were grown in a common garden near the site of the Maryland population. In the garden, plants from the Illinois population flowered, set fruit, and died earlier than those from Maryland, likely reflecting adaptations to differences in growing season length between the two populations. Although reproductive components at the flower and whole plant level differed between the two populations, reproductive output as measured by fruit and seed production was similar. Cytoplasmic genes had a subtle but pervasive effect on population differentiation; hybrids with Maryland cytoplasm were significantly differentiated from those with Illinois cytoplasm when all characters were evaluated jointly. The nuclear genetic architecture of population differentiation was evaluated with joint scaling tests. Depending on the trait, both additive and nonadditive genetic effects contributed to population differentiation. Intraspecific genetic differentiation in this wild plant species appears to reflect a complex genetic architecture that includes the contribution of additive, dominance, and epistatic components in addition to subtle cytoplasmic effects.

Artificial selection of microbial ecosystems for 3-chloroaniline biodegradation

We present a method for selecting entire microbial ecosystems for bioremediation and other practical purposes. A population of ecosystems is established in the laboratory, each ecosystem is measured for a desired property (in our case, degradation of the environmental pollutant 3-chloroaniline), and the best ecosystems are used as 'parents' to inoculate a new generation of 'offspring' ecosystems. Over many generations of variation and selection, the ecosystems become increasingly well adapted to produce the desired property. The procedure is similar to standard artificial selection experiments except that whole ecosystems, rather than single individuals, are the units of selection. The procedure can also be understood in terms of complex system theory as a way of searching a vast combinatorial space (many thousands of microbial species and many thousands of genes within species) for combinations that are especially good at producing the desired property. Ecosystem-level selection can be performed without any specific knowledge of the species that comprise the ecosystems and can select ensembles of species that would be difficult to discover with more reductionistic methods. Once a 'designer ecosystem' has been created by ecosystem-level selection, reductionistic methods can be used to identify the component species and to discover how they interact to produce the desired effect.

The changes in genetic and environmental variance with inbreeding in Drosophila melanogaster

We performed a large-scale experiment on the effects of inbreeding and population bottlenecks on the additive genetic and environmental variance for morphological traits in Drosophila melanogaster. Fifty-two inbred lines were created from the progeny of single pairs, and 90 parent-offspring families on average were measured in each of these lines for six wing size and shape traits, as well as 1945 families from the outbred population from which the lines were derived. The amount of additive genetic variance has been observed to increase after such population bottlenecks in other studies; in contrast here the mean change in additive genetic variance was in very good agreement with classical additive theory, decreasing proportionally to the inbreeding coefficient of the lines. The residual, probably environmental, variance increased on average after inbreeding. Both components of variance were highly variable among inbred lines, with increases and decreases recorded for both. The variance among lines in the residual variance provides some evidence for a genetic basis of developmental stability. Changes in the phenotypic variance of these traits are largely due to changes in the genetic variance.