Spatial environmental variation can select for evolvability

Previous studies have shown that temporally fluctuating environments can create indirect selection for modifiers of evolvability. Here, we use a simple computational model to investigate whether spatially varying environments (multiple demes with limited migration among them, and a different, static selective optimum in each) can also create indirect selection for increased evolvability. The answer is surprisingly complicated. Spatial variation in the environment can sharply reduce the survival rate of migrants, because migrants may be maladapted to their new deme, relative to incumbents. The incumbent advantage can be removed by occasional extinctions in single demes. After all incumbents in a particular deme die, incoming migrants from other demes will, on average, be similarly maladapted to the new environment. This sets off a race to adapt rapidly. Over many extinction events, and the subsequent invasions by maladapted immigrants into a new environment, indirect selection for the ability to adapt rapidly, also known as high evolvability, may result.     

Maintaining heritable variation via sex-limited temporally fluctuating selection: a phenotypic model accommodating non-Mendelian epigenetic effects

Using a phenotypic model, we show that significant heritable variation can be maintained in a population subjected to temporally fluctuating selection if only one sex is subject to selection. In fact, more variation is maintained with sex-limited selection at a given selection intensity than if both sexes are subject to half that selection intensity. This result is commensurate with existing population genetic models. However, genetic models may be inappropriate for sexually selected traits because many of them may be of non-genetic origin, such as maternal effects or – more likely –epigenetic effects. Phenotypic models obviate this problem by accommodating both genetic and epigenetic effects, as well as maternaleffects. Our phenotypic model of sex-limited temporally fluctuating selection shows that substantial heritable variation can be maintained and therebyprovides impetus to develop population epigenetic models.     

Fluctuating Environments Maintain Genetic Diversity through Neutral Fitness Effects and Balancing Selection

Genetic variation is the raw material upon which selection acts. The majority of environmental conditions change over time and therefore may result in variable selective effects. How temporally fluctuating environments impact the distribution of fitness effects and in turn population diversity is an unresolved question in evolutionary biology. Here, we employed continuous culturing using chemostats to establish environments that switch periodically between different nutrient limitations and compared the dynamics of selection to static conditions. We used the pooled Saccharomyces cerevisiae haploid gene deletion collection as a synthetic model for populations comprising thousands of unique genotypes. Using barcode sequencing, we find that static environments are uniquely characterized by a small number of high-fitness genotypes that rapidly dominate the population leading to dramatic decreases in genetic diversity. By contrast, fluctuating environments are enriched in genotypes with neutral fitness effects and an absence of extreme fitness genotypes contributing to the maintenance of genetic diversity. We also identified a unique class of genotypes whose frequencies oscillate sinusoidally with a period matching the environmental fluctuation. Oscillatory behavior corresponds to large differences in short-term fitness that are not observed across long timescales pointing to the importance of balancing selection in maintaining genetic diversity in fluctuating environments. Our results are consistent with a high degree of environmental specificity in the distribution of fitness effects and the combined effects of reduced and balancing selection in maintaining genetic diversity in the presence of variable selection.     

The role of migration in the evolution of phenotypic switching

Stochastic switching is an example of phenotypic bet hedging, where an individual can switch between different phenotypic states in a fluctuating environment. Although the evolution of stochastic switching has been studied when the environment varies temporally, there has been little theoretical work on the evolution of phenotypic switching under both spatially and temporally fluctuating selection pressures. Here, we explore the interaction of temporal and spatial change in determining the evolutionary dynamics of phenotypic switching. We find that spatial variation in selection is important; when selection pressures are similar across space, migration can decrease the rate of switching, but when selection pressures differ spatially, increasing migration between demes can facilitate the evolution of higher rates of switching. These results may help explain the diverse array of non-genetic contributions to phenotypic variability and phenotypic inheritance observed in both wild and experimental populations.     

Gene genealogies in a metapopulation

