Deciphering human ribonucleoprotein regulatory networks

Predicting functional UTR variants by integrating region-specific features

The untranslated region (UTR) of messenger ribonucleic acid (mRNA), including the 5'UTR and 3'UTR, plays a critical role in regulating gene expression and translation. Variants within the UTR can lead to changes associated with human traits and diseases; however, computational prediction of UTR variant effect is challenging. Current noncoding variant prediction mainly focuses on the promoters and enhancers, neglecting the unique sequence of the UTR and thereby limiting their predictive accuracy. In this study, using consolidated datasets of UTR variants from disease databases and large-scale experimental data, we systematically analyzed more than 50 region-specific features of UTR, including functional elements, secondary structure, sequence composition and site conservation. Our analysis reveals that certain features, such as C/G-related sequence composition in 5'UTR and A/T-related sequence composition in 3'UTR, effectively differentiate between nonfunctional and functional variant sets, unveiling potential sequence determinants of functional UTR variants. Leveraging these insights, we developed two classification models to predict functional UTR variants using machine learning, achieving an area under the curve (AUC) value of 0.94 for 5'UTR and 0.85 for 3'UTR, outperforming all existing methods. Our models will be valuable for enhancing clinical interpretation of genetic variants, facilitating the prediction and management of disease risk.

Deep Learning for Elucidating Modifications to RNA-Status and Challenges Ahead

RNA-binding proteins and chemical modifications to RNA play vital roles in the co- and post-transcriptional regulation of genes. In order to fully decipher their biological roles, it is an essential task to catalogue their precise target locations along with their preferred contexts and sequence-based determinants. Recently, deep learning approaches have significantly advanced in this field. These methods can predict the presence or absence of modification at specific genomic regions based on diverse features, particularly sequence and secondary structure, allowing us to decipher the highly non-linear sequence patterns and structures that underlie site preferences. This article provides an overview of how deep learning is being applied to this area, with a particular focus on the problem of mRNA-RBP binding, while also considering other types of chemical modification to RNA. It discusses how different types of model can handle sequence-based and/or secondary-structure-based inputs, the process of model training, including choice of negative regions and separating sets for testing and training, and offers recommendations for developing biologically relevant models. Finally, it highlights four key areas that are crucial for advancing the field.

Oligonucleotide-Based Photoaffinity Probes: Chemical Tools and Applications for Protein Labeling

A variety of proteins interact with DNA and RNA, including polymerases, histones, ribosomes, transcription factors, and repair enzymes. However, the transient non-covalent nature of these interactions poses challenges for analysis. Introducing a covalent bond between proteins and DNA via photochemical activation of a photosensitive functional group introduced onto nucleic acids offers a means to stabilize these often weak interactions without significantly altering the binding interface. Consequently, photoactivatable oligonucleotides are powerful tools for investigating nucleic acid-protein interactions involved in numerous biological and pathological processes. In this review, we provide a comprehensive overview of the chemical tools developed so far and the different strategies used for incorporating the most commonly used photoreactive reagents into oligonucleotide probes or nucleic acids. Furthermore, we illustrate their application with several examples including protein binding site mapping, identification of protein binding partners, and in cell studies.

Improved discovery of RNA-binding protein binding sites in eCLIP data using DEWSeq

