iCLIP data analysis: A complete pipeline from sequencing reads to RBP binding sites

A high-resolution map of functional miR-181 response elements in the thymus reveals the role of coding sequence targeting and an alternative seed match

MicroRNAs (miRNAs) are critical post-transcriptional regulators in many biological processes. They act by guiding RNA-induced silencing complexes to miRNA response elements (MREs) in target mRNAs, inducing translational inhibition and/or mRNA degradation. Functional MREs are expected to predominantly occur in the 3' untranslated region and involve perfect base-pairing of the miRNA seed. Here, we generate a high-resolution map of miR-181a/b-1 (miR-181) MREs to define the targeting rules of miR-181 in developing murine T cells. By combining a multi-omics approach with computational high-resolution analyses, we uncover novel miR-181 targets and demonstrate that miR-181 acts predominantly through RNA destabilization. Importantly, we discover an alternative seed match and identify a distinct set of targets with repeat elements in the coding sequence which are targeted by miR-181 and mediate translational inhibition. In conclusion, deep profiling of MREs in primary cells is critical to expand physiologically relevant targetomes and establish context-dependent miRNA targeting rules.

Arid5a uses disordered extensions of its core ARID domain for distinct DNA- and RNA-recognition and gene regulation

AT-rich interacting domain (ARID)-containing proteins, Arids, are a heterogeneous DNA-binding protein family involved in transcription regulation and chromatin processing. For the member Arid5a, no exact DNA-binding preference has been experimentally defined so far. Additionally, the protein binds to mRNA motifs for transcript stabilization, supposedly through the DNA-binding ARID domain. To date, however, no unbiased RNA motif definition and clear dissection of nucleic acid-binding through the ARID domain have been undertaken. Using NMR-centered biochemistry, we here define the Arid5a DNA preference. Further, high-throughput in vitro binding reveals a consensus RNA-binding motif engaged by the core ARID domain. Finally, transcriptome-wide binding (iCLIP2) reveals that Arid5a has a weak preference for (A)U-rich regions in pre-mRNA transcripts of factors related to RNA processing. We find that the intrinsically disordered regions flanking the ARID domain modulate the specificity and affinity of DNA binding, while they appear crucial for RNA interactions. Ultimately, our data suggest that Arid5a uses its extended ARID domain for bifunctional gene regulation and that the involvement of IDR extensions is a more general feature of Arids in interacting with different nucleic acids at the chromatin-mRNA interface.

racoon_clip-a complete pipeline for single-nucleotide analyses of iCLIP and eCLIP data

Motivation

A vast variety of biological questions connected to RNA-binding proteins can be tackled with UV crosslinking and immunoprecipitation (CLIP) experiments. However, the processing and analysis of CLIP data are rather complex. Moreover, different types of CLIP experiments like iCLIP or eCLIP are often processed in different ways, reducing comparability between multiple experiments. Therefore, we aimed to build an easy-to-use computational tool for the processing of CLIP data that can be used for both iCLIP and eCLIP data, as well as data from other truncation-based CLIP methods.

Results

Here, we introduce racoon_clip, a sustainable and fully automated pipeline for the complete processing of iCLIP and eCLIP data to extract RNA binding signal at single-nucleotide resolution. racoon_clip is easy to install and execute, with multiple pre-settings and fully customizable parameters, and outputs a conclusive summary report with visualizations and statistics for all analysis steps.

Availability and implementation

