Genome-wide mapping of spontaneous genetic alterations in diploid yeast cells

Large-scale genomic rearrangements boost SCRaMbLE in Saccharomyces cerevisiae

Synthetic Chromosome Rearrangement and Modification by LoxP-mediated Evolution (SCRaMbLE) is a promising tool to study genomic rearrangements. However, the potential of SCRaMbLE to study genomic rearrangements is currently hindered, because a strain containing all 16 synthetic chromosomes is not yet available. Here, we construct SparLox83R, a yeast strain containing 83 loxPsym sites distributed across all 16 chromosomes. SCRaMbLE of SparLox83R produces versatile genome-wide genomic rearrangements, including inter-chromosomal events. Moreover, when combined with synthetic chromosomes, SCRaMbLE of hetero-diploids with SparLox83R leads to increased diversity of genomic rearrangements and relatively faster evolution of traits compared to hetero-diploids only with wild-type chromosomes. Analysis of the SCRaMbLEd strain with increased tolerance to nocodazole demonstrates that genomic rearrangements can perturb the transcriptome and 3D genome structure and consequently impact phenotypes. In summary, a genome with sparsely distributed loxPsym sites can serve as a powerful tool for studying the consequence of genomic rearrangements and accelerating strain engineering in Saccharomyces cerevisiae.

Splitting the yeast centromere by recombination

Although fusions between the centromeres of different human chromosomes have been observed cytologically in cancer cells, since the centromeres are long arrays of satellite sequences, the details of these fusions have been difficult to investigate. We developed methods of detecting recombination within the centromeres of the yeast Saccharomyces cerevisiae (intercentromere recombination). These events occur at similar rates (about 10-8/cell division) between two active or two inactive centromeres. We mapped the breakpoints of most of the recombination events to a region of 43 base pairs of uninterrupted homology between the two centromeres. By whole-genome DNA sequencing, we showed that most (>90%) of the events occur by non-reciprocal recombination (gene conversion/break-induced replication). We also found that intercentromere recombination can involve non-homologous chromosome, generating whole-arm translocations. In addition, intercentromere recombination is associated with very frequent chromosome missegregation. These observations support the conclusion that intercentromere recombination generally has negative genetic consequences.

Novel insights into the effects of 5-hydroxymethfurural on genomic instability and phenotypic evolution using a yeast model

5-Hydroxymethfurural (5-HMF) is naturally found in a variety of foods and beverages and represents a main inhibitor in the lignocellulosic hydrolysates used for fermentation. This study investigated the impact of 5-HMF on the genomic stability and phenotypic plasticity of the yeast Saccharomyces cerevisiae. Using next-generation sequencing technology, we examined the genomic alterations of diploid S. cerevisiae isolates that were subcultured on a medium containing 1.2 g/L 5-HMF. We found that in 5-HMF-treated cells, the rates of chromosome aneuploidy, large deletions/duplications, and loss of heterozygosity were elevated compared with that in untreated cells. 5-HMF exposure had a mild impact on the rate of point mutations but altered the mutation spectrum. Contrary to what was observed in untreated cells, more monosomy than trisomy occurred in 5-HMF-treated cells. The aneuploidy mutant with monosomic chromosome IX was more resistant to 5-HMF than the diploid parent strain because of the enhanced activity of alcohol dehydrogenase. Finally, we found that overexpression of ADH6 and ZWF1 effectively stabilized the yeast genome under 5-HMF stress. Our findings not only elucidated the global effect of 5-HMF on the genomic integrity of yeast but also provided novel insights into how chromosomal instability drives the environmental adaptability of eukaryotic cells.IMPORTANCESingle-cell microorganisms are exposed to a range of stressors in both natural and industrial settings. This study investigated the effects of 5-hydroxymethfurural (5-HMF), a major inhibitor found in baked foods and lignocellulosic hydrolysates, on the chromosomal instability of yeast. We examined the mechanisms leading to the distinct patterns of 5-HMF-induced genomic alterations and discovered that chromosomal loss, typically viewed as detrimental to cell growth under most conditions, can contribute to yeast tolerance to 5-HMF. Our results increased the understanding of how specific stressors stimulate genomic plasticity and environmental adaptation in yeast.

Molecular underpinnings and environmental drivers of loss of heterozygosity in Drosophila intestinal stem cells

During development and aging, genome mutation leading to loss of heterozygosity (LOH) can uncover recessive phenotypes within tissue compartments. This phenomenon occurs in normal human tissues and is prevalent in pathological genetic conditions and cancers. While studies in yeast have defined DNA repair mechanisms that can promote LOH, the predominant pathways and environmental triggers in somatic tissues of multicellular organisms are not well understood. Here, we investigate mechanisms underlying LOH in intestinal stem cells in Drosophila. Infection with the pathogenic bacteria, Erwinia carotovora carotovora 15, but not Pseudomonas entomophila, increases LOH frequency. Using whole genome sequencing of somatic LOH events, we demonstrate that they arise primarily via mitotic recombination. Molecular features and genetic evidence argue against a break-induced replication mechanism and instead support cross-over via double Holliday junction-based repair. This study provides a mechanistic understanding of mitotic recombination, an important mediator of LOH, and its effects on stem cells in vivo.

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.

Ploidy evolution in a wild yeast is linked to an interaction between cell type and metabolism

