Quantitative analysis of fission yeast transcriptomes and proteomes in proliferating and quiescent cells.

Data on absolute molecule numbers will empower the modeling, understanding, and comparison of cellular functions and biological systems. We quantified transcriptomes and proteomes in fission yeast during cellular proliferation and quiescence. This rich resource provides the first comprehensive reference for all RNA and most protein concentrations in a eukaryote under two key physiological conditions. The integrated data set supports quantitative biology and affords unique insights into cell regulation. Although mRNAs are typically expressed in a narrow range above 1 copy/cell, most long, noncoding RNAs, except for a distinct subset, are tightly repressed below 1 copy/cell. Cell-cycle-regulated transcription tunes mRNA numbers to phase-specific requirements but can also bring about more switch-like expression. Proteins greatly exceed mRNAs in abundance and dynamic range, and concentrations are regulated to functional demands. Upon transition to quiescence, the proteome changes substantially, but, in stark contrast to mRNAs, proteins do not uniformly decrease but scale with cell volume.

Critical importance of DNA binding for CSL protein functions in fission yeast.

CSL proteins [named after the homologs CBF1 (RBP-Jκ in mice), Suppressor of Hairless and LAG-1] are conserved transcription factors found in animals and fungi. In the fission yeast Schizosaccharomyces pombe, they regulate various cellular processes, including cell cycle progression, lipid metabolism and cell adhesion. CSL proteins bind to DNA through their N-terminal Rel-like domain and central ß-trefoil domain. Here, we investigated the importance of DNA binding for CSL protein functions in fission yeast. We created CSL protein mutants with disrupted DNA binding and found that the vast majority of CSL protein functions depend on intact DNA binding. Specifically, DNA binding is crucial for the regulation of cell adhesion, lipid metabolism, cell cycle progression, long non-coding RNA expression and genome integrity maintenance. Interestingly, perturbed lipid metabolism leads to chromatin structure changes, potentially linking lipid metabolism to the diverse phenotypes associated with CSL protein functions. Our study highlights the critical role of DNA binding for CSL protein functions in fission yeast.

Perturbed fatty-acid metabolism is linked to localized chromatin hyperacetylation, increased stress-response gene expression and resistance to oxidative stress.

Oxidative stress is associated with cardiovascular and neurodegenerative diseases, diabetes, cancer, psychiatric disorders and aging. In order to counteract, eliminate and/or adapt to the sources of stress, cells possess elaborate stress-response mechanisms, which also operate at the level of regulating transcription. Interestingly, it is becoming apparent that the metabolic state of the cell and certain metabolites can directly control the epigenetic information and gene expression. In the fission yeast Schizosaccharomyces pombe, the conserved Sty1 stress-activated protein kinase cascade is the main pathway responding to most types of stresses, and regulates the transcription of hundreds of genes via the Atf1 transcription factor. Here we report that fission yeast cells defective in fatty acid synthesis (cbf11, mga2 and ACC/cut6 mutants; FAS inhibition) show increased expression of a subset of stress-response genes. This altered gene expression depends on Sty1-Atf1, the Pap1 transcription factor, and the Gcn5 and Mst1 histone acetyltransferases, is associated with increased acetylation of histone H3 at lysine 9 in the corresponding gene promoters, and results in increased cellular resistance to oxidative stress. We propose that changes in lipid metabolism can regulate the chromatin and transcription of specific stress-response genes, which in turn might help cells to maintain redox homeostasis.

Genetic effects on molecular network states explain complex traits.

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.

R-loops and regulatory changes in chronologically ageing fission yeast cells drive non-random patterns of genome rearrangements.

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.

The Nrd1-like protein Seb1 coordinates cotranscriptional 3' end processing and polyadenylation site selection.

Termination of RNA polymerase II (RNAPII) transcription is associated with RNA 3' end formation. For coding genes, termination is initiated by the cleavage/polyadenylation machinery. In contrast, a majority of noncoding transcription events in Saccharomyces cerevisiae does not rely on RNA cleavage for termination but instead terminates via a pathway that requires the Nrd1-Nab3-Sen1 (NNS) complex. Here we show that the Schizosaccharomyces pombe ortholog of Nrd1, Seb1, does not function in NNS-like termination but promotes polyadenylation site selection of coding and noncoding genes. We found that Seb1 associates with 3' end processing factors, is enriched at the 3' end of genes, and binds RNA motifs downstream from cleavage sites. Importantly, a deficiency in Seb1 resulted in widespread changes in 3' untranslated region (UTR) length as a consequence of increased alternative polyadenylation. Given that Seb1 levels affected the recruitment of conserved 3' end processing factors, our findings indicate that the conserved RNA-binding protein Seb1 cotranscriptionally controls alternative polyadenylation.

