Global donor and acceptor splicing site kinetics in human cells

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.

Gene-specific RNA homeostasis revealed by perturbation of coactivator complexes

Transcript buffering entails the reciprocal modulation of mRNA synthesis and degradation rates to maintain stable RNA levels under varying cellular conditions. Current research supports a global, non-sequence-specific connection between mRNA synthesis and degradation, but the underlying mechanisms are still unclear. In this study, we investigated changes in RNA metabolism following acute depletion of TIP60/KAT5, the acetyltransferase subunit of the NuA4 transcriptional coactivator complex, in mouse embryonic stem cells. By combining RNA sequencing of nuclear, cytoplasmic, and newly synthesised transcript fractions with biophysical modelling, we demonstrate that TIP60 predominantly enhances transcription of numerous genes, while a smaller set of genes undergoes TIP60-dependent transcriptional repression. Surprisingly, transcription changes caused by TIP60 depletion were offset by corresponding changes in RNA nuclear export and cytoplasmic stability, indicating gene-specific buffering mechanisms. Similarly, disruption of the unrelated ATAC coactivator complex also resulted in gene-specific transcript buffering. These findings reveal that transcript buffering functions at a gene-specific level and suggest that cells dynamically adjust RNA splicing, export, and degradation in response to individual RNA synthesis alterations, thereby sustaining cellular homeostasis.

Nuclear export is a limiting factor in eukaryotic mRNA metabolism

The eukaryotic mRNA life cycle includes transcription, nuclear mRNA export and degradation. To quantify all these processes simultaneously, we perform thiol-linked alkylation after metabolic labeling of RNA with 4-thiouridine (4sU), followed by sequencing of RNA (SLAM-seq) in the nuclear and cytosolic compartments of human cancer cells. We develop a model that reliably quantifies mRNA-specific synthesis, nuclear export, and nuclear and cytosolic degradation rates on a genome-wide scale. We find that nuclear degradation of polyadenylated mRNA is negligible and nuclear mRNA export is slow, while cytosolic mRNA degradation is comparatively fast. Consequently, an mRNA molecule generally spends most of its life in the nucleus. We also observe large differences in the nuclear export rates of different 3'UTR transcript isoforms. Furthermore, we identify genes whose expression is abruptly induced upon metabolic labeling. These transcripts are exported substantially faster than average mRNAs, suggesting the existence of alternative export pathways. Our results highlight nuclear mRNA export as a limiting factor in mRNA metabolism and gene regulation.

A kinetic dichotomy between mitochondrial and nuclear gene expression processes

Oxidative phosphorylation (OXPHOS) complexes, encoded by both mitochondrial and nuclear DNA, are essential producers of cellular ATP, but how nuclear and mitochondrial gene expression steps are coordinated to achieve balanced OXPHOS subunit biogenesis remains unresolved. Here, we present a parallel quantitative analysis of the human nuclear and mitochondrial messenger RNA (mt-mRNA) life cycles, including transcript production, processing, ribosome association, and degradation. The kinetic rates of nearly every stage of gene expression differed starkly across compartments. Compared with nuclear mRNAs, mt-mRNAs were produced 1,100-fold more, degraded 7-fold faster, and accumulated to 160-fold higher levels. Quantitative modeling and depletion of mitochondrial factors LRPPRC and FASTKD5 identified critical points of mitochondrial regulatory control, revealing that the mitonuclear expression disparities intrinsically arise from the highly polycistronic nature of human mitochondrial pre-mRNA. We propose that resolving these differences requires a 100-fold slower mitochondrial translation rate, illuminating the mitoribosome as a nexus of mitonuclear co-regulation.

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.

mRNA initiation and termination are spatially coordinated

The expression of a precise mRNA transcriptome is crucial for establishing cell identity and function, with dozens of alternative isoforms produced for a single gene sequence. The regulation of mRNA isoform usage occurs by the coordination of co-transcriptional mRNA processing mechanisms across a gene. Decisions involved in mRNA initiation and termination underlie the largest extent of mRNA isoform diversity, but little is known about any relationships between decisions at both ends of mRNA molecules. Here, we systematically profile the joint usage of mRNA transcription start sites (TSSs) and polyadenylation sites (PASs) across tissues and species. Using both short and long read RNA-seq data, we observe that mRNAs preferentially using upstream TSSs also tend to use upstream PASs, and congruently, the usage of downstream sites is similarly paired. This observation suggests that mRNA 5' end choice may directly influence mRNA 3' ends. Our results suggest a novel "Positional Initiation-Termination Axis" (PITA), in which the usage of alternative terminal sites are coupled based on the order in which they appear in the genome. PITA isoforms are more likely to encode alternative protein domains and use conserved sites. PITA is strongly associated with the length of genomic features, such that PITA is enriched in longer genes with more area devoted to regions that regulate alternative 5' or 3' ends. Strikingly, we found that PITA genes are more likely than non-PITA genes to have multiple, overlapping chromatin structural domains related to pairing of ordinally coupled start and end sites. In turn, PITA coupling is also associated with fast RNA Polymerase II (RNAPII) trafficking across these long gene regions. Our findings indicate that a combination of spatial and kinetic mechanisms couple transcription initiation and mRNA 3' end decisions based on ordinal position to define the expression mRNA isoforms.

