Genomic factors shape carbon and nitrogen metabolic niche breadth across Saccharomycotina yeasts.

Organisms exhibit extensive variation in ecological niche breadth, from very narrow (specialists) to very broad (generalists). Two general paradigms have been proposed to explain this variation: (i) trade-offs between performance efficiency and breadth and (ii) the joint influence of extrinsic (environmental) and intrinsic (genomic) factors. We assembled genomic, metabolic, and ecological data from nearly all known species of the ancient fungal subphylum Saccharomycotina (1154 yeast strains from 1051 species), grown in 24 different environmental conditions, to examine niche breadth evolution. We found that large differences in the breadth of carbon utilization traits between yeasts stem from intrinsic differences in genes encoding specific metabolic pathways, but we found limited evidence for trade-offs. These comprehensive data argue that intrinsic factors shape niche breadth variation in microbes.

Genomic and ecological factors shaping specialism and generalism across an entire subphylum.

Organisms exhibit extensive variation in ecological niche breadth, from very narrow (specialists) to very broad (generalists). Paradigms proposed to explain this variation either invoke trade-offs between performance efficiency and breadth or underlying intrinsic or extrinsic factors. We assembled genomic (1,154 yeast strains from 1,049 species), metabolic (quantitative measures of growth of 843 species in 24 conditions), and ecological (environmental ontology of 1,088 species) data from nearly all known species of the ancient fungal subphylum Saccharomycotina to examine niche breadth evolution. We found large interspecific differences in carbon breadth stem from intrinsic differences in genes encoding specific metabolic pathways but no evidence of trade-offs and a limited role of extrinsic ecological factors. These comprehensive data argue that intrinsic factors driving microbial niche breadth variation.

Macroevolutionary diversity of traits and genomes in the model yeast genus Saccharomyces.

Species is the fundamental unit to quantify biodiversity. In recent years, the model yeast Saccharomyces cerevisiae has seen an increased number of studies related to its geographical distribution, population structure, and phenotypic diversity. However, seven additional species from the same genus have been less thoroughly studied, which has limited our understanding of the macroevolutionary events leading to the diversification of this genus over the last 20 million years. Here, we show the geographies, hosts, substrates, and phylogenetic relationships for approximately 1,800 Saccharomyces strains, covering the complete genus with unprecedented breadth and depth. We generated and analyzed complete genome sequences of 163 strains and phenotyped 128 phylogenetically diverse strains. This dataset provides insights about genetic and phenotypic diversity within and between species and populations, quantifies reticulation and incomplete lineage sorting, and demonstrates how gene flow and selection have affected traits, such as galactose metabolism. These findings elevate the genus Saccharomyces as a model to understand biodiversity and evolution in microbial eukaryotes.

Delila-PY, a Pipeline for Utilizing the Delila Suite of Software to Identify Potential DNA Binding Motifs.

Predicting potential DNA binding motifs is a critical part of understanding gene expression across all domains of life. Here, we report the development of Delila-PY, an easy-to-use pipeline for utilizing the Delila suite of software to identify DNA binding motifs.

Repeated horizontal gene transfer of GALactose metabolism genes violates Dollo's law of irreversible loss.

Dollo's law posits that evolutionary losses are irreversible, thereby narrowing the potential paths of evolutionary change. While phenotypic reversals to ancestral states have been observed, little is known about their underlying genetic causes. The genomes of budding yeasts have been shaped by extensive reductive evolution, such as reduced genome sizes and the losses of metabolic capabilities. However, the extent and mechanisms of trait reacquisition after gene loss in yeasts have not been thoroughly studied. Here, through phylogenomic analyses, we reconstructed the evolutionary history of the yeast galactose utilization pathway and observed widespread and repeated losses of the ability to utilize galactose, which occurred concurrently with the losses of GALactose (GAL) utilization genes. Unexpectedly, we detected multiple galactose-utilizing lineages that were deeply embedded within clades that underwent ancient losses of galactose utilization. We show that at least two, and possibly three, lineages reacquired the GAL pathway via yeast-to-yeast horizontal gene transfer. Our results show how trait reacquisition can occur tens of millions of years after an initial loss via horizontal gene transfer from distant relatives. These findings demonstrate that the losses of complex traits and even whole pathways are not always evolutionary dead-ends, highlighting how reversals to ancestral states can occur.

