Whole Genome Analysis of 132 Clinical Saccharomyces cerevisiae Strains Reveals Extensive Ploidy Variation

Comparative modeling reveals the molecular determinants of aneuploidy fitness cost in a wild yeast model

Although implicated as deleterious in many organisms, aneuploidy can underlie rapid phenotypic evolution. However, aneuploidy will be maintained only if the benefit outweighs the cost, which remains incompletely understood. To quantify this cost and the molecular determinants behind it, we generated a panel of chromosome duplications in Saccharomyces cerevisiae and applied comparative modeling and molecular validation to understand aneuploidy toxicity. We show that 74%-94% of the variance in aneuploid strains' growth rates is explained by the cumulative cost of genes on each chromosome, measured for single-gene duplications using a genomic library, along with the deleterious contribution of small nucleolar RNAs (snoRNAs) and beneficial effects of tRNAs. Machine learning to identify properties of detrimental gene duplicates provided no support for the balance hypothesis of aneuploidy toxicity and instead identified gene length as the best predictor of toxicity. Our results present a generalized framework for the cost of aneuploidy with implications for disease biology and evolution.

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.

Natural proteome diversity links aneuploidy tolerance to protein turnover

Accessing the natural genetic diversity of species unveils hidden genetic traits, clarifies gene functions and allows the generalizability of laboratory findings to be assessed. One notable discovery made in natural isolates of Saccharomyces cerevisiae is that aneuploidy-an imbalance in chromosome copy numbers-is frequent1,2 (around 20%), which seems to contradict the substantial fitness costs and transient nature of aneuploidy when it is engineered in the laboratory3-5. Here we generate a proteomic resource and merge it with genomic1 and transcriptomic6 data for 796 euploid and aneuploid natural isolates. We find that natural and lab-generated aneuploids differ specifically at the proteome. In lab-generated aneuploids, some proteins-especially subunits of protein complexes-show reduced expression, but the overall protein levels correspond to the aneuploid gene dosage. By contrast, in natural isolates, more than 70% of proteins encoded on aneuploid chromosomes are dosage compensated, and average protein levels are shifted towards the euploid state chromosome-wide. At the molecular level, we detect an induction of structural components of the proteasome, increased levels of ubiquitination, and reveal an interdependency of protein turnover rates and attenuation. Our study thus highlights the role of protein turnover in mediating aneuploidy tolerance, and shows the utility of exploiting the natural diversity of species to attain generalizable molecular insights into complex biological processes.

Advanced genetic techniques in fungal pathogen research

Although fungi have been important model organisms for solving genetic, molecular, and ecological problems, recently, they are also becoming an important source of infectious disease. Despite their high medical burden, fungal pathogens are understudied, and relative to other pathogenic microbes, less is known about how their gene functions contribute to disease. This is due, in part, to a lack of powerful genetic tools to study these organisms. In turn, this has resulted in inappropriate treatments and diagnostics and poor disease management. There are a variety of reasons genetic studies were challenging in pathogenic fungi, but in recent years, most of them have been overcome or advances have been made to circumvent these barriers. In this minireview, we highlight how recent advances in genetic studies in fungal pathogens have resulted in the discovery of important biology and potential new antifungals and have created the tools to comprehensively study these important pathogens.

Comparative modeling reveals the molecular determinants of aneuploidy fitness cost in a wild yeast model

Although implicated as deleterious in many organisms, aneuploidy can underlie rapid phenotypic evolution. However, aneuploidy will only be maintained if the benefit outweighs the cost, which remains incompletely understood. To quantify this cost and the molecular determinants behind it, we generated a panel of chromosome duplications in Saccharomyces cerevisiae and applied comparative modeling and molecular validation to understand aneuploidy toxicity. We show that 74-94% of the variance in aneuploid strains' growth rates is explained by the additive cost of genes on each chromosome, measured for single-gene duplications using a genomic library, along with the deleterious contribution of snoRNAs and beneficial effects of tRNAs. Machine learning to identify properties of detrimental gene duplicates provided no support for the balance hypothesis of aneuploidy toxicity and instead identified gene length as the best predictor of toxicity. Our results present a generalized framework for the cost of aneuploidy with implications for disease biology and evolution.

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.

Assembly collapsing versus heterozygosity oversizing: detection of homokaryotic and heterokaryotic Laccaria trichodermophora strains by hybrid genome assembly

Genome assembly and annotation using short-paired reads is challenging for eukaryotic organisms due to their large size, variable ploidy and large number of repetitive elements. However, the use of single-molecule long reads improves assembly quality (completeness and contiguity), but haplotype duplications still pose assembly challenges. To address the effect of read length on genome assembly quality, gene prediction and annotation, we compared genome assemblers and sequencing technologies with four strains of the ectomycorrhizal fungus Laccaria trichodermophora. By analysing the predicted repertoire of carbohydrate enzymes, we investigated the effects of assembly quality on functional inferences. Libraries were generated using three different sequencing platforms (Illumina Next-Seq, Mi-Seq and PacBio Sequel), and genomes were assembled using single and hybrid assemblies/libraries. Long reads or hybrid assemby resolved the collapsing of repeated regions, but the nuclear heterozygous versions remained unresolved. In dikaryotic fungi, each cell includes two nuclei and each nucleus has differences not only in allelic gene version but also in gene composition and synteny. These heterokaryotic cells produce fragmentation and size overestimation of the genome assembly of each nucleus. Hybrid assembly revealed a wider functional diversity of genomes. Here, several predicted oxidizing activities on glycosyl residues of oligosaccharides and several chitooligosaccharide acetylase activities would have passed unnoticed in short-read assemblies. Also, the size and fragmentation of the genome assembly, in combination with heterozygosity analysis, allowed us to distinguish homokaryotic and heterokaryotic strains isolated from L. trichodermophora fruit bodies.

Species-wide survey of the expressivity and complexity spectrum of traits in yeast

Assessing the complexity and expressivity of traits at the species level is an essential first step to better dissect the genotype-phenotype relationship. As trait complexity behaves dynamically, the classic dichotomy between monogenic and complex traits is too simplistic. However, no systematic assessment of this complexity spectrum has been carried out on a population scale to date. In this context, we generated a large diallel hybrid panel composed of 190 unique hybrids coming from 20 natural isolates representative of the S. cerevisiae genetic diversity. For each of these hybrids, a large progeny of 160 individuals was obtained, leading to a total of 30,400 offspring individuals. Their mitotic growth was evaluated on 38 conditions inducing various cellular stresses. We developed a classification algorithm to analyze the phenotypic distributions of offspring and assess the trait complexity. We clearly found that traits are mainly complex at the population level. On average, we found that 91.2% of cross/trait combinations exhibit high complexity, while monogenic and oligogenic cases accounted for only 4.1% and 4.7%, respectively. However, the complexity spectrum is very dynamic, trait specific and tightly related to genetic backgrounds. Overall, our study provided greater insight into trait complexity as well as the underlying genetic basis of its spectrum in a natural population.

