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.

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.

Wickerhamiella verensis f.a. sp. nov., a novel yeast species isolated from subsoil groundwater contaminated with hydrocarbons and from a human infection.

Yeast strains belonging to a novel anamorphic yeast species were isolated from subsoil groundwater contaminated with hydrocarbons in a metal working factory located in northern Spain, and from a human infection in the USA. Comparison of ITS sequences between the isolates revealed 0.2 % divergence between the Spanish isolates and 0.46 % divergence between those and the USA isolate. Phylogenetic analysis based on the D1/D2 domains of the LSU rRNA gene showed that these isolates belong to the Wickerhamiella clade with W. sorbophila and W. infanticola as their closest relatives. Sequence divergence between the new isolates and W. sorbophila and W. infanticola was 1.97 and 1.79 %, respectively. The isolates in the novel species are not fermentative and pseudohyphae were not produced. Sexual reproduction was not observed for individual isolates or in mixtures of isolates. Conjugation between the isolates in the novel species and close relatives W. sorbophila and W. infanticola was not observed. These data support the proposal of Wickerhamiella verensis as a novel species, with CECT 12028T as the holotype.

Screening for Oily Yeasts Able to Convert Hydrolysates from Biomass to Biofuels While Maintaining Industrial Process Relevance.

Research has recently intensified to discover new oleaginous yeast strains able to function quickly and efficiently in low-cost lignocellulosic hydrolysates to produce high-quality lipids for use in biodiesel and chemicals. Detailed techniques are given here for ranking candidate yeast strains based on conversion of hydrolysate sugars to lipids and then optimizing cultivation conditions for best performers in a 96-well aerobic microcultivation format. A full battery of assays applicable to high throughput of small-volume samples are described for efficiently evaluating cell biomass production, lipid accumulation, fatty acid composition, and sugar utilization. Original data is additionally presented on the validation of the microtechnique for GC analysis of lipid composition in yeast since this application involved modification of a previously published assay for microalgae.

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.

Three new species of Tremellomycetes isolated from maize and northern wild rice.

The present work studied novel basidiomycetous yeasts from maize and northern wild rice plants. From comparisons of ribosomal internal transcribed spacer region (ITS) and large subunit (LSU) (D1 and D2 domains), and subsequent phylogenetic analyses, the following species were resolved and described: Papiliotrema zeae Yurkov & Kurtzman sp. nov. (ex-type cultures DSM 104035, NRRL Y-63980, MB 827356, GenBank MH718306), Solicoccozyma zizaniae Yurkov & Kurtzman sp. nov. (ex-type cultures DSM 104031, NRRL Y-7649, MB 827354, GenBank MH718302) and Vishniacozyma kurtzmanii Yurkov sp. nov. (ex-type cultures DSM 104032, NRRL Y-63981, MB 827355, GenBank MH718303). A search among environmental sequences showed that all three yeasts were previously detected, but not reliably assigned to a genus or clade. Papiliotrema zeae from maize and S. zizaniae from northern wild rice were previously found in agricultural soils under maize and rice, respectively.

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).

Four new species of Metschnikowia and the transfer of seven Candida species to Metschnikowia and Clavispora as new combinations.

From comparisons of ITS1-5.8S-ITS2 and gene sequences for nuclear D1/D2 LSU rRNA, nuclear SSU (18S) rRNA, translation elongation factor 1-α (EF1-α) and RNA polymerase II subunit 2 (RPB2), the following four new ascosporogenous yeast species were resolved and are described as Metschnikowia anglica (NRRL Y-7298T [type strain], CBS 15342, MycoBank MB 823167), Metschnikowia leonuri (NRRL Y-6546T, CBS 15341, MB 823166), Metschnikowia peoriensis (NRRL Y-5942T, CBS 15345, MB 823164) and Metschnikowia rubicola (NRRL Y-6064T, CBS 15344, MB 823165). The following six species of Candida are members of the Metschnikowia clade and are proposed for transfer to Metschnikowia as new combinations: Candida chrysomelidarum (NRRL Y-27749T, CBS 9904, MB 823223), Candida gelsemii (NRRL Y-48212T, CBS 10509, MB 823192), Candida kofuensis (NRRL Y-27226T, CBS 8058, MB 823195), Candida picachoensis (NRRL Y-27607T, CBS 9804, MB 823197), Candida pimensis (NRRL Y-27619T, CBS 9805, MB 823205) and Candida rancensis (NRRL Y-48702T, CBS 8174, MB 823224). Candida fructus (NRRL Y-17072T, CBS 6380, MB 823206) is transferred to Clavispora as a new combination, and Candida musae is shown to be a synonym of C. fructus. Apparent multiple alleles for ITS, D1/D2, EF1-α and RPB2 were detected in strains of some species.

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.

Factors driving metabolic diversity in the budding yeast subphylum.

