Ploidy Controls the Success of Mutators and Nature of Mutations during Budding Yeast Evolution

Resolving Deleterious and Near-Neutral Effects Requires Different Pooled Fitness Assay Designs

Pooled sequencing-based fitness assays are a powerful and widely used approach to quantifying fitness of thousands of genetic variants in parallel. Despite the throughput of such assays, they are prone to biases in fitness estimates, and errors in measurements are typically larger for deleterious fitness effects, relative to neutral effects. In practice, designing pooled fitness assays involves tradeoffs between the number of timepoints, the sequencing depth, and other parameters to gain as much information as possible within a feasible experiment. Here, we combined simulations and reanalysis of an existing experimental dataset to explore how assay parameters impact measurements of near-neutral and deleterious fitness effects using a standard fitness estimator. We found that sequencing multiple timepoints at relatively modest depth improved estimates of near-neutral fitness effects, but systematically biased measurements of deleterious effects. We showed that a fixed total number of reads, deeper sequencing at fewer timepoints improved resolution of deleterious fitness effects. Our results highlight a tradeoff between measurement of deleterious and near-neutral effect sizes for a fixed amount of data and suggest that fitness assay design should be tuned for fitness effects that are relevant to the specific biological question.

Identifying Targets of Selection in Laboratory Evolution Experiments

Adaptive evolution navigates a balance between chance and determinism. The stochastic processes of mutation and drift generate phenotypic variation; however, once mutations reach an appreciable frequency in the population, their fate is governed by the deterministic action of selection, enriching for favorable genotypes and purging the less-favorable ones. The net result is that replicate populations will traverse similar-but not identical-pathways to higher fitness. This parallelism in evolutionary outcomes can be leveraged to identify the genes and pathways under selection. However, distinguishing between beneficial and neutral mutations is challenging because many beneficial mutations will be lost due to drift and clonal interference, and many neutral (and even deleterious) mutations will fix by hitchhiking. Here, we review the best practices that our laboratory uses to identify genetic targets of selection from next-generation sequencing data of evolved yeast populations. The general principles for identifying the mutations driving adaptation will apply more broadly.

A fitness trade-off explains the early fate of yeast aneuploids with chromosome gains

The early development of aneuploidy from an accidental chromosome missegregation shows contrasting effects. On the one hand, it is associated with significant cellular stress and decreased fitness. On the other hand, it often carries a beneficial effect and provides a quick (but typically transient) solution to external stress. These apparently controversial trends emerge in several experimental contexts, particularly in the presence of duplicated chromosomes. However, we lack a mathematical evolutionary modeling framework that comprehensively captures these trends from the mutational dynamics and the trade-offs involved in the early stages of aneuploidy. Here, focusing on chromosome gains, we address this point by introducing a fitness model where a fitness cost of chromosome duplications is contrasted by a fitness advantage from the dosage of specific genes. The model successfully captures the experimentally measured probability of emergence of extra chromosomes in a laboratory evolution setup. Additionally, using phenotypic data collected in rich media, we explored the fitness landscape, finding evidence supporting the existence of a per-gene cost of extra chromosomes. Finally, we show that the substitution dynamics of our model, evaluated in the empirical fitness landscape, explains the relative abundance of duplicated chromosomes observed in yeast population genomics data. These findings lay a firm framework for the understanding of the establishment of newly duplicated chromosomes, providing testable quantitative predictions for future observations.

Selecting for infectivity across metapopulations can increase virulence in the social microbe Bacillus thuringiensis

Passage experiments that sequentially infect hosts with parasites have long been used to manipulate virulence. However, for many invertebrate pathogens, passage has been applied naively without a full theoretical understanding of how best to select for increased virulence and this has led to very mixed results. Understanding the evolution of virulence is complex because selection on parasites occurs across multiple spatial scales with potentially different conflicts operating on parasites with different life histories. For example, in social microbes, strong selection on replication rate within hosts can lead to cheating and loss of virulence, because investment in public goods virulence reduces replication rate. In this study, we tested how varying mutation supply and selection for infectivity or pathogen yield (population size in hosts) affected the evolution of virulence against resistant hosts in the specialist insect pathogen Bacillus thuringiensis, aiming to optimize methods for strain improvement against a difficult to kill insect target. We show that selection for infectivity using competition between subpopulations in a metapopulation prevents social cheating, acts to retain key virulence plasmids, and facilitates increased virulence. Increased virulence was associated with reduced efficiency of sporulation, and possible loss of function in putative regulatory genes but not with altered expression of the primary virulence factors. Selection in a metapopulation provides a broadly applicable tool for improving the efficacy of biocontrol agents. Moreover, a structured host population can facilitate artificial selection on infectivity, while selection on life-history traits such as faster replication or larger population sizes can reduce virulence in social microbes.

Oncogenic BRAF induces whole-genome doubling through suppression of cytokinesis

Melanomas and other solid tumors commonly have increased ploidy, with near-tetraploid karyotypes being most frequently observed. Such karyotypes have been shown to arise through whole-genome doubling events that occur during early stages of tumor progression. The generation of tetraploid cells via whole-genome doubling is proposed to allow nascent tumor cells the ability to sample various pro-tumorigenic genomic configurations while avoiding the negative consequences that chromosomal gains or losses have in diploid cells. Whereas a high prevalence of whole-genome doubling events has been established, the means by which whole-genome doubling arises is unclear. Here, we find that BRAFV600E, the most common mutation in melanomas, can induce whole-genome doubling via cytokinesis failure in vitro and in a zebrafish melanoma model. Mechanistically, BRAFV600E causes decreased activation and localization of RhoA, a critical cytokinesis regulator. BRAFV600E activity during G1/S phases of the cell cycle is required to suppress cytokinesis. During G1/S, BRAFV600E activity causes inappropriate centriole amplification, which is linked in part to inhibition of RhoA and suppression of cytokinesis. Together these data suggest that common abnormalities of melanomas linked to tumorigenesis - amplified centrosomes and whole-genome doubling events - can be induced by oncogenic BRAF and other mutations that increase RAS/MAPK pathway activity.

Evolutionary Significance of Fungal Hypermutators: Lessons Learned from Clinical Strains and Implications for Fungal Plant Pathogens

Rapid evolution of fungal pathogens poses a serious threat to medicine and agriculture. The mutation rate determines the pace of evolution of a fungal pathogen. Hypermutator fungal strains have an elevated mutation rate owing to certain defects such as those in the DNA mismatch repair system. Studies in Saccharomyces cerevisiae show that hypermutators expedite evolution by generating beneficial alleles at a faster pace than the wild-type strains. However, an accumulation of deleterious alleles in a hypermutator may reduce its fitness. The balance between fitness cost and mutation benefit determines the prevalence of hypermutators in a population. This balance is affected by a complex interaction of ploidy, mode of reproduction, population size, and recent population history. Studies in human fungal pathogens like Aspergillus fumigatus, Candida albicans, Candida glabrata, Cryptococcus deuterogattii, and Cryptococcus neoformans have highlighted the importance of hypermutators in host adaptation and development of antifungal resistance. However, a critical examination of hypermutator biology, experimental evolution studies, and epidemiological studies suggests that hypermutators may impact evolutionary investigations. This review aims to integrate the knowledge about biology, experimental evolution, and dynamics of fungal hypermutators to critically examine the evolutionary role of hypermutators in fungal pathogen populations and project implications of hypermutators in the evolution of fungal plant pathogen populations. Understanding the factors determining the emergence and evolution of fungal hypermutators can open a novel avenue of managing rapidly evolving fungal pathogens in medicine and agriculture.

Bulk-Fitness Measurements Using Barcode Sequencing Analysis in Yeast