Cellular Stress Impact on Yeast Activity in Biotechnological Processes-A Short Overview

The importance of Saccharomyces cerevisiae yeast cells is known worldwide, as they are the most used microorganisms in biotechnology for bioethanol and biofuel production. Also, they are analyzed and studied for their similar internal biochemical processes to human cells, for a better understanding of cell aging and response to cell stressors. The special ability of S. cerevisiae cells to develop in both aerobic and anaerobic conditions makes this microorganism a viable model to study the transformations and the way in which cellular metabolism is directed to face the stress conditions due to environmental changes. Thus, this review will emphasize the effects of oxidative, ethanol, and osmotic stress and also the physiological and genetic response of stress mitigation in yeast cells.

Mechanisms of Intrinsic Postzygotic Isolation: From Traditional Genic and Chromosomal Views to Genomic and Epigenetic Perspectives

Intrinsic postzygotic isolation typically appears as reduced viability or fertility of interspecific hybrids caused by genetic incompatibilities between diverged parental genomes. Dobzhansky-Muller interactions among individual genes, and chromosomal rearrangements causing problems with chromosome synapsis and recombination in meiosis, have both long been considered as major mechanisms behind intrinsic postzygotic isolation. Recent research has, however, suggested that the genetic basis of intrinsic postzygotic isolation can be more complex and involves, for example, overall divergence of the DNA sequence or epigenetic changes. Here, we review the mechanisms of intrinsic postzygotic isolation from genic, chromosomal, genomic, and epigenetic perspectives across diverse taxa. We provide empirical evidence for these mechanisms, discuss their importance in the speciation process, and highlight questions that remain unanswered.

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.

Copy number variation alters local and global mutational tolerance

Copy number variants (CNVs), duplications and deletions of genomic sequences, contribute to evolutionary adaptation but can also confer deleterious effects and cause disease. Whereas the effects of amplifying individual genes or whole chromosomes (i.e., aneuploidy) have been studied extensively, much less is known about the genetic and functional effects of CNVs of differing sizes and structures. Here, we investigated Saccharomyces cerevisiae (yeast) strains that acquired adaptive CNVs of variable structures and copy numbers following experimental evolution in glutamine-limited chemostats. Although beneficial in the selective environment, CNVs result in decreased fitness compared with the euploid ancestor in rich media. We used transposon mutagenesis to investigate mutational tolerance and genome-wide genetic interactions in CNV strains. We find that CNVs increase mutational target size, confer increased mutational tolerance in amplified essential genes, and result in novel genetic interactions with unlinked genes. We validated a novel genetic interaction between different CNVs and BMH1 that was common to multiple strains. We also analyzed global gene expression and found that transcriptional dosage compensation does not affect most genes amplified by CNVs, although gene-specific transcriptional dosage compensation does occur for ∼12% of amplified genes. Furthermore, we find that CNV strains do not show previously described transcriptional signatures of aneuploidy. Our study reveals the extent to which local and global mutational tolerance is modified by CNVs with implications for genome evolution and CNV-associated diseases, such as cancer.

Comparative genomics of infective Saccharomyces cerevisiae strains reveals their food origin

Fungal infections are less studied than viral or bacterial infections and often more difficult to treat. Saccharomyces cerevisiae is usually identified as an innocuous human-friendly yeast; however, this yeast can be responsible for infections mainly in immunosuppressed individuals. S. cerevisiae is a relevant organism widely used in the food industry. Therefore, the study of food yeasts as the source of clinical infection is becoming a pivotal question for food safety. In this study, we demonstrate that S. cerevisiae strains cause infections to spread mostly from food environments. Phylogenetic analysis, genome structure analysis, and phenotypic characterization showed that the key sources of the infective strains are food products, such as bread and probiotic supplements. We observed that the adaptation to host infection can drive important phenotypic and genomic changes in these strains that could be good markers to determine the source of infection. These conclusions add pivotal evidence to reinforce the need for surveillance of food-related S. cerevisiae strains as potential opportunistic pathogens.

DASH/Dam1 complex mutants stabilize ploidy in histone-humanized yeast by weakening kinetochore-microtubule attachments

Forcing budding yeast to chromatinize their DNA with human histones manifests an abrupt fitness cost. We previously proposed chromosomal aneuploidy and missense mutations as two potential modes of adaptation to histone humanization. Here, we show that aneuploidy in histone-humanized yeasts is specific to a subset of chromosomes that are defined by their centromeric evolutionary origins but that these aneuploidies are not adaptive. Instead, we find that a set of missense mutations in outer kinetochore proteins drives adaptation to human histones. Furthermore, we characterize the molecular mechanism underlying adaptation in two mutants of the outer kinetochore DASH/Dam1 complex, which reduce aneuploidy by suppression of chromosome instability. Molecular modeling and biochemical experiments show that these two mutants likely disrupt a conserved oligomerization interface thereby weakening microtubule attachments. We propose a model through which weakened microtubule attachments promote increased kinetochore-microtubule turnover and thus suppress chromosome instability. In sum, our data show how a set of point mutations evolved in histone-humanized yeasts to counterbalance human histone-induced chromosomal instability through weakening microtubule interactions, eventually promoting a return to euploidy.

ploidyfrost: Reference-free estimation of ploidy level from whole genome sequencing data based on de Bruijn graphs

Polyploidy is ubiquitous and its consequences are complex and variable. A change of ploidy level generally influences genetic diversity and results in morphological, physiological and ecological differences between cells or organisms with different ploidy levels. To avoid cumbersome experiments and take advantage of the less biased information provided by the vast amounts of genome sequencing data, computational tools for ploidy estimation are urgently needed. Until now, although a few such tools have been developed, many aspects of this estimation, such as the requirement of a reference genome, the lack of informative results and objective inferences, and the influence of false positives from errors and repeats, need further improvement. We have developed ploidyfrost, a de Bruijn graph-based method, to estimate ploidy levels from whole genome sequencing data sets without a reference genome. ploidyfrost provides a visual representation of allele frequency distribution generated using the ggplot2 package as well as quantitative results using the Gaussian mixture model. In addition, it takes advantage of colouring information encoded in coloured de Bruijn graphs to analyse multiple samples simultaneously and to flexibly filter putative false positives. We evaluated the performance of ploidyfrost by analysing highly heterozygous or repetitive samples of Cyclocarya paliurus and a complex allooctoploid sample of Fragaria × ananassa. Moreover, we demonstrated that the accuracy of analysis results can be improved by constraining a threshold such as Cramér's V coefficient on variant features, which may significantly reduce the side effects of sequencing errors and annoying repeats on the graphical structure constructed.

Copy Number Variation and Allele Ratio Analysis in Candida albicans Using Whole Genome Sequencing Data

