Determinants of RNA metabolism in the Schizosaccharomyces pombe genome

Co-transcriptional gene regulation in eukaryotes and prokaryotes

Many steps of RNA processing occur during transcription by RNA polymerases. Co-transcriptional activities are deemed commonplace in prokaryotes, in which the lack of membrane barriers allows mixing of all gene expression steps, from transcription to translation. In the past decade, an extraordinary level of coordination between transcription and RNA processing has emerged in eukaryotes. In this Review, we discuss recent developments in our understanding of co-transcriptional gene regulation in both eukaryotes and prokaryotes, comparing methodologies and mechanisms, and highlight striking parallels in how RNA polymerases interact with the machineries that act on nascent RNA. The development of RNA sequencing and imaging techniques that detect transient transcription and RNA processing intermediates has facilitated discoveries of transcription coordination with splicing, 3'-end cleavage and dynamic RNA folding and revealed physical contacts between processing machineries and RNA polymerases. Such studies indicate that intron retention in a given nascent transcript can prevent 3'-end cleavage and cause transcriptional readthrough, which is a hallmark of eukaryotic cellular stress responses. We also discuss how coordination between nascent RNA biogenesis and transcription drives fundamental aspects of gene expression in both prokaryotes and eukaryotes.

Timing is everything: advances in quantifying splicing kinetics

Splicing is a highly regulated process critical for proper pre-mRNA maturation and the maintenance of a healthy cellular environment. Splicing events are impacted by ongoing transcription, neighboring splicing events, and cis and trans regulatory factors on the respective pre-mRNA transcript. Within this complex regulatory environment, splicing kinetics have the potential to influence splicing outcomes but have historically been challenging to study in vivo. In this review, we highlight recent technological advancements that have enabled measurements of global splicing kinetics and of the variability of splicing kinetics at single introns. We demonstrate how identifying features that are correlated with splicing kinetics has increased our ability to form potential models for how splicing kinetics may be regulated in vivo.

Post-transcriptional splicing can occur in a slow-moving zone around the gene

Splicing is the stepwise molecular process by which introns are removed from pre-mRNA and exons are joined together to form mature mRNA sequences. The ordering and spatial distribution of these steps remain controversial, with opposing models suggesting splicing occurs either during or after transcription. We used single-molecule RNA FISH, expansion microscopy, and live-cell imaging to reveal the spatiotemporal distribution of nascent transcripts in mammalian cells. At super-resolution levels, we found that pre-mRNA formed clouds around the transcription site. These clouds indicate the existence of a transcription-site-proximal zone through which RNA move more slowly than in the nucleoplasm. Full-length pre-mRNA undergo continuous splicing as they move through this zone following transcription, suggesting a model in which splicing can occur post-transcriptionally but still within the proximity of the transcription site, thus seeming co-transcriptional by most assays. These results may unify conflicting reports of co-transcriptional versus post-transcriptional splicing.

Species-aware DNA language models capture regulatory elements and their evolution

BACKGROUND

The rise of large-scale multi-species genome sequencing projects promises to shed new light on how genomes encode gene regulatory instructions. To this end, new algorithms are needed that can leverage conservation to capture regulatory elements while accounting for their evolution.

RESULTS

Here, we introduce species-aware DNA language models, which we trained on more than 800 species spanning over 500 million years of evolution. Investigating their ability to predict masked nucleotides from context, we show that DNA language models distinguish transcription factor and RNA-binding protein motifs from background non-coding sequence. Owing to their flexibility, DNA language models capture conserved regulatory elements over much further evolutionary distances than sequence alignment would allow. Remarkably, DNA language models reconstruct motif instances bound in vivo better than unbound ones and account for the evolution of motif sequences and their positional constraints, showing that these models capture functional high-order sequence and evolutionary context. We further show that species-aware training yields improved sequence representations for endogenous and MPRA-based gene expression prediction, as well as motif discovery.

CONCLUSIONS

Collectively, these results demonstrate that species-aware DNA language models are a powerful, flexible, and scalable tool to integrate information from large compendia of highly diverged genomes.

Simultaneous measurement of nascent transcriptome and translatome using 4-thiouridine metabolic RNA labeling and translating ribosome affinity purification

Regulation of gene expression in response to various biological processes, including extracellular stimulation and environmental adaptation requires nascent RNA synthesis and translation. Analysis of the coordinated regulation of dynamic RNA synthesis and translation is required to determine functional protein production. However, reliable methods for the simultaneous measurement of nascent RNA synthesis and translation at the gene level are limited. Here, we developed a novel method for the simultaneous assessment of nascent RNA synthesis and translation by combining 4-thiouridine (4sU) metabolic RNA labeling and translating ribosome affinity purification (TRAP) using a monoclonal antibody against evolutionarily conserved ribosomal P-stalk proteins. The P-stalk-mediated TRAP (P-TRAP) technique recovered endogenous translating ribosomes, allowing easy translatome analysis of various eukaryotes. We validated this method in mammalian cells by demonstrating that acute unfolded protein response (UPR) in the endoplasmic reticulum (ER) induces dynamic reprogramming of nascent RNA synthesis and translation. Our nascent P-TRAP (nP-TRAP) method may serve as a simple and powerful tool for analyzing the coordinated regulation of transcription and translation of individual genes in various eukaryotes.

RNA-Mediated Regulation of Meiosis in Budding Yeast

Cells change their physiological state in response to environmental cues. In the absence of nutrients, unicellular fungi such as budding yeast exit mitotic proliferation and enter the meiotic cycle, leading to the production of haploid cells that are encased within spore walls. These cell state transitions are orchestrated in a developmentally coordinated manner. Execution of the meiotic cell cycle program in budding yeast, Saccharomyces cerevisiae, is regulated by the key transcription factor, Ime1. Recent developments have uncovered the role of non-coding RNA in the regulation of Ime1 and meiosis. In this review, we summarize the role of ncRNA-mediated and RNA homeostasis-based processes in the regulation of meiosis in Saccharomyces cerevisiae.

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.

Mitotic checkpoint gene expression is tuned by codon usage bias

The mitotic checkpoint (also called spindle assembly checkpoint, SAC) is a signaling pathway that safeguards proper chromosome segregation. Correct functioning of the SAC depends on adequate protein concentrations and appropriate stoichiometries between SAC proteins. Yet very little is known about the regulation of SAC gene expression. Here, we show in the fission yeast Schizosaccharomyces pombe that a combination of short mRNA half-lives and long protein half-lives supports stable SAC protein levels. For the SAC genes mad2⁺ and mad3⁺ , their short mRNA half-lives are caused, in part, by a high frequency of nonoptimal codons. In contrast, mad1⁺ mRNA has a short half-life despite a higher frequency of optimal codons, and despite the lack of known RNA-destabilizing motifs. Hence, different SAC genes employ different strategies of expression. We further show that Mad1 homodimers form co-translationally, which may necessitate a certain codon usage pattern. Taken together, we propose that the codon usage of SAC genes is fine-tuned to ensure proper SAC function. Our work shines light on gene expression features that promote spindle assembly checkpoint function and suggests that synonymous mutations may weaken the checkpoint.