BACKGROUND: Associations between traits are prevalent in nature, occurring across a diverse range of taxa and traits. Individual traits may co-evolve with one other, and these correlations can be driven by factors intrinsic or extrinsic to an organism. However, few studies, especially in microbes, have simultaneously investigated both across a broad taxonomic range. Here we quantify pairwise associations among 48 traits across 784 diverse yeast species of the ancient budding yeast subphylum Saccharomycotina, assessing the effects of phylogenetic history, genetics, and ecology. RESULTS: We find extensive negative (traits that tend to not occur together) and positive (traits that tend to co-occur) pairwise associations among traits, as well as between traits and environments. These associations can largely be explained by the biological properties of the traits, such as overlapping biochemical pathways. The isolation environments of the yeasts explain a minor but significant component of the variance, while phylogeny (the retention of ancestral traits in descendant species) plays an even more limited role. Positive correlations are pervasive among carbon utilization traits and track with chemical structures (e.g., glucosides and sugar alcohols) and metabolic pathways, suggesting a molecular basis for the presence of suites of traits. In several cases, characterized genes from model organisms suggest that enzyme promiscuity and overlapping biochemical pathways are likely mechanisms to explain these macroevolutionary trends. Interestingly, fermentation traits are negatively correlated with the utilization of pentose sugars, which are major components of the plant biomass degraded by fungi and present major bottlenecks to the production of cellulosic biofuels. Finally, we show that mammalian pathogenic and commensal yeasts have a suite of traits that includes growth at high temperature and, surprisingly, the utilization of a narrowed panel of carbon sources. CONCLUSIONS: These results demonstrate how both intrinsic physiological factors and extrinsic ecological factors drive the distribution of traits present in diverse organisms across macroevolutionary timescales.

Importance of Resolving Fungal Nomenclature: the Case of Multiple Pathogenic Species in the Cryptococcus Genus.

Cryptococcosis is a major fungal disease caused by members of the Cryptococcus gattii and Cryptococcus neoformans species complexes. After more than 15 years of molecular genetic and phenotypic studies and much debate, a proposal for a taxonomic revision was made. The two varieties within C. neoformans were raised to species level, and the same was done for five genotypes within C. gattii. In a recent perspective (K. J. Kwon-Chung et al., mSphere 2:e00357-16, 2017, https://doi.org/10.1128/mSphere.00357-16), it was argued that this taxonomic proposal was premature and without consensus in the community. Although the authors of the perspective recognized the existence of genetic diversity, they preferred the use of the informal nomenclature "C. neoformans species complex" and "C. gattii species complex." Here we highlight the advantage of recognizing these seven species, as ignoring these species will impede deciphering further biologically and clinically relevant differences between them, which may in turn delay future clinical advances.

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.

A survey of yeast from the Yarrowia clade for lipid production in dilute acid pretreated lignocellulosic biomass hydrolysate.

Yarrowia lipolytica is an oleaginous yeast species that has attracted attention as a model organism for synthesis of single cell oil. Among over 50 isolates of Y. lipolytica identified, only a few of the strains have been studied extensively. Furthermore, 12 other yeast species were recently assigned to the Yarrowia clade, and most are not well characterized in terms of cell growth and lipid accumulation, especially in industrially relevant conditions. In the present study, we investigated biomass and lipid production by 57 yeast isolates, representing all 13 species in the Yarrowia clade, on a non-detoxified dilute acid-pretreated switchgrass hydrolysate under highly aerobic conditions. The objective was to compare yeast physiology during growth in an abundant, low-cost biomass feedstock and to expand diversity of genetically tractable, oleaginous yeasts available for lipid research. Screening of 45 Y. lipolytica isolates demonstrated considerable variation within the species in terms of lipid accumulation (min = 0.1 g/L; max = 5.1 g/L; mean = 2.3 g/L); three strains (NRRL YB-420, YB-419, and YB-392) were especially promising for cellulosic biomass conversion with average improvements of 43, 57, and 64%, respectively, in final lipid titer as compared to control strain W29. Subsequently, evaluation of strains from 13 distinct species in the Yarrowia clade identified Candida phangngensis PT1-17 as the top lipid producer with a maximum titer of 9.8 g/L lipid, which was over twofold higher than the second-best species in the clade (Candida hollandica NRRL Y-48254). A small set of the most promising strains from the screenings was further characterized to evaluate inhibitor tolerance, lipid production kinetics, and fatty acid distribution. We expect that the results of this study will pave the way for new biotechnological applications involving previously overlooked and under-characterized strains within the Yarrowia clade.

Techniques for the Evolution of Robust Pentose-fermenting Yeast for Bioconversion of Lignocellulose to Ethanol.

Lignocellulosic biomass is an abundant, renewable feedstock useful for production of fuel-grade ethanol and other bio-products. Pretreatment and enzyme saccharification processes release sugars that can be fermented by yeast. Traditional industrial yeasts do not ferment xylose (comprising up to 40% of plant sugars) and are not able to function in concentrated hydrolyzates. Concentrated hydrolyzates are needed to support economical ethanol recovery, but they are laden with toxic byproducts generated during pretreatment. While detoxification methods can render hydrolyzates fermentable, they are costly and generate waste disposal liabilities. Here, adaptive evolution and isolation techniques are described and demonstrated to yield derivatives of the native Scheffersomyces stipitis strain NRRL Y-7124 that are able to efficiently convert hydrolyzates to economically recoverable ethanol despite adverse culture conditions. Improved individuals are enriched in an evolving population using multiple selection pressures reliant on natural genetic diversity of the S. stipitis population and mutations induced by exposures to two diverse hydrolyzates, ethanol or UV radiation. Final evolution cultures are dilution plated to harvest predominant isolates, while intermediate populations, frozen in glycerol at various stages of evolution, are enriched on selective media using appropriate stress gradients to recover most promising isolates through dilution plating. Isolates are screened on various hydrolyzate types and ranked using a novel procedure involving dimensionless relative performance index (RPI) transformations of the xylose uptake rate and ethanol yield data. Using the RPI statistical parameter, an overall relative performance average is calculated to rank isolates based on multiple factors, including culture conditions (varying in nutrients and inhibitors) and kinetic characteristics. Through application of these techniques, derivatives of the parent strain had the following improved features in enzyme saccharified hydrolyzates at pH 5-6: reduced initial lag phase preceding growth, reduced diauxic lag during glucose-xylose transition, significantly enhanced fermentation rates, improved ethanol tolerance and accumulation to 40 g/L.

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.