racoon_clip is implemented as a Snakemake-powered command line tool (Snakemake version ≥7.22, Python version ≥3.9). The latest release can be downloaded from GitHub (https://github.com/ZarnackGroup/racoon_clip/tree/main) and installed via pip. A detailed documentation, including installation, usage, and customization, can be found at https://racoon-clip.readthedocs.io/en/latest/. The example datasets can be downloaded from the Short Read Archive (SRA; iCLIP: SRR5646576, SRR5646577, SRR5646578) or the ENCODE Project (eCLIP: ENCSR202BFN).

It is not just about transcription: involvement of brain RNA splicing in substance use disorders

Alternative splicing is a co-transcriptional process that significantly contributes to the molecular landscape of the cell. It plays a multifaceted role in shaping gene transcription, protein diversity, and functional adaptability in response to environmental cues. Recent studies demonstrate that drugs of abuse have a profound impact on alternative splicing patterns within different brain regions. Drugs like alcohol and cocaine modify the expression of genes responsible for encoding splicing factors, thereby influencing alternative splicing of crucial genes involved in neurotransmission, neurogenesis, and neuroinflammation. Notable examples of these alterations include alcohol-induced changes in splicing factors such as HSPA6 and PCBP1, as well as cocaine's impact on PTBP1 and SRSF11. Beyond the immediate effects of drug exposure, recent research has shed light on the role of alternative splicing in contributing to the risk of substance use disorders (SUDs). This is exemplified by exon skipping events in key genes like ELOVL7, which can elevate the risk of alcohol use disorder. Lastly, drugs of abuse can induce splicing alterations through epigenetic modifications. For example, cocaine exposure leads to alterations in levels of trimethylated lysine 36 of histone H3, which exhibits a robust association with alternative splicing and serves as a reliable predictor for exon exclusion. In summary, alternative splicing has emerged as a critical player in the complex interplay between drugs of abuse and the brain, offering insights into the molecular underpinnings of SUDs.

Regulation by the RNA-binding protein Unkempt at its effector interface

How RNA-binding proteins (RBPs) convey regulatory instructions to the core effectors of RNA processing is unclear. Here, we document the existence and functions of a multivalent RBP-effector interface. We show that the effector interface of a conserved RBP with an essential role in metazoan development, Unkempt, is mediated by a novel type of 'dual-purpose' peptide motifs that can contact two different surfaces of interacting proteins. Unexpectedly, we find that the multivalent contacts do not merely serve effector recruitment but are required for the accuracy of RNA recognition by Unkempt. Systems analyses reveal that multivalent RBP-effector contacts can repurpose the principal activity of an effector for a different function, as we demonstrate for the reuse of the central eukaryotic mRNA decay factor CCR4-NOT in translational control. Our study establishes the molecular assembly and functional principles of an RBP-effector interface.

Phosphorylation of the compartmentalized PKA substrate TAF15 regulates RNA-protein interactions

Spatiotemporal-controlled second messengers alter molecular interactions of central signaling nodes for ensuring physiological signal transmission. One prototypical second messenger molecule which modulates kinase signal transmission is the cyclic-adenosine monophosphate (cAMP). The main proteinogenic cellular effectors of cAMP are compartmentalized protein kinase A (PKA) complexes. Their cell-type specific compositions precisely coordinate substrate phosphorylation and proper signal propagation which is indispensable for numerous cell-type specific functions. Here we present evidence that TAF15, which is implicated in the etiology of amyotrophic lateral sclerosis, represents a novel nuclear PKA substrate. In cross-linking and immunoprecipitation experiments (iCLIP) we showed that TAF15 phosphorylation alters the binding to target transcripts related to mRNA maturation, splicing and protein-binding related functions. TAF15 appears to be one of multiple PKA substrates that undergo RNA-binding dynamics upon phosphorylation. We observed that the activation of the cAMP-PKA signaling axis caused a change in the composition of a collection of RNA species that interact with TAF15. This observation appears to be a broader principle in the regulation of molecular interactions, as we identified a significant enrichment of RNA-binding proteins within endogenous PKA complexes. We assume that phosphorylation of RNA-binding domains adds another layer of regulation to binary protein-RNAs interactions with consequences to RNA features including binding specificities, localization, abundance and composition.

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 ).

Arabidopsis thaliana GLYCINE RICH RNA-BINDING PROTEIN 7 interaction with its iCLIP target LHCB1.1 correlates with changes in RNA stability and circadian oscillation

The importance of RNA-binding proteins (RBPs) for plant responses to environmental stimuli and development is well documented. Insights into the portfolio of RNAs they recognize, however, clearly lack behind the understanding gathered in non-plant model organisms. Here, we characterize binding of the circadian clock-regulated Arabidopsis thaliana GLYCINE-RICH RNA-BINDING PROTEIN 7 (AtGRP7) to its target transcripts. We identified novel RNA targets from individual-nucleotide resolution UV crosslinking and immunoprecipitation (iCLIP) data using an improved bioinformatics pipeline that will be broadly applicable to plant RBP iCLIP data. 2705 transcripts with binding sites were identified in plants expressing AtGRP7-GFP that were not recovered in plants expressing an RNA-binding dead variant or GFP alone. A conserved RNA motif enriched in uridine residues was identified at the AtGRP7 binding sites. NMR titrations confirmed the preference of AtGRP7 for RNAs with a central U-rich motif. Among the bound RNAs, circadian clock-regulated transcripts were overrepresented. Peak abundance of the LHCB1.1 transcript encoding a chlorophyll-binding protein was reduced in plants overexpressing AtGRP7 whereas it was elevated in atgrp7 mutants, indicating that LHCB1.1 was regulated by AtGRP7 in a dose-dependent manner. In plants overexpressing AtGRP7, the LHCB1.1 half-life was shorter compared to wild-type plants whereas in atgrp7 mutant plants, the half-life was significantly longer. Thus, AtGRP7 modulates circadian oscillations of its in vivo binding target LHCB1.1 by affecting RNA stability.

hGRAD: A versatile "one-fits-all" system to acutely deplete RNA binding proteins from condensates

