Precise estimates of mutation rate and spectrum in yeast

Instability throughout the Saccharomyces cerevisiae genome resulting from Pms1 endonuclease deficiency

The endonuclease activity of Pms1 directs mismatch repair by generating a nick in the newly replicated DNA strand. Inactivating Pms2, the human homologue of yeast Pms1, increases the chances of colorectal and uterine cancers. Here we use whole genome sequencing to show that loss of this endonuclease activity, via the pms1-DE variant, results in strong mutator effects throughout the Saccharomyces cerevisiae genome. Mutation rates are strongly increased for mutations resulting from all types of single-base substitutions and for a wide variety of single- and multi-base indel mutations. Rates for these events are further increased in strains combining pms1-DE with mutator variants of each of the three major leading and lagging strand replicases. In all cases, mutation rates, spectra, biases, and context preferences are statistically indistinguishable from strains with equivalent polymerases but lacking initial mismatch recognition due to deletion of MSH2. This implies that, across the nuclear genome, strand discrimination via the Pms1 endonuclease is as important for MMR as is initial mismatch recognition by Msh2 heterodimers.

Irc20 modulates LOH frequency and distribution in S. cerevisiae

Loss of Heterozygosity (LOH) due to mitotic recombination is frequently associated with the development of various cancers (e.g. retinoblastoma). LOH is also an important source of genetic diversity, especially in organisms where meiosis is infrequent. Irc20 is a putative helicase, and E3 ubiquitin ligase involved in DNA double-strand break repair pathway. We analyzed genome-wide LOH events, gross chromosomal changes, small insertion-deletions and single nucleotide mutations in eleven S. cerevisiae mutation accumulation lines of irc20∆, which underwent 50 mitotic bottlenecks. LOH enhancement in irc20∆ was small (1.6 fold), but statistically significant as compared to the wild type. Short (≤ 1 kb) and long (> 10 kb) LOH tracts were significantly enhanced in irc20∆. Both interstitial and terminal LOH events were also significantly enhanced in irc20∆ compared to the wild type. LOH events in irc20∆ were more telomere proximal and away from centromeres compared to the wild type. Gross chromosomal changes, single nucleotide mutations and in-dels were comparable between irc20∆ and wild type. Locus based and genome-wide analysis of meiotic recombination showed that meiotic crossover frequencies are not altered in irc20∆. These results suggest Irc20 primarily regulates mitotic recombination and does not affect meiotic crossovers. Our results suggest that the IRC20 gene is important for regulating LOH frequency and distribution.

A Proteogenomic Pipeline for the Analysis of Protein Biosynthesis Errors in the Human Pathogen Candida albicans

Candida albicans is a diploid pathogen known for its ability to live as a commensal fungus in healthy individuals but causing both superficial infections and disseminated candidiasis in immunocompromised patients where it is associated with high morbidity and mortality. Its success in colonizing the human host is attributed to a wide range of virulence traits that modulate interactions between the host and the pathogen, such as optimal growth rate at 37 °C, the ability to switch between yeast and hyphal forms, and a remarkable genomic and phenotypic plasticity. A fascinating aspect of its biology is a prominent heterogeneous proteome that arises from frequent genomic rearrangements, high allelic variation, and high levels of amino acid misincorporations in proteins. This leads to increased morphological and physiological phenotypic diversity of high adaptive potential, but the scope of such protein mistranslation is poorly understood due to technical difficulties in detecting and quantifying amino acid misincorporation events in complex protein samples. We have developed and optimized mass spectrometry and bioinformatics pipelines capable of identifying rare amino acid misincorporation events at the proteome level. We have also analyzed the proteomic profile of an engineered C. albicans strain that exhibits high level of leucine misincorporation at protein CUG sites and employed an in vivo quantitative gain-of-function fluorescence reporter system to validate our LC-MS/MS data. C. albicans misincorporates amino acids above the background level at protein sites of diverse codons, particularly at CUG, confirming our previous data on the quantification of leucine incorporation at single CUG sites of recombinant reporter proteins, but increasing misincorporation of Leucine at these sites does not alter the translational fidelity of the other codons. These findings indicate that the C. albicans statistical proteome exceeds prior estimates, suggesting that its highly plastic phenome may also be modulated by environmental factors due to translational ambiguity.

The limit to evolutionary rescue depends on ploidy in yeast exposed to nystatin

The number of copies of each chromosome, or ploidy, of an organism is a major genomic factor affecting adaptation. We set out to determine how ploidy can impact the outcome of evolution, as well as the likelihood of evolutionary rescue, using short-term experiments with yeast (Saccharomyces cerevisiae) in a high concentration of the fungicide nystatin. In similar experiments using haploid yeast, the genetic changes underlying evolutionary rescue were highly repeatable, with all rescued lines containing a single mutation in the ergosterol biosynthetic pathway. All of these beneficial mutations were recessive, which led to the expectation that diploids would find alternative genetic routes to adaptation. To test this, we repeated the experiment using both haploid and diploid strains and found that diploid populations did not evolve resistance. Although diploids are able to adapt at the same rate as haploids to a lower, not fully inhibitory, concentration of nystatin, the present study suggests that diploids are limited in their ability to adapt to an inhibitory concentration of nystatin, while haploids may undergo evolutionary rescue. These results demonstrate that ploidy can tip the balance between adaptation and extinction when organisms face an extreme environmental change.

The cytidine deaminase APOBEC3C has unique sequence and genome feature preferences

APOBEC proteins are cytidine deaminases that restrict the replication of viruses and transposable elements. Several members of the APOBEC3 family, APOBEC3A, APOBEC3B, and APOBEC3H-I, can access the nucleus and cause what is thought to be indiscriminate deamination of the genome, resulting in mutagenesis and genome instability. Although APOBEC3C is also present in the nucleus, the full scope of its deamination target preferences is unknown. By expressing human APOBEC3C in a yeast model system, I have defined the APOBEC3C mutation signature, as well as the preferred genome features of APOBEC3C targets. The APOBEC3C mutation signature is distinct from those of the known cancer genome mutators APOBEC3A and APOBEC3B. APOBEC3C produces DNA strand-coordinated mutation clusters, and APOBEC3C mutations are enriched near the transcription start sites of active genes. Surprisingly, APOBEC3C lacks the bias for the lagging strand of DNA replication that is seen for APOBEC3A and APOBEC3B. The unique preferences of APOBEC3C constitute a mutation profile that will be useful in defining sites of APOBEC3C mutagenesis in human genomes.

Variation in the Spectrum of New Mutations among Inbred Strains of Mice

The mouse serves as a mammalian model for understanding the nature of variation from new mutations, a question that has both evolutionary and medical significance. Previous studies suggest that the rate of single-nucleotide mutations (SNMs) in mice is ∼50% of that in humans. However, information largely comes from studies involving the C57BL/6 strain, and there is little information from other mouse strains. Here, we study the mutations that accumulated in 59 mouse lines derived from four inbred strains that are commonly used in genetics and clinical research (BALB/cAnNRj, C57BL/6JRj, C3H/HeNRj, and FVB/NRj), maintained for eight to nine generations by brother-sister mating. By analyzing Illumina whole-genome sequencing data, we estimate that the average rate of new SNMs in mice is ∼μ = 6.7 × 10-9. However, there is substantial variation in the spectrum of SNMs among strains, so the burden from new mutations also varies among strains. For example, the FVB strain has a spectrum that is markedly skewed toward C→A transversions and is likely to experience a higher deleterious load than other strains, due to an increased frequency of nonsense mutations in glutamic acid codons. Finally, we observe substantial variation in the rate of new SNMs among DNA sequence contexts, CpG sites, and their adjacent nucleotides playing an important role.

