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.

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.

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.

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.

Adaptation to High Ethanol Reveals Complex Evolutionary Pathways.

Tolerance to high levels of ethanol is an ecologically and industrially relevant phenotype of microbes, but the molecular mechanisms underlying this complex trait remain largely unknown. Here, we use long-term experimental evolution of isogenic yeast populations of different initial ploidy to study adaptation to increasing levels of ethanol. Whole-genome sequencing of more than 30 evolved populations and over 100 adapted clones isolated throughout this two-year evolution experiment revealed how a complex interplay of de novo single nucleotide mutations, copy number variation, ploidy changes, mutator phenotypes, and clonal interference led to a significant increase in ethanol tolerance. Although the specific mutations differ between different evolved lineages, application of a novel computational pipeline, PheNetic, revealed that many mutations target functional modules involved in stress response, cell cycle regulation, DNA repair and respiration. Measuring the fitness effects of selected mutations introduced in non-evolved ethanol-sensitive cells revealed several adaptive mutations that had previously not been implicated in ethanol tolerance, including mutations in PRT1, VPS70 and MEX67. Interestingly, variation in VPS70 was recently identified as a QTL for ethanol tolerance in an industrial bio-ethanol strain. Taken together, our results show how, in contrast to adaptation to some other stresses, adaptation to a continuous complex and severe stress involves interplay of different evolutionary mechanisms. In addition, our study reveals functional modules involved in ethanol resistance and identifies several mutations that could help to improve the ethanol tolerance of industrial yeasts.

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.

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.

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

Overlapping binding sites of the frataxin homologue assembly factor and the heat shock protein 70 transfer factor on the Isu iron-sulfur cluster scaffold protein.

In mitochondria FeS clusters, prosthetic groups critical for the activity of many proteins, are first assembled on Isu, a 14-kDa scaffold protein, and then transferred to recipient apoproteins. The assembly process involves interaction of Isu with both Nfs1, the cysteine desulfurase serving as a sulfur donor, and the yeast frataxin homolog (Yfh1) serving as a regulator of desulfurase activity and/or iron donor. Here, based on the results of biochemical experiments with purified wild-type and variant proteins, we report that interaction of Yfh1 with both Nfs1 and Isu are required for formation of a stable tripartite assembly complex. Disruption of either Yfh1-Isu or Nfs1-Isu interactions destabilizes the complex. Cluster transfer to recipient apoprotein is known to require the interaction of Isu with the J-protein/Hsp70 molecular chaperone pair, Jac1 and Ssq1. Here we show that the Yfh1 interaction with Isu involves the PVK sequence motif, which is also the site key for the interaction of Isu with Hsp70 Ssq1. Coupled with our previous observation that Nfs1 and Jac1 binding to Isu is mutually exclusive due to partially overlapping binding sites, we propose that such mutual exclusivity of cluster assembly factor (Nfs1/Yfh1) and cluster transfer factor (Jac1/Ssq1) binding to Isu has functional consequences for the transition from the assembly process to the transfer process, and thus regulation of the biogenesis of FeS cluster proteins.

Binding of the chaperone Jac1 protein and cysteine desulfurase Nfs1 to the iron-sulfur cluster scaffold Isu protein is mutually exclusive.

Biogenesis of mitochondrial iron-sulfur (Fe/S) cluster proteins requires the interaction of multiple proteins with the highly conserved 14-kDa scaffold protein Isu, on which clusters are built prior to their transfer to recipient proteins. For example, the assembly process requires the cysteine desulfurase Nfs1, which serves as the sulfur donor for cluster assembly. The transfer process requires Jac1, a J-protein Hsp70 cochaperone. We recently identified three residues on the surface of Jac1 that form a hydrophobic patch critical for interaction with Isu. The results of molecular modeling of the Isu1-Jac1 interaction, which was guided by these experimental data and structural/biophysical information available for bacterial homologs, predicted the importance of three hydrophobic residues forming a patch on the surface of Isu1 for interaction with Jac1. Using Isu variants having alterations in residues that form the hydrophobic patch on the surface of Isu, this prediction was experimentally validated by in vitro binding assays. In addition, Nfs1 was found to require the same hydrophobic residues of Isu for binding, as does Jac1, suggesting that Jac1 and Nfs1 binding is mutually exclusive. In support of this conclusion, Jac1 and Nfs1 compete for binding to Isu. Evolutionary analysis revealed that residues involved in these interactions are conserved and that they are critical residues for the biogenesis of Fe/S cluster protein in vivo. We propose that competition between Jac1 and Nfs1 for Isu binding plays an important role in transitioning the Fe/S cluster biogenesis machinery from the cluster assembly step to the Hsp70-mediated transfer of the Fe/S cluster to recipient proteins.

Sequential duplications of an ancient member of the DnaJ-family expanded the functional chaperone network in the eukaryotic cytosol.

Across eukaryotes, Hsp70-based chaperone machineries display an underlying unity in their sequence, structure, and biochemical mechanism of action, while working in a myriad of cellular processes. In good part, this extraordinary functional versatility is derived from the ability of a single Hsp70 to interact with an array of J-protein cochaperones to form a functional chaperone network. Among J-proteins, the DnaJ-type is the most prevalent, being present in all three kingdoms and in several different compartments of eukaryotic cells. However, because these ancient DnaJ-type proteins diverged at the base of the eukaryotic phylogeny, little is understood about the evolutionary basis of their diversification and thus the functional expansion of the chaperone network. Here, we report results of evolutionary and experimental analyses of two more recent members of the cytosolic DnaJ family of Saccharomyces cerevisiae, Xdj1 and Apj1, which emerged by sequential duplications of the ancient YDJ1 in Ascomycota. Sequence comparison and molecular modeling revealed that both Xdj1 and Apj1 maintained a domain organization similar to that of multifunctional Ydj1. However, despite these similarities, both Xdj1 and Apj1 evolved highly specialized functions. Xdj1 plays a unique role in the translocation of proteins from the cytosol into mitochondria. Apj1's specialized role is related to degradation of sumolyated proteins. Together these data provide the first clear example of cochaperone duplicates that evolved specialized functions, allowing expansion of the chaperone functional network, while maintaining the overall structural organization of their parental gene.

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.

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.

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.

Co-evolution-driven switch of J-protein specificity towards an Hsp70 partner.

Molecular mechanisms by which protein-protein interactions are preserved or lost after gene duplication are not understood. Taking advantage of the well-studied yeast mtHsp70:J-protein molecular chaperone system, we considered whether changes in partner proteins accompanied specialization of gene duplicates. Here, we report that existence of the Hsp70 Ssq1, which arose by duplication of the gene encoding multifunction mtHsp70 and specializes in iron-sulphur cluster biogenesis, correlates with functional and structural changes in the J domain of its J-protein partner Jac1. All species encoding this shorter alternative version of the J domain share a common ancestry, suggesting that all short JAC1 proteins arose from a single deletion event. Construction of a variant that extended the length of the J domain of a 'short' Jac1 enhanced its ability to partner with multifunctional Hsp70. Our data provide a causal link between changes in the J protein partner and specialization of duplicate Hsp70.

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.