A network of DZF proteins controls alternative splicing regulation and fidelity

Proteins containing DZF (domain associated with zinc fingers) modules play important roles throughout gene expression, from transcription to translation. Derived from nucleotidyltransferases but lacking catalytic residues, DZF domains serve as heterodimerization surfaces between DZF protein pairs. Three DZF proteins are widely expressed in mammalian tissues, ILF2, ILF3 and ZFR, which form mutually exclusive ILF2-ILF3 and ILF2-ZFR heterodimers. Using eCLIP-Seq, we find that ZFR binds across broad intronic regions to regulate the alternative splicing of cassette and mutually exclusive exons. ZFR preferentially binds dsRNA in vitro and is enriched on introns containing conserved dsRNA elements in cells. Many splicing events are similarly altered upon depletion of any of the three DZF proteins; however, we also identify independent and opposing roles for ZFR and ILF3 in alternative splicing regulation. Along with widespread involvement in cassette exon splicing, the DZF proteins control the fidelity and regulation of over a dozen highly validated mutually exclusive splicing events. Our findings indicate that the DZF proteins form a complex regulatory network that leverages dsRNA binding by ILF3 and ZFR to modulate splicing regulation and fidelity.

Insplico: effective computational tool for studying splicing order of adjacent introns genome-wide with short and long RNA-seq reads

Although splicing occurs largely co-transcriptionally, the order by which introns are removed does not necessarily follow the order in which they are transcribed. Whereas several genomic features are known to influence whether or not an intron is spliced before its downstream neighbor, multiple questions related to adjacent introns' splicing order (AISO) remain unanswered. Here, we present Insplico, the first standalone software for quantifying AISO that works with both short and long read sequencing technologies. We first demonstrate its applicability and effectiveness using simulated reads and by recapitulating previously reported AISO patterns, which unveiled overlooked biases associated with long read sequencing. We next show that AISO around individual exons is remarkably constant across cell and tissue types and even upon major spliceosomal disruption, and it is evolutionarily conserved between human and mouse brains. We also establish a set of universal features associated with AISO patterns across various animal and plant species. Finally, we used Insplico to investigate AISO in the context of tissue-specific exons, particularly focusing on SRRM4-dependent microexons. We found that the majority of such microexons have non-canonical AISO, in which the downstream intron is spliced first, and we suggest two potential modes of SRRM4 regulation of microexons related to their AISO and various splicing-related features. Insplico is available on gitlab.com/aghr/insplico.

How Estrogen, Testosterone, and Sex Differences Influence Serum Immunoglobulin Isotype Patterns in Mice and Humans

Females often exhibit superior immune responses compared to males toward vaccines and pathogens such as influenza viruses and SARS-CoV-2. To help explain these differences, we first studied serum immunoglobulin isotype patterns in C57BL/6 male and female mice. We focused on IgG2b, an isotype that lends to virus control and that has been previously shown to be elevated in murine females compared to males. Improvements in IgG2b serum levels, and/or IgG2b ratios with other non-IgM isotypes, were observed when: (i) wildtype (WT) female mice were compared to estrogen receptor knockout mice (IgG2b, IgG2b/IgG3, IgG2b/IgG1, and IgG2b/IgA were all higher in WT mice), (ii) unmanipulated female mice were compared to ovariectomized mice (IgG2b/IgA was higher in unmanipulated animals), (iii) female mice were supplemented with estrogen in the context of an inflammatory insult (IgG2b and IgG2b/IgG3 were improved by estrogen supplementation), and (iv) male mice were supplemented with testosterone, a hormone that can convert to estrogen in vivo (IgG2b, IgG2b/IgG3, IgG2b/IgG1, and IgG2b/IgA were all improved by supplementation). We next examined data from three sets of previously described male and female human blood samples. In each case, there were higher IgG2 levels, and/or ratios of IgG2 with non-IgM isotypes, in human females compared to males. The effects of sex and sex hormones in the mouse and human studies were subtle, but frequent, suggesting that sex hormones represent only a fraction of the factors that influence isotype patterns. Examination of the gene loci suggested that upregulation of murine IgG2b or human IgG2 could be mediated by estrogen receptor binding to estrogen response elements and cytosine-adenine (CA) repeats upstream of respective Cγ genes. Given that murine IgG2b and human IgG2 lend to virus control, the isotype biases in females may be sufficient to improve outcomes following vaccination or infection. Future attention to sex hormone levels, and consequent immunoglobulin isotype patterns, in clinical trials are encouraged to support the optimization of vaccine and drug products for male and female hosts.