Genomic sequences and RNA binding proteins predict RNA splicing efficiency in various single-cell contexts

MOTIVATION

The RNA splicing efficiency is of high interest for both understanding the regulatory machinery of gene expression and estimating the RNA velocity in single cells. However, its genomic regulation and stochasticity across contexts remain poorly understood.

RESULTS

Here, by leveraging the recent RNA velocity tool, we estimated the relative splicing efficiency across a variety of single-cell RNA-Seq data sets. We further extracted large sets of genomic features and 120 RNA binding protein features and found they are highly predictive to relative RNA splicing efficiency across multiple tissues and organs on human and mouse. This predictive power brings promise to reveal the complexity of RNA processing and to enhance the analysis of single-cell transcription activities.

AVAILABILITY AND IMPLEMENTATION

In order to ensure reproducibility, all preprocessed data sets and scripts used for the prediction and figure generation are publicly available at https://doi.org/10.5281/zenodo.6513669.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Cdc48 influence on separase levels is independent of mitosis and suggests translational sensitivity of separase

Cdc48 (p97/VCP) is a AAA-ATPase that can extract ubiquitinated proteins from their binding partners and can cooperate with the proteasome for their degradation. A fission yeast cdc48 mutant (cdc48-353) shows low levels of the cohesin protease, separase, and pronounced chromosome segregation defects in mitosis. Separase initiates chromosome segregation when its binding partner securin is ubiquitinated and degraded. The low separase levels in the cdc48-353 mutant have been attributed to a failure to extract ubiquitinated securin from separase, resulting in co-degradation of separase along with securin. If true, Cdc48 would be important in mitosis. In contrast, we show here that low separase levels in the cdc48-353 mutant are independent of mitosis. Moreover, we find no evidence of enhanced separase degradation in the mutant. Instead, we suggest that the cdc48-353 mutant uncovers specific requirements for separase translation. Our results highlight a need to better understand how this key mitotic enzyme is synthesized.

Functional profiling of long intergenic non-coding RNAs in fission yeast

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.

RNA-binding protein Mub1 and the nuclear RNA exosome act to fine-tune environmental stress response

Comparative RNA interactome capture identifies potential regulators of RNA metabolism in fission yeast and reveals RNA exosome–dependent buffering of stress-responsive gene expression networks.

The nuclear RNA exosome plays a key role in controlling the levels of multiple protein-coding and non-coding RNAs. Recruitment of the exosome to specific RNA substrates is mediated by RNA-binding co-factors. The transient interaction between co-factors and the exosome as well as the rapid decay of RNA substrates make identification of exosome co-factors challenging. Here, we use comparative poly(A)+ RNA interactome capture in fission yeast expressing three different mutants of the exosome to identify proteins that interact with poly(A)+ RNA in an exosome-dependent manner. Our analyses identify multiple RNA-binding proteins whose association with RNA is altered in exosome mutants, including the zinc-finger protein Mub1. Mub1 is required to maintain the levels of a subset of exosome RNA substrates including mRNAs encoding for stress-responsive proteins. Removal of the zinc-finger domain leads to loss of RNA suppression under non-stressed conditions, altered expression of heat shock genes in response to stress, and reduced growth at elevated temperature. These findings highlight the importance of exosome-dependent mRNA degradation in buffering gene expression networks to mediate cellular adaptation to stress.

Transcription and chromatin-based surveillance mechanism controls suppression of cryptic antisense transcription

Phosphorylation of the RNA polymerase II C-terminal domain Y₁S₂P₃T₄S₅P₆S₇ consensus sequence coordinates key events during transcription, and its deregulation leads to defects in transcription and RNA processing. Here, we report that the histone deacetylase activity of the fission yeast Hos2/Set3 complex plays an important role in suppressing cryptic initiation of antisense transcription when RNA polymerase II phosphorylation is dysregulated due to the loss of Ssu72 phosphatase. Interestingly, although single Hos2 and Set3 mutants have little effect, loss of Hos2 or Set3 combined with ssu72Δ results in a synergistic increase in antisense transcription globally and correlates with elevated sensitivity to genotoxic agents. We demonstrate a key role for the Ssu72/Hos2/Set3 mechanism in the suppression of cryptic antisense transcription at the 3' end of convergent genes that are most susceptible to these defects, ensuring the fidelity of gene expression within dense genomes of simple eukaryotes.

Role of promoters in regulating alternative splicing

Alternative splicing (AS) plays a critical role in enhancing proteome complexity in higher eukaryotes. Almost all the multi intron-containing genes undergo AS in humans. Splicing mainly occurs co-transcriptionally, where RNA polymerase II (RNA pol II) plays a crucial role in coordinating transcription and pre-mRNA splicing. Aberrant AS leads to non-functional proteins causative in various pathophysiological conditions such as cancers, neurodegenerative diseases, and muscular dystrophies. Transcription and pre-mRNA splicing are deeply interconnected and can influence each other's functions. Several studies evinced that specific promoters employed by RNA pol II dictate the RNA processing decisions. Promoter-specific recruitment of certain transcriptional factors or transcriptional coactivators influences splicing, and the extent to which these factors affect splicing has not been discussed in detail. Here, in this review, various DNA-binding proteins and their influence on promoter-specific AS are extensively discussed. Besides, this review highlights how the promoter-specific epigenetic changes might regulate AS.

A Trans-Omics Comparison Reveals Common Gene Expression Strategies in Four Model Organisms and Exposes Similarities and Differences between Them

The ultimate goal of gene expression regulation is on the protein level. However, because the amounts of mRNAs and proteins are controlled by their synthesis and degradation rates, the cellular amount of a given protein can be attained by following different strategies. By studying omics data for six expression variables (mRNA and protein amounts, plus their synthesis and decay rates), we previously demonstrated the existence of common expression strategies (CESs) for functionally related genes in the yeast Saccharomyces cerevisiae. Here we extend that study to two other eukaryotes: the yeast Schizosaccharomyces pombe and cultured human HeLa cells. We also use genomic data from the model prokaryote Escherichia coli as an external reference. We show that six-variable profiles (6VPs) can be constructed for every gene and that these 6VPs are similar for genes with similar functions in all the studied organisms. The differences in 6VPs between organisms can be used to establish their phylogenetic relationships. The analysis of the correlations among the six variables supports the hypothesis that most gene expression control occurs in actively growing organisms at the transcription rate level, and that translation plays a minor role. We propose that living organisms use CESs for the genes acting on the same physiological pathways, especially for those belonging to stable macromolecular complexes, but CESs have been modeled by evolution to adapt to the specific life circumstances of each organism.