The use of DNA barcodes for determining changes in genotype frequencies has been instrumental to increase the scale at which we can phenotype strain libraries by using next-generation sequencing technologies. Here, we describe the determination of strain fitness for thousands of yeast strains simultaneously in a single assay using recent innovations that increase the precision of these measurements, such as the inclusion of unique-molecular identifiers (UMIs) and purification by solid-phase reverse immobilization (SPRI) beads.

Whole-genome doubling confers unique genetic vulnerabilities on tumour cells

Whole genome doubling (WGD) occurs early in tumorigenesis and generates genetically unstable tetraploid cells that fuel tumor development^1,2. Cells that undergo WGD (WGD⁺) must adapt to accommodate their abnormal tetraploid state; however, the nature of these adaptations, and whether they confer vulnerabilities that can subsequently be exploited therapeutically, is unclear. Using sequencing data from ~10,000 primary human cancer samples and essentiality data from ~600 cancer cell lines, we show that WGD gives rise to common genetic traits that are accompanied by unique vulnerabilities. We reveal that WGD⁺ cells are more dependent on spindle assembly checkpoint signaling, DNA replication factors, and proteasome function than WGD^– cells. We also identify KIF18A, which encodes for a mitotic kinesin, as being specifically required for the viability of WGD⁺ cells. While loss of KIF18A is largely dispensable for accurate chromosome segregation during mitosis in WGD^– cells, its loss induces dramatic mitotic errors in WGD⁺ cells, ultimately impairing cell viability. Collectively, our results reveal new strategies to specifically target WGD⁺ cancer cells while sparing the normal, non-transformed WGD^– cells that comprise human tissue.

The population genetics of ploidy change in unicellular fungi

Changes in ploidy are a significant type of genetic variation, describing the number of chromosome sets per cell. Ploidy evolves in natural populations, clinical populations, and lab experiments, particularly in fungi. Despite a long history of theoretical work on this topic, predicting how ploidy will evolve has proven difficult, as it is often unclear why one ploidy state outperforms another. Here, we review what is known about contemporary ploidy evolution in diverse fungal species through the lens of population genetics. As with typical genetic variants, ploidy evolution depends on the rate that new ploidy states arise by mutation, natural selection on alternative ploidy states, and random genetic drift. However, ploidy variation also has unique impacts on evolution, with the potential to alter chromosomal stability, the rate and patterns of point mutation, and the nature of selection on all loci in the genome. We discuss how ploidy evolution depends on these general and unique factors and highlight areas where additional experimental evidence is required to comprehensively explain the ploidy transitions observed in the field and the lab.

Trajectory and uniqueness of mutational signatures in yeast mutators

The acquisition of mutations plays critical roles in adaptation, evolution, senescence, and tumorigenesis. Massive genome sequencing has allowed extraction of specific features of many mutational landscapes but it remains difficult to retrospectively determine the mechanistic origin(s), selective forces, and trajectories of transient or persistent mutations and genome rearrangements. Here, we conducted a prospective reciprocal approach to inactivate 13 single or multiple evolutionary conserved genes involved in distinct genome maintenance processes and characterize de novo mutations in 274 diploid Saccharomyces cerevisiae mutation accumulation lines. This approach revealed the diversity, complexity, and ultimate uniqueness of mutational landscapes, differently composed of base substitutions, small insertions/deletions (InDels), structural variants, and/or ploidy variations. Several landscapes parallel the repertoire of mutational signatures in human cancers while others are either novel or composites of subsignatures resulting from distinct DNA damage lesions. Notably, the increase of base substitutions in the homologous recombination-deficient Rad51 mutant, specifically dependent on the Polζ translesion polymerase, yields COSMIC signature 3 observed in BRCA1/BRCA2-mutant breast cancer tumors. Furthermore, "mutome" analyses in highly polymorphic diploids and single-cell bottleneck lineages revealed a diverse spectrum of loss-of-heterozygosity (LOH) signatures characterized by interstitial and terminal chromosomal events resulting from interhomolog mitotic cross-overs. Following the appearance of heterozygous mutations, the strong stimulation of LOHs in the rad27/FEN1 and tsa1/PRDX1 backgrounds leads to fixation of homozygous mutations or their loss along the lineage. Overall, these mutomes and their trajectories provide a mechanistic framework to understand the origin and dynamics of genome variations that accumulate during clonal evolution.

Advances in understanding the evolution of fungal genome architecture

Diversity within the fungal kingdom is evident from the wide range of morphologies fungi display as well as the various ecological roles and industrial purposes they serve. Technological advances, particularly in long-read sequencing, coupled with the increasing efficiency and decreasing costs across sequencing platforms have enabled robust characterization of fungal genomes. These sequencing efforts continue to reveal the rampant diversity in fungi at the genome level. Here, we discuss studies that have furthered our understanding of fungal genetic diversity and genomic evolution. These studies revealed the presence of both small-scale and large-scale genomic changes. In fungi, research has recently focused on many small-scale changes, such as how hypermutation and allelic transmission impact genome evolution as well as how and why a few specific genomic regions are more susceptible to rapid evolution than others. High-throughput sequencing of a diverse set of fungal genomes has also illuminated the frequency, mechanisms, and impacts of large-scale changes, which include chromosome structural variation and changes in chromosome number, such as aneuploidy, polyploidy, and the presence of supernumerary chromosomes. The studies discussed herein have provided great insight into how the architecture of the fungal genome varies within species and across the kingdom and how modern fungi may have evolved from the last common fungal ancestor and might also pave the way for understanding how genomic diversity has evolved in all domains of life.

Tetraploidy-Associated Genetic Heterogeneity Confers Chemo-Radiotherapy Resistance to Colorectal Cancer Cells

Tetraploidy, or whole-genome duplication, is a common phenomenon in cancer and preludes chromosome instability, which strongly correlates with disease progression, metastasis, and treatment failure. Therefore, it is reasonable to hypothesize that tetraploidization confers multidrug resistance. Nevertheless, the contribution of whole-genome duplication to chemo-radiotherapy resistance remains unclear. Here, using isogenic diploid and near-tetraploid clones from three colorectal cancer cell lines and one non-transformed human epithelial cell line, we show a consistent growth impairment but a divergent tumorigenic potential of near-tetraploid cells. Next, we assessed the effects of first-line chemotherapeutic drugs, other commonly used agents and ionizing radiation, and found that whole-genome duplication promoted increased chemotherapy resistance and also conferred protection against irradiation. When testing the activation of apoptosis, we observed that tetraploid cells were less prone to caspase 3 activation after treatment with first-line chemotherapeutic agents. Furthermore, we found that pre-treatment with ataxia telangiectasia and Rad3 related (ATR) inhibitors, which targets response to replication stress, significantly enhanced the sensitivity of tetraploid cells to first-line chemotherapeutic agents as well as to ionizing radiation. Our findings provide further insight into how tetraploidy results in greater levels of tolerance to chemo-radiotherapeutic agents and, moreover, we show that ATR inhibitors can sensitize near-tetraploid cells to commonly used chemo-radiotherapy regimens.

Selection on mutators is not frequency-dependent

The evolutionary fate of mutator mutations – genetic variants that raise the genome-wide mutation rate – in asexual populations is often described as being frequency (or number) dependent. Mutators can invade a population by hitchhiking with a sweeping beneficial mutation, but motivated by earlier experiments results, it has been repeatedly suggested that mutators must be sufficiently frequent to produce such a driver mutation before non-mutators do. Here, we use stochastic, agent-based simulations to show that neither the strength nor the sign of selection on mutators depend on their initial frequency, and while the overall probability of hitchhiking increases predictably with frequency, the per-capita probability of fixation remains unchanged.

The fitness cost of mismatch repair mutators in Saccharomyces cerevisiae: partitioning the mutational load