Single nucleotide polymorphism (SNP) in the doublesex ( dsx) gene splice sites and relevance for its alternative splicing in the malaria vector Anopheles gambiae

Background: Malaria burden continues to be significant in tropical regions, and conventional vector control methods are faced with challenges such as insecticide resistance. To overcome these challenges, additional vector control interventions are vital and include modern genetic approaches as well as classical methods like the sterile insect technique (SIT). In the major human malaria vector Anopheles gambiae, a candidate gene favourable for sterility induction is the doublesex ( dsx) gene, involved in mosquitos' somatic sexually dimorphic traits determination. However, the pathways that trigger the signal of dsx gene exon skipping alternative splicing mechanism in anopheline mosquitoes are not well characterized. This study aims to screen the An. gambiae dsx gene splice site sequences for single-nucleotide polymorphisms (SNPs) that could be critical to its alternative splicing. Methods: Variant annotation data from Ag1000G project phase 2 was analysed, in order to identify splice-relevant SNPs within acceptor and donor splice sites of the An. gambiae dsx gene ( Agdsx). Results: SNPs were found in both donor and acceptor sites of the Agdsx. No splice-relevant SNPs were identified in the female-specific intron 4 acceptor site and the corresponding region in males. Two SNPs (rs48712947, rs48712962) were found in the female-specific donor site of exon 5. They were not specific to either males or females as the rs48712947 was found in female mosquitoes from Cameroon, and in both males and females from Burkina Faso. In the other splice sites, the intron 3 acceptor site carried the greatest abundance of SNPs.   Conclusions: There were no gender association between the identified SNPs and the random distribution of these SNPs in mosquito populations. The SNPs in Agdsx splice sites are not critical for the alternative splicing. Other molecular mechanisms should be considered and investigated.

Inefficient splicing of long non-coding RNAs is associated with higher transcript complexity in human and mouse

Recent reports show that long non-coding RNAs (lncRNAs) have inefficient splicing and fewer alternative splice variants than mRNAs. Here, we have explored the efficiency of lncRNAs and mRNAs in producing various splice variants, given the number of exons in humans and mice. Intriguingly, lncRNAs produce more splice variants per exon, referred to as Transcript Complexity, than mRNAs. Most lncRNA splice variants are the product of the alternative last exon and exon skipping. LncRNAs and mRNAs with higher transcript complexity have shorter intron lengths. Longer exon length and GC/AG at 5'/3' splice sites are associated with higher transcript complexity in lncRNAs. Lastly, our results indicate that inefficient splicing of lncRNAs may facilitate multiple introns splicing and, thus, more spliced products per exon.

Single nucleotide polymorphism (SNP) in the doublesex (dsx) gene splice sites and relevance for its alternative splicing in the malaria vector Anopheles gambiae

Background: Malaria burden continues to be significant in tropical regions, and conventional vector control methods are faced with challenges such as insecticide resistance. To overcome these challenges, additional vector control interventions are vital and include modern genetic approaches as well as classical methods like the sterile insect technique (SIT). In the major human malaria vector Anopheles gambiae, a candidate gene favourable for sterility induction is the doublesex (dsx) gene, involved in mosquitos’ somatic sexually dimorphic traits determination. However, the pathways that trigger the signal of dsx gene exon skipping alternative splicing mechanism in anopheline mosquitoes are not well characterized. This study aims to screen the An. gambiae dsx gene splice site sequences for single-nucleotide polymorphisms (SNPs) that could be critical to its alternative splicing. Methods: Variant annotation data from Ag1000G project phase 2 was analysed, in order to identify splice-relevant SNPs within acceptor and donor splice sites of the An. gambiae dsx gene (Agdsx). Results: SNPs were found in both donor and acceptor sites of the Agdsx. No splice-relevant SNPs were identified in the female-specific intron 4 acceptor site and the corresponding region in males. Two SNPs (rs48712947, rs48712962) were found in the female-specific donor site of exon 5. They were not specific to either males or females as the rs48712947 was found in female mosquitoes from Cameroon, and in both males and females from Burkina Faso. In the other splice sites, the intron 3 acceptor site carried the greatest abundance of SNPs.   Conclusions: There were no gender association between the identified SNPs and the random distribution of these SNPs in mosquito populations. The SNPs in Agdsx splice sites are not critical for the alternative splicing. Other molecular mechanisms should be considered and investigated.

Profiling lariat intermediates reveals genetic determinants of early and late co-transcriptional splicing