Nuclear RNA binding proteins (RBPs) are difficult to study because they often belong to large protein families and form extensive networks of auto- and crossregulation. They are highly abundant and many localize to condensates with a slow turnover, requiring long depletion times or knockouts that cannot distinguish between direct and indirect or compensatory effects. Here, we developed a system that is optimized for the rapid degradation of nuclear RBPs, called hGRAD. It comes as a "one-fits-all" plasmid, and integration into any cell line with endogenously GFP-tagged proteins allows for an inducible, rapid, and complete knockdown. We show that the nuclear RBPs SRSF3, SRSF5, SRRM2, and NONO are completely cleared from nuclear speckles and paraspeckles within 2 h. hGRAD works in various cell types, is more efficient than previous methods, and does not require the expression of exogenous ubiquitin ligases. Combining SRSF5 hGRAD degradation with Nascent-seq uncovered transient transcript changes, compensatory mechanisms, and an effect of SRSF5 on transcript stability.

Characterization of RVFV Nucleocapsid Protein Binding Sites on RNA by iCLIP-seq

The nucleocapsid protein (N) in Rift Valley fever virus is an RNA-binding protein that functions in viral transcription, replication, and packaging. In this chapter, the method for studying protein-RNA interactions in context of viral infection using individual nucleotide resolution, cross-linking, immunoprecipitation, and sequencing (iCLIP-seq) is explained. The method is useful for identifying the interactions between both host and viral RNAs with N and can identify RNA motifs that interact with the protein of interest.

K6-linked ubiquitylation marks formaldehyde-induced RNA-protein crosslinks for resolution

Reactive aldehydes are produced by normal cellular metabolism or after alcohol consumption, and they accumulate in human tissues if aldehyde clearance mechanisms are impaired. Their toxicity has been attributed to the damage they cause to genomic DNA and the subsequent inhibition of transcription and replication. However, whether interference with other cellular processes contributes to aldehyde toxicity has not been investigated. We demonstrate that formaldehyde induces RNA-protein crosslinks (RPCs) that stall the ribosome and inhibit translation in human cells. RPCs in the messenger RNA (mRNA) are recognized by the translating ribosomes, marked by atypical K6-linked ubiquitylation catalyzed by the RING-in-between-RING (RBR) E3 ligase RNF14, and subsequently resolved by the ubiquitin- and ATP-dependent unfoldase VCP. Our findings uncover an evolutionary conserved formaldehyde-induced stress response pathway that protects cells against RPC accumulation in the cytoplasm, and they suggest that RPCs contribute to the cellular and tissue toxicity of reactive aldehydes.

MePMe-seq: antibody-free simultaneous m6A and m5C mapping in mRNA by metabolic propargyl labeling and sequencing

Internal modifications of mRNA have emerged as widespread and versatile regulatory mechanism to control gene expression at the post-transcriptional level. Most of these modifications are methyl groups, making S-adenosyl-L-methionine (SAM) a central metabolic hub. Here we show that metabolic labeling with a clickable metabolic precursor of SAM, propargyl-selenohomocysteine (PSH), enables detection and identification of various methylation sites. Propargylated A, C, and G nucleosides form at detectable amounts via intracellular generation of the corresponding SAM analogue. Integration into next generation sequencing enables mapping of N6-methyladenosine (m6A) and 5-methylcytidine (m5C) sites in mRNA with single nucleotide precision (MePMe-seq). Analysis of the termination profiles can be used to distinguish m6A from 2'-O-methyladenosine (Am) and N1-methyladenosine (m1A) sites. MePMe-seq overcomes the problems of antibodies for enrichment and sequence-motifs for evaluation, which was limiting previous methodologies. Metabolic labeling via clickable SAM facilitates the joint evaluation of methylation sites in RNA and potentially DNA and proteins.

The RNA-binding protein RBP42 regulates cellular energy metabolism in mammalian-infective Trypanosoma brucei

RNA-binding proteins (RBPs) are key players in coordinated post-transcriptional regulation of functionally related genes, defined as RNA regulons. RNA regulons play particularly critical roles in parasitic trypanosomes, which exhibit unregulated co-transcription of long unrelated gene arrays. In this report, we present a systematic analysis of an essential RBP, RBP42, in the mammalian-infective bloodstream form of African trypanosome and show that RBP42 is a key regulator of parasite's central carbon and energy metabolism. Using individual-nucleotide resolution UV cross-linking and immunoprecipitation to identify genome-wide RBP42-RNA interactions, we show that RBP42 preferentially binds within the coding region of mRNAs encoding core metabolic enzymes. Global quantitative transcriptomic and proteomic analyses reveal that loss of RBP42 reduces the abundance of target mRNA-encoded proteins, but not target mRNA, suggesting a positive translational regulatory role of RBP42. Significant changes in central carbon metabolic intermediates, following loss of RBP42, further support its critical role in cellular energy metabolism. Trypanosoma brucei infection, transmitted through the bite of blood-feeding tsetse flies, causes deadly diseases in humans and livestock. This disease, if left untreated, is almost always fatal. Existing therapies are toxic and difficult to administer. During T. brucei's lifecycle in two different host environments, the parasite progresses through distinctive life stages with major morphological and metabolic changes, requiring precise alteration of parasite gene expression program. In the absence of regulated transcription, post-transcriptional processes mediated by RNA-binding proteins play critical roles in T. brucei gene regulation. In this study, we show that the RNA-binding protein RBP42 plays crucial roles in cellular energy metabolic regulation of this important human pathogen. Metabolic dysregulation observed in RBP42 knockdown cells offers a breadth of potential interest to researchers studying parasite biology and can also impact research in general eukaryotic biology.

