A challenge in solving the genotype-to-phenotype relationship is to predict a cell's metabolome, believed to correlate poorly with gene expression. Using comparative quantitative proteomics, we found that differential protein expression in 97 Saccharomyces cerevisiae kinase deletion strains is non-redundant and dominated by abundance changes in metabolic enzymes. Associating differential enzyme expression landscapes to corresponding metabolomes using network models provided reasoning for poor proteome-metabolome correlations; differential protein expression redistributes flux control between many enzymes acting in concert, a mechanism not captured by one-to-one correlation statistics. Mapping these regulatory patterns using machine learning enabled the prediction of metabolite concentrations, as well as identification of candidate genes important for the regulation of metabolism. Overall, our study reveals that a large part of metabolism regulation is explained through coordinated enzyme expression changes. Our quantitative data indicate that this mechanism explains more than half of metabolism regulation and underlies the interdependency between enzyme levels and metabolism, which renders the metabolome a predictable phenotype.

Target of Rapamycin Complex 1 (TORC1) signaling promotes growth and aging. Inhibition of TORC1 leads to reduced protein translation, which promotes longevity. TORC1-dependent post-transcriptional regulation of protein translation has been well studied, while analogous transcriptional regulation is less understood. Here we screen fission yeast mutants for resistance to Torin1, which inhibits TORC1 and cell growth. Cells lacking the GATA factor Gaf1 (gaf1Δ) grow normally even in high doses of Torin1. The gaf1Δ mutation shortens the chronological lifespan of non-dividing cells and diminishes Torin1-mediated longevity. Expression profiling and genome-wide binding experiments show that upon TORC1 inhibition, Gaf1 directly upregulates genes for small-molecule metabolic pathways and indirectly represses genes for protein translation. Surprisingly, Gaf1 binds to and downregulates the tRNA genes, so it also functions as a transcription factor for RNA polymerase III. Thus, Gaf1 controls the transcription of both protein-coding and tRNA genes to inhibit translation and growth downstream of TORC1.

Eukaryotic genomes express numerous long intergenic non-coding RNAs (lincRNAs) that do not overlap any coding genes. Some lincRNAs function in various aspects of gene regulation, but it is not clear in general to what extent lincRNAs contribute to the information flow from genotype to phenotype. To explore this question, we systematically analysed cellular roles of lincRNAs in Schizosaccharomyces pombe. Using seamless CRISPR/Cas9-based genome editing, we deleted 141 lincRNA genes to broadly phenotype these mutants, together with 238 diverse coding-gene mutants for functional context. We applied high-throughput colony-based assays to determine mutant growth and viability in benign conditions and in response to 145 different nutrient, drug, and stress conditions. These analyses uncovered phenotypes for 47.5% of the lincRNAs and 96% of the protein-coding genes. For 110 lincRNA mutants, we also performed high-throughput microscopy and flow cytometry assays, linking 37% of these lincRNAs with cell-size and/or cell-cycle control. With all assays combined, we detected phenotypes for 84 (59.6%) of all lincRNA deletion mutants tested. For complementary functional inference, we analysed colony growth of strains ectopically overexpressing 113 lincRNA genes under 47 different conditions. Of these overexpression strains, 102 (90.3%) showed altered growth under certain conditions. Clustering analyses provided further functional clues and relationships for some of the lincRNAs. These rich phenomics datasets associate lincRNA mutants with hundreds of phenotypes, indicating that most of the lincRNAs analysed exert cellular functions in specific environmental or physiological contexts. This study provides groundwork to further dissect the roles of these lincRNAs in the relevant conditions.

Aberrant repair of DNA double-strand breaks can recombine distant chromosomal breakpoints. Chromosomal rearrangements compromise genome function and are a hallmark of ageing. Rearrangements are challenging to detect in non-dividing cell populations, because they reflect individually rare, heterogeneous events. The genomic distribution of de novo rearrangements in non-dividing cells, and their dynamics during ageing, remain therefore poorly characterized. Studies of genomic instability during ageing have focussed on mitochondrial DNA, small genetic variants, or proliferating cells. To characterize genome rearrangements during cellular ageing in non-dividing cells, we interrogated a single diagnostic measure, DNA breakpoint junctions, using Schizosaccharomyces pombe as a model system. Aberrant DNA junctions that accumulated with age were associated with microhomology sequences and R-loops. Global hotspots for age-associated breakpoint formation were evident near telomeric genes and linked to remote breakpoints elsewhere in the genome, including the mitochondrial chromosome. Formation of breakpoint junctions at global hotspots was inhibited by the Sir2 histone deacetylase and might be triggered by an age-dependent de-repression of chromatin silencing. An unexpected mechanism of genomic instability may cause more local hotspots: age-associated reduction in an RNA-binding protein triggering R-loops at target loci. This result suggests that biological processes other than transcription or replication can drive genome rearrangements. Notably, we detected similar signatures of genome rearrangements that accumulated in old brain cells of humans. These findings provide insights into the unique patterns and possible mechanisms of genome rearrangements in non-dividing cells, which can be promoted by ageing-related changes in gene-regulatory proteins.

