Fission stories: using PomBase to understand Schizosaccharomyces pombe biology

Arrayed CRISPRi library to suppress genes required for Schizosaccharomyces pombe viability

The fission yeast, Schizosaccharomyces pombe, is an excellent eukaryote model organism for studying essential biological processes. Its genome contains ∼1,200 genes essential for cell viability, most of which are evolutionarily conserved. To study these essential genes, resources enabling conditional perturbation of target genes are required. Here, we constructed comprehensive arrayed libraries of plasmids and strains to knock down essential genes in S. pombe using dCas9-mediated CRISPRi. These libraries cover ∼98% of all essential genes in fission yeast. We estimate that in ∼60% of these strains, transcription of a target gene was repressed so efficiently that cell proliferation was significantly inhibited. To demonstrate the usefulness of these libraries, we performed metabolic analyses with knockdown strains and revealed flexible interaction among metabolic pathways. Libraries established in this study enable comprehensive functional analyses of essential genes in S. pombe and will facilitate the understanding of essential biological processes in eukaryotes.

The DNA Damage Repair Function of Fission Yeast CK1 Involves Targeting Arp8, a Subunit of the INO80 Chromatin Remodeling Complex

The CK1 family are conserved serine/threonine kinases with numerous substrates and cellular functions. The fission yeast CK1 orthologues Hhp1 and Hhp2 were first characterized as regulators of DNA repair, but the mechanism(s) by which CK1 activity promotes DNA repair had not been investigated. Here, we found that deleting Hhp1 and Hhp2 or inhibiting CK1 catalytic activities in yeast or in human cells increased double-strand breaks (DSBs). The primary pathways to repair DSBs, homologous recombination and nonhomologous end joining, were both less efficient in cells lacking Hhp1 and Hhp2 activity. To understand how Hhp1 and Hhp2 promote DNA damage repair, we identified new substrates of these enzymes using quantitative phosphoproteomics. We confirmed that Arp8, a component of the INO80 chromatin remodeling complex, is a bona fide substrate of Hhp1 and Hhp2 important for DNA repair. Our data suggest that Hhp1 and Hhp2 facilitate DNA repair by phosphorylating multiple substrates, including Arp8.

Ageing-associated long non-coding RNA extends lifespan and reduces translation in non-dividing cells

Genomes produce widespread long non-coding RNAs (lncRNAs) of largely unknown functions. We characterize aal1 (ageing-associated lncRNA), which is induced in quiescent fission yeast cells. Deletion of aal1 shortens the chronological lifespan of non-dividing cells, while ectopic overexpression prolongs their lifespan, indicating that aal1 acts in trans. Overexpression of aal1 represses ribosomal-protein gene expression and inhibits cell growth, and aal1 genetically interacts with coding genes functioning in protein translation. The aal1 lncRNA localizes to the cytoplasm and associates with ribosomes. Notably, aal1 overexpression decreases the cellular ribosome content and inhibits protein translation. The aal1 lncRNA binds to the rpl1901 mRNA, encoding a ribosomal protein. The rpl1901 levels are reduced ~2-fold by aal1, which is sufficient to extend lifespan. Remarkably, the expression of the aal1 lncRNA in Drosophila boosts fly lifespan. We propose that aal1 reduces the ribosome content by decreasing Rpl1901 levels, thus attenuating the translational capacity and promoting longevity. Although aal1 is not conserved, its effect in flies suggests that animals feature related mechanisms that modulate ageing, based on the conserved translational machinery.

Fission yeast Duc1 links to ER-PM contact sites and influences PM lipid composition and cytokinetic ring anchoring

Cytokinesis is the final stage of the cell cycle that results in the physical separation of daughter cells. To accomplish cytokinesis, many organisms build an actin- and myosin-based cytokinetic ring (CR) that is anchored to the plasma membrane (PM). Defects in CR-PM anchoring can arise when the PM lipid phosphatidylinositol (4,5)-bisphosphate [PI(4,5)P2] is depleted. In Schizosaccharomyces pombe, reduced PM PI(4,5)P2 results in a CR that cannot maintain a medial position and slides toward one cell end, resulting in two differently sized daughter cells. S. pombe PM PI(4,5)P2 is synthesized by the phosphatidylinositol 4-phosphate 5-kinase (PI5-kinase) Its3, but what regulates this enzyme to maintain appropriate PM PI(4,5)P2 levels in S. pombe is not known. To identify Its3 regulators, we used proximity-based biotinylation, and the uncharacterized protein Duc1 was specifically detected. We discovered that Duc1 decorates the PM except at the cell division site and that its unique localization pattern is dictated by binding to the endoplasmic reticulum (ER)-PM contact site proteins Scs2 and Scs22. Our evidence suggests that Duc1 also binds PI(4,5)P2 and helps enrich Its3 at the lateral PM, thereby promoting PM PI(4,5)P2 synthesis and robust CR-PM anchoring.

Structure-guided Evolutionary Analysis of Interactome Network Rewiring at Single Residue Resolution in Yeasts

Protein-protein interactions (PPIs) are known to rewire extensively during evolution leading to lineage-specific and species-specific changes in molecular processes. However, the detailed molecular evolutionary mechanisms underlying interactome network rewiring are not well-understood. Here, we combine high-confidence PPI data, high-resolution three-dimensional structures of protein complexes, and homology-based structural annotation transfer to construct structurally-resolved interactome networks for the two yeasts S. cerevisiae and S. pombe. We then classify PPIs according to whether they are preserved or different between the two yeast species and compare site-specific evolutionary rates of interfacial versus non-interfacial residues for these different categories of PPIs. We find that residues in PPI interfaces evolve significantly more slowly than non-interfacial residues when using lineage-specific measures of evolutionary rate, but not when using non-lineage-specific measures. Furthermore, both lineage-specific and non-lineage-specific evolutionary rate measures can distinguish interfacial residues from non-interfacial residues for preserved PPIs between the two yeasts, but only the lineage-specific measure is appropriate for rewired PPIs. Finally, both lineage-specific and non-lineage-specific evolutionary rate measures are appropriate for elucidating structural determinants of protein evolution for residues outside of PPI interfaces. Overall, our results demonstrate that unlike tertiary structures of single proteins, PPIs and PPI interfaces can be highly volatile in their evolution, thus requiring the use of lineage-specific measures when studying their evolution. These results yield insight into the evolutionary design principles of PPIs and the mechanisms by which interactions are preserved or rewired between species, improving our understanding of the molecular evolution of PPIs and PPI interfaces at the residue level.

DEAD-box ATPase Dbp2 is the key enzyme in an mRNP assembly checkpoint at the 3'-end of genes and involved in the recycling of cleavage factors

mRNA biogenesis in the eukaryotic nucleus is a highly complex process. The numerous RNA processing steps are tightly coordinated to ensure that only fully processed transcripts are released from chromatin for export from the nucleus. Here, we present the hypothesis that fission yeast Dbp2, a ribonucleoprotein complex (RNP) remodelling ATPase of the DEAD-box family, is the key enzyme in an RNP assembly checkpoint at the 3'-end of genes. We show that Dbp2 interacts with the cleavage and polyadenylation complex (CPAC) and localises to cleavage bodies, which are enriched for 3'-end processing factors and proteins involved in nuclear RNA surveillance. Upon loss of Dbp2, 3'-processed, polyadenylated RNAs accumulate on chromatin and in cleavage bodies, and CPAC components are depleted from the soluble pool. Under these conditions, cells display an increased likelihood to skip polyadenylation sites and a delayed transcription termination, suggesting that levels of free CPAC components are insufficient to maintain normal levels of 3'-end processing. Our data support a model in which Dbp2 is the active component of an mRNP remodelling checkpoint that licenses RNA export and is coupled to CPAC release.

