Selection on modifiers of genetic architecture under migration load.

Gene flow between populations adapting to differing local environmental conditions might be costly because individuals can disperse to habitats where their survival is low or because they can reproduce with locally maladapted individuals. The amount by which the mean relative population fitness is kept below one creates an opportunity for modifiers of the genetic architecture to spread due to selection. Prior work that separately considered modifiers changing dispersal, recombination rates, or altering dominance or epistasis, has typically focused on the direction of selection rather than its absolute magnitude. We here develop methods to determine the strength of selection on modifiers of the genetic architecture, including modifiers of the dispersal rate, in populations that have previously evolved local adaptation. We consider scenarios with up to five loci contributing to local adaptation and derive a new model for the deterministic spread of modifiers. We find that selection for modifiers of epistasis and dominance is stronger than selection for decreased recombination, and that selection for partial reductions in recombination are extremely weak, regardless of the number of loci contributing to local adaptation. The spread of modifiers that reduce dispersal depends on the number of loci, epistasis and extent of local adaptation in the ancestral population. We identify a novel effect, that modifiers of dominance are more strongly selected when they are unlinked to the locus that they modify. These findings help explain population differentiation and reproductive isolation and provide a benchmark to compare selection on modifiers under finite population sizes and demographic stochasticity.

Slower environmental change hinders adaptation from standing genetic variation.

Evolutionary responses to environmental change depend on the time available for adaptation before environmental degradation leads to extinction. Explicit tests of this relationship are limited to microbes where adaptation usually depends on the sequential fixation of de novo mutations, excluding standing variation for genotype-by-environment fitness interactions that should be key for most natural species. For natural species evolving from standing genetic variation, adaptation at slower rates of environmental change may be impeded since the best genotypes at the most extreme environments can be lost during evolution due to genetic drift or founder effects. To address this hypothesis, we perform experimental evolution with self-fertilizing populations of the nematode Caenorhabditis elegans and develop an inference model to describe natural selection on extant genotypes under environmental change. Under a sudden environmental change, we find that selection rapidly increases the frequency of genotypes with high fitness in the most extreme environment. In contrast, under a gradual environmental change selection first favors genotypes that are worse at the most extreme environment. We demonstrate with a second set of evolution experiments that, as a consequence of slower environmental change and thus longer periods to reach the most extreme environments, genetic drift and founder effects can lead to the loss of the most beneficial genotypes. We further find that maintenance of standing genetic variation can retard the fixation of the best genotypes in the most extreme environment because of interference between them. Taken together, these results show that slower environmental change can hamper adaptation from standing genetic variation and they support theoretical models indicating that standing variation for genotype-by-environment fitness interactions critically alters the pace and outcome of adaptation under environmental change.

How differing modes of non-genetic inheritance affect population viability in fluctuating environments.

Different modes of non-genetic inheritance are expected to affect population persistence in fluctuating environments. We here analyse Caenorhabditis elegans density-independent per capita growth rate time series on 36 populations experiencing six controlled sequences of challenging oxygen level fluctuations across 60 generations, and parameterise competing models of non-genetic inheritance in order to explain observed dynamics. Our analysis shows that phenotypic plasticity and anticipatory maternal effects are sufficient to explain growth rate dynamics, but that a carryover model where 'epigenetic' memory is imperfectly transmitted and might be reset at each generation is a better fit to the data. We further find that this epigenetic memory is asymmetric since it is kept for longer when populations are exposed to the more challenging environment. Our analysis suggests that population persistence in fluctuating environments depends on the non-genetic inheritance of phenotypes whose expression is regulated across multiple generations.

Adaptation to Temporally Fluctuating Environments by the Evolution of Maternal Effects.

All organisms live in temporally fluctuating environments. Theory predicts that the evolution of deterministic maternal effects (i.e., anticipatory maternal effects or transgenerational phenotypic plasticity) underlies adaptation to environments that fluctuate in a predictably alternating fashion over maternal-offspring generations. In contrast, randomizing maternal effects (i.e., diversifying and conservative bet-hedging), are expected to evolve in response to unpredictably fluctuating environments. Although maternal effects are common, evidence for their adaptive significance is equivocal since they can easily evolve as a correlated response to maternal selection and may or may not increase the future fitness of offspring. Using the hermaphroditic nematode Caenorhabditis elegans, we here show that the experimental evolution of maternal glycogen provisioning underlies adaptation to a fluctuating normoxia-anoxia hatching environment by increasing embryo survival under anoxia. In strictly alternating environments, we found that hermaphrodites evolved the ability to increase embryo glycogen provisioning when they experienced normoxia and to decrease embryo glycogen provisioning when they experienced anoxia. At odds with existing theory, however, populations facing irregularly fluctuating normoxia-anoxia hatching environments failed to evolve randomizing maternal effects. Instead, adaptation in these populations may have occurred through the evolution of fitness effects that percolate over multiple generations, as they maintained considerably high expected growth rates during experimental evolution despite evolving reduced fecundity and reduced embryo survival under one or two generations of anoxia. We develop theoretical models that explain why adaptation to a wide range of patterns of environmental fluctuations hinges on the existence of deterministic maternal effects, and that such deterministic maternal effects are more likely to contribute to adaptation than randomizing maternal effects.