Genome-scale phylogeny and contrasting modes of genome evolution in the fungal phylum Ascomycota.

Ascomycota, the largest and most well-studied phylum of fungi, contains three subphyla: Saccharomycotina (budding yeasts), Pezizomycotina (filamentous fungi), and Taphrinomycotina (fission yeasts). Despite its importance, we lack a comprehensive genome-scale phylogeny or understanding of the similarities and differences in the mode of genome evolution within this phylum. By examining 1107 genomes from Saccharomycotina (332), Pezizomycotina (761), and Taphrinomycotina (14) species, we inferred a robust genome-wide phylogeny that resolves several contentious relationships and estimated that the Ascomycota last common ancestor likely originated in the Ediacaran period. Comparisons of genomic properties revealed that Saccharomycotina and Pezizomycotina differ greatly in their genome properties and enabled inference of the direction of evolutionary change. The Saccharomycotina typically have smaller genomes, lower guanine-cytosine contents, lower numbers of genes, and higher rates of molecular sequence evolution compared with Pezizomycotina. These results provide a robust evolutionary framework for understanding the diversity and ecological lifestyles of the largest fungal phylum.

Ethanol exposure increases mutation rate through error-prone polymerases.

Ethanol is a ubiquitous environmental stressor that is toxic to all lifeforms. Here, we use the model eukaryote Saccharomyces cerevisiae to show that exposure to sublethal ethanol concentrations causes DNA replication stress and an increased mutation rate. Specifically, we find that ethanol slows down replication and affects localization of Mrc1, a conserved protein that helps stabilize the replisome. In addition, ethanol exposure also results in the recruitment of error-prone DNA polymerases to the replication fork. Interestingly, preventing this recruitment through mutagenesis of the PCNA/Pol30 polymerase clamp or deleting specific error-prone polymerases abolishes the mutagenic effect of ethanol. Taken together, this suggests that the mutagenic effect depends on a complex mechanism, where dysfunctional replication forks lead to recruitment of error-prone polymerases. Apart from providing a general mechanistic framework for the mutagenic effect of ethanol, our findings may also provide a route to better understand and prevent ethanol-associated carcinogenesis in higher eukaryotes.

Multiple Reinventions of Mating-type Switching during Budding Yeast Evolution.

Cell type in budding yeasts is determined by the genotype at the mating-type (MAT) locus, but yeast species differ widely in their mating compatibility systems and life cycles. Among sexual yeasts, heterothallic species are those in which haploid strains fall into two distinct and stable mating types (MATa and MATα), whereas homothallic species are those that can switch mating types or that appear not to have distinct mating types [1, 2]. The evolutionary history of these mating compatibility systems is uncertain, particularly regarding the number and direction of transitions between homothallism and heterothallism, and regarding whether the process of mating-type switching had a single origin [3-5]. Here, we inferred the mating compatibility systems of 332 budding yeast species from their genome sequences. By reference to a robust phylogenomic tree [6], we detected evolutionary transitions between heterothallism and homothallism, and among different forms of homothallism. We find that mating-type switching has arisen independently at least 11 times during yeast evolution and that transitions from heterothallism to homothallism greatly outnumber transitions in the opposite direction (31 versus 3). Although the 3-locus MAT-HML-HMR mechanism of mating-type switching as seen in Saccharomyces cerevisiae had a single evolutionary origin in budding yeasts, simpler "flip/flop" mechanisms of switching evolved separately in at least 10 other groups of yeasts. These results point to the adaptive value of homothallism and mating-type switching to unicellular fungi.

Extensive loss of cell-cycle and DNA repair genes in an ancient lineage of bipolar budding yeasts.