Mutational load is the depression in a population's mean fitness that results from the continual influx of deleterious mutations. Here, we directly estimate the mutational load in a population of haploid Saccharomyces cerevisiae that are deficient for mismatch repair. We partition the load in haploids into two components. To estimate the load due to nonlethal mutations, we measure the competitive fitness of hundreds of randomly selected clones from both mismatch-repair-deficient and -proficient populations. Computation of the mean clone fitness for the mismatch-repair-deficient strain permits an estimation of the nonlethal load, and the histogram of fitness provides an interesting visualization of a loaded population. In a separate experiment, in order to estimate the load due to lethal mutations (i.e. the lethal mutation rate), we manipulate thousands of individual pairs of mother and daughter cells and track their fates. These two approaches yield point estimates for the two contributors to load, and the addition of these estimates is nearly equal to the separately measured short-term competitive fitness deficit for the mismatch-repair-deficient strain. This correspondence suggests that there is no need to invoke direct fitness effects to explain the fitness difference between mismatch-repair-deficient and -proficient strains. Assays in diploids are consistent with deleterious mutations in diploids tending towards recessivity. These results enhance our understanding of mutational load, a central population genetics concept, and we discuss their implications for the evolution of mutation rates.

The Mutator Phenotype: Adapting Microbial Evolution to Cancer Biology

The mutator phenotype hypothesis was postulated almost 40 years ago to reconcile the observation that while cancer cells display widespread mutational burden, acquisition of mutations in non-transformed cells is a rare event. Moreover, it also suggested that cancer evolution could be fostered by increased genome instability. Given the evolutionary conservation throughout the tree of life and the genetic tractability of model organisms, yeast and bacterial species pioneered studies to dissect the functions of genes required for genome maintenance (caretaker genes) or for cell growth control (gatekeeper genes). In this review, we first provide an overview of what we learned from model organisms about the roles of these genes and the genome instability that arises as a consequence of their dysregulation. We then discuss our current understanding of how mutator phenotypes shape the evolution of bacteria and yeast species. We end by bringing clinical evidence that lessons learned from single-cell organisms can be applied to tumor evolution.

Incompatibilities in Mismatch Repair Genes MLH1-PMS1 Contribute to a Wide Range of Mutation Rates in Human Isolates of Baker's Yeast

Laboratory baker's yeast strains bearing an incompatible combination of MLH1 and PMS1 mismatch repair alleles are mutators that can adapt more rapidly to stress, but do so at the cost of long-term fitness. We identified 18 baker's yeast isolates from 1011 surveyed that contain the incompatible MLH1-PMS1 genotype in a heterozygous state. Surprisingly, the incompatible combination from two human clinical heterozygous diploid isolates, YJS5845 and YJS5885, contain the exact MLH1 (S288c-derived) and PMS1 (SK1-derived) open reading frames originally shown to confer incompatibility. While these isolates were nonmutators, their meiotic spore clone progeny displayed mutation rates in a DNA slippage assay that varied over a 340-fold range. This range was 30-fold higher than observed between compatible and incompatible combinations of laboratory strains. Genotyping analysis indicated that MLH1-PMS1 incompatibility was the major driver of mutation rate in the isolates. The variation in the mutation rate of incompatible spore clones could be due to background suppressors and enhancers, as well as aneuploidy seen in the spore clones. Our data are consistent with the observed variance in mutation rate contributing to adaptation to stress conditions (e.g., in a human host) through the acquisition of beneficial mutations, with high mutation rates leading to long-term fitness costs that are buffered by mating or eliminated through natural selection.

Genome doubling shapes the evolution and prognosis of advanced cancers

Ploidy abnormalities are a hallmark of cancer, but their impact on the evolution and outcomes of cancers is unknown. Here, we identified whole-genome doubling (WGD) in the tumors of nearly 30% of 9,692 prospectively sequenced advanced cancer patients. WGD varied by tumor lineage and molecular subtype, and arose early in carcinogenesis after an antecedent transforming driver mutation. While associated with TP53 mutations, 46% of all WGD arose in TP53-wild-type tumors and in such cases was associated with an E2F-mediated G1 arrest defect, although neither aberration was obligate in WGD tumors. The variability of WGD across cancer types can be explained in part by cancer cell proliferation rates. WGD predicted for increased morbidity across cancer types, including KRAS-mutant colorectal cancers and estrogen receptor-positive breast cancers, independently of established clinical prognostic factors. We conclude that WGD is highly common in cancer and is a macro-evolutionary event associated with poor prognosis across cancer types.

Stable aneuploid tumors cells are more sensitive to TTK inhibition than chromosomally unstable cell lines

Inhibition of the spindle assembly checkpoint kinase TTK causes chromosome mis-segregation and tumor cell death. However, high levels of TTK correlate with chromosomal instability (CIN), which can lead to aneuploidy. We show that treatment of tumor cells with the selective small molecule TTK inhibitor NTRC 0066-0 overrides the mitotic checkpoint, irrespective of cell line sensitivity. In stable aneuploid cells NTRC 0066-0 induced acute CIN, whereas in cells with high levels of pre-existing CIN there was only a small additional fraction of cells mis-segregating their chromosomes. In proliferation assays stable aneuploid cells were more sensitive than cell lines with pre-existing CIN. Tetraploids are thought to be an intermediate between diploid and unstable aneuploid cells. TTK inhibitors had the same potency on post-tetraploid and parental diploid cells, which is remarkable because the post-tetraploids are more resistant to mitotic drugs. Finally, we confirm that the reference compound reversine is a TTK inhibitor and like NTRC 0066-0, inhibits the proliferation of patient-derived colorectal cancer organoids. In contrast, treatment with TTK inhibitor did not reduce the viability of non-proliferating T cell acute lymphoblastic leukemia cells samples. Consequently, TTK inhibitor therapy is expected to spare non-dividing cells, and may be used to target stable aneuploid tumors.

Frequent ploidy changes in growing yeast cultures

Ploidy is considered a very stable cellular characteristic. Although rare, changes in ploidy play important roles in the acquisition of long-term adaptations. Since these duplications allow the subsequent loss of individual chromosomes and accumulation of mutations, changes in ploidy can also cause genomic instability, and have been found to promote cancer. Despite the importance of the subject, measuring the rate of whole-genome duplications has proven extremely challenging. We have recently measured the rate of diploidization in yeast using long-term, in-lab experiments. We found that spontaneous diploidization occurs frequently, by two different mechanisms: endoreduplication and mating type switching. Despite its common occurrence, spontaneous diploidization is usually selected against, although it can be advantageous under some stressful conditions. Our results have implications for the understanding of evolutionary processes, as well as for the use of yeast cells in biotechnological applications.

Construction of Comprehensive Dosage-Matching Core Histone Mutant Libraries for Saccharomyces cerevisiae

Saccharomyces cerevisiae contains two genes for each core histone, which are presented as pairs under the control of a divergent promoter, i.e., HHT1-HHF1, HHT2-HHF2, HTA1-HTB1 and HTA2-HTB2HHT1-HHF1, and HHT2-HHF2 encode histone H3 and H4 with identical amino acid sequences but under the control of differently regulated promoters. Previous mutagenesis studies were carried out by deleting one pair and mutating the other one. Here, we present the design and construction of three additional libraries covering HTA1-HTB1, HTA2-HTB2, and HHT1-HHF1 respectively. Together with the previously described library of HHT2-HHF2 mutants, a systematic and complete collection of mutants for each of the eight core S. cerevisiae histone genes becomes available. Each designed mutant was incorporated into the genome, generating three more corresponding libraries of yeast strains. We demonstrated that, although, under normal growth conditions, strains with single-copy integrated histone genes lacked phenotypes, in some growth conditions, growth deficiencies were observed. Specifically, we showed that addition of a second copy of the mutant histone gene could rescue the lethality in some previously known mutants that cannot survive with a single copy. This resource enables systematic studies of function of each nucleosome residue in plasmid, single-copy, and double-copy integrated formats.