Whole genome sequencing of human fungal pathogens has revolutionized the speed and accuracy in which sequence variants that cause antifungal resistance can be identified. Genome rearrangements resulting in copy number variation (CNV) are a significant source of acquired antifungal drug resistance across diverse fungal species. Some CNVs are transient in nature, while other CNVs are stable and well tolerated even in the absence of antifungal drugs. By visualizing whole genome sequencing read depth as a function of genomic location, CNVs and CNV breakpoints (genomic positions where the copy number changes occur relative to the rest of the genome) are rapidly identified. A similar analysis can be used to visualize allele ratio changes that occur across the genomes of heterozygous fungal species, both in the presence and absence of CNVs. This protocol walks through the bioinformatic analysis of CNVs and allele ratios utilizing free, open-source visualization tools. We provide code to use with an example dataset (matched antifungal drug-sensitive and drug-resistant Candida albicans isolates) and notes on how to expand this protocol to other fungal genomes.

UGDR: a generic pipeline to detect recombined regions in polyploid and complex hybrid yeast genomes

MOTIVATION

In eukaryotes, homologous recombination between the parental genomes frequently occurs during the evolutionary conserved process of meiosis, generating the genetic diversity transmitted by the gametes. The genome-wide determination of the frequency and location of the recombination events can now be efficiently performed by genotyping the offspring's polymorphic markers. However, genotyping recombination in complex hybrid genomes with existing methods remains challenging because of their strain and ploidy specificity and the degree of diversity and complexity of the parental genomes, especially in [Formula: see text] polyploids.

RESULTS

We present UGDR, a pipeline to genotype the polymorphisms of complex hybrid yeast genomes. It is based on optimal mapping strategies of NGS reads, comparative analyses of the allelic ratio variation and read depth coverage. We tested the UGDR pipeline with sequencing reads from recombined hybrid diploid yeast strains and various clinical strains exhibiting different degrees of ploidy. UGDR allows to plot the markers distribution and recombination profile per chromosome.

CONCLUSION

UGDR detects and plots recombination events in haploids and polyploid yeasts, which facilitates the discovery and understanding of the yeast genetic recombination map and identify new out-performing recombinants.

Karyon: a computational framework for the diagnosis of hybrids, aneuploids, and other nonstandard architectures in genome assemblies

BACKGROUND

Recent technological developments have made genome sequencing and assembly highly accessible and widely used. However, the presence in sequenced organisms of certain genomic features such as high heterozygosity, polyploidy, aneuploidy, heterokaryosis, or extreme compositional biases can challenge current standard assembly procedures and result in highly fragmented assemblies. Hence, we hypothesized that genome databases must contain a nonnegligible fraction of low-quality assemblies that result from such type of intrinsic genomic factors.

FINDINGS

Here we present Karyon, a Python-based toolkit that uses raw sequencing data and de novo genome assembly to assess several parameters and generate informative plots to assist in the identification of nonchanonical genomic traits. Karyon includes automated de novo genome assembly and variant calling pipelines. We tested Karyon by diagnosing 35 highly fragmented publicly available assemblies from 19 different Mucorales (Fungi) species.

CONCLUSIONS

Our results show that 10 (28.57%) of the assemblies presented signs of unusual genomic configurations, suggesting that these are common, at least for some lineages within the Fungi.

Clonality, inbreeding, and hybridization in two extremotolerant black yeasts

BACKGROUND

The great diversity of lifestyles and survival strategies observed in fungi is reflected in the many ways in which they reproduce and recombine. Although a complete absence of recombination is rare, it has been reported for some species, among them 2 extremotolerant black yeasts from Dothideomycetes: Hortaea werneckii and Aureobasidium melanogenum. Therefore, the presence of diploid strains in these species cannot be explained as the product of conventional sexual reproduction.

RESULTS

Genome sequencing revealed that the ratio of diploid to haploid strains in both H. werneckii and A. melanogenum is about 2:1. Linkage disequilibrium between pairs of polymorphic loci and a high degree of concordance between the phylogenies of different genomic regions confirmed that both species are clonal. Heterozygosity of diploid strains is high, with several hybridizing genome pairs reaching the intergenomic distances typically seen between different fungal species. The origin of diploid strains collected worldwide can be traced to a handful of hybridization events that produced diploids, which were stable over long periods of time and distributed over large geographic areas.

CONCLUSIONS

Our results, based on the genomes of over 100 strains of 2 black yeasts, show that although they are clonal, they occasionally form stable and highly heterozygous diploid intraspecific hybrids. The mechanism of these apparently rare hybridization events, which are not followed by meiosis or haploidization, remains unknown. Both extremotolerant yeasts, H. werneckii and even more so A. melanogenum, a close relative of the intensely recombining and biotechnologically relevant Aureobasidium pullulans, provide an attractive model for studying the role of clonality and ploidy in extremotolerant fungi.

Saccharomyces cerevisiae as a Model System for Eukaryotic Cell Biology, from Cell Cycle Control to DNA Damage Response

The yeast Saccharomyces cerevisiae has been used for bread making and beer brewing for thousands of years. In addition, its ease of manipulation, well-annotated genome, expansive molecular toolbox, and its strong conservation of basic eukaryotic biology also make it a prime model for eukaryotic cell biology and genetics. In this review, we discuss the characteristics that made yeast such an extensively used model organism and specifically focus on the DNA damage response pathway as a prime example of how research in S. cerevisiae helped elucidate a highly conserved biological process. In addition, we also highlight differences in the DNA damage response of S. cerevisiae and humans and discuss the challenges of using S. cerevisiae as a model system.

Patterns of gene expression in ovaries of sexual vs. asexual lineages of a freshwater snail

Why sexual reproduction is so common when asexual reproduction should be much more efficient and less costly remains an open question in evolutionary biology. Comparisons between otherwise similar sexual and asexual taxa allow us to characterize the genetic architecture underlying asexuality, which can, in turn, illuminate how this reproductive mode transition occurred and the mechanisms by which it is maintained or disrupted. Here, we used transcriptome sequencing to compare patterns of ovarian gene expression between actively reproducing obligately sexual and obligately asexual females from multiple lineages of Potamopyrgus antipodarum, a freshwater New Zealand snail characterized by frequent separate transitions to asexuality and coexistence of otherwise similar sexual and asexual lineages. We also used these sequence data to evaluate whether population history accounts for variation in patterns of gene expression. We found that source population was a major source of gene expression variation, and likely more influential than reproductive mode. This outcome for these common garden-raised snails is strikingly similar to earlier results from field-collected snails. While we did not identify a likely set of candidate genes from expression profiles that could plausibly explain how transitions to asexuality occurred, we identified around 1,000 genes with evidence of differential expression between sexual and asexual reproductive modes, and 21 genes that appear to exhibit consistent expression differences between sexuals and asexuals across genetic backgrounds. This second smaller set of genes provides a good starting point for further exploration regarding a potential role in the transition to asexual reproduction. These results mark the first effort to characterize the causes of asexuality in P. antipodarum, demonstrate the apparently high heritability of gene expression patterns in this species, and hint that for P. antipodarum, transitions to asexuality might not necessarily be strongly associated with broad changes in gene expression.

Intracellular trehalose accumulation via the Agt1 transporter promotes freeze-thaw tolerance in Saccharomyces cerevisiae

