Transcript-specific determinants of pre-mRNA splicing revealed through in vivo kinetic analyses of the 1st and 2nd chemical steps.

We generate high-precision measurements of the in vivo rates of both chemical steps of pre-mRNA splicing across the genome-wide complement of substrates in yeast by coupling metabolic labeling, multiplexed primer-extension sequencing, and kinetic modeling. We demonstrate that the rates of intron removal vary widely, splice-site sequences are primary determinants of 1st step but have little apparent impact on 2nd step rates, and the 2nd step is generally faster than the 1st step. Ribosomal protein genes (RPGs) are spliced faster than non-RPGs at each step, and RPGs share evolutionarily conserved properties that may contribute to their faster splicing. A genetic variant defective in the 1st step of the pathway reveals a genome-wide defect in the 1st step but an unexpected, transcript-specific change in the 2nd step. Our work demonstrates that extended co-transcriptional association is an important determinant of splicing rate, a conclusion at odds with recent claims of ultra-fast splicing.

Multiplexed primer extension sequencing: A targeted RNA-seq method that enables high-precision quantitation of mRNA splicing isoforms and rare pre-mRNA splicing intermediates.

The study of pre-mRNA splicing has been greatly aided by the advent of RNA sequencing (RNA-seq), which enables the genome-wide detection of discrete splice isoforms. Quantification of these splice isoforms requires analysis of splicing informative sequencing reads, those that unambiguously map to a single splice isoform, including exon-intron spanning alignments corresponding to retained introns, as well as exon-exon junction spanning alignments corresponding to either canonically- or alternatively-spliced isoforms. Because most RNA-seq experiments are designed to produce sequencing alignments that uniformly cover the entirety of transcripts, only a comparatively small number of splicing informative alignments are generated for any given splice site, leading to a decreased ability to detect and/or robustly quantify many splice isoforms. To address this problem, we have recently described a method termed Multiplexed Primer Extension sequencing, or MPE-seq, which uses pools of reverse transcription primers to target sequencing to user selected loci. By targeting reverse transcription to pre-mRNA splice junctions, this approach enables a dramatic enrichment in the fraction of splicing informative alignments generated per splicing event, yielding an increase in both the precision with which splicing efficiency can be measured, and in the detection of splice isoforms including rare splicing intermediates. Here we provide a brief review of the shortcomings associated with RNA-seq that drove our development of MPE-seq, as well as a detailed protocol for implementation of MPE-seq.

Detection of splice isoforms and rare intermediates using multiplexed primer extension sequencing.

Targeted RNA sequencing (RNA-seq) aims to focus coverage on areas of interest that are inadequately sampled in standard RNA-seq experiments. Here we present multiplexed primer extension sequencing (MPE-seq), an approach for targeted RNA-seq that uses complex pools of reverse-transcription primers to enable sequencing enrichment at user-selected locations across the genome. We targeted hundreds to thousands of pre-mRNA splice junctions and obtained high-precision detection of splice isoforms, including rare pre-mRNA splicing intermediates.

Sde2 is an intron-specific pre-mRNA splicing regulator activated by ubiquitin-like processing.

The expression of intron-containing genes in eukaryotes requires generation of protein-coding messenger RNAs (mRNAs) via RNA splicing, whereby the spliceosome removes non-coding introns from pre-mRNAs and joins exons. Spliceosomes must ensure accurate removal of highly diverse introns. We show that Sde2 is a ubiquitin-fold-containing splicing regulator that supports splicing of selected pre-mRNAs in an intron-specific manner in Schizosaccharomyces pombe Both fission yeast and human Sde2 are translated as inactive precursor proteins harbouring the ubiquitin-fold domain linked through an invariant GGKGG motif to a C-terminal domain (referred to as Sde2-C). Precursor processing after the first di-glycine motif by the ubiquitin-specific proteases Ubp5 and Ubp15 generates a short-lived activated Sde2-C fragment with an N-terminal lysine residue, which subsequently gets incorporated into spliceosomes. Absence of Sde2 or defects in Sde2 activation both result in inefficient excision of selected introns from a subset of pre-mRNAs. Sde2 facilitates spliceosomal association of Cactin/Cay1, with a functional link between Sde2 and Cactin further supported by genetic interactions and pre-mRNA splicing assays. These findings suggest that ubiquitin-like processing of Sde2 into a short-lived activated form may function as a checkpoint to ensure proper splicing of certain pre-mRNAs in fission yeast.

The conserved AU dinucleotide at the 5' end of nascent U1 snRNA is optimized for the interaction with nuclear cap-binding-complex.