Genomics of Adaptation Depends on the Rate of Environmental Change in Experimental Yeast Populations

The rate of directional environmental change may have profound consequences for evolutionary dynamics and outcomes. Yet, most evolution experiments impose a sudden large change in the environment, after which the environment is kept constant. We previously cultured replicate Saccharomyces cerevisiae populations for 500 generations in the presence of either gradually increasing or constant high concentrations of the heavy metals cadmium, nickel, and zinc. Here, we investigate how each of these treatments affected genomic evolution. Whole-genome sequencing of evolved clones revealed that adaptation occurred via a combination of SNPs, small indels, and whole-genome duplications and other large-scale structural changes. In contrast to some theoretical predictions, gradual and abrupt environmental change caused similar numbers of genomic changes. For cadmium, which is toxic already at comparatively low concentrations, mutations in the same genes were used for adaptation to both gradual and abrupt increase in concentration. Conversely, for nickel and zinc, which are toxic at high concentrations only, mutations in different genes were used for adaptation depending on the rate of change. Moreover, evolution was more repeatable following a sudden change in the environment, particularly for nickel and zinc. Our results show that the rate of environmental change and the nature of the selection pressure are important drivers of evolutionary dynamics and outcomes, which has implications for a better understanding of societal problems such as climate change and pollution.

The Influence of Polyploidy on the Evolution of Yeast Grown in a Sub-Optimal Carbon Source

Polyploidization events have occurred during the evolution of many fungi, plant, and animal species and are thought to contribute to speciation and tumorigenesis, however little is known about how ploidy level contributes to adaptation at the molecular level. Here we integrate whole genome sequencing, RNA expression analysis, and relative fitness of ∼100 evolved clones at three ploidy levels. Independent haploid, diploid, and tetraploid populations were grown in a low carbon environment for 250 generations. We demonstrate that the key adaptive mutation in the evolved clones is predicted by a gene expression signature of just five genes. All of the adaptive mutations identified encompass a narrow set of genes, however the tetraploid clones gain a broader spectrum of adaptive mutations than haploid or diploid clones. While many of the adaptive mutations occur in genes that encode proteins with known roles in glucose sensing and transport, we discover mutations in genes with no canonical role in carbon utilization (IPT1 and MOT3), as well as identify novel dominant mutations in glucose signal transducers thought to only accumulate recessive mutations in carbon limited environments (MTH1 and RGT1). We conclude that polyploid cells explore more genotypic and phenotypic space than lower ploidy cells. Our study provides strong evidence for the beneficial role of polyploidization events that occur during the evolution of many species and during tumorigenesis.

Natural mismatch repair mutations mediate phenotypic diversity and drug resistance in Cryptococcus deuterogattii

Pathogenic microbes confront an evolutionary conflict between the pressure to maintain genome stability and the need to adapt to mounting external stresses. Bacteria often respond with elevated mutation rates, but little evidence exists of stable eukaryotic hypermutators in nature. Whole genome resequencing of the human fungal pathogen Cryptococcus deuterogattii identified an outbreak lineage characterized by a nonsense mutation in the mismatch repair component MSH2. This defect results in a moderate mutation rate increase in typical genes, and a larger increase in genes containing homopolymer runs. This allows facile inactivation of genes with coding homopolymer runs including FRR1, which encodes the target of the immunosuppresive antifungal drugs FK506 and rapamycin. Our study identifies a eukaryotic hypermutator lineage spread over two continents and suggests that pathogenic eukaryotic microbes may experience similar selection pressures on mutation rate as bacterial pathogens, particularly during long periods of clonal growth or while expanding into new environments.

As humans, we often think of genetic mutations as being bad. Over the past several decades we have seen health warnings issued on a variety of environmental exposures, from cigarettes to tanning beds, and with good reason because they cause mutations. For multicellular organisms like humans, these mutations are strongly associated with cancer. But in bacteria, this is not true. In fact, the rate at which mutations occur sometimes increases to help bacteria cope with stressful environments.

Unlike bacteria, humans are eukaryotes – the name given to organisms whose cells contain different compartments separated by membranes, such as the nucleus of the cell. For years, we have assumed that eukaryotic microbes, like fungi and parasites, act more like humans than like bacteria because work in budding yeast (another eukaryote) has suggested this to be the case. However, recent work in disease-causing fungi has shown that, much like bacteria, elevated mutation rates may help them to respond to stress. This could also enable fungi to become resistant to drugs used to treat fungal infections.

Cryptococcus deuterogattii is a fungus that causes human diseases including meningoencephalitis and a lung infection called pulmonary cryptococcosis. An ongoing outbreak of the fungus began in the Pacific Northwest of Canada in the late 1990s and emerged in the United States in 2006/2007. Among isolates closely related to those fungi causing the outbreak, three were found that appear to have a specific mutation in their DNA mismatch repair pathway, meaning that they may also experience a higher mutation rate. These strains are also less able to cause disease than others.

Billmyre et al. now demonstrate experimentally that all three isolates have a specific DNA mismatch repair defect, and show that these fungi experience elevated mutation rates, resulting in what is known as a hypermutator state. Furthermore, whole genome sequencing and phylogenetic analysis showed that these hypermutator strains are derived from the outbreak-causing fungi, and that their reduced ability to cause disease is likely a result of accumulating mutations and the loss of the ability to grow at the higher temperatures found in the human body.

Fungal infections are difficult to treat, in part because there are a limited number of available drugs. Elevated mutation rates will likely increase how often and how rapidly fungi develop resistance to these drugs. Understanding how commonly fungi exhibit a hypermutator state that could impact the development of drug resistance will therefore be important for treating patients with fungal infections, which account for millions of infections and hundreds of thousands of deaths annually worldwide.

Adaptive tuning of mutation rates allows fast response to lethal stress in Escherichia coli

While specific mutations allow organisms to adapt to stressful environments, most changes in an organism's DNA negatively impact fitness. The mutation rate is therefore strictly regulated and often considered a slowly-evolving parameter. In contrast, we demonstrate an unexpected flexibility in cellular mutation rates as a response to changes in selective pressure. We show that hypermutation independently evolves when different Escherichia coli cultures adapt to high ethanol stress. Furthermore, hypermutator states are transitory and repeatedly alternate with decreases in mutation rate. Specifically, population mutation rates rise when cells experience higher stress and decline again once cells are adapted. Interestingly, we identified cellular mortality as the major force driving the quick evolution of mutation rates. Together, these findings show how organisms balance robustness and evolvability and help explain the prevalence of hypermutation in various settings, ranging from emergence of antibiotic resistance in microbes to cancer relapses upon chemotherapy.

Known mutator alleles do not markedly increase mutation rate in clinical Saccharomyces cerevisiae strains

Natural selection has the potential to act on all phenotypes, including genomic mutation rate. Classic evolutionary theory predicts that in asexual populations, mutator alleles, which cause high mutation rates, can fix due to linkage with beneficial mutations. This phenomenon has been demonstrated experimentally and may explain the frequency of mutators found in bacterial pathogens. By contrast, in sexual populations, recombination decouples mutator alleles from beneficial mutations, preventing mutator fixation. In the facultatively sexual yeast Saccharomyces cerevisiae, segregating alleles of MLH1 and PMS1 have been shown to be incompatible, causing a high mutation rate when combined. These alleles had never been found together naturally, but were recently discovered in a cluster of clinical isolates. Here we report that the incompatible mutator allele combination only marginally elevates mutation rate in these clinical strains. Genomic and phylogenetic analyses provide no evidence of a historically elevated mutation rate. We conclude that the effect of the mutator alleles is dampened by background genetic modifiers. Thus, the relationship between mutation rate and microbial pathogenicity may be more complex than once thought. Our findings provide rare observational evidence that supports evolutionary theory suggesting that sexual organisms are unlikely to harbour alleles that increase their genomic mutation rate.