A paradigm for regulation at the effector interface with RNA-binding proteins

RNA-binding proteins (RBPs) are key regulators of gene expression, but how RBPs convey regulatory instructions to the core effectors of RNA processing is unclear. Here we document the existence and functions of a multivalent RBP-effector interface. We show that the effector interface of a deeply conserved RBP with an essential role in metazoan development, Unkempt, is mediated by a novel type of 'dual-purpose' peptide motifs that can contact two different surfaces of interacting proteins. Unexpectedly, we find that the multivalent contacts do not merely serve effector recruitment but are required for the accuracy of RNA recognition by the recruiting RBP. Systems analyses reveal that multivalent RBP-effector contacts can repurpose the principal activity of an effector for a different function, as we demonstrate for reuse of the central eukaryotic mRNA decay factor CCR4-NOT in translational control. Our study establishes the molecular assembly and functional principles of an RBP-effector interface, with implications for the evolution and function of RBP-operated regulatory networks.

Global role of IGF2BP1 in controlling the expression of Wnt/β-catenin-regulated genes in colorectal cancer cells

Introduction: Wnt/β-catenin signaling controls cell division and lineage specification during embryonic development, and is crucial for stem cells maintenance and gut tissue regeneration in adults. Aberrant activation of Wnt/β-catenin signaling is also essential for the pathogenesis of a variety of malignancies. The RNA-binding protein IGF2BP1 is a transcriptional target of Wnt/β-catenin signaling, normally expressed during development and often reactivated in cancer cells, where it regulates the stability of oncogenic mRNA. Methods: In this study, we employed iCLIP and RNA sequencing techniques to investigate the role of IGF2BP1 in the post-transcriptional regulation of Wnt/β-catenin-induced genes at a global level within colorectal cancer (CRC) cells characterized by constitutively active Wnt/β-catenin signaling. Results and Discussion: In our study, we show that, in contrast to normal cells, CRC cells exhibit a much stronger dependency on IGF2BP1 expression for Wnt/β-catenin-regulated genes. We show that both untransformed and CRC cells have their unique subsets of Wnt/β-catenin-regulated genes that IGF2BP1 directly controls through binding to their mRNA. Our iCLIP analysis revealed a significant change in the IGF2BP1-binding sites throughout the target transcriptomes and a significant change in the enrichment of 6-mer motifs associated with IGF2BP1 binding in response to Wnt/β-catenin signaling. Our study also revealed a signature of IGF2BP1-regulated genes that are significantly associated with colon cancer-free survival in humans, as well as potential targets for CRC treatment. Overall, this study highlights the complex and context-dependent regulation of Wnt/β-catenin signaling target genes by IGF2BP1 in non-transformed and CRC cells and identifies potential targets for colon cancer treatment.

METTL1 promotes tumorigenesis through tRNA-derived fragment biogenesis in prostate cancer

Newly growing evidence highlights the essential role that epitranscriptomic marks play in the development of many cancers; however, little is known about the role and implications of altered epitranscriptome deposition in prostate cancer. Here, we show that the transfer RNA N7-methylguanosine (m7G) transferase METTL1 is highly expressed in primary and advanced prostate tumours. Mechanistically, we find that METTL1 depletion causes the loss of m7G tRNA methylation and promotes the biogenesis of a novel class of small non-coding RNAs derived from 5'tRNA fragments. 5'tRNA-derived small RNAs steer translation control to favour the synthesis of key regulators of tumour growth suppression, interferon pathway, and immune effectors. Knockdown of Mettl1 in prostate cancer preclinical models increases intratumoural infiltration of pro-inflammatory immune cells and enhances responses to immunotherapy. Collectively, our findings reveal a therapeutically actionable role of METTL1-directed m7G tRNA methylation in cancer cell translation control and tumour biology.

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.

Investigating the Prevalence of RNA-Binding Metabolic Enzymes in E. coli

An open research field in cellular regulation is the assumed crosstalk between RNAs, metabolic enzymes, and metabolites, also known as the REM hypothesis. High-throughput assays have produced extensive interactome data with metabolic enzymes frequently found as hits, but only a few examples have been biochemically validated, with deficits especially in prokaryotes. Therefore, we rationally selected nineteen Escherichia coli enzymes from such datasets and examined their ability to bind RNAs using two complementary methods, iCLIP and SELEX. Found interactions were validated by EMSA and other methods. For most of the candidates, we observed no RNA binding (12/19) or a rather unspecific binding (5/19). Two of the candidates, namely glutamate-5-kinase (ProB) and quinone oxidoreductase (QorA), displayed specific and previously unknown binding to distinct RNAs. We concentrated on the interaction of QorA to the mRNA of yffO, a grounded prophage gene, which could be validated by EMSA and MST. Because the physiological function of both partners is not known, the biological relevance of this interaction remains elusive. Furthermore, we found novel RNA targets for the MS2 phage coat protein that served us as control. Our results indicate that RNA binding of metabolic enzymes in procaryotes is less frequent than suggested by the results of high-throughput studies, but does occur.