Splicing is initiated by a productive interaction between the pre-mRNA and the U1 snRNP, in which a short RNA duplex is established between the 5' splice site of a pre-mRNA and the 5' end of the U1 snRNA. A long-standing puzzle has been why the AU dincucleotide at the 5'-end of the U1 snRNA is highly conserved, despite the absence of an apparent role in the formation of the duplex. To explore this conundrum, we varied this AU dinucleotide into all possible permutations and analyzed the resulting molecular consequences. This led to the unexpected findings that the AU dinucleotide dictates the optimal binding of cap-binding complex (CBC) to the 5' end of the nascent U1 snRNA, which ultimately influences the utilization of U1 snRNP in splicing. Our data also provide a structural interpretation as to why the AU dinucleotide is conserved during evolution.

Structural toggle in the RNaseH domain of Prp8 helps balance splicing fidelity and catalytic efficiency.

Pre-mRNA splicing is an essential step of eukaryotic gene expression that requires both high efficiency and high fidelity. Prp8 has long been considered the "master regulator" of the spliceosome, the molecular machine that executes pre-mRNA splicing. Cross-linking and structural studies place the RNaseH domain (RH) of Prp8 near the spliceosome's catalytic core and demonstrate that prp8 alleles that map to a 17-aa extension in RH stabilize it in one of two mutually exclusive structures, the biological relevance of which are unknown. We performed an extensive characterization of prp8 alleles that map to this extension and, using in vitro and in vivo reporter assays, show they fall into two functional classes associated with the two structures: those that promote error-prone/efficient splicing and those that promote hyperaccurate/inefficient splicing. Identification of global locations of endogenous splice-site activation by lariat sequencing confirms the fidelity effects seen in our reporter assays. Furthermore, we show that error-prone/efficient RH alleles suppress a prp2 mutant deficient at promoting the first catalytic step of splicing, whereas hyperaccurate/inefficient RH alleles exhibit synthetic sickness. Together our data indicate that prp8 RH alleles link splicing fidelity with catalytic efficiency by biasing the relative stabilities of distinct spliceosome conformations. We hypothesize that the spliceosome "toggles" between such error-prone/efficient and hyperaccurate/inefficient conformations during the splicing cycle to regulate splicing fidelity.

The power of fission: yeast as a tool for understanding complex splicing.

Pre-mRNA splicing is an essential component of eukaryotic gene expression. Many metazoans, including humans, regulate alternative splicing patterns to generate expansions of their proteome from a limited number of genes. Importantly, a considerable fraction of human disease causing mutations manifest themselves through altering the sequences that shape the splicing patterns of genes. Thus, understanding the mechanistic bases of this complex pathway will be an essential component of combating these diseases. Dating almost to the initial discovery of splicing, researchers have taken advantage of the genetic tractability of budding yeast to identify the components and decipher the mechanisms of splicing. However, budding yeast lacks the complex splicing machinery and alternative splicing patterns most relevant to humans. More recently, many researchers have turned their efforts to study the fission yeast, Schizosaccharomyces pombe, which has retained many features of complex splicing, including degenerate splice site sequences, the usage of exonic splicing enhancers, and SR proteins. Here, we review recent work using fission yeast genetics to examine pre-mRNA splicing, highlighting its promise for modeling the complex splicing seen in higher eukaryotes.

Phospho-site mutants of the RNA Polymerase II C-terminal domain alter subtelomeric gene expression and chromatin modification state in fission yeast.

Eukaryotic gene expression requires that RNA Polymerase II (RNAP II) gain access to DNA in the context of chromatin. The C-terminal domain (CTD) of RNAP II recruits chromatin modifying enzymes to promoters, allowing for transcription initiation or repression. Specific CTD phosphorylation marks facilitate recruitment of chromatin modifiers, transcriptional regulators, and RNA processing factors during the transcription cycle. However, the readable code for recruiting such factors is still not fully defined and how CTD modifications affect related families of genes or regional gene expression is not well understood. Here, we examine the effects of manipulating the Y1S2P3T4S5P6S7 heptapeptide repeat of the CTD of RNAP II in Schizosaccharomyces pombe by substituting non-phosphorylatable alanines for Ser2 and/or Ser7 and the phosphomimetic glutamic acid for Ser7. Global gene expression analyses were conducted using splicing-sensitive microarrays and validated via RT-qPCR. The CTD mutations did not affect pre-mRNA splicing or snRNA levels. Rather, the data revealed upregulation of subtelomeric genes and alteration of the repressive histone H3 lysine 9 methylation (H3K9me) landscape. The data further indicate that H3K9me and expression status are not fully correlated, suggestive of CTD-dependent subtelomeric repression mechansims that act independently of H3K9me levels.