AIM

This study is to investigate the use of a constitutively expressed trehalose transport protein to directly control intracellular trehalose levels and protect baker's yeast (Saccharomyces cerevisiae) cells against freeze-thaw stress in vivo.

METHODS AND RESULTS

We used a constitutively overexpressed Agt1 transporter system to investigate the role of trehalose in the freeze-thaw tolerance of yeast cells by regulating intracellular trehalose concentrations independently of intracellular biosynthesis. Using this method, we found that increasing intracellular trehalose in yeast cells improved cell survival rate after 8 days of freezing at -80°C and -20°C. We also observed that freeze-thaw tolerance promoted by intracellular trehalose only occurs in highly concentrated cell pellets rather than cells in liquid suspension.

CONCLUSIONS

Trehalose is sufficient to provide freeze-thaw tolerance using our Agt1 overexpression system. Freeze-thaw tolerance can be further enhanced by deletion of genes encoding intracellular trehalose degradation enzymes.

SIGNIFICANCE AND IMPACT OF STUDY

These findings are relevant to improving the freeze-thaw tolerance of baker's yeast in the frozen baked goods industry through engineering strains that can accumulate intracellular trehalose via a constitutively expressed trehalose transporter and inclusion of trehalose into the growth medium.

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.

Heme Oxygenase-1 (HMX1) Loss of Function Increases the In-Host Fitness of the Saccharomyces 'boulardii' Probiotic Yeast in a Mouse Fungemia Model

The use of yeast-containing probiotics is on the rise; however, these products occasionally cause fungal infections and possibly even fungemia among susceptible probiotic-treated patients. The incidence of such cases is probably underestimated, which is why it is important to delve deeper into the pathomechanism and the adaptive features of S. ‘boulardii’. Here in this study, the potential role of the gene heme oxygenase-1 (HMX1) in probiotic yeast bloodstream-derived infections was studied by generating marker-free HMX1 deletion mutants with CRISPR/Cas9 technology from both commercial and clinical S. ‘boulardii’ isolates. The six commercial and clinical yeasts used here represented closely related but different genetic backgrounds as revealed by comparative genomic analysis. We compared the wild-type isolates against deletion mutants for their tolerance of iron starvation, hemolytic activity, as well as kidney burden in immunosuppressed BALB/c mice after lateral tail vein injection. Our results reveal that the lack of HMX1 in S. ‘boulardii’ significantly (p < 0.0001) increases the kidney burden of the mice in most genetic backgrounds, while at the same time causes decreased growth in iron-deprived media in vitro. These findings indicate that even a single-gene loss-of-function mutation can, surprisingly, cause elevated fitness in the host during an opportunistic systemic infection. Our findings indicate that the safety assessment of S. ‘boulardii’ strains should not only take strain-to-strain variation into account, but also avoid extrapolating in vitro results to in vivo virulence factor determination.

Overdominant and partially dominant mutations drive clonal adaptation in diploid Saccharomyces cerevisiae

Identification of adaptive targets in experimental evolution typically relies on extensive replication and genetic reconstruction. An alternative approach is to directly assay all mutations in an evolved clone by generating pools of segregants that contain random combinations of evolved mutations. Here, we apply this method to six Saccharomyces cerevisiae clones isolated from four diploid populations that were clonally evolved for 2,000 generations in rich glucose medium. Each clone contains 17-26 mutations relative to the ancestor. We derived intermediate genotypes between the founder and the evolved clones by bulk mating sporulated cultures of the evolved clones to a barcoded haploid version of the ancestor. We competed the resulting barcoded diploids en masse and quantified fitness in the experimental and alternative environments by barcode sequencing. We estimated average fitness effects of evolved mutations using barcode-based fitness assays and whole genome sequencing for a subset of segregants. In contrast to our previous work with haploid evolved clones, we find that diploids carry fewer beneficial mutations, with modest fitness effects (up to 5.4%) in the environment in which they arose. In agreement with theoretical expectations, reconstruction experiments show that all mutations with a detectable fitness effect manifest some degree of dominance over the ancestral allele, and most are overdominant. Genotypes with lower fitness effects in alternative environments allowed us to identify conditions that drive adaptation in our system.

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.

Gene Amplification as a Mechanism of Yeast Adaptation to Nonsense Mutations in Release Factor Genes

Protein synthesis (translation) is one of the fundamental processes occurring in the cells of living organisms. Translation can be divided into three key steps: initiation, elongation, and termination. In the yeast Saccharomyces cerevisiae, there are two translation termination factors, eRF1 and eRF3. These factors are encoded by the SUP45 and SUP35 genes, which are essential; deletion of any of them leads to the death of yeast cells. However, viable strains with nonsense mutations in both the SUP35 and SUP45 genes were previously obtained in several groups. The survival of such mutants clearly involves feedback control of premature stop codon readthrough; however, the exact molecular basis of such feedback control remain unclear. To investigate the genetic factors supporting the viability of these SUP35 and SUP45 nonsense mutants, we performed whole-genome sequencing of strains carrying mutant sup35-n and sup45-n alleles; while no common SNPs or indels were found in these genomes, we discovered a systematic increase in the copy number of the plasmids carrying mutant sup35-n and sup45-n alleles. We used the qPCR method which confirmed the differences in the relative number of SUP35 and SUP45 gene copies between strains carrying wild-type or mutant alleles of SUP35 and SUP45 genes. Moreover, we compare the number of copies of the SUP35 and SUP45 genes in strains carrying different nonsense mutant variants of these genes as a single chromosomal copy. qPCR results indicate that the number of mutant gene copies is increased compared to the wild-type control. In case of several sup45-n alleles, this was due to a disomy of the entire chromosome II, while for the sup35-218 mutation we observed a local duplication of a segment of chromosome IV containing the SUP35 gene. Taken together, our results indicate that gene amplification is a common mechanism of adaptation to nonsense mutations in release factor genes in yeast.

Dynamic genome plasticity during unisexual reproduction in the human fungal pathogen Cryptococcus deneoformans

Genome copy number variation occurs during each mitotic and meiotic cycle and it is crucial for organisms to maintain their natural ploidy. Defects in ploidy transitions can lead to chromosome instability, which is a hallmark of cancer. Ploidy in the haploid human fungal pathogen Cryptococcus neoformans is exquisitely orchestrated and ranges from haploid to polyploid during sexual development and under various environmental and host conditions. However, the mechanisms controlling these ploidy transitions are largely unknown. During C. deneoformans (formerly C. neoformans var. neoformans, serotype D) unisexual reproduction, ploidy increases prior to the onset of meiosis, can be independent from cell-cell fusion and nuclear fusion, and likely occurs through an endoreplication pathway. To elucidate the molecular mechanisms underlying this ploidy transition, we identified twenty cell cycle-regulating genes encoding cyclins, cyclin-dependent kinases (CDK), and CDK regulators. We characterized four cyclin genes and two CDK regulator genes that were differentially expressed during unisexual reproduction and contributed to diploidization. To detect ploidy transition events, we generated a ploidy reporter, called NURAT, which can detect copy number increases via double selection for nourseothricin-resistant, uracil-prototrophic cells. Utilizing this ploidy reporter, we showed that ploidy transition from haploid to diploid can be detected during the early phases of unisexual reproduction. Interestingly, selection for the NURAT reporter revealed several instances of segmental aneuploidy of multiple chromosomes, which conferred azole resistance in some isolates. These findings provide further evidence of ploidy plasticity in fungi with significant biological and public health implications.