nf-core/clipseq - a robust Nextflow pipeline for comprehensive CLIP data analysis

Crosslinking and immunoprecipitation (CLIP) technologies have become a central component of the molecular biologists' toolkit to study protein-RNA interactions and thus to uncover core principles of RNA biology. There has been a proliferation of CLIP-based experimental protocols, as well as computational tools, especially for peak-calling. Consequently, there is an urgent need for a well-documented bioinformatic pipeline that enshrines the principles of robustness, reproducibility, scalability, portability and flexibility while embracing the diversity of experimental and computational CLIP tools. To address this, we present nf-core/clipseq - a robust Nextflow pipeline for quality control and analysis of CLIP sequencing data. It is part of the international nf-core community effort to develop and curate a best-practice, gold-standard set of pipelines for data analysis. The standards enabled by Nextflow and nf-core, including workflow management, version control, continuous integration and containerisation ensure that these key needs are met. Furthermore, multiple tools are implemented ( e.g. for peak-calling), alongside visualisation of quality control metrics to empower the user to make their own informed decisions based on their data. nf-core/clipseq remains under active development, with plans to incorporate newly released tools to ensure that pipeline remains up-to-date and relevant for the community. Engagement with users and developers is encouraged through the nf-core GitHub repository and Slack channel to promote collaboration. It is available at https://nf-co.re/clipseq.

The U1 antisense morpholino oligonucleotide (AMO) disrupts U1 snRNP structure to promote intronic PCPA modification of pre-mRNAs

Functional depletion of the U1 small nuclear ribonucleoprotein (snRNP) with a 25 nt U1 AMO (antisense morpholino oligonucleotide) may lead to intronic premature cleavage and polyadenylation of thousands of genes, a phenomenon known as U1 snRNP telescripting; however, the underlying mechanism remains elusive. In this study, we demonstrated that U1 AMO could disrupt U1 snRNP structure both in vitro and in vivo, thereby affecting the U1 snRNP-RNAP polymerase II interaction. By performing chromatin immunoprecipitation sequencing for phosphorylation of Ser2 and Ser5 of the C-terminal domain of RPB1, the largest subunit of RNAP polymerase II, we showed that transcription elongation was disturbed upon U1 AMO treatment, with a particular high phosphorylation of Ser2 signal at intronic cryptic polyadenylation sites (PASs). In addition, we showed that core 3'processing factors CPSF/CstF are involved in the processing of intronic cryptic PAS. Their recruitment accumulated toward cryptic PASs upon U1 AMO treatment, as indicated by chromatin immunoprecipitation sequencing and individual-nucleotide resolution CrossLinking and ImmunoPrecipitation sequencing analysis. Conclusively, our data suggest that disruption of U1 snRNP structure mediated by U1 AMO provides a key for understanding the U1 telescripting mechanism.

clipplotr-a comparative visualization and analysis tool for CLIP data

CLIP technologies are now widely used to study RNA-protein interactions and many data sets are now publicly available. An important first step in CLIP data exploration is the visual inspection and assessment of processed genomic data on selected genes or regions and performing comparisons: either across conditions within a particular project, or incorporating publicly available data. However, the output files produced by data processing pipelines or preprocessed files available to download from data repositories are often not suitable for direct comparison and usually need further processing. Furthermore, to derive biological insight it is usually necessary to visualize a CLIP signal alongside other data such as annotations, or orthogonal functional genomic data (e.g., RNA-seq). We have developed a simple, but powerful, command-line tool: clipplotr, which facilitates these visual comparative and integrative analyses with normalization and smoothing options for CLIP data and the ability to show these alongside reference annotation tracks and functional genomic data. These data can be supplied as input to clipplotr in a range of file formats, which will output a publication quality figure. It is written in R and can both run on a laptop computer independently or be integrated into computational workflows on a high-performance cluster. Releases, source code, and documentation are freely available at https://github.com/ulelab/clipplotr.

Depletion of the RNA-binding protein PURA triggers changes in posttranscriptional gene regulation and loss of P-bodies