Long introns with short exons in vertebrate genes are thought to require spliceosome assembly across exons (exon definition), rather than introns, thereby requiring transcription of an exon to splice an upstream intron. Here, we developed CoLa-seq (co-transcriptional lariat sequencing) to investigate the timing and determinants of co-transcriptional splicing genome wide. Unexpectedly, 90% of all introns, including long introns, can splice before transcription of a downstream exon, indicating that exon definition is not obligatory for most human introns. Still, splicing timing varies dramatically across introns, and various genetic elements determine this variation. Strong U2AF2 binding to the polypyrimidine tract predicts early splicing, explaining exon definition-independent splicing. Together, our findings question the essentiality of exon definition and reveal features beyond intron and exon length that are determinative for splicing timing.

The genetic and biochemical determinants of mRNA degradation rates in mammals

BACKGROUND

Degradation rate is a fundamental aspect of mRNA metabolism, and the factors governing it remain poorly characterized. Understanding the genetic and biochemical determinants of mRNA half-life would enable more precise identification of variants that perturb gene expression through post-transcriptional gene regulatory mechanisms.

RESULTS

We establish a compendium of 39 human and 27 mouse transcriptome-wide mRNA decay rate datasets. A meta-analysis of these data identified a prevalence of technical noise and measurement bias, induced partially by the underlying experimental strategy. Correcting for these biases allowed us to derive more precise, consensus measurements of half-life which exhibit enhanced consistency between species. We trained substantially improved statistical models based upon genetic and biochemical features to better predict half-life and characterize the factors molding it. Our state-of-the-art model, Saluki, is a hybrid convolutional and recurrent deep neural network which relies only upon an mRNA sequence annotated with coding frame and splice sites to predict half-life (r=0.77). The key novel principle learned by Saluki is that the spatial positioning of splice sites, codons, and RNA-binding motifs within an mRNA is strongly associated with mRNA half-life. Saluki predicts the impact of RNA sequences and genetic mutations therein on mRNA stability, in agreement with functional measurements derived from massively parallel reporter assays.

CONCLUSIONS

Our work produces a more robust ground truth for transcriptome-wide mRNA half-lives in mammalian cells. Using these revised measurements, we trained Saluki, a model that is over 50% more accurate in predicting half-life from sequence than existing models. Saluki succinctly captures many of the known determinants of mRNA half-life and can be rapidly deployed to predict the functional consequences of arbitrary mutations in the transcriptome.

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.

Oct4 differentially regulates chromatin opening and enhancer transcription in pluripotent stem cells

The transcription factor Oct4 is essential for the maintenance and induction of stem cell pluripotency, but its functional roles are not fully understood. Here, we investigate the functions of Oct4 by depleting and subsequently recovering it in mouse embryonic stem cells (ESCs) and conducting a time-resolved multiomics analysis. Oct4 depletion leads to an immediate loss of its binding to enhancers, accompanied by a decrease in mRNA synthesis from its target genes that are part of the transcriptional network that maintains pluripotency. Gradual decrease of Oct4 binding to enhancers does not immediately change the chromatin accessibility but reduces transcription of enhancers. Conversely, partial recovery of Oct4 expression results in a rapid increase in chromatin accessibility, whereas enhancer transcription does not fully recover. These results indicate different concentration-dependent activities of Oct4. Whereas normal ESC levels of Oct4 are required for transcription of pluripotency enhancers, low levels of Oct4 are sufficient to retain chromatin accessibility, likely together with other factors such as Sox2.

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.

Co-transcriptional splicing efficiency is a gene-specific feature that can be regulated by TGFβ

Differential splicing efficiency of specific introns is a mechanism that dramatically increases protein diversity, based on selection of alternative exons for the final mature mRNA. However, it is unclear whether splicing efficiency of introns within the same gene is coordinated and eventually regulated as a mechanism to control mature mRNA levels. Based on nascent chromatin-associated RNA-sequencing data, we now find that co-transcriptional splicing (CTS) efficiency tends to be similar between the different introns of a gene. We establish that two well-differentiated strategies for CTS efficiency exist, at the extremes of a gradient: short genes that produce high levels of pre-mRNA undergo inefficient splicing, while long genes with relatively low levels of pre-mRNA have an efficient splicing. Notably, we observe that genes with efficient CTS display a higher level of mature mRNA relative to their pre-mRNA levels. Further, we show that the TGFβ signal transduction pathway regulates the general CTS efficiency, causing changes in mature mRNA levels. Taken together, our data indicate that CTS efficiency is a gene-specific characteristic that can be regulated to control gene expression.

Co-transcriptional splicing efficiency is a gene-specific characteristic that can be regulated by TGFβ to modulate gene expression.

Single nucleotide polymorphism (SNP) in the doublesex (dsx) gene splice sites and relevance for its alternative splicing in the malaria vector Anopheles gambiae