Ploidy is an intrinsic fundamental feature of all eukaryotic organisms, and ploidy variation and maintenance are critical to the organism survival and evolution. Fungi exhibit exquisite plasticity in ploidy variation in adaptation to various environmental stresses. For example, the haploid opportunistic human fungal pathogen C. deneoformans can generate diploid blastospores during unisexual reproduction and also forms polyploid titan cells during host infection; however, the mechanisms underlying these ploidy transitions are largely unknown. In this study, we elucidated the genetic regulatory circuitry governing ploidy duplication during C. deneoformans unisexual reproduction through the identification and characterization of cell cycle regulators that are differentially expressed during unisexual reproduction. We showed that four cyclin and two cyclin-dependent kinase regulator genes function in concert to orchestrate ploidy transition during unisexual reproduction. To trace and track ploidy transition events, we also generated a ploidy reporter and revealed the formation of segmental aneuploidy in addition to diploidization, illustrating the diverse mechanisms of genome plasticity in C. deneoformans.

Patterns of Genomic Instability in Interspecific Yeast Hybrids With Diverse Ancestries

The genomes of hybrids often show substantial deviations from the features of the parent genomes, including genomic instabilities characterized by chromosomal rearrangements, gains, and losses. This plastic genomic architecture generates phenotypic diversity, potentially giving hybrids access to new ecological niches. It is however unclear if there are any generalizable patterns and predictability in the type and prevalence of genomic variation and instability across hybrids with different genetic and ecological backgrounds. Here, we analyzed the genomic architecture of 204 interspecific Saccharomyces yeast hybrids isolated from natural, industrial fermentation, clinical, and laboratory environments. Synchronous mapping to all eight putative parental species showed significant variation in read depth indicating frequent aneuploidy, affecting 44% of all hybrid genomes and particularly smaller chromosomes. Early generation hybrids with largely equal genomic content from both parent species were more likely to contain aneuploidies than introgressed genomes with an older hybridization history, which presumably stabilized the genome. Shared k-mer analysis showed that the degree of genomic diversity and variability varied among hybrids with different parent species. Interestingly, more genetically distant crosses produced more similar hybrid genomes, which may be a result of stronger negative epistasis at larger genomic divergence, putting constraints on hybridization outcomes. Mitochondrial genomes were typically inherited from the species also contributing the majority nuclear genome, but there were clear exceptions to this rule. Together, we find reliable genomic predictors of instability in hybrids, but also report interesting cross- and environment-specific idiosyncrasies. Our results are an important step in understanding the factors shaping divergent hybrid genomes and their role in adaptive evolution.

Generation of unequal nuclear genotype proportions in Rhizophagus irregularis progeny causes allelic imbalance in gene transcription

Arbuscular mycorrhizal fungi (AMF) form mutualisms with most plant species. The model AMF Rhizophagus irregularis is common in many ecosystems and naturally forms homokaryons and dikaryons. Quantitative variation in allele frequencies in clonally dikaryon offspring suggests they disproportionately inherit two distinct nuclear genotypes from their parent. This is interesting, because such progeny strongly and differentially affect plant growth. Neither the frequency and magnitude of this occurrence nor its effect on gene transcription are known. Using reduced representation genome sequencing, transcriptomics, and quantitative analysis tools, we show that progeny of homokaryons and dikaryons are qualitatively genetically identical to the parent. However, dikaryon progeny differ quantitatively due to unequal inheritance of nuclear genotypes. Allele frequencies of actively transcribed biallelic genes resembled the frequencies of the two nuclear genotypes. More biallelic genes showed transcription of both alleles than monoallelic transcription, but biallelic transcription was less likely with greater allelic divergence. Monoallelic transcription levels of biallelic genes were reduced compared with biallelic gene transcription, a finding consistent with genomic conflict. Given that genetic variation in R. irregularis is associated with plant growth, our results establish quantitative genetic variation as a future consideration when selecting AMF lines to improve plant production.

Genomic diversity, chromosomal rearrangements, and interspecies hybridization in the Ogataea polymorpha species complex

The methylotrophic yeast Ogataea polymorpha has long been a useful system for recombinant protein production, as well as a model system for methanol metabolism, peroxisome biogenesis, thermotolerance, and nitrate assimilation. It has more recently become an important model for the evolution of mating-type switching. Here, we present a population genomics analysis of 47 isolates within the O. polymorpha species complex, including representatives of the species O. polymorpha, Ogataea parapolymorpha, Ogataea haglerorum, and Ogataea angusta. We found low levels of nucleotide sequence diversity within the O. polymorpha species complex and identified chromosomal rearrangements both within and between species. In addition, we found that one isolate is an interspecies hybrid between O. polymorpha and O. parapolymorpha and present evidence for loss of heterozygosity following hybridization.

High Complexity and Degree of Genetic Variation in Brettanomyces bruxellensis Population

Genome-wide characterization of genetic variants of a large population of individuals within the same species is essential to have a deeper insight into its evolutionary history as well as the genotype–phenotype relationship. Population genomic surveys have been performed in multiple yeast species, including the two model organisms, Saccharomyces cerevisiae and Schizosaccharomyces pombe. In this context, we sought to characterize at the population level the Brettanomyces bruxellensis yeast species, which is a major cause of wine spoilage and can contribute to the specific flavor profile of some Belgium beers. We have completely sequenced the genome of 53 B. bruxellensis strains isolated worldwide. The annotation of the reference genome allowed us to define the gene content of this species. As previously suggested, our genomic data clearly highlighted that genetic diversity variation is related to ploidy level, which is variable in the B. bruxellensis species. Genomes are punctuated by multiple loss-of-heterozygosity regions, whereas aneuploidies as well as segmental duplications are uncommon. Interestingly, triploid genomes are more prone to gene copy number variation than diploids. Finally, the pangenome of the species was reconstructed and was found to be small with few accessory genes compared with S. cerevisiae. The pangenome is composed of 5,409 ORFs (open reading frames) among which 5,106 core ORFs and 303 ORFs that are variable within the population. All these results highlight the different trajectories of species evolution and consequently the interest of establishing population genomic surveys in more populations.

Species-wide transposable element repertoires retrace the evolutionary history of the Saccharomyces cerevisiae host

