A systems view of spliceosomal assembly and branchpoints with iCLIP

RBM22 regulates RNA polymerase II 5' pausing, elongation rate, and termination by coordinating 7SK-P-TEFb complex and SPT5

BACKGROUND

Splicing factors are vital for the regulation of RNA splicing, but some have also been implicated in regulating transcription. The underlying molecular mechanisms of their involvement in transcriptional processes remain poorly understood.

RESULTS

Here, we describe a direct role of splicing factor RBM22 in coordinating multiple steps of RNA Polymerase II (RNAPII) transcription in human cells. The RBM22 protein widely occupies the RNAPII-transcribed gene locus in the nucleus. Loss of RBM22 promotes RNAPII pause release, reduces elongation velocity, and provokes transcriptional readthrough genome-wide, coupled with production of transcripts containing sequences from downstream of the gene. RBM22 preferentially binds to the hyperphosphorylated, transcriptionally engaged RNAPII and coordinates its dynamics by regulating the homeostasis of the 7SK-P-TEFb complex and the association between RNAPII and SPT5 at the chromatin level.

CONCLUSIONS

Our results uncover the multifaceted role of RBM22 in orchestrating the transcriptional program of RNAPII and provide evidence implicating a splicing factor in both RNAPII elongation kinetics and termination control.

Mapping protein-RNA binding in plants with individual-nucleotide-resolution UV cross-linking and immunoprecipitation (plant iCLIP2)

