Application of simultaneous selective pressures slows adaptation.

Beneficial mutations that arise in an evolving asexual population may compete or interact in ways that alter the overall rate of adaptation through mechanisms such as clonal or functional interference. The application of multiple selective pressures simultaneously may allow for a greater number of adaptive mutations, increasing the opportunities for competition between selectively advantageous alterations, and thereby reducing the rate of adaptation. We evolved a strain of Saccharomyces cerevisiae that could not produce its own histidine or uracil for ~500 generations under one or three selective pressures: limitation of the concentration of glucose, histidine, and/or uracil in the media. The rate of adaptation was obtained by measuring evolved relative fitness using competition assays. Populations evolved under a single selective pressure showed a statistically significant increase in fitness on those pressures relative to the ancestral strain, but the populations evolved on all three pressures did not show a statistically significant increase in fitness over the ancestral strain on any single pressure. Simultaneously limiting three essential nutrients for a population of S. cerevisiae effectively slows the rate of evolution on any one of the three selective pressures applied, relative to the single selective pressure cases. We identify possible mechanisms for fitness changes seen between populations evolved on one or three limiting nutrient pressures by high-throughput sequencing. Adding multiple selective pressures to evolving disease like cancer and infectious diseases could reduce the rate of adaptation and thereby may slow disease progression, prolong drug efficacy and prevent deaths.

Migration promotes mutator alleles in subdivided populations.

Mutator alleles that elevate the genomic mutation rate may invade nonrecombining populations by hitchhiking with beneficial mutations. Mutators have been repeatedly observed to take over adapting laboratory populations and have been found at high frequencies in both microbial pathogen and cancer populations in nature. Recently, we have shown that mutators are only favored by selection in sufficiently large populations and transition to being disfavored as population size decreases. This population size-dependent sign inversion in selective effect suggests that population structure may also be an important determinant of mutation rate evolution. Although large populations may favor mutators, subdividing such populations into sufficiently small subpopulations (demes) might effectively inhibit them. On the other hand, migration between small demes that otherwise inhibit hitchhiking may promote mutator fixation in the whole metapopulation. Here, we use stochastic, agent-based simulations and evolution experiments with the yeast Saccharomyces cerevisiae to show that mutators can, indeed, be favored by selection in subdivided metapopulations composed of small demes connected by sufficient migration. In fact, we show that population structure plays a previously unsuspected role in promoting mutator success in subdivided metapopulations when migration is rare.

Sign of selection on mutation rate modifiers depends on population size.

The influence of population size (N) on natural selection acting on alleles that affect fitness has been understood for almost a century. As N declines, genetic drift overwhelms selection and alleles with direct fitness effects are rendered neutral. Often, however, alleles experience so-called indirect selection, meaning they affect not the fitness of an individual but the fitness distribution of its offspring. Some of the best-studied examples of indirect selection include alleles that modify aspects of the genetic system such as recombination and mutation rates. Here, we use analytics, simulations, and experimental populations of Saccharomyces cerevisiae to examine the influence of N on indirect selection acting on alleles that increase the genomic mutation rate (mutators). Mutators experience indirect selection via genomic associations with beneficial and deleterious mutations they generate. We show that, as N declines, indirect selection driven by linked beneficial mutations is overpowered by drift before drift can neutralize the cost of the deleterious load. As a result, mutators transition from being favored by indirect selection in large populations to being disfavored as N declines. This surprising phenomenon of sign inversion in selective effect demonstrates that indirect selection on mutators exhibits a profound and qualitatively distinct dependence on N.

Systematically Altering Bacterial SOS Activity under Stress Reveals Therapeutic Strategies for Potentiating Antibiotics.