What Kind of Maternal Effects Can Be Selected For in Fluctuating Environments?

Just as phenotypic plasticity can evolve when developing individuals get informational cues about their future adult environment, deterministic maternal effects, where offspring trait values depend on the maternal environment, can evolve when mothers gain reliable information about the environments their offspring will face. Randomizing maternal effects (a type of diversifying bet hedging), where offspring trait values are randomized, can evolve by natural selection even when information about future environments is unavailable. We investigate selection on both randomizing and deterministic maternal effects in environments that show correlated fluctuations between two environmental states. We compare the strength of selection for deterministic and randomizing maternal effects and explicitly consider maternal fitness costs of producing offspring with different phenotypes. Only a small set of environmental parameters allow randomizing maternal effects to outcompete deterministic maternal effects; not only must there be little or no information available about future environments, but the frequency of each environment must fall within a narrow range. By contrast, deterministic maternal effects can always invade an ancestral state lacking a maternal effect even if the amount of environmental information available is low. The long-term outcome may involve offspring trait value randomization but only if trait values first evolve to cause extreme differences in environment-specific fitness. Overall, deterministic maternal effects are more likely to evolve by natural selection than randomizing maternal effects.

Rapid DNA loss as a counterbalance to genome expansion through retrotransposon proliferation in plants.

Transposable elements, particularly LTR-retrotransposons, comprise the primary vehicle for genome size expansion in plants, while DNA removal through illegitimate recombination and intrastrand homologous recombination serve as the most important counteracting forces to plant genomic obesity. Despite extensive research, the relative impact of these opposing forces and hence the directionality of genome size change remains unknown. In Gossypium (cotton), the 3-fold genome size variation among diploids is due largely to copy number variation of the gypsy-like retrotransposon Gorge3. Here we combine comparative sequence analysis with a modeling approach to study the directionality of genome size change in Gossypium. We demonstrate that the rate of DNA removal in the smaller genomes is sufficient to reverse genome expansion through Gorge3 proliferation. These data indicate that rates of DNA loss can be highly variable even within a single plant genus, and that the known mechanisms of DNA loss can indeed reverse the march toward genomic obesity.

Modeling the metabolic evolution of mixotrophic phytoplankton in response to rising ocean surface temperatures.

BACKGROUND: Climate change is expected to lead to warming in ocean surface temperatures which will have unequal effects on the rates of photosynthesis and heterotrophy. As a result of this changing metabolic landscape, directional phenotypic evolution will occur, with implications that cascade up to the ecosystem level. While mixotrophic phytoplankton, organisms that combine photosynthesis and heterotrophy to meet their energetic and nutritional needs, are expected to become more heterotrophic with warmer temperatures due to heterotrophy increasing at a faster rate than photosynthesis, it is unclear how evolution will influence how these organisms respond to warmer temperatures. In this study, we used adaptive dynamics to model the consequences of temperature-mediated increases in metabolic rates for the evolution of mixotrophic phytoplankton, focusing specifically on phagotrophic mixotrophs. RESULTS: We find that mixotrophs tend to evolve to become more reliant on phagotrophy as temperatures rise, leading to reduced prey abundance through higher grazing rates. However, if prey abundance becomes too low, evolution favors greater reliance on photosynthesis. These responses depend upon the trade-off that mixotrophs experience between investing in photosynthesis and phagotrophy. Mixotrophs with a convex trade-off maintain mixotrophy over the greatest range of temperatures; evolution in these "generalist" mixotrophs was found to exacerbate carbon cycle impacts, with evolving mixotrophs exhibiting increased sensitivity to rising temperature. CONCLUSIONS: Our results show that mixotrophs may respond more strongly to climate change than predicted by phenotypic plasticity alone due to evolutionary shifts in metabolic investment. However, the type of metabolic trade-off experienced by mixotrophs as well as ecological feedback on prey abundance may ultimately limit the extent of evolutionary change along the heterotrophy-phototrophy spectrum.

The rate of multi-step evolution in Moran and Wright-Fisher populations.