Co-transcriptional splicing regulates 3' end cleavage during mammalian erythropoiesis

Pre-mRNA processing steps are tightly coordinated with transcription in many organisms. To determine how co-transcriptional splicing is integrated with transcription elongation and 3' end formation in mammalian cells, we performed long-read sequencing of individual nascent RNAs and precision run-on sequencing (PRO-seq) during mouse erythropoiesis. Splicing was not accompanied by transcriptional pausing and was detected when RNA polymerase II (Pol II) was within 75-300 nucleotides of 3' splice sites (3'SSs), often during transcription of the downstream exon. Interestingly, several hundred introns displayed abundant splicing intermediates, suggesting that splicing delays can take place between the two catalytic steps. Overall, splicing efficiencies were correlated among introns within the same transcript, and intron retention was associated with inefficient 3' end cleavage. Remarkably, a thalassemia patient-derived mutation introducing a cryptic 3'SS improved both splicing and 3' end cleavage of individual β-globin transcripts, demonstrating functional coupling between the two co-transcriptional processes as a determinant of productive gene output.

Genome-scale molecular principles of mRNA half-life regulation in yeast

Precise control of protein and messenger RNA (mRNA) degradation is essential for cellular metabolism and homeostasis. Controlled and specific degradation of both molecular species necessitates their engagements with the respective degradation machineries; this engagement involves a disordered/unstructured segment of the substrate traversing the degradation tunnel of the machinery and accessing the catalytic sites. However, while molecular factors influencing protein degradation have been extensively explored on a genome scale, and in multiple organisms, such a comprehensive understanding remains missing for mRNAs. Here, we analyzed multiple genome-scale experimental yeast mRNA half-life data in light of experimentally derived mRNA secondary structures and protein binding data, along with high-resolution X-ray crystallographic structures of the RNase machines. Results unraveled a consistent genome-scale trend that mRNAs comprising longer terminal and/or internal unstructured segments have significantly shorter half-lives; the lengths of the 5'-terminal, 3'-terminal, and internal unstructured segments that affect mRNA half-life are compatible with molecular structures of the 5' exo-, 3' exo-, and endoribonuclease machineries. Sequestration into ribonucleoprotein complexes elongates mRNA half-life, presumably by burying ribonuclease engagement sites under oligomeric interfaces. After gene duplication, differences in terminal unstructured lengths, proportions of internal unstructured segments, and oligomerization modes result in significantly altered half-lives of paralogous mRNAs. Side-by-side comparison of molecular principles underlying controlled protein and mRNA degradation in yeast unravels their remarkable mechanistic similarities and suggests how the intrinsic structural features of the two molecular species, at two different levels of the central dogma, regulate their half-lives on genome scale.

Metabolic labeling of RNA using multiple ribonucleoside analogs enables the simultaneous evaluation of RNA synthesis and degradation rates

Gene expression is determined by a balance between RNA synthesis and RNA degradation. To elucidate the underlying regulatory mechanisms and principles of this, simultaneous measurements of RNA synthesis and degradation are required. Here, we report the development of “Dyrec-seq,” which uses 4-thiouridine and 5-bromouridine to simultaneously quantify RNA synthesis and degradation rates. Dyrec-seq enabled the quantification of RNA synthesis and degradation rates of 4702 genes in HeLa cells. Functional enrichment analysis showed that the RNA synthesis and degradation rates of genes are actually determined by the genes’ biological functions. A comparison of theoretical and experimental analyses revealed that the amount of RNA is determined by the ratio of RNA synthesis to degradation rates, whereas the rapidity of responses to external stimuli is determined only by the degradation rate. This study emphasizes that not only RNA synthesis but also RNA degradation is important in shaping gene expression patterns.

Comprehensive profiling of the fission yeast transcription start site activity during stress and media response

Fission yeast, Schizosaccharomyces pombe, is an attractive model organism for transcriptional and chromatin biology research. Such research is contingent on accurate annotation of transcription start sites (TSSs). However, comprehensive genome-wide maps of TSSs and their usage across commonly applied laboratory conditions and treatments for S. pombe are lacking. To this end, we profiled TSS activity genome-wide in S. pombe cultures exposed to heat shock, nitrogen starvation, hydrogen peroxide and two commonly applied media, YES and EMM2, using Cap Analysis of Gene Expression (CAGE). CAGE-based annotation of TSSs is substantially more accurate than existing PomBase annotation; on average, CAGE TSSs fall 50–75 bp downstream of PomBase TSSs and co-localize with nucleosome boundaries. In contrast to higher eukaryotes, dispersed TSS distributions are not common in S. pombe. Our data recapitulate known S. pombe stress expression response patterns and identify stress- and media-responsive alternative TSSs. Notably, alteration of growth medium induces changes of similar magnitude as some stressors. We show a link between nucleosome occupancy and genetic variation, and that the proximal promoter region is genetically diverse between S. pombe strains. Our detailed TSS map constitutes a central resource for S. pombe gene regulation research.

On the optimal design of metabolic RNA labeling experiments

Massively parallel RNA sequencing (RNA-seq) in combination with metabolic labeling has become the de facto standard approach to study alterations in RNA transcription, processing or decay. Regardless of advances in the experimental protocols and techniques, every experimentalist needs to specify the key aspects of experimental design: For example, which protocol should be used (biochemical separation vs. nucleotide conversion) and what is the optimal labeling time? In this work, we provide approximate answers to these questions using the asymptotic theory of optimal design. Specifically, we investigate, how the variance of degradation rate estimates depends on the time and derive the optimal time for any given degradation rate. Subsequently, we show that an increase in sample numbers should be preferred over an increase in sequencing depth. Lastly, we provide some guidance on use cases when laborious biochemical separation outcompetes recent nucleotide conversion based methods (such as SLAMseq) and show, how inefficient conversion influences the precision of estimates. Code and documentation can be found at https://github.com/dieterich-lab/DesignMetabolicRNAlabeling.

Massively parallel RNA sequencing (RNA-seq) in combination with metabolic labeling has become the de facto standard approach to study alterations in RNA transcription, processing or decay. In our manuscript, we address several key aspects of experimental design: 1) The optimal labeling time, 2) the number of replicate samples over sequencing depth and 3) the choice of experimental protocol. We provide approximate answers to these questions using asymptotic theory of optimal design.

Nascent RNA and the Coordination of Splicing with Transcription