A simple genealogical process is found for samples from a metapopulation, which is a population that is subdivided into a large number of demes, each of which is subject to extinction and recolonization and receives migrants from other demes. As in the migration-only models studied previously, the genealogy of any sample includes two phases: a brief sample-size adjustment followed by a coalescent process that dominates the history. This result will hold for metapopulations that are composed of a large number of demes. It is robust to the details of population structure, as long as the number of possible source demes of migrants and colonists for each deme is large. Analytic predictions about levels of genetic variation are possible, and results for average numbers of pairwise differences within and between demes are given. Further analysis of the expected number of segregating sites in a sample from a single deme illustrates some previously known differences between migration and extinction/recolonization. The ancestral process is also amenable to computer simulation. Simulation results show that migration and extinction/recolonization have very different effects on the site-frequency distribution in a sample from a single deme. Migration can cause a U-shaped site-frequency distribution, which is qualitatively similar to the pattern reported recently for positive selection. Extinction and recolonization, in contrast, can produce a mode in the site-frequency distribution at intermediate frequencies, even in a sample from a single deme.     

Fixation probability for a beneficial allele and a mutant strategy in a linear game under weak selection in a finite island model

The effect of population structure on the probability of fixation of a newly introduced mutant under weak selection is studied using a coalescent approach. Wright's island model in a framework of a finite number of demes is assumed and two selection regimes are considered: a beneficial allele model and a linear game among offspring. A first-order approximation of the fixation probability for a single mutant with respect to the intensity of selection is deduced. The approximation requires the calculation of expected coalescence times, under neutrality, for lineages starting from two or three sampled individuals. The results are obtained in a general setting without assumptions on the number of demes, the deme size or the migration rate, which allows for simultaneous coalescence or migration events in the genealogy of the sampled individuals. Comparisons are made with limit cases as the deme size or the number of demes goes to infinity or the migration rate goes to zero for which a diffusion approximation approach is possible. Conditions for selection to favor a mutant strategy replacing a resident strategy in the context of a linear game in a finite island population are addressed.     

Evolution of altruism under group selection in large and small populations in fluctuating environments

A continuous, graded form of group selection which does not involve extinction of demes can effectively oppose selection on the individual level against an altruistic allele under fluctuating environments in infinitely large demes among which uniform mixing occurs every generation. Although group selection cannot alter the conditions necessary for the initial increase of altruistic alleles, group selection can significantly influence the stationary distribution of gene frequency which is attained once stochastic forces have allowed theirintroduction. Drift is a more effective source of variation than fluctuations in selection when the variance in selection is moderate to small. High numbers of demes promote polymorphism under both graded group selection and extinction group selection.     

Polymorphism and divergence for island-model species

Estimates of the scaled selection coefficient, gamma of Sawyer and Hartl, are shown to be remarkably robust to population subdivision. Estimates of mutation parameters and divergence times, in contrast, are very sensitive to subdivision. These results follow from an analysis of natural selection and genetic drift in the island model of subdivision in the limit of a very large number of subpopulations, or demes. In particular, a diffusion process is shown to hold for the average allele frequency among demes in which the level of subdivision sets the timescale of drift and selection and determines the dynamic equilibrium of allele frequencies among demes. This provides a framework for inference about mutation, selection, divergence, and migration when data are available from a number of unlinked nucleotide sites. The effects of subdivision on parameter estimates depend on the distribution of samples among demes. If samples are taken singly from different demes, the only effect of subdivision is in the rescaling of mutation and divergence-time parameters. If multiple samples are taken from one or more demes, high levels of within-deme relatedness lead to low levels of intraspecies polymorphism and increase the number of fixed differences between samples from two species. If subdivision is ignored, mutation parameters are underestimated and the species divergence time is overestimated, sometimes quite drastically. Estimates of the strength of selection are much less strongly affected and always in a conservative direction.     

EVOLUTION IN FLUCTUATING ENVIRONMENTS: DECOMPOSING SELECTION INTO ADDITIVE COMPONENTS OF THE ROBERTSON–PRICE EQUATION

We analyze the stochastic components of the Robertson–Price equation for the evolution of quantitative characters that enables decomposition of the selection differential into components due to demographic and environmental stochasticity. We show how these two types of stochasticity affect the evolution of multivariate quantitative characters by defining demographic and environmental variances as components of individual fitness. The exact covariance formula for selection is decomposed into three components, the deterministic mean value, as well as stochastic demographic and environmental components. We show that demographic and environmental stochasticity generate random genetic drift and fluctuating selection, respectively. This provides a common theoretical framework for linking ecological and evolutionary processes. Demographic stochasticity can cause random variation in selection differentials independent of fluctuating selection caused by environmental variation. We use this model of selection to illustrate that the effect on the expected selection differential of random variation in individual fitness is dependent on population size, and that the strength of fluctuating selection is affected by how environmental variation affects the covariance in Malthusian fitness between individuals with different phenotypes. Thus, our approach enables us to partition out the effects of fluctuating selection from the effects of selection due to random variation in individual fitness caused by demographic stochasticity.     