Several groups have recently modeled evolutionary transitions from an ancestral allele to a beneficial allele separated by one or more intervening mutants. The beneficial allele can become fixed if a succession of intermediate mutants are fixed or alternatively if successive mutants arise while the previous intermediate mutant is still segregating. This latter process has been termed stochastic tunneling. Previous work has focused on the Moran model of population genetics. I use elementary methods of analyzing stochastic processes to derive the probability of tunneling in the limit of large population size for both Moran and Wright-Fisher populations. I also show how to efficiently obtain numerical results for finite populations. These results show that the probability of stochastic tunneling is twice as large under the Wright-Fisher model as it is under the Moran model.

Testing the adaptive value of sporulation in budding yeast using experimental evolution.

Saccharomyces yeast grow through mitotic cell division, converting resources into biomass. When cells experience starvation, sporulation is initiated and meiosis produces haploid cells inside a protective ascus. The protected spore state does not acquire resources and is partially protected from desiccation, heat, and caustic chemicals. Because cells cannot both be protected and acquire resources simultaneously, committing to sporulation represents a trade-off between current and future reproduction. Recent work has suggested that passaging through insect guts selects for spore formation, as surviving insect ingestion represents a major way that yeasts are vectored to new food sources. We subjected replicate populations from five yeast strains to passaging through insects, and evolved control populations by pipette passaging. We assayed populations for their propensity to sporulate after resource depletion. We found that ancestral domesticated strains produced fewer spores, and all strains evolved increased spore production in response to passaging through flies, but domesticated strains responded less. Exposure to flies led to a more rapid shift to sporulation that was more extreme in wild-derived strains. Our results indicate that insect passaging selects for spore production and suggest that domestication led to genetic canalization of the response to cues in the environment and initiation of sporulation.

A garter snake transcriptome: pyrosequencing, de novo assembly, and sex-specific differences.

BACKGROUND: The reptiles, characterized by both diversity and unique evolutionary adaptations, provide a comprehensive system for comparative studies of metabolism, physiology, and development. However, molecular resources for ectothermic reptiles are severely limited, hampering our ability to study the genetic basis for many evolutionarily important traits such as metabolic plasticity, extreme longevity, limblessness, venom, and freeze tolerance. Here we use massively parallel sequencing (454 GS-FLX Titanium) to generate a transcriptome of the western terrestrial garter snake (Thamnophis elegans) with two goals in mind. First, we develop a molecular resource for an ectothermic reptile; and second, we use these sex-specific transcriptomes to identify differences in the presence of expressed transcripts and potential genes of evolutionary interest. RESULTS: Using sex-specific pools of RNA (one pool for females, one pool for males) representing 7 tissue types and 35 diverse individuals, we produced 1.24 million sequence reads, which averaged 366 bp in length after cleaning. Assembly of the cleaned reads from both sexes with NEWBLER and MIRA resulted in 96,379 contigs containing 87% of the cleaned reads. Over 34% of these contigs and 13% of the singletons were annotated based on homology to previously identified proteins. From these homology assignments, additional clustering, and ORF predictions, we estimate that this transcriptome contains ~13,000 unique genes that were previously identified in other species and over 66,000 transcripts from unidentified protein-coding genes. Furthermore, we use a graph-clustering method to identify contigs linked by NEWBLER-split reads that represent divergent alleles, gene duplications, and alternatively spliced transcripts. Beyond gene identification, we identified 95,295 SNPs and 31,651 INDELs. From these sex-specific transcriptomes, we identified 190 genes that were only present in the mRNA sequenced from one of the sexes (84 female-specific, 106 male-specific), and many highly variable genes of evolutionary interest. CONCLUSIONS: This is the first large-scale, multi-organ transcriptome for an ectothermic reptile. This resource provides the most comprehensive set of EST sequences available for an individual ectothermic reptile species, increasing the number of snake ESTs 50-fold. We have identified genes that appear to be under evolutionary selection and those that are sex-specific. This resource will assist studies on gene expression and comparative genomics, and will facilitate the study of evolutionarily important traits at the molecular level.

The evolutionary origins of gene regulation.