Mitochondrial respiration is required to provide amino acids during fermentative proliferation of fission yeast.

When glucose is available, many organisms repress mitochondrial respiration in favour of aerobic glycolysis, or fermentation in yeast, that suffices for ATP production. Fission yeast cells, however, rely partially on respiration for rapid proliferation under fermentative conditions. Here, we determined the limiting factors that require respiratory function during fermentation. When inhibiting the electron transport chain, supplementation with arginine was necessary and sufficient to restore rapid proliferation. Accordingly, a systematic screen for mutants growing poorly without arginine identified mutants defective in mitochondrial oxidative metabolism. Genetic or pharmacological inhibition of respiration triggered a drop in intracellular levels of arginine and amino acids derived from the Krebs cycle metabolite alpha-ketoglutarate: glutamine, lysine and glutamic acid. Conversion of arginine into these amino acids was required for rapid proliferation when blocking the respiratory chain. The respiratory block triggered an immediate gene expression response diagnostic of TOR inhibition, which was muted by arginine supplementation or without the AMPK-activating kinase Ssp1. The TOR-controlled proteins featured biased composition of amino acids reflecting their shortage after respiratory inhibition. We conclude that respiration supports rapid proliferation in fermenting fission yeast cells by boosting the supply of Krebs cycle-derived amino acids.

Cdk9 and H2Bub1 signal to Clr6-CII/Rpd3S to suppress aberrant antisense transcription.

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.

Longevity is determined by ETS transcription factors in multiple tissues and diverse species.

Ageing populations pose one of the main public health crises of our time. Reprogramming gene expression by altering the activities of sequence-specific transcription factors (TFs) can ameliorate deleterious effects of age. Here we explore how a circuit of TFs coordinates pro-longevity transcriptional outcomes, which reveals a multi-tissue and multi-species role for an entire protein family: the E-twenty-six (ETS) TFs. In Drosophila, reduced insulin/IGF signalling (IIS) extends lifespan by coordinating activation of Aop, an ETS transcriptional repressor, and Foxo, a Forkhead transcriptional activator. Aop and Foxo bind the same genomic loci, and we show that, individually, they effect similar transcriptional programmes in vivo. In combination, Aop can both moderate or synergise with Foxo, dependent on promoter context. Moreover, Foxo and Aop oppose the gene-regulatory activity of Pnt, an ETS transcriptional activator. Directly knocking down Pnt recapitulates aspects of the Aop/Foxo transcriptional programme and is sufficient to extend lifespan. The lifespan-limiting role of Pnt appears to be balanced by a requirement for metabolic regulation in young flies, in which the Aop-Pnt-Foxo circuit determines expression of metabolic genes, and Pnt regulates lipolysis and responses to nutrient stress. Molecular functions are often conserved amongst ETS TFs, prompting us to examine whether other Drosophila ETS-coding genes may also affect ageing. We show that five out of eight Drosophila ETS TFs play a role in fly ageing, acting from a range of organs and cells including the intestine, adipose and neurons. We expand the repertoire of lifespan-limiting ETS TFs in C. elegans, confirming their conserved function in ageing and revealing that the roles of ETS TFs in physiology and lifespan are conserved throughout the family, both within and between species.

An essential role for dNTP homeostasis following CDK-induced replication stress.

Replication stress is a common feature of cancer cells, and thus a potentially important therapeutic target. Here, we show that cyclin-dependent kinase (CDK)-induced replication stress, resulting from Wee1 inactivation, is synthetic lethal with mutations disrupting dNTP homeostasis in fission yeast. Wee1 inactivation leads to increased dNTP demand and replication stress through CDK-induced firing of dormant replication origins. Subsequent dNTP depletion leads to inefficient DNA replication, DNA damage and to genome instability. Cells respond to this replication stress by increasing dNTP supply through histone methyltransferase Set2-dependent MBF-induced expression of Cdc22, the catalytic subunit of ribonucleotide reductase (RNR). Disrupting dNTP synthesis following Wee1 inactivation, through abrogating Set2-dependent H3K36 tri-methylation or DNA integrity checkpoint inactivation results in critically low dNTP levels, replication collapse and cell death, which can be rescued by increasing dNTP levels. These findings support a 'dNTP supply and demand' model in which maintaining dNTP homeostasis is essential to prevent replication catastrophe in response to CDK-induced replication stress.