Mismatch Repair Incompatibilities in Diverse Yeast Populations

An elevated mutation rate can provide cells with a source of mutations to adapt to changing environments. We identified a negative epistatic interaction involving naturally occurring variants in the MLH1 and PMS1 mismatch repair (MMR) genes of Saccharomyces cerevisiae We hypothesized that this MMR incompatibility, created through mating between divergent S. cerevisiae, yields mutator progeny that can rapidly but transiently adapt to an environmental stress. Here we analyzed the MLH1 and PMS1 genes across 1010 S. cerevisiae natural isolates spanning a wide range of ecological sources (tree exudates, Drosophila, fruits, and various fermentation and clinical isolates) and geographical sources (Europe, America, Africa, and Asia). We identified one homozygous clinical isolate and 18 heterozygous isolates containing the incompatible MMR genotype. The MLH1-PMS1 gene combination isolated from the homozygous clinical isolate conferred a mutator phenotype when expressed in the S288c laboratory background. Using a novel reporter to measure mutation rates, we showed that the overall mutation rate in the homozygous incompatible background was similar to that seen in compatible strains, indicating the presence of suppressor mutations in the clinical isolate that lowered its mutation rate. This observation and the identification of 18 heterozygous isolates, which can lead to MMR incompatible genotypes in the offspring, are consistent with an elevated mutation rate rapidly but transiently facilitating adaptation. To avoid long-term fitness costs, the incompatibility is apparently buffered by mating or by acquiring suppressors. These observations highlight effective strategies in eukaryotes to avoid long-term fitness costs associated with elevated mutation rates.

Interactions between Genetic and Ecological Effects on the Evolution of Life Cycles

Sexual reproduction leads to an alternation between haploid and diploid phases, whose relative length varies widely across taxa. Previous genetical models showed that diploid or haploid life cycles may be favored, depending on dominance interactions and on effective recombination rates. By contrast, niche differentiation between haploids and diploids may favor biphasic life cycles, in which development occurs in both phases. In this article, we explore the interplay between genetical and ecological factors, assuming that deleterious mutations affect the competitivity of individuals within their ecological niche and allowing different effects of mutations in haploids and diploids (including antagonistic selection). We show that selection on a modifier gene affecting the relative length of both phases can be decomposed into a direct selection term favoring the phase with the highest mean fitness (due to either ecological differences or differential effects of mutations) and an indirect selection term favoring the phase in which selection is more efficient. When deleterious alleles occur at many loci and in the presence of ecological differentiation between haploids and diploids, evolutionary branching often occurs and leads to the stable coexistence of alleles coding for haploid and diploid cycles, while temporal variations in niche sizes may stabilize biphasic cycles.

Evolution of Mutation Rates in Rapidly Adapting Asexual Populations

Mutator and antimutator alleles often arise and spread in both natural microbial populations and laboratory evolution experiments. The evolutionary dynamics of these mutation rate modifiers are determined by indirect selection on linked beneficial and deleterious mutations. These indirect selection pressures have been the focus of much earlier theoretical and empirical work, but we still have a limited analytical understanding of how the interplay between hitchhiking and deleterious load influences the fates of modifier alleles. Our understanding is particularly limited when clonal interference is common, which is the regime of primary interest in laboratory microbial evolution experiments. Here, we calculate the fixation probability of a mutator or antimutator allele in a rapidly adapting asexual population, and we show how this quantity depends on the population size, the beneficial and deleterious mutation rates, and the strength of a typical driver mutation. In the absence of deleterious mutations, we find that clonal interference enhances the fixation probability of mutators, even as they provide a diminishing benefit to the overall rate of adaptation. When deleterious mutations are included, natural selection pushes the population toward a stable mutation rate that can be suboptimal for the adaptation of the population as a whole. The approach to this stable mutation rate is not necessarily monotonic: even in the absence of epistasis, selection can favor mutator and antimutator alleles that "overshoot" the stable mutation rate by substantial amounts.

Transcriptome analysis of tetraploid cells identifies cyclin D2 as a facilitator of adaptation to genome doubling in the presence of p53

Tetraploidization, or genome doubling, is a prominent event in tumorigenesis, primarily because cell division in polyploid cells is error-prone and produces aneuploid cells. This study investigates changes in gene expression evoked in acute and adapted tetraploid cells and their effect on cell-cycle progression. Acute polyploidy was generated by knockdown of the essential regulator of cytokinesis anillin, which resulted in cytokinesis failure and formation of binucleate cells, or by chemical inhibition of Aurora kinases, causing abnormal mitotic exit with formation of single cells with aberrant nuclear morphology. Transcriptome analysis of these acute tetraploid cells revealed common signatures of activation of the tumor-suppressor protein p53. Suppression of proliferation in these cells was dependent on p53 and its transcriptional target, CDK inhibitor p21. Rare proliferating tetraploid cells can emerge from acute polyploid populations. Gene expression analysis of single cell-derived, adapted tetraploid clones showed up-regulation of several p53 target genes and cyclin D2, the activator of CDK4/6/2. Overexpression of cyclin D2 in diploid cells strongly potentiated the ability to proliferate with increased DNA content despite the presence of functional p53. These results indicate that p53-mediated suppression of proliferation of polyploid cells can be averted by increased levels of oncogenes such as cyclin D2, elucidating a possible route for tetraploidy-mediated genomic instability in carcinogenesis.

Experimental evolution in fungi: An untapped resource

Historically, evolutionary biology has been considered an observational science. Examining populations and inferring evolutionary histories mold evolutionary theories. In contrast, laboratory evolution experiments make use of the amenability of traditional model organisms to study fundamental processes underlying evolution in real time in simple, but well-controlled, environments. With advances in high-throughput biology and next generation sequencing, it is now possible to propagate hundreds of parallel populations over thousands of generations and to quantify precisely the frequencies of various mutations over time. Experimental evolution combines the ability to simultaneously monitor replicate populations with the power to vary individual parameters to test specific evolutionary hypotheses, something that is impractical or infeasible in natural populations. Many labs are now conducting laboratory evolution experiments in nearly all model systems including viruses, bacteria, yeast, nematodes, and fruit flies. Among these systems, fungi occupy a unique niche: with a short generation time, small compact genomes, and sexual cycles, fungi are a particularly valuable and largely untapped resource for propelling future growth in the field of experimental evolution. Here, we describe the current state of fungal experimental evolution and why fungi are uniquely positioned to answer many of the outstanding questions in the field. We also review which fungal species are most well suited for experimental evolution.

Genomic confirmation of nutrient-dependent mutability of mutators in Escherichia coli

Mutators with increased mutation rates are prevalent in various environments and have important roles in accelerating adaptive evolution. Previous studies on mutator strains of microorganisms have shown that some mutators have constant mutation rates, whereas others exhibit switchable mutation rates depending on nutritional conditions. This suggests that the contributions of mutators on evolution vary with fluctuating nutritional conditions. However, such conditional mutability has been unclear at the genomic level. In addition, it is still unknown why mutation rates change with nutritional condition. Here, we used two mutator strains of Escherichia coli to explore the nutrient dependence of mutation rates at the genomic level. These strains were transferred repeatedly under different nutritional conditions for hundreds of generations to accumulate mutations. Whole-genome sequencing of the offspring showed that the nutrient dependence of the mutation rates was pervasive at the genomic scale. Neutrality in the mutation accumulation processes and constancy in the mutational bias suggested that nutrient dependence was not derived from conditional selective purges or from shifts of mutational bias. Some mutators could simply switch their mutation rates for both transitions and transversions in response to nutritional shifts.