Cell-cycle checkpoints and DNA repair processes protect organisms from potentially lethal mutational damage. Compared to other budding yeasts in the subphylum Saccharomycotina, we noticed that a lineage in the genus Hanseniaspora exhibited very high evolutionary rates, low Guanine-Cytosine (GC) content, small genome sizes, and lower gene numbers. To better understand Hanseniaspora evolution, we analyzed 25 genomes, including 11 newly sequenced, representing 18/21 known species in the genus. Our phylogenomic analyses identify two Hanseniaspora lineages, a faster-evolving lineage (FEL), which began diversifying approximately 87 million years ago (mya), and a slower-evolving lineage (SEL), which began diversifying approximately 54 mya. Remarkably, both lineages lost genes associated with the cell cycle and genome integrity, but these losses were greater in the FEL. E.g., all species lost the cell-cycle regulator WHIskey 5 (WHI5), and the FEL lost components of the spindle checkpoint pathway (e.g., Mitotic Arrest-Deficient 1 [MAD1], Mitotic Arrest-Deficient 2 [MAD2]) and DNA-damage-checkpoint pathway (e.g., Mitosis Entry Checkpoint 3 [MEC3], RADiation sensitive 9 [RAD9]). Similarly, both lineages lost genes involved in DNA repair pathways, including the DNA glycosylase gene 3-MethylAdenine DNA Glycosylase 1 (MAG1), which is part of the base-excision repair pathway, and the DNA photolyase gene PHotoreactivation Repair deficient 1 (PHR1), which is involved in pyrimidine dimer repair. Strikingly, the FEL lost 33 additional genes, including polymerases (i.e., POLymerase 4 [POL4] and POL32) and telomere-associated genes (e.g., Repressor/activator site binding protein-Interacting Factor 1 [RIF1], Replication Factor A 3 [RFA3], Cell Division Cycle 13 [CDC13], Pbp1p Binding Protein [PBP2]). Echoing these losses, molecular evolutionary analyses reveal that, compared to the SEL, the FEL stem lineage underwent a burst of accelerated evolution, which resulted in greater mutational loads, homopolymer instabilities, and higher fractions of mutations associated with the common endogenously damaged base, 8-oxoguanine. We conclude that Hanseniaspora is an ancient lineage that has diversified and thrived, despite lacking many otherwise highly conserved cell-cycle and genome integrity genes and pathways, and may represent a novel, to our knowledge, system for studying cellular life without them.

Eukaryotic Acquisition of a Bacterial Operon.

Operons are a hallmark of bacterial genomes, where they allow concerted expression of functionally related genes as single polycistronic transcripts. They are rare in eukaryotes, where each gene usually drives expression of its own independent messenger RNAs. Here, we report the horizontal operon transfer of a siderophore biosynthesis pathway from relatives of Escherichia coli into a group of budding yeast taxa. We further show that the co-linearly arranged secondary metabolism genes are expressed, exhibit eukaryotic transcriptional features, and enable the sequestration and uptake of iron. After transfer, several genetic changes occurred during subsequent evolution, including the gain of new transcription start sites that were sometimes within protein-coding sequences, acquisition of polyadenylation sites, structural rearrangements, and integration of eukaryotic genes into the cluster. We conclude that the genes were likely acquired as a unit, modified for eukaryotic gene expression, and maintained by selection to adapt to the highly competitive, iron-limited environment.

Tempo and Mode of Genome Evolution in the Budding Yeast Subphylum.

Budding yeasts (subphylum Saccharomycotina) are found in every biome and are as genetically diverse as plants or animals. To understand budding yeast evolution, we analyzed the genomes of 332 yeast species, including 220 newly sequenced ones, which represent nearly one-third of all known budding yeast diversity. Here, we establish a robust genus-level phylogeny comprising 12 major clades, infer the timescale of diversification from the Devonian period to the present, quantify horizontal gene transfer (HGT), and reconstruct the evolution of 45 metabolic traits and the metabolic toolkit of the budding yeast common ancestor (BYCA). We infer that BYCA was metabolically complex and chronicle the tempo and mode of genomic and phenotypic evolution across the subphylum, which is characterized by very low HGT levels and widespread losses of traits and the genes that control them. More generally, our results argue that reductive evolution is a major mode of evolutionary diversification.