Ploidy is an evolutionarily labile trait, and its variation across the tree of life has profound impacts on evolutionary trajectories and life histories. The immediate consequences and molecular causes of ploidy variation on organismal fitness are frequently less clear, although extreme mating type skews in some fungi hint at links between cell type and adaptive traits. Here, we report an unusual recurrent ploidy reduction in replicate populations of the budding yeast Saccharomyces eubayanus experimentally evolved for improvement of a key metabolic trait, the ability to use maltose as a carbon source. We find that haploids have a substantial, but conditional, fitness advantage in the absence of other genetic variation. Using engineered genotypes that decouple the effects of ploidy and cell type, we show that increased fitness is primarily due to the distinct transcriptional program deployed by haploid-like cell types, with a significant but smaller contribution from absolute ploidy. The link between cell-type specification and the carbon metabolism adaptation can be traced to the noncanonical regulation of a maltose transporter by a haploid-specific gene. This study provides novel mechanistic insight into the molecular basis of an environment-cell type fitness interaction and illustrates how selection on traits unexpectedly linked to ploidy states or cell types can drive karyotypic evolution in fungi.

Origin of fungal hybrids with pathogenic potential from warm seawater environments

Hybridisation is a common event in yeasts often leading to genomic variability and adaptation. The yeast Candida orthopsilosis is a human-associated opportunistic pathogen belonging to the Candida parapsilosis species complex. Most C. orthopsilosis clinical isolates are hybrids resulting from at least four independent crosses between two parental lineages, of which only one has been identified. The rare presence or total absence of parentals amongst clinical isolates is hypothesised to be a consequence of a reduced pathogenicity with respect to their hybrids. Here, we sequence and analyse the genomes of environmental C. orthopsilosis strains isolated from warm marine ecosystems. We find that a majority of environmental isolates are hybrids, phylogenetically closely related to hybrid clinical isolates. Furthermore, we identify the missing parental lineage, thus providing a more complete overview of the genomic evolution of this species. Additionally, we discover phenotypic differences between the two parental lineages, as well as between parents and hybrids, under conditions relevant for pathogenesis. Our results suggest a marine origin of C. orthopsilosis hybrids, with intrinsic pathogenic potential, and pave the way to identify pre-existing environmental adaptations that rendered hybrids more prone than parental lineages to colonise and infect the mammalian host.

Contribution of Spontaneous Mutations to Quantitative and Molecular Variation at the Highly Repetitive rDNA Locus in Yeast

The ribosomal DNA array in Saccharomyces cerevisiae consists of many tandem repeats whose copy number is believed to be functionally important but highly labile. Regulatory mechanisms have evolved to maintain copy number by directed mutation, but how spontaneous variation at this locus is generated and selected has not been well characterized. We applied a mutation accumulation approach to quantify the impacts of mutation and selection on this unique genomic feature across hundreds of mutant strains. We find that mutational variance for this trait is relatively high, and that unselected mutations elsewhere in the genome can disrupt copy number maintenance. In consequence, copy number generally declines gradually, consistent with a previously proposed model of rDNA maintenance where a downward mutational bias is normally compensated by mechanisms that increase copy number when it is low. This pattern holds across ploidy levels and strains in the standard lab environment but differs under some stressful conditions. We identify several alleles, gene categories, and genomic features that likely affect copy number, including aneuploidy for chromosome XII. Copy number change is associated with reduced growth in diploids, consistent with stabilizing selection. Levels of standing variation in copy number are well predicted by a balance between mutation and stabilizing selection, suggesting this trait is not subject to strong diversifying selection in the wild. The rate and spectrum of point mutations within the rDNA locus itself are distinct from the rest of the genome and predictive of polymorphism locations. Our findings help differentiate the roles of mutation and selection and indicate that spontaneous mutation patterns shape several aspects of ribosomal DNA evolution.

High quality de novo genome assembly of the non-conventional yeast Kazachstania bulderi describes a potential low pH production host for biorefineries

Kazachstania bulderi is a non-conventional yeast species able to grow efficiently on glucose and δ-gluconolactone at low pH. These unique traits make K. bulderi an ideal candidate for use in sustainable biotechnology processes including low pH fermentations and the production of green chemicals including organic acids. To accelerate strain development with this species, detailed information of its genetics is needed. Here, by employing long read sequencing we report a high-quality phased genome assembly for three strains of K. bulderi species, including the type strain. The sequences were assembled into 12 chromosomes with a total length of 14 Mb, and the genome was fully annotated at structural and functional levels, including allelic and structural variants, ribosomal array and mating type locus. This high-quality reference genome provides a resource to advance our fundamental knowledge of biotechnologically relevant non-conventional yeasts and to support the development of genetic tools for manipulating such strains towards their use as production hosts in biotechnological processes.

Nonlethal Furfural Exposure Causes Genomic Alterations and Adaptability Evolution in Saccharomyces cerevisiae