Enhanced crosslinking and immunoprecipitation (eCLIP) sequencing is a method for transcriptome-wide detection of binding sites of RNA-binding proteins (RBPs). However, identified crosslink sites can deviate from experimentally established functional elements of even well-studied RBPs. Current peak-calling strategies result in low replication and high false positive rates. Here, we present the R/Bioconductor package DEWSeq that makes use of replicate information and size-matched input controls. We benchmarked DEWSeq on 107 RBPs for which both eCLIP data and RNA sequence motifs are available and were able to more than double the number of motif-containing binding regions relative to standard eCLIP processing. The improvement not only relates to the number of binding sites (3.1-fold with known motifs for RBFOX2), but also their subcellular localization (1.9-fold of mitochondrial genes for FASTKD2) and structural targets (2.2-fold increase of stem-loop regions for SLBP. On several orthogonal CLIP-seq datasets, DEWSeq recovers a larger number of motif-containing binding sites (3.3-fold). DEWSeq is a well-documented R/Bioconductor package, scalable to adequate numbers of replicates, and tends to substantially increase the proportion and total number of RBP binding sites containing biologically relevant features.

A systematic benchmark of machine learning methods for protein-RNA interaction prediction

RNA-binding proteins (RBPs) are central actors of RNA post-transcriptional regulation. Experiments to profile-binding sites of RBPs in vivo are limited to transcripts expressed in the experimental cell type, creating the need for computational methods to infer missing binding information. While numerous machine-learning based methods have been developed for this task, their use of heterogeneous training and evaluation datasets across different sets of RBPs and CLIP-seq protocols makes a direct comparison of their performance difficult. Here, we compile a set of 37 machine learning (primarily deep learning) methods for in vivo RBP-RNA interaction prediction and systematically benchmark a subset of 11 representative methods across hundreds of CLIP-seq datasets and RBPs. Using homogenized sample pre-processing and two negative-class sample generation strategies, we evaluate methods in terms of predictive performance and assess the impact of neural network architectures and input modalities on model performance. We believe that this study will not only enable researchers to choose the optimal prediction method for their tasks at hand, but also aid method developers in developing novel, high-performing methods by introducing a standardized framework for their evaluation.

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.

A computational map of the human-SARS-CoV-2 protein-RNA interactome predicted at single-nucleotide resolution

RNA-binding proteins (RBPs) are critical host factors for viral infection, however, large scale experimental investigation of the binding landscape of human RBPs to viral RNAs is costly and further complicated due to sequence variation between viral strains. To fill this gap, we investigated the role of RBPs in the context of SARS-CoV-2 by constructing the first in silico map of human RBP-viral RNA interactions at nucleotide-resolution using two deep learning methods (pysster and DeepRiPe) trained on data from CLIP-seq experiments on more than 100 human RBPs. We evaluated conservation of RBP binding between six other human pathogenic coronaviruses and identified sites of conserved and differential binding in the UTRs of SARS-CoV-1, SARS-CoV-2 and MERS. We scored the impact of mutations from 11 variants of concern on protein-RNA interaction, identifying a set of gain- and loss-of-binding events, as well as predicted the regulatory impact of putative future mutations. Lastly, we linked RBPs to functional, OMICs and COVID-19 patient data from other studies, and identified MBNL1, FTO and FXR2 RBPs as potential clinical biomarkers. Our results contribute towards a deeper understanding of how viruses hijack host cellular pathways and open new avenues for therapeutic intervention.

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.

RBP Image Database: A resource for the systematic characterization of the subcellular distribution properties of human RNA binding proteins

RNA binding proteins (RBPs) are central regulators of gene expression implicated in all facets of RNA metabolism. As such, they play key roles in cellular physiology and disease etiology. Since different steps of post-transcriptional gene expression tend to occur in specific regions of the cell, including nuclear or cytoplasmic locations, defining the subcellular distribution properties of RBPs is an important step in assessing their potential functions. Here, we present the RBP Image Database, a resource that details the subcellular localization features of 301 RBPs in the human HepG2 and HeLa cell lines, based on the results of systematic immuno-fluorescence studies conducted using a highly validated collection of RBP antibodies and a panel of 12 markers for specific organelles and subcellular structures. The unique features of the RBP Image Database include: (i) hosting of comprehensive representative images for each RBP-marker pair, with ∼250,000 microscopy images; (ii) a manually curated controlled vocabulary of annotation terms detailing the localization features of each factor; and (iii) a user-friendly interface allowing the rapid querying of the data by target or annotation. The RBP Image Database is freely available at https://rnabiology.ircm.qc.ca/RBPImage/.

Identifying interpretable gene-biomarker associations with functionally informed kernel-based tests in 190,000 exomes

Here we present an exome-wide rare genetic variant association study for 30 blood biomarkers in 191,971 individuals in the UK Biobank. We compare gene-based association tests for separate functional variant categories to increase interpretability and identify 193 significant gene-biomarker associations. Genes associated with biomarkers were ~ 4.5-fold enriched for conferring Mendelian disorders. In addition to performing weighted gene-based variant collapsing tests, we design and apply variant-category-specific kernel-based tests that integrate quantitative functional variant effect predictions for missense variants, splicing and the binding of RNA-binding proteins. For these tests, we present a computationally efficient combination of the likelihood-ratio and score tests that found 36% more associations than the score test alone while also controlling the type-1 error. Kernel-based tests identified 13% more associations than their gene-based collapsing counterparts and had advantages in the presence of gain of function missense variants. We introduce local collapsing by amino acid position for missense variants and use it to interpret associations and identify potential novel gain of function variants in PIEZO1. Our results show the benefits of investigating different functional mechanisms when performing rare-variant association tests, and demonstrate pervasive rare-variant contribution to biomarker variability.

Tunable Cellular Localization and Extensive Cytoskeleton-Interplay of Reflectins

Reflectin proteins are natural copolymers consisting of repeated canonical domains. They are located in a biophotonic system called Bragg lamellae and manipulate the dynamic structural coloration of iridocytes. Their biological functions are intriguing, but the underlying mechanism is not fully understood. Reflectin A1, A2, B1, and C were found to present distinguished cyto-/nucleoplasmic localization preferences in the work. Comparable intracellular localization was reproduced by truncated reflectin variants, suggesting a conceivable evolutionary order among reflectin proteins. The size-dependent access of reflectin variants into the nucleus demonstrated a potential model of how reflectins get into Bragg lamellae. Moreover, RfA1 was found to extensively interact with the cytoskeleton, including its binding to actin and enrichment at the microtubule organizing center. This implied that the cytoskeleton system plays a fundamental role during the organization and transportation of reflectin proteins. The findings presented here provide evidence to get an in-depth insight into the evolutionary processes and working mechanisms of reflectins, as well as novel molecular tools to achieve tunable intracellular transportation.

RiPCA: An Assay for the Detection of RNA-Protein Interactions in Live Cells

Increasing interest in studying and modulating the interactions between RNAs and their RNA-binding proteins has indicated the need for enabling technologies. Existing means of detecting RNA-protein interactions (RPIs) are often limited to biochemical or post-lysis methods or cell-based methods that require the addition of an RNA-based affinity tag, such as the MS2 hairpin, precluding them from use in detecting small or highly processed RNAs. Taking advantage of bioorthogonal chemistry- and split-luciferase-based technologies, we developed an assay for the detection of RPIs in live cells. This article details the protocol and design considerations for RiPCA, or RNA interaction with Protein-mediated Complementation Assay. © 2022 Wiley Periodicals LLC.

Targeted RNA editing: novel tools to study post-transcriptional regulation

RNA binding proteins (RBPs) regulate nearly all post-transcriptional processes within cells. To fully understand RBP function, it is essential to identify their in vivo targets. Standard techniques for profiling RBP targets, such as crosslinking immunoprecipitation (CLIP) and its variants, are limited or suboptimal in some situations, e.g. when compatible antibodies are not available and when dealing with small cell populations such as neuronal subtypes and primary stem cells. This review summarizes and compares several genetic approaches recently designed to identify RBP targets in such circumstances. TRIBE (targets of RNA binding proteins identified by editing), RNA tagging, and STAMP (surveying targets by APOBEC-mediated profiling) are new genetic tools useful for the study of post-transcriptional regulation and RBP identification. We describe the underlying RNA base editing technology, recent applications, and therapeutic implications.

A global analysis on the differential regulation of RNA binding proteins (RBPs) by TNF–α as potential modulators of metabolic syndromes

Metabolic syndrome (MetS) is associated with a group of conditions, which enhances the risk of diabetes, heart diseases and stroke in the affected individuals. Earlier reports from our lab have shown that Tumor necrosis factor-α (TNF-α) significantly modulates the expression of 56 genes at the alternative splicing level which are involved in various signaling and metabolic pathways (MetS genes) connected to MetS. These MetS genes were predicted to interact with various RNA-binding proteins (RBPs) when exposed to TNF-α, resulting changes in their alternative splicing patterns. Here we are presenting data of an RNA-Seq analysis, which identified 1218 unique, and significantly regulated genes by TNF-α, 15% of which are RBPs . Among the 1218 genes, 204 genes have been identified as MetS genes by the ingenuity pathway analysis, and 10% of the MetS genes are found as RBPs. Our results also show that TNF-α changes the phosphorylation status of certain RBPs such as SR proteins, crucial players in alternative splicing, possibly via changing the activation status of certain upstream signaling molecules which also act as upstream kinases for these proteins. Taken together, these findings suggest that TNF-α influences the regulation of the RBPs at the various levels for their expression, which may lead to the alteration of the splicing pattern of the MetS genes. MetS genes acting as RBPs and are modulated by TNF-α, predict the existence of highly interconnected mechanisms which require further analysis to understand their dual roles on the onset of these diseases.

Robust T cell activation requires an eIF3-driven burst in T cell receptor translation

Activation of T cells requires a rapid surge in cellular protein synthesis. However, the role of translation initiation in the early induction of specific genes remains unclear. Here, we show human translation initiation factor eIF3 interacts with select immune system related mRNAs including those encoding the T cell receptor (TCR) subunits TCRA and TCRB. Binding of eIF3 to the TCRA and TCRB mRNA 3'-untranslated regions (3'-UTRs) depends on CD28 coreceptor signaling and regulates a burst in TCR translation required for robust T cell activation. Use of the TCRA or TCRB 3'-UTRs to control expression of an anti-CD19 chimeric antigen receptor (CAR) improves the ability of CAR-T cells to kill tumor cells in vitro. These results identify a new mechanism of eIF3-mediated translation control that can aid T cell engineering for immunotherapy applications.

