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.

Gene loss and cis-regulatory novelty shaped core histone gene evolution in the apiculate yeast Hanseniaspora uvarum.

Core histone genes display a remarkable diversity of cis-regulatory mechanisms despite their protein sequence conservation. However, the dynamics and significance of this regulatory turnover are not well understood. Here, we describe the evolutionary history of core histone gene regulation across 400 million years in budding yeasts. We find that canonical mode of core histone regulation-mediated by the trans-regulator Spt10-is ancient, likely emerging between 320 and 380 million years ago and is fixed in the majority of extant species. Unexpectedly, we uncovered the emergence of a novel core histone regulatory mode in the Hanseniaspora genus, from its fast-evolving lineage, which coincided with the loss of 1 copy of its paralogous core histone genes. We show that the ancestral Spt10 histone regulatory mode was replaced, via cis-regulatory changes in the histone control regions, by a derived Mcm1 histone regulatory mode and that this rewiring event occurred with no changes to the trans-regulator, Mcm1, itself. Finally, we studied the growth dynamics of the cell cycle and histone synthesis in genetically modified Hanseniaspora uvarum. We find that H. uvarum divides rapidly, with most cells completing a cell cycle within 60 minutes. Interestingly, we observed that the regulatory coupling between histone and DNA synthesis was lost in H. uvarum. Our results demonstrate that core histone gene regulation was fixed anciently in budding yeasts, however it has greatly diverged in the Hanseniaspora fast-evolving lineage.

Humanization reveals pervasive incompatibility of yeast and human kinetochore components.

Kinetochores assemble on centromeres to drive chromosome segregation in eukaryotic cells. Humans and budding yeast share most of the structural subunits of the kinetochore, whereas protein sequences have diverged considerably. The conserved centromeric histone H3 variant, CenH3 (CENP-A in humans and Cse4 in budding yeast), marks the site for kinetochore assembly in most species. A previous effort to complement Cse4 in yeast with human CENP-A was unsuccessful; however, co-complementation with the human core nucleosome was not attempted. Previously, our lab successfully humanized the core nucleosome in yeast; however, this severely affected cellular growth. We hypothesized that yeast Cse4 is incompatible with humanized nucleosomes and that the kinetochore represented a limiting factor for efficient histone humanization. Thus, we argued that including the human CENP-A or a Cse4-CENP-A chimera might improve histone humanization and facilitate kinetochore function in humanized yeast. The opposite was true: CENP-A expression reduced histone humanization efficiency, was toxic to yeast, and disrupted cell cycle progression and kinetochore function in wild-type (WT) cells. Suppressors of CENP-A toxicity included gene deletions of subunits of 3 conserved chromatin remodeling complexes, highlighting their role in CenH3 chromatin positioning. Finally, we attempted to complement the subunits of the NDC80 kinetochore complex, individually and in combination, without success, in contrast to a previous study indicating complementation by the human NDC80/HEC1 gene. Our results suggest that limited protein sequence similarity between yeast and human components in this very complex structure leads to failure of complementation.

Manipulating the 3D organization of the largest synthetic yeast chromosome.

Whether synthetic genomes can power life has attracted broad interest in the synthetic biology field. Here, we report de novo synthesis of the largest eukaryotic chromosome thus far, synIV, a 1,454,621-bp yeast chromosome resulting from extensive genome streamlining and modification. We developed megachunk assembly combined with a hierarchical integration strategy, which significantly increased the accuracy and flexibility of synthetic chromosome construction. Besides the drastic sequence changes, we further manipulated the 3D structure of synIV to explore spatial gene regulation. Surprisingly, we found few gene expression changes, suggesting that positioning inside the yeast nucleoplasm plays a minor role in gene regulation. Lastly, we tethered synIV to the inner nuclear membrane via its hundreds of loxPsym sites and observed transcriptional repression of the entire chromosome, demonstrating chromosome-wide transcription manipulation without changing the DNA sequences. Our manipulation of the spatial structure of synIV sheds light on higher-order architectural design of the synthetic genomes.

Context-dependent neocentromere activity in synthetic yeast chromosome VIII.