A Genetic Incompatibility Accelerates Adaptation in Yeast

During mismatch repair (MMR) MSH proteins bind to mismatches that form as the result of DNA replication errors and recruit MLH factors such as Mlh1-Pms1 to initiate excision and repair steps. Previously, we identified a negative epistatic interaction involving naturally occurring polymorphisms in the MLH1 and PMS1 genes of baker’s yeast. Here we hypothesize that a mutagenic state resulting from this negative epistatic interaction increases the likelihood of obtaining beneficial mutations that can promote adaptation to stress conditions. We tested this by stressing yeast strains bearing mutagenic (incompatible) and non-mutagenic (compatible) mismatch repair genotypes. Our data show that incompatible populations adapted more rapidly and without an apparent fitness cost to high salt stress. The fitness advantage of incompatible populations was rapid but disappeared over time. The fitness gains in both compatible and incompatible strains were due primarily to mutations in PMR1 that appeared earlier in incompatible evolving populations. These data demonstrate a rapid and reversible role (by mating) for genetic incompatibilities in accelerating adaptation in eukaryotes. They also provide an approach to link experimental studies to observational population genomics.

In nature, bacterial populations with high mutation rates can adapt faster to new environments by acquiring beneficial mutations. However, such populations also accumulate harmful mutations that reduce their fitness. We show that the model eukaryote baker’s yeast can use a similar mutator strategy to adapt to new environments. The mutator state that we observed resulted from an incompatibility involving two genes, MLH1 and PMS1, that work together to remove DNA replication errors through a spellchecking mismatch repair mechanism. This incompatibility can occur through mating between baker’s yeast from different genetic backgrounds, yielding mutator offspring containing an MLH1-PMS1 combination not present in either parent. Interestingly, these offspring adapted more rapidly to stress, compared to the parental strains, and did so without an overall loss in fitness. DNA sequencing analyses of baker’s yeast strains from across the globe support the presence of incompatible hybrid yeast strains in nature. These observations provide a powerful model to understand how the segregation of defects in DNA mismatch repair can serve as an effective strategy to enable eukaryotes to adapt to changing environments.

Real-Time Evolution of a Subtelomeric Gene Family in Candida albicans

Subtelomeric regions of the genome are notable for high rates of sequence evolution and rapid gene turnover. Evidence of subtelomeric evolution has relied heavily on comparisons of historical evolutionary patterns to infer trends and frequencies of these events. Here, we describe evolution of the subtelomeric TLO gene family in Candida albicans during laboratory passaging for over 4000 generations. C. albicans is a commensal and opportunistic pathogen of humans and the TLO gene family encodes a subunit of the Mediator complex that regulates transcription and affects a range of virulence factors. We identified 16 distinct subtelomeric recombination events that altered the TLO repertoire. Ectopic recombination between subtelomeres on different chromosome ends occurred approximately once per 5000 generations and was often followed by loss of heterozygosity, resulting in the complete loss of one TLO gene sequence with expansion of another. In one case, recombination within TLO genes produced a novel TLO gene sequence. TLO copy number changes were biased, with some TLOs preferentially being copied to novel chromosome arms and other TLO genes being frequently lost. The majority of these nonreciprocal recombination events occurred either within the 3′ end of the TLO coding sequence or within a conserved 50-bp sequence element centromere-proximal to TLO coding sequence. Thus, subtelomeric recombination is a rapid mechanism of generating genotypic diversity through alterations in the number and sequence of related gene family members.

dNTP pool levels modulate mutator phenotypes of error-prone DNA polymerase ε variants

Mutator phenotypes create genetic diversity that fuels tumor evolution. DNA polymerase (Pol) ε mediates leading strand DNA replication. Proofreading defects in this enzyme drive a number of human malignancies. Here, using budding yeast, we show that mutator variants of Pol ε depend on damage uninducible (Dun)1, an S-phase checkpoint kinase that maintains dNTP levels during a normal cell cycle and up-regulates dNTP synthesis upon checkpoint activation. Deletion of DUN1 (dun1Δ) suppresses the mutator phenotype of pol2-4 (encoding Pol ε proofreading deficiency) and is synthetically lethal with pol2-M644G (encoding altered Pol ε base selectivity). Although pol2-4 cells cycle normally, pol2-M644G cells progress slowly through S-phase. The pol2-M644G cells tolerate deletions of mediator of the replication checkpoint (MRC) 1 (mrc1Δ) and radiation sensitive (Rad) 9 (rad9Δ), which encode mediators of checkpoint responses to replication stress and DNA damage, respectively. The pol2-M644G mutator phenotype is partially suppressed by mrc1Δ but not rad9Δ; neither deletion suppresses the pol2-4 mutator phenotype. Thus, checkpoint activation augments the Dun1 effect on replication fidelity but is not required for it. Deletions of genes encoding key Dun1 targets that negatively regulate dNTP synthesis, suppress the dun1Δ pol2-M644G synthetic lethality and restore the mutator phenotype of pol2-4 in dun1Δ cells. DUN1 pol2-M644G cells have constitutively high dNTP levels, consistent with checkpoint activation. In contrast, pol2-4 and POL2 cells have similar dNTP levels, which decline in the absence of Dun1 and rise in the absence of the negative regulators of dNTP synthesis. Thus, dNTP pool levels correlate with Pol ε mutator severity, suggesting that treatments targeting dNTP pools could modulate mutator phenotypes for therapy.

The evolution of drug resistance in clinical isolates of Candida albicans

Candida albicans is both a member of the healthy human microbiome and a major pathogen in immunocompromised individuals. Infections are typically treated with azole inhibitors of ergosterol biosynthesis often leading to drug resistance. Studies in clinical isolates have implicated multiple mechanisms in resistance, but have focused on large-scale aberrations or candidate genes, and do not comprehensively chart the genetic basis of adaptation. Here, we leveraged next-generation sequencing to analyze 43 isolates from 11 oral candidiasis patients. We detected newly selected mutations, including single-nucleotide polymorphisms (SNPs), copy-number variations and loss-of-heterozygosity (LOH) events. LOH events were commonly associated with acquired resistance, and SNPs in 240 genes may be related to host adaptation. Conversely, most aneuploidies were transient and did not correlate with drug resistance. Our analysis also shows that isolates also varied in adherence, filamentation, and virulence. Our work reveals new molecular mechanisms underlying the evolution of drug resistance and host adaptation.

DOI:http://dx.doi.org/10.7554/eLife.00662.001

Nearly all humans are infected with the fungus Candida albicans. In most people, the infection does not produce any symptoms because their immune system is able to counteract the fungus' attempts to spread around the body. However, if the balance between fungal attack and body defence fails, the fungus is able to spread, which can lead to serious disease that is fatal in 42% of cases.

How does C. albicans outcompete the body's defences to cause disease? This is a pertinent question because the most effective antifungal medicines—including the drug fluconazole—do not kill the fungus; they only stop it from growing. This gives the fungus time to develop resistance to the drug by becoming able to quickly replace the fungal proteins the drug destroys, or to efficiently remove the drug from its cells.

In this study, Ford et al. studied the changes that occur in the DNA of C. albicans over time in patients who are being treated with fluconazole. Ford et al. took 43 samples of C. albicans from 11 patients with weakened immune systems. The experiments show that the fungus samples collected early on were more sensitive to the drug than the samples collected later.

In most cases, the genetic data suggest that the infections begin with a single fungal cell; the cells in the later samples are its offspring. Despite this, there is a lot of genetic variation between samples from the same patient, which indicates that the fungus is under pressure to become more resistant to the drug. There were 240 genes—including those that can alter the surface on the fungus cells to make it better at evading the host immune system—in which small changes occurred over time in three or more patients. Laboratory tests revealed that many of these genes are likely important for the fungus to survive in an animal host in the presence of the drug.