RBPSpot: Learning on appropriate contextual information for RBP binding sites discovery

Identifying the factors determining the RBP-RNA interactions remains a big challenge. It involves sparse binding motifs and a suitable sequence context for binding. The present work describes an approach to detect RBP binding sites in RNAs using an ultra-fast inexact k-mers search for statistically significant seeds. The seeds work as an anchor to evaluate the context and binding potential using flanking region information while leveraging from Deep Feed-forward Neural Network. The developed models also received support from MD-simulation studies. The implemented software, RBPSpot, scored consistently high for all the performance metrics including average accuracy of ∼90% across a large number of validated datasets. It outperformed the compared tools, including some with much complex deep-learning models, during a comprehensive benchmarking process. RBPSpot can identify RBP binding sites in the human system and can also be used to develop new models, making it a valuable resource in the area of regulatory system studies.

•

Efficient motif anchoring helps to get good quality contextual information on binding

•

Realistic and high granularity datasets ensure better performance of the classifiers

•

DNN models on the contextual features outperform more complex deep learning tools

•

RBPSpot algorithm may be used to develop RBP binding models for other species also

Multi-resBind: a residual network-based multi-label classifier for in vivo RNA binding prediction and preference visualization

BACKGROUND

Protein-RNA interactions play key roles in many processes regulating gene expression. To understand the underlying binding preference, ultraviolet cross-linking and immunoprecipitation (CLIP)-based methods have been used to identify the binding sites for hundreds of RNA-binding proteins (RBPs) in vivo. Using these large-scale experimental data to infer RNA binding preference and predict missing binding sites has become a great challenge. Some existing deep-learning models have demonstrated high prediction accuracy for individual RBPs. However, it remains difficult to avoid significant bias due to the experimental protocol. The DeepRiPe method was recently developed to solve this problem via introducing multi-task or multi-label learning into this field. However, this method has not reached an ideal level of prediction power due to the weak neural network architecture.