Functional and evolutionary characterization of a secondary metabolite gene cluster in budding yeasts.

Secondary metabolites are key in how organisms from all domains of life interact with their environment and each other. The iron-binding molecule pulcherrimin was described a century ago, but the genes responsible for its production in budding yeasts have remained uncharacterized. Here, we used phylogenomic footprinting on 90 genomes across the budding yeast subphylum Saccharomycotina to identify the gene cluster associated with pulcherrimin production. Using targeted gene replacements in Kluyveromyces lactis, we characterized the four genes that make up the cluster, which likely encode two pulcherriminic acid biosynthesis enzymes, a pulcherrimin transporter, and a transcription factor involved in both biosynthesis and transport. The requirement of a functional putative transporter to utilize extracellular pulcherrimin-complexed iron demonstrates that pulcherriminic acid is a siderophore, a chelator that binds iron outside the cell for subsequent uptake. Surprisingly, we identified homologs of the putative transporter and transcription factor genes in multiple yeast genera that lacked the biosynthesis genes and could not make pulcherrimin, including the model yeast Saccharomyces cerevisiae We deleted these previously uncharacterized genes and showed they are also required for pulcherrimin utilization in S. cerevisiae, raising the possibility that other genes of unknown function are linked to secondary metabolism. Phylogenetic analyses of this gene cluster suggest that pulcherrimin biosynthesis and utilization were ancestral to budding yeasts, but the biosynthesis genes and, subsequently, the utilization genes, were lost in many lineages, mirroring other microbial public goods systems that lead to the rise of cheater organisms.

Evolutionary instability of CUG-Leu in the genetic code of budding yeasts.

The genetic code used in nuclear genes is almost universal, but here we report that it changed three times in parallel during the evolution of budding yeasts. All three changes were reassignments of the codon CUG, which is translated as serine (in 2 yeast clades), alanine (1 clade), or the 'universal' leucine (2 clades). The newly discovered Ser2 clade is in the final stages of a genetic code transition. Most species in this clade have genes for both a novel tRNASer(CAG) and an ancestral tRNALeu(CAG) to read CUG, but only tRNASer(CAG) is used in standard growth conditions. The coexistence of these alloacceptor tRNA genes indicates that the genetic code transition occurred via an ambiguous translation phase. We propose that the three parallel reassignments of CUG were not driven by natural selection in favor of their effects on the proteome, but by selection to eliminate the ancestral tRNALeu(CAG).

Repeated Cis-Regulatory Tuning of a Metabolic Bottleneck Gene during Evolution.

Repeated evolutionary events imply underlying genetic constraints that can make evolutionary mechanisms predictable. Morphological traits are thought to evolve frequently through cis-regulatory changes because these mechanisms bypass constraints in pleiotropic genes that are reused during development. In contrast, the constraints acting on metabolic traits during evolution are less well studied. Here we show how a metabolic bottleneck gene has repeatedly adopted similar cis-regulatory solutions during evolution, likely due to its pleiotropic role integrating flux from multiple metabolic pathways. Specifically, the genes encoding phosphoglucomutase activity (PGM1/PGM2), which connect GALactose catabolism to glycolysis, have gained and lost direct regulation by the transcription factor Gal4 several times during yeast evolution. Through targeted mutations of predicted Gal4-binding sites in yeast genomes, we show this galactose-mediated regulation of PGM1/2 supports vigorous growth on galactose in multiple yeast species, including Saccharomyces uvarum and Lachancea kluyveri. Furthermore, the addition of galactose-inducible PGM1 alone is sufficient to improve the growth on galactose of multiple species that lack this regulation, including Saccharomyces cerevisiae. The strong association between regulation of PGM1/2 by Gal4 even enables remarkably accurate predictions of galactose growth phenotypes between closely related species. This repeated mode of evolution suggests that this specific cis-regulatory connection is a common way that diverse yeasts can govern flux through the pathway, likely due to the constraints imposed by this pleiotropic bottleneck gene. Since metabolic pathways are highly interconnected, we argue that cis-regulatory evolution might be widespread at pleiotropic genes that control metabolic bottlenecks and intersections.