Cells balance glycolysis with respiration to support their metabolic needs in different environmental or physiological contexts. With abundant glucose, many cells prefer to grow by aerobic glycolysis or fermentation. Using 161 natural isolates of fission yeast, we investigated the genetic basis and phenotypic effects of the fermentation-respiration balance. The laboratory and a few other strains depended more on respiration. This trait was associated with a single nucleotide polymorphism in a conserved region of Pyk1, the sole pyruvate kinase in fission yeast. This variant reduced Pyk1 activity and glycolytic flux. Replacing the "low-activity" pyk1 allele in the laboratory strain with the "high-activity" allele was sufficient to increase fermentation and decrease respiration. This metabolic rebalancing triggered systems-level adjustments in the transcriptome and proteome and in cellular traits, including increased growth and chronological lifespan but decreased resistance to oxidative stress. Thus, low Pyk1 activity does not lead to a growth advantage but to stress tolerance. The genetic tuning of glycolytic flux may reflect an adaptive trade-off in a species lacking pyruvate kinase isoforms.

Genetically identical cells are known to differ in many physiological parameters such as growth rate and drug tolerance. Metabolic specialization is believed to be a cause of such phenotypic heterogeneity, but detection of metabolically divergent subpopulations remains technically challenging. We developed a proteomics-based technology, termed differential isotope labelling by amino acids (DILAC), that can detect producer and consumer subpopulations of a particular amino acid within an isogenic cell population by monitoring peptides with multiple occurrences of the amino acid. We reveal that young, morphologically undifferentiated yeast colonies contain subpopulations of lysine producers and consumers that emerge due to nutrient gradients. Deconvoluting their proteomes using DILAC, we find evidence for in situ cross-feeding where rapidly growing cells ferment and provide the more slowly growing, respiring cells with ethanol. Finally, by combining DILAC with fluorescence-activated cell sorting, we show that the metabolic subpopulations diverge phenotypically, as exemplified by a different tolerance to the antifungal drug amphotericin B. Overall, DILAC captures previously unnoticed metabolic heterogeneity and provides experimental evidence for the role of metabolic specialization and cross-feeding interactions as a source of phenotypic heterogeneity in isogenic cell populations.

Genetic factors underlying aging are remarkably conserved from yeast to human. The fission yeast Schizosaccharomyces pombe is an emerging genetic model to analyze cellular aging. Chronological lifespan (CLS) has been studied in stationary-phase yeast cells depleted for glucose, which only survive for a few days. Here, we analyzed CLS in quiescent S. pombe cells deprived of nitrogen, which arrest in a differentiated, G0-like state and survive for more than 2 months. We applied parallel mutant phenotyping by barcode sequencing (Bar-seq) to assay pooled haploid deletion mutants as they aged together during long-term quiescence. As expected, mutants with defects in autophagy or quiescence were under-represented or not detected. Lifespan scores could be calculated for 1199 mutants. We focus the discussion on the 48 most long-lived mutants, including both known aging genes in other model systems and genes not previously implicated in aging. Genes encoding membrane proteins were particularly prominent as pro-aging factors. We independently verified the extended CLS in individual assays for 30 selected mutants, showing the efficacy of the screen. We also applied Bar-seq to profile all pooled deletion mutants for proliferation under a standard growth condition. Unlike for stationary-phase cells, no inverse correlation between growth and CLS of quiescent cells was evident. These screens provide a rich resource for further studies, and they suggest that the quiescence model can provide unique, complementary insights into cellular aging.