Understanding the molecular mechanisms of human diseases: the benefits of fission yeasts

The role of model organisms such as yeasts in life science research is crucial. Although the baker's yeast (Saccharomyces cerevisiae) is the most popular model among yeasts, the contribution of the fission yeasts (Schizosaccharomyces) to life science is also indisputable. Since both types of yeasts share several thousands of common orthologous genes with humans, they provide a simple research platform to investigate many fundamental molecular mechanisms and functions, thereby contributing to the understanding of the background of human diseases. In this review, we would like to highlight the many advantages of fission yeasts over budding yeasts. The usefulness of fission yeasts in virus research is shown as an example, presenting the most important research results related to the Human Immunodeficiency Virus Type 1 (HIV-1) Vpr protein. Besides, the potential role of fission yeasts in the study of prion biology is also discussed. Furthermore, we are keen to promote the uprising model yeast Schizosaccharomyces japonicus, which is a dimorphic species in the fission yeast genus. We propose the hyphal growth of S. japonicus as an unusual opportunity as a model to study the invadopodia of human cancer cells since the two seemingly different cell types can be compared along fundamental features. Here we also collect the latest laboratory protocols and bioinformatics tools for the fission yeasts to highlight the many possibilities available to the research community. In addition, we present several limiting factors that everyone should be aware of when working with yeast models.

Schizosaccharomyces pombe as a fundamental model for research on mitochondrial gene expression: Progress, achievements and outlooks

Schizosaccharomyces pombe (fission yeast) is an attractive model for mitochondrial research. The organism resembles human cells in terms of mitochondrial inheritance, mitochondrial transport, sugar metabolism, mitogenome structure and dependence of viability on the mitogenome (the petite-negative phenotype). Transcriptions of these genomes produce only a few polycistronic transcripts, which then undergo processing as per the tRNA punctuation model. In general, the machinery for mitochondrial gene expression is structurally and functionally conserved between fission yeast and humans. Furthermore, molecular research on S. pombe is supported by a considerable number of experimental techniques and database resources. Owing to these advantages, fission yeast has significantly contributed to biomedical and fundamental research. Here, we review the current state of knowledge regarding S. pombe mitochondrial gene expression, and emphasise the pertinence of fission yeast as both a model and tool, especially for studies on mitochondrial translation.

Electrophoretic mobility shift assays (EMSAs) for in vitro detection of protein-nucleic acid interactions

Protein-nucleic acid interactions drive some of the most important physiological events in cells. Here, we present a protocol for detecting protein-DNA or protein-RNA interactions in vitro. We describe steps for labeling nucleic acid species and electrophoretic mobility shift assays (EMSAs). This protocol can be used to confirm suspected in vivo interactions using recombinantly expressed/purified proteins of interest and a nucleic acid substrate. It can further be used to investigate mutations that can disrupt interaction or compensatory mutations that restore it. For complete details on the use and execution of this protocol, please refer to Mansouri-Noori et al.1.

In Vivo Proximity Labeling Identifies a New Function for the Lifespan and Autophagy-regulating Kinase Pef1, an Ortholog of Human Cdk5

Cdk5 is a highly-conserved, noncanonical cell division kinase important to the terminal differentiation of mammalian cells in multiple organ systems. We previously identified Pef1, the Schizosaccharomyces pombe ortholog of cdk5, as regulator of chronological lifespan. To reveal the processes impacted by Pef1, we developed APEX2-biotin phenol-mediated proximity labeling in S. pombe. Efficient labeling required a short period of cell wall digestion and eliminating glucose and nitrogen sources from the medium. We identified 255 high-confidence Pef1 neighbors in growing cells and a novel Pef1-interacting partner, the DNA damage response protein Rad24. The Pef1-Rad24 interaction was validated by reciprocal proximity labeling and co-immunoprecipitation. Eliminating Pef1 partially rescued the DNA damage sensitivity of cells lacking Rad24. To monitor how Pef1 neighbors change under different conditions, cells induced for autophagy were labeled and 177 high-confidence Pef1 neighbors were identified. Gene ontology (GO) analysis of the Pef1 neighbors identified proteins participating in processes required for autophagosome expansion including regulation of actin dynamics and vesicle-mediated transport. Some of these proteins were identified in both exponentially growing and autophagic cells. Pef1-APEX2 proximity labeling therefore identified a new Pef1 function in modulating the DNA damage response and candidate processes that Pef1 and other cdk5 orthologs may regulate.

Laboratory horror stories: Poison in the agars

The fission yeast Schizosaccharomyces pombe is a single-celled eukaryote that can be cultured as a haploid or as a diploid. Scientists employ mating, meiosis, and the plating of ascospores and cells to generate strains with novel genotypes and to discover biological processes. Our two laboratories encountered independently sudden-onset, major impediments to such research. Spore suspensions and vegetative cells no longer plated effectively on minimal media. By systematically analyzing multiple different media components from multiple different suppliers, we identified the source of the problem. Specific lots of agar, from different suppliers, were toxic. Interestingly, the inhibitory effect was attenuated on rich media. Consequently, quality control checks that use only rich media can provide false assurances on the quality of the agar. Lastly, we describe likely sources of the toxicity and we provide specific guidance for quality control measures that should be applied by all vendors as preconditions for their sale of agar.

Fungal microtubule organizing centers are evolutionarily unstable structures

For most Eukaryotic species the requirements of cilia formation dictate the structure of microtubule organizing centers (MTOCs). In this study we find that loss of cilia corresponds to loss of evolutionary stability for fungal MTOCs. We used iterative search algorithms to identify proteins homologous to those found in Saccharomyces cerevisiae, and Schizosaccharomyces pombe MTOCs, and calculated site-specific rates of change for those proteins that were broadly phylogenetically distributed. Our results indicate that both the protein composition of MTOCs as well as the sequence of MTOC proteins are poorly conserved throughout the fungal kingdom. To begin to reconcile this rapid evolutionary change with the rigid structure and essential function of the S. cerevisiae MTOC we further analyzed how structural interfaces among proteins influence the rates of change for specific residues within a protein. We find that a more stable protein may stabilize portions of an interacting partner where the two proteins are in contact. In summary, while the protein composition and sequences of the MTOC may be rapidly changing the proteins within the structure have a stabilizing effect on one another. Further exploration of fungal MTOCs will expand our understanding of how changes in the functional needs of a cell have affected physical structures, proteomes, and protein sequences throughout fungal evolution.

In Silico Analysis of Protein-Protein Interactions of Putative Endoplasmic Reticulum Metallopeptidase 1 in Schizosaccharomyces pombe

Ermp1 is a putative metalloprotease from Schizosaccharomyces pombe and a member of the Fxna peptidases. Although their function is unknown, orthologous proteins from rats and humans have been associated with the maturation of ovarian follicles and increased ER stress. This study focuses on proposing the first prediction of PPI by comparison of the interologues between humans and yeasts, as well as the molecular docking and dynamics of the M28 domain of Ermp1 with possible target proteins. As results, 45 proteins are proposed that could interact with the metalloprotease. Most of these proteins are related to the transport of Ca2+ and the metabolism of amino acids and proteins. Docking and molecular dynamics suggest that the M28 domain of Ermp1 could hydrolyze leucine and methionine residues of Amk2, Ypt5 and Pex12. These results could support future experimental investigations of other Fxna peptidases, such as human ERMP1.