Evidence for loss and reacquisition of alcoholic fermentation in a fructophilic yeast lineage.

Fructophily is a rare trait that consists of the preference for fructose over other carbon sources. Here, we show that in a yeast lineage (the Wickerhamiella/Starmerella, W/S clade) comprised of fructophilic species thriving in the high-sugar floral niche, the acquisition of fructophily is concurrent with a wider remodeling of central carbon metabolism. Coupling comparative genomics with biochemical and genetic approaches, we gathered ample evidence for the loss of alcoholic fermentation in an ancestor of the W/S clade and subsequent reinstatement through either horizontal acquisition of homologous bacterial genes or modification of a pre-existing yeast gene. An enzyme required for sucrose assimilation was also acquired from bacteria, suggesting that the genetic novelties identified in the W/S clade may be related to adaptation to the high-sugar environment. This work shows how even central carbon metabolism can be remodeled by a surge of HGT events.

Variable repeats in the eukaryotic polyubiquitin gene ubi4 modulate proteostasis and stress survival.

Ubiquitin conjugation signals for selective protein degradation by the proteasome. In eukaryotes, ubiquitin is encoded both as a monomeric ubiquitin unit fused to a ribosomal gene and as multiple ubiquitin units in tandem. The polyubiquitin gene is a unique, highly conserved open reading frame composed solely of tandem repeats, yet it is still unclear why cells utilize this unusual gene structure. Using the Saccharomyces cerevisiae UBI4 gene, we show that this multi-unit structure allows cells to rapidly produce large amounts of ubiquitin needed to respond to sudden stress. The number of ubiquitin units encoded by UBI4 influences cellular survival and the rate of ubiquitin-proteasome system (UPS)-mediated proteolysis following heat stress. Interestingly, the optimal number of repeats varies under different types of stress indicating that natural variation in repeat numbers may optimize the chance for survival. Our results demonstrate how a variable polycistronic transcript provides an evolutionary alternative for gene copy number variation.Eukaryotic cells rely on the ubiquitin-proteasome system for selective degradation of proteins, a process vital to organismal fitness. Here the authors show that the number of repeats in the polyubiquitin gene is evolutionarily unstable within and between yeast species, and that this variability may tune the cell's capacity to respond to sudden environmental perturbations.

Spathaspora boniae sp. nov., a D-xylose-fermenting species in the Candida albicans/Lodderomyces clade.

Two yeast isolates producing asci-containing elongate ascospores with curved ends typical of the genus Spathaspora were isolated from rotting wood samples collected in an Atlantic rainforest ecosystem in Brazil. Phylogenetic analysis of the LSU rRNA gene D1/D2 domain sequences demonstrated that the strains represent a new species and placed it next to Candida blackwellae, in a clade that also contains Candida albicans and Candida dubliniensis. Other sequences of the ribosomal gene cluster supported same placementin the same clade, and a phylogenomic analysis placed this new species in an early emerging position relative to the larger C. albicans/Lodderomyces clade. One interpretation is that the genus Spathaspora is, in fact, paraphyletic. In conformity with this view, we propose the novel species Spathaspora boniae sp. nov. to accommodate the isolates. The type strain of Spathaspora boniae sp. nov. is UFMG-CM-Y306T (=CBS 13262T). The MycoBank number is MB 821297. A detailed analysis of xylose metabolism was conducted for the new species.

Genome sequence and physiological analysis of Yamadazyma laniorum f.a. sp. nov. and a reevaluation of the apocryphal xylose fermentation of its sister species, Candida tenuis.