The bacterial SOS response is a DNA damage repair network that is strongly implicated in both survival and acquired drug resistance under antimicrobial stress. The two SOS regulators, LexA and RecA, have therefore emerged as potential targets for adjuvant therapies aimed at combating resistance, although many open questions remain. For example, it is not well understood whether SOS hyperactivation is a viable therapeutic approach or whether LexA or RecA is a better target. Furthermore, it is important to determine which antimicrobials could serve as the best treatment partners with SOS-targeting adjuvants. Here we derived Escherichia coli strains that have mutations in either lexA or recA genes in order to cover the full spectrum of possible SOS activity levels. We then systematically analyzed a wide range of antimicrobials by comparing the mean inhibitory concentrations (MICs) and induced mutation rates for each drug-strain combination. We first show that significant changes in MICs are largely confined to DNA-damaging antibiotics, with strains containing a constitutively repressed SOS response impacted to a greater extent than hyperactivated strains. Second, antibiotic-induced mutation rates were suppressed when SOS activity was reduced, and this trend was observed across a wider spectrum of antibiotics. Finally, perturbing either LexA or RecA proved to be equally viable strategies for targeting the SOS response. Our work provides support for multiple adjuvant strategies, while also suggesting that the combination of an SOS inhibitor with a DNA-damaging antibiotic could offer the best potential for lowering MICs and decreasing acquired drug resistance. IMPORTANCE Our antibiotic arsenal is becoming depleted, in part, because bacteria have the ability to rapidly adapt and acquire resistance to our best agents. The SOS pathway, a widely conserved DNA damage stress response in bacteria, is activated by many antibiotics and has been shown to play central role in promoting survival and the evolution of resistance under antibiotic stress. As a result, targeting the SOS response has been proposed as an adjuvant strategy to revitalize our current antibiotic arsenal. However, the optimal molecular targets and partner antibiotics for such an approach remain unclear. In this study, focusing on the two key regulators of the SOS response, LexA and RecA, we provide the first comprehensive assessment of how to target the SOS response in order to increase bacterial susceptibility and reduce mutagenesis under antibiotic treatment.

Natural selection promotes antigenic evolvability.

The hypothesis that evolvability - the capacity to evolve by natural selection - is itself the object of natural selection is highly intriguing but remains controversial due in large part to a paucity of direct experimental evidence. The antigenic variation mechanisms of microbial pathogens provide an experimentally tractable system to test whether natural selection has favored mechanisms that increase evolvability. Many antigenic variation systems consist of paralogous unexpressed 'cassettes' that recombine into an expression site to rapidly alter the expressed protein. Importantly, the magnitude of antigenic change is a function of the genetic diversity among the unexpressed cassettes. Thus, evidence that selection favors among-cassette diversity is direct evidence that natural selection promotes antigenic evolvability. We used the Lyme disease bacterium, Borrelia burgdorferi, as a model to test the prediction that natural selection favors amino acid diversity among unexpressed vls cassettes and thereby promotes evolvability in a primary surface antigen, VlsE. The hypothesis that diversity among vls cassettes is favored by natural selection was supported in each B. burgdorferi strain analyzed using both classical (dN/dS ratios) and Bayesian population genetic analyses of genetic sequence data. This hypothesis was also supported by the conservation of highly mutable tandem-repeat structures across B. burgdorferi strains despite a near complete absence of sequence conservation. Diversification among vls cassettes due to natural selection and mutable repeat structures promotes long-term antigenic evolvability of VlsE. These findings provide a direct demonstration that molecular mechanisms that enhance evolvability of surface antigens are an evolutionary adaptation. The molecular evolutionary processes identified here can serve as a model for the evolution of antigenic evolvability in many pathogens which utilize similar strategies to establish chronic infections.

Genomic mutation rates that neutralize adaptive evolution and natural selection.

When mutation rates are low, natural selection remains effective, and increasing the mutation rate can give rise to an increase in adaptation rate. When mutation rates are high to begin with, however, increasing the mutation rate may have a detrimental effect because of the overwhelming presence of deleterious mutations. Indeed, if mutation rates are high enough: (i) adaptive evolution may be neutralized, resulting in a zero (or negative) adaptation rate despite the continued availability of adaptive and/or compensatory mutations, or (ii) natural selection may be neutralized, because the fitness of lineages bearing adaptive and/or compensatory mutations--whether established or newly arising--is eroded by excessive mutation, causing such lineages to decline in frequency. We apply these two criteria to a standard model of asexual adaptive evolution and derive mathematical expressions--some new, some old in new guise--delineating the mutation rates under which either adaptive evolution or natural selection is neutralized. The expressions are simple and require no a priori knowledge of organism- and/or environment-specific parameters. Our discussion connects these results to each other and to previous theory, showing convergence or equivalence of the different results in most cases.