How antisense transcripts can evolve to encode novel proteins

Protein coding features can emerge de novo in non coding transcripts, resulting in emergence of new protein coding genes. Studies across many species show that a large fraction of evolutionarily novel non-coding RNAs have an antisense overlap with protein coding genes. The open reading frames (ORFs) in these antisense RNAs could also overlap with existing ORFs. In this study, we investigate how the evolution an ORF could be constrained by its overlap with an existing ORF in three different reading frames. Using a combination of mathematical modeling and genome/transcriptome data analysis in two different model organisms, we show that antisense overlap can increase the likelihood of ORF emergence and reduce the likelihood of ORF loss, especially in one of the three reading frames. In addition to rationalising the repeatedly reported prevalence of de novo emerged genes in antisense transcripts, our work also provides a generic modeling and an analytical framework that can be used to understand evolution of antisense genes.

Modeling Length Changes in De Novo Open Reading Frames during Neutral Evolution

For protein coding genes to emerge de novo from a non-genic DNA, the DNA sequence must gain an open reading frame (ORF) and the ability to be transcribed. The newborn de novo gene can further evolve to accumulate changes in its sequence. Consequently, it can also elongate or shrink with time. Existing literature shows that older de novo genes have longer ORF, but it is not clear if they elongated with time or remained of the same length since their inception. To address this question we developed a mathematical model of ORF elongation as a Markov-jump process, and show that ORFs tend to keep their length in short evolutionary timescales. We also show that if change occurs it is likely to be a truncation. Our genomics and transcriptomics data analyses of seven Drosophila melanogaster populations are also in agreement with the model's prediction. We conclude that selection could facilitate ORF length extension that may explain why longer ORFs were observed in old de novo genes in studies analysing longer evolutionary time scales. Alternatively, shorter ORFs may be purged because they may be less likely to yield functional proteins.

Spontaneous single-nucleotide substitutions and microsatellite mutations have distinct distributions of fitness effects

The fitness effects of new mutations determine key properties of evolutionary processes. Beneficial mutations drive evolution, yet selection is also shaped by the frequency of small-effect deleterious mutations, whose combined effect can burden otherwise adaptive lineages and alter evolutionary trajectories and outcomes in clonally evolving organisms such as viruses, microbes, and tumors. The small effect sizes of these important mutations have made accurate measurements of their rates difficult. In microbes, assessing the effect of mutations on growth can be especially instructive, as this complex phenotype is closely linked to fitness in clonally evolving organisms. Here, we perform high-throughput time-lapse microscopy on cells from mutation-accumulation strains to precisely infer the distribution of mutational effects on growth rate in the budding yeast, Saccharomyces cerevisiae. We show that mutational effects on growth rate are overwhelmingly negative, highly skewed towards very small effect sizes, and frequent enough to suggest that deleterious hitchhikers may impose a significant burden on evolving lineages. By using lines that accumulated mutations in either wild-type or slippage repair-defective backgrounds, we further disentangle the effects of two common types of mutations, single-nucleotide substitutions and simple sequence repeat indels, and show that they have distinct effects on yeast growth rate. Although the average effect of a simple sequence repeat mutation is very small (~0.3%), many do alter growth rate, implying that this class of frequent mutations has an important evolutionary impact.

Genomic factors shaping codon usage across the Saccharomycotina subphylum

Codon usage bias, or the unequal use of synonymous codons, is observed across genes, genomes, and between species. The biased use of synonymous codons has been implicated in many cellular functions, such as translation dynamics and transcript stability, but can also be shaped by neutral forces. The Saccharomycotina, the fungal subphylum containing the yeasts Saccharomyces cerevisiae and Candida albicans , has been a model system for studying codon usage. We characterized codon usage across 1,154 strains from 1,051 species to gain insight into the biases, molecular mechanisms, evolution, and genomic features contributing to codon usage patterns across the subphylum. We found evidence of a general preference for A/T-ending codons and correlations between codon usage bias, GC content, and tRNA-ome size. Codon usage bias is also distinct between the 12 orders within the subphylum to such a degree that yeasts can be classified into orders with an accuracy greater than 90% using a machine learning algorithm trained on codon usage. We also characterized the degree to which codon usage bias is impacted by translational selection. Interestingly, the degree of translational selection was influenced by a combination of genome features and assembly metrics that included the number of coding sequences, BUSCO count, and genome length. Our analysis also revealed an extreme bias in codon usage in the Saccharomycodales associated with a lack of predicted arginine tRNAs. The order contains 24 species, and 23 are computationally predicted to lack tRNAs that decode CGN codons, leaving only the AGN codons to encode arginine. Analysis of Saccharomycodales gene expression, tRNA sequences, and codon evolution suggests that extreme avoidance of the CGN codons is associated with a decline in arginine tRNA function. Codon usage bias within the Saccharomycotina is generally consistent with previous investigations in fungi, which show a role for both genomic features and GC bias in shaping codon usage. However, we find cases of extreme codon usage preference and avoidance along yeast lineages, suggesting additional forces may be shaping the evolution of specific codons.

Mobile genetic element-based gene editing and genome engineering: Recent advances and applications

Genome engineering has revolutionized several scientific fields, ranging from biochemistry and fundamental research to therapeutic uses and crop development. Diverse engineering toolkits have been developed and used to effectively modify the genome sequences of organisms. However, there is a lack of extensive reviews on genome engineering technologies based on mobile genetic elements (MGEs), which induce genetic diversity within host cells by changing their locations in the genome. This review provides a comprehensive update on the versatility of MGEs as powerful genome engineering tools that offers efficient solutions to challenges associated with genome engineering. MGEs, including DNA transposons, retrotransposons, retrons, and CRISPR-associated transposons, offer various advantages, such as a broad host range, genome-wide mutagenesis, efficient large-size DNA integration, multiplexing capabilities, and in situ single-stranded DNA generation. We focused on the components, mechanisms, and features of each MGE-based tool to highlight their cellular applications. Finally, we discussed the current challenges of MGE-based genome engineering and provided insights into the evolving landscape of this transformative technology. In conclusion, the combination of genome engineering with MGE demonstrates remarkable potential for addressing various challenges and advancing the field of genetic manipulation, and promises to revolutionize our ability to engineer and understand the genomes of diverse organisms.

Aneuploidy Can Be an Evolutionary Diversion on the Path to Adaptation