Pioneering advances in genome engineering, and specifically in genome writing, have revolutionized the field of synthetic biology, propelling us toward the creation of synthetic genomes. The Sc2.0 project aims to build the first fully synthetic eukaryotic organism by assembling the genome of Saccharomyces cerevisiae. With the completion of synthetic chromosome VIII (synVIII) described here, this goal is within reach. In addition to writing the yeast genome, we sought to manipulate an essential functional element: the point centromere. By relocating the native centromere sequence to various positions along chromosome VIII, we discovered that the minimal 118-bp CEN8 sequence is insufficient for conferring chromosomal stability at ectopic locations. Expanding the transplanted sequence to include a small segment (∼500 bp) of the CDEIII-proximal pericentromere improved chromosome stability, demonstrating that minimal centromeres display context-dependent functionality.

Taxogenomic analysis of a novel yeast species isolated from soil, Pichia galeolata sp. nov.

A novel budding yeast species was isolated from a soil sample collected in the United States of America. Phylogenetic analyses of multiple loci and phylogenomic analyses conclusively placed the species within the genus Pichia. Strain yHMH446 falls within a clade that includes Pichia norvegensis, Pichia pseudocactophila, Candida inconspicua, and Pichia cactophila. Whole genome sequence data were analyzed for the presence of genes known to be important for carbon and nitrogen metabolism, and the phenotypic data from the novel species were compared to all Pichia species with publicly available genomes. Across the genus, including the novel species candidate, we found that the inability to use many carbon and nitrogen sources correlated with the absence of metabolic genes. Based on these results, Pichia galeolata sp. nov. is proposed to accommodate yHMH446T (=NRRL Y-64187 = CBS 16864). This study shows how integrated taxogenomic analysis can add mechanistic insight to species descriptions.

Engineered transcription-associated Cas9 targeting in eukaryotic cells.

DNA targeting Class 2 CRISPR-Cas effector nucleases, including the well-studied Cas9 proteins, evolved protospacer-adjacent motif (PAM) and guide RNA interactions that sequentially license their binding and cleavage activities at protospacer target sites. Both interactions are nucleic acid sequence specific but function constitutively; thus, they provide intrinsic spatial control over DNA targeting activities but naturally lack temporal control. Here we show that engineered Cas9 fusion proteins which bind to nascent RNAs near a protospacer can facilitate spatiotemporal coupling between transcription and DNA targeting at that protospacer: Transcription-associated Cas9 Targeting (TraCT). Engineered TraCT is enabled when suboptimal PAM interactions limit basal activity in vivo and when one or more nascent RNA substrates are still tethered to the actively transcribing target DNA in cis. We further show that this phenomenon can be exploited for selective editing at one of two identical targets in distinct gene loci, or, in diploid allelic loci that are differentially transcribed. Our work demonstrates that temporal control over Cas9's targeting activity at specific DNA sites may be engineered without modifying Cas9's core domains and guide RNA components or their expression levels. More broadly, it establishes RNA binding in cis as a mechanism that can conditionally stimulate CRISPR-Cas DNA targeting in 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.

Human macroH2A1 drives nucleosome dephasing and genome instability in histone-humanized yeast.

In addition to replicative histones, eukaryotic genomes encode a repertoire of non-replicative variant histones providing additional layers of structural and epigenetic regulation. Here, we systematically replaced individual replicative human histones with non-replicative human variant histones using a histone replacement system in yeast. Variants H2A.J, TsH2B, and H3.5 complemented for their respective replicative counterparts. However, macroH2A1 failed to complement and its expression was toxic in yeast, negatively interacting with native yeast histones and kinetochore genes. To isolate yeast with "macroH2A1 chromatin" we decoupled the effects of its macro and histone fold domains, which revealed that both domains sufficed to override native yeast nucleosome positioning. Furthermore, both modified constructs of macroH2A1 exhibited lower nucleosome occupancy that correlated with decreased short-range chromatin interactions (<20 Kb), disrupted centromeric clustering, and increased chromosome instability. While supporting viability, macroH2A1 dramatically alters chromatin organization in yeast, leading to genome instability and massive fitness defects.

Protein-protein interactions in the Mla lipid transport system probed by computational structure prediction and deep mutational scanning.