RESULTS

Compared to the DeepRiPe approach, our Multi-resBind method demonstrated substantial improvements using the same large-scale PAR-CLIP dataset with respect to an increase in the area under the receiver operating characteristic curve and average precision. We conducted extensive experiments to evaluate the impact of various types of input data on the final prediction accuracy. The same approach was used to evaluate the effect of loss functions. Finally, a modified integrated gradient was employed to generate attribution maps. The patterns disentangled from relative contributions according to context offer biological insights into the underlying mechanism of protein-RNA interactions.

CONCLUSIONS

Here, we propose Multi-resBind as a new multi-label deep-learning approach to infer protein-RNA binding preferences and predict novel interactions. The results clearly demonstrate that Multi-resBind is a promising tool to predict unknown binding sites in vivo and gain biology insights into why the neural network makes a given prediction.

Genome-wide functional screen of 3'UTR variants uncovers causal variants for human disease and evolution

3' untranslated region (3'UTR) variants are strongly associated with human traits and diseases, yet few have been causally identified. We developed the massively parallel reporter assay for 3'UTRs (MPRAu) to sensitively assay 12,173 3'UTR variants. We applied MPRAu to six human cell lines, focusing on genetic variants associated with genome-wide association studies (GWAS) and human evolutionary adaptation. MPRAu expands our understanding of 3'UTR function, suggesting that simple sequences predominately explain 3'UTR regulatory activity. We adapt MPRAu to uncover diverse molecular mechanisms at base pair resolution, including an adenylate-uridylate (AU)-rich element of LEPR linked to potential metabolic evolutionary adaptations in East Asians. We nominate hundreds of 3'UTR causal variants with genetically fine-mapped phenotype associations. Using endogenous allelic replacements, we characterize one variant that disrupts a miRNA site regulating the viral defense gene TRIM14 and one that alters PILRB abundance, nominating a causal variant underlying transcriptional changes in age-related macular degeneration.

Matrix-screening reveals a vast potential for direct protein-protein interactions among RNA binding proteins

RNA-binding proteins (RBPs) are crucial factors of post-transcriptional gene regulation and their modes of action are intensely investigated. At the center of attention are RNA motifs that guide where RBPs bind. However, sequence motifs are often poor predictors of RBP-RNA interactions in vivo. It is hence believed that many RBPs recognize RNAs as complexes, to increase specificity and regulatory possibilities. To probe the potential for complex formation among RBPs, we assembled a library of 978 mammalian RBPs and used rec-Y2H matrix screening to detect direct interactions between RBPs, sampling > 600 K interactions. We discovered 1994 new interactions and demonstrate that interacting RBPs bind RNAs adjacently in vivo. We further find that the mRNA binding region and motif preferences of RBPs deviate, depending on their adjacently binding interaction partners. Finally, we reveal novel RBP interaction networks among major RNA processing steps and show that splicing impairing RBP mutations observed in cancer rewire spliceosomal interaction networks. The dataset we provide will be a valuable resource for understanding the combinatorial interactions of RBPs with RNAs and the resulting regulatory outcomes.

Thermodynamic modeling reveals widespread multivalent binding by RNA-binding proteins

MOTIVATION

Understanding how proteins recognize their RNA targets is essential to elucidate regulatory processes in the cell. Many RNA-binding proteins (RBPs) form complexes or have multiple domains that allow them to bind to RNA in a multivalent, cooperative manner. They can thereby achieve higher specificity and affinity than proteins with a single RNA-binding domain. However, current approaches to de novo discovery of RNA binding motifs do not take multivalent binding into account.

RESULTS