Background: The malaria burden continues to be significant in tropical regions, and conventional vector control methods are faced with challenges such as insecticide resistance. To overcome these challenges, additional vector control interventions are vital and include modern genetic approaches as well as classical methods like the sterile insect technique (SIT). In the major human malaria vector Anopheles gambiae, a candidate gene favourable for sterility induction is the doublesex (dsx) gene, encoding somatic sexually dimorphic traits in mosquitoes. However, the mechanism that regulates the expression of this gene in anopheline mosquitoes is poorly understood. This study aimed to screen the An. gambiae dsx gene splice site sequences for single nucleotide polymorphisms (SNPs) that could be critical to its alternative splicing. Methods: Variant annotation data from Ag1000G project phase 2 was analysed, in order to identify splice-relevant SNPs within acceptor and donor splice sites of the An. gambiae dsx gene (Agdsx). Results: SNPs were found in both donor and acceptor sites of the Agdsx. No splice-relevant SNPs were identified in the female-specific intron 4 acceptor site and the corresponding region in males. Two SNPs (rs48712947, rs48712962) were found in the female-specific donor site of exon 5. They were not specific to either males or females as the rs48712947 was found in female mosquitoes from Cameroon, and in both males and females from Burkina Faso. In the other splice sites, the intron 3 acceptor site carried the greatest abundance of SNPs.   Conclusions: There were no gender association between the identified SNPs and the random distribution of these SNPs in mosquito populations. The SNPs in Agdsx splice sites are not critical for the alternative splicing. Other molecular mechanisms should be considered and investigated.

Splice site proximity influences alternative exon definition

Alternative splicing enables higher eukaryotes to expand mRNA diversity from a finite number of genes through highly combinatorial splice site selection mechanisms that are influenced by the sequence of competing splice sites, cis-regulatory elements binding trans-acting factors, the length of exons and introns harbouring alternative splice sites and RNA secondary structures at putative splice junctions. To test the hypothesis that the intron definition or exon definition modes of splice site recognition direct the selection of alternative splice patterns, we created a database of alternative splice site usage (ALTssDB). When alternative splice sites are embedded within short introns (intron definition), the 5' and 3' splice sites closest to each other across the intron preferentially pair, consistent with previous observations. However, when alternative splice sites are embedded within large flanking introns (exon definition), the 5' and 3' splice sites closest to each other across the exon are preferentially selected. Thus, alternative splicing decisions are influenced by the intron and exon definition modes of splice site recognition. The results demonstrate that the spliceosome pairs splice sites that are closest in proximity within the unit of initial splice site selection.

Evolution of the Early Spliceosomal Complex-From Constitutive to Regulated Splicing

Pre-mRNA splicing is a major process in the regulated expression of genes in eukaryotes, and alternative splicing is used to generate different proteins from the same coding gene. Splicing is a catalytic process that removes introns and ligates exons to create the RNA sequence that codifies the final protein. While this is achieved in an autocatalytic process in ancestral group II introns in prokaryotes, the spliceosome has evolved during eukaryogenesis to assist in this process and to finally provide the opportunity for intron-specific splicing. In the early stage of splicing, the RNA 5' and 3' splice sites must be brought within proximity to correctly assemble the active spliceosome and perform the excision and ligation reactions. The assembly of this first complex, termed E-complex, is currently the least understood process. We focused in this review on the formation of the E-complex and compared its composition and function in three different organisms. We highlight the common ancestral mechanisms in S. cerevisiae, S. pombe, and mammals and conclude with a unifying model for intron definition in constitutive and regulated co-transcriptional splicing.

What's Wrong in a Jump? Prediction and Validation of Splice Site Variants

Alternative splicing (AS) is a crucial process to enhance gene expression driving organism development. Interestingly, more than 95% of human genes undergo AS, producing multiple protein isoforms from the same transcript. Any alteration (e.g., nucleotide substitutions, insertions, and deletions) involving consensus splicing regulatory sequences in a specific gene may result in the production of aberrant and not properly working proteins. In this review, we introduce the key steps of splicing mechanism and describe all different types of genomic variants affecting this process (splicing variants in acceptor/donor sites or branch point or polypyrimidine tract, exonic, and deep intronic changes). Then, we provide an updated approach to improve splice variants detection. First, we review the main computational tools, including the recent Machine Learning-based algorithms, for the prediction of splice site variants, in order to characterize how a genomic variant interferes with splicing process. Next, we report the experimental methods to validate the predictive analyses are defined, distinguishing between methods testing RNA (transcriptomics analysis) or proteins (proteomics experiments). For both prediction and validation steps, benefits and weaknesses of each tool/procedure are accurately reported, as well as suggestions on which approaches are more suitable in diagnostic rather than in clinical research.

Aging, inflammation and DNA damage in the somatic testicular niche with idiopathic germ cell aplasia

Molecular mechanisms associated with human germ cell aplasia in infertile men remain undefined. Here we perform single-cell transcriptome profiling to highlight differentially expressed genes and pathways in each somatic cell type in testes of men with idiopathic germ cell aplasia. We identify immaturity of Leydig cells, chronic tissue inflammation, fibrosis, and senescence phenotype of the somatic cells, as well markers of chronic inflammation in the blood. We find that deregulated expression of parentally imprinted genes in myoid and immature Leydig cells, with relevant changes in the ratio of Lamin A/C transcripts and an active DNA damage response in Leydig and peritubular myoid cells are also indicative of senescence of the testicular niche. This study offers molecular insights into the pathogenesis of idiopathic germ cell aplasia.

