A genomic compendium of cultivated human gut fungi characterizes the gut mycobiome and its relevance to common diseases.

The gut fungal community represents an essential element of human health, yet its functional and metabolic potential remains insufficiently elucidated, largely due to the limited availability of reference genomes. To address this gap, we presented the cultivated gut fungi (CGF) catalog, encompassing 760 fungal genomes derived from the feces of healthy individuals. This catalog comprises 206 species spanning 48 families, including 69 species previously unidentified. We explored the functional and metabolic attributes of the CGF species and utilized this catalog to construct a phylogenetic representation of the gut mycobiome by analyzing over 11,000 fecal metagenomes from Chinese and non-Chinese populations. Moreover, we identified significant common disease-related variations in gut mycobiome composition and corroborated the associations between fungal signatures and inflammatory bowel disease (IBD) through animal experimentation. These resources and findings substantially enrich our understanding of the biological diversity and disease relevance of the human gut mycobiome.

Debugging and consolidating multiple synthetic chromosomes reveals combinatorial genetic interactions.

The Sc2.0 project is building a eukaryotic synthetic genome from scratch. A major milestone has been achieved with all individual Sc2.0 chromosomes assembled. Here, we describe the consolidation of multiple synthetic chromosomes using advanced endoreduplication intercrossing with tRNA expression cassettes to generate a strain with 6.5 synthetic chromosomes. The 3D chromosome organization and transcript isoform profiles were evaluated using Hi-C and long-read direct RNA sequencing. We developed CRISPR Directed Biallelic URA3-assisted Genome Scan, or "CRISPR D-BUGS," to map phenotypic variants caused by specific designer modifications, known as "bugs." We first fine-mapped a bug in synthetic chromosome II (synII) and then discovered a combinatorial interaction associated with synIII and synX, revealing an unexpected genetic interaction that links transcriptional regulation, inositol metabolism, and tRNASerCGA abundance. Finally, to expedite consolidation, we employed chromosome substitution to incorporate the largest chromosome (synIV), thereby consolidating >50% of the Sc2.0 genome in one strain.

Design, construction, and functional characterization of a tRNA neochromosome in yeast.

Here, we report the design, construction, and characterization of a tRNA neochromosome, a designer chromosome that functions as an additional, de novo counterpart to the native complement of Saccharomyces cerevisiae. Intending to address one of the central design principles of the Sc2.0 project, the ∼190-kb tRNA neochromosome houses all 275 relocated nuclear tRNA genes. To maximize stability, the design incorporates orthogonal genetic elements from non-S. cerevisiae yeast species. Furthermore, the presence of 283 rox recombination sites enables an orthogonal tRNA SCRaMbLE system. Following construction in yeast, we obtained evidence of a potent selective force, manifesting as a spontaneous doubling in cell ploidy. Furthermore, tRNA sequencing, transcriptomics, proteomics, nucleosome mapping, replication profiling, FISH, and Hi-C were undertaken to investigate questions of tRNA neochromosome behavior and function. Its construction demonstrates the remarkable tractability of the yeast model and opens up opportunities to directly test hypotheses surrounding these essential non-coding RNAs.

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.

Scalable, Continuous Evolution of Genes at Mutation Rates above Genomic Error Thresholds.

Directed evolution is a powerful approach for engineering biomolecules and understanding adaptation. However, experimental strategies for directed evolution are notoriously labor intensive and low throughput, limiting access to demanding functions, multiple functions in parallel, and the study of molecular evolution in replicate. We report OrthoRep, an orthogonal DNA polymerase-plasmid pair in yeast that stably mutates ∼100,000-fold faster than the host genome in vivo, exceeding the error threshold of genomic replication that causes single-generation extinction. User-defined genes in OrthoRep continuously and rapidly evolve through serial passaging, a highly straightforward and scalable process. Using OrthoRep, we evolved drug-resistant malarial dihydrofolate reductases (DHFRs) in 90 independent replicates. We uncovered a more complex fitness landscape than previously realized, including common adaptive trajectories constrained by epistasis, rare outcomes that avoid a frequent early adaptive mutation, and a suboptimal fitness peak that occasionally traps evolving populations. OrthoRep enables a new paradigm of routine, high-throughput evolution of biomolecular and cellular function.

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.

Development of a Comprehensive Genotype-to-Fitness Map of Adaptation-Driving Mutations in Yeast.

