Saccharomyces Genome Database: the genomics resource of budding yeast

Interpreting protein abundance in Saccharomyces cerevisiae through relational learning

MOTIVATION

Proteomic profiles reflect the functional readout of the physiological state of an organism. An increased understanding of what controls and defines protein abundances is of high scientific interest. Saccharomyces cerevisiae is a well-studied model organism, and there is a large amount of structured knowledge on yeast systems biology in databases such as the Saccharomyces Genome Database, and highly curated genome-scale metabolic models like Yeast8. These datasets, the result of decades of experiments, are abundant in information, and adhere to semantically meaningful ontologies.

RESULTS

By representing this knowledge in an expressive Datalog database we generated data descriptors using relational learning that, when combined with supervised machine learning, enables us to predict protein abundances in an explainable manner. We learnt predictive relationships between protein abundances, function and phenotype; such as α-amino acid accumulations and deviations in chronological lifespan. We further demonstrate the power of this methodology on the proteins His4 and Ilv2, connecting qualitative biological concepts to quantified abundances.

AVAILABILITY AND IMPLEMENTATION

All data and processing scripts are available at the following Github repository: https://github.com/DanielBrunnsaker/ProtPredict.

CryptoCEN: A Co-Expression Network for Cryptococcus neoformans reveals novel proteins involved in DNA damage repair

Elucidating gene function is a major goal in biology, especially among non-model organisms. However, doing so is complicated by the fact that molecular conservation does not always mirror functional conservation, and that complex relationships among genes are responsible for encoding pathways and higher-order biological processes. Co-expression, a promising approach for predicting gene function, relies on the general principal that genes with similar expression patterns across multiple conditions will likely be involved in the same biological process. For Cryptococcus neoformans, a prevalent human fungal pathogen greatly diverged from model yeasts, approximately 60% of the predicted genes in the genome lack functional annotations. Here, we leveraged a large amount of publicly available transcriptomic data to generate a C. neoformans Co-Expression Network (CryptoCEN), successfully recapitulating known protein networks, predicting gene function, and enabling insights into the principles influencing co-expression. With 100% predictive accuracy, we used CryptoCEN to identify 13 new DNA damage response genes, underscoring the utility of guilt-by-association for determining gene function. Overall, co-expression is a powerful tool for uncovering gene function, and decreases the experimental tests needed to identify functions for currently under-annotated genes.

Calcineurin-dependent contributions to fitness in the opportunistic pathogen Candida glabrata

The protein phosphatase calcineurin is vital for the virulence of the opportunistic fungal pathogen Candida glabrata. The host-induced stresses that activate calcineurin signaling are unknown, as are the targets of calcineurin relevant to virulence. To potentially shed light on these processes, millions of transposon insertion mutants throughout the genome of C. glabrata were profiled en masse for fitness defects in the presence of FK506, a specific inhibitor of calcineurin. Eighty-seven specific gene deficiencies depended on calcineurin signaling for full viability in vitro both in wild-type and pdr1∆ null strains lacking pleiotropic drug resistance. Three genes involved in cell wall biosynthesis (FKS1, DCW1, FLC1) possess co-essential paralogs whose expression depended on calcineurin and Crz1 in response to micafungin, a clinical antifungal that interferes with cell wall biogenesis. Interestingly, 80% of the FK506-sensitive mutants were deficient in different aspects of vesicular trafficking, such as endocytosis, exocytosis, sorting, and biogenesis of secretory proteins in the endoplasmic reticulum (ER). In response to the experimental antifungal manogepix that blocks GPI-anchor biosynthesis in the ER, calcineurin signaling increased and strongly prevented cell death independent of Crz1, one of its major targets. Comparisons between manogepix, micafungin, and the ER-stressing tunicamycin reveal a correlation between the degree of calcineurin signaling and the degree of cell survival. These findings suggest that calcineurin plays major roles in mitigating stresses of vesicular trafficking. Such stresses may arise during host infection and in response to antifungal therapies.IMPORTANCECalcineurin plays critical roles in the virulence of most pathogenic fungi. This study sheds light on those roles in the opportunistic pathogen Candida glabrata using a genome-wide analysis in vitro. The findings could lead to antifungal developments that also avoid immunosuppression.

Selected cachaça yeast strains share a genomic profile related to traits relevant to industrial fermentation processes

The isolation and selection of yeast strains to improve the quality of the cachaça-Brazilian Spirit-have been studied in our research group. Our strategy considers Saccharomyces cerevisiae as the predominant species involved in sugarcane juice fermentation and the presence of different stressors (osmolarity, temperature, ethanol content, and competition with other microorganisms). It also considers producing balanced concentrations of volatile compounds (higher alcohols and acetate and/or ethyl esters), flocculation capacity, and ethanol production. Since the genetic bases behind these traits of interest are not fully established, the whole genome sequencing of 11 different Saccharomyces cerevisiae strains isolated and selected from different places was analyzed to identify the presence of a specific genetic variation common to cachaça yeast strains. We have identified 20,128 single-nucleotide variants shared by all genomes. Of these shared variants, 37 were new variants (being six missenses), and 4,451 were identified as missenses. We performed a detailed functional annotation (using enrichment analysis, protein-protein interaction network analysis, and database and in-depth literature searches) of these new and missense variants. Many genes carrying these variations were involved in the phenotypes of flocculation, tolerance to fermentative stresses, and production of volatile compounds and ethanol. These results demonstrate the existence of a genetic profile shared by the 11 strains under study that could be associated with the applied selective strategy. Thus, this study points out genes and variants that may be used as molecular markers for selecting strains well suited to the fermentation process, including genetic improvement by genome editing, ultimately producing high-quality beverages and adding value.IMPORTANCEThis work demonstrates the existence of new genetic markers related to different phenotypes used to select yeast strains and mutations in genes directly involved in producing flavoring compounds and ethanol, and others related to flocculation and stress resistance.

Pseudouridylation-mediated gene expression modulation

RNA-guided pseudouridylation, a widespread post-transcriptional RNA modification, has recently gained recognition for its role in cellular processes such as pre-mRNA splicing and the modulation of premature termination codon (PTC) readthrough. This review provides insights into its mechanisms, functions, and potential therapeutic applications. It examines the mechanisms governing RNA-guided pseudouridylation, emphasizing the roles of guide RNAs and pseudouridine synthases in catalyzing uridine-to-pseudouridine conversion. A key focus is the impact of RNA-guided pseudouridylation of U2 small nuclear RNA on pre-mRNA splicing, encompassing its influence on branch site recognition and spliceosome assembly. Additionally, the review discusses the emerging role of RNA-guided pseudouridylation in regulating PTC readthrough, impacting translation termination and genetic disorders. Finally, it explores the therapeutic potential of pseudouridine modifications, offering insights into potential treatments for genetic diseases and cancer and the development of mRNA vaccine.

Data integration strategies for whole-cell modeling

Data makes the world go round-and high quality data is a prerequisite for precise models, especially for whole-cell models (WCM). Data for WCM must be reusable, contain information about the exact experimental background, and should-in its entirety-cover all relevant processes in the cell. Here, we review basic requirements to data for WCM and strategies how to combine them. As a species-specific resource, we introduce the Yeast Cell Model Data Base (YCMDB) to illustrate requirements and solutions. We discuss recent standards for data as well as for computational models including the modeling process as data to be reported. We outline strategies for constructions of WCM despite their inherent complexity.

The Monarch Initiative in 2024: an analytic platform integrating phenotypes, genes and diseases across species

Bridging the gap between genetic variations, environmental determinants, and phenotypic outcomes is critical for supporting clinical diagnosis and understanding mechanisms of diseases. It requires integrating open data at a global scale. The Monarch Initiative advances these goals by developing open ontologies, semantic data models, and knowledge graphs for translational research. The Monarch App is an integrated platform combining data about genes, phenotypes, and diseases across species. Monarch's APIs enable access to carefully curated datasets and advanced analysis tools that support the understanding and diagnosis of disease for diverse applications such as variant prioritization, deep phenotyping, and patient profile-matching. We have migrated our system into a scalable, cloud-based infrastructure; simplified Monarch's data ingestion and knowledge graph integration systems; enhanced data mapping and integration standards; and developed a new user interface with novel search and graph navigation features. Furthermore, we advanced Monarch's analytic tools by developing a customized plugin for OpenAI's ChatGPT to increase the reliability of its responses about phenotypic data, allowing us to interrogate the knowledge in the Monarch graph using state-of-the-art Large Language Models. The resources of the Monarch Initiative can be found at monarchinitiative.org and its corresponding code repository at github.com/monarch-initiative/monarch-app.