C. albicans cells usually have two genetically distinct copies of every gene. Ford et al. found that for some genes—including some that make surface components or are involved in expelling drugs from cells—the loss of genetic information from one copy, so that both copies become identical, is linked to resistance to fluconazole. However, the gain of whole or partial chromosomes—which contain large numbers of genes—is not linked to resistance, but may provide additional genetic material for generating diversity in the yeast population that may help the cells to evolve resistance in the future.

These experiments have identified many new candidate genes that are important for drug resistance and evading the host immune system, and which could be used to guide the development of new therapeutics to treat these life-threatening infections.

DOI:http://dx.doi.org/10.7554/eLife.00662.002

Nutrient-dependent growth defects and mutability of mutators in Escherichia coli

So-called mutators emerge when mismatch repair and proofreading mechanisms are defective. Mutators not only accelerate the accumulation of mutations that are beneficial for adaptation but also cause a large number of deleterious mutations that are disadvantageous for cell growth. However, such growth defects may be compensated by nutrient availability. How the growth burden is associated with high mutability in relation to nutritional variation is an intriguing question. To address this question, we constructed a variety of Escherichia coli mutator strains through combinatorial deletions of mismatch repair and proofreading genes and quantitatively evaluated their growth and mutation rates under different nutritional conditions. Growth defects caused by high mutation rates were commonly observed in all mutators, and these defects were alleviated by nutrient supplementation in most mutators. In addition, the mutation rates of the mutators fluctuated greatly in response to nutritional conditions, in contrast to the nearly constant mutation rate of the wild-type strain under varying nutritional conditions. The results showed conditional growth defects and nutrition-sensitive mutability as general features of mutators. This study indicates the importance of modulating mutability in response to changing nutrient conditions to minimize the risk of extinction due to genetic load.

Cancer: evolution within a lifetime

Subclonal cancer populations change spatially and temporally during the disease course. Studies are revealing branched evolutionary cancer growth with low-frequency driver events present in subpopulations of cells, providing escape mechanisms for targeted therapeutic approaches. Despite such complexity, evidence is emerging for parallel evolution of subclones, mediated through distinct somatic events converging on the same gene, signal transduction pathway, or protein complex in different subclones within the same tumor. Tumors may follow gradualist paths (microevolution) as well as major shifts in evolutionary trajectories (macroevolution). Although macroevolution has been subject to considerable controversy in post-Darwinian evolutionary theory, we review evidence that such nongradual, saltatory leaps, driven through chromosomal rearrangements or genome doubling, may be particularly relevant to tumor evolution. Adapting cancer care to the challenges imposed by tumor micro- and macroevolution and developing deeper insight into parallel evolutionary events may prove central to improving outcome and reducing drug development costs.

Size Does Matter: Why Polyploid Tumor Cells are Critical Drug Targets in the War on Cancer

Tumor evolution presents a formidable obstacle that currently prevents the development of truly curative treatments for cancer. In this perspective, we advocate for the hypothesis that tumor cells with significantly elevated genomic content (polyploid tumor cells) facilitate rapid tumor evolution and the acquisition of therapy resistance in multiple incurable cancers. We appeal to studies conducted in yeast, cancer models, and cancer patients, which all converge on the hypothesis that polyploidy enables large phenotypic leaps, providing access to many different therapy-resistant phenotypes. We develop a flow-cytometry based method for quantifying the prevalence of polyploid tumor cells, and show the frequency of these cells in patient tumors may be higher than is generally appreciated. We then present recent studies identifying promising new therapeutic strategies that could be used to specifically target polyploid tumor cells in cancer patients. We argue that these therapeutic approaches should be incorporated into new treatment strategies aimed at blocking tumor evolution by killing the highly evolvable, therapy-resistant polyploid cell subpopulations, thus helping to maintain patient tumors in a drug sensitive state.

Experimental evolution and the dynamics of genomic mutation rate modifiers

Because genes that affect mutation rates are themselves subject to mutation, mutation rates can be influenced by natural selection and other evolutionary forces. The population genetics of mutation rate modifier alleles has been a subject of theoretical interest for many decades. Here, we review experimental contributions to our understanding of mutation rate modifier dynamics. Numerous evolution experiments have shown that mutator alleles (modifiers that elevate the genomic mutation rate) can readily rise to high frequencies via genetic hitchhiking in non-recombining microbial populations. Whereas these results certainly provide an explanatory framework for observations of sporadically high mutation rates in pathogenic microbes and in cancer lineages, it is nonetheless true that most natural populations have very low mutation rates. This raises the interesting question of how mutator hitchhiking is suppressed or its phenotypic effect reversed in natural populations. Very little experimental work has addressed this question; with this in mind, we identify some promising areas for future experimental investigation.

DNA replication error-induced extinction of diploid yeast

Genetic defects in DNA polymerase accuracy, proofreading, or mismatch repair (MMR) induce mutator phenotypes that accelerate adaptation of microbes and tumor cells. Certain combinations of mutator alleles synergistically increase mutation rates to levels that drive extinction of haploid cells. The maximum tolerated mutation rate of diploid cells is unknown. Here, we define the threshold for replication error-induced extinction (EEX) of diploid Saccharomyces cerevisiae. Double-mutant pol3 alleles that carry mutations for defective DNA polymerase-δ proofreading (pol3-01) and accuracy (pol3-L612M or pol3-L612G) induce strong mutator phenotypes in heterozygous diploids (POL3/pol3-01,L612M or POL3/pol3-01,L612G). Both pol3-01,L612M and pol3-01,L612G alleles are lethal in the homozygous state; cells with pol3-01,L612M divide up to 10 times before arresting at random stages in the cell cycle. Antimutator eex mutations in the pol3 alleles suppress this lethality (pol3-01,L612M,eex or pol3-01,L612G,eex). MMR defects synergize with pol3-01,L612M,eex and pol3-01,L612G,eex alleles, increasing mutation rates and impairing growth. Conversely, inactivation of the Dun1 S-phase checkpoint kinase suppresses strong pol3-01,L612M,eex and pol3-01,L612G,eex mutator phenotypes as well as the lethal pol3-01,L612M phenotype. Our results reveal that the lethal error threshold in diploids is 10 times higher than in haploids and likely determined by homozygous inactivation of essential genes. Pronounced loss of fitness occurs at mutation rates well below the lethal threshold, suggesting that mutator-driven cancers may be susceptible to drugs that exacerbate replication errors.

The effect of population bottlenecks on mutation rate evolution in asexual populations

In the absence of recombination, a mutator allele can spread through a population by hitchhiking with beneficial mutations that appear in its genetic background. Theoretical studies over the past decade have shown that the survival and fixation probability of beneficial mutations can be severely reduced by population size bottlenecks. Here, we use computational modelling and evolution experiments with the yeast S. cerevisiae to examine whether population bottlenecks can affect mutator dynamics in adapting asexual populations. In simulation, we show that population bottlenecks can inhibit mutator hitchhiking with beneficial mutations and are most effective at lower beneficial mutation supply rates. We then subjected experimental populations of yeast propagated at the same effective population size to three different bottleneck regimes and observed that the speed of mutator hitchhiking was significantly slower at smaller bottlenecks, consistent with our theoretical expectations. Our results, thus, suggest that bottlenecks can be an important factor in mutation rate evolution and can in certain circumstances act to stabilize or, at least, delay the progressive elevation of mutation rates in asexual populations. Additionally, our findings provide the first experimental support for the theoretically postulated effect of population bottlenecks on beneficial mutations and demonstrate the usefulness of studying mutator frequency dynamics for understanding the underlying dynamics of fitness-affecting mutations.