Adaptive evolution plays a large role in generating the phenotypic diversity observed in nature, yet current methods are impractical for characterizing the molecular basis and fitness effects of large numbers of individual adaptive mutations. Here, we used a DNA barcoding approach to generate the genotype-to-fitness map for adaptation-driving mutations from a Saccharomyces cerevisiae population experimentally evolved by serial transfer under limiting glucose. We isolated and measured the fitness of thousands of independent adaptive clones and sequenced the genomes of hundreds of clones. We found only two major classes of adaptive mutations: self-diploidization and mutations in the nutrient-responsive Ras/PKA and TOR/Sch9 pathways. Our large sample size and precision of measurement allowed us to determine that there are significant differences in fitness between mutations in different genes, between different paralogs, and even between different classes of mutations within the same gene.

Domestication and Divergence of Saccharomyces cerevisiae Beer Yeasts.

Whereas domestication of livestock, pets, and crops is well documented, it is still unclear to what extent microbes associated with the production of food have also undergone human selection and where the plethora of industrial strains originates from. Here, we present the genomes and phenomes of 157 industrial Saccharomyces cerevisiae yeasts. Our analyses reveal that today's industrial yeasts can be divided into five sublineages that are genetically and phenotypically separated from wild strains and originate from only a few ancestors through complex patterns of domestication and local divergence. Large-scale phenotyping and genome analysis further show strong industry-specific selection for stress tolerance, sugar utilization, and flavor production, while the sexual cycle and other phenotypes related to survival in nature show decay, particularly in beer yeasts. Together, these results shed light on the origins, evolutionary history, and phenotypic diversity of industrial yeasts and provide a resource for further selection of superior strains. PAPERCLIP.

Genomic Nucleosome Organization Reconstituted with Pure Proteins.

Chromatin remodelers regulate genes by organizing nucleosomes around promoters, but their individual contributions are obfuscated by the complex in vivo milieu of factor redundancy and indirect effects. Genome-wide reconstitution of promoter nucleosome organization with purified proteins resolves this problem and is therefore a critical goal. Here, we reconstitute four stages of nucleosome architecture using purified components: yeast genomic DNA, histones, sequence-specific Abf1/Reb1, and remodelers RSC, ISW2, INO80, and ISW1a. We identify direct, specific, and sufficient contributions that in vivo observations validate. First, RSC clears promoters by translating poly(dA:dT) into directional nucleosome removal. Second, partial redundancy is recapitulated where INO80 alone, or ISW2 at Abf1/Reb1sites, positions +1 nucleosomes. Third, INO80 and ISW2 each align downstream nucleosomal arrays. Fourth, ISW1a tightens the spacing to canonical repeat lengths. Such a minimal set of rules and proteins establishes core mechanisms by which promoter chromatin architecture arises through a blend of redundancy and specialization.

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.

Global analysis of mRNA isoform half-lives reveals stabilizing and destabilizing elements in yeast.

We measured half-lives of 21,248 mRNA 3' isoforms in yeast by rapidly depleting RNA polymerase II from the nucleus and performing direct RNA sequencing throughout the decay process. Interestingly, half-lives of mRNA isoforms from the same gene, including nearly identical isoforms, often vary widely. Based on clusters of isoforms with different half-lives, we identify hundreds of sequences conferring stabilization or destabilization upon mRNAs terminating downstream. One class of stabilizing element is a polyU sequence that can interact with poly(A) tails, inhibit the association of poly(A)-binding protein, and confer increased stability upon introduction into ectopic transcripts. More generally, destabilizing and stabilizing elements are linked to the propensity of the poly(A) tail to engage in double-stranded structures. Isoforms engineered to fold into 3' stem-loop structures not involving the poly(A) tail exhibit even longer half-lives. We suggest that double-stranded structures at 3' ends are a major determinant of mRNA stability.

Subnucleosomal structures and nucleosome asymmetry across a genome.

Genes are packaged into nucleosomal arrays, each nucleosome typically having two copies of histones H2A, H2B, H3, and H4. Histones have distinct posttranslational modifications, variant isoforms, and dynamics. Whether each histone copy within a nucleosome has distinct properties, particularly in relation to the direction of transcription, is unknown. Here we use chromatin immunoprecipitation-exonuclease (ChIP-exo) to resolve the organization of individual histones on a genomic scale. We detect widespread subnucleosomal structures in dynamic chromatin, including what appear to be half-nucleosomes consisting of one copy of each histone. We also detect interactions of H3 tails with linker DNA between nucleosomes, which may be negatively regulated by methylation of H3K36. Histone variant H2A.Z is enriched on the promoter-distal half of the +1 nucleosome, whereas H2BK123 ubiquitylation and H3K9 acetylation are enriched on the promoter-proximal half in a transcription-linked manner. Subnucleosome asymmetries might serve as molecular beacons that guide transcription.