Dissection of the role of a Src homology 3 domain in the evolution of binding preference of paralogous proteins

Protein-protein interactions (PPIs) drive many cellular processes. Some interactions are directed by Src homology 3 (SH3) domains that bind proline-rich motifs on other proteins. The evolution of the binding specificity of SH3 domains is not completely understood, particularly following gene duplication. Paralogous genes accumulate mutations that can modify protein functions and, for SH3 domains, their binding preferences. Here, we examined how the binding of the SH3 domains of 2 paralogous yeast type I myosins, Myo3 and Myo5, evolved following duplication. We found that the paralogs have subtly different SH3-dependent interaction profiles. However, by swapping SH3 domains between the paralogs and characterizing the SH3 domains freed from their protein context, we find that very few of the differences in interactions, if any, depend on the SH3 domains themselves. We used ancestral sequence reconstruction to resurrect the preduplication SH3 domains and examined, moving back in time, how the binding preference changed. Although the most recent ancestor of the 2 domains had a very similar binding preference as the extant ones, older ancestral domains displayed a gradual loss of interaction with the modern interaction partners when inserted in the extant paralogs. Molecular docking and experimental characterization of the free ancestral domains showed that their affinity with the proline motifs is likely not the cause for this loss of binding. Taken together, our results suggest that a SH3 and its host protein could create intramolecular or allosteric interactions essential for the SH3-dependent PPIs, making domains not functionally equivalent even when they have the same binding specificity.

Meta-analysis of dispensable essential genes and their interactions with bypass suppressors

Genes have been historically classified as essential or non-essential based on their requirement for viability. However, genomic mutations can sometimes bypass the requirement for an essential gene, challenging the binary classification of gene essentiality. Such dispensable essential genes represent a valuable model for understanding the incomplete penetrance of loss-of-function mutations often observed in natural populations. Here, we compiled data from multiple studies on essential gene dispensability in Saccharomyces cerevisiae to comprehensively characterize these genes. In analyses spanning different evolutionary timescales, dispensable essential genes exhibited distinct phylogenetic properties compared with other essential and non-essential genes. Integration of interactions with suppressor genes that can bypass the gene essentiality revealed the high functional modularity of the bypass suppression network. Furthermore, dispensable essential and bypass suppressor gene pairs reflected simultaneous changes in the mutational landscape of S. cerevisiae strains. Importantly, species in which dispensable essential genes were non-essential tended to carry bypass suppressor mutations in their genomes. Overall, our study offers a comprehensive view of dispensable essential genes and illustrates how their interactions with bypass suppressors reflect evolutionary outcomes.

The genetic and molecular basis of haploinsufficiency in flowering plants

In diploid organisms, haploinsufficiency can be defined as the requirement for more than one fully functional copy of a gene. In contrast to most genes, whose loss-of-function alleles are recessive, loss-of-function alleles of haploinsufficient genes are dominant. However, forward and reverse genetic screens are biased toward obtaining recessive, loss-of-function mutations, and therefore, dominant mutations of all types are underrepresented in mutant collections. Despite this underrepresentation, haploinsufficient loci have intriguing implications for studies of genome evolution, gene dosage, stability of protein complexes, genetic redundancy, and gene expression. Here we review examples of haploinsufficiency in flowering plants and describe the underlying molecular mechanisms and evolutionary forces driving haploinsufficiency. Finally, we discuss the masking of haploinsufficiency by genetic redundancy, a widespread phenomenon among angiosperms.

A Light in the Dark: Uncovering Wolbachia-Host Interactions Using Fluorescence Imaging

The success of microbial endosymbionts, which reside naturally within a eukaryotic "host" organism, requires effective microbial interaction with, and manipulation of, the host cells. Fluorescence microscopy has played a key role in elucidating the molecular mechanisms of endosymbiosis. For 30 years, fluorescence analyses have been a cornerstone in studies of endosymbiotic Wolbachia bacteria, focused on host colonization, maternal transmission, reproductive parasitism, horizontal gene transfer, viral suppression, and metabolic interactions in arthropods and nematodes. Fluorescence-based studies stand to continue informing Wolbachia-host interactions in increasingly detailed and innovative ways.

Development of a fumonisin-sensitive Saccharomyces cerevisiae indicator strain and utilization for activity testing of candidate detoxification genes

IMPORTANCE

Fumonisins can cause diseases in animals and humans consuming Fusarium-contaminated food or feed. The search for microbes capable of fumonisin degradation, or for enzymes that can detoxify fumonisins, currently relies primarily on chemical detection methods. Our constructed fumonisin B1-sensitive yeast strain can be used to phenotypically detect detoxification activity and should be useful in screening for novel fumonisin resistance genes and to elucidate fumonisin metabolism and resistance mechanisms in fungi and plants, and thereby, in the long term, help to mitigate the threat of fumonisins in feed and food.

Building synthetic chromosomes from natural DNA

De novo chromosome synthesis is costly and time-consuming, limiting its use in research and biotechnology. Building synthetic chromosomes from natural components is an unexplored alternative with many potential applications. In this paper, we report CReATiNG (Cloning, Reprogramming, and Assembling Tiled Natural Genomic DNA), a method for constructing synthetic chromosomes from natural components in yeast. CReATiNG entails cloning segments of natural chromosomes and then programmably assembling them into synthetic chromosomes that can replace the native chromosomes in cells. We use CReATiNG to synthetically recombine chromosomes between strains and species, to modify chromosome structure, and to delete many linked, non-adjacent regions totaling 39% of a chromosome. The multiplex deletion experiment reveals that CReATiNG also enables recovery from flaws in synthetic chromosome design via recombination between a synthetic chromosome and its native counterpart. CReATiNG facilitates the application of chromosome synthesis to diverse biological problems.

Characterizing Different Modes of Interplay Between Rap1 and H3 Using Inducible H3-depletion Yeast

Histones and transcription factors (TFs) are two important DNA-binding proteins that interact, compete, and together regulate transcriptional processes in response to diverse internal and external stimuli. Condition-specific depletion of histones in Saccharomyces cerevisiae using a galactose-inducible H3 promoter provides a suitable framework for examining transcriptional alteration resulting from reduced nucleosome content. However, the effect on DNA binding activities of TFs is yet to be fully explored. In this work, we combine ChIP-seq of H3 with RNA-seq to elucidate the genome-scale relationships between H3 occupancy patterns and transcriptional dynamics before and after global H3 depletion. ChIP-seq of Rap1 is also conducted in the H3-depletion and control treatments, to investigate the interplay between this master regulator TF and nucleosomal H3, and to explore the impact on diverse transcriptional responses of different groups of target genes and functions. Ultimately, we propose a working model and testable hypotheses regarding the impact of global and local H3 depletion on transcriptional changes. We also demonstrate different potential modes of interaction between Rap1 and H3, which sheds light on the potential multifunctional regulatory capabilities of Rap1 and potentially other pioneer factors.

Potency of CRISPR-Cas Antifungals Is Enhanced by Cotargeting DNA Repair and Growth Regulatory Machinery at the Genetic Level

The emergence of virulent, resistant, and rapidly evolving fungal pathogens poses a significant threat to public health, agriculture, and the environment. Targeting cellular processes with standard small-molecule intervention may be effective but requires long development times and is prone to antibiotic resistance. To overcome the current limitations of antibiotic development and treatment, this study harnesses CRISPR-Cas systems as antifungals by capitalizing on their adaptability, specificity, and efficiency in target design. The conventional design of CRISPR-Cas antimicrobials, based on induction of DNA double-strand breaks (DSBs), is potentially less effective in fungi due to robust eukaryotic DNA repair machinery. Here, we report a novel design principle to formulate more effective CRISPR-Cas antifungals by cotargeting essential genes with DNA repair defensive genes that remove the fungi's ability to repair the DSB sites of essential genes. By evaluating this design on the model fungus Saccharomyces cerevisiae, we demonstrated that essential and defensive gene cotargeting is more effective than either essential or defensive gene targeting alone. The top-performing CRISPR-Cas antifungals performed as effectively as the antibiotic Geneticin. A gene cotargeting interaction analysis revealed that cotargeting essential genes with RAD52 involved in homologous recombination (HR) was the most synergistic combination. Fast growth kinetics of S. cerevisiae induced resistance to CRISPR-Cas antifungals, where genetic mutations mostly occurred in defensive genes and guide RNA sequences.

