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.

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.

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.

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.

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.

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.

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.

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.

Migration-selection balance at multiple Loci and selection on dominance and recombination.

A steady influx of a single deleterious multilocus genotype will impose genetic load on the resident population and leave multiple descendants carrying various numbers of the foreign alleles. Provided that the foreign types are rare at equilibrium, and all immigrant genes are eventually eliminated by selection, the population structure can be inferred explicitly from the branching process taking place within a single immigrant lineage. Unless the migration and recombination rates were high, this novel method gives a close approximation to the simulation with all possible multilocus genotypes considered. Once the load and the foreign genotypes frequencies are known, it becomes possible to estimate selection acting on the invading modifiers of (i) dominance and (ii) recombination rate on the foreign gene block. We found that the modifiers of the (i) type are able to invade faster than the type (ii) modifier, however, this result only applies in the strong selection/low migration/low recombination scenario. Varying the number of genes in the immigrant genotype can have a non-monotonic effect on the migration load and the modifier's invasion rate: although blocks carrying more genes can give rise to longer lineages, they also experience stronger selection pressure. The heaviest load is therefore imposed by the genotypes carrying moderate numbers of genes.

Evolution of transcription networks in response to temporal fluctuations.

Organisms respond to changes in their environment over a wide range of biological and temporal scales. Such phenotypic plasticity can involve developmental, behavioral, physiological, and genetic shifts. The adaptive value of a plastic response is known to depend on the nature of the information that is available to the organism as well as the direct and indirect costs of the plastic response. We modeled the dynamic process of simple gene regulatory networks as they responded to temporal fluctuations in environmental conditions. We simulated the evolution of networks to determine when genes that function solely as transcription factors, with no direct function of their own, are beneficial to the function of the network. When there is perfect information about the environment and there is no timing information to be extracted then there is no advantage to adding pure transcription factor genes to the network. In contrast, when there is either timing information that can be extracted or only indirect information about the current state of the environment then additional transcription factor genes improve the evolved network fitness.

Multiple routes to subfunctionalization and gene duplicate specialization.

Gene duplication is arguably the most significant source of new functional genetic material. A better understanding of the processes that lead to the stable incorporation of gene duplications into the genome is important both because it relates to interspecific differences in genome composition and because it can shed light on why some classes of gene are more prone to duplication than others. Typically, models of gene duplication consider the periods before duplication, during the spread and fixation of a new duplicate, and following duplication as distinct phases without a common underlying selective environment. I consider a scenario where a gene that is initially expressed in multiple contexts can undergo mutations that alter its expression profile or its functional coding sequence. The selective regime that acts on the functional output of the allele copies carried by an individual is constant. If there is a potential selective benefit to having different coding sequences expressed in each context, then, regardless of the constraints on functional variation at the single-locus gene, the waiting time until a gene duplication is incorporated goes down as population size increases.

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.

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.

Sex allocation based on relative and absolute condition.

Traditional models predict that organisms should allocate to sex based on their condition relative to the condition of their competitors, tracking shifts in mean condition in fluctuating environments, and maintaining an equilibrium sex ratio. In contrast, when individuals are constrained to define their condition absolutely, environmental fluctuations induce fluctuating sex ratios and the evolutionary loss of condition-dependent sex allocation in short-lived organisms. Here, we present a simulation model of temperature-dependent sex determination (TSD) in fluctuating environments that specifically examines the importance of relativity in defining individual condition. When relativity in condition is allowed to evolve, short-lived organisms evolve switchlike TSD reaction norms and define their condition relative to the annual temperature distribution, thus preventing biased cohort sex ratios in extreme years. Long-lived organisms also evolve switchlike reaction norms, but define condition less relatively and experience biased cohort sex ratios. The predictions are supported by data from painted turtles, where TSD reaction norms exhibit pivotal temperatures of sex determination that partially track mean annual temperature. Examining relativity in amniotic vertebrates provides a conceptual framework for multifactorial sex determination and suggests new ways of exploring adaptive hypotheses of sex allocation by focusing on the importance of frequency-dependent selection on sex.

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.

Dissecting selection on female mating preferences during secondary contact.

Population-specific preferences involved in premating isolation may be based on several different types of mating cues. Here, we compare the rates of spread of 12 different mating preferences that reflect preferences for local adaptation, male condition, and reinforcement. We introduce methods to dissect the components of the rate of spread to determine why certain mating preferences spread more quickly than others. We confirm the result that female preferences based on population-specific markers alone always spread faster than female preferences based only on a single local adaptation locus, regardless of the strength of natural selection on hybrid incompatibility. However, we find that this occurs for different reasons depending on the strength of selection against hybrids. Female preferences based on total male condition also achieved high rates of spread, suggesting that preferences for condition-dependent male displays may evolve under reinforcement scenarios.

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.