The outer membrane (OM) of Gram-negative bacteria is an asymmetric bilayer that protects the cell from external stressors, such as antibiotics. The Mla transport system is implicated in the Maintenance of OM Lipid Asymmetry by mediating retrograde phospholipid transport across the cell envelope. Mla uses a shuttle-like mechanism to move lipids between the MlaFEDB inner membrane complex and the MlaA-OmpF/C OM complex, via a periplasmic lipid-binding protein, MlaC. MlaC binds to MlaD and MlaA, but the underlying protein-protein interactions that facilitate lipid transfer are not well understood. Here, we take an unbiased deep mutational scanning approach to map the fitness landscape of MlaC from Escherichia coli, which provides insights into important functional sites. Combining this analysis with AlphaFold2 structure predictions and binding experiments, we map the MlaC-MlaA and MlaC-MlaD protein-protein interfaces. Our results suggest that the MlaD and MlaA binding surfaces on MlaC overlap to a large extent, leading to a model in which MlaC can only bind one of these proteins at a time. Low-resolution cryo-electron microscopy (cryo-EM) maps of MlaC bound to MlaFEDB suggest that at least two MlaC molecules can bind to MlaD at once, in a conformation consistent with AlphaFold2 predictions. These data lead us to a model for MlaC interaction with its binding partners and insights into lipid transfer steps that underlie phospholipid transport between the bacterial inner and OMs.

Two differentially stable rDNA loci coexist on the same chromosome and form a single nucleolus.

The nucleolus is the most prominent membraneless compartment within the nucleus-dedicated to the metabolism of ribosomal RNA. Nucleoli are composed of hundreds of ribosomal DNA (rDNA) repeated genes that form large chromosomal clusters, whose high recombination rates can cause nucleolar dysfunction and promote genome instability. Intriguingly, the evolving architecture of eukaryotic genomes appears to have favored two strategic rDNA locations-where a single locus per chromosome is situated either near the centromere (CEN) or the telomere. Here, we deployed an innovative genome engineering approach to cut and paste to an ectopic chromosomal location-the ~1.5 mega-base rDNA locus in a single step using CRISPR technology. This "megablock" rDNA engineering was performed in a fused-karyotype strain of Saccharomyces cerevisiae. The strategic repositioning of this locus within the megachromosome allowed experimentally mimicking and monitoring the outcome of an rDNA migratory event, in which twin rDNA loci coexist on the same chromosomal arm. We showed that the twin-rDNA yeast readily adapts, exhibiting wild-type growth and maintaining rRNA homeostasis, and that the twin loci form a single nucleolus throughout the cell cycle. Unexpectedly, the size of each rDNA array appears to depend on its position relative to the CEN, in that the locus that is CEN-distal undergoes size reduction at a higher frequency compared to the CEN-proximal counterpart. Finally, we provided molecular evidence supporting a mechanism called paralogous cis-rDNA interference, which potentially explains why placing two identical repeated arrays on the same chromosome may negatively affect their function and structural stability.

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

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

Substrate, temperature, and geographical patterns among nearly 2000 natural yeast isolates.

Yeasts have broad importance as industrially and clinically relevant microbes and as powerful models for fundamental research, but we are only beginning to understand the roles yeasts play in natural ecosystems. Yeast ecology is often more difficult to study compared to other, more abundant microbes, but growing collections of natural yeast isolates are beginning to shed light on fundamental ecological questions. Here, we used environmental sampling and isolation to assemble a dataset of 1962 isolates collected from throughout the contiguous United States of America (USA) and Alaska, which were then used to uncover geographic patterns, along with substrate and temperature associations among yeast taxa. We found some taxa, including the common yeasts Torulaspora delbrueckii and Saccharomyces paradoxus, to be repeatedly isolated from multiple sampled regions of the USA, and we classify these as broadly distributed cosmopolitan yeasts. A number of yeast taxon-substrate associations were identified, some of which were novel and some of which support previously reported associations. Further, we found a strong effect of isolation temperature on the phyla of yeasts recovered, as well as for many species. We speculate that substrate and isolation temperature associations reflect the ecological diversity of and niche partitioning by yeast taxa.

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.

Population genomics of the pathogenic yeast Candida tropicalis identifies hybrid isolates in environmental samples.

