Rapid,nonparallel genomic evolution of Brassica rapa (field mustard) under experimental drought

While we know that climate change can potentially cause rapid phenotypic evolution, our understanding of the genetic basis and degree of genetic parallelism of rapid evolutionary responses to climate change is limited. In this study, we combined the resurrection approach with an evolve-and-resequence design to examine genome-wide evolutionary changes following drought. We exposed genetically similar replicate populations of the annual plant Brassica rapa derived from a field population in southern California to four generations of experimental drought or watered conditions in a greenhouse. Genome-wide sequencing of ancestral and descendant population pools identified hundreds of SNPs that showed evidence of rapidly evolving in response to drought. Several of these were in stress response genes, and two were identified in a prior study of drought response in this species. However, almost all genetic changes were unique among experimental populations, indicating that the evolutionary changes were largely nonparallel, despite the fact that genetically similar replicates of the same founder population had experienced controlled and consistent selection regimes. This nonparallelism of evolution at the genetic level is potentially because of polygenetic adaptation allowing for multiple different genetic routes to similar phenotypic outcomes. Our findings help to elucidate the relationship between rapid phenotypic and genomic evolution and shed light on the degree of parallelism and predictability of genomic evolution to environmental change.     

The probability of genetic parallelism and convergence in natural populations

Genomic and genetic methods allow investigation of how frequently the same genes are used by different populations during adaptive evolution, yielding insights into the predictability of evolution at the genetic level. We estimated the probability of gene reuse in parallel and convergent phenotypic evolution in nature using data from published studies. The estimates are surprisingly high, with mean probabilities of 0.32 for genetic mapping studies and 0.55 for candidate gene studies. The probability declines with increasing age of the common ancestor of compared taxa, from about 0.8 for young nodes to 0.1–0.4 for the oldest nodes in our study. Probability of gene reuse is higher when populations begin from the same ancestor (genetic parallelism) than when they begin from divergent ancestors (genetic convergence). Our estimates are broadly consistent with genomic estimates of gene reuse during repeated adaptation to similar environments, but most genomic studies lack data on phenotypic traits affected. Frequent reuse of the same genes during repeated phenotypic evolution suggests that strong biases and constraints affect adaptive evolution, resulting in changes at a relatively small subset of available genes. Declines in the probability of gene reuse with increasing age suggest that these biases diverge with time.     

HOMOGENEITY OF THE GENETIC VARIANCE-COVARIANCE MATRIX FOR ANTIPREDATOR TRAITS IN TWO NATURAL POPULATIONS OF THE GARTER SNAKE THAMNOPHIS ORDINOIDES

Quantitative genetic models of evolution rely on the genetic variance-covariance matrix to predict the phenotypic response to selection. Both prospective and retrospective studies of phenotypic evolution across generations rely on assumptions about the constancy of patterns of genetic covariance through time. In the absence of robust theoretical predictions about the stability of genetic covariances, this assumption must be tested with empirical comparisons of genetic parameters among populations and species. Genetic variance-covariance matrices were estimated for a suite of antipredator traits in two populations of the northwestern garter snake, Thamnophis ordinoides. The characters studied include color pattern and antipredator behaviors that interact to facilitate escape from predators. Significant heritabilities for all traits were detected in both populations. Genetic correlations and covariances were found among behaviors in both populations and between color pattern and behavior in one of the populations. Phenotypic means differed among populations, but pairwise comparisons revealed no heterogeneity of genetic parameters between the populations. The structure of the genetic variance-covariance matrix has apparently not changed significantly during the divergence of these two populations.     

DOES PHENOTYPIC PLASTICITY FOR ADULT SIZE VERSUS FOOD LEVEL IN DROSOPHILA MELANOGASTER EVOLVE IN RESPONSE TO ADAPTATION TO DIFFERENT REARING DENSITIES?

Recent studies with Drosophila have suggested that there is extensive genetic variability for phenotypic plasticity of body size versus food level. If true, we expect that the outcome of evolution at very different food levels should yield genotypes whose adult size show different patterns of phenotypic plasticity. We have tested this prediction with six independent populations of Drosophila melanogaster kept at extreme densities for 125 generations. We found that the phenotypic plasticity of body size versus food level is not affected by selection or the presence of competitors of a different genotype. However, we document increasing among population variation in phenotypic plasticity due to random genetic drift. Several reasons are explored to explain these results including the possibility that the use of highly inbred lines to make inferences about the evolution of genetically variable populations may be misleading.     

From nature to the laboratory: the impact of founder effects on adaptation