Furfural is a major inhibitor found in lignocellulosic hydrolysate, a promising feedstock for the biofermentation industry. In this study, we aimed to investigate the potential impact of this furan-derived chemical on yeast genome integrity and phenotypic evolution by using genetic screening systems and high-throughput analyses. Our results showed that the rates of aneuploidy, chromosomal rearrangements (including large deletions and duplications), and loss of heterozygosity (LOH) increased by 50-fold, 23-fold, and 4-fold, respectively, when yeast cells were cultured in medium containing a nonlethal dose of furfural (0.6 g/L). We observed significantly different ratios of genetic events between untreated and furfural-exposed cells, indicating that furfural exposure induced a unique pattern of genomic instability. Furfural exposure also increased the proportion of CG-to-TA and CG-to-AT base substitutions among point mutations, which was correlated with DNA oxidative damage. Interestingly, although monosomy of chromosomes often results in the slower growth of yeast under spontaneous conditions, we found that monosomic chromosome IX contributed to the enhanced furfural tolerance. Additionally, terminal LOH events on the right arm of chromosome IV, which led to homozygosity of the SSD1 allele, were associated with furfural resistance. This study sheds light on the mechanisms underlying the influence of furfural on yeast genome integrity and adaptability evolution. IMPORTANCE Industrial microorganisms are often exposed to multiple environmental stressors and inhibitors during their application. This study demonstrates that nonlethal concentrations of furfural in the culture medium can significantly induce genome instability in the yeast Saccharomyces cerevisiae. Notably, furfural-exposed yeast cells displayed frequent chromosome aberrations, indicating the potent teratogenicity of this inhibitor. We identified specific genomic alterations, including monosomic chromosome IX and loss of heterozygosity of the right arm of chromosome IV, that confer furfural tolerance to a diploid S. cerevisiae strain. These findings enhance our understanding of how microorganisms evolve and adapt to stressful environments and offer insights for developing strategies to improve their performance in industrial applications.

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.

Genome-wide contributions of the MutSα- and MutSβ-dependent DNA mismatch repair pathways to the maintenance of genetic stability in S. cerevisiae

The DNA mismatch repair (MMR) system is a major DNA repair system that suppresses inherited and sporadic cancers in humans. In eukaryotes the MutSα-dependent and MutSβ-dependent MMR pathways correct DNA polymerase errors. Here, we investigated these two pathways on a whole-genome level in S. cerevisiae. We found that inactivation of MutSα-dependent MMR by deletion of the MSH6 gene increases the genome-wide mutation rate by ∼17-fold, and loss of MutSβ-dependent MMR via deletion of MSH3 elevates the genome-wide mutation rate by ∼4-fold. We also found that MutSα-dependent MMR does not show a preference for protecting coding or noncoding DNA from mutations, whereas MutSβ-dependent MMR preferentially protects noncoding DNA from mutations. The most frequent mutations in the msh6Δ strain are C>T transitions, whereas 1-6-bp deletions are the most common genetic alterations in the msh3Δ strain. Strikingly, MutSα-dependent MMR is more important than MutSβ-dependent MMR for protection from 1-bp insertions, while MutSβ-dependent MMR has a more critical role in the defense against 1-bp deletions and 2-6-bp indels. We also determined that a mutational signature of yeast MSH6 loss is similar to mutational signatures of human MMR deficiency. Furthermore, our analysis showed that compared to other 5'-NCN-3' trinucleotides, 5'-GCA-3' trinucleotides are at the highest risk of accumulating C>T transitions at the central position in the msh6Δ cells and that the presence of a G/A base at the -1 position is important for the efficient MutSα-dependent suppression of C>T transitions. Our results highlight key differences between the roles of the MutSα-dependent and MutSβ-dependent MMR pathways.

Shuffling the yeast genome using CRISPR/Cas9-generated DSBs that target the transposable Ty1 elements

Although homologous recombination between transposable elements can drive genomic evolution in yeast by facilitating chromosomal rearrangements, the details of the underlying mechanisms are not fully clarified. In the genome of the yeast Saccharomyces cerevisiae, the most common class of transposon is the retrotransposon Ty1. Here, we explored how Cas9-induced double-strand breaks (DSBs) directed to Ty1 elements produce genomic alterations in this yeast species. Following Cas9 induction, we observed a significant elevation of chromosome rearrangements such as deletions, duplications and translocations. In addition, we found elevated rates of mitotic recombination, resulting in loss of heterozygosity. Using Southern analysis coupled with short- and long-read DNA sequencing, we revealed important features of recombination induced in retrotransposons. Almost all of the chromosomal rearrangements reflect the repair of DSBs at Ty1 elements by non-allelic homologous recombination; clustered Ty elements were hotspots for chromosome rearrangements. In contrast, a large proportion (about three-fourths) of the allelic mitotic recombination events have breakpoints in unique sequences. Our analysis suggests that some of the latter events reflect extensive processing of the broken ends produced in the Ty element that extend into unique sequences resulting in break-induced replication. Finally, we found that haploid and diploid strain have different preferences for the pathways used to repair double-stranded DNA breaks. Our findings demonstrate the importance of DNA lesions in retrotransposons in driving genome evolution.

Rates and spectra of de novo structural mutations in Chlamydomonas reinhardtii

Genetic variation originates from several types of spontaneous mutation, including single-nucleotide substitutions, short insertions and deletions (indels), and larger structural changes. Structural mutations (SMs) drive genome evolution and are thought to play major roles in evolutionary adaptation, speciation, and genetic disease, including cancers. Sequencing of mutation accumulation (MA) lines has provided estimates of rates and spectra of single-nucleotide and indel mutations in many species, yet the rate of new SMs is largely unknown. Here, we use long-read sequencing to determine the full mutation spectrum in MA lines derived from two strains (CC-1952 and CC-2931) of the green alga Chlamydomonas reinhardtii The SM rate is highly variable between strains and between MA lines, and SMs represent a substantial proportion of all mutations in both strains (CC-1952 6%; CC-2931 12%). The SM spectra differ considerably between the two strains, with almost all inversions and translocations occurring in CC-2931 MA lines. This variation is associated with heterogeneity in the number and type of active transposable elements (TEs), which comprise major proportions of SMs in both strains (CC-1952 22%; CC-2931 38%). In CC-2931, a Crypton and a previously undescribed type of DNA element have caused 71% of chromosomal rearrangements, whereas in CC-1952, a Dualen LINE is associated with 87% of duplications. Other SMs, notably large duplications in CC-2931, are likely products of various double-strand break repair pathways. Our results show that diverse types of SMs occur at substantial rates, and support prominent roles for SMs and TEs in evolution.