We present Bipartite Motif Finder (BMF), which is based on a thermodynamic model of RBPs with two cooperatively binding RNA-binding domains. We show that bivalent binding is a common strategy among RBPs, yielding higher affinity and sequence specificity. We furthermore illustrate that the spatial geometry between the binding sites can be learned from bound RNA sequences. These discovered bipartite motifs are consistent with previously known motifs and binding behaviors. Our results demonstrate the importance of multivalent binding for RNA-binding proteins and highlight the value of bipartite motif models in representing the multivalency of protein-RNA interactions.

AVAILABILITY AND IMPLEMENTATION

BMF source code is available at https://github.com/soedinglab/bipartite_motif_finder under a GPL license. The BMF web server is accessible at https://bmf.soedinglab.org.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

DDX3 depletion represses translation of mRNAs with complex 5' UTRs

DDX3 is an RNA chaperone of the DEAD-box family that regulates translation. Ded1, the yeast ortholog of DDX3, is a global regulator of translation, whereas DDX3 is thought to preferentially affect a subset of mRNAs. However, the set of mRNAs that are regulated by DDX3 are unknown, along with the relationship between DDX3 binding and activity. Here, we use ribosome profiling, RNA-seq, and PAR-CLIP to define the set of mRNAs that are regulated by DDX3 in human cells. We find that while DDX3 binds highly expressed mRNAs, depletion of DDX3 particularly affects the translation of a small subset of the transcriptome. We further find that DDX3 binds a site on helix 16 of the human ribosomal rRNA, placing it immediately adjacent to the mRNA entry channel. Translation changes caused by depleting DDX3 levels or expressing an inactive point mutation are different, consistent with different association of these genetic variant types with disease. Taken together, this work defines the subset of the transcriptome that is responsive to DDX3 inhibition, with relevance for basic biology and disease states where DDX3 is altered.

Transcriptional control of CBX5 by the RNA binding proteins RBMX and RBMXL1 maintains chromatin state in myeloid leukemia

RNA binding proteins (RBPs) are key arbiters of post-transcriptional regulation and are found to be found dysregulated in hematological malignancies. Here, we identify the RBP RBMX and its retrogene RBMXL1 to be required for murine and human myeloid leukemogenesis. RBMX/L1 are overexpressed in acute myeloid leukemia (AML) primary patients compared to healthy individuals, and RBMX/L1 loss delayed leukemia development. RBMX/L1 loss lead to significant changes in chromatin accessibility, as well as chromosomal breaks and gaps. We found that RBMX/L1 directly bind to mRNAs, affect transcription of multiple loci, including CBX5 (HP1α), and control the nascent transcription of the CBX5 locus. Forced CBX5 expression rescued the RBMX/L1 depletion effects on cell growth and apoptosis. Overall, we determine that RBMX/L1 control leukemia cell survival by regulating chromatin state through their downstream target CBX5. These findings identify a mechanism for RBPs directly promoting transcription and suggest RBMX/L1, as well as CBX5, as potential therapeutic targets in myeloid malignancies.

Zooming in on protein-RNA interactions: a multi-level workflow to identify interaction partners

Interactions between proteins and RNA are at the base of numerous cellular regulatory and functional phenomena. The investigation of the biological relevance of non-coding RNAs has led to the identification of numerous novel RNA-binding proteins (RBPs). However, defining the RNA sequences and structures that are selectively recognised by an RBP remains challenging, since these interactions can be transient and highly dynamic, and may be mediated by unstructured regions in the protein, as in the case of many non-canonical RBPs. Numerous experimental and computational methodologies have been developed to predict, identify and verify the binding between a given RBP and potential RNA partners, but navigating across the vast ocean of data can be frustrating and misleading. In this mini-review, we propose a workflow for the identification of the RNA binding partners of putative, newly identified RBPs. The large pool of potential binders selected by in-cell experiments can be enriched by in silico tools such as catRAPID, which is able to predict the RNA sequences more likely to interact with specific RBP regions with high accuracy. The RNA candidates with the highest potential can then be analysed in vitro to determine the binding strength and to precisely identify the binding sites. The results thus obtained can furthermore validate the computational predictions, offering an all-round solution to the issue of finding the most likely RNA binding partners for a newly identified potential RBP.

RNA-binding proteins regulate aldosterone homeostasis in human steroidogenic cells

Angiotensin II (AngII) stimulates adrenocortical cells to produce aldosterone, a master regulator of blood pressure. Despite extensive characterization of the transcriptional and enzymatic control of adrenocortical steroidogenesis, there are still major gaps in the precise regulation of AII-induced gene expression kinetics. Specifically, we do not know the regulatory contribution of RNA-binding proteins (RBPs) and RNA decay, which can control the timing of stimulus-induced gene expression. To investigate this question, we performed a high-resolution RNA-seq time course of the AngII stimulation response and 4-thiouridine pulse labeling in a steroidogenic human cell line (H295R). We identified twelve temporally distinct gene expression responses that contained mRNA encoding proteins known to be important for various steps of aldosterone production, such as cAMP signaling components and steroidogenic enzymes. AngII response kinetics for many of these mRNAs revealed a coordinated increase in both synthesis and decay. These findings were validated in primary human adrenocortical cells stimulated ex vivo with AngII. Using a candidate screen, we identified a subset of RNA-binding protein and RNA decay factors that activate or repress AngII-stimulated aldosterone production. Among the repressors of aldosterone were BTG2, which promotes deadenylation and global RNA decay. BTG2 was induced in response to AngII stimulation and promoted the repression of mRNAs encoding pro-steroidogenic factors indicating the existence of an incoherent feedforward loop controlling aldosterone homeostasis. These data support a model in which coordinated increases in transcription and decay facilitate the major transcriptomic changes required to implement a pro-steroidogenic expression program that actively resolved to prevent aldosterone overproduction.