Recombination Modification in a Fluctuating Environment

This paper examines the theory of the evolution of increased recombination between two loci subjected to interactive selection in a temporally fluctuating environment. Both cyclical and stochastic environments are considered. It is shown that temporal variation in the linkage disequilibrium coefficient for the pair of selected loci, due to fluctuations in the selective values of the genotypes at these loci, can give rise to selection in favor of modifier genes increasing recombination. The equilibrium level of recombination established in a given population depends on several factors; it is highest for intermediate values of the environmental periodicity or autocorrelation, for cases when the modifier genes are themselves linked to the selected loci, and for high levels of environmental variation. In general, it seems that the rate of modification of recombination values by this process will be low except when the modifiers are tightly linked to the selected loci. The possible evolutionary significance of this process is discussed in relation to observations on genetic systems of plants and animals.     

The Effects of Background and Interference Selection on Patterns of Genetic Variation in Subdivided Populations

It is well known that most new mutations that affect fitness exert deleterious effects and that natural populations are often composed of subpopulations (demes) connected by gene flow. To gain a better understanding of the joint effects of purifying selection and population structure, we focus on a scenario where an ancestral population splits into multiple demes and study neutral diversity patterns in regions linked to selected sites. In the background selection regime of strong selection, we first derive analytic equations for pairwise coalescent times and F_ST as a function of time after the ancestral population splits into two demes and then construct a flexible coalescent simulator that can generate samples under complex models such as those involving multiple demes or nonconservative migration. We have carried out extensive forward simulations to show that the new methods can accurately predict diversity patterns both in the nonequilibrium phase following the split of the ancestral population and in the equilibrium between mutation, migration, drift, and selection. In the interference selection regime of many tightly linked selected sites, forward simulations provide evidence that neutral diversity patterns obtained from both the nonequilibrium and equilibrium phases may be virtually indistinguishable for models that have identical variance in fitness, but are nonetheless different with respect to the number of selected sites and the strength of purifying selection. This equivalence in neutral diversity patterns suggests that data collected from subdivided populations may have limited power for differentiating among the selective pressures to which closely linked selected sites are subject.     

Selection in a fluctuating environment leads to decreased genetic variation and facilitates the evolution of phenotypic plasticity

Changes in the environment are expected to induce changes in the quantitative genetic variation, which influences the ability of a population to adapt to environmental change. Furthermore, environmental changes are not constant in time, but fluctuate. Here, we investigate the effect of rapid, continuous and/or fluctuating temperature changes in the seed beetle Callosobruchus maculatus, using an evolution experiment followed by a split-brood experiment. In line with expectations, individuals responded in a plastic way and had an overall higher potential to respond to selection after a rapid change in the environment. After selection in an environment with increasing temperature, plasticity remained unchanged (or decreased) and environmental variation decreased, especially when fluctuations were added; these results were unexpected. As expected, the genetic variation decreased after fluctuating selection. Our results suggest that fluctuations in the environment have major impact on the response of a population to environmental change; in a highly variable environment with low predictability, a plastic response might not be beneficial and the response is genetically and environmentally canalized resulting in a low potential to respond to selection and low environmental sensitivity. Interestingly, we found greater variation for phenotypic plasticity after selection, suggesting that the potential for plasticity to evolve is facilitated after exposure to environmental fluctuations. Our study highlights that environmental fluctuations should be considered when investigating the response of a population to environmental change.     

The many-demes limit for selection and drift in a subdivided population

A diffusion approximation is obtained for the frequency of a selected allele in a population comprised of many subpopulations or demes. The form of the diffusion is equivalent to that for an unstructured population, except that it occurs on a longer time scale when migration among demes is restricted. This many-demes diffusion limit relies on the collection of demes always being in statistical equilibrium with respect to migration and drift for a given allele frequency in the total population. Selection is assumed to be weak, in inverse proportion to the number of demes, and the results hold for any deme sizes and migration rates greater than zero. The distribution of allele frequencies among demes is also described.     