John N, Zaidi F, Doyle F, Fasullo M. High-throughput screening of the Saccharomyces cerevisiae genome for 2-amino-3-methylimidazo [4,5-f] quinoline resistance identifies colon cancer-associated genes

Heterocyclic aromatic amines (HAAs) are potent carcinogenic agents found in charred meats and cigarette smoke. However, few eukaryotic resistance genes have been identified. We used Saccharomyces cerevisiae (budding yeast) to identify genes that confer resistance to 2-amino-3-methylimidazo[4,5-f] quinoline (IQ). CYP1A2 and NAT2 activate IQ to become a mutagenic nitrenium compound. Deletion libraries expressing human CYP1A2 and NAT2 or no human genes were exposed to either 400 or 800 µM IQ for 5 or 10 generations. DNA barcodes were sequenced using the Illumina HiSeq 2500 platform and statistical significance was determined for exactly matched barcodes. We identified 424 ORFs, including 337 genes of known function, in duplicate screens of the "humanized" collection for IQ resistance; resistance was further validated for a select group of 51 genes by growth curves, competitive growth, or trypan blue assays. Screens of the library not expressing human genes identified 143 ORFs conferring resistance to IQ per se. Ribosomal protein and protein modification genes were identified as IQ resistance genes in both the original and "humanized" libraries, while nitrogen metabolism, DNA repair, and growth control genes were also prominent in the "humanized" library. Protein complexes identified included the casein kinase 2 (CK2) and histone chaperone (HIR) complex. Among DNA Repair and checkpoint genes, we identified those that function in postreplication repair (RAD18, UBC13, REV7), base excision repair (NTG1), and checkpoint signaling (CHK1, PSY2). These studies underscore the role of ribosomal protein genes in conferring IQ resistance, and illuminate DNA repair pathways for conferring resistance to activated IQ.

Quantifying turgor pressure in budding and fission yeasts based upon osmotic properties

Walled cells, such as plants, fungi, and bacteria cells, possess a high internal hydrostatic pressure, termed turgor pressure, that drives volume growth and contributes to cell shape determination. Rigorous measurement of turgor pressure, however, remains challenging, and reliable quantitative measurements, even in budding yeast are still lacking. Here, we present a simple and robust experimental approach to access turgor pressure in yeasts based upon the determination of isotonic concentration using protoplasts as osmometers. We propose three methods to identify the isotonic condition - three-dimensional cell volume, cytoplasmic fluorophore intensity, and mobility of a cytGEMs nano-rheology probe - that all yield consistent values. Our results provide turgor pressure estimates of 1.0 ± 0.1 MPa for Schizosaccharomyces pombe, 0.49 ± 0.01 MPa for Schizosaccharomyces japonicus, 0.5 ± 0.1 MPa for Saccharomyces cerevisiae W303a and 0.31 ± 0.03 MPa for Saccharomyces cerevisiae BY4741. Large differences in turgor pressure and nano-rheology measurements between the Saccharomyces cerevisiae strains demonstrate how fundamental biophysical parameters can vary even among wild-type strains of the same species. These side-by-side measurements of turgor pressure in multiple yeast species provide critical values for quantitative studies on cellular mechanics and comparative evolution.

Biological factors and statistical limitations prevent detection of most noncanonical proteins by mass spectrometry

Ribosome profiling experiments indicate pervasive translation of short open reading frames (ORFs) outside of annotated protein-coding genes. However, shotgun mass spectrometry (MS) experiments typically detect only a small fraction of the predicted protein products of this noncanonical translation. The rarity of detection could indicate that most predicted noncanonical proteins are rapidly degraded and not present in the cell; alternatively, it could reflect technical limitations. Here, we leveraged recent advances in ribosome profiling and MS to investigate the factors limiting detection of noncanonical proteins in yeast. We show that the low detection rate of noncanonical ORF products can largely be explained by small size and low translation levels and does not indicate that they are unstable or biologically insignificant. In particular, proteins encoded by evolutionarily young genes, including those with well-characterized biological roles, are too short and too lowly expressed to be detected by shotgun MS at current detection sensitivities. Additionally, we find that decoy biases can give misleading estimates of noncanonical protein false discovery rates, potentially leading to false detections. After accounting for these issues, we found strong evidence for 4 noncanonical proteins in MS data, which were also supported by evolution and translation data. These results illustrate the power of MS to validate unannotated genes predicted by ribosome profiling, but also its substantial limitations in finding many biologically relevant lowly expressed proteins.

The interactome of the UapA transporter reveals putative new players in anterograde membrane cargo trafficking

Neosynthesized plasma membrane (PM) proteins co-translationally translocate to the ER, concentrate at regions called ER-exit sites (ERes) and pack into COPII secretory vesicles which are sorted to the early-Golgi through membrane fusion. Following Golgi maturation, membrane cargoes reach the late-Golgi, from where they exit in clathrin-coated vesicles destined to the PM, directly or through endosomes. Post-Golgi membrane cargo trafficking also involves the cytoskeleton and the exocyst. The Golgi-dependent secretory pathway is thought to be responsible for the trafficking of all major membrane proteins. However, our recent findings in Aspergillus nidulans showed that several plasma membrane cargoes, such as transporters and receptors, follow a sorting route that seems to bypass Golgi functioning. To gain insight on membrane trafficking and specifically Golgi-bypass, here we used proximity dependent biotinylation (PDB) coupled with data-independent acquisition mass spectrometry (DIA-MS) for identifying transient interactors of the UapA transporter. Our assays, which included proteomes of wild-type and mutant strains affecting ER-exit or endocytosis, identified both expected and novel interactions that might be physiologically relevant to UapA trafficking. Among those, we validated, using reverse genetics and fluorescence microscopy, that COPI coatomer is essential for ER-exit and anterograde trafficking of UapA and other membrane cargoes. We also showed that ArfAArf1 GTPase activating protein (GAP) Glo3 contributes to UapA trafficking at increased temperature. This is the first report addressing the identification of transient interactions during membrane cargo biogenesis using PDB and proteomics coupled with fungal genetics. Our work provides a basis for dissecting dynamic membrane cargo trafficking via PDB assays.

Genomic analysis of the marine yeast Rhodotorula sphaerocarpa ETNP2018 reveals adaptation to the open ocean

BACKGROUND

Despite a rising interest in the diversity and ecology of fungi in marine environments, there are few published genomes of fungi isolated from the ocean. The basidiomycetous yeast (unicellular fungus) genus Rhodotorula are prevalent and abundant in the open ocean, and they have been isolated from a wide range of other environments. Many of these environments are nutrient poor, such as the Antarctica and the Atacama deserts, raising the question as to how Rhodotorula yeasts may have adapted their metabolic strategies to optimize survival under low nutrient conditions. In order to understand their adaptive strategies in the ocean, the genome of R. sphaerocarpa ETNP2018 was compared to that of fourteen representative Rhodotorula yeasts, isolated from a variety of environments.

RESULTS

Rhodotorula sphaerocarpa ETNP2018, a strain isolated from the oligotrophic part of the eastern tropical North Pacific (ETNP) oxygen minimum zone (OMZ), hosts the smallest of the fifteen genomes and yet the number of protein-coding genes it possesses is on par with the other strains. Its genome exhibits a distinct reduction in genes dedicated to Major Facilitator Superfamily transporters as well as biosynthetic enzymes. However, its core metabolic pathways are fully conserved. Our research indicates that the selective pressures of the ETNP OMZ favor a streamlined genome with reduced overall biosynthetic potential balanced by a stable set of core metabolisms and an expansion of mechanisms for nutrient acquisition.