At each active protein-encoding gene, nascent RNA is tethered to the DNA axis by elongating RNA polymerase II (Pol II) and is continuously altered by splicing and other processing events during its synthesis. This review discusses the development of three major methods that enable us to track the conversion of precursor messenger RNA (pre-mRNA) to messenger RNA (mRNA) products in vivo: live-cell imaging, metabolic labeling of RNA, and RNA-seq of purified nascent RNA. These approaches are complementary, addressing distinct issues of transcription rates and intron lifetimes alongside spatial information regarding the gene position of Pol II at which spliceosomes act. The findings will be placed in the context of active transcription units, each of which-because of the presence of nascent RNA, Pol II, and features of the chromatin environment-will recruit a potentially gene-specific constellation of RNA binding proteins and processing machineries.

Global donor and acceptor splicing site kinetics in human cells

RNA splicing is an essential part of eukaryotic gene expression. Although the mechanism of splicing has been extensively studied in vitro, in vivo kinetics for the two-step splicing reaction remain poorly understood. Here, we combine transient transcriptome sequencing (TT-seq) and mathematical modeling to quantify RNA metabolic rates at donor and acceptor splice sites across the human genome. Splicing occurs in the range of minutes and is limited by the speed of RNA polymerase elongation. Splicing kinetics strongly depends on the position and nature of nucleotides flanking splice sites, and on structural interactions between unspliced RNA and small nuclear RNAs in spliceosomal intermediates. Finally, we introduce the 'yield' of splicing as the efficiency of converting unspliced to spliced RNA and show that it is highest for mRNAs and independent of splicing kinetics. These results lead to quantitative models describing how splicing rates and yield are encoded in the human genome.

RNA-Seq approach for accurate characterization of splicing efficiency of yeast introns

Introns in different genes, or even different introns within the same gene, often have different splice sites and differ in splicing efficiency (SE). One expects mass-transcribed genes to have introns with higher SE than weakly transcribed genes. However, such a simple expectation cannot be tested directly because variable SE for these genes is often not measured. Mechanistically, SE should depend on signal strength at key splice sites (SS) such as 5'SS, 3'SS and branchpoint site (BPS), i.e., SE = F(5'SS, 3'SS, BPS). However, without SE, we again cannot model how these splice sites contribute to SE. Here I present an RNA-Seq approach to quantify SE for each of the 304 introns in yeast (Saccharomyces cerevisiae) genes, including 24 in the 5'UTR, by measuring 1) number of reads mapped to exon-exon junctions (N_EE) as a proxy for the abundance of spliced form, and 2) number of reads mapped to exon-intron junction (N_EI5 and N_EI3 at 5' and 3' ends of intron) as a proxy for the abundance of unspliced form. The total mRNA is N_Total = N_EE + p * N_EI5 + (1-p) * N_EI3, with the simplest p = 0.5 but statistical methods were presented to estimate p from data. An estimated p is needed because N_EI5 is expected to be smaller than N_EI3 due to 1) step 1 splicing occurs before step 2 so EI5 is broken before EI3, 2) enrichment of poly(A) mRNA by oligo-dT, and 3) 5' degradation. SE is defined as the proportion (N_EE/N_Total). Application of the method shows that ribosomal protein messages are efficiently and mostly cotranscriptionally spliced. Yeast genes with long introns are also spliced efficiently. HAC1/YFL031W is poorly spliced partly because its splicing involves a nonspliceosome mechanism and partly because Ire1p, which participate in splicing HAC1, is hardly expressed. Many putative yeast genes have low SE, and some splice sites are incorrectly annotated.

Quantification and discovery of sequence determinants of protein-per-mRNA amount in 29 human tissues

Despite their importance in determining protein abundance, a comprehensive catalogue of sequence features controlling protein-to-mRNA (PTR) ratios and a quantification of their effects are still lacking. Here, we quantified PTR ratios for 11,575 proteins across 29 human tissues using matched transcriptomes and proteomes. We estimated by regression the contribution of known sequence determinants of protein synthesis and degradation in addition to 45 mRNA and 3 protein sequence motifs that we found by association testing. While PTR ratios span more than 2 orders of magnitude, our integrative model predicts PTR ratios at a median precision of 3.2-fold. A reporter assay provided functional support for two novel UTR motifs, and an immobilized mRNA affinity competition-binding assay identified motif-specific bound proteins for one motif. Moreover, our integrative model led to a new metric of codon optimality that captures the effects of codon frequency on protein synthesis and degradation. Altogether, this study shows that a large fraction of PTR ratio variation in human tissues can be predicted from sequence, and it identifies many new candidate post-transcriptional regulatory elements.

Saccharomyces cerevisiae Metabolic Labeling with 4-thiouracil and the Quantification of Newly Synthesized mRNA As a Proxy for RNA Polymerase II Activity

Global defects in RNA polymerase II transcription might be overlooked by transcriptomic studies analyzing steady-state RNA. Indeed, the global decrease in mRNA synthesis has been shown to be compensated by a simultaneous decrease in mRNA degradation to restore normal steady-state levels. Hence, the genome-wide quantification of mRNA synthesis, independently from mRNA decay, is the best direct reflection of RNA polymerase II transcriptional activity. Here, we discuss a method using non-perturbing metabolic labeling of nascent RNAs in Saccharomyces cerevisiae (S. cerevisiae). Specifically, the cells are cultured for 6 min with a uracil analog, 4-thiouracil, and the labeled newly transcribed RNAs are purified and quantified to determine the synthesis rates of all individual mRNA. Moreover, using labeled Schizosaccharomyces pombe cells as internal standard allows comparing mRNA synthesis in different S. cerevisiae strains. Using this protocol and fitting the data with a dynamic kinetic model, the corresponding mRNA decay rates can be determined.

Elongation/Termination Factor Exchange Mediated by PP1 Phosphatase Orchestrates Transcription Termination

Termination of RNA polymerase II (Pol II) transcription is a key step that is important for 3' end formation of functional mRNA, mRNA release, and Pol II recycling. Even so, the underlying termination mechanism is not yet understood. Here, we demonstrate that the conserved and essential termination factor Seb1 is found on Pol II near the end of the RNA exit channel and the Rpb4/7 stalk. Furthermore, the Seb1 interaction surface with Pol II largely overlaps with that of the elongation factor Spt5. Notably, Seb1 co-transcriptional recruitment is dependent on Spt5 dephosphorylation by the conserved PP1 phosphatase Dis2, which also dephosphorylates threonine 4 within the Pol II heptad repeated C-terminal domain. We propose that Dis2 orchestrates the transition from elongation to termination phase during the transcription cycle by mediating elongation to termination factor exchange and dephosphorylation of Pol II C-terminal domain.

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.

Schizosaccharomyces pombe cardiolipin synthase is part of a mitochondrial fusion protein regulated by intron retention