Cyclin-dependent kinase 9 (Cdk9) promotes elongation by RNA polymerase II (RNAPII), mRNA processing, and co-transcriptional histone modification. Cdk9 phosphorylates multiple targets, including the conserved RNAPII elongation factor Spt5 and RNAPII itself, but how these different modifications mediate Cdk9 functions is not known. Here we describe two Cdk9-dependent pathways in the fission yeast Schizosaccharomyces pombe that involve distinct targets and elicit distinct biological outcomes. Phosphorylation of Spt5 by Cdk9 creates a direct binding site for Prf1/Rtf1, a transcription regulator with functional and physical links to the Polymerase Associated Factor (PAF) complex. PAF association with chromatin is also dependent on Cdk9 but involves alternate phosphoacceptor targets. Prf1 and PAF are biochemically separate in cell extracts, and genetic analyses show that Prf1 and PAF are functionally distinct and exert opposing effects on the RNAPII elongation complex. We propose that this opposition constitutes a Cdk9 auto-regulatory mechanism, such that a positive effect on elongation, driven by the PAF pathway, is kept in check by a negative effect of Prf1/Rtf1 and downstream mono-ubiquitylation of histone H2B. Thus, optimal RNAPII elongation may require balanced action of functionally distinct Cdk9 pathways.

Heterochromatin is essential for chromosome segregation, gene silencing and genome integrity. The fission yeast Schizosaccharomyces pombe contains heterochromatin at centromeres, subtelomeres, and mating type genes, as well as at small islands of meiotic genes dispersed across the genome. This heterochromatin is generated by partially redundant mechanisms, including the production of small interfering RNAs (siRNAs) that are incorporated into the RITS protein complex (RNAi-Induced Transcriptional Silencing). The assembly of heterochromatin islands requires the function of the RNA-binding protein Mmi1, which recruits RITS to its mRNA targets and to heterochromatin islands. In addition, Mmi1 directs its targets to an exosome-dependent RNA elimination pathway.

Both single and multicellular organisms depend on anti-stress mechanisms that enable them to deal with sudden changes in the environment, including exposure to heat and oxidants. Central to the stress response are dynamic changes in metabolism, such as the transition from the glycolysis to the pentose phosphate pathway-a conserved first-line response to oxidative insults1,2. Here we report a second metabolic adaptation that protects microbial cells in stress situations. The role of the yeast polyamine transporter Tpo1p3-5 in maintaining oxidant resistance is unknown6. However, a proteomic time-course experiment suggests a link to lysine metabolism. We reveal a connection between polyamine and lysine metabolism during stress situations, in the form of a promiscuous enzymatic reaction in which the first enzyme of the polyamine pathway, Spe1p, decarboxylates lysine and forms an alternative polyamine, cadaverine. The reaction proceeds in the presence of extracellular lysine, which is taken up by cells to reach concentrations up to one hundred times higher than those required for growth. Such extensive harvest is not observed for the other amino acids, is dependent on the polyamine pathway and triggers a reprogramming of redox metabolism. As a result, NADPH-which would otherwise be required for lysine biosynthesis-is channelled into glutathione metabolism, leading to a large increase in glutathione concentrations, lower levels of reactive oxygen species and increased oxidant tolerance. Our results show that nutrient uptake occurs not only to enable cell growth, but when the nutrient availability is favourable it also enables cells to reconfigure their metabolism to preventatively mount stress protection.

During meiosis, tethering of parental mitochondria to opposite cell poles inhibits the mixing of mitochondria with different genomes and ensures uniparental inheritance in thestandard laboratory strain of fission yeast. We here investigate mitochondrial inheritance in crosses between natural isolates using tetrad dissection and next-generation sequencing. We find that colonies grown from single spores can sometimes carry a mix of mitochondrial genotypes, that mitochondrial genomes can recombine during meiosis, that in some cases tetrads do not follow the 2:2 segregation pattern, and that certain crosses may feature a weak bias towards one of the parents. Together, these findings paint a more nuanced picture of mitochondrial inheritance in the wild.

Functional genomic strategies have become fundamental for annotating gene function and regulatory networks. Here, we combined functional genomics with proteomics by quantifying protein abundances in a genome-scale knockout library in Saccharomyces cerevisiae, using data-independent acquisition mass spectrometry. We find that global protein expression is driven by a complex interplay of (1) general biological properties, including translation rate, protein turnover, the formation of protein complexes, growth rate, and genome architecture, followed by (2) functional properties, such as the connectivity of a protein in genetic, metabolic, and physical interaction networks. Moreover, we show that functional proteomics complements current gene annotation strategies through the assessment of proteome profile similarity, protein covariation, and reverse proteome profiling. Thus, our study reveals principles that govern protein expression and provides a genome-spanning resource for functional annotation.