CONCLUSIONS

In summary, this study offers insights into the adaptation of fungi to the oligotrophic ocean and provides valuable information for understanding the ecological roles of fungi in the ocean.

Persistent Acetylation of Histone H3 Lysine 56 Compromises the Activity of DNA Replication Origins

In Saccharomyces cerevisiae, newly synthesized histones H3 are acetylated on lysine 56 (H3 K56ac) by the Rtt109 acetyltransferase prior to their deposition on nascent DNA behind replication forks. Two deacetylases of the sirtuin family, Hst3 and Hst4, remove H3 K56ac from chromatin after S phase. hst3Δ hst4Δ cells present constitutive H3 K56ac, which sensitizes cells to replicative stress via unclear mechanisms. A chemogenomic screen revealed that DBF4 heterozygosity sensitizes cells to NAM-induced inhibition of sirtuins. DBF4 and CDC7 encode subunits of the Dbf4-dependent kinase (DDK), which activates origins of DNA replication during S phase. We show that (i) cells harboring the dbf4-1 or cdc7-4 hypomorphic alleles are sensitized to NAM, and that (ii) the sirtuins Sir2, Hst1, Hst3, and Hst4 promote DNA replication in cdc7-4 cells. We further demonstrate that Rif1, an inhibitor of DDK-dependent activation of origins, causes DNA damage and replication defects in NAM-treated cells and hst3Δ hst4Δ mutants. cdc7-4 hst3Δ hst4Δ cells are shown to display delayed initiation of DNA replication, which is not due to intra-S checkpoint activation but requires Rtt109-dependent H3 K56ac. Our results suggest that constitutive H3 K56ac sensitizes cells to replicative stress in part by negatively influencing the activation of origins of DNA replication.

JL. Additional principles that govern the release of pre-ribosomes from the nucleolus into the nucleoplasm in yeast

During eukaryotic ribosome biogenesis, pre-ribosomes travel from the nucleolus, where assembly is initiated, to the nucleoplasm and then are exported to the cytoplasm, where assembly concludes. Although nuclear export of pre-ribosomes has been extensively investigated, the release of pre-ribosomes from the nucleolus is an understudied phenomenon. Initial data indicate that unfolded rRNA interacts in trans with nucleolar components and that, when rRNA folds due to ribosomal protein (RP) binding, the number of trans interactions drops below the threshold necessary for nucleolar retention. To validate and expand on this idea, we performed a bioinformatic analysis of the protein components of the Saccharomyces cerevisiae ribosome assembly pathway. We found that ribosome biogenesis factors (RiBi factors) contain significantly more predicted trans interacting regions than RPs. We also analyzed cryo-EM structures of ribosome assembly intermediates to determine how nucleolar pre-ribosomes differ from post-nucleolar pre-ribosomes, specifically the capacity of RPs, RiBi factors, and rRNA components to interact in trans. We observed a significant decrease in the theoretical trans-interacting capability of pre-ribosomes between nucleolar and post-nucleolar stages of assembly due to the release of RiBi factors from particles and the folding of rRNA. Here, we provide a mechanism for the release of pre-ribosomes from the nucleolus.

Model-based inference of a plant-specific dual role for HOPS in regulating guard cell vacuole fusion

Stomata are the pores on a leaf surface that regulate gas exchange. Each stoma consists of two guard cells whose movements regulate pore opening and thereby control CO2 fixation and water loss. Guard cell movements depend in part on the remodeling of vacuoles, which have been observed to change from a highly fragmented state to a fused morphology during stomata opening. This change in morphology requires a membrane fusion mechanism that responds rapidly to environmental signals, allowing plants to respond to diurnal and stress cues. With guard cell vacuoles being both large and responsive to external signals, stomata represent a unique system in which to delineate mechanisms of membrane fusion. Fusion of vacuole membranes is a highly conserved process in eukaryotes, with key roles played by two multi-subunit complexes: HOPS (homotypic fusion and vacuolar protein sorting) and SNARE (soluble NSF attachment protein receptor). HOPS is a vacuole tethering factor that is thought to chaperone SNAREs from apposing vacuole membranes into a fusion-competent complex capable of rearranging membranes. To resolve a counter-intuitive observation regarding the role of HOPS in regulating plant vacuole morphology, we derived a quantitative model of vacuole fusion dynamics and used it to generate testable predictions about HOPS-SNARE interactions. We derived our model by applying simulation-based inference to integrate prior knowledge about molecular interactions with limited, qualitative observations of emergent vacuole phenotypes. By constraining the model parameters to yield the emergent outcomes observed for stoma opening - as induced by two distinct chemical treatments - we predicted a dual role for HOPS and identified a stalled form of the SNARE complex that differs from phenomena reported in yeast. We predict that HOPS has contradictory actions at different points in the fusion signaling pathway, promoting the formation of SNARE complexes, but limiting their activity.

Synthetic yeast chromosome XI design provides a testbed for the study of extrachromosomal circular DNA dynamics

We describe construction of the synthetic yeast chromosome XI (synXI) and reveal the effects of redesign at non-coding DNA elements. The 660-kb synthetic yeast genome project (Sc2.0) chromosome was assembled from synthesized DNA fragments before CRISPR-based methods were used in a process of bug discovery, redesign, and chromosome repair, including precise compaction of 200 kb of repeat sequence. Repaired defects were related to poor centromere function and mitochondrial health and were associated with modifications to non-coding regions. As part of the Sc2.0 design, loxPsym sequences for Cre-mediated recombination are inserted between most genes. Using the GAP1 locus from chromosome XI, we show that these sites can facilitate induced extrachromosomal circular DNA (eccDNA) formation, allowing direct study of the effects and propagation of these important molecules. Construction and characterization of synXI contributes to our understanding of non-coding DNA elements, provides a useful tool for eccDNA study, and will inform future synthetic genome design.

Hsp90 shapes adaptation by controlling the fitness consequences of regulatory variation

The essential stress-responsive chaperone Hsp90 impacts development and adaptation from microbes to humans. Yet despite evidence of its role in evolution, pathogenesis, and oncogenic transformation, the molecular mechanisms by which Hsp90 alters the consequences of mutations remain vigorously debated. Here we exploit the power of nucleotide-resolution genetic mapping in Saccharomyces cerevisiae to uncover more than 1,000 natural variant-to-phenotype associations governed by this molecular chaperone. Strikingly, Hsp90 more frequently modified the phenotypic effects of cis-regulatory variation than variants that altered protein sequence. Moreover, these interactions made the largest contribution to Hsp90-dependent heredity. Nearly all interacting variants-both regulatory and protein-coding-fell within clients of Hsp90 or targets of its direct binding partners. Hsp90 activity affected mutations in evolutionarily young genes, segmental deletions, and heterozygotes, highlighting its influence on variation central to evolutionary novelty. Reconciling the diverse epistatic effects of this chaperone, synthetic transcriptional regulation and reconstructions of natural alleles by genome editing revealed a central role for Hsp90 in regulating the fundamental relationship between activity and phenotype. Our findings establish that non-coding variation is a core driver of Hsp90's influence on heredity, offering a mechanistic explanation for the chaperone's strong effects on evolution and development across species.

Only One of Three Bcs1 Homologs in Aspergillus fumigatus Confers Respiratory Growth

The mitochondrial translocase Bcs1 is required for the correct assembly of complex III of the mitochondrial respiratory chain. Because of its importance, Bcs1 was recently proposed as a target for antifungal agents. The function of this AAA (ATPase Associated with diverse cellular Activities) protein has been extensively characterized in Saccharomyces cerevisiae. This yeast as well as previously studied mammals each encode only one homolog. In contrast, the pathogenic mold Aspergillus fumigatus encodes three putative Bcs1 homologs, none of which have been characterized to date. To study the role of these three homologs in A. fumigatus, conditional and deletion mutants of the respective genes AFUA_3G13000 (bcs1A), AFUA_4G01260 (bcs1B), and AFUA_2G14760 (bcs1C) were generated. A deletion or downregulation of bcs1A resulted in drastically reduced growth and sporulation rates and in a significantly altered susceptibility to azole antifungals. In contrast, mutants lacking Bcs1B or Bcs1C did not show any phenotypes differing from the wild type. Salicylhydroxamic acid-an inhibitor of the alternative oxidase that allows the respiratory chain to bypass complex III in some species-caused a complete growth arrest of the bcs1A deletion mutant. In a Galleria mellonella infection model, the deletion of bcs1A resulted in significantly decreased virulence. Only Bcs1A was able to partially complement a deletion of BCS1 in S. cerevisiae. The subcellular localization of Bcs1B and Bcs1C outside of mitochondria suggests that these Bcs1 homologs exert cellular functions different from that of Bcs1. Our data demonstrate that Bcs1A is the sole Bcs1 ortholog in A. fumigatus.