Substrate displacement of CK1 C-termini regulates kinase specificity

CK1 kinases participate in many signaling pathways, and their regulation is of meaningful biological consequence. CK1s autophosphorylate their C-terminal noncatalytic tails, and eliminating these tails increases substrate phosphorylation in vitro, suggesting that the autophosphorylated C-termini act as inhibitory pseudosubstrates. To test this prediction, we comprehensively identified the autophosphorylation sites on Schizosaccharomyces pombe Hhp1 and human CK1ε. Phosphoablating mutations increased Hhp1 and CK1ε activity toward substrates. Peptides corresponding to the C-termini interacted with the kinase domains only when phosphorylated, and substrates competitively inhibited binding of the autophosphorylated tails to the substrate binding grooves. Tail autophosphorylation influenced the catalytic efficiency with which CK1s targeted different substrates, and truncating the tail of CK1δ broadened its linear peptide substrate motif, indicating that tails contribute to substrate specificity as well. Considering autophosphorylation of both T220 in the catalytic domain and C-terminal sites, we propose a displacement specificity model to describe how autophosphorylation modulates substrate specificity for the CK1 family.

PomBase: a Global Core Biodata Resource-growth, collaboration, and sustainability

PomBase (https://www.pombase.org), the model organism database (MOD) for fission yeast, was recently awarded Global Core Biodata Resource (GCBR) status by the Global Biodata Coalition (GBC; https://globalbiodata.org/) after a rigorous selection process. In this MOD review, we present PomBase's continuing growth and improvement over the last 2 years. We describe these improvements in the context of the qualitative GCBR indicators related to scientific quality, comprehensivity, accelerating science, user stories, and collaborations with other biodata resources. This review also showcases the depth of existing connections both within the biocuration ecosystem and between PomBase and its user community.

Regulation of sexual differentiation initiation in Schizosaccharomyces pombe

The fission yeast Schizosaccharomyces pombe is an excellent model organism to explore cellular events owing to rich tools in genetics, molecular biology, cellular biology, and biochemistry. Schizosaccharomyces pombe proliferates continuously when nutrients are abundant but arrests in G1 phase upon depletion of nutrients such as nitrogen and glucose. When cells of opposite mating types are present, cells conjugate, fuse, undergo meiosis, and finally form 4 spores. This sexual differentiation process in S. pombe has been studied extensively. To execute sexual differentiation, the glucose-sensing cAMP-PKA (cyclic adenosine monophosphate-protein kinase A) pathway, nitrogen-sensing TOR (target of rapamycin) pathway, and SAPK (stress-activating protein kinase) pathway are crucial, and the MAPK (mitogen-activating protein kinase) cascade is essential for pheromone sensing. These signals regulate ste11 at the transcriptional and translational levels, and Ste11 is modified in multiple ways. This review summarizes the initiation of sexual differentiation in S. pombe based on results I have helped to obtain, including the work of many excellent researchers.

The kinase Rio1 and a ribosome collision-dependent decay pathway survey the integrity of 18S rRNA cleavage

The 18S rRNA sequence is highly conserved, particularly at its 3'-end, which is formed by the endonuclease Nob1. How Nob1 identifies its target sequence is not known, and in vitro experiments have shown Nob1 to be error-prone. Moreover, the sequence around the 3'-end is degenerate with similar sites nearby. Here, we used yeast genetics, biochemistry, and next-generation sequencing to investigate a role for the ATPase Rio1 in monitoring the accuracy of the 18S rRNA 3'-end. We demonstrate that Nob1 can miscleave its rRNA substrate and that miscleaved rRNA accumulates upon bypassing the Rio1-mediated quality control (QC) step, but not in healthy cells with intact QC mechanisms. Mechanistically, we show that Rio1 binding to miscleaved rRNA is weaker than its binding to accurately processed 18S rRNA. Accordingly, excess Rio1 results in accumulation of miscleaved rRNA. Ribosomes containing miscleaved rRNA can translate, albeit more slowly, thereby inviting collisions with trailing ribosomes. These collisions result in degradation of the defective ribosomes utilizing parts of the machinery for mRNA QC. Altogether, the data support a model in which Rio1 inspects the 3'-end of the nascent 18S rRNA to prevent miscleaved 18S rRNA-containing ribosomes from erroneously engaging in translation, where they induce ribosome collisions. The data also demonstrate how ribosome collisions purify cells of altered ribosomes with different functionalities, with important implications for the concept of ribosome heterogeneity.

Generation and characterization of temperature-sensitive alleles of the glucanosyltransferase Gas1 in Schizosaccharomyces pombe

The Schizosaccharomyces pombe Gas family of β-1,3-glucanosyltransferases modify the cell wall by elongating β-1,3-glucan chains. While gas1Δ cells are inviable under standard laboratory growth conditions, they are viable in the presence of an osmotic stabilizer. Even under these conditions however, gas1Δ cells are slow-growing and display cell separation and morphology defects. Here, we isolated and characterized two gas1 temperature-sensitive alleles. Our data support that Gas1 is the primary S. pombe β-1,3-glucanosyltransferase important for cell separation and cell viability and provide useful tools for further analysis of S. pombe cell wall formation.

Schizosaccharomyces versatilis represents a distinct evolutionary lineage of fission yeast

The fission yeast species Schizosaccharomyces japonicus is currently divided into two varieties-S. japonicus var. japonicus and S. japonicus var. versatilis. Here we examine the var. versatilis isolate CBS5679. The CBS5679 genome shows 88% identity to the reference genome of S. japonicus var. japonicus at the coding sequence level, with phylogenetic analyses suggesting that it has split from the S. japonicus lineage 25 million years ago. The CBS5679 genome contains a reciprocal translocation between chromosomes 1 and 2, together with several large inversions. The products of genes linked to the major translocation are associated with 'metabolism' and 'cellular assembly' ontology terms. We further show that CBS5679 does not generate viable progeny with the reference strain of S. japonicus. Although CBS5679 shares closer similarity to the 'type' strain of var. versatilis as compared to S. japonicus, it is not identical to the type strain, suggesting population structure within var. versatilis. We recommend that the taxonomic status of S. japonicus var. versatilis is raised, with it being treated as a separate species, Schizosaccharomyces versatilis.

A clustering procedure for three-way RNA sequencing data using data transformations and matrix-variate Gaussian mixture models

RNA sequencing of time-course experiments results in three-way count data where the dimensions are the genes, the time points and the biological units. Clustering RNA-seq data allows to extract groups of co-expressed genes over time. After standardisation, the normalised counts of individual genes across time points and biological units have similar properties as compositional data. We propose the following procedure to suitably cluster three-way RNA-seq data: (1) pre-process the RNA-seq data by calculating the normalised expression profiles, (2) transform the data using the additive log ratio transform to map the composition in the D-part Aitchison simplex to a D - 1 -dimensional Euclidean vector, (3) cluster the transformed RNA-seq data using matrix-variate Gaussian mixture models and (4) assess the quality of the overall cluster solution and of individual clusters based on cluster separation in the transformed space using density-based silhouette information and on compactness of the cluster in the original space using cluster maps as a suitable visualisation. The proposed procedure is illustrated on RNA-seq data from fission yeast and results are also compared to an analogous two-way approach after flattening out the biological units.

POMBOX: A Fission Yeast Cloning Toolkit for Molecular and Synthetic Biology

The fission yeast Schizosaccharomyces pombe is a popular model organism in molecular biology and cell physiology. With its ease of genetic manipulation and growth, supported by in-depth functional annotations in the PomBase database and genome-wide metabolic models,S. pombe is an attractive option for synthetic biology applications. However,S. pombe currently lacks modular tools for generating genetic circuits with more than 1 transcriptional unit. We developed a toolkit to address this gap. Adapted from the MoClo-YTK plasmid kit for Saccharomyces cerevisiae and using the same modular cloning grammar, our POMBOX toolkit is designed to facilitate fast, efficient, and modular construction of genetic circuits inS. pombe. It allows for interoperability when working with DNA sequences that are functional in bothS. cerevisiae and S. pombe (e.g., protein tags, antibiotic resistance cassettes, and coding sequences). Moreover, POMBOX enables the modular assembly of multigene pathways and increases the possible pathway length from 6 to 12 transcriptional units. We also adapted the stable integration vector homology arms to Golden Gate assembly and tested the genomic integration success rates depending on different sequence sizes, from 4 to 24 kb. We included 14 S. pombe promoters that we characterized using two fluorescent proteins, in both minimally defined (EMM2─Edinburgh minimal media) and complex (YES─yeast extract with supplements) media. Then, we examined the efficacy of 6 S. cerevisiae and 6 synthetic terminators in S. pombe. Finally, we used the POMBOX kit for a synthetic biology application in metabolic engineering and expressed plant enzymes in S. pombe to produce specialized metabolite precursors, namely, methylxanthine, amorpha-4,11-diene, and cinnamic acid from the purine, mevalonate, and aromatic amino acid pathways.

The Cdc14 phosphatase, Clp1, does not affect genome expression

Schizosaccharomyces pombe Clp1 is a Cdc14-family phosphatase that reverses mitotic Cdk1 phosphorylation. Despite evolutionary conservation, Clp1 's mammalian orthologs do not share this function. Rather, higher eukaryotic Cdc14 enzymes act in DNA repair, ciliogenesis, and gene regulation. To examine if Clp1 regulates gene expression, we compared the transcriptional profiles of cells lacking Clp1 function to that of wildtype. Because clp1∆ cells are sensitive to the actin depolymerizing drug, LatrunculinA, we also investigated whether a transcriptional response was involved. Our results indicate that Clp1 does not detectably affect gene expression and highlight the organism-specific functions of this conserved phosphatase family.

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.

A3D Model Organism Database (A3D-MODB): a database for proteome aggregation predictions in model organisms

Protein aggregation has been associated with aging and different pathologies and represents a bottleneck in the industrial production of biotherapeutics. Numerous past studies performed in Escherichia coli and other model organisms have allowed to dissect the biophysical principles underlying this process. This knowledge fuelled the development of computational tools, such as Aggrescan 3D (A3D) to forecast and re-design protein aggregation. Here, we present the A3D Model Organism Database (A3D-MODB) http://biocomp.chem.uw.edu.pl/A3D2/MODB, a comprehensive resource for the study of structural protein aggregation in the proteomes of 12 key model species spanning distant biological clades. In addition to A3D predictions, this resource incorporates information useful for contextualizing protein aggregation, including membrane protein topology and structural model confidence, as an indirect reporter of protein disorder. The database is openly accessible without any need for registration. We foresee A3D-MOBD evolving into a central hub for conducting comprehensive, multi-species analyses of protein aggregation, fostering the development of protein-based solutions for medical, biotechnological, agricultural and industrial applications.

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.

Fission Yeast TORC1 Promotes Cell Proliferation through Sfp1, a Transcription Factor Involved in Ribosome Biogenesis

Target of rapamycin complex 1 (TORC1) is activated in response to nutrient availability and growth factors, promoting cellular anabolism and proliferation. To explore the mechanism of TORC1-mediated proliferation control, we performed a genetic screen in fission yeast and identified Sfp1, a zinc-finger transcription factor, as a multicopy suppressor of temperature-sensitive TORC1 mutants. Our observations suggest that TORC1 phosphorylates Sfp1 and protects Sfp1 from proteasomal degradation. Transcription analysis revealed that Sfp1 positively regulates genes involved in ribosome production together with two additional transcription factors, Ifh1/Crf1 and Fhl1. Ifh1 physically interacts with Fhl1, and the nuclear localization of Ifh1 is regulated in response to nutrient levels in a manner dependent on TORC1 and Sfp1. Taken together, our data suggest that the transcriptional regulation of the genes involved in ribosome biosynthesis by Sfp1, Ifh1, and Fhl1 is one of the key pathways through which nutrient-activated TORC1 promotes cell proliferation.

Functional Prokaryotic-Like Deoxycytidine Triphosphate Deaminases and Thymidylate Synthase in Eukaryotic Social Amoebae: Vertical, Endosymbiotic, or Horizontal Gene Transfer?

The de novo synthesis of deoxythymidine triphosphate uses several pathways: gram-negative bacteria use deoxycytidine triphosphate deaminase to convert deoxycytidine triphosphate into deoxyuridine triphosphate, whereas eukaryotes and gram-positive bacteria instead use deoxycytidine monophosphate deaminase to transform deoxycytidine monophosphate to deoxyuridine monophosphate. It is then unusual that in addition to deoxycytidine monophosphate deaminases, the eukaryote Dictyostelium discoideum has 2 deoxycytidine triphosphate deaminases (Dcd1Dicty and Dcd2Dicty). Expression of either DcdDicty can fully rescue the slow growth of an Escherichia coli dcd knockout. Both DcdDicty mitigate the hydroxyurea sensitivity of a Schizosaccharomyces pombe deoxycytidine monophosphate deaminase knockout. Phylogenies show that Dcd1Dicty homologs may have entered the common ancestor of the eukaryotic groups of Amoebozoa, Obazoa, Metamonada, and Discoba through an ancient horizontal gene transfer from a prokaryote or an ancient endosymbiotic gene transfer from a mitochondrion, followed by horizontal gene transfer from Amoebozoa to several other unrelated groups of eukaryotes. In contrast, the Dcd2Dicty homologs were a separate horizontal gene transfer from a prokaryote or a virus into either Amoebozoa or Rhizaria, followed by a horizontal gene transfer between them. ThyXDicty, the D. discoideum thymidylate synthase, another enzyme of the deoxythymidine triphosphate biosynthesis pathway, was suggested previously to be acquired from the ancestral mitochondria or by horizontal gene transfer from alpha-proteobacteria. ThyXDicty can fully rescue the E. coli thymidylate synthase knockout, and we establish that it was obtained by the common ancestor of social amoebae not from mitochondria but from a bacterium. We propose horizontal gene transfer and endosymbiotic gene transfer contributed to the enzyme diversity of the deoxythymidine triphosphate synthesis pathway in most social amoebae, many Amoebozoa, and other eukaryotes.

Revised fission yeast gene and allele nomenclature guidelines for machine readability

Standardized nomenclature for genes, gene products, and isoforms is crucial to prevent ambiguity and enable clear communication of scientific data, facilitating efficient biocuration and data sharing. Standardized genotype nomenclature, which describes alleles present in a specific strain that differ from those in the wild-type reference strain, is equally essential to maximize research impact and ensure that results linking genotypes to phenotypes are Findable, Accessible, Interoperable, and Reusable (FAIR). In this publication, we extend the fission yeast clade gene nomenclature guidelines to support the curation efforts at PomBase (www.pombase.org), the Schizosaccharomyces pombe Model Organism Database. This update introduces nomenclature guidelines for noncoding RNA genes, following those set forth by the Human Genome Organisation Gene Nomenclature Committee. Additionally, we provide a significant update to the allele and genotype nomenclature guidelines originally published in 1987, to standardize the diverse range of genetic modifications enabled by the fission yeast genetic toolbox. These updated guidelines reflect a community consensus between numerous fission yeast researchers. Adoption of these rules will improve consistency in gene and genotype nomenclature, and facilitate machine-readability and automated entity recognition of fission yeast genes and alleles in publications or datasets. In conclusion, our updated guidelines provide a valuable resource for the fission yeast research community, promoting consistency, clarity, and FAIRness in genetic data sharing and interpretation.

Characterization of Pik1 function in fission yeast reveals its conserved role in lipid synthesis and not cytokinesis

Phosphatidylinositol (PI)-4-phosphate (PI4P) is a lipid found at the plasma membrane (PM) and Golgi in cells from yeast to humans. PI4P is generated from PI by PI4-kinases and can be converted into PI-4,5-bisphosphate [PI(4,5)P2]. Schizosaccharomyces pombe have two essential PI4-kinases - Stt4 and Pik1. Stt4 localizes to the PM, and its loss from the PM results in a decrease of PM PI4P and PI(4,5)P2. As a result, cells divide non-medially due to disrupted cytokinetic ring-PM anchoring. However, the localization and function of S. pombe Pik1 has not been thoroughly examined. Here, we found that Pik1 localizes exclusively to the trans-Golgi and is required for Golgi PI4P production. We determined that Ncs1 regulates Pik1, but unlike in other organisms, it is not required for Pik1 Golgi localization. When Pik1 function was disrupted, PM PI4P but not PI(4,5)P2 levels were reduced, a major difference compared with Stt4. We conclude that Stt4 is the chief enzyme responsible for producing the PI4P that generates PI(4,5)P2. Also, that cells with disrupted Pik1 do not divide asymmetrically highlights the specific importance of PM PI(4,5)P2 for cytokinetic ring-PM anchoring.

The LARP1 homolog Slr1p controls the stability and expression of proto-5'TOP mRNAs in fission yeast

Messenger RNAs (mRNAs) in higher eukaryotes that encode proteins important for the assembly of the translational apparatus (e.g., ribosomal proteins) often harbor a pyrimidine-rich motif at the extreme 5' end known as a 5' terminal oligopyrimidine (5'TOP) sequence. Members of the La-related protein 1 (LARP1) family control 5'TOP expression through a conserved DM15 motif, but the mechanism is not well understood. 5'TOP motifs have not been described in many lower organisms, and fission yeast harbors a LARP1 homolog that also lacks a DM15 motif. In this work, we show that the fission yeast LARP1 homolog, Slr1p, controls the translation and stability of mRNAs encoding proteins analogous to 5'TOP mRNAs in higher eukaryotes, which we thus refer to as proto-5'TOPs. Our data suggest that the LARP1 DM15 motif and the mRNA 5'TOP motif may be features that were scaffolded over a more fundamental mechanism of LARP1-associated control of gene expression.

Delineating yeast cleavage and polyadenylation signals using deep learning

3'-end cleavage and polyadenylation is an essential process for eukaryotic mRNA maturation. In yeast species, the polyadenylation signals that recruit the processing machinery are degenerate and remain poorly characterized compared to well-defined regulatory elements in mammals. Especially, recent deep sequencing experiments showed extensive cleavage heterogeneity for some mRNAs in Saccharomyces cerevisiae and uncovered the polyA motif differences between S. cerevisiae vs. Schizosaccharomyces pombe . The findings raised the fundamental question of how polyadenylation signals are formed in yeast. Here we addressed this question by developing deep learning models to deconvolute degenerate cis -regulatory elements and quantify their positional importance in mediating yeast polyA site formation, cleavage heterogeneity, and strength. In S. cerevisiae , cleavage heterogeneity is promoted by the depletion of U-rich elements around polyA sites as well as multiple occurrences of upstream UA-rich elements. Sites with high cleavage heterogeneity show overall lower strength. The site strength and tandem site distances modulate alternative polyadenylation (APA) under the diauxic stress. Finally, we developed a deep learning model to reveal the distinct motif configuration of S. pombe polyA sites which show more precise cleavage than S. cerevisiae . Altogether, our deep learning models provide unprecedented insights into polyA site formation across yeast species.

An ESCRT grommet cooperates with a diffusion barrier to maintain nuclear integrity

The molecular mechanisms by which the endosomal sorting complexes required for transport (ESCRT) proteins contribute to the integrity of the nuclear envelope (NE) barrier are not fully defined. We leveraged the single NE hole generated by mitotic extrusion of the Schizosaccharomyces pombe spindle pole body to reveal two modes of ESCRT function executed by distinct complements of ESCRT-III proteins, both dependent on CHMP7/Cmp7. A grommet-like function is required to restrict the NE hole in anaphase B, whereas replacement of Cmp7 by a sealing module ultimately closes the NE in interphase. Without Cmp7, nucleocytoplasmic compartmentalization remains intact despite NE discontinuities of up to 540 nm, suggesting mechanisms to limit diffusion through these holes. We implicate spindle pole body proteins as key components of a diffusion barrier acting with Cmp7 in anaphase B. Thus, NE remodelling mechanisms cooperate with proteinaceous diffusion barriers beyond nuclear pore complexes to maintain the nuclear compartment.

Elevated levels of sphingolipid MIPC in the plasma membrane disrupt the coordination of cell growth with cell wall formation in fission yeast

Coupling cell wall expansion with cell growth is a universal challenge faced by walled organisms. Mutations in Schizosaccharomyces pombe css1, which encodes a PM inositol phosphosphingolipid phospholipase C, prevent cell wall expansion but not synthesis of cell wall material. To probe how Css1 modulates cell wall formation we used classical and chemical genetics coupled with quantitative mass spectrometry. We found that elevated levels of the sphingolipid biosynthetic pathway's final product, mannosylinositol phosphorylceramide (MIPC), specifically correlated with the css1-3 phenotype. We also found that an apparent indicator of sphingolipids and a sterol biosensor accumulated at the cytosolic face of the PM at cell tips and the division site of css1-3 cells and, in accord, the PM in css1-3 was less dynamic than in wildtype cells. Interestingly, disrupting the protein glycosylation machinery recapitulated the css1-3 phenotype and led us to investigate Ghs2, a glycosylated PM protein predicted to modify cell wall material. Disrupting Ghs2 function led to aberrant cell wall material accumulation suggesting Ghs2 is dysfunctional in css1-3. We conclude that preventing an excess of MIPC in the S. pombe PM is critical to the function of key PM-localized proteins necessary for coupling growth with cell wall formation.

Meiotic recombination is confirmed to be unusually high in the fission yeast Schizosaccharomyces pombe

In most eukaryotes, meiotic crossovers (COs) are limited to 1-3 per chromosome, and are prevented from occurring close to one another by CO interference. The fission yeast Schizosaccharomyces pombe, an exception to these general rules, was reported to have the highest CO number per chromosome and no or weak interference. However, global CO frequency was indirectly estimated, calling for confirmation. Here, we used an innovative strategy to determine COs genome-wide in S. pombe. We confirmed weak CO interference, acting at physical distances compatible with the patterning of recombination precursors. We revealed a slight co-variation in CO number between chromosomes, suggesting that a limiting pro-CO factor varies between meiocytes. CO number per chromosome varies proportionally with chromosome size, with the three chromosomes having, on average, 15.9, 12.5, and 7.0 COs, respectively. This reinforces S. pombe's status as the eukaryote with the highest CO number per chromosome described to date.

The minimal intrinsic stochasticity of constitutively expressed eukaryotic genes is sub-Poissonian

Gene expression inherently gives rise to stochastic variation ("noise") in the production of gene products. Minimizing noise is crucial for ensuring reliable cellular functions. However, noise cannot be suppressed below a certain intrinsic limit. For constitutively expressed genes, this limit is typically assumed to be Poissonian noise, wherein the variance in mRNA numbers is equal to their mean. Here, we demonstrate that several cell division genes in fission yeast exhibit mRNA variances significantly below this limit. The reduced variance can be explained by a gene expression model incorporating multiple transcription and mRNA degradation steps. Notably, in this sub-Poissonian regime, distinct from Poissonian or super-Poissonian regimes, cytoplasmic noise is effectively suppressed through a higher mRNA export rate. Our findings redefine the lower limit of eukaryotic gene expression noise and uncover molecular requirements for achieving ultralow noise, which is expected to be important for vital cellular functions.

Characterization of Pik1 function in fission yeast reveals its conserved role in lipid synthesis and not cytokinesis

Phosphatidylinositol (PI)-4-phosphate (PI4P) is a lipid found at the plasma membrane (PM) and Golgi in cells from yeast to humans. PI4P is generated from PI by PI4-kinases and can be converted to PI-4,5-bisphosphate [PI(4,5)P 2 ]. Schizosaccharomyces pombe have 2 essential PI4-kinases: Stt4 and Pik1. Stt4 localizes to the PM and its loss from the PM results in a decrease of PM PI4P and PI(4,5)P 2 . As a result, cells divide non-medially due to disrupted cytokinetic ring-PM anchoring. However, the localization and function of S. pombe Pik1 has not been thoroughly examined. Here, we found that Pik1 localizes exclusively to the trans-Golgi and is required for Golgi PI4P production. We determined that Ncs1 regulates Pik1, but unlike in other organisms, it is not required for Pik1 Golgi localization. When Pik1 function was disrupted, PM PI4P but not PI(4,5)P 2 levels were reduced, a major difference with Stt4. We conclude that Stt4 is the chief enzyme responsible for producing the PI4P that generates PI(4,5)P 2 . Also, that cells with disrupted Pik1 do not divide asymmetrically highlights the specific importance of PM PI(4,5)P 2 for cytokinetic ring-PM anchoring.

Summary statement

Fission yeast Pik1 localizes exclusively to the trans-Golgi independently of Ncs1, where it contributes to PI4P but not PI(4,5)P 2 synthesis. Pik1 does not affect cytokinesis.

Status of genome function annotation in model organisms and crops

Since the entry into genome-enabled biology several decades ago, much progress has been made in determining, describing, and disseminating the functions of genes and their products. Yet, this information is still difficult to access for many scientists and for most genomes. To provide easy access and a graphical summary of the status of genome function annotation for model organisms and bioenergy and food crop species, we created a web application (https://genomeannotation.rheelab.org) to visualize, search, and download genome annotation data for 28 species. The summary graphics and data tables will be updated semi-annually, and snapshots will be archived to provide a historical record of the progress of genome function annotation efforts. Clear and simple visualization of up-to-date genome function annotation status, including the extent of what is unknown, will help address the grand challenge of elucidating the functions of all genes in organisms.

Substrate displacement of CK1 C-termini regulates kinase specificity

CK1 kinases participate in many signaling pathways; how these enzymes are regulated is therefore of significant biological consequence. CK1s autophosphorylate their C-terminal non-catalytic tails, and eliminating these modifications increases substrate phosphorylation in vitro, suggesting that the autophosphorylated C-termini act as inhibitory pseudosubstrates. To test this prediction, we comprehensively identified the autophosphorylation sites on Schizosaccharomyces pombe Hhp1 and human CK1ε. Peptides corresponding to the C-termini interacted with the kinase domains only when phosphorylated, and phosphoablating mutations increased Hhp1 and CK1ε activity towards substrates. Interestingly, substrates competitively inhibited binding of the autophosphorylated tails to the substrate binding grooves. The presence or absence of tail autophosphorylation influenced the catalytic efficiency with which CK1s targeted different substrates, indicating that tails contribute to substrate specificity. Combining this mechanism with autophosphorylation of the T220 site in the catalytic domain, we propose a displacement specificity model to describe how autophosphorylation regulates substrate specificity for the CK1 family.

Stn1-Ten1 and Taz1 independently promote replication of subtelomeric fragile sequences in fission yeast

Efficient replication of terminal DNA is crucial to maintain telomere stability. In fission yeast, Taz1 and the Stn1-Ten1 (ST) complex play prominent roles in DNA-ends replication. However, their function remains elusive. Here, we have analyzed genome-wide replication and show that ST does not affect genome-wide replication but is crucial for the efficient replication of a subtelomeric region called STE3-2. We further show that, when ST function is compromised, a homologous recombination (HR)-based fork restart mechanism becomes necessary for STE3-2 stability. While both Taz1 and Stn1 bind to STE3-2, we find that the STE3-2 replication function of ST is independent of Taz1 but relies on its association with the shelterin proteins Pot1-Tpz1-Poz1. Finally, we demonstrate that the firing of an origin normally inhibited by Rif1 can circumvent the replication defect of subtelomeres when ST function is compromised. Our results help illuminate why fission yeast telomeres are terminal fragile sites.

Optimization of energy production and central carbon metabolism in a non-respiring eukaryote

Most eukaryotes respire oxygen, using it to generate biomass and energy. However, a few organisms have lost the capacity to respire. Understanding how they manage biomass and energy production may illuminate the critical points at which respiration feeds into central carbon metabolism and explain possible routes to its optimization. Here, we use two related fission yeasts, Schizosaccharomyces pombe and Schizosaccharomyces japonicus, as a comparative model system. We show that although S. japonicus does not respire oxygen, unlike S. pombe, it is capable of efficient NADH oxidation, amino acid synthesis, and ATP generation. We probe possible optimization strategies through the use of stable isotope tracing metabolomics, mass isotopologue distribution analysis, genetics, and physiological experiments. S. japonicus appears to have optimized cytosolic NADH oxidation via glycerol-3-phosphate synthesis. It runs a fully bifurcated TCA pathway, sustaining amino acid production. Finally, we propose that it has optimized glycolysis to maintain high ATP/ADP ratio, in part by using the pentose phosphate pathway as a glycolytic shunt, reducing allosteric inhibition of glycolysis and supporting biomass generation. By comparing two related organisms with vastly different metabolic strategies, our work highlights the versatility and plasticity of central carbon metabolism in eukaryotes, illuminating critical adaptations supporting the preferential use of glycolysis over oxidative phosphorylation.

Comprehensive mutational analysis of the checkpoint signaling function of Rpa1/Ssb1 in fission yeast

Replication protein A (RPA) is a heterotrimeric complex and the major single-strand DNA (ssDNA) binding protein in eukaryotes. It plays important roles in DNA replication, repair, recombination, telomere maintenance, and checkpoint signaling. Because RPA is essential for cell survival, understanding its checkpoint signaling function in cells has been challenging. Several RPA mutants have been reported previously in fission yeast. None of them, however, has a defined checkpoint defect. A separation-of-function mutant of RPA, if identified, would provide significant insights into the checkpoint initiation mechanisms. We have explored this possibility and carried out an extensive genetic screen for Rpa1/Ssb1, the large subunit of RPA in fission yeast, looking for mutants with defects in checkpoint signaling. This screen has identified twenty-five primary mutants that are sensitive to genotoxins. Among these mutants, two have been confirmed partially defective in checkpoint signaling primarily at the replication fork, not the DNA damage site. The remaining mutants are likely defective in other functions such as DNA repair or telomere maintenance. Our screened mutants, therefore, provide a valuable tool for future dissection of the multiple functions of RPA in fission yeast.

Telomere-to-Telomere genome assemblies of human-infecting Encephalitozoon species

BACKGROUND

Microsporidia are diverse spore forming, fungal-related obligate intracellular pathogens infecting a wide range of hosts. This diversity is reflected at the genome level with sizes varying by an order of magnitude, ranging from less than 3 Mb in Encephalitozoon species (the smallest known in eukaryotes) to more than 50 Mb in Edhazardia spp. As a paradigm of genome reduction in eukaryotes, the small Encephalitozoon genomes have attracted much attention with investigations revealing gene dense, repeat- and intron-poor genomes characterized by a thorough pruning of molecular functions no longer relevant to their obligate intracellular lifestyle. However, because no Encephalitozoon genome has been sequenced from telomere-to-telomere and since no methylation data is available for these species, our understanding of their overall genetic and epigenetic architectures is incomplete.

METHODS

In this study, we sequenced the complete genomes from telomere-to-telomere of three human-infecting Encephalitozoon spp. -E. intestinalis ATCC 50506, E. hellem ATCC 50604 and E. cuniculi ATCC 50602- using short and long read platforms and leveraged the data generated as part of the sequencing process to investigate the presence of epigenetic markers in these genomes. We also used a mixture of sequence- and structure-based computational approaches, including protein structure prediction, to help identify which Encephalitozoon proteins are involved in telomere maintenance, epigenetic regulation, and heterochromatin formation.

RESULTS

The Encephalitozoon chromosomes were found capped by TTAGG 5-mer telomeric repeats followed by telomere associated repeat elements (TAREs) flanking hypermethylated ribosomal RNA (rRNA) gene loci featuring 5-methylcytosines (5mC) and 5-hemimethylcytosines (5hmC), themselves followed by lesser methylated subtelomeres and hypomethylated chromosome cores. Strong nucleotide biases were identified between the telomeres/subtelomeres and chromosome cores with significant changes in GC/AT, GT/AC and GA/CT contents. The presence of several genes coding for proteins essential to telomere maintenance, epigenetic regulation, and heterochromatin formation was further confirmed in the Encephalitozoon genomes.

CONCLUSION

Altogether, our results strongly support the subtelomeres as sites of heterochromatin formation in Encephalitozoon genomes and further suggest that these species might shutdown their energy-consuming ribosomal machinery while dormant as spores by silencing of the rRNA genes using both 5mC/5hmC methylation and facultative heterochromatin formation at these loci.

The Gene Ontology knowledgebase in 2023

The Gene Ontology (GO) knowledgebase (http://geneontology.org) is a comprehensive resource concerning the functions of genes and gene products (proteins and noncoding RNAs). GO annotations cover genes from organisms across the tree of life as well as viruses, though most gene function knowledge currently derives from experiments carried out in a relatively small number of model organisms. Here, we provide an updated overview of the GO knowledgebase, as well as the efforts of the broad, international consortium of scientists that develops, maintains, and updates the GO knowledgebase. The GO knowledgebase consists of three components: (1) the GO-a computational knowledge structure describing the functional characteristics of genes; (2) GO annotations-evidence-supported statements asserting that a specific gene product has a particular functional characteristic; and (3) GO Causal Activity Models (GO-CAMs)-mechanistic models of molecular "pathways" (GO biological processes) created by linking multiple GO annotations using defined relations. Each of these components is continually expanded, revised, and updated in response to newly published discoveries and receives extensive QA checks, reviews, and user feedback. For each of these components, we provide a description of the current contents, recent developments to keep the knowledgebase up to date with new discoveries, and guidance on how users can best make use of the data that we provide. We conclude with future directions for the project.

Fission yeast CK1 promotes DNA double-strand break repair through both homologous recombination and non-homologous end joining

The CK1 family are conserved serine/threonine kinases with numerous substrates and cellular functions. The fission yeast CK1 orthologues Hhp1 and Hhp2 were first characterized as regulators of DNA repair, but the mechanism(s) by which CK1 activity promotes DNA repair had not been investigated. Here, we found that deleting Hhp1 and Hhp2 or inhibiting CK1 catalytic activities in yeast or in human cells activated the DNA damage checkpoint due to persistent double-strand breaks (DSBs). The primary pathways to repair DSBs, homologous recombination and non-homologous end joining, were both less efficient in cells lacking Hhp1 and Hhp2 activity. In order to understand how Hhp1 and Hhp2 promote DSB repair, we identified new substrates using quantitative phosphoproteomics. We confirmed that Arp8, a component of the INO80 chromatin remodeling complex, is a bona fide substrate of Hhp1 and Hhp2 that is important for DSB repair. Our data suggest that Hhp1 and Hhp2 facilitate DSB repair by phosphorylating multiple substrates, including Arp8.

A high-quality reference genome for the fission yeast Schizosaccharomyces osmophilus

Fission yeasts are an ancient group of fungal species that diverged from each other from tens to hundreds of million years ago. Among them is the preeminent model organism Schizosaccharomyces pombe, which has significantly contributed to our understandings of molecular mechanisms underlying fundamental cellular processes. The availability of the genomes of S. pombe and 3 other fission yeast species S. japonicus, S. octosporus, and S. cryophilus has enabled cross-species comparisons that provide insights into the evolution of genes, pathways, and genomes. Here, we performed genome sequencing on the type strain of the recently identified fission yeast species S. osmophilus and obtained a complete mitochondrial genome and a nuclear genome assembly with gaps only at rRNA gene arrays. A total of 5,098 protein-coding nuclear genes were annotated and orthologs for more than 95% of them were identified. Genome-based phylogenetic analysis showed that S. osmophilus is most closely related to S. octosporus and these 2 species diverged around 16 million years ago. To demonstrate the utility of this S. osmophilus reference genome, we conducted cross-species comparative analyses of centromeres, telomeres, transposons, the mating-type region, Cbp1 family proteins, and mitochondrial genomes. These analyses revealed conservation of repeat arrangements and sequence motifs in centromere cores, identified telomeric sequences composed of 2 types of repeats, delineated relationships among Tf1/sushi group retrotransposons, characterized the evolutionary origins and trajectories of Cbp1 family domesticated transposases, and discovered signs of interspecific transfer of 2 types of mitochondrial selfish elements.

Protein S-palmitoylation regulates different stages of meiosis in Schizosaccharomyces pombe

Posttranslational protein S-palmitoylation regulates the localization and function of its target proteins involved in diverse cellular processes including meiosis. In this study, we demonstrate that S-palmitoylation mediated by Erf2-Erf4 and Akr1 palmitoylacyltransferases is required at multiple meiotic stages in the fission yeast Schizosaccharomyces pombe We find that S-palmitoylation by Erf2-Erf4 is required for Ras1 localization at the cell periphery to enrich at the cell conjugation site for mating pheromone response. In the absence of Erf2 or Erf4, mutant cells are sterile. A role of Akr1 S-palmitoylating the nuclear fusion protein Tht1 to function in karyogamy is identified. We demonstrate that S-palmitoylation stabilizes and localizes Tht1 to ER, interacting with Sey1 ER fusion GTPase for proper meiotic nuclear fusion. In akr1, tht1, or sey1 mutant, meiotic cells, haploid nuclei are unfused with subsequent chromosome segregation defects. Erf2-Erf4 has an additional substrate of the spore coat protein Isp3. In the absence of Erf2, Isp3 is mislocalized from the spore coat. Together, these results highlight the versatility of the cellular processes in which protein S-palmitoylation participates.

Multiple polarity kinases inhibit phase separation of F-BAR protein Cdc15 and antagonize cytokinetic ring assembly in fission yeast

The F-BAR protein Cdc15 is essential for cytokinesis in Schizosaccharomyces pombe and plays a key role in attaching the cytokinetic ring (CR) to the plasma membrane (PM). Cdc15's abilities to bind to the membrane and oligomerize via its F-BAR domain are inhibited by phosphorylation of its intrinsically disordered region (IDR). Multiple cell polarity kinases regulate Cdc15 IDR phosphostate, and of these the DYRK kinase Pom1 phosphorylation sites on Cdc15 have been shown in vivo to prevent CR formation at cell tips. Here, we compared the ability of Pom1 to control Cdc15 phosphostate and cortical localization to that of other Cdc15 kinases: Kin1, Pck1, and Shk1. We identified distinct but overlapping cohorts of Cdc15 phosphorylation sites targeted by each kinase, and the number of sites correlated with each kinases' abilities to influence Cdc15 PM localization. Coarse-grained simulations predicted that cumulative IDR phosphorylation moves the IDRs of a dimer apart and toward the F-BAR tips. Further, simulations indicated that the overall negative charge of phosphorylation masks positively charged amino acids necessary for F-BAR oligomerization and membrane interaction. Finally, simulations suggested that dephosphorylated Cdc15 undergoes phase separation driven by IDR interactions. Indeed, dephosphorylated but not phosphorylated Cdc15 undergoes liquid-liquid phase separation to form droplets in vitro that recruit Cdc15 binding partners. In cells, Cdc15 phosphomutants also formed PM-bound condensates that recruit other CR components. Together, we propose that a threshold of Cdc15 phosphorylation by assorted kinases prevents Cdc15 condensation on the PM and antagonizes CR assembly.

Translation-complex profiling of fission yeast cells reveals dynamic rearrangements of scanning ribosomal subunits upon nutritional stress

Control of mRNA translation is key for stress responses. Translation initiation is usually rate-limiting and, in eukaryotes, involves mRNA scanning by the small ribosomal subunit. Despite its importance, many aspects of translation in vivo have not been explored fully, especially at the transcriptome-wide level. A recent method termed translation-complex profiling (TCP-seq) allows transcriptome-wide views of scanning ribosomal subunits. We applied TCP-seq to nutritional stress in the fission yeast Schizosaccharomyces pombe. At initiation sites, we observed multiple complexes resembling those of mammals, and consistent with queuing of scanning subunits. In 5' UTRs, small subunit accumulations were common and may reflect impediments to scanning. A key mediator of stress responses in S. pombe is the Fil1 transcription factor, which is regulated translationally by a poorly-understood mechanism involving upstream Open Reading Frames (uORFs). TCP-seq data of fil1 shows that stress allows scanning subunits to by-pass specific uORFs and reach the fil1 coding sequence. The integration of these observations with reporter assays revealed that fil1 translational control is mediated by a combination of scanning reinitiation-repressive and permissive uORFs, and establishes fil1 as a model for uORF-mediated translational control. Altogether, our transcriptome-wide study reveals general and gene-specific features of translation in a model eukaryote.

Fission Yeast Rho1p-GEFs: From Polarity and Cell Wall Synthesis to Genome Stability

Rho1p is a membrane-associated protein that belongs to the Rho family of small GTPases. These proteins coordinate processes such as actin remodelling and polarised secretion to maintain the shape and homeostasis of yeast cells. In response to extracellular stimuli, Rho1p undergoes conformational switching between a guanosine triphosphate (GTP)-bound active state and a guanosine diphosphate (GDP)-bound inactive state. Cycling is improved with guanine nucleotide exchange factor (GEF) activity necessary to activate signalling and GTPase activating protein (GAP) activity required for subsequent signal depletion. This review focuses on fission yeast Rho1p GEFs, Rgf1p, Rgf2p, and Rgf3p that belong to the family of DH-PH domain-containing Dbl-related GEFs. They are multi-domain proteins that detect biological signals that induce or inhibit their catalytic activity over Rho1p. Each of them activates Rho1p in different places and times. Rgf1p acts preferentially during polarised growth. Rgf2p is required for sporulation, and Rgf3p plays an essential function in septum synthesis. In addition, we outline the noncanonical roles of Rho1p-GEFs in genomic instability.

YEASTRACT+: a portal for the exploitation of global transcription regulation and metabolic model data in yeast biotechnology and pathogenesis

YEASTRACT+ (http://yeastract-plus.org/) is a tool for the analysis, prediction and modelling of transcription regulatory data at the gene and genomic levels in yeasts. It incorporates three integrated databases: YEASTRACT (http://yeastract-plus.org/yeastract/), PathoYeastract (http://yeastract-plus.org/pathoyeastract/) and NCYeastract (http://yeastract-plus.org/ncyeastract/), focused on Saccharomyces cerevisiae, pathogenic yeasts of the Candida genus, and non-conventional yeasts of biotechnological relevance. In this release, YEASTRACT+ offers upgraded information on transcription regulation for the ten previously incorporated yeast species, while extending the database to another pathogenic yeast, Candida auris. Since the last release of YEASTRACT+ (January 2020), a fourth database has been integrated. CommunityYeastract (http://yeastract-plus.org/community/) offers a platform for the creation, use, and future update of YEASTRACT-like databases for any yeast of the users' choice. CommunityYeastract currently provides information for two Saccharomyces boulardii strains, Rhodotorula toruloides NP11 oleaginous yeast, and Schizosaccharomyces pombe 972h-. In addition, YEASTRACT+ portal currently gathers 304 547 documented regulatory associations between transcription factors (TF) and target genes and 480 DNA binding sites, considering 2771 TFs from 11 yeast species. A new set of tools, currently implemented for S. cerevisiae and C. albicans, is further offered, combining regulatory information with genome-scale metabolic models to provide predictions on the most promising transcription factors to be exploited in cell factory optimisation or to be used as novel drug targets. The expansion of these new tools to the remaining YEASTRACT+ species is ongoing.