Danny A. Bitton

The genomic and phenotypic diversity of Schizosaccharomyces pombe

Natural 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 usefulness of this model, we surveyed the genomic and phenotypic variation in 161 natural isolates. We sequenced the genomes of all strains, finding moderate genetic diversity (π = 3 × 10−3 substitutions/site) and weak global population structure. We estimate that dispersal of S. pombe began during human antiquity (∼340 BCE), and ancestors of these strains reached the Americas at ∼1623 CE. We quantified 74 traits, finding substantial heritable phenotypic diversity. We conducted 223 genome-wide association studies, with 89 traits showing at least one association. The most significant variant for each trait explained 22% of the phenotypic variance on average, with indels having larger effects than SNPs. This analysis represents a rich resource to examine genotype-phenotype relationships in a tractable model.

LaSSO, a strategy for genome-wide mapping of intronic lariats and branch points using RNA-seq

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.

Individual letters of the RNA polymerase II CTD code govern distinct gene expression programs in fission yeast.

Significance The carboxyl-terminal domain (CTD) code encrypted within the YSPTSPS heptad repeats of RNA polymerase II is deeply rooted in eukaryal biology. This paper exploits fission yeast CTD mutants in which single coding cues (or combinations of cues) are subtracted from all consensus heptads to gauge the roles of the Tyr1, Ser2, Thr4, and Ser7 phosphoacceptors in gene expression. Transcriptome profiling revealed distinctive signatures for the CTD phosphoacceptors that were verified by cell-based activity assays. The results highlight (i) the role of Tyr1 and Ser2 in repressing the iron uptake regulon and (ii) the distinctive roles of Thr4 and Ser7 in the phosphate starvation response, whereby Thr4 facilitates and Ser7 represses the induction of Pho1 acid phosphatase. The primary structure and phosphorylation pattern of the tandem Y1S2P3T4S5P6S7 repeats of the RNA polymerase II carboxyl-terminal domain (CTD) comprise an informational code that coordinates transcription, chromatin modification, and RNA processing. To gauge the contributions of individual CTD coding “letters” to gene expression, we analyzed the poly(A)+ transcriptomes of fission yeast mutants that lack each of the four inessential CTD phosphoacceptors: Tyr1, Ser2, Thr4, and Ser7. There was a hierarchy of CTD mutational effects with respect to the number of dysregulated protein-coding RNAs, with S2A (n = 227) >> Y1F (n = 71) > S7A (n = 58) >> T4A (n = 7). The majority of the protein-coding RNAs affected in Y1F cells were coordinately affected by S2A, suggesting that Tyr1-Ser2 constitutes a two-letter code “word.” Y1F and S2A elicited increased expression of genes encoding proteins involved in iron uptake (Frp1, Fip1, Fio1, Str3, Str1, Sib1), without affecting the expression of the genes that repress the iron regulon, implying that Tyr1-Ser2 transduces a repressive signal. Y1F and S2A cells had increased levels of ferric reductase activity and were hypersensitive to phleomycin, indicative of elevated intracellular iron. The T4A and S7A mutations had opposing effects on the phosphate response pathway. T4A reduced the expression of two genes encoding proteins involved in phosphate acquisition (the Pho1 acid phosphatase and the phosphate transporter SPBC8E4.01c), without affecting the expression of known genes that regulate the phosphate response pathway, whereas S7A increased pho1+ expression. These results highlight specific cellular gene expression programs that are responsive to distinct CTD cues.

AnGeLi: A Tool for the Analysis of Gene Lists from Fission Yeast.

Genome-wide assays and screens typically result in large lists of genes or proteins. Enrichments of functional or other biological properties within such lists can provide valuable insights and testable hypotheses. To systematically detect these enrichments can be challenging and time-consuming, because relevant data to compare against query gene lists are spread over many different sources. We have developed AnGeLi (Analysis of Gene Lists), an intuitive, integrated web-tool for comprehensive and customized interrogation of gene lists from the fission yeast, Schizosaccharomyces pombe. AnGeLi searches for significant enrichments among multiple qualitative and quantitative information sources, including gene and phenotype ontologies, genetic and protein interactions, numerous features of genes, transcripts, translation, and proteins such as copy numbers, chromosomal positions, genetic diversity, RNA polymerase II and ribosome occupancy, localization, conservation, half-lives, domains, and molecular weight among others, as well as diverse sets of genes that are co-regulated or lead to the same phenotypes when mutated. AnGeLi uses robust statistics which can be tailored to specific needs. It also provides the option to upload user-defined gene sets to compare against the query list. Through an integrated data submission form, AnGeLi encourages the community to contribute additional curated gene lists to further increase the usefulness of this resource and to get the most from the ever increasing large-scale experiments. AnGeLi offers a rigorous yet flexible statistical analysis platform for rich insights into functional enrichments and biological context for query gene lists, thus providing a powerful exploratory tool through which S. pombe researchers can uncover fresh perspectives and unexpected connections from genomic data. AnGeLi is freely available at: www.bahlerlab.info/AnGeLi