The RNA-binding protein PURA has been implicated in the rare, monogenetic, neurodevelopmental disorder PURA Syndrome. PURA binds both DNA and RNA and has been associated with various cellular functions. Only little is known about its main cellular roles and the molecular pathways affected upon PURA depletion. Here, we show that PURA is predominantly located in the cytoplasm, where it binds to thousands of mRNAs. Many of these transcripts change abundance in response to PURA depletion. The encoded proteins suggest a role for PURA in immune responses, mitochondrial function, autophagy and processing (P)-body activity. Intriguingly, reduced PURA levels decrease the expression of the integral P-body components LSM14A and DDX6 and strongly affect P-body formation in human cells. Furthermore, PURA knockdown results in stabilization of P-body-enriched transcripts, whereas other mRNAs are not affected. Hence, reduced PURA levels, as reported in patients with PURA Syndrome, influence the formation and composition of this phase-separated RNA processing machinery. Our study proposes PURA Syndrome as a new model to study the tight connection between P-body-associated RNA regulation and neurodevelopmental disorders.

Poison cassette exon splicing of SRSF6 regulates nuclear speckle dispersal and the response to hypoxia

Hypoxia induces massive changes in alternative splicing (AS) to adapt cells to the lack of oxygen. Here, we identify the splicing factor SRSF6 as a key factor in the AS response to hypoxia. The SRSF6 level is strongly reduced in acute hypoxia, which serves a dual purpose: it allows for exon skipping and triggers the dispersal of nuclear speckles. Our data suggest that cells use dispersal of nuclear speckles to reprogram their gene expression during hypoxic adaptation and that SRSF6 plays an important role in cohesion of nuclear speckles. Down-regulation of SRSF6 is achieved through inclusion of a poison cassette exon (PCE) promoted by SRSF4. Removing the PCE 3' splice site using CRISPR/Cas9 abolishes SRSF6 reduction in hypoxia. Aberrantly high SRSF6 levels in hypoxia attenuate hypoxia-mediated AS and impair dispersal of nuclear speckles. As a consequence, proliferation and genomic instability are increased, while the stress response is suppressed. The SRSF4-PCE-SRSF6 hypoxia axis is active in different cancer types, and high SRSF6 expression in hypoxic tumors correlates with a poor prognosis. We propose that the ultra-conserved PCE of SRSF6 acts as a tumor suppressor and that its inclusion in hypoxia is crucial to reduce SRSF6 levels. This may prevent tumor cells from entering the metastatic route of hypoxia adaptation.

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.

Splicing factor BUD31 promotes ovarian cancer progression through sustaining the expression of anti-apoptotic BCL2L12

Dysregulated expression of splicing factors has important roles in cancer development and progression. However, it remains a challenge to identify the cancer-specific splicing variants. Here we demonstrate that spliceosome component BUD31 is increased in ovarian cancer, and its higher expression predicts worse prognosis. We characterize the BUD31-binding motif and find that BUD31 preferentially binds exon-intron regions near splicing sites. Further analysis reveals that BUD31 inhibition results in extensive exon skipping and a reduced production of long isoforms containing full coding sequence. In particular, we identify BCL2L12, an anti-apoptotic BCL2 family member, as one of the functional splicing targets of BUD31. BUD31 stimulates the inclusion of exon 3 to generate full-length BCL2L12 and promotes ovarian cancer progression. Knockdown of BUD31 or splice-switching antisense oligonucleotide treatment promotes exon 3 skipping and results in a truncated isoform of BCL2L12 that undergoes nonsense-mediated mRNA decay, and the cells subsequently undergo apoptosis. Our findings reveal BUD31-regulated exon inclusion as a critical factor for ovarian cancer cell survival and cancer progression.

High-throughput mutagenesis identifies mutations and RNA-binding proteins controlling CD19 splicing and CART-19 therapy resistance

Following CART-19 immunotherapy for B-cell acute lymphoblastic leukaemia (B-ALL), many patients relapse due to loss of the cognate CD19 epitope. Since epitope loss can be caused by aberrant CD19 exon 2 processing, we herein investigate the regulatory code that controls CD19 splicing. We combine high-throughput mutagenesis with mathematical modelling to quantitatively disentangle the effects of all mutations in the region comprising CD19 exons 1-3. Thereupon, we identify ~200 single point mutations that alter CD19 splicing and thus could predispose B-ALL patients to developing CART-19 resistance. Furthermore, we report almost 100 previously unknown splice isoforms that emerge from cryptic splice sites and likely encode non-functional CD19 proteins. We further identify cis-regulatory elements and trans-acting RNA-binding proteins that control CD19 splicing (e.g., PTBP1 and SF3B4) and validate that loss of these factors leads to pervasive CD19 mis-splicing. Our dataset represents a comprehensive resource for identifying predictive biomarkers for CART-19 therapy.

Multiple alternative splicing events in CD19 mRNA have been associated with resistance/relapse to CD19 CAR-T therapy in patients with B cell malignancies. Here, by combining patient data and a high-throughput mutagenesis screen, the authors identify single point mutations and RNA-binding proteins that can control CD19 splicing and be associated with CD19 CAR-T therapy resistance.

The function of Wtap in N6-adenosine methylation of mRNAs controls T cell receptor signaling and survival of T cells