The complexity of many cellular and organismal traits results from the integration of genetic and environmental factors via molecular networks. Network structure and effect propagation are best understood at the level of functional modules, but so far, no concept has been established to include the global network state. Here, we show when and how genetic perturbations lead to molecular changes that are confined to small parts of a network versus when they lead to modulation of network states. Integrating multi-omics profiling of genetically heterogeneous budding and fission yeast strains with an array of cellular traits identified a central state transition of the yeast molecular network that is related to PKA and TOR (PT) signaling. Genetic variants affecting this PT state globally shifted the molecular network along a single-dimensional axis, thereby modulating processes including energy and amino acid metabolism, transcription, translation, cell cycle control, and cellular stress response. We propose that genetic effects can propagate through large parts of molecular networks because of the functional requirement to centrally coordinate the activity of fundamental cellular processes.

Many proteins remain poorly characterized even in well-studied organisms, presenting a bottleneck for research. We applied phenomics and machine-learning approaches with Schizosaccharomyces pombe for broad cues on protein functions. We assayed colony-growth phenotypes to measure the fitness of deletion mutants for 3509 non-essential genes in 131 conditions with different nutrients, drugs, and stresses. These analyses exposed phenotypes for 3492 mutants, including 124 mutants of 'priority unstudied' proteins conserved in humans, providing varied functional clues. For example, over 900 proteins were newly implicated in the resistance to oxidative stress. Phenotype-correlation networks suggested roles for poorly characterized proteins through 'guilt by association' with known proteins. For complementary functional insights, we predicted Gene Ontology (GO) terms using machine learning methods exploiting protein-network and protein-homology data (NET-FF). We obtained 56,594 high-scoring GO predictions, of which 22,060 also featured high information content. Our phenotype-correlation data and NET-FF predictions showed a strong concordance with existing PomBase GO annotations and protein networks, with integrated analyses revealing 1675 novel GO predictions for 783 genes, including 47 predictions for 23 priority unstudied proteins. Experimental validation identified new proteins involved in cellular aging, showing that these predictions and phenomics data provide a rich resource to uncover new protein functions.

Our current understanding of how natural genetic variation affects gene expression beyond well-annotated coding genes is still limited. The use of deep sequencing technologies for the study of expression quantitative trait loci (eQTLs) has the potential to close this gap. Here, we generated the first recombinant strain library for fission yeast and conducted an RNA-seq-based QTL study of the coding, non-coding, and antisense transcriptomes. We show that the frequency of distal effects (trans-eQTLs) greatly exceeds the number of local effects (cis-eQTLs) and that non-coding RNAs are as likely to be affected by eQTLs as protein-coding RNAs. We identified a genetic variation of swc5 that modifies the levels of 871 RNAs, with effects on both sense and antisense transcription, and show that this effect most likely goes through a compromised deposition of the histone variant H2A.Z. The strains, methods, and datasets generated here provide a rich resource for future studies.

In the fission yeast Schizosaccharomyces pombe the prevailing approach for gene manipulations is based on homologous recombination of a PCR product that contains genomic target sequences and a selectable marker. The CRISPR/Cas9 system has recently been implemented in fission yeast, which allows for seamless genome editing without integration of a selection marker or leaving any other genomic 'scars'. The published method involves manual design of the single guide RNA (sgRNA), and digestion of a large plasmid with a problematic restriction enzyme to clone the sgRNA. To increase the efficiency of this approach, we have established and optimized a PCR-based system to clone the sgRNA without restriction enzymes into a plasmid with a dominant natMX6 (nourseothricin) selection marker. We also provide a web-tool, CRISPR4P, to support the design of the sgRNAs and the primers required for the entire process of seamless DNA deletion. Moreover, we report the preparation of G1-synchronized and cryopreserved S. pombe cells, which greatly increases the efficiency and speed for transformations, and may also facilitate standard gene manipulations. Applying this optimized CRISPR/Cas9-based approach, we have successfully deleted over 80 different non-coding RNA genes, which are generally lowly expressed, and have inserted 7 point mutations in 4 different genomic regions.

The function of transcription in dynamic gene expression programs has been extensively studied, but little is known about how it is integrated with RNA turnover at the genome-wide level. We investigated these questions using the meiotic gene expression program of Schizosaccharomyces pombe. We identified over 80 transcripts that co-purify with the meiotic-specific Meu5p RNA-binding protein. Their levels and half-lives were reduced in meu5 mutants, demonstrating that Meu5p stabilizes its targets. Most Meu5p-bound RNAs were also targets of the Mei4p transcription factor, which induces the transient expression of approximately 500 meiotic genes. Although many Mei4p targets showed sharp expression peaks, Meu5p targets had broad expression profiles. In the absence of meu5, all Mei4p targets were expressed with similar kinetics, indicating that Meu5p alters the global features of the gene expression program. As Mei4p activates meu5 transcription, Mei4p, Meu5p and their common targets form a feed-forward loop, a motif common in transcriptional networks but not studied in the context of mRNA decay. Our data provide insight into the topology of regulatory networks integrating transcriptional and posttranscriptional controls.