Molecular mechanisms associated with human germ cell aplasia in infertile men remain undefined. Here the authors perform single-cell transcriptome profiling to highlight differentially expressed genes and pathways in each somatic cell type in testes of men with idiopathic germ cell aplasia.

The upstream 5' splice site remains associated to the transcription machinery during intron synthesis

In the earliest step of spliceosome assembly, the two splice sites flanking an intron are brought into proximity by U1 snRNP and U2AF along with other proteins. The mechanism that facilitates this intron looping is poorly understood. Using a CRISPR interference-based approach to halt RNA polymerase II transcription in the middle of introns in human cells, we discovered that the nascent 5′ splice site base pairs with a U1 snRNA that is tethered to RNA polymerase II during intron synthesis. This association functionally corresponds with splicing outcome, involves bona fide 5′ splice sites and cryptic intronic sites, and occurs transcriptome-wide. Overall, our findings reveal that the upstream 5′ splice sites remain attached to the transcriptional machinery during intron synthesis and are thus brought into proximity of the 3′ splice sites; potentially mediating the rapid splicing of long introns.

We know that most splicing reactions take place co-transcriptionally, but how the transcription machinery facilitate splicing of introns is unknown. Here the authors show that the 5′ splice site remains associated with the transcription machinery during intron synthesis through U1 snRNP, providing a basis for the rapid splicing reaction of introns.

Dynamics of transcriptional and post-transcriptional regulation

Despite gene expression programs being notoriously complex, RNA abundance is usually assumed as a proxy for transcriptional activity. Recently developed approaches, able to disentangle transcriptional and post-transcriptional regulatory processes, have revealed a more complex scenario. It is now possible to work out how synthesis, processing and degradation kinetic rates collectively determine the abundance of each gene's RNA. It has become clear that the same transcriptional output can correspond to different combinations of the kinetic rates. This underscores the fact that markedly different modes of gene expression regulation exist, each with profound effects on a gene's ability to modulate its own expression. This review describes the development of the experimental and computational approaches, including RNA metabolic labeling and mathematical modeling, that have been disclosing the mechanisms underlying complex transcriptional programs. Current limitations and future perspectives in the field are also discussed.

Single-Molecule Techniques to Study Chromatin

Besides the basic organization in nucleosome core particles (NCPs), eukaryotic chromatin is further packed through interactions with numerous protein complexes including transcription factors, chromatin remodeling and modifying enzymes. This nucleoprotein complex provides the template for many important biological processes, such as DNA replication, transcription, and DNA repair. Thus, to understand the molecular basis of these DNA transactions, it is critical to define individual changes of the chromatin structure at precise genomic regions where these machineries assemble and drive biological reactions. Single-molecule approaches provide the only possible solution to overcome the heterogenous nature of chromatin and monitor the behavior of individual chromatin transactions in real-time. In this review, we will give an overview of currently available single-molecule methods to obtain mechanistic insights into nucleosome positioning, histone modifications and DNA replication and transcription analysis-previously unattainable with population-based assays.

Deletion of a pseudogene within a fragile site triggers the oncogenic expression of the mitotic CCSER1 gene

The oncogenic role of common fragile sites (CFS), focal and pervasive gaps in the cancer genome arising from replicative stress, remains controversial. Exploiting the TCGA dataset, we found that in most CFS the genes residing within the associated focal deletions are down-regulated, including proteins involved in tumour immune recognition. In a subset of CFS, however, the residing genes are surprisingly overexpressed. Within the most frequent CFS in this group, FRA4F, which is deleted in up to 18% of cancer cases and harbours the CCSER1 gene, we identified a region which includes an intronic, antisense pseudogene, TMSB4XP8. TMSB4XP8 focal ablation or transcriptional silencing elicits the overexpression of CCSER1, through a cis-acting mechanism. CCSER1 overexpression increases proliferation and triggers centrosome amplifications, multinuclearity, and aberrant mitoses. Accordingly, FRA4F is associated in patient samples to mitotic genes deregulation and genomic instability. As a result, cells overexpressing CCSER1 become sensitive to the treatment with aurora kinase inhibitors. Our findings point to a novel tumourigenic mechanism where focal deletions increase the expression of a new class of "dormant" oncogenes.

Efficient RNA polymerase II pause release requires U2 snRNP function