Despite crucial roles of RNA-binding proteins (RBPs) in plant physiology and development, methods for determining their transcriptome-wide binding landscape are less developed than those used in other model organisms. Cross-linking and immunoprecipitation (CLIP) methods (based on UV-mediated generation of covalent bonds between RNAs and cognate RBPs in vivo, purification of the cross-linked complexes and identification of the co-purified RNAs by high-throughput sequencing) have been applied mainly in mammalian cells growing in monolayers or in translucent tissue. We have developed plant iCLIP2, an efficient protocol for performing individual-nucleotide-resolution CLIP (iCLIP) in plants, tailored to overcome the experimental hurdles posed by plant tissue. We optimized the UV dosage to efficiently cross-link RNA and proteins in plants and expressed epitope-tagged RBPs under the control of their native promoters in loss-of-function mutants. We select epitopes for which nanobodies are available, allowing stringent conditions for immunopurification of the RNA-protein complexes to be established. To overcome the inherently high RNase content of plant cells, RNase inhibitors are added and the limited RNA fragmentation step is modified. We combine the optimized isolation of RBP-bound RNAs with iCLIP2, a streamlined protocol that greatly enhances the efficiency of library preparation for high-throughput sequencing. Plant researchers with experience in molecular biology and handling of RNA can complete this iCLIP2 protocol in ~5 d. Finally, we describe a bioinformatics workflow to determine targets of Arabidopsis RBPs from iCLIP data, covering all steps from downloading sequencing reads to identifying cross-linking events ( https://github.com/malewins/Plant-iCLIPseq ), and present the R/Bioconductor package BindingSiteFinder to extract reproducible binding sites ( https://bioconductor.org/packages/release/bioc/html/BindingSiteFinder.html ).

Variations of intronic branchpoint motif: identification and functional implications in splicing and disease

The branchpoint (BP) motif is an essential intronic element for spliceosomal pre-mRNA splicing. In mammals, its sequence composition, distance to the downstream exon, and number of BPs per 3´ splice site are highly variable, unlike the GT/AG dinucleotides at the intron ends. These variations appear to provide evolutionary advantages for fostering alternative splicing, satisfying more diverse cellular contexts, and promoting resilience to genetic changes, thus contributing to an extra layer of complexity for gene regulation. Importantly, variants in the BP motif itself or in genes encoding BP-interacting factors cause human genetic diseases or cancers, highlighting the critical function of BP motif and the need to precisely identify functional BPs for faithful interpretation of their roles in splicing. In this perspective, we will succinctly summarize the major findings related to BP motif variations, discuss the relevant issues/challenges, and provide our insights.

FUBP1 is a general splicing factor facilitating 3' splice site recognition and splicing of long introns

Splicing of pre-mRNAs critically contributes to gene regulation and proteome expansion in eukaryotes, but our understanding of the recognition and pairing of splice sites during spliceosome assembly lacks detail. Here, we identify the multidomain RNA-binding protein FUBP1 as a key splicing factor that binds to a hitherto unknown cis-regulatory motif. By collecting NMR, structural, and in vivo interaction data, we demonstrate that FUBP1 stabilizes U2AF2 and SF1, key components at the 3' splice site, through multivalent binding interfaces located within its disordered regions. Transcriptional profiling and kinetic modeling reveal that FUBP1 is required for efficient splicing of long introns, which is impaired in cancer patients harboring FUBP1 mutations. Notably, FUBP1 interacts with numerous U1 snRNP-associated proteins, suggesting a unique role for FUBP1 in splice site bridging for long introns. We propose a compelling model for 3' splice site recognition of long introns, which represent 80% of all human introns.

Toward a systems view on RNA-binding proteins and associated RNAs in plants: Guilt by association

RNA-binding proteins (RBPs) have a broad impact on most biochemical, physiological, and developmental processes in a plant's life. RBPs engage in an on-off relationship with their RNA partners, accompanying virtually every stage in RNA processing and function. While the function of a plethora of RBPs in plant development and stress responses has been described, we are lacking a systems-level understanding of components in RNA-based regulation. Novel techniques have substantially enlarged the compendium of proteins with experimental evidence for binding to RNAs in the cell, the RNA-binding proteome. Furthermore, ribonomics methods have been adapted for use in plants to profile the in vivo binding repertoire of RBPs genome-wide. Here, we discuss how recent technological achievements have provided novel insights into the mode of action of plant RBPs at a genome-wide scale. Furthermore, we touch upon two emerging topics, the connection of RBPs to phase separation in the cell and to extracellular RNAs. Finally, we define open questions to be addressed to move toward an integrated understanding of RBP function.

Temporal-iCLIP captures co-transcriptional RNA-protein interactions

Dynamic RNA-protein interactions govern the co-transcriptional packaging of RNA polymerase II (RNAPII)-derived transcripts. Yet, our current understanding of this process in vivo primarily stems from steady state analysis. To remedy this, we here conduct temporal-iCLIP (tiCLIP), combining RNAPII transcriptional synchronisation with UV cross-linking of RNA-protein complexes at serial timepoints. We apply tiCLIP to the RNA export adaptor, ALYREF; a component of the Nuclear Exosome Targeting (NEXT) complex, RBM7; and the nuclear cap binding complex (CBC). Regardless of function, all tested factors interact with nascent RNA as it exits RNAPII. Moreover, we demonstrate that the two transesterification steps of pre-mRNA splicing temporally separate ALYREF and RBM7 binding to splicing intermediates, and that exon-exon junction density drives RNA 5'end binding of ALYREF. Finally, we identify underappreciated steps in snoRNA 3'end processing performed by RBM7. Altogether, our data provide a temporal view of RNA-protein interactions during the early phases of transcription.

C-terminal determinants for RNA binding motif 7 protein stability and RNA recognition

The 7SK ribonucleoprotein (RNP) is a critical regulator of eukaryotic transcription. Recently, RNA binding motif 7 (RBM7) containing an RNA recognition motif (RRM) was reported to associate with 7SK RNA and core 7SK RNP protein components in response to DNA damage. However, little is known about the mode of RBM7-7SK RNA recognition. Here, we found that RRM constructs containing extended C-termini have increased solubility compared to a minimal RRM construct, although these constructs aggregate in a temperature and concentration-dependent manner. Using solution NMR dynamics experiments, we identified additional structural features observed previously in crystal but not in solution structures. To identify potential RBM7-7SK RNA binding sites, we analyzed deposited data from in cellulo crosslinking experiments and found that RBM7 primarily crosslinks to the distal region of 7SK stem-loop 3 (SL3). Electrophoretic mobility shift assays and NMR chemical shift perturbation experiments showed weak binding to 7SK SL3 constructs in vitro. Together, these results provide new insights into RBM7 RRM folding and recognition of 7SK RNA.

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.

Metadensity: a background-aware python pipeline for summarizing CLIP signals on various transcriptomic sites

Motivation

Cross-linking and immunoprecipitation (CLIP) is a technology to map the binding sites of RNA-binding proteins (RBPs). The region where an RBP binds within RNA is often indicative of its molecular function in RNA processing. As an example, the binding sites of splicing factors are found within or proximal to alternatively spliced exons. To better reveal the function of RBPs, we developed a tool to visualize the distribution of CLIP signals around various transcript features.

Results

Here, we present Metadensity (https://github.com/YeoLab/Metadensity), a software that allows users to generate metagene plots. Metadensity allows users to input features such as branchpoints and preserves the near-nucleotide resolution of CLIP technologies by not scaling the features by length. Metadensity normalizes immunoprecipitated libraries with background controls, such as size-matched inputs, then windowing in various user-defined features. Finally, the signals are averaged across a provided set of transcripts.

Availability and implementation

Metadensity is available at https://github.com/YeoLab/Metadensity, with example notebooks at https://metadensity.readthedocs.io/en/latest/tutorial.html.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Genome-wide detection of human variants that disrupt intronic branchpoints

Pre-messenger RNA splicing is initiated with the recognition of a single-nucleotide intronic branchpoint (BP) within a BP motif by spliceosome elements. Forty-eight rare variants in 43 human genes have been reported to alter splicing and cause disease by disrupting BP. However, until now, no computational approach was available to efficiently detect such variants in massively parallel sequencing data. We established a comprehensive human genome-wide BP database by integrating existing BP data and generating new BP data from RNA sequencing of lariat debranching enzyme DBR1-mutated patients and from machine-learning predictions. We characterized multiple features of BP in major and minor introns and found that BP and BP-2 (two nucleotides upstream of BP) positions exhibit a lower rate of variation in human populations and higher evolutionary conservation than the intronic background, while being comparable to the exonic background. We developed BPHunter as a genome-wide computational approach to systematically and efficiently detect intronic variants that may disrupt BP recognition. BPHunter retrospectively identified 40 of the 48 known pathogenic BP variants, in which we summarized a strategy for prioritizing BP variant candidates. The remaining eight variants all create AG-dinucleotides between the BP and acceptor site, which is the likely reason for missplicing. We demonstrated the practical utility of BPHunter prospectively by using it to identify a novel germline heterozygous BP variant of STAT2 in a patient with critical COVID-19 pneumonia and a novel somatic intronic 59-nucleotide deletion of ITPKB in a lymphoma patient, both of which were validated experimentally. BPHunter is publicly available from https://hgidsoft.rockefeller.edu/BPHunter and https://github.com/casanova-lab/BPHunter .

CDK11 regulates pre-mRNA splicing by phosphorylation of SF3B1

RNA splicing, the process of intron removal from pre-mRNA, is essential for the regulation of gene expression. It is controlled by the spliceosome, a megadalton RNA-protein complex that assembles de novo on each pre-mRNA intron through an ordered assembly of intermediate complexes^1,2. Spliceosome activation is a major control step that requires substantial protein and RNA rearrangements leading to a catalytically active complex^1-5. Splicing factor 3B subunit 1 (SF3B1) protein-a subunit of the U2 small nuclear ribonucleoprotein⁶-is phosphorylated during spliceosome activation^7-10, but the kinase that is responsible has not been identified. Here we show that cyclin-dependent kinase 11 (CDK11) associates with SF3B1 and phosphorylates threonine residues at its N terminus during spliceosome activation. The phosphorylation is important for the association between SF3B1 and U5 and U6 snRNAs in the activated spliceosome, termed the B^act complex, and the phosphorylation can be blocked by OTS964, a potent and selective inhibitor of CDK11. Inhibition of CDK11 prevents spliceosomal transition from the precatalytic complex B to the activated complex B^act and leads to widespread intron retention and accumulation of non-functional spliceosomes on pre-mRNAs and chromatin. We demonstrate a central role of CDK11 in spliceosome assembly and splicing regulation and characterize OTS964 as a highly selective CDK11 inhibitor that suppresses spliceosome activation and splicing.

Rapid factor depletion highlights intricacies of nucleoplasmic RNA degradation

Turnover of nucleoplasmic transcripts by the mammalian multi-subunit RNA exosome is mediated by two adaptors: the Nuclear EXosome Targeting (NEXT) complex and the Poly(A) tail eXosome Targeting (PAXT) connection. Functional analyses of NEXT and PAXT have largely utilized long-term factor depletion strategies, facilitating the appearance of indirect phenotypes. Here, we rapidly deplete NEXT, PAXT and core exosome components, uncovering the direct consequences of their acute losses. Generally, proteome changes are sparse and largely dominated by co-depletion of other exosome and adaptor subunits, reflecting possible subcomplex compositions. While parallel high-resolution 3′ end sequencing of newly synthesized RNA confirms previously established factor specificities, it concomitantly demonstrates an inflation of long-term depletion datasets by secondary effects. Most strikingly, a general intron degradation phenotype, observed in long-term NEXT depletion samples, is undetectable upon short-term depletion, which instead emphasizes NEXT targeting of snoRNA-hosting introns. Further analysis of these introns uncovers an unusual mode of core exosome-independent RNA decay. Our study highlights the accumulation of RNAs as an indirect result of long-term decay factor depletion, which we speculate is, at least partly, due to the exhaustion of alternative RNA decay pathways.

Spatial correlation statistics enable transcriptome-wide characterization of RNA structure binding

Molecular interactions at identical transcriptomic locations or at proximal but non-overlapping sites can mediate RNA modification and regulation, necessitating tools to uncover these spatial relationships. We present nearBynding, a flexible algorithm and software pipeline that models spatial correlation between transcriptome-wide tracks from diverse data types. nearBynding can process and correlate interval as well as continuous data and incorporate experimentally derived or in silico predicted transcriptomic tracks. nearBynding offers visualization functions for its statistics to identify colocalizations and adjacent features. We demonstrate the application of nearBynding to correlate RNA-binding protein (RBP) binding preferences with other RBPs, RNA structure, or RNA modification. By cross-correlating RBP binding and RNA structure data, we demonstrate that nearBynding recapitulates known RBP binding to structural motifs and provides biological insights into RBP binding preference of G-quadruplexes. nearBynding is available as an R/Bioconductor package and can run on a personal computer, making correlation of transcriptomic features broadly accessible.

Biology of the mRNA Splicing Machinery and Its Dysregulation in Cancer Providing Therapeutic Opportunities

Dysregulation of messenger RNA (mRNA) processing—in particular mRNA splicing—is a hallmark of cancer. Compared to normal cells, cancer cells frequently present aberrant mRNA splicing, which promotes cancer progression and treatment resistance. This hallmark provides opportunities for developing new targeted cancer treatments. Splicing of precursor mRNA into mature mRNA is executed by a dynamic complex of proteins and small RNAs called the spliceosome. Spliceosomes are part of the supraspliceosome, a macromolecular structure where all co-transcriptional mRNA processing activities in the cell nucleus are coordinated. Here we review the biology of the mRNA splicing machinery in the context of other mRNA processing activities in the supraspliceosome and present current knowledge of its dysregulation in lung cancer. In addition, we review investigations to discover therapeutic targets in the spliceosome and give an overview of inhibitors and modulators of the mRNA splicing process identified so far. Together, this provides insight into the value of targeting the spliceosome as a possible new treatment for lung cancer.

Consequences for Pancreatic β-Cell Identity and Function of Unregulated Transcript Processing

Mounting evidence suggests a role for alternative splicing (AS) of transcripts in the normal physiology and pathophysiology of the pancreatic β-cell. In the apparent absence of RNA repair systems, RNA decay pathways are likely to play an important role in controlling the stability, distribution and diversity of transcript isoforms in these cells. Around 35% of alternatively spliced transcripts in human cells contain premature termination codons (PTCs) and are targeted for degradation via nonsense-mediated decay (NMD), a vital quality control process. Inflammatory cytokines, whose levels are increased in both type 1 (T1D) and type 2 (T2D) diabetes, stimulate alternative splicing events and the expression of NMD components, and may or may not be associated with the activation of the NMD pathway. It is, however, now possible to infer that NMD plays a crucial role in regulating transcript processing in normal and stress conditions in pancreatic β-cells. In this review, we describe the possible role of Regulated Unproductive Splicing and Translation (RUST), a molecular mechanism embracing NMD activity in relationship to AS and translation of damaged transcript isoforms in these cells. This process substantially reduces the abundance of non-functional transcript isoforms, and its dysregulation may be involved in pancreatic β-cell failure in diabetes.

psiCLIP reveals dynamic RNA binding by DEAH-box helicases before and after exon ligation

RNA helicases remodel the spliceosome to enable pre-mRNA splicing, but their binding and mechanism of action remain poorly understood. To define helicase-RNA contacts in specific spliceosomal states, we develop purified spliceosome iCLIP (psiCLIP), which reveals dynamic helicase-RNA contacts during splicing catalysis. The helicase Prp16 binds along the entire available single-stranded RNA region between the branchpoint and 3'-splice site, while Prp22 binds diffusely downstream of the branchpoint before exon ligation, but then switches to more narrow binding in the downstream exon after exon ligation, arguing against a mechanism of processive translocation. Depletion of the exon-ligation factor Prp18 destabilizes Prp22 binding to the pre-mRNA, suggesting that proofreading by Prp22 may sense the stability of the spliceosome during exon ligation. Thus, psiCLIP complements structural studies by providing key insights into the binding and proofreading activity of spliceosomal RNA helicases.

Massively parallel Cas13 screens reveal principles for guide RNA design

Type VI CRISPR enzymes are RNA-targeting proteins with nuclease activity that enable specific and robust target gene knockdown without altering the genome. To define rules for the design of Cas13d guide RNAs (gRNAs), we conducted massively parallel screens targeting messenger RNAs (mRNAs) of a green fluorescent protein transgene, and CD46, CD55 and CD71 cell-surface proteins in human cells. In total, we measured the activity of 24,460 gRNAs with and without mismatches relative to the target sequences. Knockdown efficacy is driven by gRNA-specific features and target site context. Single mismatches generally reduce knockdown to a modest degree, but spacer nucleotides 15-21 are largely intolerant of target site mismatches. We developed a computational model to identify optimal gRNAs and confirm their generalizability, testing 3,979 guides targeting mRNAs of 48 endogenous genes. We show that Cas13 can be used in forward transcriptomic pooled screens and, using our model, predict optimized Cas13 gRNAs for all protein-coding transcripts in the human genome.

Principles of RNA processing from analysis of enhanced CLIP maps for 150 RNA binding proteins

BACKGROUND

A critical step in uncovering rules of RNA processing is to study the in vivo regulatory networks of RNA binding proteins (RBPs). Crosslinking and immunoprecipitation (CLIP) methods enable mapping RBP targets transcriptome-wide, but methodological differences present challenges to large-scale analysis across datasets. The development of enhanced CLIP (eCLIP) enabled the mapping of targets for 150 RBPs in K562 and HepG2, creating a unique resource of RBP interactomes profiled with a standardized methodology in the same cell types.

RESULTS

Our analysis of 223 eCLIP datasets reveals a range of binding modalities, including highly resolved positioning around splicing signals and mRNA untranslated regions that associate with distinct RBP functions. Quantification of enrichment for repetitive and abundant multicopy elements reveals 70% of RBPs have enrichment for non-mRNA element classes, enables identification of novel ribosomal RNA processing factors and sites, and suggests that association with retrotransposable elements reflects multiple RBP mechanisms of action. Analysis of spliceosomal RBPs indicates that eCLIP resolves AQR association after intronic lariat formation, enabling identification of branch points with single-nucleotide resolution, and provides genome-wide validation for a branch point-based scanning model for 3' splice site recognition. Finally, we show that eCLIP peak co-occurrences across RBPs enable the discovery of novel co-interacting RBPs.

CONCLUSIONS

This work reveals novel insights into RNA biology by integrated analysis of eCLIP profiling of 150 RBPs with distinct functions. Further, our quantification of both mRNA and other element association will enable further research to identify novel roles of RBPs in regulating RNA processing.

Alternative Splicing Regulatory Networks: Functions, Mechanisms, and Evolution. Mol Cell. 2019;76:329-345. [PMID: 31626751 DOI: 10.1016/j.molcel.2019.09.017]

High-throughput sequencing-based methods and their applications in the study of transcriptomes have revolutionized our understanding of alternative splicing. Networks of functionally coordinated and biologically important alternative splicing events continue to be discovered in an ever-increasing diversity of cell types in the context of physiologically normal and disease states. These studies have been complemented by efforts directed at defining sequence codes governing splicing and their cognate trans-acting factors, which have illuminated important combinatorial principles of regulation. Additional studies have revealed critical roles of position-dependent, multivalent protein-RNA interactions that direct splicing outcomes. Investigations of evolutionary changes in RNA binding proteins, splice variants, and associated cis elements have further shed light on the emergence, mechanisms, and functions of splicing networks. Progress in these areas has emphasized the need for a coordinated, community-based effort to systematically address the functions of individual splice variants associated with normal and disease biology.

Assessment of branch point prediction tools to predict physiological branch points and their alteration by variants

BACKGROUND

Branch points (BPs) map within short motifs upstream of acceptor splice sites (3'ss) and are essential for splicing of pre-mature mRNA. Several BP-dedicated bioinformatics tools, including HSF, SVM-BPfinder, BPP, Branchpointer, LaBranchoR and RNABPS were developed during the last decade. Here, we evaluated their capability to detect the position of BPs, and also to predict the impact on splicing of variants occurring upstream of 3'ss.

RESULTS

We used a large set of constitutive and alternative human 3'ss collected from Ensembl (n = 264,787 3'ss) and from in-house RNAseq experiments (n = 51,986 3'ss). We also gathered an unprecedented collection of functional splicing data for 120 variants (62 unpublished) occurring in BP areas of disease-causing genes. Branchpointer showed the best performance to detect the relevant BPs upstream of constitutive and alternative 3'ss (99.48 and 65.84% accuracies, respectively). For variants occurring in a BP area, BPP emerged as having the best performance to predict effects on mRNA splicing, with an accuracy of 89.17%.

CONCLUSIONS

Our investigations revealed that Branchpointer was optimal to detect BPs upstream of 3'ss, and that BPP was most relevant to predict splicing alteration due to variants in the BP area.