Aneuploidy is common in eukaryotes, often leading to decreased fitness. However, evidence from fungi and human tumur cells suggests that specific aneuploidies can be beneficial under stressful conditions and facilitate adaptation. In a previous evolutionary experiment with yeast, populations evolving under heat stress became aneuploid, only to later revert to euploidy after beneficial mutations accumulated. It was therefore suggested that aneuploidy is a "stepping stone" on the path to adaptation. Here, we test this hypothesis. We use Bayesian inference to fit an evolutionary model with both aneuploidy and mutation to the experimental results. We then predict the genotype frequency dynamics during the experiment, demonstrating that most of the evolved euploid population likely did not descend from aneuploid cells, but rather from the euploid wild-type population. Our model shows how the beneficial mutation supply-the product of population size and beneficial mutation rate-determines the evolutionary dynamics: with low supply, much of the evolved population descends from aneuploid cells; but with high supply, beneficial mutations are generated fast enough to outcompete aneuploidy due to its inherent fitness cost. Our results suggest that despite its potential fitness benefits under stress, aneuploidy can be an evolutionary "diversion" rather than a "stepping stone": it can delay, rather than facilitate, the adaptation of the population, and cells that become aneuploid may leave less descendants compared to cells that remain diploid.

Candida albicans' inorganic phosphate transport and evolutionary adaptation to phosphate scarcity

Phosphorus is essential in all cells' structural, metabolic and regulatory functions. For fungal cells that import inorganic phosphate (Pi) up a steep concentration gradient, surface Pi transporters are critical capacitators of growth. Fungi must deploy Pi transporters that enable optimal Pi uptake in pH and Pi concentration ranges prevalent in their environments. Single, triple and quadruple mutants were used to characterize the four Pi transporters we identified for the human fungal pathogen Candida albicans, which must adapt to alkaline conditions during invasion of the host bloodstream and deep organs. A high-affinity Pi transporter, Pho84, was most efficient across the widest pH range while another, Pho89, showed high-affinity characteristics only within one pH unit of neutral. Two low-affinity Pi transporters, Pho87 and Fgr2, were active only in acidic conditions. Only Pho84 among the Pi transporters was clearly required in previously identified Pi-related functions including Target of Rapamycin Complex 1 signaling and hyphal growth. We used in vitro evolution and whole genome sequencing as an unbiased forward genetic approach to probe adaptation to prolonged Pi scarcity of two quadruple mutant lineages lacking all 4 Pi transporters. Lineage-specific genomic changes corresponded to divergent success of the two lineages in fitness recovery during Pi limitation. In this process, initial, large-scale genomic alterations like aneuploidies and loss of heterozygosity were eventually lost as populations presumably gained small-scale mutations. Severity of some phenotypes linked to Pi starvation, like cell wall stress hypersensitivity, decreased in parallel to evolving populations' fitness recovery in Pi scarcity, while that of others like membrane stress responses diverged from these fitness phenotypes. C. albicans therefore has diverse options to reconfigure Pi management during prolonged scarcity. Since Pi homeostasis differs substantially between fungi and humans, adaptive processes to Pi deprivation may harbor small-molecule targets that impact fungal growth and virulence.

Upper Bound on the Mutational Burden Imposed by a CRISPR-Cas9 Gene-Drive Element

CRISPR-Cas9 gene drives (CCGDs) are powerful tools for genetic control of wild populations, useful for eradication of disease vectors, conservation of endangered species and other applications. However, Cas9 alone and in a complex with gRNA can cause double-stranded DNA breaks at off-target sites, which could increase the mutational load and lead to loss of heterozygosity (LOH). These undesired effects raise potential concerns about the long-term evolutionary safety of CCGDs, but the magnitude of these effects is unknown. To estimate how the presence of a CCGD or a Cas9 alone in the genome affects the rates of LOH events and de novo mutations, we carried out a mutation accumulation experiment in yeast Saccharomyces cerevisiae. Despite its substantial statistical power, our experiment revealed no detectable effect of CCGD or Cas9 alone on the genome-wide rates of mutations or LOH events, suggesting that these rates are affected by less than 30%. Nevertheless, we found that Cas9 caused a slight but significant shift towards more interstitial and fewer terminal LOH events, and the CCGD caused a significant difference in the distribution of LOH events on Chromosome V. Taken together, our results show that these genetic elements impose a weak and likely localized additional mutational burden in the yeast model. Although the mutagenic effects of CCGDs need to be further evaluated in other systems, our results suggest that the effect of CCGDs on off-target mutation rates and genetic diversity may be acceptable.

Gene birth in a model of non-genic adaptation

BACKGROUND

Over evolutionary timescales, genomic loci can switch between functional and non-functional states through processes such as pseudogenization and de novo gene birth. Particularly, de novo gene birth is a widespread process, and many examples continue to be discovered across diverse evolutionary lineages. However, the general mechanisms that lead to functionalization are poorly understood, and estimated rates of de novo gene birth remain contentious. Here, we address this problem within a model that takes into account mutations and structural variation, allowing us to estimate the likelihood of emergence of new functions at non-functional loci.

RESULTS

Assuming biologically reasonable mutation rates and mutational effects, we find that functionalization of non-genic loci requires the realization of strict conditions. This is in line with the observation that most de novo genes are localized to the vicinity of established genes. Our model also provides an explanation for the empirical observation that emerging proto-genes are often lost despite showing signs of adaptation.

CONCLUSIONS

Our work elucidates the properties of non-genic loci that make them fertile for adaptation, and our results offer mechanistic insights into the process of de novo gene birth.

Evolution of the Mutation Spectrum Across a Mammalian Phylogeny

Although evolutionary biologists have long theorized that variation in DNA repair efficacy might explain some of the diversity of lifespan and cancer incidence across species, we have little data on the variability of normal germline mutagenesis outside of humans. Here, we shed light on the spectrum and etiology of mutagenesis across mammals by quantifying mutational sequence context biases using polymorphism data from thirteen species of mice, apes, bears, wolves, and cetaceans. After normalizing the mutation spectrum for reference genome accessibility and k-mer content, we use the Mantel test to deduce that mutation spectrum divergence is highly correlated with genetic divergence between species, whereas life history traits like reproductive age are weaker predictors of mutation spectrum divergence. Potential bioinformatic confounders are only weakly related to a small set of mutation spectrum features. We find that clock-like mutational signatures previously inferred from human cancers cannot explain the phylogenetic signal exhibited by the mammalian mutation spectrum, despite the ability of these signatures to fit each species' 3-mer spectrum with high cosine similarity. In contrast, parental aging signatures inferred from human de novo mutation data appear to explain much of the 1-mer spectrum's phylogenetic signal in combination with a novel mutational signature. We posit that future models purporting to explain the etiology of mammalian mutagenesis need to capture the fact that more closely related species have more similar mutation spectra; a model that fits each marginal spectrum with high cosine similarity is not guaranteed to capture this hierarchy of mutation spectrum variation among species.

Genome-wide probing of eukaryotic nascent RNA structure elucidates cotranscriptional folding and its antimutagenic effect

The transcriptional intermediates of RNAs fold into secondary structures with multiple regulatory roles, yet the details of such cotranscriptional RNA folding are largely unresolved in eukaryotes. Here, we present eSPET-seq (Structural Probing of Elongating Transcripts in eukaryotes), a method to assess the cotranscriptional RNA folding in Saccharomyces cerevisiae. Our study reveals pervasive structural transitions during cotranscriptional folding and overall structural similarities between nascent and mature RNAs. Furthermore, a combined analysis with genome-wide R-loop and mutation rate approximations provides quantitative evidence for the antimutator effect of nascent RNA folding through competitive inhibition of the R-loops, known to facilitate transcription-associated mutagenesis. Taken together, we present an experimental evaluation of cotranscriptional folding in eukaryotes and demonstrate the antimutator effect of nascent RNA folding. These results suggest genome-wide coupling between the processing and transmission of genetic information through RNA folding.