Computational analysis on the dissemination of non-B DNA structural motifs in promoter regions of 1180 cellular genomes

The promoter regions of gene regulation are under evolutionary constraints and earlier studies uncovered that they are characterized by enrichment of functional non-B DNA structural signatures like curved DNA, cruciform DNA, G-quadruplex, triple-helical DNA, slipped DNA structures, and Z-DNA. However, these studies are restricted to a few model organisms, single non-B DNA motif types, or whole genomic sequences, and their comparative accumulation in promoter regions of different domains of life has not been reported comprehensively. In this study, for the first time, we investigated the preponderance of non-B DNA-prone motifs in promoter regions in 1180 genomes belonging to 28 taxonomic groups using the non-B DNA Motif Search Tool (nBMST). The trends suggest that they are predominant in promoters compared to the upstream and downstream regions of all three domains of life and variably linked to taxonomic groups. Cruciform DNA motif is the most abundant form of non-B DNA, spanning from archaea to lower eukaryotes. Curved DNA motifs are prominent in host-associated bacteria, and suppressed in mammals. Triplex-DNA and slipped DNA structure repeats are discretely dispersed in all lineages. G-quadruplex motifs are significantly enriched in mammals. We also observed that the unique enrichment of non-B DNA in promoters is strongly linked to genome GC, size, evolutionary time divergence, and ecological adaptations. Overall, our work systematically reports the unique non-B DNA structural landscape of cellular organisms from the perspective of the cis-regulatory code of genomes.

Comparative analysis of mitochondrion-related organelles in anaerobic amoebozoans

Archamoebae comprises free-living or endobiotic amoebiform protists that inhabit anaerobic or microaerophilic environments and possess mitochondrion-related organelles (MROs) adapted to function anaerobically. We compared in silico reconstructed MRO proteomes of eight species (six genera) and found that the common ancestor of Archamoebae possessed very few typical components of the protein translocation machinery, electron transport chain and tricarboxylic acid cycle. On the other hand, it contained a sulphate activation pathway and bacterial iron-sulphur (Fe-S) assembly system of MIS-type. The metabolic capacity of the MROs, however, varies markedly within this clade. The glycine cleavage system is widely conserved among Archamoebae, except in Entamoeba, probably owing to its role in catabolic function or one-carbon metabolism. MRO-based pyruvate metabolism was dispensed within subgroups Entamoebidae and Rhizomastixidae, whereas sulphate activation could have been lost in isolated cases of Rhizomastix libera, Mastigamoeba abducta and Endolimax sp. The MIS (Fe-S) assembly system was duplicated in the common ancestor of Mastigamoebidae and Pelomyxidae, and one of the copies took over Fe-S assembly in their MRO. In Entamoebidae and Rhizomastixidae, we hypothesize that Fe-S cluster assembly in both compartments may be facilitated by dual localization of the single system. We could not find evidence for changes in metabolic functions of the MRO in response to changes in habitat; it appears that such environmental drivers do not strongly affect MRO reduction in this group of eukaryotes.

Biological Factors and Statistical Limitations Prevent Detection of Most Noncanonical Proteins by Mass Spectrometry

Ribosome profiling experiments indicate pervasive translation of short open reading frames (ORFs) outside of annotated protein-coding genes. However, shotgun mass spectrometry experiments typically detect only a small fraction of the predicted protein products of this noncanonical translation. The rarity of detection could indicate that most predicted noncanonical proteins are rapidly degraded and not present in the cell; alternatively, it could reflect technical limitations. Here we leveraged recent advances in ribosome profiling and mass spectrometry to investigate the factors limiting detection of noncanonical proteins in yeast. We show that the low detection rate of noncanonical ORF products can largely be explained by small size and low translation levels and does not indicate that they are unstable or biologically insignificant. In particular, proteins encoded by evolutionarily young genes, including those with well-characterized biological roles, are too short and too lowly-expressed to be detected by shotgun mass spectrometry at current detection sensitivities. Additionally, we find that decoy biases can give misleading estimates of noncanonical protein false discovery rates, potentially leading to false detections. After accounting for these issues, we found strong evidence for four noncanonical proteins in mass spectrometry data, which were also supported by evolution and translation data. These results illustrate the power of mass spectrometry to validate unannotated genes predicted by ribosome profiling, but also its substantial limitations in finding many biologically relevant lowly-expressed proteins.

Yeast metagenomics: analytical challenges in the analysis of the eukaryotic microbiome

Even if their impact is often underestimated, yeasts and yeast-like fungi represent the most prevalent eukaryotic members of microbial communities on Earth. They play numerous roles in natural ecosystems and in association with their hosts. They are involved in the food industry and pharmaceutical production, but they can also cause diseases in other organisms, making the understanding of their biology mandatory. The ongoing loss of biodiversity due to overexploitation of environmental resources is a growing concern in many countries. Therefore, it becomes crucial to understand the ecology and evolutionary history of these organisms to systematically classify them. To achieve this, it is essential that our knowledge of the mycobiota reaches a level similar to that of the bacterial communities. To overcome the existing challenges in the study of fungal communities, the first step should be the establishment of standardized techniques for the correct identification of species, even from complex matrices, both in wet lab practices and in bioinformatic tools.

Comparative analyses of disease-linked missense mutations in the RNA exosome modeled in budding yeast reveal distinct functional consequences in translation

The RNA exosome is an evolutionarily conserved exoribonuclease complex that consists of a 3-subunit cap, a 6-subunit barrel-shaped core, and a catalytic base subunit. Missense mutations in genes encoding structural subunits of the RNA exosome cause a growing family of diseases with diverse pathologies, collectively termed RNA exosomopathies. The disease symptoms vary and can manifest as neurological defects or developmental disorders. The diversity of the RNA exosomopathy pathologies suggests that the different missense mutations in structural genes result in distinct in vivo consequences. To investigate these functional consequences and distinguish whether they are unique to each RNA exosomopathy mutation, we generated a collection of in vivo models using budding yeast by introducing pathogenic missense mutations in orthologous S. cerevisiae genes. We then performed a comparative RNA-seq analysis to assess broad transcriptomic changes in each mutant model. Three of the mutant models rrp4-G226D, rrp40-W195R and rrp46-L191H, which model mutations in the genes encoding structural subunits of the RNA exosome, EXOSC2, EXOSC3 and EXOSC5 showed the largest transcriptomic differences. Further analyses revealed shared increased transcripts enriched in translation or ribosomal RNA modification/processing pathways across the three mutant models. Studies of the impact of the mutations on translation revealed shared defects in ribosome biogenesis but distinct impacts on translation. Collectively, our results provide the first comparative analysis of several RNA exosomopathy mutant models and suggest that different RNA exosomopathy mutations result in in vivo consequences that are both unique and shared across each variant, providing more insight into the biology underlying each distinct pathology.

Effects of HSP70 chaperones Ssa1 and Ssa2 on Ste5 scaffold and the mating mitogen-activated protein kinase (MAPK) pathway in Saccharomyces cerevisiae