Interconnections Between RNA-Processing Pathways Revealed by a Sequencing-Based Genetic Screen for Pre-mRNA Splicing Mutants in Fission Yeast.

Pre-mRNA splicing is an essential component of eukaryotic gene expression and is highly conserved from unicellular yeasts to humans. Here, we present the development and implementation of a sequencing-based reverse genetic screen designed to identify nonessential genes that impact pre-mRNA splicing in the fission yeast Schizosaccharomyces pombe, an organism that shares many of the complex features of splicing in higher eukaryotes. Using a custom-designed barcoding scheme, we simultaneously queried ∼3000 mutant strains for their impact on the splicing efficiency of two endogenous pre-mRNAs. A total of 61 nonessential genes were identified whose deletions resulted in defects in pre-mRNA splicing; enriched among these were factors encoding known or predicted components of the spliceosome. Included among the candidates identified here are genes with well-characterized roles in other RNA-processing pathways, including heterochromatic silencing and 3' end processing. Splicing-sensitive microarrays confirm broad splicing defects for many of these factors, revealing novel functional connections between these pathways.

Mating-Induced Transcriptome Changes in the Reproductive Tract of Female Aedes aegypti.

The Aedes aegypti mosquito is a significant public health threat, as it is the main vector of dengue and chikungunya viruses. Disease control efforts could be enhanced through reproductive manipulation of these vectors. Previous work has revealed a relationship between male seminal fluid proteins transferred to females during mating and female post-mating physiology and behavior. To better understand this interplay, we used short-read RNA sequencing to identify gene expression changes in the lower reproductive tract of females in response to mating. We characterized mRNA expression in virgin and mated females at 0, 6 and 24 hours post-mating (hpm) and identified 364 differentially abundant transcripts between mating status groups. Surprisingly, 60 transcripts were more abundant at 0 hpm compared to virgin females, suggesting transfer from males. Twenty of these encode known Ae. aegypti seminal fluid proteins. Transfer and detection of male accessory gland-derived mRNA in females at 0 hpm was confirmed by measurement of eGFP mRNA in females mated to eGFP-expressing males. In addition, 150 transcripts were up-regulated at 6 hpm and 24 hpm, while 130 transcripts were down-regulated at 6 hpm and 24 hpm. Gene Ontology (GO) enrichment analysis revealed that proteases, a protein class broadly known to play important roles in reproduction, were among the most enriched protein classes. RNAs associated with immune system and antimicrobial function were also up-regulated at 24 hpm. Collectively, our results suggest that copulation initiates broad transcriptome changes across the mosquito female reproductive tract, "priming" her for important subsequent processes of blood feeding, egg development and immune defense. Our transcriptome analysis provides a vital foundation for future studies of the consequences of mating on female biology and will aid studies seeking to identify specific gene families, molecules and pathways that support key reproductive processes in the female mosquito.

A Proteome-wide Fission Yeast Interactome Reveals Network Evolution Principles from Yeasts to Human.

Here, we present FissionNet, a proteome-wide binary protein interactome for S. pombe, comprising 2,278 high-quality interactions, of which ∼ 50% were previously not reported in any species. FissionNet unravels previously unreported interactions implicated in processes such as gene silencing and pre-mRNA splicing. We developed a rigorous network comparison framework that accounts for assay sensitivity and specificity, revealing extensive species-specific network rewiring between fission yeast, budding yeast, and human. Surprisingly, although genes are better conserved between the yeasts, S. pombe interactions are significantly better conserved in human than in S. cerevisiae. Our framework also reveals that different modes of gene duplication influence the extent to which paralogous proteins are functionally repurposed. Finally, cross-species interactome mapping demonstrates that coevolution of interacting proteins is remarkably prevalent, a result with important implications for studying human disease in model organisms. Overall, FissionNet is a valuable resource for understanding protein functions and their evolution.

Widespread alternative and aberrant splicing revealed by lariat sequencing.