The Dynamic Fungal Genome: Polyploidy, Aneuploidy and Copy Number Variation in Response to Stress

Fungal species have dynamic genomes and often exhibit genomic plasticity in response to stress. This genome plasticity often comes with phenotypic consequences that affect fitness and resistance to stress. Fungal pathogens exhibit genome plasticity in both clinical and agricultural settings and often during adaptation to antifungal drugs, posing significant challenges to human health. Therefore, it is important to understand the rates, mechanisms, and impact of large genomic changes. This review addresses the prevalence of polyploidy, aneuploidy, and copy number variation across diverse fungal species, with special attention to prominent fungal pathogens and model species. We also explore the relationship between environmental stress and rates of genomic changes and highlight the mechanisms underlying genotypic and phenotypic changes. A comprehensive understanding of these dynamic fungal genomes is needed to identify novel solutions for the increase in antifungal drug resistance.

Environmental variance in male mating success modulates the positive versus negative impacts of sexual selection on genetic load

Sexual selection on males is predicted to increase population fitness, and delay population extinction, when mating success negatively covaries with genetic load across individuals. However, such benefits of sexual selection could be counteracted by simultaneous increases in genome-wide drift resulting from reduced effective population size caused by increased variance in fitness. Resulting fixation of deleterious mutations could be greatest in small populations, and when environmental variation in mating traits partially decouples sexual selection from underlying genetic variation. The net consequences of sexual selection for genetic load and population persistence are therefore likely to be context dependent, but such variation has not been examined. We use a genetically explicit individual-based model to show that weak sexual selection can increase population persistence time compared to random mating. However, for stronger sexual selection such positive effects can be overturned by the detrimental effects of increased genome-wide drift. Furthermore, the relative strengths of mutation-purging and drift critically depend on the environmental variance in the male mating trait. Specifically, increasing environmental variance caused stronger sexual selection to elevate deleterious mutation fixation rate and mean selection coefficient, driving rapid accumulation of drift load and decreasing population persistence times. These results highlight an intricate balance between conflicting positive and negative consequences of sexual selection on genetic load, even in the absence of sexually antagonistic selection. They imply that environmental variances in key mating traits, and intrinsic genetic drift, should be properly factored into future theoretical and empirical studies of the evolution of population fitness under sexual selection.

Genotyping and susceptibility testing uncovers large azole-resistant Candida tropicalis clade in Alexandria, Egypt

OBJECTIVES

Candida tropicalis is an emerging medically relevant Candida species. The yeast primarily causes opportunistic infections in intensive care units and is highly prevalent in tropical countries. The genetic diversity within this species is high, and nosocomial transmission has been reported. C. tropicalis genotyping of isolates from low- and middle-income countries is underrepresented when compared with that from high-income countries. Also, in Egypt, only limited genotyping has been conducted for C. tropicalis isolates, while antifungal resistance seems to increase, especially against azoles.

METHODS

Antifungal susceptibility testing was performed on 64 C. tropicalis isolates from ICU patients collected from multiple hospitals in Alexandria, Egypt. Genotyping by means of short tandem repeat (STR) and whole genome sequencing (WGS) single nucleotide polymorphism (SNP) analysis was performed.

RESULTS

Using antifungal susceptibility testing, fluconazole resistance was observed in 24 isolates (38%), of which 23 harboured an ERG11 G464S substitution, previously shown to cause resistance in Candida albicans. STR genotyping showed that these 23 isolates were related, forming a distinct resistant clade. WGS SNP analysis subsequently confirmed this genetic relationship, although isolates within this clade differed in at least 429 SNPs, suggesting that these were independently introduced.

CONCLUSION

Overall, STR and WGS SNP analysis of this collection indicates limited C. tropicalis nosocomial transmission in Alexandria, while the presence of this large azole-resistant C. tropicalis clade within this city hampers the treatment of intensive care unit patients.

Phenotype switching of the mutation rate facilitates adaptive evolution

The mutation rate plays an important role in adaptive evolution. It can be modified by mutator and anti-mutator alleles. Recent empirical evidence hints that the mutation rate may vary among genetically identical individuals: evidence from bacteria suggests that the mutation rate can be affected by expression noise of a DNA repair protein and potentially also by translation errors in various proteins. Importantly, this non-genetic variation may be heritable via a transgenerational epigenetic mode of inheritance, giving rise to a mutator phenotype that is independent from mutator alleles. Here, we investigate mathematically how the rate of adaptive evolution is affected by the rate of mutation rate phenotype switching. We model an asexual population with two mutation rate phenotypes, non-mutator and mutator. An offspring may switch from its parental phenotype to the other phenotype. We find that switching rates that correspond to so-far empirically described non-genetic systems of inheritance of the mutation rate lead to higher rates of adaptation on both artificial and natural fitness landscapes. These switching rates can maintain within the same individuals both a mutator phenotype and intermediary mutations, a combination that facilitates adaptation. Moreover, non-genetic inheritance increases the proportion of mutators in the population, which in turn increases the probability of hitchhiking of the mutator phenotype with adaptive mutations. This in turns facilitates the acquisition of additional adaptive mutations. Our results rationalize recently observed noise in the expression of proteins that affect the mutation rate and suggest that non-genetic inheritance of this phenotype may facilitate evolutionary adaptive processes.

Dynamics of Podospora anserina Genome Evolution in a Long-Term Experiment

The Podospora anserina long-term evolution experiment (PaLTEE) is the only running filamentous fungus study, which is still going on. The aim of our work is to trace the evolutionary dynamics of the accumulation of mutations in the genomes of eight haploid populations of P. anserina. The results of the genome-wide analysis of all of the lineages, performed 8 years after the start of the PaLTEE, are presented. Data analysis detected 312 single nucleotide polymorphisms (SNPs) and 39 short insertion-deletion mutations (indels) in total. There was a clear trend towards a linear increase in the number of SNPs depending on the experiment duration. Among 312 SNPs, 153 were fixed in the coding regions of P. anserina genome. Relatively few synonymous mutations were found, exactly 38; 42 were classified as nonsense mutations; 72 were assigned to missense mutations. In addition, 21 out of 39 indels identified were also localized in coding regions. Here, we also report the detection of parallel evolution at the paralog level in the P. anserina model system. Parallelism in evolution at the level of protein functions also occurs. The latter is especially true for various transcription factors, which may indicate selection leading to optimization of the wide range of cellular processes under experimental conditions.

The use of whole-genome sequencing and development of bioinformatics to monitor overlapping outbreaks of Candida auris in southern Nevada

A Candida auris outbreak has been ongoing in Southern Nevada since August 2021. In this manuscript we describe the sequencing of over 200 C. auris isolates from patients at several facilities. Genetically distinct subgroups of C. auris were detected from Clade I (3 distinct lineages) and III (1 lineage). Open-source bioinformatic tools were developed and implemented to aid in the epidemiological investigation. The work herein compares three methods for C. auris whole genome analysis: Nullarbor, MycoSNP and a new pipeline TheiaEuk. We also describe a novel analysis method focused on elucidating phylogenetic linkages between isolates within an ongoing outbreak. Moreover, this study places the ongoing outbreaks in a global context utilizing existing sequences provided worldwide. Lastly, we describe how the generated results were communicated to the epidemiologists and infection control to generate public health interventions.