Ste5 is a prototype of scaffold proteins that regulate activation of mitogen-activated protein kinase (MAPK) cascades in all eukaryotes. Ste5 associates with many proteins including Gβγ (Ste4), Ste11 MAPKKK, Ste7 MAPKK, Fus3 and Kss1 MAPKs, Bem1, Cdc24. Here we show that Ste5 also associates with heat shock protein 70 chaperone (Hsp70) Ssa1 and that Ssa1 and its ortholog Ssa2 are together important for Ste5 function and efficient mating responses. The majority of purified overexpressed Ste5 associates with Ssa1. Loss of Ssa1 and Ssa2 has deleterious effects on Ste5 abundance, integrity, and localization particularly when Ste5 is expressed at native levels. The status of Ssa1 and Ssa2 influences Ste5 electrophoresis mobility and formation of high molecular weight species thought to be phosphorylated, ubiquitinylated and aggregated and lower molecular weight fragments. A Ste5 VWA domain mutant with greater propensity to form punctate foci has reduced predicted propensity to bind Ssa1 near the mutation sites and forms more punctate foci when Ssa1 Is overexpressed, supporting a dynamic protein quality control relationship between Ste5 and Ssa1. Loss of Ssa1 and Ssa2 reduces activation of Fus3 and Kss1 MAPKs and FUS1 gene expression and impairs mating shmoo morphogenesis. Surprisingly, ssa1, ssa2, ssa3 and ssa4 single, double and triple mutants can still mate, suggesting compensatory mechanisms exist for folding. Additional analysis suggests Ssa1 is the major Hsp70 chaperone for the mating and invasive growth pathways and reveals several Hsp70-Hsp90 chaperone-network proteins required for mating morphogenesis.

Cellular Stress Impact on Yeast Activity in Biotechnological Processes-A Short Overview

The importance of Saccharomyces cerevisiae yeast cells is known worldwide, as they are the most used microorganisms in biotechnology for bioethanol and biofuel production. Also, they are analyzed and studied for their similar internal biochemical processes to human cells, for a better understanding of cell aging and response to cell stressors. The special ability of S. cerevisiae cells to develop in both aerobic and anaerobic conditions makes this microorganism a viable model to study the transformations and the way in which cellular metabolism is directed to face the stress conditions due to environmental changes. Thus, this review will emphasize the effects of oxidative, ethanol, and osmotic stress and also the physiological and genetic response of stress mitigation in yeast cells.

Contribution of Spontaneous Mutations to Quantitative and Molecular Variation at the Highly Repetitive rDNA Locus in Yeast

The ribosomal DNA array in Saccharomyces cerevisiae consists of many tandem repeats whose copy number is believed to be functionally important but highly labile. Regulatory mechanisms have evolved to maintain copy number by directed mutation, but how spontaneous variation at this locus is generated and selected has not been well characterized. We applied a mutation accumulation approach to quantify the impacts of mutation and selection on this unique genomic feature across hundreds of mutant strains. We find that mutational variance for this trait is relatively high, and that unselected mutations elsewhere in the genome can disrupt copy number maintenance. In consequence, copy number generally declines gradually, consistent with a previously proposed model of rDNA maintenance where a downward mutational bias is normally compensated by mechanisms that increase copy number when it is low. This pattern holds across ploidy levels and strains in the standard lab environment but differs under some stressful conditions. We identify several alleles, gene categories, and genomic features that likely affect copy number, including aneuploidy for chromosome XII. Copy number change is associated with reduced growth in diploids, consistent with stabilizing selection. Levels of standing variation in copy number are well predicted by a balance between mutation and stabilizing selection, suggesting this trait is not subject to strong diversifying selection in the wild. The rate and spectrum of point mutations within the rDNA locus itself are distinct from the rest of the genome and predictive of polymorphism locations. Our findings help differentiate the roles of mutation and selection and indicate that spontaneous mutation patterns shape several aspects of ribosomal DNA evolution.

Quantifying turgor pressure in budding and fission yeasts based upon osmotic properties

Walled cells, such as plants, fungi, and bacteria cells, possess a high internal hydrostatic pressure, termed turgor pressure, that drives volume growth and contributes to cell shape determination. Rigorous measurement of turgor pressure, however, remains challenging, and reliable quantitative measurements, even in budding yeast are still lacking. Here, we present a simple and robust experimental approach to access turgor pressure in yeasts based upon the determination of isotonic concentration using protoplasts as osmometers. We propose three methods to identify the isotonic condition - 3D cell volume, cytoplasmic fluorophore intensity, and mobility of a cytGEMs nano-rheology probe - that all yield consistent values. Our results provide turgor pressure estimates of 1.0 ± 0.1 MPa for S. pombe, 0.49 ± 0.01 MPa for S. japonicus, 0.5 ± 0.1 MPa for S. cerevisiae W303a and 0.31 ± 0.03 MPa for S. cerevisiae BY4741. Large differences in turgor pressure and nano-rheology measurements between the S. cerevisiae strains demonstrate how fundamental biophysical parameters can vary even among wildtype strains of the same species. These side-by-side measurements of turgor pressure in multiple yeast species provide critical values for quantitative studies on cellular mechanics and comparative evolution.

Visioning synthetic futures for yeast research within the context of current global techno-political trends

Yeast research is entering into a new period of scholarship, with new scientific tools, new questions to ask and new issues to consider. The politics of emerging and critical technology can no longer be separated from the pursuit of basic science in fields, such as synthetic biology and engineering biology. Given the intensifying race for technological leadership, yeast research is likely to attract significant investment from government, and that it offers huge opportunities to the curious minded from a basic research standpoint. This article provides an overview of new directions in yeast research with a focus on Saccharomyces cerevisiae, and places these trends in their geopolitical context. At the highest level, yeast research is situated within the ongoing convergence of the life sciences with the information sciences. This convergent effect is most strongly pronounced in areas of AI-enabled tools for the life sciences, and the creation of synthetic genomes, minimal genomes, pan-genomes, neochromosomes and metagenomes using computer-assisted design tools and methodologies. Synthetic yeast futures encompass basic and applied science questions that will be of intense interest to government and nongovernment funding sources. It is essential for the yeast research community to map and understand the context of their research to ensure their collaborations turn global challenges into research opportunities.

Systematic Profiling of Ale Yeast Protein Dynamics across Fermentation and Repitching