Transcription by RNA polymerase II (Pol II) is coupled to pre-mRNA splicing, but the underlying mechanisms remain poorly understood. Co-transcriptional splicing requires assembly of a functional spliceosome on nascent pre-mRNA, but whether and how this influences Pol II transcription remains unclear. Here we show that inhibition of pre-mRNA branch site recognition by the spliceosome component U2 snRNP leads to a widespread and strong decrease in new RNA synthesis from human genes. Multiomics analysis reveals that inhibition of U2 snRNP function increases the duration of Pol II pausing in the promoter-proximal region, impairs recruitment of the pause release factor P-TEFb, and reduces Pol II elongation velocity at the beginning of genes. Our results indicate that efficient release of paused Pol II into active transcription elongation requires the formation of functional spliceosomes and that eukaryotic mRNA biogenesis relies on positive feedback from the splicing machinery to the transcription machinery.

Revealing nascent RNA processing dynamics with nano-COP

During maturation, eukaryotic precursor RNAs undergo processing events including intron splicing, 3'-end cleavage, and polyadenylation. Here we describe nanopore analysis of co-transcriptional processing (nano-COP), a method for probing the timing and patterns of RNA processing. An extension of native elongating transcript sequencing, which quantifies transcription genome-wide through short-read sequencing of nascent RNA 3' ends, nano-COP uses long-read nascent RNA sequencing to observe global patterns of RNA processing. First, nascent RNA is stringently purified through a combination of 4-thiouridine metabolic labeling and cellular fractionation. In contrast to cDNA or short-read-based approaches relying on reverse transcription or amplification, the sample is sequenced directly through nanopores to reveal the native context of nascent RNA. nano-COP identifies both active transcription sites and splice isoforms of single RNA molecules during synthesis, providing insight into patterns of intron removal and the physical coupling between transcription and splicing. The nano-COP protocol yields data within 3 d.

Characterizing RNA stability genome-wide through combined analysis of PRO-seq and RNA-seq data

BACKGROUND

The concentrations of distinct types of RNA in cells result from a dynamic equilibrium between RNA synthesis and decay. Despite the critical importance of RNA decay rates, current approaches for measuring them are generally labor-intensive, limited in sensitivity, and/or disruptive to normal cellular processes. Here, we introduce a simple method for estimating relative RNA half-lives that is based on two standard and widely available high-throughput assays: Precision Run-On sequencing (PRO-seq) and RNA sequencing (RNA-seq).

RESULTS

Our method treats PRO-seq as a measure of transcription rate and RNA-seq as a measure of RNA concentration, and estimates the rate of RNA decay required for a steady-state equilibrium. We show that this approach can be used to assay relative RNA half-lives genome-wide, with good accuracy and sensitivity for both coding and noncoding transcription units. Using a structural equation model (SEM), we test several features of transcription units, nearby DNA sequences, and nearby epigenomic marks for associations with RNA stability after controlling for their effects on transcription. We find that RNA splicing-related features are positively correlated with RNA stability, whereas features related to miRNA binding and DNA methylation are negatively correlated with RNA stability. Furthermore, we find that a measure based on U1 binding and polyadenylation sites distinguishes between unstable noncoding and stable coding transcripts but is not predictive of relative stability within the mRNA or lincRNA classes. We also identify several histone modifications that are associated with RNA stability.

CONCLUSION

We introduce an approach for estimating the relative half-lives of individual RNAs. Together, our estimation method and systematic analysis shed light on the pervasive impacts of RNA stability on cellular RNA concentrations.

Transcription activation depends on the length of the RNA polymerase II C-terminal domain

Eukaryotic RNA polymerase II (Pol II) contains a tail‐like, intrinsically disordered carboxy‐terminal domain (CTD) comprised of heptad‐repeats, that functions in coordination of the transcription cycle and in coupling transcription to co‐transcriptional processes. The CTD repeat number varies between species and generally increases with genome size, but the reasons for this are unclear. Here, we show that shortening the CTD in human cells to half of its length does not generally change pre‐mRNA synthesis or processing in cells. However, CTD shortening decreases the duration of promoter‐proximal Pol II pausing, alters transcription of putative enhancer elements, and delays transcription activation after stimulation of the MAP kinase pathway. We suggest that a long CTD is required for efficient enhancer‐dependent recruitment of Pol II to target genes for their rapid activation.

Transcription and splicing: A two-way street

RNA synthesis by RNA polymerase II and RNA processing are closely coupled during the transcription cycle of protein-coding genes. This coupling affords opportunities for quality control and regulation of gene expression and the effects can go in both directions. For example, polymerase speed can affect splice site selection and splicing can increase transcription and affect the chromatin landscape. Here we review the many ways that transcription and splicing influence one another, including how splicing "talks back" to transcription. We will also place the connections between transcription and splicing in the context of other RNA processing events that define the exons that will make up the final mRNA. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Processing > Splicing Regulation/Alternative Splicing.

Intrinsic Regulatory Role of RNA Structural Arrangement in Alternative Splicing Control