Cardiolipin (CL) is a unique lipid component of mitochondria in all eukaryotes. It is important for the architecture of mitochondrial membranes and for mitochondrial dynamics. CL also creates a highly specific microenvironment of mitochondrial protein machineries. CL biosynthetic pathway is, however, only partially characterized in the fission yeast Schizosaccharomyces pombe. Here we show that CL synthase is an essential protein in S. pombe. It is encoded by the ORF SPAC22A12.08c as a C terminal part of a tandem fusion protein together with a mitochondrial hydrolase of unknown function. Expression of S. pombe CL synthase is able to complement deletion of the CRD1 gene of Saccharomyces cerevisiae and, vice versa, S. cerevisiae CRD1 gene complements deletion of S. pombe SPAC22A12.08c. The proper expression of CL synthase and its partner in the tandem protein, the mitochondrial hydrolase, is regulated at the level of alternate intron splicing. The first part of the SPAC22A12.08c fusion protein could be translated from both major SPAC22A12.08c derived mRNAs, with and without intron IV. Functional CL synthase, however, is produced only from the minor SPAC22A12.08c derived mRNA that has intron IV retained. Thus, intron retention is a novel mechanism for the differential expression of two proteins that evolved as a fusion protein and are under the control of the same promoter.

Histone deacetylation promotes transcriptional silencing at facultative heterochromatin

It is important to accurately regulate the expression of genes involved in development and environmental response. In the fission yeast Schizosaccharomyces pombe, meiotic genes are tightly repressed during vegetative growth. Despite being embedded in heterochromatin these genes are transcribed and believed to be repressed primarily at the level of RNA. However, the mechanism of facultative heterochromatin formation and the interplay with transcription regulation is not understood. We show genome-wide that HDAC-dependent histone deacetylation is a major determinant in transcriptional silencing of facultative heterochromatin domains. Indeed, mutation of class I/II HDACs leads to increased transcription of meiotic genes and accumulation of their mRNAs. Mechanistic dissection of the pho1 gene where, in response to phosphate, transient facultative heterochromatin is established by overlapping lncRNA transcription shows that the Clr3 HDAC contributes to silencing independently of SHREC, but in an lncRNA-dependent manner. We propose that HDACs promote facultative heterochromatin by establishing alternative transcriptional silencing.

Long-read sequencing of nascent RNA reveals coupling among RNA processing events

Pre-mRNA splicing is accomplished by the spliceosome, a megadalton complex that assembles de novo on each intron. Because spliceosome assembly and catalysis occur cotranscriptionally, we hypothesized that introns are removed in the order of their transcription in genomes dominated by constitutive splicing. Remarkably little is known about splicing order and the regulatory potential of nascent transcript remodeling by splicing, due to the limitations of existing methods that focus on analysis of mature splicing products (mRNAs) rather than substrates and intermediates. Here, we overcome this obstacle through long-read RNA sequencing of nascent, multi-intron transcripts in the fission yeast Schizosaccharomyces pombe. Most multi-intron transcripts were fully spliced, consistent with rapid cotranscriptional splicing. However, an unexpectedly high proportion of transcripts were either fully spliced or fully unspliced, suggesting that splicing of any given intron is dependent on the splicing status of other introns in the transcript. Supporting this, mild inhibition of splicing by a temperature-sensitive mutation in prp2, the homolog of vertebrate U2AF65, increased the frequency of fully unspliced transcripts. Importantly, fully unspliced transcripts displayed transcriptional read-through at the polyA site and were degraded cotranscriptionally by the nuclear exosome. Finally, we show that cellular mRNA levels were reduced in genes with a high number of unspliced nascent transcripts during caffeine treatment, showing regulatory significance of cotranscriptional splicing. Therefore, overall splicing of individual nascent transcripts, 3′ end formation, and mRNA half-life depend on the splicing status of neighboring introns, suggesting crosstalk among spliceosomes and the polyA cleavage machinery during transcription elongation.

Landscape of the spliced leader trans-splicing mechanism in Schistosoma mansoni

Spliced leader dependent trans-splicing (SLTS) has been described as an important RNA regulatory process that occurs in different organisms, including the trematode Schistosoma mansoni. We identified more than seven thousand putative SLTS sites in the parasite, comprising genes with a wide spectrum of functional classes, which underlines the SLTS as a ubiquitous mechanism in the parasite. Also, SLTS gene expression levels span several orders of magnitude, showing that SLTS frequency is not determined by the expression level of the target gene, but by the presence of particular gene features facilitating or hindering the trans-splicing mechanism. Our in-depth investigation of SLTS events demonstrates widespread alternative trans-splicing (ATS) acceptor sites occurring in different regions along the entire gene body, highlighting another important role of SLTS generating alternative RNA isoforms in the parasite, besides the polycistron resolution. Particularly for introns where SLTS directly competes for the same acceptor substrate with cis-splicing, we identified for the first time additional and important features that might determine the type of splicing. Our study substantially extends the current knowledge of RNA processing by SLTS in S. mansoni, and provide basis for future studies on the trans-splicing mechanism in other eukaryotes.

Circadian clock-dependent and -independent posttranscriptional regulation underlies temporal mRNA accumulation in mouse liver

Rhythms in gene expression propelled by the circadian clock and environmental signals are ubiquitous across cells and tissues. In particular, in mouse tissues, thousands of transcripts show oscillations with a period of 24 hours. Keys question are how such rhythms propagate and eventually exert functions, but also how these are generated. Here, we developed a mathematical model based on total RNA-seq to classify genes according to the respective contributions of transcriptional and posttranscriptional regulation toward mRNA expression profiles. We found that about one-third of rhythmically accumulating mRNA are under posttranscriptional regulation. Such regulation is only partially dependent on the circadian clock, showing that systemic pathways and feeding patterns contribute important posttranscriptional control of gene expression in liver.

The mammalian circadian clock coordinates physiology with environmental cycles through the regulation of daily oscillations of gene expression. Thousands of transcripts exhibit rhythmic accumulations across mouse tissues, as determined by the balance of their synthesis and degradation. While diurnally rhythmic transcription regulation is well studied and often thought to be the main factor generating rhythmic mRNA accumulation, the extent of rhythmic posttranscriptional regulation is debated, and the kinetic parameters (e.g., half-lives), as well as the underlying regulators (e.g., mRNA-binding proteins) are relatively unexplored. Here, we developed a quantitative model for cyclic accumulations of pre-mRNA and mRNA from total RNA-seq data, and applied it to mouse liver. This allowed us to identify that about 20% of mRNA rhythms were driven by rhythmic mRNA degradation, and another 15% of mRNAs regulated by both rhythmic transcription and mRNA degradation. The method could also estimate mRNA half-lives and processing times in intact mouse liver. We then showed that, depending on mRNA half-life, rhythmic mRNA degradation can either amplify or tune phases of mRNA rhythms. By comparing mRNA rhythms in wild-type and Bmal1^−/− animals, we found that the rhythmic degradation of many transcripts did not depend on a functional BMAL1. Interestingly clock-dependent and -independent degradation rhythms peaked at distinct times of day. We further predicted mRNA-binding proteins (mRBPs) that were implicated in the posttranscriptional regulation of mRNAs, either through stabilizing or destabilizing activities. Together, our results demonstrate how posttranscriptional regulation temporally shapes rhythmic mRNA accumulation in mouse liver.