Most founding events entail a reduction in population size, which in turn leads to genetic drift effects that can deplete alleles. Besides reducing neutral genetic variability, founder effects can in principle shift additive genetic variance for phenotypes that underlie fitness. This could then lead to different rates of adaptation among populations that have undergone a population size bottleneck as well as an environmental change, even when these populations have a common evolutionary history. Thus, theory suggests that there should be an association between observable genetic variability for both neutral markers and phenotypes related to fitness. Here, we test this scenario by monitoring the early evolutionary dynamics of six laboratory foundations derived from founders taken from the same source natural population of Drosophila subobscura. Each foundation was in turn three‐fold replicated. During their first few generations, these six foundations showed an abrupt increase in their genetic differentiation, within and between foundations. The eighteen populations that were monitored also differed in their patterns of phenotypic adaptation according to their immediately ancestral founding sample. Differences in early genetic variability and in effective population size were found to predict differences in the rate of adaptation during the first 21 generations of laboratory evolution. We show that evolution in a novel environment is strongly contingent not only on the initial composition of a newly founded population but also on the stochastic changes that occur during the first generations of colonization. Such effects make laboratory populations poor guides to the evolutionary genetic properties of their ancestral wild populations.     

Does plasticity enhance or dampen phenotypic parallelism? A test with three lake–stream stickleback pairs

Parallel (and convergent) phenotypic variation is most often studied in the wild, where it is difficult to disentangle genetic vs. environmentally induced effects. As a result, the potential contributions of phenotypic plasticity to parallelism (and nonparallelism) are rarely evaluated in a formal sense. Phenotypic parallelism could be enhanced by plasticity that causes stronger parallelism across populations in the wild than would be expected from genetic differences alone. Phenotypic parallelism could be dampened if site‐specific plasticity induced differences between otherwise genetically parallel populations. We used a common‐garden study of three independent lake–stream stickleback population pairs to evaluate the extent to which adaptive divergence has a genetic or plastic basis, and to investigate the enhancing vs. dampening effects of plasticity on phenotypic parallelism. We found that lake–stream differences in most traits had a genetic basis, but that several traits also showed contributions from plasticity. Moreover, plasticity was much more prevalent in one watershed than in the other two. In most cases, plasticity enhanced phenotypic parallelism, whereas in a few cases, plasticity had a dampening effect. Genetic and plastic contributions to divergence seem to play a complimentary, likely adaptive, role in phenotypic parallelism of lake–stream stickleback. These findings highlight the value of formally comparing wild‐caught and laboratory‐reared individuals in the study of phenotypic parallelism.     

Locally adapted traits maintained in the face of high gene flow

Gene flow between phenotypically divergent populations can disrupt local adaptation or, alternatively, may stimulate adaptive evolution by increasing genetic variation. We capitalised on historical Trinidadian guppy transplant experiments to test the phenotypic effects of increased gene flow caused by replicated introductions of adaptively divergent guppies, which were translocated from high‐ to low‐predation environments. We sampled two native populations prior to the onset of gene flow, six historic introduction sites, introduction sources and multiple downstream points in each basin. Extensive gene flow from introductions occurred in all streams, yet adaptive phenotypic divergence across a gradient in predation level was maintained. Descendants of guppies from a high‐predation source site showed high phenotypic similarity with native low‐predation guppies in as few as ~12 generations after gene flow, likely through a combination of adaptive evolution and phenotypic plasticity. Our results demonstrate that locally adapted phenotypes can be maintained despite extensive gene flow from divergent populations.     

Studying phenotypic evolution using multivariate quantitative genetics

Quantitative genetics provides a powerful framework for studying phenotypic evolution and the evolution of adaptive genetic variation. Central to the approach is G, the matrix of additive genetic variances and covariances. G summarizes the genetic basis of the traits and can be used to predict the phenotypic response to multivariate selection or to drift. Recent analytical and computational advances have improved both the power and the accessibility of the necessary multivariate statistics. It is now possible to study the relationships between G and other evolutionary parameters, such as those describing the mutational input, the shape and orientation of the adaptive landscape, and the phenotypic divergence among populations. At the same time, we are moving towards a greater understanding of how the genetic variation summarized by G evolves. Computer simulations of the evolution of G, innovations in matrix comparison methods, and rapid development of powerful molecular genetic tools have all opened the way for dissecting the interaction between allelic variation and evolutionary process. Here I discuss some current uses of G, problems with the application of these approaches, and identify avenues for future research.     

Epigenetics in natural animal populations