Genome evolution due to allopolyploidization in wheat

The wheat group has evolved through allopolyploidization, namely, through hybridization among species from the plant genera Aegilops and Triticum followed by genome doubling. This speciation process has been associated with ecogeographical expansion and with domestication. In the past few decades, we have searched for explanations for this impressive success. Our studies attempted to probe the bases for the wide genetic variation characterizing these species, which accounts for their great adaptability and colonizing ability. Central to our work was the investigation of how allopolyploidization alters genome structure and expression. We found in wheat that allopolyploidy accelerated genome evolution in two ways: (1) it triggered rapid genome alterations through the instantaneous generation of a variety of cardinal genetic and epigenetic changes (which we termed "revolutionary" changes), and (2) it facilitated sporadic genomic changes throughout the species' evolution (i.e., evolutionary changes), which are not attainable at the diploid level. Our major findings in natural and synthetic allopolyploid wheat indicate that these alterations have led to the cytological and genetic diploidization of the allopolyploids. These genetic and epigenetic changes reflect the dynamic structural and functional plasticity of the allopolyploid wheat genome. The significance of this plasticity for the successful establishment of wheat allopolyploids, in nature and under domestication, is discussed.

Emergence of DNA polymerase ε antimutators that escape error-induced extinction in yeast

DNA polymerases (Pols) ε and δ perform the bulk of yeast leading- and lagging-strand DNA synthesis. Both Pols possess intrinsic proofreading exonucleases that edit errors during polymerization. Rare errors that elude proofreading are extended into duplex DNA and excised by the mismatch repair (MMR) system. Strains that lack Pol proofreading or MMR exhibit a 10- to 100-fold increase in spontaneous mutation rate (mutator phenotype), and inactivation of both Pol δ proofreading (pol3-01) and MMR is lethal due to replication error-induced extinction (EEX). It is unclear whether a similar synthetic lethal relationship exists between defects in Pol ε proofreading (pol2-4) and MMR. Using a plasmid-shuffling strategy in haploid Saccharomyces cerevisiae, we observed synthetic lethality of pol2-4 with alleles that completely abrogate MMR (msh2Δ, mlh1Δ, msh3Δ msh6Δ, or pms1Δ mlh3Δ) but not with partial MMR loss (msh3Δ, msh6Δ, pms1Δ, or mlh3Δ), indicating that high levels of unrepaired Pol ε errors drive extinction. However, variants that escape this error-induced extinction (eex mutants) frequently emerged. Five percent of pol2-4 msh2Δ eex mutants encoded second-site changes in Pol ε that reduced the pol2-4 mutator phenotype between 3- and 23-fold. The remaining eex alleles were extragenic to pol2-4. The locations of antimutator amino-acid changes in Pol ε and their effects on mutation spectra suggest multiple mechanisms of mutator suppression. Our data indicate that unrepaired leading- and lagging-strand polymerase errors drive extinction within a few cell divisions and suggest that there are polymerase-specific pathways of mutator suppression. The prevalence of suppressors extragenic to the Pol ε gene suggests that factors in addition to proofreading and MMR influence leading-strand DNA replication fidelity.

The yield of experimental yeast populations declines during selection

The trade-off between growth rate and yield can limit population productivity. Here we tested for this life-history trade-off in replicate haploid and diploid populations of Saccharomyces cerevisiae propagated in glucose-limited medium in batch cultures for 5000 generations. The yield of single clones isolated from the haploid lineages, measured as both optical and population density at the end of a growth cycle, declined during selection and was negatively correlated with growth rate. Initially, diploid populations did not pay this cost of adaptation but haploidized after about 1000-3000 generations of selection, and this ploidy transition was associated with a decline in yield caused by reduced cell size. These results demonstrate the experimental evolution of a trade-off between growth rate and yield, caused by antagonistic pleiotropy, during adaptation in haploids and after an adaptive transition from diploidy to haploidy.

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.

Hot spots for allosteric regulation on protein surfaces

Recent work indicates a general architecture for proteins in which sparse networks of physically contiguous and coevolving amino acids underlie basic aspects of structure and function. These networks, termed sectors, are spatially organized such that active sites are linked to many surface sites distributed throughout the structure. Using the metabolic enzyme dihydrofolate reductase as a model system, we show that: (1) the sector is strongly correlated to a network of residues undergoing millisecond conformational fluctuations associated with enzyme catalysis, and (2) sector-connected surface sites are statistically preferred locations for the emergence of allosteric control in vivo. Thus, sectors represent an evolutionarily conserved "wiring" mechanism that can enable perturbations at specific surface positions to rapidly initiate conformational control over protein function. These findings suggest that sectors enable the evolution of intermolecular communication and regulation.

Genome-wide analytical approaches for reverse metabolic engineering of industrially relevant phenotypes in yeast

Successful reverse engineering of mutants that have been obtained by nontargeted strain improvement has long presented a major challenge in yeast biotechnology. This paper reviews the use of genome-wide approaches for analysis of Saccharomyces cerevisiae strains originating from evolutionary engineering or random mutagenesis. On the basis of an evaluation of the strengths and weaknesses of different methods, we conclude that for the initial identification of relevant genetic changes, whole genome sequencing is superior to other analytical techniques, such as transcriptome, metabolome, proteome, or array-based genome analysis. Key advantages of this technique over gene expression analysis include the independency of genome sequences on experimental context and the possibility to directly and precisely reproduce the identified changes in naive strains. The predictive value of genome-wide analysis of strains with industrially relevant characteristics can be further improved by classical genetics or simultaneous analysis of strains derived from parallel, independent strain improvement lineages.

Mutator suppression and escape from replication error-induced extinction in yeast

Cells rely on a network of conserved pathways to govern DNA replication fidelity. Loss of polymerase proofreading or mismatch repair elevates spontaneous mutation and facilitates cellular adaptation. However, double mutants are inviable, suggesting that extreme mutation rates exceed an error threshold. Here we combine alleles that affect DNA polymerase δ (Pol δ) proofreading and mismatch repair to define the maximal error rate in haploid yeast and to characterize genetic suppressors of mutator phenotypes. We show that populations tolerate mutation rates 1,000-fold above wild-type levels but collapse when the rate exceeds 10⁻³ inactivating mutations per gene per cell division. Variants that escape this error-induced extinction (eex) rapidly emerge from mutator clones. One-third of the escape mutants result from second-site changes in Pol δ that suppress the proofreading-deficient phenotype, while two-thirds are extragenic. The structural locations of the Pol δ changes suggest multiple antimutator mechanisms. Our studies reveal the transient nature of eukaryotic mutators and show that mutator phenotypes are readily suppressed by genetic adaptation. This has implications for the role of mutator phenotypes in cancer.

Organisms strike a balance between genetic continuity and change. Most cells are well adapted to their niches and therefore invest heavily in mechanisms that maintain accurate DNA replication. When cell populations are confronted with changing environmental conditions, “mutator” clones with high mutation rates emerge and readily adapt to the new conditions by rapidly acquiring beneficial mutations. However, deleterious mutations also accumulate, raising the question: what level of mutational burden can cell populations sustain before collapsing? Here we experimentally determine the maximal mutation rate in haploid yeast. We observe that yeast can withstand a 1,000-fold increase in mutation rate without losing colony forming capacity. Yet no strains survive a 10,000-fold increase in mutation rate. Escape mutants with an “anti-mutator” phenotype frequently emerge from cell populations undergoing this error-induced extinction. The diversity of antimutator changes suggests that strong mutator phenotypes in nature may be inherently transient, ensuring that rapid adaptation is followed by genetic attenuation which preserves the beneficial, adaptive mutations. These observations are relevant to microbial populations during infection as well as the somatic evolution of cancer cells.