Candida tropicalis is a human pathogen that primarily infects the immunocompromised. Whereas the genome of one isolate, C. tropicalis MYA-3404, was originally sequenced in 2009, there have been no large-scale, multi-isolate studies of the genetic and phenotypic diversity of this species. Here, we used whole genome sequencing and phenotyping to characterize 77 isolates of C. tropicalis from clinical and environmental sources from a variety of locations. We show that most C. tropicalis isolates are diploids with approximately 2-6 heterozygous variants per kilobase. The genomes are relatively stable, with few aneuploidies. However, we identified one highly homozygous isolate and six isolates of C. tropicalis with much higher heterozygosity levels ranging from 36-49 heterozygous variants per kilobase. Our analyses show that the heterozygous isolates represent two different hybrid lineages, where the hybrids share one parent (A) with most other C. tropicalis isolates, but the second parent (B or C) differs by at least 4% at the genome level. Four of the sequenced isolates descend from an AB hybridization, and two from an AC hybridization. The hybrids are MTLa/α heterozygotes. Hybridization, or mating, between different parents is therefore common in the evolutionary history of C. tropicalis. The new hybrids were predominantly found in environmental niches, including from soil. Hybridization is therefore unlikely to be associated with virulence. In addition, we used genotype-phenotype correlation and CRISPR-Cas9 editing to identify a genome variant that results in the inability of one isolate to utilize certain branched-chain amino acids as a sole nitrogen source.

Microbial Diversity Associated with the Pollen Stores of Captive-Bred Bumble Bee Colonies.

The pollen stores of bumble bees host diverse microbiota that influence overall colony fitness. Yet, the taxonomic identity of these symbiotic microbes is relatively unknown. In this descriptive study, we characterized the microbial community of pollen provisions within captive-bred bumble bee hives obtained from two commercial suppliers located in North America. Findings from 16S rRNA and ITS gene-based analyses revealed that pollen provisions from the captive-bred hives shared several microbial taxa that have been previously detected among wild populations. While diverse microbes across phyla Firmicutes, Proteobacteria, Bacteroidetes, Actinobacteria, and Ascomycota were detected in all commercial hives, significant differences were detected at finer-scale taxonomic resolution based on the supplier source. The causative agent of chalkbrood disease in honey bees, Ascosphaera apis, was detected in all hives obtained from one supplier source, although none of the hives showed symptoms of infection. The shared core microbiota across both commercial supplier sources consisted of two ubiquitous bee-associated groups, Lactobacillus and Wickerhamiella/Starmerella clade yeasts that potentially contribute to the beneficial function of the microbiome of bumble bee pollen provisions.

Superloser: A Plasmid Shuffling Vector for Saccharomyces cerevisiae with Exceedingly Low Background.

Here we report a new plasmid shuffle vector for forcing budding yeast (Saccharomyces cerevisiae) to incorporate a new genetic pathway in place of a native pathway - even an essential one - while maintaining low false positive rates (less than 1 in 108 per cell). This plasmid, dubbed "Superloser," was designed with reduced sequence similarity to commonly used yeast plasmids (i.e., pRS400 series) to limit recombination, a process that in our experience leads to retention of the yeast gene(s) instead of the desired gene(s). In addition, Superloser utilizes two orthogonal copies of the counter-selectable marker URA3 to reduce spontaneous 5-fluoroorotic acid resistance. Finally, the CEN/ARS sequence is fused to the GAL1-10 promoter, which disrupts plasmid segregation in the presence of the sugar galactose, causing Superloser to rapidly be removed from a population of cells. We show one proof-of-concept shuffling experiment: swapping yeast's core histones out for their human counterparts. Superloser is especially useful for forcing yeast to use highly unfavorable genes, such as human histones, as it enables plating a large number of cells (1.4x109) on a single 10 cm petri dish while maintaining a very low background. Therefore, Superloser is a useful tool for yeast geneticists to effectively shuffle low viability genes and/or pathways in yeast that may arise in as few as 1 in 108 cells.

Pathogenic budding yeasts isolated outside of clinical settings.

Budding yeasts are distributed across a wide range of habitats, including as human commensals. However, under some conditions, these commensals can cause superficial, invasive, and even lethal infections. Despite their importance to human health, little is known about the ecology of these opportunistic pathogens, aside from their associations with mammals and clinical environments. During a survey of approximately 1000 non-clinical samples across the United States of America, we isolated 54 strains of budding yeast species considered opportunistic pathogens, including Candida albicans and Candida (Nakaseomyces) glabrata. We found that, as a group, pathogenic yeasts were positively associated with fruits and soil environments, whereas the species Pichia kudriavzevii (syn. Candida krusei syn. Issatchenkia orientalis) had a significant association with plants. Of the four species that cause 95% of candidiasis, we found a positive association with soil. These results suggest that pathogenic yeast ecology is more complex and diverse than is currently appreciated and raises the possibility that these additional environments could be a point of contact for human infections.

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.

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.