The effects of mutations on gene expression and alternative splicing

Understanding the relationship between mutations and their genomic and phenotypic consequences has been a longstanding goal of evolutionary biology. However, few studies have investigated the impact of mutations on gene expression and alternative splicing on the genome-wide scale. In this study, we aim to bridge this knowledge gap by utilizing whole-genome sequencing data and RNA sequencing data from 16 obligately parthenogenetic Daphnia mutant lines to investigate the effects of ethyl methanesulfonate-induced mutations on gene expression and alternative splicing. Using rigorous analyses of mutations, expression changes and alternative splicing, we show that trans-effects are the major contributor to the variance in gene expression and alternative splicing between the wild-type and mutant lines, whereas cis mutations only affected a limited number of genes and do not always alter gene expression. Moreover, we show that there is a significant association between differentially expressed genes and exonic mutations, indicating that exonic mutations are an important driver of altered gene expression.

Experimental evolution of Saccharomyces uvarum at high temperature yields elevation of maximal growth temperature and loss of the mitochondrial genome

An organism's upper thermal tolerance is a major driver of its ecology and is a complex, polygenic trait. Given the wide variance in this critical phenotype across the tree of life, it is quite striking that this trait has not proven very evolutionarily labile in experimental evolution studies of microbes. In stark contrast to recent studies, William Henry Dallinger in the 1880s reported increasing the upper thermal limit of microbes he experimentally evolved by >40°C using a very gradual temperature ramping strategy. Using a selection scheme inspired by Dallinger, we sought to increase the upper thermal limit of Saccharomyces uvarum . This species has a maximum growth temperature of 34-35°C, considerably lower than S. cerevisiae . After 136 passages on solid plates at progressively higher temperatures, we recovered a clone that can grow at 36°C, a gain of ~1.5°C. Additionally, the evolved clone lost its mitochondrial genome and cannot respire. In contrast, an induced rho 0 derivative of the ancestor shows a decrease in thermotolerance. Also, incubation of the ancestor at 34°C for 5 days increased the frequency of petite mutants drastically compared to 22°C, supporting the notion that mutation pressure rather than selection drove loss of mtDNA in the evolved clone. These results demonstrate that S. uvarum 's upper thermal limit can be elevated slightly via experimental evolution and corroborate past observations in S. cerevisiae that high temperature selection schemes can inadvertently lead to production of the potentially undesirable respiratory incompetent phenotype in yeasts.

Loss of Heterozygosity and Its Importance in Evolution

Loss of heterozygosity (LOH) is a mitotic recombination event that converts heterozygous loci to homozygous loci. This mutation event is widespread in organisms that have asexual reproduction like budding yeasts, and is also an important and frequent mutation event in tumorigenesis. Mutation accumulation studies have demonstrated that LOH occurs at a rate higher than the point mutation rate, and can impact large portions of the genome. Laboratory evolution experiments of heterozygous yeasts have revealed that LOH often unmasks beneficial recessive alleles that can confer large fitness advantages. Here, I highlight advances in understanding dominance, fitness, and phenotypes in laboratory evolved heterozygous yeast strains. I discuss best practices for detecting LOH in intraspecific and interspecific evolved clones and populations. Utilizing heterozygous strain backgrounds in laboratory evolution experiments offers an opportunity to advance our understanding of this important mutation type in shaping adaptation and genome evolution in wild, domesticated, and clinical populations.

"Evolution of the mutation spectrum across a mammalian phylogeny"

Little is known about how the spectrum and etiology of germline mutagenesis might vary among mammalian species. To shed light on this mystery, we quantify variation in mutational sequence context biases using polymorphism data from thirteen species of mice, apes, bears, wolves, and cetaceans. After normalizing the mutation spectrum for reference genome accessibility and k -mer content, we use the Mantel test to deduce that mutation spectrum divergence is highly correlated with genetic divergence between species, whereas life history traits like reproductive age are weaker predictors of mutation spectrum divergence. Potential bioinformatic confounders are only weakly related to a small set of mutation spectrum features. We find that clocklike mutational signatures previously inferred from human cancers cannot explain the phylogenetic signal exhibited by the mammalian mutation spectrum, despite the ability of these clocklike signatures to fit each species' 3-mer spectrum with high cosine similarity. In contrast, parental aging signatures inferred from human de novo mutation data appear to explain much of the mutation spectrum's phylogenetic signal when fit to non-context-dependent mutation spectrum data in combination with a novel mutational signature. We posit that future models purporting to explain the etiology of mammalian mutagenesis need to capture the fact that more closely related species have more similar mutation spectra; a model that fits each marginal spectrum with high cosine similarity is not guaranteed to capture this hierarchy of mutation spectrum variation among species.

Low Spontaneous Mutation Rate in Complex Multicellular Eukaryotes with a Haploid-Diploid Life Cycle

The spontaneous mutation rate µ is a crucial parameter to understand evolution and biodiversity. Mutation rates are highly variable across species, suggesting that µ is susceptible to selection and drift and that species life cycle and life history may impact its evolution. In particular, asexual reproduction and haploid selection are expected to affect the mutation rate, but very little empirical data are available to test this expectation. Here, we sequence 30 genomes of a parent-offspring pedigree in the model brown alga Ectocarpus sp.7, and 137 genomes of an interspecific cross of the closely related brown alga Scytosiphon to have access to the spontaneous mutation rate of representative organisms of a complex multicellular eukaryotic lineage outside animals and plants, and to evaluate the potential impact of life cycle on the mutation rate. Brown algae alternate between a haploid and a diploid stage, both multicellular and free living, and utilize both sexual and asexual reproduction. They are, therefore, excellent models to empirically test expectations of the effect of asexual reproduction and haploid selection on mutation rate evolution. We estimate that Ectocarpus has a base substitution rate of µbs = 4.07 × 10-10 per site per generation, whereas the Scytosiphon interspecific cross had µbs = 1.22 × 10-9. Overall, our estimations suggest that these brown algae, despite being multicellular complex eukaryotes, have unusually low mutation rates. In Ectocarpus, effective population size (Ne) could not entirely explain the low µbs. We propose that the haploid-diploid life cycle, combined with extensive asexual reproduction, may be additional key drivers of the mutation rate in these organisms.

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.

Mutation bias and the predictability of evolution

Predicting evolutionary outcomes is an important research goal in a diversity of contexts. The focus of evolutionary forecasting is usually on adaptive processes, and efforts to improve prediction typically focus on selection. However, adaptive processes often rely on new mutations, which can be strongly influenced by predictable biases in mutation. Here, we provide an overview of existing theory and evidence for such mutation-biased adaptation and consider the implications of these results for the problem of prediction, in regard to topics such as the evolution of infectious diseases, resistance to biochemical agents, as well as cancer and other kinds of somatic evolution. We argue that empirical knowledge of mutational biases is likely to improve in the near future, and that this knowledge is readily applicable to the challenges of short-term prediction. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

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.