Xylose fermentation is a rare trait that is immensely important to the cellulosic biofuel industry, and Candida tenuis is one of the few yeasts that has been reported with this trait. Here we report the isolation of two strains representing a candidate sister species to C. tenuis. Integrated analysis of genome sequence and physiology suggested the genetic basis of a number of traits, including variation between the novel species and C. tenuis in lactose metabolism due to the loss of genes encoding lactose permease and ß-galactosidase in the former. Surprisingly, physiological characterization revealed that neither the type strain of C. tenuis nor this novel species fermented xylose in traditional assays. We reexamined three xylose-fermenting strains previously identified as C. tenuis and found that these strains belong to the genus Scheffersomyces and are not C. tenuis. We propose Yamadazyma laniorum f.a. sp. nov. to accommodate our new strains and designate its type strain as yHMH7 (=CBS 14780 = NRRL Y-63967T). Furthermore, we propose the transfer of Candida tenuis to the genus Yamadazyma as Yamadazyma tenuis comb. nov. This approach provides a roadmap for how integrated genome sequence and physiological analysis can yield insight into the mechanisms that generate yeast biodiversity.

In Silico Whole Genome Sequencer and Analyzer (iWGS): a Computational Pipeline to Guide the Design and Analysis of de novo Genome Sequencing Studies.

The availability of genomes across the tree of life is highly biased toward vertebrates, pathogens, human disease models, and organisms with relatively small and simple genomes. Recent progress in genomics has enabled the de novo decoding of the genome of virtually any organism, greatly expanding its potential for understanding the biology and evolution of the full spectrum of biodiversity. The increasing diversity of sequencing technologies, assays, and de novo assembly algorithms have augmented the complexity of de novo genome sequencing projects in nonmodel organisms. To reduce the costs and challenges in de novo genome sequencing projects and streamline their experimental design and analysis, we developed iWGS ( in silicoWhole Genome Sequencer and Analyzer), an automated pipeline for guiding the choice of appropriate sequencing strategy and assembly protocols. iWGS seamlessly integrates the four key steps of a de novo genome sequencing project: data generation (through simulation), data quality control, de novo assembly, and assembly evaluation and validation. The last three steps can also be applied to the analysis of real data. iWGS is designed to enable the user to have great flexibility in testing the range of experimental designs available for genome sequencing projects, and supports all major sequencing technologies and popular assembly tools. Three case studies illustrate how iWGS can guide the design of de novo genome sequencing projects, and evaluate the performance of a wide variety of user-specified sequencing strategies and assembly protocols on genomes of differing architectures. iWGS, along with a detailed documentation, is freely available at https://github.com/zhouxiaofan1983/iWGS.

Reconstructing the Backbone of the Saccharomycotina Yeast Phylogeny Using Genome-Scale Data.

Understanding the phylogenetic relationships among the yeasts of the subphylum Saccharomycotina is a prerequisite for understanding the evolution of their metabolisms and ecological lifestyles. In the last two decades, the use of rDNA and multilocus data sets has greatly advanced our understanding of the yeast phylogeny, but many deep relationships remain unsupported. In contrast, phylogenomic analyses have involved relatively few taxa and lineages that were often selected with limited considerations for covering the breadth of yeast biodiversity. Here we used genome sequence data from 86 publicly available yeast genomes representing nine of the 11 known major lineages and 10 nonyeast fungal outgroups to generate a 1233-gene, 96-taxon data matrix. Species phylogenies reconstructed using two different methods (concatenation and coalescence) and two data matrices (amino acids or the first two codon positions) yielded identical and highly supported relationships between the nine major lineages. Aside from the lineage comprised by the family Pichiaceae, all other lineages were monophyletic. Most interrelationships among yeast species were robust across the two methods and data matrices. However, eight of the 93 internodes conflicted between analyses or data sets, including the placements of: the clade defined by species that have reassigned the CUG codon to encode serine, instead of leucine; the clade defined by a whole genome duplication; and the species Ascoidea rubescens These phylogenomic analyses provide a robust roadmap for future comparative work across the yeast subphylum in the disciplines of taxonomy, molecular genetics, evolutionary biology, ecology, and biotechnology. To further this end, we have also provided a BLAST server to query the 86 Saccharomycotina genomes, which can be found at http://y1000plus.org/blast.