Adaptive topography of fluctuating selection in a Mendelian population

An adaptive topography is derived for a large randomly mating diploid population under weak density-independent selection in a fluctuating environment. Assuming a stationary distribution of environmental states with no temporal autocorrelation, a diffusion approximation for population size and allele frequency, p, reveals that the expected change in p involves the gradient with respect to p of the stochastic intrinsic rate of increase (the density-independent long-run growth rate), r = r - sigma 2 e/2, where r is the mean Malthusian fitness in the average environment and is the environmental variance in population growth rate. The expected relative fitness of a genotype is its Malthusian fitness in the average environment minus the covariance of its fitness with population growth rate. The influence of fitness correlation between genotypes is illustrated by an analysis of the Haldane-Jayakar model of fluctuating selection on a single diallelic locus, and on two loci with additive effects on a quantitative character.     

Plasticity and heritability of morphological variation within and between parapatric stickleback demes

The threespine stickleback (Gasterosteus aculeatus) has emerged as an important model organism in evolutionary ecology, largely due to the repeated, parallel evolution of divergent morphotypes found in populations having colonized freshwater habitats. However, morphological divergence following colonization is not a universal phenomenon. We explore this in a large-scale estuarine ecosystem inhabited by two parapatric stickleback demes, each physiologically adapted to divergent osmoregulatory environments (fresh vs. saline waters). Using geometric morphometric analyses of wild-caught individuals, we detected significant differences between demes, in addition to sexual dimorphism, in body shape. However, rearing full-sib families from each deme under controlled, reciprocal salinity conditions revealed no differences between genotypes and highly significant environmental effects. It is also noteworthy that fish from both demes were fully plated, whether found in the wild or reared under reciprocal salinity conditions. Although we found significant heritability for body shape, we also noted significant direct environmental effects for many latent shape variables. Moreover, we found little evidence for diversifying selection acting on body size and shape (Q(ST) ). Nevertheless, uniform compressive variation did exceed neutral expectations, yet despite evidence of both allometry and genetic correlation with body length, we detected no correlated signatures of selection. Taken together, these results suggest that much of the morphological divergence observed in this system is the result of plastic responses to environmental variation rather than adaptive differentiation.     

Stable two-allele polymorphisms maintained by fluctuating fitnesses and seed banks: protecting the blues in Linanthus parryae

Motivated by data demonstrating fluctuating relative and absolute fitnesses for white- versus blue-flowered morphs of the desert annual Linanthus parryae, we present conditions under which temporally fluctuating selection and fluctuating contributions to a persistent seed bank will maintain a stable single-locus polymorphism. In L. parryae, blue flower color is determined by a single dominant allele. To disentangle the underlying diversity-maintaining mechanism from the mathematical complications associated with departures from Hardy-Weinberg genotype frequencies and dominance, we successively analyze a haploid model, a diploid model with three distinguishable genotypes, and a diploid model with complete dominance. For each model, we present conditions for the maintenance of a stable polymorphism, then use a diffusion approximation to describe the long-term fluctuations associated with these polymorphisms. Our protected polymorphism analyses show that a genotype whose arithmetic and geometric mean relative fitnesses are both less than one can persist if its relative fitness exceeds one in years that produce the most offspring. This condition is met by data from a population of L. parryae whose white morph has higher fitness (seed set) only in years of relatively heavy rain fall. The data suggest that the observed polymorphism may be explained by fluctuating selection. However, the yearly variation in flower color frequencies cannot be fully explained by our simple models, which ignore age structure and possible selection in the seed bank. We address two additional questions--one mathematical, the other biological--concerning the applicability of diffusion approximations to intense selection and the applicability of long-term predictions to datasets spanning decades for populations with long-lived seed banks.     

Eco-evolutionary maintenance of diversity in fluctuating environments