The rate of spontaneous mutations in yeast deficient for MutSβ function

Mutations in simple sequence repeat loci underlie many inherited disorders in humans, and are increasingly recognized as important determinants of natural phenotypic variation. In eukaryotes, mutations in these sequences are primarily repaired by the MutSβ mismatch repair complex. To better understand the role of this complex in mismatch repair and the determinants of simple sequence repeat mutation predisposition, we performed mutation accumulation in yeast strains with abrogated MutSβ function. We demonstrate that mutations in simple sequence repeat loci in the absence of mismatch repair are primarily deletions. We also show that mutations accumulate at drastically different rates in short (<8 bp) and longer repeat loci. These data lend support to a model in which the mismatch repair complex is responsible for repair primarily in longer simple sequence repeats.

Spontaneous Mutation Rates and Spectra of Respiratory-Deficient Yeast

The yeast petite mutant was first discovered in the yeast Saccharomyces cerevisiae, which shows growth stress due to defects in genes encoding the respiratory chain. In a previous study, we described that deletion of the nuclear-encoded gene MRPL25 leads to mitochondrial genome (mtDNA) loss and the petite phenotype, which can be rescued by acquiring ATP3 mutations. The mrpl25Δ strain showed an elevated SNV (single nucleotide variant) rate, suggesting genome instability occurred during the crisis of mtDNA loss. However, the genome-wide mutation landscape and mutational signatures of mitochondrial dysfunction are unknown. In this study we profiled the mutation spectra in yeast strains with the genotype combination of MRPL25 and ATP3 in their wildtype and mutated status, along with the wildtype and cytoplasmic petite rho0 strains as controls. In addition to the previously described elevated SNV rate, we found the INDEL (insertion/deletion) rate also increased in the mrpl25Δ strain, reinforcing the occurrence of genome instability. Notably, although both are petites, the mrpl25Δ and rho0 strains exhibited different INDEL rates and transition/transversion ratios, suggesting differences in the mutational signatures underlying these two types of petites. Interestingly, the petite-related mutagenesis effect disappeared when ATP3 suppressor mutations were acquired, suggesting a cost-effective mechanism for restoring both fitness and genome stability. Taken together, we present an unbiased genome-wide characterization of the mutation rates and spectra of yeast strains with respiratory deficiency, which provides valuable insights into the impact of respiratory deficiency on genome instability.

Population genetic analysis of the microsporidium Ordospora colligata reveals the role of natural selection and phylogeography on its extremely compact and reduced genome

The determinants of variation in a species' genome-wide nucleotide diversity include historical, environmental, and stochastic aspects. This diversity can inform us about the species' past and present evolutionary dynamics. In parasites, the mode of transmission and the interactions with the host might supersede the effects of these aspects in shaping parasite genomic diversity. We used genomic samples from 10 populations of the microsporidian parasite Ordospora colligata to investigate present genomic diversity and how it was shaped by evolutionary processes, specifically, the role of phylogeography, co-phylogeography (with the host), natural selection, and transmission mode. Although very closely related microsporidia cause diseases in humans, O. colligata is specific to the freshwater crustacean Daphnia magna and has one of the smallest known eukaryotic genomes. We found an overlapping phylogeography between O. colligata and its host highlighting the long-term, intimate relationship between them. The observed geographic distribution reflects previous findings that O. colligata exhibits adaptations to colder habitats, which differentiates it from other microsporidian gut parasites of D. magna predominantly found in warmer areas. The co-phylogeography allowed us to calibrate the O. colligata phylogeny and thus estimate its mutation rate. We identified several genetic regions under potential selection. Our whole-genome study provides insights into the evolution of one of the most reduced eukaryotic genomes and shows how different processes shape genomic diversity of an obligate parasite.

Phenotypic and genetic alterations of Burkholderia pseudomallei in patients during relapse and persistent infections

The bacterium Burkholderia pseudomallei is the causative agent of melioidosis, a severe tropical disease associated with high mortality and relapse and persistent infections. Treatment of melioidosis requires prolonged antibiotic therapy; however, little is known about relapse and persistent infections, particularly the phenotypic and genetic alterations of B. pseudomallei in patients. In this study, we performed pulsed-field gel electrophoresis (PFGE) to compare the bacterial genotype between the initial isolate and the subsequent isolate from each of 23 suspected recurrent and persistent melioidosis patients in Northeast Thailand. We used whole-genome sequencing (WGS) to investigate multilocus sequence types and genetic alterations of within-host strain pairs. We also investigated the bacterial phenotypes associated with relapse and persistent infections, including multinucleated giant cell (MNGC) formation efficiency and intracellular multiplication. We first identified 13 (1.2%) relapse, 7 (0.7%) persistent, and 3 (0.3%) reinfection patients from 1,046 survivors. Each of the 20 within-host strain pairs from patients with relapse and persistent infections shared the same genotype, suggesting that the subsequent isolates arise from the infecting isolate. Logistic regression analysis of clinical data revealed regimen and duration of oral antibiotic therapies as risk factors associated with relapse and persistent infections. WGS analysis demonstrated 17 within-host genetic alteration events in 6 of 20 paired isolates, including a relatively large deletion and 16 single-nucleotide polymorphism (stocktickerSNP) mutations distributed across 12 genes. In 1 of 20 paired isolates, we observed significantly increased cell-to-cell fusion and intracellular replication in the second isolate compared with the initial isolate from a patient with persistent infection. WGS analysis suggested that a non-synonymous mutation in the tssB-5 gene, which encoded an essential component of the type VI secretion system, may be associated with the increased intracellular replication and MNGC formation efficiency of the second isolate of the patient. This information provides insights into genetic and phenotypic alterations in B. pseudomallei in human melioidosis, which may represent a bacterial strategy for persistent and relapse infections.

Very Low Rates of Spontaneous Gene Deletions and Gene Duplications in Dictyostelium discoideum

The study of spontaneous mutation rates has revealed a wide range of heritable point mutation rates across species, but there are comparatively few estimates for large-scale deletion and duplication rates. The handful of studies that have directly calculated spontaneous rates of deletion and duplication using mutation accumulation lines have estimated that genes are duplicated and deleted at orders of magnitude greater rates than the spontaneous point mutation rate. In our study, we tested whether spontaneous gene deletion and gene duplication rates are also high in Dictyostelium discoideum, a eukaryote with among the lowest point mutation rates (2.5 × 10-11 per site per generation) and an AT-rich genome (GC content of 22%). We calculated mutation rates of gene deletions and duplications using whole-genome sequencing data originating from a mutation accumulation experiment and determined the association between the copy number mutations and GC content. Overall, we estimated an average of 3.93 × 10-8 gene deletions and 1.18 × 10-8 gene duplications per gene per generation. While orders of magnitude greater than their point mutation rate, these rates are much lower compared to gene deletion and duplication rates estimated from mutation accumulation lines in other organisms (that are on the order of ~ 10-6 per gene/generation). The deletions and duplications were enriched in regions that were AT-rich even compared to the genomic background, in contrast to our expectations if low GC content was contributing to low mutation rates. The low deletion and duplication mutation rates in D. discoideum compared to other eukaryotes mirror their low point mutation rates, supporting previous work suggesting that this organism has high replication fidelity and effective molecular machinery to avoid the accumulation of mutations in their genome.