Native elongating transcript sequencing reveals global anti-correlation between sense and antisense nascent transcription in fission yeast

Antisense transcription can regulate sense gene expression. However, previous annotations of antisense transcription units have been based on detection of mature antisense long noncoding (aslnc)RNAs by RNA-seq and/or microarrays, only giving a partial view of the antisense transcription landscape and incomplete molecular bases for antisense-mediated regulation. Here, we used native elongating transcript sequencing to map genome-wide nascent antisense transcription in fission yeast. Strikingly, antisense transcription was detected for most protein-coding genes, correlating with low sense transcription, especially when overlapping the mRNA start site. RNA profiling revealed that the resulting aslncRNAs mainly correspond to cryptic Xrn1/Exo2-sensitive transcripts (XUTs). ChIP-seq analyses showed that antisense (as)XUT's expression is associated with specific histone modification patterns. Finally, we showed that asXUTs are controlled by the histone chaperone Spt6 and respond to meiosis induction, in both cases anti-correlating with levels of the paired-sense mRNAs, supporting physiological significance to antisense-mediated gene attenuation. Our work highlights that antisense transcription is much more extended than anticipated and might constitute an additional nonpromoter determinant of gene regulation complexity.

The kinetics of pre-mRNA splicing in the Drosophila genome and the influence of gene architecture

Production of most eukaryotic mRNAs requires splicing of introns from pre-mRNA. The splicing reaction requires definition of splice sites, which are initially recognized in either intron-spanning (‘intron definition’) or exon-spanning (‘exon definition’) pairs. To understand how exon and intron length and splice site recognition mode impact splicing, we measured splicing rates genome-wide in Drosophila, using metabolic labeling/RNA sequencing and new mathematical models to estimate rates. We found that the modal intron length range of 60–70 nt represents a local maximum of splicing rates, but that much longer exon-defined introns are spliced even faster and more accurately. We observed unexpectedly low variation in splicing rates across introns in the same gene, suggesting the presence of gene-level influences, and we identified multiple gene level variables associated with splicing rate. Together our data suggest that developmental and stress response genes may have preferentially evolved exon definition in order to enhance the rate or accuracy of splicing.

A Massively Parallel Reporter Assay of 3' UTR Sequences Identifies In Vivo Rules for mRNA Degradation

The stability of mRNAs is regulated by signals within their sequences, but a systematic and predictive understanding of the underlying sequence rules remains elusive. Here we introduce UTR-seq, a combination of massively parallel reporter assays and regression models, to survey the dynamics of tens of thousands of 3' UTR sequences during early zebrafish embryogenesis. UTR-seq revealed two temporal degradation programs: a maternally encoded early-onset program and a late-onset program that accelerated degradation after zygotic genome activation. Three signals regulated early-onset rates: stabilizing poly-U and UUAG sequences and destabilizing GC-rich signals. Three signals explained late-onset degradation: miR-430 seeds, AU-rich sequences, and Pumilio recognition sites. Sequence-based regression models translated 3' UTRs into their unique decay patterns and predicted the in vivo effect of sequence signals on mRNA stability. Their application led to the successful design of artificial 3' UTRs that conferred specific mRNA dynamics. UTR-seq provides a general strategy to uncover the rules of RNA cis regulation.

Cis-regulatory elements explain most of the mRNA stability variation across genes in yeast

The stability of mRNA is one of the major determinants of gene expression. Although a wealth of sequence elements regulating mRNA stability has been described, their quantitative contributions to half-life are unknown. Here, we built a quantitative model for Saccharomyces cerevisiae based on functional mRNA sequence features that explains 59% of the half-life variation between genes and predicts half-life at a median relative error of 30%. The model revealed a new destabilizing 3' UTR motif, ATATTC, which we functionally validated. Codon usage proves to be the major determinant of mRNA stability. Nonetheless, single-nucleotide variations have the largest effect when occurring on 3' UTR motifs or upstream AUGs. Analyzing mRNA half-life data of 34 knockout strains showed that the effect of codon usage not only requires functional decapping and deadenylation, but also the 5'-to-3' exonuclease Xrn1, the nonsense-mediated decay genes, but not no-go decay. Altogether, this study quantitatively delineates the contributions of mRNA sequence features on stability in yeast, reveals their functional dependencies on degradation pathways, and allows accurate prediction of half-life from mRNA sequence.

Phosphorylation of the RNA-binding protein Zfs1 modulates sexual differentiation in fission yeast

Sexual differentiation in the fission yeast Schizosaccharomyces pombe promotes cell cycle arrest and extensive changes in gene expression, resulting in cell-to-cell fusion, the exchange of hereditary material and specialized cell division. These events are detrimental to the cell if they are triggered in inappropriate conditions, and therefore the decision to differentiate must be precisely controlled. Here, we investigated the role of the RNA-binding protein Zfs1 in this process by identifying its targets and characterizing novel post-translational regulatory mechanisms. We found that Zfs1 negatively regulates the G1 cyclin Puc1, and deregulated Puc1 levels inhibit differentiation in the zfs1Δ mutant. We also found that Zfs1 undergoes phosphorylation, which is stimulated upon nitrogen depletion or inhibition of the TOR pathway. Phosphorylation of Zfs1 modulates accumulation of Puc1 and plays an important role in the response of the cell to sexual differentiation signals. We propose that Zfs1 functions as an integrator of nutrient information to modulate sexual differentiation, contributing to the establishment of the differentiation-activating threshold.

Emerging Properties and Functional Consequences of Noncoding Transcription

Eukaryotic genomes are rich in transcription units encoding "long noncoding RNAs" (lncRNAs). The purpose of all this transcription is unclear since most lncRNAs are quickly targeted for destruction during synthesis or shortly thereafter. As debates continue over the functional significance of many specific lncRNAs, support grows for the notion that the act of transcription rather than the RNA product itself is functionally important in many cases. Indeed, this alternative mechanism might better explain how low-abundance lncRNAs transcribed from noncoding DNA function in organisms. Here, we highlight some of the recently emerging features that distinguish coding from noncoding transcription and discuss how these differences might have important implications for the functional consequences of noncoding transcription.