T cell antigen-receptor (TCR) signaling controls the development, activation and survival of T cells by involving several layers and numerous mechanisms of gene regulation. N⁶-methyladenosine (m⁶A) is the most prevalent messenger RNA modification affecting splicing, translation and stability of transcripts. In the present study, we describe the Wtap protein as essential for m⁶A methyltransferase complex function and reveal its crucial role in TCR signaling in mouse T cells. Wtap and m⁶A methyltransferase functions were required for the differentiation of thymocytes, control of activation-induced death of peripheral T cells and prevention of colitis by enabling gut RORγt⁺ regulatory T cell function. Transcriptome and epitranscriptomic analyses reveal that m⁶A modification destabilizes Orai1 and Ripk1 mRNAs. Lack of post-transcriptional repression of the encoded proteins correlated with increased store-operated calcium entry activity and diminished survival of T cells with conditional genetic inactivation of Wtap. These findings uncover how m⁶A modification impacts on TCR signal transduction and determines activation and survival of T cells.

A web server for identifying circRNA-RBP variable-length binding sites based on stacked generalization ensemble deep learning network

Circular RNA (circRNA) can exert biological functions by interacting with RNA-binding protein (RBP), and some deep learning-based methods have been developed to predict RBP binding sites on circRNA. However, most of these methods identify circRNA-RBP binding sites are only based on single data resource and cannot provide exact binding sites, only providing the probability value of a sequence fragment. To solve these problems, we propose a binding sites localization algorithm that fuses binding sites from multiple databases, and further design a stacked generalization ensemble deep learning model named CirRBP to identify RBP binding sites on circRNA. The CirRBP is trained by combining the binding sites from multiple databases and makes predictions by weighted aggregating the predictions of each sub-model. The results show that the CirRBP outperforms any sub-model and existing online prediction model. For better access to our research results, we develop an open-source web application called CRWS (CircRNA-RBP Web Server). Its back-end learning model of the CRWS is a stacked generalization ensemble learning model CirRBP based on different deep learning frameworks. Given a full-length circRNA or fragment sequence and a target RBP, the CRWS can analyze and provide the exact potential binding sites of the target RBP on the given sequence through the binding sites localization algorithm, and visualize it. In addition, the CRWS can discover the most widely distributed motif in each RBP dataset. Up to now, CRWS is the first significant online tool that uses multi-source data to train models and predict exact binding sites. CRWS is now publicly and freely available without login requirement at: http://www.bioinformatics.team.

The XRN1-regulated RNA helicase activity of YTHDC2 ensures mouse fertility independently of m6A recognition

The functional consequence of N⁶-methyladenosine (m⁶A) RNA modification is mediated by "reader" proteins of the YTH family. YTH domain-containing 2 (YTHDC2) is essential for mammalian fertility, but its molecular function is poorly understood. Here, we identify U-rich motifs as binding sites of YTHDC2 on 3' UTRs of mouse testicular RNA targets. Although its YTH domain is an m⁶A-binder in vitro, the YTH point mutant mice are fertile. Significantly, the loss of its 3'→5' RNA helicase activity causes mouse infertility, with the catalytic-dead mutation being dominant negative. Biochemical studies reveal that the weak helicase activity of YTHDC2 is enhanced by its interaction with the 5'→3' exoribonuclease XRN1. Single-cell transcriptomics indicate that Ythdc2 mutant mitotic germ cells transition into meiosis but accumulate a transcriptome with mixed mitotic/meiotic identity that fail to progress further into meiosis. Finally, our demonstration that ythdc2 mutant zebrafish are infertile highlights its conserved role in animal germ cell development.

Rett syndrome linked to defects in forming the MeCP2/Rbfox/LASR complex in mouse models

Rett syndrome (RTT) is a severe neurological disorder and a leading cause of intellectual disability in young females. RTT is mainly caused by mutations found in the X-linked gene encoding methyl-CpG binding protein 2 (MeCP2). Despite extensive studies, the molecular mechanism underlying RTT pathogenesis is still poorly understood. Here, we report MeCP2 as a key subunit of a higher-order multiunit protein complex Rbfox/LASR. Defective MeCP2 in RTT mouse models disrupts the assembly of the MeCP2/Rbfox/LASR complex, leading to reduced binding of Rbfox proteins to target pre-mRNAs and aberrant splicing of Nrxns and Nlgn1 critical for synaptic plasticity. We further show that MeCP2 disease mutants display defective condensate properties and fail to promote phase-separated condensates with Rbfox proteins in vitro and in cultured cells. These data link an impaired function of MeCP2 with disease mutation in splicing control to its defective properties in mediating the higher-order assembly of the MeCP2/Rbfox/LASR complex.

MeCP2 mutations can cause Rett syndrome, a severe childhood neurological disorder. Here the authors show that MeCP2 mediates the higher-order assembly of a large splicing complex Rbfox/LASR, which is disrupted in the mouse models of Rett syndrome.

Deep and accurate detection of m6A RNA modifications using miCLIP2 and m6Aboost machine learning