RBM24 exacerbates bladder cancer progression by forming a Runx1t1/TCF4/miR-625-5p feedback loop

RNA–binding motif protein 24 (RBM24) acts as a multifunctional determinant of cell fate, proliferation, apoptosis, and differentiation during development by regulating premRNA splicing and mRNA stability. It is also implicated in carcinogenesis, but the functions of RBM24 in bladder cancer (BC) remain unclear. In the present study, we revealed that RBM24 was upregulated in BC tissues. Importantly, we found that a higher level of RBM24 was correlated with poor prognosis in BC patients. Overexpression of RBM24 promoted BC cell proliferation, while depletion of RBM24 inhibited BC cell proliferation in vivo and in vitro. Mechanistically, RBM24 positively regulated Runx1t1 expression in BC cells by binding to and enhancing Runx1t1 mRNA stability. Furthermore, Runx1t1 in turn promoted RBM24 expression by interacting with the transcription factor TCF4 and suppressing the transcription of miR-625-5p, which directly targets RBM24 and suppresses RBM24 expression. RBM24-regulated BC cell proliferation was moderated via the Runx1t1/TCF4/miR-625-5p feedback loop. These results indicate that the RBM24/Runx1t1/TCF4/miR-625-5p positive feedback loop participates in BC progression. Disruption of this pathway may be a potential therapeutic strategy for BC treatment.

A protein called RBM24 promotes progression of bladder cancer (BC) by forming a positive feedback loop with a specific transcription factor, driving cancer cell proliferation. Survival rates for BC are low, and the current imperfect understanding of the underlying mechanisms makes it difficult to treat. Ping-Ying Guo at Hebei Medical University in Shijiazhuang, China, and co-workers investigated the role of RBM24, known to be involved in other cancers, and found increased levels in BC tissues. Higher levels were associated with a poor prognosis. Further investigation revealed that RBM24 boosts levels of the transcription factor, which suppresses a molecule that in turn suppresses RBM24, forming a positive feedback loop promoting BC cell proliferation. Interrupting the feedback loop decreased tumor size in a mouse model. These results may help identify better treatments for BC.

ELAVL1 primarily couples mRNA stability with the 3' UTRs of interferon-stimulated genes

Upon pathogen detection, the innate immune system triggers signaling events leading to upregulation of pro-inflammatory and anti-microbial mRNA transcripts. RNA-binding proteins (RBPs) interact with these critical mRNAs and regulate their fates at the post-transcriptional level. One such RBP is ELAVL1. Although significant progress has been made in understanding how embryonic lethal vision-like protein 1 (ELAVL1) regulates mRNAs, its target repertoire and binding distribution within an immunological context remain poorly understood. We overlap four high-throughput approaches to define its context-dependent targets and determine its regulatory impact during immune activation. ELAVL1 transitions from binding overwhelmingly intronic sites to 3′ UTR sites upon immune stimulation of cells, binding previously and newly expressed mRNAs. We find that ELAVL1 mediates the RNA stability of genes that regulate pathways essential to pathogen sensing and cytokine production. Our findings reveal the importance of examining RBP regulatory impact under dynamic transcriptomic events to understand their post-transcriptional regulatory roles within specific biological circuitries.

Rothamel et al. show that upon immune activation, the RNA-binding protein ELAVL1 accumulates in the cytoplasm and redistributes from introns to mRNA 3′ UTRs. 3′ UTR binding confers enrichment and transcript stability. Many top-ranking transcripts are interferon-stimulated genes (ISGs), indicating that ELAVL1 is a positive regulator of an innate immune response.

ELAV/Hu RNA binding proteins determine multiple programs of neural alternative splicing