The African Swine Fever Virus Transcriptome.

African swine fever virus (ASFV) causes hemorrhagic fever in domestic pigs, presenting the biggest global threat to animal farming in recorded history. Despite the importance of ASFV, little is known about the mechanisms and regulation of ASFV transcription. Using RNA sequencing methods, we have determined total RNA abundance, transcription start sites, and transcription termination sites at single-nucleotide resolution. This allowed us to characterize DNA consensus motifs of early and late ASFV core promoters, as well as a polythymidylate sequence determinant for transcription termination. Our results demonstrate that ASFV utilizes alternative transcription start sites between early and late stages of infection and that ASFV RNA polymerase (RNAP) undergoes promoter-proximal transcript slippage at 5' ends of transcription units, adding quasitemplated AU- and AUAU-5' extensions to mRNAs. Here, we present the first much-needed genome-wide transcriptome study that provides unique insight into ASFV transcription and serves as a resource to aid future functional analyses of ASFV genes which are essential to combat this devastating disease.IMPORTANCE African swine fever virus (ASFV) causes incurable and often lethal hemorrhagic fever in domestic pigs. In 2020, ASF presents an acute and global animal health emergency that has the potential to devastate entire national economies as effective vaccines or antiviral drugs are not currently available (according to the Food and Agriculture Organization of the United Nations). With major outbreaks ongoing in Eastern Europe and Asia, urgent action is needed to advance our knowledge about the fundamental biology of ASFV, including the mechanisms and temporal control of gene expression. A thorough understanding of RNAP and transcription factor function, and of the sequence context of their promoter motifs, as well as accurate knowledge of which genes are expressed when and the amino acid sequence of the encoded proteins, is direly needed for the development of antiviral drugs and vaccines.

Pyruvate kinase variant of fission yeast tunes carbon metabolism, cell regulation, growth and stress resistance.

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.

PomBase 2018: user-driven reimplementation of the fission yeast database provides rapid and intuitive access to diverse, interconnected information.

PomBase (www.pombase.org), the model organism database for the fission yeast Schizosaccharomyces pombe, has undergone a complete redevelopment, resulting in a more fully integrated, better-performing service. The new infrastructure supports daily data updates as well as fast, efficient querying and smoother navigation within and between pages. New pages for publications and genotypes provide routes to all data curated from a single source and to all phenotypes associated with a specific genotype, respectively. For ontology-based annotations, improved displays balance comprehensive data coverage with ease of use. The default view now uses ontology structure to provide a concise, non-redundant summary that can be expanded to reveal underlying details and metadata. The phenotype annotation display also offers filtering options to allow users to focus on specific areas of interest. An instance of the JBrowse genome browser has been integrated, facilitating loading of and intuitive access to, genome-scale datasets. Taken together, the new data and pages, along with improvements in annotation display and querying, allow users to probe connections among different types of data to form a comprehensive view of fission yeast biology. The new PomBase implementation also provides a rich set of modular, reusable tools that can be deployed to create new, or enhance existing, organism-specific databases.

General amino acid control in fission yeast is regulated by a nonconserved transcription factor, with functions analogous to Gcn4/Atf4.

Eukaryotes respond to amino acid starvation by enhancing the translation of mRNAs encoding b-ZIP family transcription factors (GCN4 in Saccharomyces cerevisiae and ATF4 in mammals), which launch transcriptional programs to counter this stress. This pathway involves phosphorylation of the eIF2 translation factor by Gcn2-protein kinases and is regulated by upstream ORFs (uORFs) in the GCN4/ATF4 5' leaders. Here, we present evidence that the transcription factors that mediate this response are not evolutionarily conserved. Although cells of the fission yeast Schizosaccharomyces pombe respond transcriptionally to amino acid starvation, they lack clear Gcn4 and Atf4 orthologs. We used ribosome profiling to identify mediators of this response in S. pombe, looking for transcription factors that behave like GCN4 We discovered a transcription factor (Fil1) translationally induced by amino acid starvation in a 5' leader and Gcn2-dependent manner. Like Gcn4, Fil1 is required for the transcriptional response to amino acid starvation, and Gcn4 and Fil1 regulate similar genes. Despite their similarities in regulation, function, and targets, Fil1 and Gcn4 belong to different transcription factor families (GATA and b-ZIP, respectively). Thus, the same functions are performed by nonorthologous proteins under similar regulation. These results highlight the plasticity of transcriptional networks, which maintain conserved principles with nonconserved regulators.