Genome-wide mapping of human DNA replication by optical replication mapping supports a stochastic model of eukaryotic replication.

The heterogeneous nature of eukaryotic replication kinetics and the low efficiency of individual initiation sites make mapping the location and timing of replication initiation in human cells difficult. To address this challenge, we have developed optical replication mapping (ORM), a high-throughput single-molecule approach, and used it to map early-initiation events in human cells. The single-molecule nature of our data and a total of >2,500-fold coverage of the human genome on 27 million fibers averaging ∼300 kb in length allow us to identify initiation sites and their firing probability with high confidence. We find that the distribution of human replication initiation is consistent with inefficient, stochastic activation of heterogeneously distributed potential initiation complexes enriched in accessible chromatin. These observations are consistent with stochastic models of initiation-timing regulation and suggest that stochastic regulation of replication kinetics is a fundamental feature of eukaryotic replication, conserved from yeast to humans.

Esc2 orchestrates substrate-specific sumoylation by acting as a SUMO E2 cofactor in genome maintenance.

SUMO modification regulates diverse cellular processes by targeting hundreds of proteins. However, the limited number of sumoylation enzymes raises the question of how such a large number of substrates are efficiently modified. Specifically, how genome maintenance factors are dynamically sumoylated at DNA replication and repair sites to modulate their functions is poorly understood. Here, we demonstrate a role for the conserved yeast Esc2 protein in this process by acting as a SUMO E2 cofactor. Esc2 is required for genome stability and binds to Holliday junctions and replication fork structures. Our targeted screen found that Esc2 promotes the sumoylation of a Holliday junction dissolution complex and specific replisome proteins. Esc2 does not elicit these effects via stable interactions with substrates or their common SUMO E3. Rather, we show that a SUMO-like domain of Esc2 stimulates sumoylation by exploiting a noncovalent SUMO binding site on the E2 enzyme. This role of Esc2 in sumoylation is required for Holliday junction clearance and genome stability. Our findings thus suggest that Esc2 acts as a SUMO E2 cofactor at distinct DNA structures to promote the sumoylation of specific substrates and genome maintenance.

Genome-scale chromatin binding dynamics of RNA Polymerase II general transcription machinery components.

A great deal of work has revealed, in structural detail, the components of the preinitiation complex (PIC) machinery required for initiation of mRNA gene transcription by RNA polymerase II (Pol II). However, less-well understood are the in vivo PIC assembly pathways and their kinetics, an understanding of which is vital for determining how rates of in vivo RNA synthesis are established. We used competition ChIP in budding yeast to obtain genome-scale estimates of the residence times for five general transcription factors (GTFs): TBP, TFIIA, TFIIB, TFIIE and TFIIF. While many GTF-chromatin interactions were short-lived ( < 1 min), there were numerous interactions with residence times in the range of several minutes. Sets of genes with a shared function also shared similar patterns of GTF kinetic behavior. TFIIE, a GTF that enters the PIC late in the assembly process, had residence times correlated with RNA synthesis rates. The datasets and results reported here provide kinetic information for most of the Pol II-driven genes in this organism, offering a rich resource for exploring the mechanistic relationships between PIC assembly, gene regulation, and transcription.

The Dynamic Fungal Genome: Polyploidy, Aneuploidy and Copy Number Variation in Response to Stress.

Fungal species have dynamic genomes and often exhibit genomic plasticity in response to stress. This genome plasticity often comes with phenotypic consequences that affect fitness and resistance to stress. Fungal pathogens exhibit genome plasticity in both clinical and agricultural settings and often during adaptation to antifungal drugs, posing significant challenges to human health. Therefore, it is important to understand the rates, mechanisms, and impact of large genomic changes. This review addresses the prevalence of polyploidy, aneuploidy, and copy number variation across diverse fungal species, with special attention to prominent fungal pathogens and model species. We also explore the relationship between environmental stress and rates of genomic changes and highlight the mechanisms underlying genotypic and phenotypic changes. A comprehensive understanding of these dynamic fungal genomes is needed to identify novel solutions for the increase in antifungal drug resistance.