Cancer in light of experimental evolution.

Cancer initiation, progression, and the emergence of therapeutic resistance are evolutionary phenomena of clonal somatic cell populations. Studies in microbial experimental evolution and the theoretical work inspired by such studies are yielding deep insights into the evolutionary dynamics of clonal populations, yet there has been little explicit consideration of the relevance of this rapidly growing field to cancer biology. Here, we examine how the understanding of mutation, selection, and spatial structure in clonal populations that is emerging from experimental evolution may be applicable to cancer. Along the way, we discuss some significant ways in which cancer differs from the model systems used in experimental evolution. Despite these differences, we argue that enhanced prediction and control of cancer may be possible using ideas developed in the context of experimental evolution, and we point out some prospects for future research at the interface between these traditionally separate areas.

Contrasting dynamics of a mutator allele in asexual populations of differing size.

Mutators have been shown to hitchhike in asexual populations when the anticipated beneficial mutation supply rate of the mutator subpopulation, NU(b) (for subpopulation of size N and beneficial mutation rate U(b)) exceeds that of the wild-type subpopulation. Here, we examine the effect of total population size on mutator dynamics in asexual experimental populations of Saccharomyces cerevisiae. Although mutators quickly hitchhike to fixation in smaller populations, mutator fixation requires more and more time as population size increases; this observed delay in mutator hitchhiking is consistent with the expected effect of clonal interference. Interestingly, despite their higher beneficial mutation supply rate, mutators are supplanted by the wild type in very large populations. We postulate that this striking reversal in mutator dynamics is caused by an interaction between clonal interference, the fitness cost of the mutator allele, and infrequent large-effect beneficial mutations in our experimental populations. Our work thus identifies a potential set of circumstances under which mutator hitchhiking can be inhibited in natural asexual populations, despite recent theoretical predictions that such populations should have a net tendency to evolve ever-higher genomic mutation rates.

Real time forecasting of near-future evolution.

A metaphor for adaptation that informs much evolutionary thinking today is that of mountain climbing, where horizontal displacement represents change in genotype, and vertical displacement represents change in fitness. If it were known a priori what the 'fitness landscape' looked like, that is, how the myriad possible genotypes mapped onto fitness, then the possible paths up the fitness mountain could each be assigned a probability, thus providing a dynamical theory with long-term predictive power. Such detailed genotype-fitness data, however, are rarely available and are subject to change with each change in the organism or in the environment. Here, we take a very different approach that depends only on fitness or phenotype-fitness data obtained in real time and requires no a priori information about the fitness landscape. Our general statistical model of adaptive evolution builds on classical theory and gives reasonable predictions of fitness and phenotype evolution many generations into the future.

Mutator dynamics in sexual and asexual experimental populations of yeast.

BACKGROUND: In asexual populations, mutators may be expected to hitchhike with associated beneficial mutations. In sexual populations, recombination is predicted to erode such associations, inhibiting mutator hitchhiking. To investigate the effect of recombination on mutators experimentally, we compared the frequency dynamics of a mutator allele (msh2Δ) in sexual and asexual populations of Saccharomyces cerevisiae. RESULTS: Mutator strains increased in frequency at the expense of wild-type strains in all asexual diploid populations, with some approaching fixation in 150 generations of propagation. Over the same period of time, mutators declined toward loss in all corresponding sexual diploid populations as well as in haploid populations propagated asexually. CONCLUSIONS: We report the first experimental investigation of mutator dynamics in sexual populations. We show that a strong mutator quickly declines in sexual populations while hitchhiking to high frequency in asexual diploid populations, as predicted by theory. We also show that the msh2Δ mutator has a high and immediate realized cost that is alone sufficient to explain its decline in sexual populations. We postulate that this cost is indirect; namely, that it is due to a very high rate of recessive lethal or strongly deleterious mutation. However, we cannot rule out the possibility that msh2Δ also has unknown directly deleterious effects on fitness, and that these effects may differ between haploid asexual and sexual populations. Despite these reservations, our results prompt us to speculate that the short-term cost of highly deleterious recessive mutations can be as important as recombination in preventing mutator hitchhiking in sexual populations.