N6-methyladenosine (m6A) is the most abundant internal RNA modification in eukaryotic mRNAs and influences many aspects of RNA processing. miCLIP (m6A individual-nucleotide resolution UV crosslinking and immunoprecipitation) is an antibody-based approach to map m6A sites with single-nucleotide resolution. However, due to broad antibody reactivity, reliable identification of m6A sites from miCLIP data remains challenging. Here, we present miCLIP2 in combination with machine learning to significantly improve m6A detection. The optimized miCLIP2 results in high-complexity libraries from less input material. Importantly, we established a robust computational pipeline to tackle the inherent issue of false positives in antibody-based m6A detection. The analyses were calibrated with Mettl3 knockout cells to learn the characteristics of m6A deposition, including m6A sites outside of DRACH motifs. To make our results universally applicable, we trained a machine learning model, m6Aboost, based on the experimental and RNA sequence features. Importantly, m6Aboost allows prediction of genuine m6A sites in miCLIP2 data without filtering for DRACH motifs or the need for Mettl3 depletion. Using m6Aboost, we identify thousands of high-confidence m6A sites in different murine and human cell lines, which provide a rich resource for future analysis. Collectively, our combined experimental and computational methodology greatly improves m6A identification.

Ultraplex: A rapid, flexible, all-in-one fastq demultiplexer

Background: The first step of virtually all next generation sequencing analysis involves the splitting of the raw sequencing data into separate files using sample-specific barcodes, a process known as "demultiplexing". However, we found that existing software for this purpose was either too inflexible or too computationally intensive for fast, streamlined processing of raw, single end fastq files containing combinatorial barcodes. Results: Here, we introduce a fast and uniquely flexible demultiplexer, named Ultraplex, which splits a raw FASTQ file containing barcodes either at a single end or at both 5' and 3' ends of reads, trims the sequencing adaptors and low-quality bases, and moves unique molecular identifiers (UMIs) into the read header, allowing subsequent removal of PCR duplicates. Ultraplex is able to perform such single or combinatorial demultiplexing on both single- and paired-end sequencing data, and can process an entire Illumina HiSeq lane, consisting of nearly 500 million reads, in less than 20 minutes. Conclusions: Ultraplex greatly reduces computational burden and pipeline complexity for the demultiplexing of complex sequencing libraries, such as those produced by various CLIP and ribosome profiling protocols, and is also very user friendly, enabling streamlined, robust data processing. Ultraplex is available on PyPi and Conda and via Github.

The RNA-binding profile of the splicing factor SRSF6 in immortalized human pancreatic β-cells

In pancreatic β-cells, the expression of the splicing factor SRSF6 is regulated by GLIS3, a transcription factor encoded by a diabetes susceptibility gene. SRSF6 down-regulation promotes β-cell demise through splicing dysregulation of central genes for β-cells function and survival, but how RNAs are targeted by SRSF6 remains poorly understood. Here, we define the SRSF6 binding landscape in the human pancreatic β-cell line EndoC-βH1 by integrating individual-nucleotide resolution UV cross-linking and immunoprecipitation (iCLIP) under basal conditions with RNA sequencing after SRSF6 knockdown. We detect thousands of SRSF6 bindings sites in coding sequences. Motif analyses suggest that SRSF6 specifically recognizes a purine-rich consensus motif consisting of GAA triplets and that the number of contiguous GAA triplets correlates with increasing binding site strength. The SRSF6 positioning determines the splicing fate. In line with its role in β-cell function, we identify SRSF6 binding sites on regulated exons in several diabetes susceptibility genes. In a proof-of-principle, the splicing of the susceptibility gene LMO7 is modulated by antisense oligonucleotides. Our present study unveils the splicing regulatory landscape of SRSF6 in immortalized human pancreatic β-cells.

I-KCKT allows dissection-free RNA profiling of adult Drosophila intestinal progenitor cells

The adult Drosophila intestinal epithelium is a model system for stem cell biology, but its utility is limited by current biochemical methods that lack cell type resolution. Here, we describe a new proximity-based profiling method that relies upon a GAL4 driver, termed intestinal-kickout-GAL4 (I-KCKT-GAL4), that is exclusively expressed in intestinal progenitor cells. This method uses UV crosslinked whole animal frozen powder as its starting material to immunoprecipitate the RNA cargoes of transgenic epitope-tagged RNA binding proteins driven by I-KCKT-GAL4 When applied to the general mRNA-binder, poly(A)-binding protein, the RNA profile obtained by this method identifies 98.8% of transcripts found after progenitor cell sorting, and has low background noise despite being derived from whole animal lysate. We also mapped the targets of the more selective RNA binder, Fragile X mental retardation protein (FMRP), using enhanced crosslinking and immunoprecipitation (eCLIP), and report for the first time its binding motif in Drosophila cells. This method will therefore enable the RNA profiling of wild-type and mutant intestinal progenitor cells from intact flies exposed to normal and altered environments, as well as the identification of RNA-protein interactions crucial for stem cell function.