Alternative splicing is a highly sophisticated process, playing a significant role in posttranscriptional gene expression and underlying the diversity and complexity of organisms. Its regulation is multilayered, including an intrinsic role of RNA structural arrangement which undergoes time- and tissue-specific alterations. In this review, we describe the principles of RNA structural arrangement and briefly decipher its cis- and trans-acting cellular modulators which serve as crucial determinants of biological functionality of the RNA structure. Subsequently, we engage in a discussion about the RNA structure-mediated mechanisms of alternative splicing regulation. On one hand, the impairment of formation of optimal RNA structures may have critical consequences for the splicing outcome and further contribute to understanding the pathomechanism of severe disorders. On the other hand, the structural aspects of RNA became significant features taken into consideration in the endeavor of finding potential therapeutic treatments. Both aspects have been addressed by us emphasizing the importance of ongoing studies in both fields.

Nascent Transcript Folding Plays a Major Role in Determining RNA Polymerase Elongation Rates

Transcription elongation rates influence RNA processing, but sequence-specific regulation is poorly understood. We addressed this in vivo, analyzing RNAPI in S. cerevisiae. Mapping RNAPI by Miller chromatin spreads or UV crosslinking revealed 5' enrichment and strikingly uneven local polymerase occupancy along the rDNA, indicating substantial variation in transcription speed. Two features of the nascent transcript correlated with RNAPI distribution: folding energy and GC content in the transcription bubble. In vitro experiments confirmed that strong RNA structures close to the polymerase promote forward translocation and limit backtracking, whereas high GC in the transcription bubble slows elongation. A mathematical model for RNAPI elongation confirmed the importance of nascent RNA folding in transcription. RNAPI from S. pombe was similarly sensitive to transcript folding, as were S. cerevisiae RNAPII and RNAPIII. For RNAPII, unstructured RNA, which favors slowed elongation, was associated with faster cotranscriptional splicing and proximal splice site use, indicating regulatory significance for transcript folding.

Direct RNA Sequencing for the Study of Synthesis, Processing, and Degradation of Modified Transcripts

It has been known for a few decades that transcripts can be marked by dozens of different modifications. Yet, we are just at the beginning of charting these marks and understanding their functional impact. High-quality methods were developed for the profiling of some of these marks, and approaches to finely study their impact on specific phases of the RNA life-cycle are available, including RNA metabolic labeling. Thanks to these improvements, the most abundant marks, including N6-methyladenosine, are emerging as important determinants of the fate of marked RNAs. However, we still lack approaches to directly study how the set of marks for a given RNA molecule shape its fate. In this perspective, we first review current leading approaches in the field. Then, we propose an experimental and computational setup, based on direct RNA sequencing and mathematical modeling, to decipher the functional consequences of RNA modifications on the fate of individual RNA molecules and isoforms.

Splicing Kinetics and Coordination Revealed by Direct Nascent RNA Sequencing through Nanopores

Understanding how splicing events are coordinated across numerous introns in metazoan RNA transcripts requires quantitative analyses of transient RNA processing events in living cells. We developed nanopore analysis of co-transcriptional processing (nano-COP), in which nascent RNAs are directly sequenced through nanopores, exposing the dynamics and patterns of RNA splicing without biases introduced by amplification. Long nano-COP reads reveal that, in human and Drosophila cells, splicing occurs after RNA polymerase II transcribes several kilobases of pre-mRNA, suggesting that metazoan splicing transpires distally from the transcription machinery. Inhibition of the branch-site recognition complex SF3B rapidly diminished global co-transcriptional splicing. We found that splicing order does not strictly follow the order of transcription and is associated with cis-acting elements, alternative splicing, and RNA-binding factors. Further, neighboring introns in human cells tend to be spliced concurrently, implying that splicing of these introns occurs cooperatively. Thus, nano-COP unveils the organizational complexity of RNA processing.

Using TTchem-seq for profiling nascent transcription and measuring transcript elongation

The dynamics of transcription can be studied genome wide by high-throughput sequencing of nascent and newly synthesized RNA. 4-thiouridine (4SU) labeling in vivo enables the specific capture of such new transcripts, with 4SU residues being tagged by biotin linkers and captured using streptavidin beads before library production and high-throughput sequencing. To achieve high-resolution profiles of transcribed regions, an RNA fragmentation step before biotin tagging was introduced, in an approach known as transient transcriptome sequencing (TT-seq). We recently introduced a chemical approach for RNA fragmentation that we refer to as TTchem-seq. We describe how TTchem-seq can be used in combination with transient inhibition of early elongation using the reversible CDK9 inhibitor, 5,6-dichlorobenzimidazole 1-β-D-ribofuranoside (DRB), to measure RNA polymerase II (RNAPII) elongation rates in vivo, a technique we call DRB/TTchem-seq. Here, we provide detailed protocols for carrying out TTchem-seq and DRB/TTchem-seq, including computational analysis. Experiments and data analysis can be performed over a period of 10-13 d and require molecular biology and bioinformatics skills.

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.