Fitness Landscape of the Fission Yeast Genome.

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.

Long noncoding RNA repertoire and targeting by nuclear exosome, cytoplasmic exonuclease, and RNAi in fission yeast.

Long noncoding RNAs (lncRNAs), which are longer than 200 nucleotides but often unstable, contribute a substantial and diverse portion to pervasive noncoding transcriptomes. Most lncRNAs are poorly annotated and understood, although several play important roles in gene regulation and diseases. Here we systematically uncover and analyze lncRNAs in Schizosaccharomyces pombe. Based on RNA-seq data from twelve RNA-processing mutants and nine physiological conditions, we identify 5775 novel lncRNAs, nearly 4× the previously annotated lncRNAs. The expression of most lncRNAs becomes strongly induced under the genetic and physiological perturbations, most notably during late meiosis. Most lncRNAs are cryptic and suppressed by three RNA-processing pathways: the nuclear exosome, cytoplasmic exonuclease, and RNAi. Double-mutant analyses reveal substantial coordination and redundancy among these pathways. We classify lncRNAs by their dominant pathway into cryptic unstable transcripts (CUTs), Xrn1-sensitive unstable transcripts (XUTs), and Dicer-sensitive unstable transcripts (DUTs). XUTs and DUTs are enriched for antisense lncRNAs, while CUTs are often bidirectional and actively translated. The cytoplasmic exonuclease, along with RNAi, dampens the expression of thousands of lncRNAs and mRNAs that become induced during meiosis. Antisense lncRNA expression mostly negatively correlates with sense mRNA expression in the physiological, but not the genetic conditions. Intergenic and bidirectional lncRNAs emerge from nucleosome-depleted regions, upstream of positioned nucleosomes. Our results highlight both similarities and differences to lncRNA regulation in budding yeast. This broad survey of the lncRNA repertoire and characteristics in S. pombe, and the interwoven regulatory pathways that target lncRNAs, provides a rich framework for their further functional analyses.

A Pre-mRNA degradation pathway that selectively targets intron-containing genes requires the nuclear poly(A)-binding protein.

General discard pathways eliminate unprocessed and irregular pre-mRNAs to control the quality of gene expression. In contrast to such general pre-mRNA decay, we describe here a nuclear pre-mRNA degradation pathway that controls the expression of select intron-containing genes. We show that the fission yeast nuclear poly(A)-binding protein, Pab2, and the nuclear exosome subunit, Rrp6, are the main factors involved in this polyadenylation-dependent pre-mRNA degradation pathway. Transcriptome analysis and intron swapping experiments revealed that inefficient splicing is important to dictate susceptibility to Pab2-dependent pre-mRNA decay. We also show that negative splicing regulation can promote the poor splicing efficiency required for this pre-mRNA decay pathway, and in doing so, we identified a mechanism of cross-regulation between paralogous ribosomal proteins through nuclear pre-mRNA decay. Our findings unveil a layer of regulation in the nucleus in which the turnover of specific pre-mRNAs, besides the turnover of mature mRNAs, is used to control gene expression.

Elimination of a specific histone H3K14 acetyltransferase complex bypasses the RNAi pathway to regulate pericentric heterochromatin functions.

In Schizosaccharomyces pombe, the RNAi pathway is required for the formation of pericentric heterochromatin, proper chromosome segregation, and repression of pericentric meiotic recombination. Here we demonstrate that, when the activity of the histone H3 Lys 14 (H3K14) acetyltransferase Mst2 is eliminated, the RNAi machinery is no longer required for pericentric heterochromatin functions. We further reveal that reducing RNA polymerase II recruitment to pericentric regions is essential for maintaining heterochromatin in the absence of RNAi.

Cell-based screens and phenomics with fission yeast.

Next-generation sequencing approaches have considerably advanced our understanding of genome function and regulation. However, the knowledge of gene function and complex cellular processes remains a challenge and bottleneck in biological research. Phenomics is a rapidly emerging area, which seeks to rigorously characterize all phenotypes associated with genes or gene variants. Such high-throughput phenotyping under different conditions can be a potent approach toward gene function. The fission yeast Schizosaccharomyces pombe (S. pombe) is a proven eukaryotic model organism that is increasingly used for genomewide screens and phenomic assays. In this review, we highlight current large-scale, cell-based approaches used with S. pombe, including computational colony-growth measurements, genetic interaction screens, parallel profiling using barcodes, microscopy-based cell profiling, metabolomic methods and transposon mutagenesis. These diverse methods are starting to offer rich insights into the relationship between genotypes and phenotypes.