Studying the genetic and molecular characteristics of brewing yeast strains is crucial for understanding their domestication history and adaptations accumulated over time in fermentation environments, and for guiding optimizations to the brewing process itself. Saccharomyces cerevisiae (brewing yeast) is amongst the most profiled organisms on the planet, yet the temporal molecular changes that underlie industrial fermentation and beer brewing remain understudied. Here, we characterized the genomic makeup of a Saccharomyces cerevisiae ale yeast widely used in the production of Hefeweizen beers, and applied shotgun mass spectrometry to systematically measure the proteomic changes throughout two fermentation cycles which were separated by 14 rounds of serial repitching. The resulting brewing yeast proteomics resource includes 64,740 protein abundance measurements. We found that this strain possesses typical genetic characteristics of Saccharomyces cerevisiae ale strains and displayed progressive shifts in molecular processes during fermentation based on protein abundance changes. We observed protein abundance differences between early fermentation batches compared to those separated by 14 rounds of serial repitching. The observed abundance differences occurred mainly in proteins involved in the metabolism of ergosterol and isobutyraldehyde. Our systematic profiling serves as a starting point for deeper characterization of how the yeast proteome changes during commercial fermentations and additionally serves as a resource to guide fermentation protocols, strain handling, and engineering practices in commercial brewing and fermentation environments. Finally, we created a web interface (https://brewing-yeast-proteomics.ccbb.utexas.edu/) to serve as a valuable resource for yeast geneticists, brewers, and biochemists to provide insights into the global trends underlying commercial beer production.

Human SHQ1 variants R335C and A426V lead to severe ribosome biogenesis defects when expressed in yeast

SHQ1 is an essential chaperone that binds the pseudouridine synthase dyskerin in the cytoplasm and escorts the enzyme to the nucleus, where dyskerin is assembled into small nucleolar RNPs (snoRNPs) of the H/ACA class. These particles carry out pseudouridine formation in ribosomal RNAs (rRNAs) and participate in maturation of rRNA precursors (pre-rRNAs). Variants of human SHQ1 have been linked to neurodevelopmental deficiencies; here we focused on two compound heterozygous mutations identified in a child showing a severe neurological disorder comprising cerebellar degeneration. To investigate the molecular defects caused by mutations R335C and A426V we used a conditional yeast strain that can be depleted of the endogenous Shq1 protein while constitutively expressing human SHQ1 (wild-type or variants). Although wild-type SHQ1 complemented the Shq1-depleted strain, cells expressing variant R335C could not support growth, and cells expressing variant A426V were temperature-sensitive. When shifted to restrictive conditions, yeast cells progressively lost H/ACA snoRNAs and accumulated unprocessed pre-rRNAs, which led to reduced production of ribosomes. Levels of Cbf5 (yeast homologue of dyskerin) were decreased in yeast cells expressing SHQ1 variants under restrictive conditions. Immunoprecipitation experiments revealed that interaction of Cbf5 with SHQ1 variants was weakened but not abolished, and yeast two-hybrid assays showed that mutation R335C is more deleterious than mutation A426V. Our data provide additional evidence for the critical role of SHQ1 in chaperoning the pseudouridine synthase dyskerin, and how its inadequate function has detrimental consequences on the production of H/ACA snoRNPs and ribosomes.

The architecture of binding cooperativity between densely bound transcription factors

The binding of transcription factors (TFs) along genomes is restricted to a subset of sites containing their preferred motifs. TF-binding specificity is often attributed to the co-binding of interacting TFs; however, apart from specific examples, this model remains untested. Here, we define dependencies among budding yeast TFs that localize to overlapping promoters by profiling the genome-wide consequences of co-depleting multiple TFs. We describe unidirectional interactions, revealing Msn2 as a central factor allowing TF binding at its target promoters. By contrast, no case of mutual cooperation was observed. Particularly, Msn2 retained binding at its preferred promoters upon co-depletion of fourteen similarly bound TFs. Overall, the consequences of TF co-depletions were moderate, limited to a subset of promoters, and failed to explain the role of regions outside the DNA-binding domain in directing TF-binding preferences. Our results call for re-evaluating the role of cooperative interactions in directing TF-binding preferences.

A3DyDB: exploring structural aggregation propensities in the yeast proteome

BACKGROUND

The budding yeast Saccharomyces cerevisiae (S. cerevisiae) is a well-established model system for studying protein aggregation due to the conservation of essential cellular structures and pathways found across eukaryotes. However, limited structural knowledge of its proteome has prevented a deeper understanding of yeast functionalities, interactions, and aggregation.

RESULTS

In this study, we introduce the A3D yeast database (A3DyDB), which offers an extensive catalog of aggregation propensity predictions for the S. cerevisiae proteome. We used Aggrescan 3D (A3D) and the newly released protein models from AlphaFold2 (AF2) to compute the structure-based aggregation predictions for 6039 yeast proteins. The A3D algorithm exploits the information from 3D protein structures to calculate their intrinsic aggregation propensities. To facilitate simple and intuitive data analysis, A3DyDB provides a user-friendly interface for querying, browsing, and visualizing information on aggregation predictions from yeast protein structures. The A3DyDB also allows for the evaluation of the influence of natural or engineered mutations on protein stability and solubility. The A3DyDB is freely available at http://biocomp.chem.uw.edu.pl/A3D2/yeast .

CONCLUSION

The A3DyDB addresses a gap in yeast resources by facilitating the exploration of correlations between structural aggregation propensity and diverse protein properties at the proteome level. We anticipate that this comprehensive database will become a standard tool in the modeling of protein aggregation and its implications in budding yeast.

A dynamic actin cytoskeleton is required to prevent constitutive VDAC-dependent MAPK signalling and aberrant lipid homeostasis

The dynamic nature of the actin cytoskeleton is required to coordinate many cellular processes, and a loss of its plasticity has been linked to accelerated cell aging and attenuation of adaptive response mechanisms. Cofilin is an actin-binding protein that controls actin dynamics and has been linked to mitochondrial signaling pathways that control drug resistance and cell death. Here we show that cofilin-driven chronic depolarization of the actin cytoskeleton activates cell wall integrity mitogen-activated protein kinase (MAPK) signalling and disrupts lipid homeostasis in a voltage-dependent anion channel (VDAC)-dependent manner. Expression of the cof1-5 mutation, which reduces the dynamic nature of actin, triggers loss of cell wall integrity, vacuole fragmentation, disruption of lipid homeostasis, lipid droplet (LD) accumulation, and the promotion of cell death. The integrity of the actin cytoskeleton is therefore essential to maintain the fidelity of MAPK signaling, lipid homeostasis, and cell health in S. cerevisiae.

Nucleosome retention by histone chaperones and remodelers occludes pervasive DNA-protein binding

DNA packaging within chromatin depends on histone chaperones and remodelers that form and position nucleosomes. Cells express multiple such chromatin regulators with overlapping in-vitro activities. Defining specific in-vivo activities requires monitoring histone dynamics during regulator depletion, which has been technically challenging. We have recently generated histone-exchange sensors in Saccharomyces cerevisiae, which we now use to define the contributions of 15 regulators to histone dynamics genome-wide. While replication-independent exchange in unperturbed cells maps to promoters, regulator depletions primarily affected gene bodies. Depletion of Spt6, Spt16 or Chd1 sharply increased nucleosome replacement sequentially at the beginning, middle or end of highly expressed gene bodies. They further triggered re-localization of chaperones to affected gene body regions, which compensated for nucleosome loss during transcription complex passage, but concurred with extensive TF binding in gene bodies. We provide a unified quantitative screen highlighting regulator roles in retaining nucleosome binding during transcription and preserving genomic packaging.

C-terminal sequence stability profiling in Saccharomyces cerevisiae reveals protective protein quality control pathways

Protein quality control (PQC) mechanisms are essential for degradation of misfolded or dysfunctional proteins. An essential part of protein homeostasis is recognition of defective proteins by PQC components and their elimination by the ubiquitin-proteasome system, often concentrating on protein termini as indicators of protein integrity. Changes in amino acid composition of C-terminal ends arise through protein disintegration, alternative splicing, or during the translation step of protein synthesis from premature termination or translational stop-codon read-through. We characterized reporter protein stability using light-controlled exposure of the random C-terminal peptide collection (CtPC) in budding yeast revealing stabilizing and destabilizing features of amino acids at positions -5 to -1 of the C terminus. The (de)stabilization properties of CtPC-degrons depend on amino acid identity, position, as well as composition of the C-terminal sequence and are transferable. Evolutionary pressure toward stable proteins in yeast is evidenced by amino acid residues under-represented in cytosolic and nuclear proteins at corresponding C-terminal positions, but over-represented in unstable CtPC-degrons, and vice versa. Furthermore, analysis of translational stop-codon read-through peptides suggested that such extended proteins have destabilizing C termini. PQC pathways targeting CtPC-degrons involved the ubiquitin-protein ligase Doa10 and the cullin-RING E3 ligase SCFDas1 (Skp1-Cullin-F-box protein). Overall, our data suggest a proteome protection mechanism that targets proteins with unnatural C termini by recognizing a surprisingly large number of C-terminal sequence variants.

Inferring sparse structure in genotype-phenotype maps

Correlation among multiple phenotypes across related individuals may reflect some pattern of shared genetic architecture: individual genetic loci affect multiple phenotypes (an effect known as pleiotropy), creating observable relationships between phenotypes. A natural hypothesis is that pleiotropic effects reflect a relatively small set of common "core" cellular processes: each genetic locus affects one or a few core processes, and these core processes in turn determine the observed phenotypes. Here, we propose a method to infer such structure in genotype-phenotype data. Our approach, sparse structure discovery (SSD) is based on a penalized matrix decomposition designed to identify latent structure that is low-dimensional (many fewer core processes than phenotypes and genetic loci), locus-sparse (each locus affects few core processes), and/or phenotype-sparse (each phenotype is influenced by few core processes). Our use of sparsity as a guide in the matrix decomposition is motivated by the results of a novel empirical test indicating evidence of sparse structure in several recent genotype-phenotype datasets. First, we use synthetic data to show that our SSD approach can accurately recover core processes if each genetic locus affects few core processes or if each phenotype is affected by few core processes. Next, we apply the method to three datasets spanning adaptive mutations in yeast, genotoxin robustness assay in human cell lines, and genetic loci identified from a yeast cross, and evaluate the biological plausibility of the core process identified. More generally, we propose sparsity as a guiding prior for resolving latent structure in empirical genotype-phenotype maps.

Identification and characterization of protein interactions with the major Niemann-Pick type C disease protein in yeast reveals pathways of therapeutic potential

Niemann-Pick type C (NP-C) disease is a rare lysosomal storage disease caused by mutations in NPC1 (95% cases) or NPC2 (5% cases). These proteins function together in cholesterol egress from the lysosome, whereby upon mutation, cholesterol and other lipids accumulate causing major pathologies. However, it is not fully understood how cholesterol is transported from NPC1 residing at the lysosomal membrane to the endoplasmic reticulum (ER) and plasma membrane. The yeast ortholog of NPC1, Niemann-Pick type C-related protein-1 (Ncr1), functions similarly to NPC1; when transfected into a mammalian cell lacking NPC1, Ncr1 rescues the diagnostic hallmarks of cholesterol and sphingolipid accumulation. Here, we aimed to identify and characterize protein-protein interactions (PPIs) with the yeast Ncr1 protein. A genome-wide split-ubiquitin membrane yeast two-hybrid (MYTH) protein interaction screen identified 11 ER membrane-localized, full-length proteins interacting with Ncr1 at the lysosomal/vacuolar membrane. These highlight the importance of ER-vacuole membrane interface and include PPIs with the Cyb5/Cbr1 electron transfer system, the ceramide synthase complex, and the Sec61/Sbh1 protein translocation complex. These PPIs were not detected in a sterol auxotrophy condition and thus depend on normal sterol metabolism. To provide biological context for the Ncr1-Cyb5 PPI, a yeast strain lacking this PPI (via gene deletions) exhibited altered levels of sterols and sphingolipids including increased levels of glucosylceramide that mimic NP-C disease. Overall, the results herein provide new physical and genetic interaction models to further use the yeast model of NP-C disease to better understand human NP-C disease.

Fungal COP9 signalosome assembly requires connection of two trimeric intermediates for integration of intrinsic deneddylase

The conserved eight-subunit COP9 signalosome (CSN) is required for multicellular fungal development. The CSN deneddylase cooperates with the Cand1 exchange factor to control replacements of E3 ubiquitin cullin RING ligase receptors, providing specificity to eukaryotic protein degradation. Aspergillus nidulans CSN assembles through a heptameric pre-CSN, which is activated by integration of the catalytic CsnE deneddylase. Combined genetic and biochemical approaches provided the assembly choreography within a eukaryotic cell for native fungal CSN. Interactomes of functional GFP-Csn subunit fusions in pre-CSN deficient fungal strains were compared by affinity purifications and mass spectrometry. Two distinct heterotrimeric CSN subcomplexes were identified as pre-CSN assembly intermediates. CsnA-C-H and CsnD-F-G form independently of CsnB, which connects the heterotrimers to a heptamer and enables subsequent integration of CsnE to form the enzymatically active CSN complex. Surveillance mechanisms control accurate Csn subunit amounts and correct cellular localization for sequential assembly since deprivation of Csn subunits changes the abundance and location of remaining Csn subunits.

The genomic region of the 3' untranslated region (3'UTR) of PHO84, rather than the antisense RNA, promotes gene repression