Emerging role of HuR in inflammatory response in kidney diseases

Human antigen R (HuR) is a member of the embryonic lethal abnormal vision (ELAV) family which can bind to the A/U rich elements in 3' un-translated region of mRNA and regulate mRNA splicing, transportation, and stability. Unlike other members of the ELAV family, HuR is ubiquitously expressed. Early studies mainly focused on HuR function in malignant diseases. As researches proceed, more and more proofs demonstrate its relationship with inflammation. Since most kidney diseases involve pathological changes of inflammation, HuR is now suggested to play a pivotal role in glomerular nephropathy, tubular ischemia-reperfusion damage, renal fibrosis and even renal tumors. By regulating the mRNAs of target genes, HuR is causally linked to the onset and progression of kidney diseases. Reports on this topic are steadily increasing, however, the detailed function and mechanism of action of HuR are still not well understood. The aim of this review article is to summarize the present understanding of the role of HuR in inflammation in kidney diseases, and we anticipate that future research will ultimately elucidate the therapeutic value of this novel target.

Architecture of the RNA polymerase II-Paf1C-TFIIS transcription elongation complex

The conserved polymerase-associated factor 1 complex (Paf1C) plays multiple roles in chromatin transcription and genomic regulation. Paf1C comprises the five subunits Paf1, Leo1, Ctr9, Cdc73 and Rtf1, and binds to the RNA polymerase II (Pol II) transcription elongation complex (EC). Here we report the reconstitution of Paf1C from Saccharomyces cerevisiae, and a structural analysis of Paf1C bound to a Pol II EC containing the elongation factor TFIIS. Cryo-electron microscopy and crosslinking data reveal that Paf1C is highly mobile and extends over the outer Pol II surface from the Rpb2 to the Rpb3 subunit. The Paf1-Leo1 heterodimer and Cdc73 form opposite ends of Paf1C, whereas Ctr9 bridges between them. Consistent with the structural observations, the initiation factor TFIIF impairs Paf1C binding to Pol II, whereas the elongation factor TFIIS enhances it. We further show that Paf1C is globally required for normal mRNA transcription in yeast. These results provide a three-dimensional framework for further analysis of Paf1C function in transcription through chromatin.

Size-Dependent Expression of the Mitotic Activator Cdc25 Suggests a Mechanism of Size Control in Fission Yeast

Proper cell size is essential for cellular function. Nonetheless, despite more than 100 years of work on the subject, the mechanisms that maintain cell-size homeostasis are largely mysterious [1]. Cells in growing populations maintain cell size within a narrow range by coordinating growth and division. Bacterial and eukaryotic cells both demonstrate homeostatic size control, which maintains population-level variation in cell size within a certain range and returns the population average to that range if it is perturbed [1, 2]. Recent work has proposed two different strategies for size control: budding yeast has been proposed to use an inhibitor-dilution strategy to regulate size at the G1/S transition [3], whereas bacteria appear to use an adder strategy, in which a fixed amount of growth each generation causes cell size to converge on a stable average [4-6]. Here we present evidence that cell size in the fission yeast Schizosaccharomyces pombe is regulated by a third strategy: the size-dependent expression of the mitotic activator Cdc25. cdc25 transcript levels are regulated such that smaller cells express less Cdc25 and larger cells express more Cdc25, creating an increasing concentration of Cdc25 as cells grow and providing a mechanism for cells to trigger cell division when they reach a threshold concentration of Cdc25. Because regulation of mitotic entry by Cdc25 is well conserved, this mechanism may provide a widespread solution to the problem of size control in eukaryotes.

Extremely fast and incredibly close: cotranscriptional splicing in budding yeast

RNA splicing, an essential part of eukaryotic pre-messenger RNA processing, can be simultaneous with transcription by RNA polymerase II. Here, we compare and review independent next-generation sequencing methods that jointly quantify transcription and splicing in budding yeast. For many yeast transcripts, splicing is fast, taking place within seconds of intron transcription, while polymerase is within a few dozens of nucleotides of the 3' splice site. Ribosomal protein transcripts are spliced particularly fast and cotranscriptionally. However, some transcripts are spliced inefficiently or mainly post-transcriptionally. Intron-mediated regulation of some genes is likely to be cotranscriptional. We suggest that intermediates of the splicing reaction, missing from current data sets, may hold key information about splicing kinetics.

The conserved protein Seb1 drives transcription termination by binding RNA polymerase II and nascent RNA

Termination of RNA polymerase II (Pol II) transcription is an important step in the transcription cycle, which involves the dislodgement of polymerase from DNA, leading to release of a functional transcript. Recent studies have identified the key players required for this process and showed that a common feature of these proteins is a conserved domain that interacts with the phosphorylated C-terminus of Pol II (CTD-interacting domain, CID). However, the mechanism by which transcription termination is achieved is not understood. Using genome-wide methods, here we show that the fission yeast CID-protein Seb1 is essential for termination of protein-coding and non-coding genes through interaction with S2-phosphorylated Pol II and nascent RNA. Furthermore, we present the crystal structures of the Seb1 CTD- and RNA-binding modules. Unexpectedly, the latter reveals an intertwined two-domain arrangement of a canonical RRM and second domain. These results provide important insights into the mechanism underlying eukaryotic transcription termination.

Spt5 Plays Vital Roles in the Control of Sense and Antisense Transcription Elongation

Spt5 is an essential and conserved factor that functions in transcription and co-transcriptional processes. However, many aspects of the requirement for Spt5 in transcription are poorly understood. We have analyzed the consequences of Spt5 depletion in Schizosaccharomyces pombe using four genome-wide approaches. Our results demonstrate that Spt5 is crucial for a normal rate of RNA synthesis and distribution of RNAPII over transcription units. In the absence of Spt5, RNAPII localization changes dramatically, with reduced levels and a relative accumulation over the first ∼500 bp, suggesting that Spt5 is required for transcription past a barrier. Spt5 depletion also results in widespread antisense transcription initiating within this barrier region. Deletions of this region alter the distribution of RNAPII on the sense strand, suggesting that the barrier observed after Spt5 depletion is normally a site at which Spt5 stimulates elongation. Our results reveal a global requirement for Spt5 in transcription elongation.

Perfect timing: splicing and transcription rates in living cells