Stress combined with loss of the Candida albicans SUMO protease Ulp2 triggers selection of aneuploidy via a two-step process

A delicate balance between genome stability and instability ensures genome integrity while generating genetic diversity, a critical step for evolution. Indeed, while excessive genome instability is harmful, moderated genome instability can drive adaptation to novel environments by maximising genetic variation. Candida albicans, a human fungal pathogen that colonises different parts of the human body, adapts rapidly and frequently to different hostile host microenvironments. In this organism, the ability to generate large-scale genomic variation is a key adaptative mechanism triggering dangerous infections even in the presence of antifungal drugs. Understanding how fitter novel karyotypes are selected is key to determining how C. albicans and other microbial pathogens establish infections. Here, we identified the SUMO protease Ulp2 as a regulator of C. albicans genome integrity through genetic screening. Deletion of ULP2 leads to increased genome instability, enhanced genome variation and reduced fitness in the absence of additional stress. The combined stress caused by the lack of ULP2 and antifungal drug treatment leads to the selection of adaptive segmental aneuploidies that partially rescue the fitness defects of ulp2Δ/Δ cells. Short and long-read genomic sequencing demonstrates that these novel genotypes are selected via a two-step process leading to the formation of novel chromosomal fragments with breakpoints at microhomology regions and DNA repeats.

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.

An apple somatic mutation of delayed fruit maturation date is primarily caused by a retrotransposon insertion-associated large deletion

Somatic mutations may alter important traits in tree fruits, such as fruit color, size and maturation date. Autumn Gala (AGala), a somatic mutation from apple cultivar Gala, matures 4 weeks later than Gala. To understand the mechanisms underlying the delayed maturation, RNA-seq analyses were conducted with fruit sampled at 13 (Gala) and 16 (AGala) time-points during their growth and development. Weighted gene co-expression network analysis (WGCNA) of 23 372 differentially expressed genes resulted in 25 WGCNA modules. Of these, modules 1 (r = -0.98, P = 2E-21) and 2 (r = -0.52, P = 0.004), which were suppressed in AGala, were correlated with fruit maturation date. Surprisingly, 77 of the 152 member genes in module 1 were harbored in a 2.8-Mb genomic region on chromosome 6 that was deleted and replaced by a 10.7-kb gypsy-like retrotransposon (Gy-36) from chromosome 7 in AGala. Among the 77 member genes, MdACT7 was the most suppressed (by 10.5-fold) in AGala due to a disruptive 2.5-kb insertion in coding sequence. Moreover, MdACT7 is the exclusive apple counterpart of Arabidopsis ACT7 known of essential roles in plant development, and the functional allele MdACT7, which was lost to the deletion in AGala, was associated with early fruit maturation in 268 apple accessions. Overexpressing alleles MdACT7 and Mdact7 in an Arabidopsis act7 line showed that MdACT7 largely rescued its stunted growth and delayed initial flowering while Mdact7 did not. Therefore, the 2.8-Mb hemizygous deletion is largely genetically causal for fruit maturation delay in AGala, and the total loss of MdACT7 might have contributed to the phenotype.

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.

MULTI-SCULPT: Multiplex Integration via Selective, CRISPR-Mediated, Ultralong Pathway Transformation in Yeast for Plant Natural Product Synthesis

Yeast has been a versatile model host for complex and valuable natural product biosynthesis via the reconstruction of heterologous biosynthetic pathways. Recent advances in natural product pathway elucidation have uncovered many large and complicated plant pathways that contain 10-30 genes for the biosynthesis of structurally complex, valuable natural products. However, the ability to reconstruct ultralong pathways efficiently in yeast does not match the increasing demand for valuable plant natural product biomanufacturing. Here, we developed a one-pot, multigene pathway integration method in yeast, named MULTI-SCULPT for multiplex integration via selective, CRISPR-mediated, ultralong pathway transformation. Leveraging multilocus genomic disruption via CRISPR/Cas9, newly developed native and synthetic genetic parts, and fine-tuned gene integration and characterization methods, we managed to integrate 21 DNA inserts that contain a 12-gene plant isoflavone biosynthetic pathway into yeast with a 90-100% success rate in 12 days. This method enables fast and efficient ultralong biosynthetic pathway integration and can allow for the fast iterative integration of even longer pathways in the future. Ultimately, this method will accelerate combinatorial optimization of elucidated plant natural product pathways and accelerate putative pathway characterization heterologously.

Bursts of Genomic Instability Potentiate Phenotypic and Genomic Diversification in Saccharomyces cerevisiae