Growing evidence suggests that temporally fluctuating environments are important in maintaining variation both within and between species. To date, however, studies of genetic variation within a population have been largely conducted by evolutionary biologists (particularly population geneticists), while population and community ecologists have concentrated more on diversity at the species level. Despite considerable conceptual overlap, the commonalities and differences of these two alternative paradigms have yet to come under close scrutiny. Here, we review theoretical and empirical studies in population genetics and community ecology focusing on the ‘temporal storage effect’ and synthesise theories of diversity maintenance across different levels of biological organisation. Drawing on Chesson's coexistence theory, we explain how temporally fluctuating environments promote the maintenance of genetic variation and species diversity. We propose a further synthesis of the two disciplines by comparing models employing traditional frequency-dependent dynamics and those adopting density-dependent dynamics. We then address how temporal fluctuations promote genetic and species diversity simultaneously via rapid evolution and eco-evolutionary dynamics. Comparing and synthesising ecological and evolutionary approaches will accelerate our understanding of diversity maintenance in nature.     

Spatial heterogeneity in the strength of selection against deleterious alleles and the mutation load

According to current estimates of genomic deleterious mutation rates (which are often of the order 0.1-1) the mutation load (defined as a reduction in the average fitness of a population due to the presence of deleterious alleles) may be important in many populations. In this paper, I use multilocus simulations to explore the effect of spatial heterogeneity in the strength of selection against deleterious alleles on the mutation load (for example, it has been suggested that stressful environments may increase the strength of selection). These simulations show contrasted results: in some situations, spatial heterogeneity may greatly reduce the mutation load, due to the fact that migrants coming from demes under stronger selection carry relatively few deleterious alleles, and benefit from a strong advantage within demes under weaker selection (where individuals carry many more deleterious alleles); in other situations, however, deleterious alleles accumulate within demes under stronger selection, due to migration pressure from demes under weaker selection, leading to fitness erosion within those demes. This second situation is more frequent when the productivity of the different demes is proportional to their mean fitness. The effect of spatial heterogeneity is greatly reduced, however, when the response to environmental differences is inconsistent across loci.     

PHASE III OF WRIGHT'S SHIFTING BALANCE PROCESS AND THE VARIANCE AMONG DEMES IN MIGRATION RATE

Interdemic selection by the differential migration of individuals out from demes of high fitness and into demes of low fitness (Phase III) is one of the most controversial aspects of Wright's Shifting Balance Theory. I derive a relationship between Phase III migration and the interdemic selection differential, S, and show its potential effect on F_ST. The relationship reveals a diversifying effect of interdemic selection by Phase III migration on the genetic structure of a metapopulation. Using experimental metapopulations, I explored the effect of Phase III migration on F_ST by comparing the genetic variance among demes for two different patterns of migration: (1) island model migration and (2) Wright's Phase III migration. Although mean migration rates were the same, I found that the variance among demes in migration rate was significantly higher with Phase III than with island model migration. As a result, F_ST for the frequency of a neutral marker locus was higher with Phase III than it was with island model migration. By increasing F_ST, Phase III enhanced the genetic differentiation among demes for traits not subject to interdemic selection. This feature makes Wright's process different from individual selection which, by reducing effective population size, decreases the genetic variance within demes for all other traits. I discussed this finding in relation to the efficacy of Phase III and random migration for effecting peak shifts, and the contribution of genes with indirect effects to among‐deme variation.     

EVOLUTION OF PHENOTYPE–ENVIRONMENT ASSOCIATIONS BY GENETIC RESPONSES TO SELECTION AND PHENOTYPIC PLASTICITY IN A TEMPORALLY AUTOCORRELATED ENVIRONMENT

Covariation between population‐mean phenotypes and environmental variables, sometimes termed a “phenotype–environment association” (PEA), can result from phenotypic plasticity, genetic responses to natural selection, or both. PEAs can potentially provide information on the evolutionary dynamics of a particular set of populations, but this requires a full theoretical characterization of PEAs and their evolution. Here, we derive formulas for the expected PEA in a temporally fluctuating environment for a quantitative trait with a linear reaction norm. We compare several biologically relevant scenarios, including constant versus evolving plasticity, and the situation in which an environment affects both development and selection but at different time periods. We find that PEAs are determined not only by biological factors (e.g., magnitude of plasticity, genetic variation), but also environmental factors, such as the association between the environments of development and of selection, and in some cases the level of temporal autocorrelation. We also describe how a PEA can be used to estimate the relationship between an optimum phenotype and an environmental variable (i.e., the environmental sensitivity of selection), an important parameter for determining the extinction risk of populations experiencing environmental change. We illustrate this ability using published data on the predator‐induced morphological responses of tadpoles to predation risk.