Phenotypic plasticity is an important mechanism for populations to buffer themselves from environmental change. While it has long been appreciated that natural populations possess genetic variation in the extent of plasticity, a surge of recent evidence suggests that epigenetic variation could also play an important role in shaping phenotypic responses. Compared with genetic variation, epigenetic variation is more likely to have higher spontaneous rates of mutation and a more sensitive reaction to environmental inputs. In our review, we first provide an overview of recent studies on epigenetically encoded thermal plasticity in animals to illustrate environmentally‐mediated epigenetic effects within and across generations. Second, we discuss the role of epigenetic effects during adaptation by exploring population epigenetics in natural animal populations. Finally, we evaluate the evolutionary potential of epigenetic variation depending on its autonomy from genetic variation and its transgenerational stability. Although many of the causal links between epigenetic variation and phenotypic plasticity remain elusive, new data has explored the role of epigenetic variation in facilitating evolution in natural populations. This recent progress in ecological epigenetics will be helpful for generating predictive models of the capacity of organisms to adapt to changing climates.     

Molecular Parallelism Underlies Convergent Highland Adaptation of Maize Landraces

Convergent phenotypic evolution provides some of the strongest evidence for adaptation. However, the extent to which recurrent phenotypic adaptation has arisen via parallelism at the molecular level remains unresolved, as does the evolutionary origin of alleles underlying such adaptation. Here, we investigate genetic mechanisms of convergent highland adaptation in maize landrace populations and evaluate the genetic sources of recurrently selected alleles. Population branch excess statistics reveal substantial evidence of parallel adaptation at the level of individual single-nucleotide polymorphism (SNPs), genes, and pathways in four independent highland maize populations. The majority of convergently selected SNPs originated via migration from a single population, most likely in the Mesoamerican highlands, while standing variation introduced by ancient gene flow was also a contributor. Polygenic adaptation analyses of quantitative traits reveal that alleles affecting flowering time are significantly associated with elevation, indicating the flowering time pathway was targeted by highland adaptation. In addition, repeatedly selected genes were significantly enriched in the flowering time pathway, indicating their significance in adapting to highland conditions. Overall, our study system represents a promising model to study convergent evolution in plants with potential applications to crop adaptation across environmental gradients.     

Resolving the paradox of stasis: models with stabilizing selection explain evolutionary divergence on all timescales

We tested the ability of six quantitative genetic models to explain the evolution of phenotypic means using an extensive database compiled by Gingerich. Our approach differs from past efforts in that we use explicit models of evolutionary process, with parameters estimated from contemporary populations, to analyze a large sample of divergence data on many different timescales. We show that one quantitative genetic model yields a good fit to data on phenotypic divergence across timescales ranging from a few generations to 10 million generations. The key feature of this model is a fitness optimum that moves within fixed limits. Conversely, a model of neutral evolution, models with a stationary optimum that undergoes Brownian or white noise motion, a model with a moving optimum, and a peak shift model all fail to account for the data on most or all timescales. We discuss our results within the framework of Simpson's concept of adaptive landscapes and zones. Our analysis suggests that the underlying process causing phenotypic stasis is adaptation to an optimum that moves within an adaptive zone with stable boundaries. We discuss the implication of our results for comparative studies and phylogeny inference based on phenotypic characters.     

The repeatability of adaptive radiation during long-term experimental evolution of Escherichia coli in a multiple nutrient environment

Adaptive radiations occur when a species diversifies into different ecological specialists due to competition for resources and trade-offs associated with the specialization. The evolutionary outcome of an instance of adaptive radiation cannot generally be predicted because chance (stochastic events) and necessity (deterministic events) contribute to the evolution of diversity. With increasing contributions of chance, the degree of parallelism among different instances of adaptive radiations and the predictability of an outcome will decrease. To assess the relative contributions of chance and necessity during adaptive radiation, we performed a selection experiment by evolving twelve independent microcosms of Escherichia coli for 1000 generations in an environment that contained two distinct resources. Specialization to either of these resources involves strong trade-offs in the ability to use the other resource. After selection, we measured three phenotypic traits: 1) fitness, 2) mean colony size, and 3) colony size diversity. We used fitness relative to the ancestor as a measure of adaptation to the selective environment; changes in colony size as a measure of the evolution of new resource specialists because colony size has been shown to correlate with resource specialization; and colony size diversity as a measure of the evolved ecological diversity. Resource competition led to the rapid evolution of phenotypic diversity within microcosms. Measurements of fitness, colony size, and colony size diversity within and among microcosms showed that the repeatability of adaptive radiation was high, despite the evolution of genetic variation within microcosms. Consistent with the observation of parallel evolution, we show that the relative contributions of chance are far smaller and less important than effects due to adaptation for the traits investigated. The two-resource environment imposed similar selection pressures in independent populations and promoted parallel phenotypic adaptive radiations in all independently evolved microcosms.     

Insight into the genetic bases of climatic adaptation in experimentally evolving wheat populations