One way that organisms cope with constantly changing physical and biological conditions is by regulating the expression of genes and thereby altering protein production. Clearly, altering the protein production to match the environmental demands can be adaptive, but there may be evolutionary barriers to the transition from constitutive expression to regulated expression. In particular, down-regulating a gene when it is not needed means that there will necessarily be a delay in protein production when the protein is up-regulated in the future. We develop a model of simple gene regulation in response to randomly changing environmental conditions. We calculate the long-term behavior of gene expression and determine the fitness consequences of changes in the gene regulation. We then embed this model into a population genetic framework in order to determine the conditions that allow populations to evolve environment-specific transcription rates. The population genetic model follows the evolutionary transition from constitutive expression to regulated expression. There are three distinct possible evolutionary outcomes. The gene may be stuck in the always on position, the gene may first evolve to an intermediate constitutive expression level and then evolve regulation, or regulation can evolve directly from the ancestral state in a smooth fashion. Regulation is most likely to evolve when the costs of mis-expression are low and the transcript decay rate is high. This suggests that genes that have less severe reductions in fitness when mis-expressed are more likely to initially evolve regulation.

Gene expression dynamics in randomly varying environments.

A simple model of gene regulation in response to stochastically changing environmental conditions is developed and analyzed. The model consists of a differential equation driven by a continuous time 2-state Markov process. The density function of the resulting process converges to a beta distribution. We show that the moments converge to their stationary values exponentially in time. Simulations of a two-stage process where protein production depends on mRNA concentrations are also presented demonstrating that protein concentration tracks the environment whenever the rate of protein turnover is larger than the rate of environmental change. Single-celled organisms are therefore expected to have relatively high mRNA and protein turnover rates for genes that respond to environmental fluctuations.

Partial Selfing Can Reduce Genetic Loads While Maintaining Diversity During Experimental Evolution.

Partial selfing, whereby self- and cross- fertilization occur in populations at intermediate frequencies, is generally thought to be evolutionarily unstable. Yet, it is found in natural populations. This could be explained if populations with partial selfing are able to reduce genetic loads and the possibility for inbreeding depression while keeping genetic diversity that may be important for future adaptation. To address this hypothesis, we compare the experimental evolution of Caenorhabditis elegans populations under partial selfing, exclusive selfing or predominant outcrossing, while they adapt to osmotically challenging conditions. We find that the ancestral genetic load, as measured by the risk of extinction upon inbreeding by selfing, is maintained as long as outcrossing is the main reproductive mode, but becomes reduced otherwise. Analysis of genome-wide single-nucleotide polymorphisms (SNPs) during experimental evolution and among the inbred lines that survived enforced inbreeding indicates that populations with predominant outcrossing or partial selfing maintained more genetic diversity than expected with neutrality or purifying selection. We discuss the conditions under which this could be explained by the presence of recessive deleterious alleles and/or overdominant loci. Taken together, our observations suggest that populations evolving under partial selfing can gain some of the benefits of eliminating unlinked deleterious recessive alleles and also the benefits of maintaining genetic diversity at partially dominant or overdominant loci that become associated due to variance of inbreeding levels.

Mutual information reveals variation in temperature-dependent sex determination in response to environmental fluctuation, lifespan and selection.

Quantifying the degree to which sex determination depends on the environment can yield insight into the evolution, ecological dynamics, and functional aspects of sex determination. In temperature-dependent sex determination (TSD), theory often predicts a complete dependence of sex on temperature, with a switch-like reaction norm. However, empirical data suggest more shallow relationships between sex and temperature. Here, we demonstrate the usefulness of an index, mutual information (MI), to reflect the degree of temperature dependence in sex. MI depends on both the shape of a reaction norm and the natural temperature variation, thus providing a measure of TSD that is ecologically dependent. We demonstrate that increased lifespan and decreased environmental fluctuation predict reaction norms with high MI (switch-like). However, mutation and weaker selection on sex-specific performance reduce average MI in a population, suggesting that mutation-selection balance can resolve some of the conflict between theoretical predictions of individual-based optimality and population-based empirical results. The MI index allows clear comparison of TSD across life histories and habitats and reveals functional similarities between reaction norms that may appear different. The model provides testable predictions for TSD across populations, namely that MI should increase with lifespan and decrease with historical environmental fluctuations.

Quasi-species evolution in subdivided populations favours maximally deleterious mutations.

Most models of quasi-species evolution predict that populations will evolve to occupy areas of sequence space with the greatest concentration of neutral sequences, thus minimizing the deleterious mutation rate and creating mutationally 'robust' genomes. In contrast, empirical studies of the principal model of quasi-species evolution, RNA viruses, suggest that the effects of deleterious mutations are more severe than in similar DNA-based microbes. We demonstrate that populations divided into discrete patches connected by dispersal may favour genotypes where the deleterious effect of non-neutral mutations is maximized. This effect is especially strong in the absence of back mutation and when the amount of time spent in hosts prior to dispersal is intermediate. Our results indicate that RNA viruses that produce acute infections initiated by a small number of virions are expected to evolve fragile genetic architectures when compared with other RNA viruses.