How microbial cells leverage their phenotypic potential to survive in a changing environment is a complex biological problem, with important implications for pathogenesis and species evolution. Stochastic phenotype switching, a particularly fascinating adaptive approach observed in numerous species across the tree of life, introduces phenotypic diversity into a population through mechanisms which have remained difficult to define. Here we describe our investigations into the mechanistic basis of colony morphology phenotype switching which occurs in populations of a pathogenic isolate of Saccharomyces cerevisiae, YJM311. We observed that clonal populations of YJM311 cells produce variant colonies that display altered morphologies and, using whole genome sequence analysis, discovered that these variant clones harbored an exceptional collection of karyotypes newly altered by de novo structural genomic variations (SVs). Overall, our analyses indicate that copy number alterations, more often than changes in allelic identity, provide the causative basis of this phenotypic variation. Individual variants carried between 1 and 16 de novo copy number variations, most of which were whole chromosomal aneuploidies. Notably, we found that the inherent stability of the diploid YJM311 genome is comparable to that of domesticated laboratory strains, indicating that the collections of SVs harbored by variant clones did not arise by a chronic chromosomal instability (CIN) mechanism. Rather, our data indicate that these variant clones acquired such complex karyotypic configurations simultaneously, during stochastic and transient episodes of punctuated systemic genomic instability (PSGI). Surprisingly, we found that the majority of these highly altered variant karyotypes were propagated with perfect fidelity in long-term passaging experiments, demonstrating that high aneuploidy burdens can often be conducive with prolonged genomic integrity. Together, our results demonstrate that colony morphology switching in YJM311 is driven by a stochastic process in which genome stability and plasticity are integrally coupled to phenotypic heterogeneity. Consequently, this system simultaneously introduces both phenotypic and genomic variation into a population of cells, which can, in turn perpetuate population diversity for many generations thereafter.

Genomic Variation-Mediating Fluconazole Resistance in Yeast

Fungal infections pose a serious and growing threat to public health. These infections can be treated with antifungal drugs by killing hazardous fungi in the body. However, the resistance can develop over time when fungi are exposed to antifungal drugs by generating genomic variations, including mutation, aneuploidy, and loss of heterozygosity. The variations could reduce the binding affinity of a drug to its target or block the pathway through which drugs exert their activity. Here, we review genomic variation-mediating fluconazole resistance in the yeast Candida, with the hope of highlighting the functional consequences of genomic variations for the antifungal resistance.

Neural networks enable efficient and accurate simulation-based inference of evolutionary parameters from adaptation dynamics

The rate of adaptive evolution depends on the rate at which beneficial mutations are introduced into a population and the fitness effects of those mutations. The rate of beneficial mutations and their expected fitness effects is often difficult to empirically quantify. As these 2 parameters determine the pace of evolutionary change in a population, the dynamics of adaptive evolution may enable inference of their values. Copy number variants (CNVs) are a pervasive source of heritable variation that can facilitate rapid adaptive evolution. Previously, we developed a locus-specific fluorescent CNV reporter to quantify CNV dynamics in evolving populations maintained in nutrient-limiting conditions using chemostats. Here, we use CNV adaptation dynamics to estimate the rate at which beneficial CNVs are introduced through de novo mutation and their fitness effects using simulation-based likelihood-free inference approaches. We tested the suitability of 2 evolutionary models: a standard Wright-Fisher model and a chemostat model. We evaluated 2 likelihood-free inference algorithms: the well-established Approximate Bayesian Computation with Sequential Monte Carlo (ABC-SMC) algorithm, and the recently developed Neural Posterior Estimation (NPE) algorithm, which applies an artificial neural network to directly estimate the posterior distribution. By systematically evaluating the suitability of different inference methods and models, we show that NPE has several advantages over ABC-SMC and that a Wright-Fisher evolutionary model suffices in most cases. Using our validated inference framework, we estimate the CNV formation rate at the GAP1 locus in the yeast Saccharomyces cerevisiae to be 10-4.7 to 10-4 CNVs per cell division and a fitness coefficient of 0.04 to 0.1 per generation for GAP1 CNVs in glutamine-limited chemostats. We experimentally validated our inference-based estimates using 2 distinct experimental methods-barcode lineage tracking and pairwise fitness assays-which provide independent confirmation of the accuracy of our approach. Our results are consistent with a beneficial CNV supply rate that is 10-fold greater than the estimated rates of beneficial single-nucleotide mutations, explaining the outsized importance of CNVs in rapid adaptive evolution. More generally, our study demonstrates the utility of novel neural network-based likelihood-free inference methods for inferring the rates and effects of evolutionary processes from empirical data with possible applications ranging from tumor to viral evolution.

Global genomic instability caused by reduced expression of DNA polymerase ε in yeast

SignificanceAlthough most studies of the genetic regulation of genome stability involve an analysis of mutations within the coding sequences of genes required for DNA replication or DNA repair, recent studies in yeast show that reduced levels of wild-type enzymes can also produce a mutator phenotype. By whole-genome sequencing and other methods, we find that reduced levels of the wild-type DNA polymerase ε in yeast greatly increase the rates of mitotic recombination, aneuploidy, and single-base mutations. The observed pattern of genome instability is different from those observed in yeast strains with reduced levels of the other replicative DNA polymerases, Pol α and Pol δ. These observations are relevant to our understanding of cancer and other diseases associated with genetic instability.

Genome-scale patterns in the loss of heterozygosity incidence in Saccharomyces cerevisiae