An important step toward understanding gene regulation is the elucidation of the time necessary for the completion of individual steps. Measurement of reaction rates can reveal potential nodes for regulation. For example, measurements of in vivo transcription elongation rates reveal regulation by DNA sequence, gene architecture, and chromatin. Pre-mRNA splicing is regulated by transcription elongation rates and vice versa, yet the rates of RNA processing reactions remain largely elusive. Since the 1980s, numerous model systems and approaches have been used to determine the precise timing of splicing in vivo. Because splicing can be co-transcriptional, the position of Pol II when splicing is detected has been used as a proxy for time by some investigators. In addition to these 'distance-based' measurements, 'time-based' measurements have been possible through live cell imaging, metabolic labeling of RNA, and gene induction. Yet splicing rates can be convolved by the time it takes for transcription, spliceosome assembly and spliceosome disassembly. The variety of assays and systems used has, perhaps not surprisingly, led to reports of widely differing splicing rates in vivo. Recently, single molecule RNA-seq has indicated that splicing occurs more quickly than previously deduced. Here we comprehensively review these findings and discuss evidence that splicing and transcription rates are closely coordinated, facilitating the efficiency of gene expression. On the other hand, introduction of splicing delays through as yet unknown mechanisms provide opportunity for regulation. More work is needed to understand how cells optimize the rates of gene expression for a range of biological conditions. WIREs RNA 2017, 8:e1401. doi: 10.1002/wrna.1401 For further resources related to this article, please visit the WIREs website.

The dynamic landscape of fission yeast meiosis alternative-splice isoforms

Alternative splicing increases the diversity of transcriptomes and proteomes in metazoans. The extent to which alternative splicing is active and functional in unicellular organisms is less understood. Here, we exploit a single-molecule long-read sequencing technique and develop an open-source software program called SpliceHunter to characterize the transcriptome in the meiosis of fission yeast. We reveal 14,353 alternative splicing events in 17,669 novel isoforms at different stages of meiosis, including antisense and read-through transcripts. Intron retention is the major type of alternative splicing, followed by alternate “intron in exon.” Seven hundred seventy novel transcription units are detected; 53 of the predicted proteins show homology in other species and form theoretical stable structures. We report the complexity of alternative splicing along isoforms, including 683 intra-molecularly co-associated intron pairs. We compare the dynamics of novel isoforms based on the number of supporting full-length reads with those of annotated isoforms and explore the translational capacity and quality of novel isoforms. The evaluation of these factors indicates that the majority of novel isoforms are unlikely to be both condition-specific and translatable but consistent with the possibility of biologically functional novel isoforms. Moreover, the co-option of these unusual transcripts into newly born genes seems likely. Together, the results of this study highlight the diversity and dynamics at the isoform level in the sexual development of fission yeast.

Measures of RNA metabolism rates: Toward a definition at the level of single bonds

We give an overview of experimental and computational methods to estimate RNA metabolism rates genome-wide. We then advocate a local definition of RNA metabolism rate at the level of individual phosphodiester bonds. Rates of formation and disappearance of individual bonds are unambiguously defined, in contrast to rates of complete transcripts. We show that over previous approaches, the recently developed transient transcriptome sequencing (TT-seq) protocol allows for estimation of metabolism rates of individual bonds with least positional bias.

Analysis of the association between codon optimality and mRNA stability in Schizosaccharomyces pombe

Background

Recent experiments have shown that codon optimality is a major determinant of mRNA stability in Saccharomyces cerevisiae and that this phenomenon may be conserved in Escherichia coli and some metazoans, although work in Neurospora crassa is not consistent with this model.

Results

We examined the association between codon optimality and mRNA stability in the fission yeast Schizosaccharomyces pombe. Our analysis revealed the following points. First, we observe a genome-wide association between codon optimality and mRNA stability also in S. pombe, suggesting evolutionary conservation of the phenomenon. Second, in both S. pombe and S. cerevisiae, mRNA synthesis rates are also correlated at the genome-wide analysis with codon optimality, suggesting that the long-appreciated association between codon optimality and mRNA abundance is due to regulation of both mRNA synthesis and degradation. However, when we examined correlation of codon optimality and either mRNA half-lives or synthesis rates controlling for mRNA abundance, codon optimality was still positively correlated with mRNA half-lives in S. cerevisiae, but the association was no longer significant for mRNA half-lives in S. pombe or for synthesis rates in either organism. This illustrates how only the pairwise analysis of multiple correlating variables may limit these types of analyses. Finally, in S. pombe, codon optimality is associated with known DNA/RNA sequence motifs that are associated with mRNA production/stability, suggesting these two features have been under similar selective pressures for optimal gene expression.

Conclusions

Consistent with the emerging body of studies, this study suggests that the association between codon optimality and mRNA stability may be a broadly conserved phenomenon. It also suggests that the association can be explained at least in part by independent adaptations of codon optimality and other transcript features for elevated expression during evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3237-6) contains supplementary material, which is available to authorized users.

A conserved role of the RSC chromatin remodeler in the establishment of nucleosome-depleted regions

The occupancy of nucleosomes governs access to the eukaryotic genomes and results from a combination of biophysical features and the effect of ATP-dependent remodelling complexes. Most promoter regions show a conserved pattern characterized by a nucleosome-depleted region (NDR) flanked by nucleosomal arrays. The conserved RSC remodeler was reported to be critical to establish NDR in vivo in budding yeast but other evidences suggested that this activity may not be conserved in fission yeast. By reanalysing and expanding previously published data, we propose that NDR formation requires, at least partially, RSC in both yeast species. We also discuss the most prominent biological role of RSC and the possibility that non-essential subunits do not define alternate versions of the complex.

Statistical modeling of isoform splicing dynamics from RNA-seq time series data

MOTIVATION

Isoform quantification is an important goal of RNA-seq experiments, yet it remains problematic for genes with low expression or several isoforms. These difficulties may in principle be ameliorated by exploiting correlated experimental designs, such as time series or dosage response experiments. Time series RNA-seq experiments, in particular, are becoming increasingly popular, yet there are no methods that explicitly leverage the experimental design to improve isoform quantification.

RESULTS

Here, we present DICEseq, the first isoform quantification method tailored to correlated RNA-seq experiments. DICEseq explicitly models the correlations between different RNA-seq experiments to aid the quantification of isoforms across experiments. Numerical experiments on simulated datasets show that DICEseq yields more accurate results than state-of-the-art methods, an advantage that can become considerable at low coverage levels. On real datasets, our results show that DICEseq provides substantially more reproducible and robust quantifications, increasing the correlation of estimates from replicate datasets by up to 10% on genes with low or moderate expression levels (bottom third of all genes). Furthermore, DICEseq permits to quantify the trade-off between temporal sampling of RNA and depth of sequencing, frequently an important choice when planning experiments. Our results have strong implications for the design of RNA-seq experiments, and offer a novel tool for improved analysis of such datasets.

AVAILABILITY AND IMPLEMENTATION

Python code is freely available at http://diceseq.sf.net

CONTACT

G.Sanguinetti@ed.ac.uk

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.