Competition between high- and higher-mutating strains of Escherichia coli.

Experimental studies have shown that a mutator allele can readily hitchhike to fixation with beneficial mutations in an asexual population having a low, wild-type mutation rate. Here, we show that a genotype bearing two mutator alleles can supplant a population already fixed for one mutator allele. Our results provide experimental support for recent theory predicting that mutator alleles will tend to accumulate in asexual populations by hitchhiking with beneficial mutations, causing an ever-higher genomic mutation rate.

Heterochronic evolution reveals modular timing changes in budding yeast transcriptomes.

BACKGROUND: Gene expression is a dynamic trait, and the evolution of gene regulation can dramatically alter the timing of gene expression without greatly affecting mean expression levels. Moreover, modules of co-regulated genes may exhibit coordinated shifts in expression timing patterns during evolutionary divergence. Here, we examined transcriptome evolution in the dynamical context of the budding yeast cell-division cycle, to investigate the extent of divergence in expression timing and the regulatory architecture underlying timing evolution. RESULTS: Using a custom microarray platform, we obtained 378 measurements for 6,263 genes over 18 timepoints of the cell-division cycle in nine strains of S. cerevisiae and one strain of S. paradoxus. Most genes show significant divergence in expression dynamics at all scales of transcriptome organization, suggesting broad potential for timing changes. A model test comparing expression level evolution versus timing evolution revealed a better fit with timing evolution for 82% of genes. Analysis of shared patterns of timing evolution suggests the existence of seven dynamically-autonomous modules, each of which shows coherent evolutionary timing changes. Analysis of transcription factors associated with these gene modules suggests a modular pleiotropic source of divergence in expression timing. CONCLUSIONS: We propose that transcriptome evolution may generally entail changes in timing (heterochrony) rather than changes in levels (heterometry) of expression. Evolution of gene expression dynamics may involve modular changes in timing control mediated by module-specific transcription factors. We hypothesize that genome-wide gene regulation may utilize a general architecture comprised of multiple semi-autonomous event timelines, whose superposition could produce combinatorial complexity in timing control patterns.

Beneficial mutations and the dynamics of adaptation in asexual populations.

We discuss the dynamics of adaptive evolution in asexual (clonal) populations. The classical 'periodic selection' model of clonal evolution assumed that beneficial mutations are very rare and therefore substitute unfettered into populations as occasional, isolated events. Newer models allow for the possibility that beneficial mutations are sufficiently common to coexist and compete for fixation within populations. Experimental evolution studies in microbes provide empirical support for stochastic models in which both selection and mutation are strong effects and clones compete for fixation; however, the relative importance of competition among clones bearing mutations of different selective effects versus competition among clones bearing multiple mutations remains unresolved. We provide some new theoretical results, moreover, suggesting that population dynamics consistent with the periodic selection model can arise even in a deterministic model that can accommodate a very high beneficial mutation rate.

Identifying mutator phenotypes among fluoroquinolone-resistant strains of Streptococcus pneumoniae using fluctuation analysis.

The occurrence of mutator phenotypes among laboratory-generated and clinical levofloxacin-resistant strains of Streptococcus pneumoniae was determined using fluctuation analysis. The in vitro selection for levofloxacin-resistant mutants of strain D39, each with point mutations in both gyrA and parC or parE, was not associated with a significant change in the mutation rate. Two of eight clinical isolates resistant to levofloxacin (MIC, >8 microg/ml) had estimated mutation rates of 1.2 x 10(-7) and 9.4 x 10(-8) mutations per cell division, indicating potential mutator phenotypes, compared to strain D39, which had an estimated mutation rate of 1.4 x 10(-8) mutations per cell division. The levofloxacin-resistant isolates with the highest mutation rates showed evidence of dysfunctional mismatch repair and contained missense mutations in mut genes at otherwise highly conserved sites. The association of hypermutability in levofloxacin-resistant S. pneumoniae clinical isolates with mutations in DNA mismatch repair genes provides further evidence that mismatch repair mutants may have a selective advantage in the setting of antibiotic pressure, facilitating the development of further antibiotic resistance.