Transposable elements (TE) are an important source of genetic variation with a dynamic and content that greatly differ in a wide range of species. The origin of the intraspecific content variation is not always clear and little is known about the precise nature of it. Here, we surveyed the species-wide content of the Ty LTR-retrotransposons in a broad collection of 1,011 Saccharomyces cerevisiae natural isolates to understand what can stand behind the variation of the repertoire that is the type and number of Ty elements. We have compiled an exhaustive catalog of all the TE sequence variants present in the S. cerevisiae species by identifying a large set of new sequence variants. The characterization of the TE content in each isolate clearly highlighted that each subpopulation exhibits a unique and specific repertoire, retracing the evolutionary history of the species. Most interestingly, we have shown that ancient interspecific hybridization events had a major impact in the birth of new sequence variants and therefore in the shaping of the TE repertoires. We also investigated the transpositional activity of these elements in a large set of natural isolates, and we found a broad variability related to the level of ploidy as well as the genetic background. Overall, our results pointed out that the evolution of the Ty content is deeply impacted by clade-specific events such as introgressions and therefore follows the population structure. In addition, our study lays the foundation for future investigations to better understand the transpositional regulation and more broadly the TE–host interactions.

The neutral rate of whole-genome duplication varies among yeast species and their hybrids

Hybridization and polyploidization are powerful mechanisms of speciation. Hybrid speciation often coincides with whole-genome duplication (WGD) in eukaryotes. This suggests that WGD may allow hybrids to thrive by increasing fitness, restoring fertility and/or increasing access to adaptive mutations. Alternatively, it has been suggested that hybridization itself may trigger WGD. Testing these models requires quantifying the rate of WGD in hybrids without the confounding effect of natural selection. Here we show, by measuring the spontaneous rate of WGD of more than 1300 yeast crosses evolved under relaxed selection, that some genotypes or combinations of genotypes are more prone to WGD, including some hybrids between closely related species. We also find that higher WGD rate correlates with higher genomic instability and that WGD increases fertility and genetic variability. These results provide evidence that hybridization itself can promote WGD, which in turn facilitates the evolution of hybrids.

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.

How to characterize a strain? Clonal heterogeneity in industrial Saccharomyces influences both phenotypes and heterogeneity in phenotypes

Populations of microbes are constantly evolving heterogeneity that selection acts upon, yet heterogeneity is nontrivial to assess methodologically. The necessary practice of isolating single-cell colonies and thus subclone lineages for establishing, transferring, and using a strain results in single-cell bottlenecks with a generally neglected effect on the characteristics of the strain itself. Here, we present evidence that various subclone lineages for industrial yeasts sequenced for recent genomic studies show considerable differences, ranging from loss of heterozygosity to aneuploidies. Subsequently, we assessed whether phenotypic heterogeneity is also observable in industrial yeast, by individually testing subclone lineages obtained from products. Phenotyping of industrial yeast samples and their newly isolated subclones showed that single-cell bottlenecks during isolation can indeed considerably influence the observable phenotype. Next, we decoupled fitness distributions on the level of individual cells from clonal interference by plating single-cell colonies and quantifying colony area distributions. We describe and apply an approach using statistical modeling to compare the heterogeneity in phenotypes across samples and subclone lineages. One strain was further used to show how individual subclonal lineages are remarkably different not just in phenotype but also in the level of heterogeneity in phenotype. With these observations, we call attention to the fact that choosing an initial clonal lineage from an industrial yeast strain may vastly influence downstream performances and observations on karyotype, on phenotype, and also on heterogeneity.

Genetic variation in aneuploidy prevalence and tolerance across Saccharomyces cerevisiae lineages

Individuals carrying an aberrant number of chromosomes can vary widely in their expression of aneuploidy phenotypes. A major unanswered question is the degree to which an individual's genetic makeup influences its tolerance of karyotypic imbalance. Here we investigated within-species variation in aneuploidy prevalence and tolerance, using Saccharomyces cerevisiae as a model for eukaryotic biology. We analyzed genotypic and phenotypic variation recently published for over 1,000 S. cerevisiae strains spanning dozens of genetically defined clades and ecological associations. Our results show that the prevalence of chromosome gain and loss varies by clade and can be better explained by differences in genetic background than ecology. The relationships between lineages with high aneuploidy frequencies suggest that increased aneuploidy prevalence emerged multiple times in S. cerevisiae evolution. Separate from aneuploidy prevalence, analyzing growth phenotypes revealed that some genetic backgrounds-such as the European Wine lineage-show fitness costs in aneuploids compared to euploids, whereas other clades with high aneuploidy frequencies show little evidence of major deleterious effects. Our analysis confirms that chromosome gain can produce phenotypic benefits, which could influence evolutionary trajectories. These results have important implications for understanding genetic variation in aneuploidy prevalence in health, disease, and evolution.

The population genetics of ploidy change in unicellular fungi

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

Genome-Wide Association Study in New York Phytophthora capsici Isolates Reveals Loci Involved in Mating Type and Mefenoxam Sensitivity

Phytophthora capsici is a soilborne oomycete plant pathogen that causes severe vegetable crop losses in New York (NY) state and worldwide. This pathogen is difficult to manage, in part due to its production of long-lasting sexual spores and its tendency to quickly evolve fungicide resistance. We single nucleotide polymorphism (SNP) genotyped 252 P. capsici isolates, predominantly from NY, in order to conduct a genome-wide association study for mating type and mefenoxam sensitivity. The population structure and extent of chromosomal copy number variation in this collection of isolates were also characterized. Population structure analyses showed isolates largely clustered by the field site where they were collected, with values of F_ST between pairs of fields ranging from 0.10 to 0.31. Thirty-three isolates were putative aneuploids, demonstrating evidence for having up to four linkage groups present in more than two copies, and an additional two isolates appeared to be genome-wide triploids. Mating type was mapped to a region on scaffold 4, consistent with previous findings, and mefenoxam sensitivity was associated with several SNP markers at a novel locus on scaffold 62. We identified several candidate genes for mefenoxam sensitivity, including a homolog of yeast ribosome synthesis factor Rrp5, but failed to locate near the scaffold 62 locus any subunits of RNA polymerase I, the hypothesized target site of phenylamide fungicides in oomycetes. This work expands our knowledge of the population biology of P. capsici and provides a foundation for functional validation of candidate genes associated with epidemiologically important phenotypes.

Genomic instability in an interspecific hybrid of the genus Saccharomyces: a matter of adaptability

Ancient events of polyploidy have been linked to huge evolutionary leaps in the tree of life, while increasing evidence shows that newly established polyploids have adaptive advantages in certain stress conditions compared to their relatives with a lower ploidy. The genus Saccharomyces is a good model for studying such events, as it contains an ancient whole-genome duplication event and many sequenced Saccharomyces cerevisiae are, evolutionary speaking, newly formed polyploids. Many polyploids have unstable genomes and go through large genome erosions; however, it is still unknown what mechanisms govern this reduction. Here, we sequenced and studied the natural S. cerevisiae × Saccharomyces kudriavzevii hybrid strain, VIN7, which was selected for its commercial use in the wine industry. The most singular observation is that its nuclear genome is highly unstable and drastic genomic alterations were observed in only a few generations, leading to a widening of its phenotypic landscape. To better understand what leads to the loss of certain chromosomes in the VIN7 cell population, we looked for genetic features of the genes, such as physical interactions, complex formation, epistatic interactions and stress responding genes, which could have beneficial or detrimental effects on the cell if their dosage is altered by a chromosomal copy number variation. The three chromosomes lost in our VIN7 population showed different patterns, indicating that multiple factors could explain the mechanisms behind the chromosomal loss. However, one common feature for two out of the three chromosomes is that they are among the smallest ones. We hypothesize that small chromosomes alter their copy numbers more frequently as a low number of genes is affected, meaning that it is a by-product of genome instability, which might be the chief driving force of the adaptability and genome architecture of this hybrid.