Former studies have established that loss of heterozygosity can be a key driver of sequence evolution in unicellular eukaryotes and tissues of metazoans. However, little is known about whether the distribution of loss of heterozygosity events is largely random or forms discernible patterns across genomes. To initiate our experiments, we introduced selectable markers to both arms of all chromosomes of the budding yeast. Subsequent extensive assays, repeated over several genetic backgrounds and environments, provided a wealth of information on the genetic and environmental determinants of loss of heterozygosity. Three findings stand out. First, the number of loss of heterozygosity events per unit time was more than 25 times higher for growing than starving cells. Second, loss of heterozygosity was most frequent when regions of homology around a recombination site were identical, about a half-% sequence divergence was sufficient to reduce its incidence. Finally, the density of loss of heterozygosity events was highly dependent on the genome's physical architecture. It was several-fold higher on short chromosomal arms than on long ones. Comparably large differences were seen within a single arm where regions close to a centromere were visibly less affected than regions close, though usually not strictly adjacent, to a telomere. We suggest that the observed uneven distribution of loss of heterozygosity events could have been caused not only by an uneven density of initial DNA damages. Location-depended differences in the mode of DNA repair, or its effect on fitness, were likely to operate as well.

Characteristics and possible mechanisms of formation of microinversions distinguishing human and chimpanzee genomes

Genomic inversions come in various sizes. While long inversions are relatively easy to identify by aligning high-quality genome sequences, unambiguous identification of microinversions is more problematic. Here, using a set of extra stringent criteria to distinguish microinversions from other mutational events, we describe microinversions that occurred after the divergence of humans and chimpanzees. In total, we found 59 definite microinversions that range from 17 to 33 nucleotides in length. In majority of them, human genome sequences matched exactly the reverse-complemented chimpanzee genome sequences, implying that the inverted DNA segment was copied precisely. All these microinversions were flanked by perfect or nearly perfect inverted repeats pointing to their key role in their formation. Template switching at inverted repeats during DNA replication was previously discussed as a possible mechanism for the microinversion formation. However, many of definite microinversions found by us cannot be easily explained via template switching owing to the combination of the short length and imperfect nature of their flanking inverted repeats. We propose a novel, alternative mechanism that involves repair of a double-stranded break within the inverting segment via microhomology-mediated break-induced replication, which can consistently explain all definite microinversion events.

Clonal evolution in serially passaged Cryptococcus neoformans × deneoformans hybrids reveals a heterogenous landscape of genomic change

Cryptococcus neoformans × deneoformans hybrids (also known as serotype AD hybrids) are basidiomycete yeasts that are common in a clinical setting. Like many hybrids, the AD hybrids are largely locked at the F1 stage and are mostly unable to undergo normal meiotic reproduction. However, these F1 hybrids, which display a high (∼10%) sequence divergence are known to genetically diversify through mitotic recombination and aneuploidy, and this diversification may be adaptive. In this study, we evolved a single AD hybrid genotype in six diverse environments by serial passaging and then used genome resequencing of evolved clones to determine evolutionary mechanisms of adaptation. The evolved clones generally increased fitness after passaging, accompanied by an average of 3.3 point mutations, 2.9 loss of heterozygosity (LOH) events, and 0.7 trisomic chromosomes per clone. LOH occurred through nondisjunction of chromosomes, crossing over consistent with break-induced replication, and gene conversion, in that order of prevalence. The breakpoints of these recombination events were significantly associated with regions of the genome with lower sequence divergence between the parents and clustered in sub-telomeric regions, notably in regions that had undergone introgression between the two parental species. Parallel evolution was observed, particularly through repeated homozygosity via nondisjunction, yet there was little evidence of environment-specific parallel change for either LOH, aneuploidy, or mutations. These data show that AD hybrids have both a remarkable genomic plasticity and yet are challenged in the ability to recombine through sequence divergence and chromosomal rearrangements, a scenario likely limiting the precision of adaptive evolution to novel environments.

Mitotic recombination in yeast: what we know and what we don't know

Saccharomyces cerevisiae is at the forefront of defining the major recombination mechanisms/models that repair targeted double-strand breaks during mitosis. Each of these models predicts specific molecular intermediates as well as genetic outcomes. Recent use of single-nucleotide polymorphisms to track the exchange of sequences in recombination products has provided an unprecedented level of detail about the corresponding intermediates and the extents to which different mechanisms are utilized. This approach also has revealed complexities that are not predicted by canonical models, suggesting that modifications to these models are needed. Current data are consistent with the initiation of most inter-homolog spontaneous mitotic recombination events by a double-strand break. In addition, the sister chromatid is preferred over the homolog as a repair template.

Uncovering bleomycin-induced genomic alterations and underlying mechanisms in yeast

Bleomycin (BLM) is a widely used chemotherapeutic drug. BLM-treated cells showed an elevated rate of mutations, but the underlying mechanisms remained unclear. In this study, the global genomic alterations in BLM-treated cells were explored in the yeast Saccharomyces cerevisiae. Using genetic assay and whole-genome sequencing, we found that the mutation rate could be greatly elevated in S. cerevisiae cells that underwent Zeocin^TM (a BLM member) treatment. One-base deletion and T to G substitution at the 5'-GT-3' motif was the most striking signature of Zeocin^TM-induced mutations. This was mainly the result of translesion DNA synthesis involving Rev1 and polymerase ζ. Zeocin^TM treatment led to the frequent loss of heterozygosity and chromosomal rearrangements in the diploid strains. The breakpoints of recombination events were significantly associated with certain chromosomal elements. Lastly, we identified multiple genomic alterations that contributed to BLM resistance in the Zeocin^TM-treated mutants. Overall, this study provides new insights into the genotoxicity and evolutional effects of BLM. Importance Bleomycin is an antitumor antibiotic that can mutate genomic DNA. Using yeast models in combination with genome sequencing, the mutational signatures of Zeocin^TM (a member of the bleomycin family) are disclosed. Translesion-synthesis polymerases are crucial for the viability of Zeocin^TM-treated yeast cells at the sacrifice of a higher mutation rate. We also confirmed that multiple genomic alterations were associated with the improved resistance to Zeocin^TM, providing novel insights into how bleomycin resistance is developed in cells.