Experimental populations evolving under natural selection represent an interesting tool to study genetic bases of adaptation. Evolution of genes possibly involved in adaptive response can be followed together with the corresponding phenotypic traits. Using experimental populations of hexaploid wheat, we studied the evolution of flowering time, a major adaptive trait that synchronizes the initiation of reproduction and the occurrence of favourable environmental conditions. During 12 generations, three populations were grown in contrasted environments (Vervins North France, Le Moulon near Paris, Toulouse South France) under the influence of natural selection, drift, mutation and recombination. Evolution of diversity at the major gene VRN-1 involved in wheat vernalization response has been analysed jointly with earliness estimated in controlled conditions. Whatever the population, rapid phenotypic changes as well as parallel genotypic variations were observed in the first seven generations, probably as the result of selection acting on this major gene which explains 80% of the trait variation overall. Different allelic combinations at physically unlinked copies of VRN-1 located on distinct genomes (A, B and D) were selected between populations. As theoretically expected, due to population differentiation, a high level of genetic diversity was maintained overall in generation 12. Surprisingly, in two populations out of three, the emergence of new alleles by mutation or migration, coupled with temporal variable selection or frequency-dependent selection, allowed to maintain within-population diversity despite local genetic drift and natural selection. This result may plead for an evolutionary approach of wheat genetic resource conservation.     

Rapid evolution of larval development time during mass-rearing in the melon fly, Bactrocera cucurbitae

A quantitative genetic analysis of rapid evolution of a life history trait has been conducted on the first 24 generations of mass-rearing in the melon fly Bactrocera cucurbitae (Coquillett) (Diptera: Tephritidae). The phenotypic variance of larval development time in each generation was divided into genetic and residual components. Mean and phenotypic coefficients of variation of larval development time decreased gradually as generations proceeded as a result of artificial selection for shorter larval period in the mass-rearing procedure. There was a trend that additive genetic coefficients of variation in larval development time decreased with generations. These changes are entirely attributed to genetic responses to laboratory selection under the mass-rearing environment because the population was maintained at a very large size so as to exclude random genetic drift and inbreeding depression, which would be other factors responsible for the observed genetic changes. The residual coefficients of variation in larval development time did not change with generations. Realized heritability of larval development time was low. The heritabilities for larval development time estimated from parent–offspring regression at generations 60 and 70, when the evolutionary plateau was asymptotically reached, were not significantly larger than 0. Received: April 22, 1999 / Accepted: September 20, 1999     

Altered regulation of the OmpF porin by Fis in Escherichia coli during an evolution experiment and between B and K-12 strains

The phenotypic plasticity of global regulatory networks provides bacteria with rapid acclimation to a wide range of environmental conditions, while genetic changes in those networks provide additional flexibility as bacteria evolve across long time scales. We previously identified mutations in the global regulator-encoding gene fis that enhanced organismal fitness during a long-term evolution experiment with Escherichia coli. To gain insight into the effects of these mutations, we produced two-dimensional protein gels with strains carrying different fis alleles, including a beneficial evolved allele and one with an in-frame deletion. We found that Fis controls the expression of the major porin-encoding gene ompF in the E. coli B-derived ancestral strain used in the evolution experiment, a relationship that has not been described before. We further showed that this regulatory connection evolved over two different time scales, perhaps explaining why it was not observed before. On the longer time scale, we showed that this regulation of ompF by Fis is absent from the more widely studied K-12 strain and thus is specific to the B strain. On a shorter time scale, this regulatory linkage was lost during 20,000 generations of experimental evolution of the B strain. Finally, we mapped the Fis binding sites in the ompF regulatory region, and we present a hypothetical model of ompF expression that includes its other known regulators.     

Reverse evolution of fitness in Drosophila melanogaster

Abstract The evolution of fitness is central to evolutionary theory, yet few experimental systems allow us to track its evolution in genetically and environmentally relevant contexts. Reverse evolution experiments allow the study of the evolutionary return to ancestral phenotypic states, including fitness. This in turn permits well‐defined tests for the dependence of adaptation on evolutionary history and environmental conditions. In the experiments described here, 20 populations of heterogeneous evolutionary histories were returned to their common ancestral environment for 50 generations, and were then compared with both their immediate differentiated ancestors and populations which had remained in the ancestral environment. One measure of fitness returned to ancestral levels to a greater extent than other characters did. The phenotypic effects of reverse evolution were also contingent on previous selective history. Moreover, convergence to the ancestral state was highly sensitive to environmental conditions. The phenotypic plasticity of fecundity, a character directly selected for, evolved during the experimental time frame. Reverse evolution appears to force multiple, diverged populations to converge on a common fitness state through different life‐history and genetic changes.