ELAV/Hu factors are conserved RNA binding proteins (RBPs) that play diverse roles in mRNA processing and regulation. The founding member, Drosophila Elav, was recognized as a vital neural factor 35 years ago. Nevertheless, little was known about its impacts on the transcriptome, and potential functional overlap with its paralogs. Building on our recent findings that neural-specific lengthened 3' UTR isoforms are co-determined by ELAV/Hu factors, we address their impacts on splicing. While only a few splicing targets of Drosophila are known, ectopic expression of each of the three family members (Elav, Fne and Rbp9) alters hundreds of cassette exon and alternative last exon (ALE) splicing choices. Reciprocally, double mutants of elav/fne, but not elav alone, exhibit opposite effects on both classes of regulated mRNA processing events in larval CNS. While manipulation of Drosophila ELAV/Hu RBPs induces both exon skipping and inclusion, characteristic ELAV/Hu motifs are enriched only within introns flanking exons that are suppressed by ELAV/Hu factors. Moreover, the roles of ELAV/Hu factors in global promotion of distal ALE splicing are mechanistically linked to terminal 3' UTR extensions in neurons, since both processes involve bypass of proximal polyadenylation signals linked to ELAV/Hu motifs downstream of cleavage sites. We corroborate the direct action of Elav in diverse modes of mRNA processing using RRM-dependent Elav-CLIP data from S2 cells. Finally, we provide evidence for conservation in mammalian neurons, which undergo broad programs of distal ALE and APA lengthening, linked to ELAV/Hu motifs downstream of regulated polyadenylation sites. Overall, ELAV/Hu RBPs orchestrate multiple broad programs of neuronal mRNA processing and isoform diversification in Drosophila and mammalian neurons.

Translational remodeling by RNA-binding proteins and noncoding RNAs

Responsible for generating the proteome that controls phenotype, translation is the ultimate convergence point for myriad upstream signals that influence gene expression. System-wide adaptive translational reprogramming has recently emerged as a pillar of cellular adaptation. As classic regulators of mRNA stability and translation efficiency, foundational studies established the concept of collaboration and competition between RNA-binding proteins (RBPs) and noncoding RNAs (ncRNAs) on individual mRNAs. Fresh conceptual innovations now highlight stress-activated, evolutionarily conserved RBP networks and ncRNAs that increase the translation efficiency of populations of transcripts encoding proteins that participate in a common cellular process. The discovery of post-transcriptional functions for long noncoding RNAs (lncRNAs) was particularly intriguing given their cell-type-specificity and historical definition as nuclear-functioning epigenetic regulators. The convergence of RBPs, lncRNAs, and microRNAs on functionally related mRNAs to enable adaptive protein synthesis is a newer biological paradigm that highlights their role as "translatome (protein output) remodelers" and reinvigorates the paradigm of "RNA operons." Together, these concepts modernize our understanding of cellular stress adaptation and strategies for therapeutic development. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications Translation > Translation Regulation Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs.

Global analysis of RNA-protein interactions in TNF-α induced alternative splicing in metabolic disorders

In this report, using the database of RNA-binding protein specificities (RBPDB) and our previously published RNA-seq data, we analyzed the interactions between RNA and RNA-binding proteins to decipher the role of alternative splicing in metabolic disorders induced by TNF-α. We identified 13 395 unique RNA-RBP interactions, including 385 unique RNA motifs and 35 RBPs, some of which (including MBNL-1 and 3, ZFP36, ZRANB2, and SNRPA) are transcriptionally regulated by TNF-α. In addition to some previously reported RBPs, such as RBMX and HuR/ELAVL1, we found a few novel RBPs, such as ZRANB2 and SNRPA, to be involved in the regulation of metabolic syndrome-associated genes that contain an enrichment of tetrameric RNA sequences (AUUU). Taken together, this study paves the way for novel RNA-protein interaction-based therapeutics for treating metabolic syndromes.

On the Way to Understanding the Interplay between the RNA Structure and Functions in Cells: A Genome-Wide Perspective

RNAs adopt specific structures in order to perform their biological activities. The structure of RNA is an important layer of gene expression regulation, and can impact a plethora of cellular processes, starting with transcription, RNA processing, and translation, and ending with RNA turnover. The development of high-throughput technologies has enabled a deeper insight into the sophisticated interplay between the structure of the cellular transcriptome and the living cells environment. In this review, we present the current view on the RNA structure in vivo resulting from the most recent transcriptome-wide studies in different organisms, including mammalians, yeast, plants, and bacteria. We focus on the relationship between the mRNA structure and translation, mRNA stability and degradation, protein binding, and RNA posttranscriptional modifications.

High throughput approaches to study RNA-protein interactions in vitro

To understand the regulation of gene expression it is critical to determine how proteins interact with and discriminate between different RNAs. In this review, we discuss experimental techniques that utilize high throughput approaches to characterize the interactions of proteins with large numbers of RNAs in vitro. We describe the underlying principles for the main methods, briefly discuss their scope and limitations, and outline how insight from the techniques contributes to our understanding of specificity for RNA-protein interactions.

Characterization of RNP Networks of PUM1 and PUM2 Post-Transcriptional Regulators in TCam-2 Cells, a Human Male Germ Cell Model