Parallel Profiling of Fission Yeast Deletion Mutants for Proliferation and for Lifespan During Long-Term Quiescence

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.

Structural and Functional Characterization of the N Terminus of Schizosaccharomyces pombe Cwf10

ABSTRACT The spliceosome is a dynamic macromolecular machine that catalyzes the removal of introns from pre-mRNA, yielding mature message. Schizosaccharomyces pombe Cwf10 (homolog of Saccharomyces cerevisiae Snu114 and human U5-116K), an integral member of the U5 snRNP, is a GTPase that has multiple roles within the splicing cycle. Cwf10/Snu114 family members are highly homologous to eukaryotic translation elongation factor EF2, and they contain a conserved N-terminal extension (NTE) to the EF2-like portion, predicted to be an intrinsically unfolded domain. Using S. pombe as a model system, we show that the NTE is not essential, but cells lacking this domain are defective in pre-mRNA splicing. Genetic interactions between cwf10-ΔNTE and other pre-mRNA splicing mutants are consistent with a role for the NTE in spliceosome activation and second-step catalysis. Characterization of Cwf10-NTE by various biophysical techniques shows that in solution the NTE contains regions of both structure and disorder. The first 23 highly conserved amino acids of the NTE are essential for its role in splicing but when overexpressed are not sufficient to restore pre-mRNA splicing to wild-type levels in cwf10-ΔNTE cells. When the entire NTE is overexpressed in the cwf10-ΔNTE background, it can complement the truncated Cwf10 protein in trans, and it immunoprecipitates a complex similar in composition to the late-stage U5.U2/U6 spliceosome. These data show that the structurally flexible NTE is capable of independently incorporating into the spliceosome and improving splicing function, possibly indicating a role for the NTE in stabilizing conformational rearrangements during a splice cycle.

Identification of new players in cell division, DNA damage response, and morphogenesis through construction of Schizosaccharomyces pombe deletion strains.

Many fundamental biological processes are studied using the fission yeast, Schizosaccharomyces pombe. Here we report the construction of a set of 281 haploid gene deletion strains covering many previously uncharacterized genes. This collection of strains was tested for growth under a variety of different stress conditions. We identified new genes involved in DNA metabolism, completion of the cell cycle, and morphogenesis. This subset of nonessential gene deletions will add to the toolkits available for the study of biological processes in S. pombe.

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.

Functional and Regulatory Profiling of Energy Metabolism in Fission Yeast

BackgroundThe control of energy metabolism is fundamental for cell growth and function and anomalies in it are implicated in complex diseases and ageing. Metabolism in yeast cells can be manipulated by supplying different carbon sources: yeast grown on glucose rapidly proliferates by fermentation, analogous to tumour cells growing by aerobic glycolysis, whereas on non-fermentable carbon sources metabolism shifts towards respiration.ResultsWe screened deletion libraries of fission yeast to identify over 200 genes required for respiratory growth. Growth media and auxotrophic mutants strongly influenced respiratory metabolism. Most genes uncovered in the mutant screens have not been implicated in respiration in budding yeast. We applied gene-expression profiling approaches to compare steady-state fermentative and respiratory growth and to analyse the dynamic adaptation to respiratory growth. The transcript levels of most genes functioning in energy metabolism pathways are coherently tuned, reflecting anticipated differences in metabolic flows between fermenting and respiring cells. We show that acetyl-CoA synthase, rather than citrate lyase, is essential for acetyl-CoA synthesis in fission yeast. We also investigated the transcriptional response to mitochondrial damage by genetic or chemical perturbations, defining a retrograde response that involves the concerted regulation of distinct groups of nuclear genes that may avert harm from mitochondrial malfunction.ConclusionsThis study provides a rich framework of the genetic and regulatory basis of energy metabolism in fission yeast and beyond, and it pinpoints weaknesses of commonly used auxotroph mutants for investigating metabolism. As a model for cellular energy regulation, fission yeast provides an attractive and complementary system to budding yeast.

Fitness Landscape of the Fission Yeast Genome

Background Non-protein-coding regions of eukaryotic genomes remain poorly understood. Diversity studies, comparative genomics and biochemical outputs of genomic sites can be indicators of functional elements, but none produce fine-scale genome-wide descriptions of all functional elements. Results Towards the generation of a comprehensive description of functional elements in the haploid Schizosaccharomyces pombe genome, we generated transposon mutagenesis libraries to a density of one insertion per 13 nucleotides of the genome. We applied a five-state hidden Markov model (HMM) to characterise insertion-depleted regions at nucleotide-level resolution. HMM-defined functional constraint was consistent with genetic diversity, comparative genomics, gene-expression data and genome annotation. Conclusions We infer that transposon insertions lead to fitness consequences in 90% of the genome, including 80% of the non-protein-coding regions, reflecting the presence of numerous non-coding elements in this compact genome that have functional roles. Display of this data in genome browsers provides fine-scale views of structure-function relationships within specific genes.