Alternative splicing is an important and ancient feature of eukaryotic gene structure, the existence of which has likely facilitated eukaryotic proteome expansions. Here, we have used intron lariat sequencing to generate a comprehensive profile of splicing events in Schizosaccharomyces pombe, amongst the simplest organisms that possess mammalian-like splice site degeneracy. We reveal an unprecedented level of alternative splicing, including alternative splice site selection for over half of all annotated introns, hundreds of novel exon-skipping events, and thousands of novel introns. Moreover, the frequency of these events is far higher than previous estimates, with alternative splice sites on average activated at ∼3% the rate of canonical sites. Although a subset of alternative sites are conserved in related species, implying functional potential, the majority are not detectably conserved. Interestingly, the rate of aberrant splicing is inversely related to expression level, with lowly expressed genes more prone to erroneous splicing. Although we validate many events with RNAseq, the proportion of alternative splicing discovered with lariat sequencing is far greater, a difference we attribute to preferential decay of aberrantly spliced transcripts. Together, these data suggest the spliceosome possesses far lower fidelity than previously appreciated, highlighting the potential contributions of alternative splicing in generating novel gene structures.

Mdm1/Snx13 is a novel ER-endolysosomal interorganelle tethering protein.

Although endolysosomal trafficking is well defined, how it is regulated and coordinates with cellular metabolism is unclear. To identify genes governing endolysosomal dynamics, we conducted a global fluorescence-based screen to reveal endomembrane effector genes. Screening implicated Phox (PX) domain-containing protein Mdm1 in endomembrane dynamics. Surprisingly, we demonstrate that Mdm1 is a novel interorganelle tethering protein that localizes to endoplasmic reticulum (ER)-vacuole/lysosome membrane contact sites (MCSs). We show that Mdm1 is ER anchored and contacts the vacuole surface in trans via its lipid-binding PX domain. Strikingly, overexpression of Mdm1 induced ER-vacuole hypertethering, underscoring its role as an interorganelle tether. We also show that Mdm1 and its paralogue Ydr179w-a (named Nvj3 in this study) localize to ER-vacuole MCSs independently of established tether Nvj1. Finally, we find that Mdm1 truncations analogous to neurological disease-associated SNX14 alleles fail to tether the ER and vacuole and perturb sphingolipid metabolism. Our work suggests that human Mdm1 homologues may play previously unappreciated roles in interorganelle communication and lipid metabolism.

Cytoplasmic polyadenylation is a major mRNA regulator during oogenesis and egg activation in Drosophila.

The GLD-2 class of poly(A) polymerases regulate the timing of translation of stored transcripts by elongating the poly(A) tails of target mRNAs in the cytoplasm. WISPY is a GLD-2 enzyme that acts in the Drosophila female germline and is required for the completion of the egg-to-embryo transition. Though a handful of WISPY target mRNAs have been identified during both oogenesis and early embryogenesis, it was unknown whether WISP simply regulated a small pool of patterning or cell cycle genes, or whether, instead, cytoplasmic polyadenylation was widespread during this developmental transition. To identify the full range of WISPY targets, we carried out microarray analysis to look for maternal mRNAs whose poly(A) tails fail to elongate in the absence of WISP function. We examined the polyadenylated portion of the maternal transcriptome in both stage 14 (mature) oocytes and in early embryos that had completed egg activation. Our analysis shows that the poly(A) tails of thousands of maternal mRNAs fail to elongate in wisp-deficient oocytes and embryos. Furthermore, we have identified specific classes of genes that are highly regulated in this manner at each stage. Our study shows that cytoplasmic polyadenylation is a major regulatory mechanism during oocyte maturation and egg activation.

Lariat sequencing in a unicellular yeast identifies regulated alternative splicing of exons that are evolutionarily conserved with humans.

Alternative splicing is a potent regulator of gene expression that vastly increases proteomic diversity in multicellular eukaryotes and is associated with organismal complexity. Although alternative splicing is widespread in vertebrates, little is known about the evolutionary origins of this process, in part because of the absence of phylogenetically conserved events that cross major eukaryotic clades. Here we describe a lariat-sequencing approach, which offers high sensitivity for detecting splicing events, and its application to the unicellular fungus, Schizosaccharomyces pombe, an organism that shares many of the hallmarks of alternative splicing in mammalian systems but for which no previous examples of exon-skipping had been demonstrated. Over 200 previously unannotated splicing events were identified, including examples of regulated alternative splicing. Remarkably, an evolutionary analysis of four of the exons identified here as subject to skipping in S. pombe reveals high sequence conservation and perfect length conservation with their homologs in scores of plants, animals, and fungi. Moreover, alternative splicing of two of these exons have been documented in multiple vertebrate organisms, making these the first demonstrations of identical alternative-splicing patterns in species that are separated by over 1 billion y of evolution.

Cross-species protein interactome mapping reveals species-specific wiring of stress response pathways.

