Jacek Kominek

Genomics and the making of yeast biodiversity

Yeasts are unicellular fungi that do not form fruiting bodies. Although the yeast lifestyle has evolved multiple times, most known species belong to the subphylum Saccharomycotina (syn. Hemiascomycota, hereafter yeasts). This diverse group includes the premier eukaryotic model system, Saccharomyces cerevisiae; the common human commensal and opportunistic pathogen, Candida albicans; and over 1000 other known species (with more continuing to be discovered). Yeasts are found in every biome and continent and are more genetically diverse than angiosperms or chordates. Ease of culture, simple life cycles, and small genomes (∼10-20Mbp) have made yeasts exceptional models for molecular genetics, biotechnology, and evolutionary genomics. Here we discuss recent developments in understanding the genomic underpinnings of the making of yeast biodiversity, comparing and contrasting natural and human-associated evolutionary processes. Only a tiny fraction of yeast biodiversity and metabolic capabilities has been tapped by industry and science. Expanding the taxonomic breadth of deep genomic investigations will further illuminate how genome function evolves to encode their diverse metabolisms and ecologies.

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.

Genomic analysis and D-xylose fermentation of three novel Spathaspora species: Spathaspora girioi sp. nov., Spathaspora hagerdaliae f. a., sp. nov. and Spathaspora gorwiae f. a., sp. nov.

Three novel D-xylose-fermenting yeast species of Spathaspora clade were recovered from rotting wood in regions of the Atlantic Rainforest ecosystem in Brazil. Differentiation of new species was based on analyses of the gene encoding the D1/D2 sequences of large subunit of rRNA and on 642 conserved, single-copy, orthologous genes from genome sequence assemblies from the newly described species and 15 closely-related Debaryomycetaceae/Metschnikowiaceae species. Spathaspora girioi sp. nov. produced unconjugated asci with a single elongated ascospore with curved ends; ascospore formation was not observed for the other two species. The three novel species ferment D-xylose with different efficiencies. Spathaspora hagerdaliae sp. nov. and Sp. girioi sp. nov. showed xylose reductase (XR) activity strictly dependent on NADPH, whereas Sp. gorwiae sp. nov. had XR activity that used both NADH and NADPH as co-factors. The genes that encode enzymes involved in D-xylose metabolism (XR, xylitol dehydrogenase and xylulokinase) were also identified for these novel species. The type strains are Sp. girioi sp. nov. UFMG-CM-Y302(T) (=CBS 13476), Sp. hagerdaliae f.a., sp. nov. UFMG-CM-Y303(T) (=CBS 13475) and Sp. gorwiae f.a., sp. nov. UFMG-CM-Y312(T) (=CBS 13472).

in silico Whole Genome Sequencer & 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 silico Whole 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.

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.

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.

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

&NA; 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 &bgr;‐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. Graphical Abstract Figure. Here we report the discovery, genome sequence, physiological characterization and taxonomic description of a new species of yeast that is sister to a species previously believed to ferment xylose. Figure. No caption abailable.

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).The genetic code for amino acids is nearly universal, and among eukaryotic nuclear genomes the only known reassignments are of codon CUG in yeasts. Here, the authors identify a third independent CUG transition in budding yeasts that is still ongoing with alternative tRNAs present in the genome.

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

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

Eukaryotic acquisition of a bacterial operon

Operons are a hallmark of bacterial genomes, where they allow concerted expression of multiple 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 catecholate-class siderophore biosynthesis pathway from Enterobacteriaceae into a group of closely related yeast taxa. We further show that the co-linearly arranged secondary metabolism genes are actively expressed, exhibit mainly eukaryotic transcriptional features, and enable the sequestration and uptake of iron. After transfer to the eukaryotic host, several genetic changes occurred, including the acquisition of polyadenylation sites, structural rearrangements, integration of eukaryotic genes, and secondary loss in some lineages. We conclude that the operon genes were likely captured in the shared insect gut habitat, modified for eukaryotic gene expression, and maintained by selection to adapt to the highly-competitive, iron-limited environment.