SummaryData 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.

Author contributions DCJ coordinated all analyses, isolated DNA for sequencing, analysed and filtered SNP calls, conducted diversity analysis and GWAS and drafted the manuscript. CR produced phenotype data for growth on various solid media and growth rates in liquid media. AR conducted analysis of dating using mitochondrial data. DS conducted GWAS. MP analysed all phenotype data. TM identified LTR transposon insertions and analysed transposon insertion data. FXM conducted crosses for analysis of spore viability ZI produced indel calls with Cortex. WL conducted analysis of recombination rate, linkage disequilibrium decay and PCA for distance between strains. TMKC assisted with phenotype and population analysis. RP analysed Cortex and GATK indel calls. MM conducted amino acid profiling. JLDL and AC produced automated measures of cell morphology. SB aligned reads and produced GATK SNP calls. GH analysed population structure using fineSTRUCTURE. BO'F estimated the TMRCA from the nuclear genome using ACG. TK identified LTR transposon insertions JTS produced de novo assemblies. LB developed the custom Workspace workflow Spotsizer. BT assisted with sequence analysis. DAB assisted with analysis of novel genes. TS assisted with strain verification. SC produced images of wild strains and assisted with strain verification. JEEUH assisted with SNP validation. LvT and MT assisted with LTR validation. LJ and JL assisted with manual measures of cell morphology and FACS. SA produced gene expression data. MF, KM and ND assisted with sequencing. WB initiated and assisted with strain collection. JH coordinated manual measures of cell morphology and FACS. RECS coordinated automated measures of cell morphology. MR coordinated amino acid profiling. NM conducted analysis of recombination, linkage disequilibrium and advised on aspects of diversity and GWAS. DJB advised on GWAS. RD facilitated sequencing. JB contributed to the initiation and development of the project and financed the JB laboratory. AccessionsSequence data are archived in the European Nucleotide Archive (www.ebi.ac.uk/ena/), Study Accessions PRJEB2733 and PRJEB6284 (Supplementary Table 7). All SNPs and indels were submitted to NCBI dbSNP (www.ncbi.nlm.nih.gov/SNP/). Accessions are 974514578-974688138 (SNPs) and 974702618-974688139 (indels). Europe PMC Funders Group AbstractNatural variation within species reveals aspects of genome evolution and function. The fission yeast Schizosaccharomyces pombe is an important model for eukaryotic biology, but researchers typically use one standard laboratory strain. To extend the utility of this model, we surveyed the genomic and phenotypic variation in 161 natural isolates. We sequenced the genomes of all strains, revealing moderate genetic diversity (π = 3 ×10 −3 ) and weak global population structure. We estimate that dispersal of S. pombe began within human antiquity (~340 BCE), and ancestors of these strains reached the Americas at ~1623 CE. We quantified 74 traits, revealing substantial heritable phenotypic diversity. We cond...

Target of rapamycin complex 1 (TORC1) is implicated in growth control and aging from yeast to humans. Fission yeast is emerging as a popular model organism to study TOR signaling, although rapamycin has been thought to not affect cell growth in this organism. Here, we analyzed the effects of rapamycin and caffeine, singly and combined, on multiple cellular processes in fission yeast. The two drugs led to diverse and specific phenotypes that depended on TORC1 inhibition, including prolonged chronological lifespan, inhibition of global translation, inhibition of cell growth and division, and reprograming of global gene expression mimicking nitrogen starvation. Rapamycin and caffeine differentially affected these various TORC1-dependent processes. Combined drug treatment augmented most phenotypes and effectively blocked cell growth. Rapamycin showed a much more subtle effect on global translation than did caffeine, while both drugs were effective in prolonging chronological lifespan. Rapamycin and caffeine did not affect the lifespan via the pH of the growth media. Rapamycin prolonged the lifespan of nongrowing cells only when applied during the growth phase but not when applied after cells had stopped proliferation. The doses of rapamycin and caffeine strongly correlated with growth inhibition and with lifespan extension. This comprehensive analysis will inform future studies into TORC1 function and cellular aging in fission yeast and beyond.

DNA double-strand break (DSB) repair is a highly regulated process performed predominantly by non-homologous end joining (NHEJ) or homologous recombination (HR) pathways. How these pathways are coordinated in the context of chromatin is unclear. Here we uncover a role for histone H3K36 modification in regulating DSB repair pathway choice in fission yeast. We find Set2-dependent H3K36 methylation reduces chromatin accessibility, reduces resection and promotes NHEJ, while antagonistic Gcn5-dependent H3K36 acetylation increases chromatin accessibility, increases resection and promotes HR. Accordingly, loss of Set2 increases H3K36Ac, chromatin accessibility and resection, while Gcn5 loss results in the opposite phenotypes following DSB induction. Further, H3K36 modification is cell cycle regulated with Set2-dependent H3K36 methylation peaking in G1 when NHEJ occurs, while Gcn5-dependent H3K36 acetylation peaks in S/G2 when HR prevails. These findings support an H3K36 chromatin switch in regulating DSB repair pathway choice.

No abstract

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.

Both canonical and alternative splicing of RNAs are governed by intronic sequence elements and produce transient lariat structures fastened by branch points within introns. To map precisely the location of branch points on a genomic scale, we developed LaSSO (Lariat Sequence Site Origin), a data-driven algorithm which utilizes RNA-seq data. Using fission yeast cells lacking the debranching enzyme Dbr1, LaSSO not only accurately identified canonical splicing events, but also pinpointed novel, but rare, exon-skipping events, which may reflect aberrantly spliced transcripts. Compromised intron turnover perturbed gene regulation at multiple levels, including splicing and protein translation. Notably, Dbr1 function was also critical for the expression of mitochondrial genes and for the processing of self-spliced mitochondrial introns. LaSSO showed better sensitivity and accuracy than algorithms used for computational branch-point prediction or for empirical branch-point determination. Even when applied to a human data set acquired in the presence of debranching activity, LaSSO identified both canonical and exon-skipping branch points. LaSSO thus provides an effective approach for defining high-resolution maps of branch-site sequences and intronic elements on a genomic scale. LaSSO should be useful to validate introns and uncover branch-point sequences in any eukaryote, and it could be integrated into RNA-seq pipelines. [Supplemental material is available for this article.]Introns and exons refer to noncoding and coding sequences, respectively, that constitute protein-coding genes (Gilbert 1978). To create a functional messenger RNA (mRNA), introns are excised via a highly conserved and accurate process called splicing that culminates in concatenation of exon sequences into translatable transcripts. Splicing entails two transesterification reactions catalyzed by the spliceosome, a large RNA-protein complex (Wahl et al. 2009). The first reaction involves a nucleophilic attack of an adenosine (branch point) on the 59-splice donor, resulting in a lariat structure fixed by a 29-59 phosphodiester bond; the intron remains only attached to the downstream exon (Fig. 1A,B ;Padgett et al. 1985). The second reaction involves an attack of the detached upstream exon on the 39-splice acceptor, resulting in intron lariat release and exon ligation (Fig. 1C). The intron is then processed by exonucleolytic cleavage of the 39-lariat tail and linearization by the debranching enzyme Dbr1 (Fig.

No abstract