Allelic divergence precedes and promotes gene duplication.

One of the striking observations from recent whole-genome comparisons is that changes in the number of specialized genes in existing gene families, as opposed to novel taxon-specific gene families, are responsible for the majority of the difference in genome composition between major taxa. Previous models of duplicate gene evolution focused primarily on the role that neutral processes can play in evolutionary divergence after the duplicates are already fixed in the population. By instead including the entire cycle of duplication and divergence, we show that specialized functions are most likely to evolve through strong selection acting on segregating alleles at a single locus, even before the duplicate arises. We show that the fitness relationships that allow divergent alleles to evolve at a single locus largely overlap with the conditions that allow divergence of previously duplicated genes. Thus, a solution to the paradox of the origin of organismal complexity via the expansion of gene families exists in the form of the deterministic spread of novel duplicates via natural selection.

The opportunity for canalization and the evolution of genetic networks.

There has been a recent revival of interest in how genetic interactions evolve, spurred on by an increase in our knowledge of genetic interactions at the molecular level. Empirical work on genetic networks has revealed a surprising amount of robustness to perturbations, suggesting that robustness is an evolved feature of genetic networks. Here, we derive a general model for the evolution of canalization that can incorporate any form of perturbation. We establish an upper bound to the strength of selection on canalization that is approximately equal to the fitness load in the system. This method makes it possible to compare different forms of perturbation, including genetic, developmental, and environmental effects. In general, load that arises from mutational processes is low because the mutation rate is itself low. Mutation load can create selection for canalization in a small network that can be achieved through dominance evolution or gene duplication, and in each case selection for canalization is weak at best. In larger genetic networks, selection on genetic canalization can be reasonably strong because larger networks have higher mutational load. Because load induced through migration, segregation, developmental noise, and environmental variance is not mutation limited, each can cause strong selection for canalization.

Matings Systems and the Evolution of Niche Breadth.

Several theoretical studies of niche breadth evolution have shown that niche breadth expansion can only occur under a limited set of conditions. These studies have assumed that mating pairs form at random within spatial subpopulations. I show here that nonrandom mating within spatial subpopulations can greatly alter niche breadth evolution. Niche breadth expansion is restricted by assortative mating but is permitted by sexual selection. When new populations are formed by a single pulse of immigrants, the ability to persist in novel but similar environments is enhanced by sexual selection. When new populations receive a constant stream of immigrants, sexual selection allows the evolution of increased niche breadth even when the novel environment is very different from the ancestral environment. These results suggest that species lineages that exhibit sexual selection will have broader niches than species lineages that mate randomly or assortatively.

A general theory for the evolutionary dynamics of virulence.

Most theory on the evolution of virulence is based on a game-theoretic approach. One potential shortcoming of this approach is that it does not allow the prediction of the evolutionary dynamics of virulence. Such dynamics are of interest for several reasons: for experimental tests of theory, for the development of useful virulence management protocols, and for understanding virulence evolution in situations where the epidemiological dynamics never reach equilibrium and/or when evolutionary change occurs on a timescale comparable to that of the epidemiological dynamics. Here we present a general theory similar to that of quantitative genetics in evolutionary biology that allows for the easy construction of models that include both within-host mutation as well as superinfection and that is capable of predicting both the short- and long-term evolution of virulence. We illustrate the generality and intuitive appeal of the theory through a series of examples showing how it can lead to transparent interpretations of the selective forces governing virulence evolution. It also leads to novel predictions that are not possible using the game-theoretic approach. The general theory can be used to model the evolution of other pathogen traits as well.

Older males signal more reliably.

The hypothesis that females prefer older males because they have higher mean fitness than younger males has been the centre of recent controversy. These discussions have focused on the success of a female who prefers males of a particular age class when age cues, but not quality cues, are available. Thus, if the distribution of male quality changes with age, such that older males have on average genotypes with higher fitness than younger males, then a female who mates with older males has fitter offspring, which allows the female preference to spread through a genetic correlation. We develop a general model for male display in a species with multiple reproductive bouts that allows us to identify the conditions that promote reliable signalling within an age class. Because males have opportunities for future reproduction, they will reduce their levels of advertising compared with a semelparous species. In addition, because higher-quality males have more future reproduction, they will reduce their advertising more than low-quality males. Thus, the conditions for reliable signalling in a semelparous organism are generally not sufficient to produce reliable signalling in species with multiple reproductive bouts. This result is due to the possibility of future reproduction so that, as individuals age and the opportunities for future reproduction fade, signalling becomes more reliable. This provides a novel rationale for female preference for older mates; older males reveal more information in their sexual displays.