A mapping platform for mitotic crossover by single-cell multi-omics

Mitotic crossovers have the potential to cause large-scale genome rearrangements. Here, we describe high-throughput, single-cell, whole-genome sequencing methods for mapping crossovers genome-wide at scale. The methods are generalizable to various eukaryotes and to other end points requiring high-throughput, high-coverage single cell sequencing.

Epigenetic Response of Yarrowia lipolytica to Stress: Tracking Methylation Level and Search for Methylation Patterns via Whole-Genome Sequencing

DNA methylation is a common, but not universal, epigenetic modification that plays an important role in multiple cellular processes. While definitely settled for numerous plant, mammalian, and bacterial species, the genome methylation in different fungal species, including widely studied and industrially-relevant yeast species, Yarrowia lipolytica, is still a matter of debate. In this paper, we report a differential DNA methylation level in the genome of Y. lipolytica subjected to sequential subculturing and to heat stress conditions. To this end, we adopted repeated batch bioreactor cultivations of Y. lipolytica subjected to thermal stress in specific time intervals. To analyze the variation in DNA methylation between stressed and control cultures, we (a) quantified the global DNA methylation status using an immuno-assay, and (b) studied DNA methylation patterns through whole-genome sequencing. Primarily, we demonstrated that 5 mC modification can be detected using a commercial immuno-assay, and that the modifications are present in Y. lipolytica’s genome at ~0.5% 5 mC frequency. On the other hand, we did not observe any changes in the epigenetic response of Y. lipolytica to heat shock (HS) treatment. Interestingly, we identified a general phenomenon of decreased 5 mC level in Y. lipolytica’s genome in the stationary phase of growth, when compared to a late-exponential epigenome. While this study provides an insight into the subculturing stress response and adaptation to the stress at epigenetic level by Y. lipolytica, it also leaves an open question of inability to detect any genomic DNA methylation level (either in CpG context or context-less) through whole-genome sequencing. The results of ONT sequencing, suggesting that 5 mC modification is either rare or non-existent in Y. lipolytica genome, are contradicted with the results of the immunoassay.

Overdominant Mutations Restrict Adaptive Loss of Heterozygosity at Linked Loci

Loss of heterozygosity is a common mode of adaptation in asexual diploid populations. Because mitotic recombination frequently extends the full length of a chromosome arm, the selective benefit of loss of heterozygosity may be constrained by linked heterozygous mutations. In a previous laboratory evolution experiment with diploid yeast, we frequently observed homozygous mutations in the WHI2 gene on the right arm of Chromosome XV. However, when heterozygous mutations arose in the STE4 gene, another common target on Chromosome XV, loss of heterozygosity at WHI2 was not observed. Here, we show that mutations at WHI2 are partially dominant and that mutations at STE4 are overdominant. We test whether beneficial heterozygous mutations at these two loci interfere with one another by measuring loss of heterozygosity at WHI2 over 1,000 generations for ∼300 populations that differed initially only at STE4 and WHI2. We show that the presence of an overdominant mutation in STE4 reduces, but does not eliminate, loss of heterozygosity at WHI2. By sequencing 40 evolved clones, we show that populations with linked overdominant and partially dominant mutations show less parallelism at the gene level, more varied evolutionary outcomes, and increased rates of aneuploidy. Our results show that the degree of dominance and the phasing of heterozygous beneficial mutations can constrain loss of heterozygosity along a chromosome arm, and that conflicts between partially dominant and overdominant mutations can affect evolutionary outcomes.

The fidelity of DNA replication, particularly on GC-rich templates, is reduced by defects of the Fe-S cluster in DNA polymerase δ

Iron-sulfur clusters (4Fe–4S) exist in many enzymes concerned with DNA replication and repair. The contribution of these clusters to enzymatic activity is not fully understood. We identified the MET18 (MMS19) gene of Saccharomyces cerevisiae as a strong mutator on GC-rich genes. Met18p is required for the efficient insertion of iron-sulfur clusters into various proteins. met18 mutants have an elevated rate of deletions between short flanking repeats, consistent with increased DNA polymerase slippage. This phenotype is very similar to that observed in mutants of POL3 (encoding the catalytic subunit of Pol δ) that weaken binding of the iron-sulfur cluster. Comparable mutants of POL2 (Pol ϵ) do not elevate deletions. Further support for the conclusion that met18 strains result in impaired DNA synthesis by Pol δ are the observations that Pol δ isolated from met18 strains has less bound iron and is less processive in vitro than the wild-type holoenzyme.

Loss of heterozygosity results in rapid but variable genome homogenization across yeast genetic backgrounds