PHO84 is a budding yeast gene reported to be negatively regulated by its cognate antisense transcripts both in cis and in trans. In this study, we performed Transient-transcriptome sequencing (TT-seq) to investigate the correlation of sense/antisense pairs in a dbp2Δ strain and found over 700 sense/antisense pairs, including PHO84, to be positively correlated, contrasting the prevailing model. To define what mechanism regulates the PHO84 gene and how this regulation could have been originally attributed to repression by the antisense transcript, we conducted a series of molecular biology and genetics experiments. We now report that the 3' untranslated region (3'UTR) of PHO84 plays a repressive role in sense expression, an activity not linked to the antisense transcripts. Moreover, we provide results of a genetic screen for 3'UTR-dependent repression of PHO84 and show that the vast majority of identified factors are linked to negative regulation. Finally, we show that the PHO84 promoter and terminator form gene loops which correlate with transcriptional repression, and that the RNA-binding protein, Tho1, increases this looping and the 3'UTR-dependent repression. Our results negate the current model for antisense non-coding transcripts of PHO84 and suggest that many of these transcripts are byproducts of open chromatin.

HmbC, a Protein of the HMG Family, Participates in the Regulation of Carotenoid Biosynthesis in Fusarium fujikuroi

In the fungus Fusarium fujikuroi, carotenoid production is up-regulated by light and down-regulated by the CarS RING finger protein, which modulates the mRNA levels of carotenoid pathway genes (car genes). To identify new potential regulators of car genes, we used a biotin-mediated pull-down procedure to detect proteins capable of binding to their promoters. We focused our attention on one of the proteins found in the screening, belonging to the High-Mobility Group (HMG) family that was named HmbC. The deletion of the hmbC gene resulted in increased carotenoid production due to higher mRNA levels of car biosynthetic genes. In addition, the deletion resulted in reduced carS mRNA levels, which could also explain the partial deregulation of the carotenoid pathway. The mutants exhibited other phenotypic traits, such as alterations in development under certain stress conditions, or reduced sensitivity to cell wall degrading enzymes, revealed by less efficient protoplast formation, indicating that HmbC is also involved in other cellular processes. In conclusion, we identified a protein of the HMG family that participates in the regulation of carotenoid biosynthesis. This is probably achieved through an epigenetic mechanism related to chromatin structure, as is frequent in this class of proteins.

CryptoCEN: A Co-Expression Network for Cryptococcus neoformans reveals novel proteins involved in DNA damage repair

Elucidating gene function is a major goal in biology, especially among non-model organisms. However, doing so is complicated by the fact that molecular conservation does not always mirror functional conservation, and that complex relationships among genes are responsible for encoding pathways and higher-order biological processes. Co-expression, a promising approach for predicting gene function, relies on the general principal that genes with similar expression patterns across multiple conditions will likely be involved in the same biological process. For Cryptococcus neoformans, a prevalent human fungal pathogen greatly diverged from model yeasts, approximately 60% of the predicted genes in the genome lack functional annotations. Here, we leveraged a large amount of publicly available transcriptomic data to generate a C. neoformans Co-Expression Network (CryptoCEN), successfully recapitulating known protein networks, predicting gene function, and enabling insights into the principles influencing co-expression. With 100% predictive accuracy, we used CryptoCEN to identify 13 new DNA damage response genes, underscoring the utility of guilt-by-association for determining gene function. Overall, co-expression is a powerful tool for uncovering gene function, and decreases the experimental tests needed to identify functions for currently under-annotated genes.

In vivo characterization of the critical interaction between the RNA exosome and the essential RNA helicase Mtr4 in Saccharomyces cerevisiae

The RNA exosome is a conserved molecular machine that processes/degrades numerous coding and non-coding RNAs. The 10-subunit complex is composed of three S1/KH cap subunits (human EXOSC2/3/1; yeast Rrp4/40/Csl4), a lower ring of six PH-like subunits (human EXOSC4/7/8/9/5/6; yeast Rrp41/42/43/45/46/Mtr3), and a singular 3'-5' exo/endonuclease DIS3/Rrp44. Recently, several disease-linked missense mutations have been identified in structural cap and core RNA exosome genes. In this study, we characterize a rare multiple myeloma patient missense mutation that was identified in the cap subunit gene EXOSC2. This missense mutation results in a single amino acid substitution, p.Met40Thr, in a highly conserved domain of EXOSC2. Structural studies suggest that this Met40 residue makes direct contact with the essential RNA helicase, MTR4, and may help stabilize the critical interaction between the RNA exosome complex and this cofactor. To assess this interaction in vivo, we utilized the Saccharomyces cerevisiae system and modeled the EXOSC2 patient mutation into the orthologous yeast gene RRP4, generating the variant rrp4-M68T. The rrp4-M68T cells show accumulation of certain RNA exosome target RNAs and show sensitivity to drugs that impact RNA processing. We also identified robust negative genetic interactions between rrp4-M68T and specific mtr4 mutants. A complementary biochemical approach revealed that Rrp4 M68T shows decreased interaction with Mtr4, consistent with these genetic results. This study suggests that the EXOSC2 mutation identified in a multiple myeloma patient impacts the function of the RNA exosome and provides functional insight into a critical interface between the RNA exosome and Mtr4.