Fungal evolution: cellular, genomic and metabolic complexity

The question of how phenotypic and genomic complexity are inter-related and how they are shaped through evolution is a central question in biology that historically has been approached from the perspective of animals and plants. In recent years, however, fungi have emerged as a promising alternative system to address such questions. Key to their ecological success, fungi present a broad and diverse range of phenotypic traits. Fungal cells can adopt many different shapes, often within a single species, providing them with great adaptive potential. Fungal cellular organizations span from unicellular forms to complex, macroscopic multicellularity, with multiple transitions to higher or lower levels of cellular complexity occurring throughout the evolutionary history of fungi. Similarly, fungal genomes are very diverse in their architecture. Deep changes in genome organization can occur very quickly, and these phenomena are known to mediate rapid adaptations to environmental changes. Finally, the biochemical complexity of fungi is huge, particularly with regard to their secondary metabolites, chemical products that mediate many aspects of fungal biology, including ecological interactions. Herein, we explore how the interplay of these cellular, genomic and metabolic traits mediates the emergence of complex phenotypes, and how this complexity is shaped throughout the evolutionary history of Fungi.

Defining and Disrupting Species Boundaries in Saccharomyces

The genus Saccharomyces is an evolutionary paradox. On the one hand, it is composed of at least eight clearly phylogenetically delineated species; these species are reproductively isolated from each other, and hybrids usually cannot complete their sexual life cycles. On the other hand, Saccharomyces species have a long evolutionary history of hybridization, which has phenotypic consequences for adaptation and domestication. A variety of cellular, ecological, and evolutionary mechanisms are responsible for this partial reproductive isolation among Saccharomyces species. These mechanisms have caused the evolution of diverse Saccharomyces species and hybrids, which occupy a variety of wild and domesticated habitats. In this article, we introduce readers to the mechanisms isolating Saccharomyces species, the circumstances in which reproductive isolation mechanisms are effective and ineffective, and the evolutionary consequences of partial reproductive isolation. We discuss both the evolutionary history of the genus Saccharomyces and the human history of taxonomists and biologists struggling with species concepts in this fascinating genus.

Virulence phenotypes result from interactions between pathogen ploidy and genetic background

Studying fungal virulence is often challenging and frequently depends on many contexts, including host immune status and pathogen genetic background. However, the role of ploidy has often been overlooked when studying virulence in eukaryotic pathogens. Since fungal pathogens, including the human opportunistic pathogen Candida albicans, can display extensive ploidy variation, assessing how ploidy impacts virulence has important clinical relevance. As an opportunistic pathogen, C. albicans causes nonlethal, superficial infections in healthy individuals, but life-threatening bloodstream infections in individuals with compromised immune function. Here, we determined how both ploidy and genetic background of C. albicans impacts virulence phenotypes in healthy and immunocompromised nematode hosts by characterizing virulence phenotypes in four near-isogenic diploid and tetraploid pairs of strains, which included both laboratory and clinical genetic backgrounds. We found that C. albicans infections decreased host survival and negatively impacted host reproduction, and we leveraged these two measures to survey both lethal and nonlethal virulence phenotypes across the multiple C. albicans strains. In this study, we found that regardless of pathogen ploidy or genetic background, immunocompromised hosts were susceptible to fungal infection compared to healthy hosts. Furthermore, for each host context, we found a significant interaction between C. albicans genetic background and ploidy on virulence phenotypes, but no global differences between diploid and tetraploid pathogens were observed.

Heat shock drives genomic instability and phenotypic variations in yeast

High temperature causes ubiquitous environmental stress to microorganisms, but studies have not fully explained whether and to what extent heat shock would affect genome stability. Hence, this study explored heat-shock-induced genomic alterations in the yeast Saccharomyces cerevisiae. Using genetic screening systems and customized single nucleotide polymorphism (SNP) microarrays, we found that heat shock (52 °C) for several minutes could heighten mitotic recombination by at least one order of magnitude. More than half of heat-shock-induced mitotic recombinations were likely to be initiated by DNA breaks in the S/G₂ phase of the cell cycle. Chromosomal aberration, mainly trisomy, was elevated hundreds of times in heat-shock-treated cells than in untreated cells. Distinct chromosomal instability patterns were also observed between heat-treated and carbendazim-treated yeast cells. Finally, we demonstrated that heat shock stimulates fast phenotypic evolutions (such as tolerance to ethanol, vanillin, fluconazole, and tunicamycin) in the yeast population. This study not only provided novel insights into the effect of temperature fluctuations on genomic integrity but also developed a simple protocol to generate an aneuploidy mutant of yeast.

Candida albicans Genetic Background Influences Mean and Heterogeneity of Drug Responses and Genome Stability during Evolution in Fluconazole

Antimicrobial resistance is an evolutionary phenomenon with clinical implications. We tested how replicates from diverse strains of Candida albicans, a prevalent human fungal pathogen, evolve in the commonly prescribed antifungal drug fluconazole. Replicates on average increased in fitness in the level of drug they were evolved to, with the least fit parental strains improving the most. Very few replicates increased resistance above the drug level they were evolved in. Notably, many replicates increased in genome size and changed in drug tolerance (a drug response where a subpopulation of cells grow slowly in high levels of drug), and variability among replicates in fitness, tolerance, and genome size was higher in strains that initially were more sensitive to the drug. Genetic background influenced the average degree of adaptation and the evolved variability of many phenotypes, highlighting that different strains from the same species may respond and adapt very differently during adaptation.

The importance of within-species diversity in determining the evolutionary potential of a population to evolve drug resistance or tolerance is not well understood, including in eukaryotic pathogens. To examine the influence of genetic background, we evolved replicates of 20 different clinical isolates of Candida albicans, a human fungal pathogen, in fluconazole, the commonly used antifungal drug. The isolates hailed from the major C. albicans clades and had different initial levels of drug resistance and tolerance to the drug. The majority of replicates rapidly increased in fitness in the evolutionary environment, with the degree of improvement inversely correlated with parental strain fitness in the drug. Improvement was largely restricted to up to the evolutionary level of drug: only 4% of the evolved replicates increased resistance (MIC) above the evolutionary level of drug. Prevalent changes were altered levels of drug tolerance (slow growth of a subpopulation of cells at drug concentrations above the MIC) and increased diversity of genome size. The prevalence and predominant direction of these changes differed in a strain-specific manner, but neither correlated directly with parental fitness or improvement in fitness. Rather, low parental strain fitness was correlated with high levels of heterogeneity in fitness, tolerance, and genome size among evolved replicates. Thus, parental strain background is an important determinant in mean improvement to the evolutionary environment as well as the diversity of evolved phenotypes, and the range of possible responses of a pathogen to an antimicrobial drug cannot be captured by in-depth study of a single strain background.