Stalled spliceosomes are a signal for RNAi-mediated genome defense.

Using the yeast Cryptococcus neoformans, we describe a mechanism by which transposons are initially targeted for RNAi-mediated genome defense. We show that intron-containing mRNA precursors template siRNA synthesis. We identify a Spliceosome-Coupled And Nuclear RNAi (SCANR) complex required for siRNA synthesis and demonstrate that it physically associates with the spliceosome. We find that RNAi target transcripts are distinguished by suboptimal introns and abnormally high occupancy on spliceosomes. Functional investigations demonstrate that the stalling of mRNA precursors on spliceosomes is required for siRNA accumulation. Lariat debranching enzyme is also necessary for siRNA production, suggesting a requirement for processing of stalled splicing intermediates. We propose that recognition of mRNA precursors by the SCANR complex is in kinetic competition with splicing, thereby promoting siRNA production from transposon transcripts stalled on spliceosomes. Disparity in the strength of expression signals encoded by transposons versus host genes offers an avenue for the evolution of genome defense.

High-resolution mapping reveals a conserved, widespread, dynamic mRNA methylation program in yeast meiosis.

N(6)-methyladenosine (m(6)A) is the most ubiquitous mRNA base modification, but little is known about its precise location, temporal dynamics, and regulation. Here, we generated genomic maps of m(6)A sites in meiotic yeast transcripts at nearly single-nucleotide resolution, identifying 1,308 putatively methylated sites within 1,183 transcripts. We validated eight out of eight methylation sites in different genes with direct genetic analysis, demonstrated that methylated sites are significantly conserved in a related species, and built a model that predicts methylated sites directly from sequence. Sites vary in their methylation profiles along a dense meiotic time course and are regulated both locally, via predictable methylatability of each site, and globally, through the core meiotic circuitry. The methyltransferase complex components localize to the yeast nucleolus, and this localization is essential for mRNA methylation. Our data illuminate a conserved, dynamically regulated methylation program in yeast meiosis and provide an important resource for studying the function of this epitranscriptomic modification.

Divergent genomic trajectories predate the origin of animals and fungi.

Animals and fungi have radically distinct morphologies, yet both evolved within the same eukaryotic supergroup: Opisthokonta1,2. Here we reconstructed the trajectory of genetic changes that accompanied the origin of Metazoa and Fungi since the divergence of Opisthokonta with a dataset that includes four novel genomes from crucial positions in the Opisthokonta phylogeny. We show that animals arose only after the accumulation of genes functionally important for their multicellularity, a tendency that began in the pre-metazoan ancestors and later accelerated in the metazoan root. By contrast, the pre-fungal ancestors experienced net losses of most functional categories, including those gained in the path to Metazoa. On a broad-scale functional level, fungal genomes contain a higher proportion of metabolic genes and diverged less from the last common ancestor of Opisthokonta than did the gene repertoires of Metazoa. Metazoa and Fungi also show differences regarding gene gain mechanisms. Gene fusions are more prevalent in Metazoa, whereas a larger fraction of gene gains were detected as horizontal gene transfers in Fungi and protists, in agreement with the long-standing idea that transfers would be less relevant in Metazoa due to germline isolation3-5. Together, our results indicate that animals and fungi evolved under two contrasting trajectories of genetic change that predated the origin of both groups. The gradual establishment of two clearly differentiated genomic contexts thus set the stage for the emergence of Metazoa and Fungi.

Disome and Trisome Profiling Reveal Genome-wide Targets of Ribosome Quality Control.

The ribosome-associated protein quality control (RQC) system that resolves stalled translation events is activated when ribosomes collide and form disome, trisome, or higher-order complexes. However, it is unclear whether this system distinguishes collision complexes formed on defective mRNAs from those with functional roles on endogenous transcripts. Here, we performed disome and trisome footprint profiling in yeast and found collisions were enriched on diverse sequence motifs known to slow translation. When 60S recycling was inhibited, disomes accumulated at stop codons and could move into the 3' UTR to reinitiate translation. The ubiquitin ligase and RQC factor Hel2/ZNF598 generally recognized collisions but did not induce degradation of endogenous transcripts. However, loss of Hel2 triggered the integrated stress response, via phosphorylation of eIF2α, thus linking these pathways. Our results suggest that Hel2 has a role in sensing ribosome collisions on endogenous mRNAs, and such events may be important for cellular homeostasis.