The relationship between DNA sequence, biochemical function, and molecular evolution is relatively well-described for protein-coding regions of genomes, but far less clear in noncoding regions, particularly, in eukaryote genomes. In part, this is because we lack a complete description of the essential noncoding elements in a eukaryote genome. To contribute to this challenge, we used saturating transposon mutagenesis to interrogate the Schizosaccharomyces pombe genome. We generated 31 million transposon insertions, a theoretical coverage of 2.4 insertions per genomic site. We applied a five-state hidden Markov model (HMM) to distinguish insertion-depleted regions from insertion biases. Both raw insertion-density and HMM-defined fitness estimates showed significant quantitative relationships to gene knockout fitness, genetic diversity, divergence, and expected functional regions based on transcription and gene annotations. Through several analyses, we conclude that transposon insertions produced fitness effects in 66-90% of the genome, including substantial portions of the noncoding regions. Based on the HMM, we estimate that 10% of the insertion depleted sites in the genome showed no signal of conservation between species and were weakly transcribed, demonstrating limitations of comparative genomics and transcriptomics to detect functional units. In this species, 3'- and 5'-untranslated regions were the most prominent insertion-depleted regions that were not represented in measures of constraint from comparative genomics. We conclude that the combination of transposon mutagenesis, evolutionary, and biochemical data can provide new insights into the relationship between genome function and molecular evolution.

Mono-ubiquitylation of histone H2B (H2Bub1) and phosphorylation of elongation factor Spt5 by cyclin-dependent kinase 9 (Cdk9) occur during transcription by RNA polymerase II (RNAPII), and are mutually dependent in fission yeast. It remained unclear whether Cdk9 and H2Bub1 cooperate to regulate the expression of individual genes. Here, we show that Cdk9 inhibition or H2Bub1 loss induces intragenic antisense transcription of ∼10% of fission yeast genes, with each perturbation affecting largely distinct subsets; ablation of both pathways de-represses antisense transcription of over half the genome. H2Bub1 and phospho-Spt5 have similar genome-wide distributions; both modifications are enriched, and directly proportional to each other, in coding regions, and decrease abruptly around the cleavage and polyadenylation signal (CPS). Cdk9-dependence of antisense suppression at specific genes correlates with high H2Bub1 occupancy, and with promoter-proximal RNAPII pausing. Genetic interactions link Cdk9, H2Bub1 and the histone deacetylase Clr6-CII, while combined Cdk9 inhibition and H2Bub1 loss impair Clr6-CII recruitment to chromatin and lead to decreased occupancy and increased acetylation of histones within gene coding regions. These results uncover novel interactions between co-transcriptional histone modification pathways, which link regulation of RNAPII transcription elongation to suppression of aberrant initiation.

Ageing-related processes are largely conserved, with simple organisms remaining the main platform to discover and dissect new ageing-associated genes. Yeasts provide potent model systems to study cellular ageing owing their amenability to systematic functional assays under controlled conditions. Even with yeast cells, however, ageing assays can be laborious and resource-intensive. Here we present improved experimental and computational methods to study chronological lifespan in Schizosaccharomyces pombe. We decoded the barcodes for 3206 mutants of the latest gene-deletion library, enabling the parallel profiling of ~700 additional mutants compared to previous screens. We then applied a refined method of barcode sequencing (Bar-seq), addressing technical and statistical issues raised by persisting DNA in dead cells and sampling bottlenecks in aged cultures, to screen for mutants showing altered lifespan during stationary phase. This screen identified 341 long-lived mutants and 1246 short-lived mutants which point to many previously unknown ageing-associated genes, including 46 conserved but entirely uncharacterized genes. The ageing-associated genes showed coherent enrichments in processes also associated with human ageing, particularly with respect to ageing in non-proliferative brain cells. We also developed an automated colony-forming unit assay to facilitate medium- to high-throughput chronological-lifespan studies by saving time and resources compared to the traditional assay. Results from the Bar-seq screen showed good agreement with this new assay. This study provides an effective methodological platform and identifies many new ageing-associated genes as a framework for analysing cellular ageing in yeast and beyond.