Complete genetic linkage can subvert natural selection.

The intricate adjustment of organisms to their environment demonstrates the effectiveness of natural selection. But Darwin himself recognized that certain biological features could limit this effectiveness, features that generally reduce the efficiency of natural selection or yield suboptimal adaptation. Genetic linkage is known to be one such feature, and here we show theoretically that it can introduce a more sinister flaw: when there is complete linkage between loci affecting fitness and loci affecting mutation rate, positive natural selection and recurrent mutation can drive mutation rates in an adapting population to intolerable levels. We discuss potential implications of this finding for the early establishment of recombination, the evolutionary fate of asexual populations, and immunological clearance of clonal pathogens.

Allopatric divergence, secondary contact, and genetic isolation in wild yeast populations.

In plants and animals, new biological species clearly have arisen as a byproduct of genetic divergence in allopatry. However, our understanding of the processes that generate new microbial species remains limited [1] despite the large contribution of microbes to the world's biodiversity. A recent hypothesis claims that microbes lack biogeographical divergence because their population sizes are large and their migration rates are presumably high [2, 3]. In recapitulating the classic microbial-ecology dictum that "everything is everywhere, and the environment selects"[4, 5], this hypothesis casts doubt on whether geographic divergence promotes speciation in microbes. To date, its predictions have been tested primarily with data from eubacteria and archaebacteria [6-8]. However, this hypothesis's most important implication is in sexual eukaryotic microbes, where migration and genetic admixture are specifically predicted to inhibit allopatric divergence and speciation [9]. Here, we use nuclear-sequence data from globally distributed natural populations of the yeast Saccharomyces paradoxus to investigate the role of geography in generating diversity in sexual eukaryotic microbes. We show that these populations have undergone allopatric divergence and then secondary contact without genetic admixture. Our data thus support the occurrence of evolutionary processes necessary for allopatric speciation in sexual microbes.

Mate choice assays and mating propensity differences in natural yeast populations.

In sexual microbes, mating occurs by fusion of individual cells. This complete fitness investment suggests that cell behaviour could potentially mediate prezygotic isolation between microbial species, a topic about which very little is known. To investigate this possibility, we conducted individual cell mate choice trials and mass-culture mating propensity assays with isolates from sympatric natural populations of the closely related yeasts Saccharomyces cerevisiae and Saccharomyces paradoxus. Although we found no evidence for active species recognition in mate choice, we observed a marked difference in mating propensity between these two species. We briefly discuss the possibility that this mating propensity difference may contribute to reproductive isolation between S. cerevisiae and S. paradoxus in nature.

Population genetic variation in gene expression is associated with phenotypic variation in Saccharomyces cerevisiae.

BACKGROUND: The relationship between genetic variation in gene expression and phenotypic variation observable in nature is not well understood. Identifying how many phenotypes are associated with differences in gene expression and how many gene-expression differences are associated with a phenotype is important to understanding the molecular basis and evolution of complex traits. RESULTS: We compared levels of gene expression among nine natural isolates of Saccharomyces cerevisiae grown either in the presence or absence of copper sulfate. Of the nine strains, two show a reduced growth rate and two others are rust colored in the presence of copper sulfate. We identified 633 genes that show significant differences in expression among strains. Of these genes, 20 were correlated with resistance to copper sulfate and 24 were correlated with rust coloration. The function of these genes in combination with their expression pattern suggests the presence of both correlative and causative expression differences. But the majority of differentially expressed genes were not correlated with either phenotype and showed the same expression pattern both in the presence and absence of copper sulfate. To determine whether these expression differences may contribute to phenotypic variation under other environmental conditions, we examined one phenotype, freeze tolerance, predicted by the differential expression of the aquaporin gene AQY2. We found freeze tolerance is associated with the expression of AQY2. CONCLUSIONS: Gene expression differences provide substantial insight into the molecular basis of naturally occurring traits and can be used to predict environment dependent phenotypic variation.