The fission yeast Schizosaccharomyces pombe has more metazoan-like features than the budding yeast Saccharomyces cerevisiae, yet it has similarly facile genetics. We present a large-scale verified binary protein-protein interactome network, "StressNet," based on high-throughput yeast two-hybrid screens of interacting proteins classified as part of stress response and signal transduction pathways in S. pombe. We performed systematic, cross-species interactome mapping using StressNet and a protein interactome network of orthologous proteins in S. cerevisiae. With cross-species comparative network studies, we detected a previously unidentified component (Snr1) of the S. pombe mitogen-activated protein kinase Sty1 pathway. Coimmunoprecipitation experiments showed that Snr1 interacted with Sty1 and that deletion of snr1 increased the sensitivity of S. pombe cells to stress. Comparison of StressNet with the interactome network of orthologous proteins in S. cerevisiae showed that most of the interactions among these stress response and signaling proteins are not conserved between species but are "rewired"; orthologous proteins have different binding partners in both species. In particular, transient interactions connecting proteins in different functional modules were more likely to be rewired than conserved. By directly testing interactions between proteins in one yeast species and their corresponding binding partners in the other yeast species with yeast two-hybrid assays, we found that about half of the interactions that are traditionally considered "conserved" form modified interaction interfaces that may potentially accommodate novel functions.

Systematic dissection of roles for chromatin regulators in a yeast stress response.

Packaging of eukaryotic genomes into chromatin has wide-ranging effects on gene transcription. Curiously, it is commonly observed that deletion of a global chromatin regulator affects expression of only a limited subset of genes bound to or modified by the regulator in question. However, in many single-gene studies it has become clear that chromatin regulators often do not affect steady-state transcription, but instead are required for normal transcriptional reprogramming by environmental cues. We therefore have systematically investigated the effects of 83 histone mutants, and 119 gene deletion mutants, on induction/repression dynamics of 170 transcripts in response to diamide stress in yeast. Importantly, we find that chromatin regulators play far more pronounced roles during gene induction/repression than they do in steady-state expression. Furthermore, by jointly analyzing the substrates (histone mutants) and enzymes (chromatin modifier deletions) we identify specific interactions between histone modifications and their regulators. Combining these functional results with genome-wide mapping of several histone marks in the same time course, we systematically investigated the correspondence between histone modification occurrence and function. We followed up on one pathway, finding that Set1-dependent H3K4 methylation primarily acts as a gene repressor during multiple stresses, specifically at genes involved in ribosome biosynthesis. Set1-dependent repression of ribosomal genes occurs via distinct pathways for ribosomal protein genes and ribosomal biogenesis genes, which can be separated based on genetic requirements for repression and based on chromatin changes during gene repression. Together, our dynamic studies provide a rich resource for investigating chromatin regulation, and identify a significant role for the "activating" mark H3K4me3 in gene repression.

A quantitative, high-throughput reverse genetic screen reveals novel connections between Pre-mRNA splicing and 5' and 3' end transcript determinants.

Here we present the development and implementation of a genome-wide reverse genetic screen in the budding yeast, Saccharomyces cerevisiae, that couples high-throughput strain growth, robotic RNA isolation and cDNA synthesis, and quantitative PCR to allow for a robust determination of the level of nearly any cellular RNA in the background of ~5,500 different mutants. As an initial test of this approach, we sought to identify the full complement of factors that impact pre-mRNA splicing. Increasing lines of evidence suggest a relationship between pre-mRNA splicing and other cellular pathways including chromatin remodeling, transcription, and 3' end processing, yet in many cases the specific proteins responsible for functionally connecting these pathways remain unclear. Moreover, it is unclear whether all pathways that are coupled to splicing have been identified. As expected, our approach sensitively detects pre-mRNA accumulation in the vast majority of strains containing mutations in known splicing factors. Remarkably, however, several additional candidates were found to cause increases in pre-mRNA levels similar to that seen for canonical splicing mutants, none of which had previously been implicated in the splicing pathway. Instead, several of these factors have been previously implicated to play roles in chromatin remodeling, 3' end processing, and other novel categories. Further analysis of these factors using splicing-sensitive microarrays confirms that deletion of Bdf1, a factor that links transcription initiation and chromatin remodeling, leads to a global splicing defect, providing evidence for a novel connection between pre-mRNA splicing and this component of the SWR1 complex. By contrast, mutations in 3' end processing factors such as Cft2 and Yth1 also result in pre-mRNA splicing defects, although only for a subset of transcripts, suggesting that spliceosome assembly in S. cerevisiae may more closely resemble mammalian models of exon-definition. More broadly, our work demonstrates the capacity of this approach to identify novel regulators of various cellular RNAs.