RatesTools: a Nextflow pipeline for detecting de novo germline mutations in pedigree sequence data

SUMMARY

Here, we introduce RatesTools, an automated pipeline to infer de novo mutation rates from parent-offspring trio data of diploid organisms. By providing a reference genome and high-coverage, whole-genome resequencing data of a minimum of three individuals (sire, dam and offspring), RatesTools provides a list of candidate de novo mutations and calculates a putative mutation rate. RatesTools uses several quality filtering steps, such as discarding sites with low mappability and highly repetitive regions, as well as sites with low genotype and mapping qualities to find potential de novo mutations. In addition, RatesTools implements several optional filters based on post hoc assumptions of the heterozygosity and mutation rate of the organism. Filters are highly customizable to user specifications in order to maximize utility across a wide range of applications.

AVAILABILITY AND IMPLEMENTATION

RatesTools is freely available at https://github.com/campanam/RatesTools under a Creative Commons Zero (CC0) license. The pipeline is implemented in Nextflow (Di Tommaso et al., 2017), Ruby (http://www.ruby-lang.org), Bash (https://www.gnu.org/software/bash/) and R (R Core Team, 2020) with reliance upon several other freely available tools. RatesTools is compatible with macOS and Linux operating systems.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Mutualism-enhancing mutations dominate early adaptation in a two-species microbial community

Species interactions drive evolution while evolution shapes these interactions. The resulting eco-evolutionary dynamics and their repeatability depend on how adaptive mutations available to community members affect fitness and ecologically relevant traits. However, the diversity of adaptive mutations is not well characterized, and we do not know how this diversity is affected by the ecological milieu. Here we use barcode lineage tracking to address this question in a community of yeast Saccharomyces cerevisiae and alga Chlamydomonas reinhardtii that have a net commensal relationship that results from a balance between competitive and mutualistic interactions. We find that yeast has access to many adaptive mutations with diverse ecological consequences, in particular those that increase and reduce the yields of both species. The presence of the alga does not change which mutations are adaptive in yeast (that is, there is no fitness trade-off for yeast between growing alone or with alga), but rather shifts selection to favour yeast mutants that increase the yields of both species and make the mutualism stronger. Thus, in the presence of the alga, adaptative mutations contending for fixation in yeast are more likely to enhance the mutualism, even though cooperativity is not directly favoured by natural selection in our system. Our results demonstrate that ecological interactions not only alter the trajectory of evolution but also dictate its repeatability; in particular, weak mutualisms can repeatably evolve to become stronger.

Can changes in 3D genome architecture create new regulatory landscapes that contribute to phenotypic evolution?

Animal genomes are compartmentalized into insulated regulatory units named topology-associated domains (TADs). TADs insulate gene promoters from enhancers that occupy neighboring TADs. Chromosomal rearrangements that disrupt TAD structure can generate new regulatory interactions between enhancers and promoters that were once separated into different TADs, which might lead to new gene expression patterns. On the one hand, TAD rearrangements are known to cause deleterious phenotypes, but, on the other hand, rearrangements can also create novel expression patterns that may be selected during evolution because they generate advantageous phenotypes. Here, we review recent studies that explore the effects of chromosomal rearrangements and genetic perturbations on TAD structure and gene regulation in the context of development and evolution. We discuss the possible contribution of evolutionary breakpoints (EBRs) that affect TAD structure to the evolution of gene regulation and the phenotype.

Architectural groups of a subtelomeric gene family evolve along distinct paths in Candida albicans

Subtelomeres are dynamic genomic regions shaped by elevated rates of recombination, mutation, and gene birth/death. These processes contribute to formation of lineage-specific gene family expansions that commonly occupy subtelomeres across eukaryotes. Investigating the evolution of subtelomeric gene families is complicated by the presence of repetitive DNA and high sequence similarity among gene family members that prevents accurate assembly from whole genome sequences. Here, we investigated the evolution of the telomere-associated (TLO) gene family in Candida albicans using 189 complete coding sequences retrieved from 23 genetically diverse strains across the species. Tlo genes conformed to the 3 major architectural groups (α/β/γ) previously defined in the genome reference strain but significantly differed in the degree of within-group diversity. One group, Tloβ, was always found at the same chromosome arm with strong sequence similarity among all strains. In contrast, diverse Tloα sequences have proliferated among chromosome arms. Tloγ genes formed 7 primary clades that included each of the previously identified Tloγ genes from the genome reference strain with 3 Tloγ genes always found on the same chromosome arm among strains. Architectural groups displayed regions of high conservation that resolved newly identified functional motifs, providing insight into potential regulatory mechanisms that distinguish groups. Thus, by resolving intraspecies subtelomeric gene variation, it is possible to identify previously unknown gene family complexity that may underpin adaptive functional variation.

Eco-evolutionary extinction and recolonization dynamics reduce genetic load and increase time to extinction in highly inbred populations

Understanding how genetic and ecological effects can interact to shape genetic loads within and across local populations is key to understanding ongoing persistence of systems that should otherwise be susceptible to extinction through mutational meltdown. Classic theory predicts short persistence times for metapopulations comprising small local populations with low connectivity, due to accumulation of deleterious mutations. Yet, some such systems have persisted over evolutionary time, implying the existence of mechanisms that allow metapopulations to avoid mutational meltdown. We first hypothesize a mechanism by which the combination of stochasticity in the numbers and types of mutations arising locally (genetic stochasticity), resulting local extinction, and recolonization through evolving dispersal facilitates metapopulation persistence. We then test this mechanism using a spatially and genetically explicit individual-based model. We show that genetic stochasticity in highly structured metapopulations can result in local extinctions, which can favor increased dispersal, thus allowing recolonization of empty habitat patches. This causes fluctuations in metapopulation size and transient gene flow, which reduces genetic load and increases metapopulation persistence over evolutionary time. Our suggested mechanism and simulation results provide an explanation for the conundrum presented by the continued persistence of highly structured populations with inbreeding mating systems that occur in diverse taxa.

Intrinsic base substitution patterns in diverse species reveal links to cancer and metabolism

Analyses of large-scale cancer sequencing data have revealed that mutagenic processes can create distinctive patterns of base substitutions, called mutational signatures. Interestingly, mutational patterns resembling some of these signatures can also be observed in normal cells. To determine whether similar patterns exist more generally, we analyzed large data sets of genetic variation, including mutations from 7 model species and single nucleotide polymorphisms in 42 species, totaling >1.9 billion variants. We found that base substitution patterns for most species closely match single base substitution (SBS) mutational signature 5 in the Catalog of Somatic Mutations in Cancer (COSMIC) database. SBS5 is ubiquitous in cancers and also present in normal human cells, suggesting that similar patterns of genetic variation across so many species are likely due to conserved biochemistry. We investigated the mechanistic origins of the SBS5-like mutational pattern in Saccharomyces cerevisiae, and show that translesion DNA synthesis and sugar metabolism are directly linked to this form of mutagenesis. We propose that conserved metabolic processes in cells are coupled to continuous generation of genetic variants, which can be acted upon by selection to drive the evolution of biological entities.

Recent insights into the evolution of mutation rates in yeast