The dynamics and diversity of the appearance of genetic variants play an essential role in the evolution of the genome and the shaping of biodiversity. Recent population-wide genome sequencing surveys have highlighted the importance of loss of heterozygosity (LOH) events and have shown that they are a neglected part of the genetic diversity landscape. To assess the extent, variability, and spectrum, we explored the accumulation of LOH events in 169 heterozygous diploid Saccharomyces cerevisiae mutation accumulation lines across nine genetic backgrounds. In total, we detected a large set of 22,828 LOH events across distinct genetic backgrounds with a heterozygous level ranging from 0.1% to 1%. LOH events are very frequent with a rate consistently much higher than the mutation rate, showing their importance for genome evolution. We observed that the interstitial LOH (I-LOH) events, resulting in internal short LOH tracts, were much frequent (n = 19,660) than the terminal LOH (T-LOH) events, that is, tracts extending to the end of the chromosome (n = 3168). However, the spectrum, the rate, and the fraction of the genome under LOH vary across genetic backgrounds. Interestingly, we observed that the more the ancestors were heterozygous, the more they accumulated T-LOH events. In addition, frequent short I-LOH tracts are a signature of the lines derived from hybrids with low spore fertility. Finally, we found lines showing almost complete homozygotization during vegetative progression. Overall, our results highlight that the variable dynamics of the LOH accumulation across distinct genetic backgrounds might lead to rapid differential genome evolution during vegetative growth.

Mitotic recombination between homologous chromosomes drives genomic diversity in diatoms

Diatoms, an evolutionarily successful group of microalgae, display high levels of intraspecific genetic variability in natural populations. However, the contribution of various mechanisms generating such diversity is unknown. Here we estimated the genetic micro-diversity within a natural diatom population and mapped the genomic changes arising within clonally propagated diatom cell cultures. Through quantification of haplotype diversity by next-generation sequencing and amplicon re-sequencing of selected loci, we documented a rapid accumulation of multiple haplotypes accompanied by the appearance of novel protein variants in cell cultures initiated from a single founder cell. Comparison of the genomic changes between mother and daughter cells revealed copy number variation and copy-neutral loss of heterozygosity leading to the fixation of alleles within individual daughter cells. The loss of heterozygosity can be accomplished by recombination between homologous chromosomes. To test this hypothesis, we established an endogenous readout system and estimated that the frequency of interhomolog mitotic recombination was under standard growth conditions 4.2 events per 100 cell divisions. This frequency is increased under environmental stress conditions, including treatment with hydrogen peroxide and cadmium. These data demonstrate that copy number variation and mitotic recombination between homologous chromosomes underlie clonal variability in diatom populations. We discuss the potential adaptive evolutionary benefits of the plastic response in the interhomolog mitotic recombination rate, and we propose that this may have contributed to the ecological success of diatoms.

Stability across the Whole Nuclear Genome in the Presence and Absence of DNA Mismatch Repair

We describe the contribution of DNA mismatch repair (MMR) to the stability of the eukaryotic nuclear genome as determined by whole-genome sequencing. To date, wild-type nuclear genome mutation rates are known for over 40 eukaryotic species, while measurements in mismatch repair-defective organisms are fewer in number and are concentrated on Saccharomyces cerevisiae and human tumors. Well-studied organisms include Drosophila melanogaster and Mus musculus, while less genetically tractable species include great apes and long-lived trees. A variety of techniques have been developed to gather mutation rates, either per generation or per cell division. Generational rates are described through whole-organism mutation accumulation experiments and through offspring–parent sequencing, or they have been identified by descent. Rates per somatic cell division have been estimated from cell line mutation accumulation experiments, from systemic variant allele frequencies, and from widely spaced samples with known cell divisions per unit of tissue growth. The latter methods are also used to estimate generational mutation rates for large organisms that lack dedicated germlines, such as trees and hyphal fungi. Mechanistic studies involving genetic manipulation of MMR genes prior to mutation rate determination are thus far confined to yeast, Arabidopsis thaliana, Caenorhabditis elegans, and one chicken cell line. A great deal of work in wild-type organisms has begun to establish a sound baseline, but far more work is needed to uncover the variety of MMR across eukaryotes. Nonetheless, the few MMR studies reported to date indicate that MMR contributes 100-fold or more to genome stability, and they have uncovered insights that would have been impossible to obtain using reporter gene assays.

The Role of Structural Variation in Adaptation and Evolution of Yeast and Other Fungi

Mutations in DNA can be limited to one or a few nucleotides, or encompass larger deletions, insertions, duplications, inversions and translocations that span long stretches of DNA or even full chromosomes. These so-called structural variations (SVs) can alter the gene copy number, modify open reading frames, change regulatory sequences or chromatin structure and thus result in major phenotypic changes. As some of the best-known examples of SV are linked to severe genetic disorders, this type of mutation has traditionally been regarded as negative and of little importance for adaptive evolution. However, the advent of genomic technologies uncovered the ubiquity of SVs even in healthy organisms. Moreover, experimental evolution studies suggest that SV is an important driver of evolution and adaptation to new environments. Here, we provide an overview of the causes and consequences of SV and their role in adaptation, with specific emphasis on fungi since these have proven to be excellent models to study SV.

Origin, Regulation, and Fitness Effect of Chromosomal Rearrangements in the Yeast Saccharomyces cerevisiae

Chromosomal rearrangements comprise unbalanced structural variations resulting in gain or loss of DNA copy numbers, as well as balanced events including translocation and inversion that are copy number neutral, both of which contribute to phenotypic evolution in organisms. The exquisite genetic assay and gene editing tools available for the model organism Saccharomyces cerevisiae facilitate deep exploration of the mechanisms underlying chromosomal rearrangements. We discuss here the pathways and influential factors of chromosomal rearrangements in S. cerevisiae. Several methods have been developed to generate on-demand chromosomal rearrangements and map the breakpoints of rearrangement events. Finally, we highlight the contributions of chromosomal rearrangements to drive phenotypic evolution in various S. cerevisiae strains. Given the evolutionary conservation of DNA replication and recombination in organisms, the knowledge gathered in the small genome of yeast can be extended to the genomes of higher eukaryotes.