IMPORTANCE Antimicrobial resistance is an evolutionary phenomenon with clinical implications. We tested how replicates from diverse strains of Candida albicans, a prevalent human fungal pathogen, evolve in the commonly prescribed antifungal drug fluconazole. Replicates on average increased in fitness in the level of drug they were evolved to, with the least fit parental strains improving the most. Very few replicates increased resistance above the drug level they were evolved in. Notably, many replicates increased in genome size and changed in drug tolerance (a drug response where a subpopulation of cells grow slowly in high levels of drug), and variability among replicates in fitness, tolerance, and genome size was higher in strains that initially were more sensitive to the drug. Genetic background influenced the average degree of adaptation and the evolved variability of many phenotypes, highlighting that different strains from the same species may respond and adapt very differently during adaptation.

Aneuploidy Enables Cross-Adaptation to Unrelated Drugs

Aneuploidy is common both in tumor cells responding to chemotherapeutic agents and in fungal cells adapting to antifungal drugs. Because aneuploidy simultaneously affects many genes, it has the potential to confer multiple phenotypes to the same cells. Here, we analyzed the mechanisms by which Candida albicans, the most prevalent human fungal pathogen, acquires the ability to survive both chemotherapeutic agents and antifungal drugs. Strikingly, adaptation to both types of drugs was accompanied by the acquisition of specific whole-chromosome aneuploidies, with some aneuploid karyotypes recovered independently and repeatedly from very different drug conditions. Specifically, strains selected for survival in hydroxyurea, an anticancer drug, acquired cross-adaptation to caspofungin, a first-line antifungal drug, and both acquired traits were attributable to trisomy of the same chromosome: loss of trisomy was accompanied by loss of adaptation to both drugs. Mechanistically, aneuploidy simultaneously altered the copy number of most genes on chromosome 2, yet survival in hydroxyurea or caspofungin required different genes and stress response pathways. Similarly, chromosome 5 monosomy conferred increased tolerance to both fluconazole and to caspofungin, antifungals with different mechanisms of action. Thus, the potential for cross-adaptation is not a feature of aneuploidy per se; rather, it is dependent on specific genes harbored on given aneuploid chromosomes. Furthermore, pre-exposure to hydroxyurea increased the frequency of appearance of caspofungin survivors, and hydroxyurea-adapted C. albicans cells were refractory to antifungal drug treatment in a mouse model of systemic candidiasis. This highlights the potential clinical consequences for the management of cancer chemotherapy patients at risk of fungal infections.

Molecular signatures of aneuploidy-driven adaptive evolution

Alteration of normal ploidy (aneuploidy) can have a number of opposing effects, such as unbalancing protein abundances and inhibiting cell growth but also accelerating genetic diversification and rapid adaptation. The interplay of these detrimental and beneficial effects remains puzzling. Here, to understand how cells develop tolerance to aneuploidy, we subject disomic (i.e. with an extra chromosome copy) strains of yeast to long-term experimental evolution under strong selection, by forcing disomy maintenance and daily population dilution. We characterize mutations, karyotype alterations and gene expression changes, and dissect the associated molecular strategies. Cells with different extra chromosomes accumulated mutations at distinct rates and displayed diverse adaptive events. They tended to evolve towards normal ploidy through chromosomal DNA loss and gene expression changes. We identify genes with recurrent mutations and altered expression in multiple lines, revealing a variant that improves growth under genotoxic stresses. These findings support rapid evolvability of disomic strains that can be used to characterize fitness effects of mutations under different stress conditions.

Aneuploidy (abnormal chromosome number) can enable rapid adaptation to stress conditions, but it also entails fitness costs from gene imbalance. Here, the authors experimentally evolve yeast while forcing maintenance of aneuploidy to identify the mechanisms that promote tolerance of aneuploidy.

A Saccharomyces paradox: chromosomes from different species are incompatible because of anti-recombination, not because of differences in number or arrangement

Many species are able to hybridize, but the sterility of these hybrids effectively prevents gene flow between the species, reproductively isolating them and allowing them to evolve independently. Yeast hybrids formed by Saccharomyces cerevisiae and Saccharomyces paradoxus parents are viable and able to grow by mitosis, but they are sexually sterile because most of the gametes they make by meiosis are inviable. The genomes of these two species are so diverged that they cannot recombine properly during meiosis, so they fail to segregate efficiently. Thus most hybrid gametes are inviable because they lack essential chromosomes. Recent work shows that chromosome mis-segregation explains nearly all observed hybrid sterility—genetic incompatibilities have only a small sterilising effect, and there are no significant sterilising incompatibilities in chromosome arrangement or number between the species. It is interesting that chromosomes from these species have diverged so much in sequence without changing in configuration, even though large chromosomal changes occur quite frequently, and sometimes beneficially, in evolving yeast populations.

Baker's Yeast Clinical Isolates Provide a Model for How Pathogenic Yeasts Adapt to Stress

Global outbreaks of drug-resistant fungi such as Candida auris are thought to be due at least in part to excessive use of antifungal drugs. Baker's yeast Saccharomyces cerevisiae has gained importance as an emerging opportunistic fungal pathogen that can cause infections in immunocompromised patients. Analyses of over 1000 S. cerevisiae isolates are providing rich resources to better understand how fungi can grow in human environments. A large percentage of clinical S. cerevisiae isolates are heterozygous across many nucleotide sites, and a significant proportion are of mixed ancestry and/or are aneuploid or polyploid. Such features potentially facilitate adaptation to new environments. These observations provide strong impetus for expanding genomic and molecular studies on clinical and wild isolates to understand the prevalence of genetic diversity and instability-generating mechanisms, and how they are selected for and maintained. Such work can also lead to the identification of new targets for antifungal drugs.

The life and times of yeasts in traditional food fermentations

Yeasts are eukaryotic microorganisms which have a long history in the biotechnology of food production, as they have been used since centuries in bread-making or in the production of alcoholic beverages such as wines or beers. Relative to this importance, a lot of research has been devoted to the study of yeasts involved in making these important products. The role of yeasts in other fermentations in association with other microorganisms - mainly lactic acid bacteria - has been relatively less studied, and often it is not clear if yeasts occurring in such fermentations are contaminants with no role in the fermentation, spoilage microorganisms or whether they actually serve a technological or functional purpose. Some knowledge is available for yeasts used as starter cultures in fermented raw sausages or in the production of acid curd cheeses. This review aimed to summarize the current knowledge on the taxonomy, the presence and potential functional or technological roles of yeasts in traditional fermented plant, dairy, fish and meat fermentations.