Mutation is the origin of all genetic variation, good and bad. The mutation process can evolve in response to mutations, positive or negative selection, and genetic drift, but how these forces contribute to mutation-rate variation is an unsolved problem at the heart of genetics research. Mutations can be challenging to measure, but genome sequencing and other tools have allowed for the collection of larger and more detailed datasets, particularly in the yeast-model system. We review key hypotheses for the evolution of mutation rates and describe recent advances in understanding variation in mutational properties within and among yeast species. The multidimensional spectrum of mutations is increasingly recognized as holding valuable clues about how this important process evolves.

Selection in a growing colony biases results of mutation accumulation experiments

Mutations provide the raw material for natural selection to act. Therefore, understanding the variety and relative frequency of different type of mutations is critical to understanding the nature of genetic diversity in a population. Mutation accumulation (MA) experiments have been used in this context to estimate parameters defining mutation rates, distribution of fitness effects (DFE), and spectrum of mutations. MA experiments can be performed with different effective population sizes. In MA experiments with bacteria, a single founder is grown to a size of a colony (~ 108). It is assumed that natural selection plays a minimal role in dictating the dynamics of colony growth. In this work, we simulate colony growth via a mathematical model, and use our model to mimic an MA experiment. We demonstrate that selection ensures that, in an MA experiment, fraction of all mutations that are beneficial is over-represented by a factor of almost two, and that the distribution of fitness effects of beneficial and deleterious mutations are inaccurately captured in an MA experiment. Given this, the estimate of mutation rates from MA experiments is non-trivial. We then perform an MA experiment with 160 lines of E. coli, and show that due to the effect of selection in a growing colony, the size and sector of a colony from which the experiment is propagated impacts the results. Overall, we demonstrate that the results of MA experiments need to be revisited taking into account the action of selection in a growing colony.

Is the mutation rate lower in genomic regions of stronger selective constraints?

A study of the plant Arabidopsis thaliana detected lower mutation rates in genomic regions where mutations are more likely to be deleterious, challenging the principle that mutagenesis is blind to its consequence. To examine the generality of this finding, we analyze large mutational data from baker's yeast and humans. The yeast data do not exhibit this trend, whereas the human data show an opposite trend that disappears upon the control of potential confounders. We find that the Arabidopsis study identified substantially more mutations than reported in the original data-generating studies and expected from Arabidopsis' mutation rate. These extra mutations are enriched in polynucleotide tracts and have relatively low sequencing qualities so are likely sequencing errors. Furthermore, the polynucleotide “mutations” can produce the purported mutational trend in Arabidopsis. Together, our results do not support lower mutagenesis of genomic regions of stronger selective constraints in the plant, fungal, and animal models examined.

NIH4215: A mutation-prone thiamine auxotrophic clinical Aspergillus fumigatus isolate

Aspergillus fumigatus is the main cause of life-threatening invasive aspergillosis. Despite the availability of various antifungals, therapy remains challenging and requires further studies. Accordingly, the clinical A. fumigatus isolate NIH4215 deriving from a fatal case of human pulmonary aspergillosis has frequently been used in drug efficacy studies. Unexpectedly, our initial attempts to generate a bioluminescent reporter of strain NIH4215 for in vivo drug efficacy studies failed, as NIH4215 was unable to grow on defined minimal medium. Subsequent analyses discovered a previously undescribed thiamine auxotrophy of strain NIH4215 and transformation with thiamine biosynthesis genes from A. fumigatus strain Af293 identified the nmt1 gene as cause of the thiamine auxotrophy. Sequencing of the defective nmt1 gene revealed the loss of a cysteine codon within an essential iron-binding motif. Subsequently, the wild-type nmt1 gene was successfully used to generate a bioluminescent reporter strain in NIH4215 by simultaneously deleting the akuB locus. The resulting bioluminescent ΔakuB strains showed a high frequency of homologous integration as confirmed by generation of pyrG and niaD deletion mutants. When tested in a Galleria mellonella infection model, neither thiamine auxotrophy nor the deletion of the akuB locus had a significant effect on virulence. However, besides thiamine auxotrophy, sectors with altered morphology and albino mutants frequently arose on colony edges of strain NIH4215 and its derivatives, and stable albino mutants were successfully isolated. A proposed increased mutation rate of NIH4215 was confirmed by screening for spontaneous occurrence of fluoorotic acid resistant mutants. Independent mutations in the pyrG and pyrE gene were identified in the fluoroorotic acid resistant NIH4215 isolates and the frequency of mutation was by at least one order of magnitude higher than that observed for the clinical A. fumigatus isolate CBS144.89. In summary, despite its virulence in animal models, strain NIH4215 is a thiamine auxotroph and prone to accumulate mutations. Our results suggest that thiamine biosynthesis is dispensable for host infection and mutation-prone strains such as NIH4215 could potentially facilitate the evolution of azole resistant strains as increasingly observed in the environment.

Genetically controlled mtDNA deletions prevent ROS damage by arresting oxidative phosphorylation

Deletion of mitochondrial DNA in eukaryotes is currently attributed to rare accidental events associated with mitochondrial replication or repair of double-strand breaks. We report the discovery that yeast cells arrest harmful intramitochondrial superoxide production by shutting down respiration through genetically controlled deletion of mitochondrial oxidative phosphorylation genes. We show that this process critically involves the antioxidant enzyme superoxide dismutase 2 and two-way mitochondrial-nuclear communication through Rtg2 and Rtg3. While mitochondrial DNA homeostasis is rapidly restored after cessation of a short-term superoxide stress, long-term stress causes maladaptive persistence of the deletion process, leading to complete annihilation of the cellular pool of intact mitochondrial genomes and irrevocable loss of respiratory ability. This shows that oxidative stress-induced mitochondrial impairment may be under strict regulatory control. If the results extend to human cells, the results may prove to be of etiological as well as therapeutic importance with regard to age-related mitochondrial impairment and disease.

Ribodysgenesis: sudden genome instability in the yeast Saccharomyces cerevisiae arising from RNase H2 cleavage at genomic-embedded ribonucleotides

Ribonucleotides can be incorporated into DNA during replication by the replicative DNA polymerases. These aberrant DNA subunits are efficiently recognized and removed by Ribonucleotide Excision Repair, which is initiated by the heterotrimeric enzyme RNase H2. While RNase H2 is essential in higher eukaryotes, the yeast Saccharomyces cerevisiae can survive without RNase H2 enzyme, although the genome undergoes mutation, recombination and other genome instability events at an increased rate. Although RNase H2 can be considered as a protector of the genome from the deleterious events that can ensue from recognition and removal of embedded ribonucleotides, under conditions of high ribonucleotide incorporation and retention in the genome in a RNase H2-negative strain, sudden introduction of active RNase H2 causes massive DNA breaks and genome instability in a condition which we term 'ribodysgenesis'. The DNA breaks and genome instability arise solely from RNase H2 cleavage directed to the ribonucleotide-containing genome. Survivors of ribodysgenesis have massive loss of heterozygosity events stemming from recombinogenic lesions on the ribonucleotide-containing DNA, with increases of over 1000X from wild-type. DNA breaks are produced over one to two divisions and subsequently cells adapt to RNase H2 and ribonucleotides in the genome and grow with normal levels of genome instability.