Mammalian Pumilio (PUM) proteins are sequence-specific, RNA-binding proteins (RBPs) with wide-ranging roles. They are involved in germ cell development, which has functional implications in development and fertility. Although human PUM1 and PUM2 are closely related to each other and recognize the same RNA binding motif, there is some evidence for functional diversity. To address that problem, first we used RIP-Seq and RNA-Seq approaches, and identified mRNA pools regulated by PUM1 and PUM2 proteins in the TCam-2 cell line, a human male germ cell model. Second, applying global mass spectrometry-based profiling, we identified distinct PUM1- and PUM2-interacting putative protein cofactors, most of them involved in RNA processing. Third, combinatorial analysis of RIP and RNA-Seq, mass spectrometry, and RNA motif enrichment analysis revealed that PUM1 and PUM2 form partially varied RNP-regulatory networks (RNA regulons), which indicate different roles in human reproduction and testicular tumorigenesis. Altogether, this work proposes that protein paralogues with very similar and evolutionary highly conserved functional domains may play divergent roles in the cell by combining with different sets of protein cofactors. Our findings highlight the versatility of PUM paralogue-based post-transcriptional regulation, offering insight into the mechanisms underlying their diverse biological roles and diseases resulting from their dysfunction.

Deep neural networks for interpreting RNA-binding protein target preferences

Deep learning has become a powerful paradigm to analyze the binding sites of regulatory factors including RNA-binding proteins (RBPs), owing to its strength to learn complex features from possibly multiple sources of raw data. However, the interpretability of these models, which is crucial to improve our understanding of RBP binding preferences and functions, has not yet been investigated in significant detail. We have designed a multitask and multimodal deep neural network for characterizing in vivo RBP targets. The model incorporates not only the sequence but also the region type of the binding sites as input, which helps the model to boost the prediction performance. To interpret the model, we quantified the contribution of the input features to the predictive score of each RBP. Learning across multiple RBPs at once, we are able to avoid experimental biases and to identify the RNA sequence motifs and transcript context patterns that are the most important for the predictions of each individual RBP. Our findings are consistent with known motifs and binding behaviors and can provide new insights about the regulatory functions of RBPs.

The G3BP1-Family-USP10 Deubiquitinase Complex Rescues Ubiquitinated 40S Subunits of Ribosomes Stalled in Translation from Lysosomal Degradation

Ribosome-associated quality control (RQC) purges aberrant mRNAs and nascent polypeptides in a multi-step molecular process initiated by the E3 ligase ZNF598 through sensing of ribosomes collided at aberrant mRNAs and monoubiquitination of distinct small ribosomal subunit proteins. We show that G3BP1-family-USP10 complexes are required for deubiquitination of RPS2, RPS3, and RPS10 to rescue modified 40S subunits from programmed degradation. Knockout of USP10 or G3BP1 family proteins increased lysosomal ribosomal degradation and perturbed ribosomal subunit stoichiometry, both of which were rescued by a single K214R substitution of RPS3. While the majority of RPS2 and RPS3 monoubiquitination resulted from ZNF598-dependent sensing of ribosome collisions initiating RQC, another minor pathway contributed to their monoubiquitination. G3BP1 family proteins have long been considered RNA-binding proteins, however, our results identified 40S subunits and associated mRNAs as their predominant targets, a feature shared by stress granules to which G3BP1 family proteins localize under stress.

Approaches for measuring the dynamics of RNA-protein interactions

RNA-protein interactions are pivotal for the regulation of gene expression from bacteria to human. RNA-protein interactions are dynamic; they change over biologically relevant timescales. Understanding the regulation of gene expression at the RNA level therefore requires knowledge of the dynamics of RNA-protein interactions. Here, we discuss the main experimental approaches to measure dynamic aspects of RNA-protein interactions. We cover techniques that assess dynamics of cellular RNA-protein interactions that accompany biological processes over timescales of hours or longer and techniques measuring the kinetic dynamics of RNA-protein interactions in vitro. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes RNA Evolution and Genomics > Ribonomics.

The Architecture of the Human RNA-Binding Protein Regulatory Network

RNA-binding proteins (RBPs) are key players of post-transcriptional regulation of gene expression, relying on competitive and cooperative interactions to fine-tune their action. Several studies have described individual interactions of RBPs with RBP mRNAs. Here we present a systematic network investigation of fifty thousand interactions between RBPs and the UTRs of RBP mRNAs. We identified two structural features in this network. RBP clusters are groups of densely interconnected RBPs co-binding their targets, suggesting a tight control of cooperative and competitive behaviors. RBP chains are hierarchical structures connecting RBP clusters and driven by evolutionarily ancient RBPs. These features suggest that RBP chains may coordinate the different cell programs controlled by RBP clusters. Under this model, the regulatory signal flows through chains from one cluster to another, implementing elaborate regulatory plans. This work thus suggests RBP-RBP interactions as a backbone driving post-transcriptional regulation of gene expression to control RBPs action on their targets.

•

The RBP-RBP network is a robust and efficient hierarchical structure

•

The RBP-RBP network is formed by RBP clusters and chains

•

RBP-RBP interactions cluster to cooperate and compete on common target mRNAs

•

RBP chains are master regulatory units of the cell led by ancient RBPs