iCLIP predicts the dual splicing effects of TIA-RNA interactions

Decoding protein-RNA interactions using CLIP-based methodologies

Protein-RNA interactions are central to all RNA processing events, with pivotal roles in the regulation of gene expression and cellular functions. Dysregulation of these interactions has been increasingly linked to the pathogenesis of human diseases. High-throughput approaches to identify RNA-binding proteins and their binding sites on RNA - in particular, ultraviolet crosslinking followed by immunoprecipitation (CLIP) - have helped to map the RNA interactome, yielding transcriptome-wide protein-RNA atlases that have contributed to key mechanistic insights into gene expression and gene-regulatory networks. Here, we review these recent advances, explore the effects of cellular context on RNA binding, and discuss how these insights are shaping our understanding of cellular biology. We also review the potential therapeutic applications arising from new knowledge of protein-RNA interactions.

The miR-30-5p/TIA-1 axis directs cellular senescence by regulating mitochondrial dynamics

Senescent cells exhibit a diverse spectrum of changes in their morphology, proliferative capacity, senescence-associated secretory phenotype (SASP) production, and mitochondrial homeostasis. These cells often manifest with elongated mitochondria, a hallmark of cellular senescence. However, the precise regulatory mechanisms orchestrating this phenomenon remain predominantly unexplored. In this study, we provide compelling evidence for decreases in TIA-1, a pivotal regulator of mitochondrial dynamics, in models of both replicative senescence and ionizing radiation (IR)-induced senescence. The downregulation of TIA-1 was determined to trigger mitochondrial elongation and enhance the expression of senescence-associated β-galactosidase, a marker of cellular senescence, in human foreskin fibroblast HS27 cells and human keratinocyte HaCaT cells. Conversely, the overexpression of TIA-1 mitigated IR-induced cellular senescence. Notably, we identified the miR-30-5p family as a novel factor regulating TIA-1 expression. Augmented expression of the miR-30-5p family was responsible for driving mitochondrial elongation and promoting cellular senescence in response to IR. Taken together, our findings underscore the significance of the miR-30-5p/TIA-1 axis in governing mitochondrial dynamics and cellular senescence.

Untacking small RNA profiling and RNA fragment footprinting: Approaches and challenges in library construction

Small RNAs (sRNAs) with sizes ranging from 15 to 50 nucleotides (nt) are critical regulators of gene expression control. Prior studies have shown that sRNAs are involved in a broad range of biological processes, such as organ development, tumorigenesis, and epigenomic regulation; however, emerging evidence unveils a hidden layer of diversity and complexity of endogenously encoded sRNAs profile in eukaryotic organisms, including novel types of sRNAs and the previously unknown post-transcriptional RNA modifications. This underscores the importance for accurate, unbiased detection of sRNAs in various cellular contexts. A multitude of high-throughput methods based on next-generation sequencing (NGS) are developed to decipher the sRNA expression and their modifications. Nonetheless, distinct from mRNA sequencing, the data from sRNA sequencing suffer frequent inconsistencies and high variations emanating from the adapter contaminations and RNA modifications, which overall skew the sRNA libraries. Here, we summarize the sRNA-sequencing approaches, and discuss the considerations and challenges for the strategies and methods of sRNA library construction. The pros and cons of sRNA sequencing have significant implications for implementing RNA fragment footprinting approaches, including CLIP-seq and Ribo-seq. We envision that this review can inspire novel improvements in small RNA sequencing and RNA fragment footprinting in future. This article is categorized under: RNA Evolution and Genomics > Computational Analyses of RNA RNA Processing > Processing of Small RNAs Regulatory RNAs/RNAi/Riboswitches > Biogenesis of Effector Small RNAs.

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

Bibliometric Overview on T-Cell Intracellular Antigens and Their Pathological Implications

T-cell intracellular antigen 1 (TIA1) and TIA1-like/related protein (TIAL1/TIAR) are two members of the classical family of RNA binding proteins. Through their selective interactions with distinct RNAs and proteins, these multifunctional regulators are involved in chromatin remodeling, RNA splicing and processing and translation regulation, linking them to a wide range of diseases including neuronal disorders, cancer and other pathologies. From their discovery to the present day, many studies have focused on the behavior of these proteins in order to understand their impact on molecular and cellular processes and to understand their relationship to human pathologies. The volume of research on these proteins in various fields, including molecular biology, biochemistry, cell biology, immunology and cancer, has steadily increased, indicating a growing interest in these gene expression regulators among researchers. This information can be used to know the most productive institutions working in the field, understand the focus of research, identify key areas of involvement, delve deeper into their relationship and impact on different diseases, and to establish the level of study associated with them.

Multiple roles for AU-rich RNA binding proteins in the development of haematologic malignancies and their resistance to chemotherapy

Post-transcriptional regulation by RNA binding proteins can determine gene expression levels and drive changes in cancer cell proteomes. Identifying mechanisms of protein-RNA binding, including preferred sequence motifs bound in vivo, provides insights into protein-RNA networks and how they impact mRNA structure, function, and stability. In this review, we will focus on proteins that bind to AU-rich elements (AREs) in nascent or mature mRNA where they play roles in response to stresses encountered by cancer cells. ARE-binding proteins (ARE-BPs) specifically impact alternative splicing, stability, decay and translation, and formation of RNA-rich biomolecular condensates like cytoplasmic stress granules (SGs). For example, recent findings highlight the role of ARE-BPs - like TIAR and HUR - in chemotherapy resistance and in translational regulation of mRNAs encoding pro-inflammatory cytokines. We will discuss emerging evidence that different modes of ARE-BP activity impact leukaemia and lymphoma development, progression, adaptation to microenvironment and chemotherapy resistance.

Subcytoplasmic location of translation controls protein output

The cytoplasm is highly compartmentalized, but the extent and consequences of subcytoplasmic mRNA localization in non-polarized cells are largely unknown. We determined mRNA enrichment in TIS granules (TGs) and the rough endoplasmic reticulum (ER) through particle sorting and isolated cytosolic mRNAs by digitonin extraction. When focusing on genes that encode non-membrane proteins, we observed that 52% have transcripts enriched in specific compartments. Compartment enrichment correlates with a combinatorial code based on mRNA length, exon length, and 3' UTR-bound RNA-binding proteins. Compartment-biased mRNAs differ in the functional classes of their encoded proteins: TG-enriched mRNAs encode low-abundance proteins with strong enrichment of transcription factors, whereas ER-enriched mRNAs encode large and highly expressed proteins. Compartment localization is an important determinant of mRNA and protein abundance, which is supported by reporter experiments showing that redirecting cytosolic mRNAs to the ER increases their protein expression. In summary, the cytoplasm is functionally compartmentalized by local translation environments.

Functional characterization of the disease-associated CCL2 rs1024611G-rs13900T haplotype: The role of the RNA-binding protein HuR

CC-chemokine ligand 2 (CCL2) is involved in the pathogenesis of several diseases associated with monocyte/macrophage recruitment, such as HIV-associated neurocognitive disorder (HAND), tuberculosis, and atherosclerosis. The rs1024611 (alleles:A>G; G is the risk allele) polymorphism in the CCL2 cis-regulatory region is associated with increased CCL2 expression in vitro and ex vivo, leukocyte mobilization in vivo, and deleterious disease outcomes. However, the molecular basis for the rs1024611-associated differential CCL2 expression remains poorly characterized. It is conceivable that genetic variant(s) in linkage disequilibrium (LD) with rs1024611 could mediate such effects. Previously, we used rs13900 (alleles:_C>T) in the CCL2 3' untranslated region (3' UTR) that is in perfect LD with rs1024611 to demonstrate allelic expression imbalance (AEI) of CCL2 in heterozygous individuals. Here we tested the hypothesis that the rs13900 could modulate CCL2 expression by altering mRNA turnover and/or translatability. The rs13900 T allele conferred greater stability to the CCL2 transcript when compared to the rs13900 C allele. The rs13900 T allele also had increased binding to Human Antigen R (HuR), an RNA-binding protein, in vitro and ex vivo. The rs13900 alleles imparted differential activity to reporter vectors and influenced the translatability of the reporter transcript. We further demonstrated a role for HuR in mediating allele-specific effects on CCL2 expression in overexpression and silencing studies. The presence of the rs1024611G-rs13900T conferred a distinct transcriptomic signature related to inflammation and immunity. Our studies suggest that the differential interactions of HuR with rs13900 could modulate CCL2 expression and explain the interindividual differences in CCL2-mediated disease susceptibility.

Cooperation and Competition of RNA Secondary Structure and RNA-Protein Interactions in the Regulation of Alternative Splicing

The regulation of alternative splicing in eukaryotic cells is carried out through the coordinated action of a large number of factors, including RNA-binding proteins and RNA structure. The RNA structure influences alternative splicing by blocking cis-regulatory elements, or bringing them closer or farther apart. In combination with RNA-binding proteins, it generates transcript conformations that help to achieve the necessary splicing outcome. However, the binding of regulatory proteins depends on RNA structure and, vice versa, the formation of RNA structure depends on the interaction with regulators. Therefore, RNA structure and RNA-binding proteins are inseparable components of common regulatory mechanisms. This review highlights examples of alternative splicing regulation by RNA-binding proteins, the regulation through local and long-range RNA structures, as well as how these elements work together, cooperate, and compete.

Sensing nucleotide composition in virus RNA

Nucleotide composition plays a crucial role in the structure, function and recognition of RNA molecules. During infection, virus RNA is exposed to multiple endogenous proteins that detect local or global compositional biases and interfere with virus replication. Recent advancements in RNA:protein mapping technologies have enabled the identification of general RNA-binding preferences in the human proteome at basal level and in the context of virus infection. In this review, we explore how cellular proteins recognise nucleotide composition in virus RNA and the impact these interactions have on virus replication. Protein-binding G-rich and C-rich sequences are common examples of how host factors detect and limit infection, and, in contrast, viruses may have evolved to purge their genomes from such motifs. We also give examples of how human RNA-binding proteins inhibit virus replication, not only by destabilising virus RNA, but also by interfering with viral protein translation and genome encapsidation. Understanding the interplay between cellular proteins and virus RNA composition can provide insights into host-virus interactions and uncover potential targets for antiviral strategies.

The RNA binding proteins TIA1 and TIAL1 promote Mcl1 mRNA translation to protect germinal center responses from apoptosis

Germinal centers (GCs) are essential for the establishment of long-lasting antibody responses. GC B cells rely on post-transcriptional RNA mechanisms to translate activation-associated transcriptional programs into functional changes in the cell proteome. However, the critical proteins driving these key mechanisms are still unknown. Here, we show that the RNA binding proteins TIA1 and TIAL1 are required for the generation of long-lasting GC responses. TIA1- and TIAL1-deficient GC B cells fail to undergo antigen-mediated positive selection, expansion and differentiation into B-cell clones producing high-affinity antibodies. Mechanistically, TIA1 and TIAL1 control the transcriptional identity of dark- and light-zone GC B cells and enable timely expression of the prosurvival molecule MCL1. Thus, we demonstrate here that TIA1 and TIAL1 are key players in the post-transcriptional program that selects high-affinity antigen-specific GC B cells.

ORC1 binds to cis-transcribed RNAs for efficient activation of replication origins

Cells must coordinate the activation of thousands of replication origins dispersed throughout their genome. Active transcription is known to favor the formation of mammalian origins, although the role that RNA plays in this process remains unclear. We show that the ORC1 subunit of the human Origin Recognition Complex interacts with RNAs transcribed from genes with origins in their transcription start sites (TSSs), displaying a positive correlation between RNA binding and origin activity. RNA depletion, or the use of ORC1 RNA-binding mutant, result in inefficient activation of proximal origins, linked to impaired ORC1 chromatin release. ORC1 RNA binding activity resides in its intrinsically disordered region, involved in intra- and inter-molecular interactions, regulation by phosphorylation, and phase-separation. We show that RNA binding favors ORC1 chromatin release, by regulating its phosphorylation and subsequent degradation. Our results unveil a non-coding function of RNA as a dynamic component of the chromatin, orchestrating the activation of replication origins.

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.

In vivo PAR-CLIP (viP-CLIP) of liver TIAL1 unveils targets regulating cholesterol synthesis and secretion

System-wide cross-linking and immunoprecipitation (CLIP) approaches have unveiled regulatory mechanisms of RNA-binding proteins (RBPs) mainly in cultured cells due to limitations in the cross-linking efficiency of tissues. Here, we describe viP-CLIP (in vivo PAR-CLIP), a method capable of identifying RBP targets in mammalian tissues, thereby facilitating the functional analysis of RBP-regulatory networks in vivo. We applied viP-CLIP to mouse livers and identified Insig2 and ApoB as prominent TIAL1 target transcripts, indicating an important role of TIAL1 in cholesterol synthesis and secretion. The functional relevance of these targets was confirmed by showing that TIAL1 influences their translation in hepatocytes. Mutant Tial1 mice exhibit altered cholesterol synthesis, APOB secretion and plasma cholesterol levels. Our results demonstrate that viP-CLIP can identify physiologically relevant RBP targets by finding a factor implicated in the negative feedback regulation of cholesterol biosynthesis.

Capture RIC-seq reveals positional rules of PTBP1-associated RNA loops in splicing regulation

RNA-binding proteins (RBPs) bind at different positions of the pre-mRNA molecules to promote or reduce the usage of a particular exon. Seeking to understand the working principle of these positional effects, we develop a capture RIC-seq (CRIC-seq) method to enrich specific RBP-associated in situ proximal RNA-RNA fragments for deep sequencing. We determine hnRNPA1-, SRSF1-, and PTBP1-associated proximal RNA-RNA contacts and regulatory mechanisms in HeLa cells. Unexpectedly, the 3D RNA map analysis shows that PTBP1-associated loops in individual introns preferentially promote cassette exon splicing by accelerating asymmetric intron removal, whereas the loops spanning across cassette exon primarily repress splicing. These "positional rules" can faithfully predict PTBP1-regulated splicing outcomes. We further demonstrate that cancer-related splicing quantitative trait loci can disrupt RNA loops by reducing PTBP1 binding on pre-mRNAs to cause aberrant splicing in tumors. Our study presents a powerful method for exploring the functions of RBP-associated RNA-RNA proximal contacts in gene regulation and disease.

Position-dependent effects of RNA-binding proteins in the context of co-transcriptional splicing

Alternative splicing is an important step in eukaryotic mRNA pre-processing which increases the complexity of gene expression programs, but is frequently altered in disease. Previous work on the regulation of alternative splicing has demonstrated that splicing is controlled by RNA-binding proteins (RBPs) and by epigenetic DNA/histone modifications which affect splicing by changing the speed of polymerase-mediated pre-mRNA transcription. The interplay of these different layers of splicing regulation is poorly understood. In this paper, we derived mathematical models describing how splicing decisions in a three-exon gene are made by combinatorial spliceosome binding to splice sites during ongoing transcription. We additionally take into account the effect of a regulatory RBP and find that the RBP binding position within the sequence is a key determinant of how RNA polymerase velocity affects splicing. Based on these results, we explain paradoxical observations in the experimental literature and further derive rules explaining why the same RBP can act as inhibitor or activator of cassette exon inclusion depending on its binding position. Finally, we derive a stochastic description of co-transcriptional splicing regulation at the single-cell level and show that splicing outcomes show little noise and follow a binomial distribution despite complex regulation by a multitude of factors. Taken together, our simulations demonstrate the robustness of splicing outcomes and reveal that quantitative insights into kinetic competition of co-transcriptional events are required to fully understand this important mechanism of gene expression diversity.

PAR-CLIP Assay in Ferroptosis

Ferroptosis is a regulatory cell death process that is accompanied by large amounts of iron ion accumulation and lipid peroxidation. Photoactivated ribonucleoside-enhanced crosslinking and immunoprecipitation (PAR-CLIP) is a method used to identify the binding sites of RNA-binding proteins (RBPs) on target RNAs with high resolution at the nucleotide level. By inserting photosensitive ribonucleoside analogs into new RNA transcripts of living cells, characteristic mutations can be generated during reverse transcription and be used to accurately locate the crosslinking position of RNAs and RBPs. The use of PAR-CLIP to detect interactions and determine precise crosslinking sites between RNAs and RBPs, or to search for RNAs upstream or downstream of ferroptosis pathways genes through known proteins, can help to clarify and verify the occurrence and regulation mechanisms of the various signaling pathways of ferroptosis. Furthermore, it may reveal new targets for ferroptosis detection and improve the treatment efficiency of ferroptosis-related diseases such as cancer and neurodegenerative diseases. Here, we introduce a specific PAR-CLIP protocol for monitoring the ferroptosis process.

The splicing regulators TIA1 and TIAL1 are required for the expression of the DNA damage repair machinery during B cell lymphopoiesis

B cell lymphopoiesis requires dynamic modulation of the B cell transcriptome for timely coordination of somatic mutagenesis and DNA repair in progenitor B (pro-B) cells. Here, we show that, in pro-B cells, the RNA-binding proteins T cell intracellular antigen 1 (TIA1) and TIA1-like protein (TIAL1) act redundantly to enable developmental progression. They are global splicing regulators that control the expression of hundreds of mRNAs, including those involved in DNA damage repair. Mechanistically, TIA1 and TIAL1 bind to 5' splice sites for exon definition, splicing, and expression of DNA damage sensors, such as Chek2 and Rif1. In their absence, pro-B cells show exacerbated DNA damage, altered P53 expression, and increased cell death. Our study uncovers the importance of tight regulation of RNA splicing by TIA1 and TIAL1 for the expression of integrative transcriptional programs that control DNA damage sensing and repair during B cell development.

Positional motif analysis reveals the extent of specificity of protein-RNA interactions observed by CLIP

BACKGROUND

Crosslinking and immunoprecipitation (CLIP) is a method used to identify in vivo RNA-protein binding sites on a transcriptome-wide scale. With the increasing amounts of available data for RNA-binding proteins (RBPs), it is important to understand to what degree the enriched motifs specify the RNA-binding profiles of RBPs in cells.

RESULTS

We develop positionally enriched k-mer analysis (PEKA), a computational tool for efficient analysis of enriched motifs from individual CLIP datasets, which minimizes the impact of technical and regional genomic biases by internal data normalization. We cross-validate PEKA with mCross and show that the use of input control for background correction is not required to yield high specificity of enriched motifs. We identify motif classes with common enrichment patterns across eCLIP datasets and across RNA regions, while also observing variations in the specificity and the extent of motif enrichment across eCLIP datasets, between variant CLIP protocols, and between CLIP and in vitro binding data. Thereby, we gain insights into the contributions of technical and regional genomic biases to the enriched motifs, and find how motif enrichment features relate to the domain composition and low-complexity regions of the studied proteins.

CONCLUSIONS

Our study provides insights into the overall contributions of regional binding preferences, protein domains, and low-complexity regions to the specificity of protein-RNA interactions, and shows the value of cross-motif and cross-RBP comparison for data interpretation. Our results are presented for exploratory analysis via an online platform in an RBP-centric and motif-centric manner ( https://imaps.goodwright.com/apps/peka/ ).

T-Cell Intracellular Antigen 1-Like Protein in Physiology and Pathology

T-cell intracellular antigen 1 (TIA1)-related/like (TIAR/TIAL1) protein is a multifunctional RNA-binding protein (RBP) involved in regulating many aspects of gene expression, independently or in combination with its paralog TIA1. TIAR was first described in 1992 by Paul Anderson’s lab in relation to the development of a cell death phenotype in immune system cells, as it possesses nucleolytic activity against cytotoxic lymphocyte target cells. Similar to TIA1, it is characterized by a subcellular nucleo-cytoplasmic localization and ubiquitous expression in the cells of different tissues of higher organisms. In this paper, we review the relevant structural and functional information available about TIAR from a triple perspective (molecular, cellular and pathophysiological), paying special attention to its expression and regulation in cellular events and processes linked to human pathophysiology.

Integrative genome-wide analysis reveals EIF3A as a key downstream regulator of translational repressor protein Musashi 2 (MSI2)

Musashi 2 (MSI2) is an RNA binding protein (RBP) that regulates asymmetric cell division and cell fate decisions in normal and cancer stem cells. MSI2 appears to repress translation by binding to 3′ untranslated regions (3′UTRs) of mRNA, but the identity of functional targets remains unknown. Here, we used individual nucleotide resolution cross-linking and immunoprecipitation (iCLIP) to identify direct RNA binding partners of MSI2 and integrated these data with polysome profiling to obtain insights into MSI2 function. iCLIP revealed specific MSI2 binding to thousands of mRNAs largely in 3′UTRs, but translational differences were restricted to a small fraction of these transcripts, indicating that MSI2 regulation is not triggered by simple binding. Instead, the functional targets identified here were bound at higher density and contain more ‘UAG’ motifs compared to targets bound nonproductively. To further distinguish direct and indirect targets, MSI2 was acutely depleted. Surprisingly, only 50 transcripts were found to undergo translational induction on acute loss. Using complementary approaches, we determined eukaryotic translation initiation factor 3A (EIF3A) to be an immediate, direct target. We propose that MSI2 downregulation of EIF3A amplifies these effects on translation. Our results also underscore the challenges in defining functional targets of RBPs since mere binding does not imply a discernible functional interaction.

Next-generation sequencing: A new avenue to understand viral RNA-protein interactions

The genomes of RNA viruses present an astonishing source of both sequence and structural diversity. From intracellular viral RNA-host interfaces to interactions between the RNA genome and structural proteins in virus particles themselves, almost the entire viral lifecycle is accompanied by a myriad of RNA-protein interactions that are required to fulfill their replicative potential. It is therefore important to characterize such rich and dynamic collections of viral RNA-protein interactions to understand virus evolution and their adaptation to their hosts and environment. Recent advances in next-generation sequencing technologies have allowed the characterization of viral RNA-protein interactions, including both transient and conserved interactions, where molecular and structural approaches have fallen short. In this review, we will provide a methodological overview of the high-throughput techniques used to study viral RNA-protein interactions, their biochemical mechanisms, and how they evolved from classical methods as well as one another. We will discuss how different techniques have fueled virus research to characterize how viral RNA and proteins interact, both locally and on a global scale. Finally, we will present examples on how these techniques influence the studies of clinically important pathogens such as HIV-1 and SARS-CoV-2.

Prediction of RNA-protein interactions using a nucleotide language model

Motivation

The accumulation of sequencing data has enabled researchers to predict the interactions between RNA sequences and RNA-binding proteins (RBPs) using novel machine learning techniques. However, existing models are often difficult to interpret and require additional information to sequences. Bidirectional encoder representations from transformer (BERT) is a language-based deep learning model that is highly interpretable. Therefore, a model based on BERT architecture can potentially overcome such limitations.

Results

Here, we propose BERT-RBP as a model to predict RNA-RBP interactions by adapting the BERT architecture pretrained on a human reference genome. Our model outperformed state-of-the-art prediction models using the eCLIP-seq data of 154 RBPs. The detailed analysis further revealed that BERT-RBP could recognize both the transcript region type and RNA secondary structure only based on sequence information. Overall, the results provide insights into the fine-tuning mechanism of BERT in biological contexts and provide evidence of the applicability of the model to other RNA-related problems.

Availability and implementation

Python source codes are freely available at https://github.com/kkyamada/bert-rbp. The datasets underlying this article were derived from sources in the public domain: [RBPsuite (http://www.csbio.sjtu.edu.cn/bioinf/RBPsuite/), Ensembl Biomart (http://asia.ensembl.org/biomart/martview/)].

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Dynamics of T-Cell Intracellular Antigen 1-Dependent Stress Granules in Proteostasis and Welander Distal Myopathy under Oxidative Stress

T-cell intracellular antigen 1 (TIA1) is an RNA-binding protein that is primarily involved in the post-transcriptional regulation of cellular RNAs. Furthermore, it is a key component of stress granules (SGs), RNA, and protein aggregates that are formed in response to stressful stimuli to reduce cellular activity as a survival mechanism. TIA1 p.E384K mutation is the genetic cause of Welander distal myopathy (WDM), a late-onset muscular dystrophy whose pathogenesis has been related to modifying SG dynamics. In this study, we present the results obtained by analyzing two specific aspects: (i) SGs properties and dynamics depending on the amino acid at position 384 of TIA1; and (ii) the formation/disassembly time-course of TIA1WT/WDM-dependent SGs under oxidative stress. The generation of TIA1 variants—in which the amino acid mutated in WDM and the adjacent ones were replaced by lysines, glutamic acids, or alanines—allowed us to verify that the inclusion of a single lysine is necessary and sufficient to alter SGs dynamics. Moreover, time-lapse microscopy analysis allowed us to establish in vivo the dynamics of TIA1WT/WDM-dependent SG formation and disassembly, after the elimination of the oxidizing agent, for 1 and 3 h, respectively. Our observations show distinct dynamics between the formation and disassembly of TIA1WT/WDM-dependent SGs. Taken together, this study has allowed us to expand the existing knowledge on the role of TIA1 and the WDM mutation in SG formation.

In Situ Peroxidase Labeling Followed by Mass-Spectrometry Reveals TIA1 Interactome

TIA1 is a broadly expressed DNA/RNA binding protein that regulates multiple aspects of RNA metabolism. It is best known for its role in stress granule assembly during the cellular stress response. Three RNA recognition motifs mediate TIA1 functions along with a prion-like domain that supports multivalent protein-protein interactions that are yet poorly characterized. Here, by fusing the enhanced ascorbate peroxidase 2 (APEX2) biotin-labeling enzyme to TIA1 combined with mass spectrometry, the proteins in the immediate vicinity of TIA1 were defined in situ. Eighty-six and 203 protein partners, mostly associated with ribonucleoprotein complexes, were identified in unstressed control and acute stress conditions, respectively. Remarkably, the repertoire of TIA1 protein partners was highly dissimilar between the two cellular states. Under unstressed control conditions, the biological processes associated with the TIA1 interactome were enriched for cytosolic ontologies related to mRNA metabolism, such as translation initiation, nucleocytoplasmic transport, and RNA catabolism, while the protein identities were primarily represented by RNA binding proteins, ribosomal subunits, and eicosanoid regulators. Under acute stress, TIA1-labeled partners displayed a broader subcellular distribution that included the chromosomes and mitochondria. The enriched biological processes included splicing, translation, and protein synthesis regulation, while the molecular function of the proteins was enriched for RNA binding activity, ribosomal subunits, DNA double-strand break repair, and amide metabolism. Altogether, these data highlight the TIA1 spatial environment with its different partners in diverse cellular states and pave the way to dissect TIA1 role in these processes.

The Multifunctional Faces of T-Cell Intracellular Antigen 1 in Health and Disease

T-cell intracellular antigen 1 (TIA1) is an RNA-binding protein that is expressed in many tissues and in the vast majority of species, although it was first discovered as a component of human cytotoxic T lymphocytes. TIA1 has a dual localization in the nucleus and cytoplasm, where it plays an important role as a regulator of gene-expression flux. As a multifunctional master modulator, TIA1 controls biological processes relevant to the physiological functioning of the organism and the development and/or progression of several human pathologies. This review summarizes our current knowledge of the molecular aspects and cellular processes involving TIA1, with relevance for human pathophysiology.

Deficiency of T-Cell Intracellular Antigen 1 in Murine Embryonic Fibroblasts Is Associated with Changes in Mitochondrial Morphology and Respiration

T-cell intracellular antigen 1 (TIA1) is a multifunctional RNA-binding protein involved in regulating gene expression and splicing during development and in response to environmental stress, to maintain cell homeostasis and promote survival. Herein, we used TIA1-deficient murine embryonic fibroblasts (MEFs) to study their role in mitochondria homeostasis. We found that the loss of TIA1 was associated with changes in mitochondrial morphology, promoting the appearance of elongated mitochondria with heterogeneous cristae density and size. The proteomic patterns of TIA1-deficient MEFs were consistent with expression changes in molecular components related to mitochondrial dynamics/organization and respiration. Bioenergetics analysis illustrated that TIA1 deficiency enhances mitochondrial respiration. Overall, our findings shed light on the role of TIA1 in mitochondrial dynamics and highlight a point of crosstalk between potential pro-survival and pro-senescence pathways.

DHX15-independent roles for TFIP11 in U6 snRNA modification, U4/U6.U5 tri-snRNP assembly and pre-mRNA splicing fidelity

The U6 snRNA, the core catalytic component of the spliceosome, is extensively modified post-transcriptionally, with 2'-O-methylation being most common. However, how U6 2'-O-methylation is regulated remains largely unknown. Here we report that TFIP11, the human homolog of the yeast spliceosome disassembly factor Ntr1, localizes to nucleoli and Cajal Bodies and is essential for the 2'-O-methylation of U6. Mechanistically, we demonstrate that TFIP11 knockdown reduces the association of U6 snRNA with fibrillarin and associated snoRNAs, therefore altering U6 2'-O-methylation. We show U6 snRNA hypomethylation is associated with changes in assembly of the U4/U6.U5 tri-snRNP leading to defects in spliceosome assembly and alterations in splicing fidelity. Strikingly, this function of TFIP11 is independent of the RNA helicase DHX15, its known partner in yeast. In sum, our study demonstrates an unrecognized function for TFIP11 in U6 snRNP modification and U4/U6.U5 tri-snRNP assembly, identifying TFIP11 as a critical spliceosome assembly regulator.

Improved Predicting of The Sequence Specificities of RNA Binding Proteins by Deep Learning

RNA-binding proteins (RBPs) have a significant role in various regulatory tasks. However, the mechanism by which RBPs identify the subsequence target RNAs is still not clear. In recent years, several machine and deep learning-based computational models have been proposed for understanding the binding preferences of RBPs. These methods required integrating multiple features with raw RNA sequences such as secondary structure and their performances can be further improved. In this paper, we propose an efficient and simple convolution neural network, RBPCNN, that relies on the combination of the raw RNA sequence and evolutionary information. We show that conservation scores (evolutionary information) for the RNA sequences can significantly improve the overall performance of the proposed predictor. In addition, the automatic extraction of the binding sequence motifs can enhance our understanding of the binding specificities of RBPs. The experimental results show that RBPCNN outperforms significantly the current state-of-the-art methods. More specifically, the average area under the receiver operator curve was improved by 2.67 percent and the mean average precision was improved by 8.03 percent. The datasets and results can be downloaded from https://home.jbnu.ac.kr/NSCL/RBPCNN.htm.

Identification of TIA1 mRNA targets during human neuronal development

Background

Neuronal development is a tightly controlled process involving multi-layered regulatory mechanisms. While transcriptional pathways regulating neurodevelopment are well characterized, post-transcriptional programs are still poorly understood. TIA1 is an RNA-binding protein that can regulate splicing, stability, or translation of target mRNAs, and has been shown to play critical roles in stress response and neurodevelopment. However, the identity of mRNAs regulated by TIA1 during neurodevelopment under unstressed conditions is still unknown.

Methods and Results

To identify the mRNAs targeted by TIA1 during the first stages of human neurodevelopment, we performed RNA immunoprecipitation-sequencing (RIP-seq) on human embryonic stem cells (hESCs) and derived neural progenitor cells (NPCs), and cortical neurons under unstressed conditions. While there was no change in TIA1 protein levels, the number of TIA1 targeted mRNAs decreased from pluripotent cells to neurons. We identified 2400, 845, and 330 TIA1 mRNA targets in hESCs, NPC, and neurons, respectively. The vast majority of mRNA targets in hESC were genes associated with neurodevelopment and included autism spectrum disorder-risk genes that were not bound in neurons. Additionally, we found that most TIA1 mRNA targets have reduced ribosomal engagement levels.

Conclusion

Our results reveal TIA1 mRNA targets in hESCs and during human neurodevelopment, indicate that translation repression is a key process targeted by TIA1 binding and implicate TIA1 function in neuronal differentiation.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11033-021-06634-0.

The SRSF4-GAS5-Glucocorticoid Receptor Axis Regulates Ventricular Hypertrophy

Supplemental Digital Content is available in the text.

RBPs (RNA-binding proteins) play critical roles in human biology and disease. Aberrant RBP expression affects various steps in RNA processing, altering the function of the target RNAs. The RBP SRSF4 (serine/arginine-rich splicing factor 4) has been linked to neuropathies and cancer. However, its role in the heart is completely unknown.

Dynamics of alternative splicing during somatic cell reprogramming reveals functions for RNA-binding proteins CPSF3, hnRNP UL1, and TIA1

BACKGROUND

Somatic cell reprogramming is the process that allows differentiated cells to revert to a pluripotent state. In contrast to the extensively studied rewiring of epigenetic and transcriptional programs required for reprogramming, the dynamics of post-transcriptional changes and their associated regulatory mechanisms remain poorly understood. Here we study the dynamics of alternative splicing changes occurring during efficient reprogramming of mouse B cells into induced pluripotent stem (iPS) cells and compare them to those occurring during reprogramming of mouse embryonic fibroblasts.

RESULTS

We observe a significant overlap between alternative splicing changes detected in the two reprogramming systems, which are generally uncoupled from changes in transcriptional levels. Correlation between gene expression of potential regulators and specific clusters of alternative splicing changes enables the identification and subsequent validation of CPSF3 and hnRNP UL1 as facilitators, and TIA1 as repressor of mouse embryonic fibroblasts reprogramming. We further find that these RNA-binding proteins control partially overlapping programs of splicing regulation, involving genes relevant for developmental and morphogenetic processes.

CONCLUSIONS

Our results reveal common programs of splicing regulation during reprogramming of different cell types and identify three novel regulators of this process and their targets.

Transite: A Computational Motif-Based Analysis Platform That Identifies RNA-Binding Proteins Modulating Changes in Gene Expression

RNA-binding proteins (RBPs) play critical roles in regulating gene expression by modulating splicing, RNA stability, and protein translation. Stimulus-induced alterations in RBP function contribute to global changes in gene expression, but identifying which RBPs are responsible for the observed changes remains an unmet need. Here, we present Transite, a computational approach that systematically infers RBPs influencing gene expression through changes in RNA stability and degradation. As a proof of principle, we apply Transite to RNA expression data from human patients with non-small-cell lung cancer whose tumors were sampled at diagnosis or after recurrence following treatment with platinum-based chemotherapy. Transite implicates known RBP regulators of the DNA damage response and identifies hnRNPC as a new modulator of chemotherapeutic resistance, which we subsequently validated experimentally. Transite serves as a framework for the identification of RBPs that drive cell-state transitions and adds additional value to the vast collection of publicly available gene expression datasets.

Krismer et al. present a computational approach to identify RNA-binding proteins (RBPs) that modulate post-transcriptional control of gene expression using RNA expression data as inputs. By applying this approach to publicly available patient datasets, they identify and experimentally confirm that the RBP hnRNPC contributes to chemotherapy resistance in lung cancer.

SRSF3 and SRSF7 modulate 3'UTR length through suppression or activation of proximal polyadenylation sites and regulation of CFIm levels

Background

Alternative polyadenylation (APA) refers to the regulated selection of polyadenylation sites (PASs) in transcripts, which determines the length of their 3′ untranslated regions (3′UTRs). We have recently shown that SRSF3 and SRSF7, two closely related SR proteins, connect APA with mRNA export. The mechanism underlying APA regulation by SRSF3 and SRSF7 remained unknown.

Results

Here we combine iCLIP and 3′-end sequencing and find that SRSF3 and SRSF7 bind upstream of proximal PASs (pPASs), but they exert opposite effects on 3′UTR length. SRSF7 enhances pPAS usage in a concentration-dependent but splicing-independent manner by recruiting the cleavage factor FIP1, generating short 3′UTRs. Protein domains unique to SRSF7, which are absent from SRSF3, contribute to FIP1 recruitment. In contrast, SRSF3 promotes distal PAS (dPAS) usage and hence long 3′UTRs directly by counteracting SRSF7, but also indirectly by maintaining high levels of cleavage factor Im (CFIm) via alternative splicing. Upon SRSF3 depletion, CFIm levels decrease and 3′UTRs are shortened. The indirect SRSF3 targets are particularly sensitive to low CFIm levels, because here CFIm serves a dual function; it enhances dPAS and inhibits pPAS usage by binding immediately downstream and assembling unproductive cleavage complexes, which together promotes long 3′UTRs.

Conclusions

We demonstrate that SRSF3 and SRSF7 are direct modulators of pPAS usage and show how small differences in the domain architecture of SR proteins can confer opposite effects on pPAS regulation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-021-02298-y.

DeepCLIP: predicting the effect of mutations on protein-RNA binding with deep learning

Nucleotide variants can cause functional changes by altering protein-RNA binding in various ways that are not easy to predict. This can affect processes such as splicing, nuclear shuttling, and stability of the transcript. Therefore, correct modeling of protein-RNA binding is critical when predicting the effects of sequence variations. Many RNA-binding proteins recognize a diverse set of motifs and binding is typically also dependent on the genomic context, making this task particularly challenging. Here, we present DeepCLIP, the first method for context-aware modeling and predicting protein binding to RNA nucleic acids using exclusively sequence data as input. We show that DeepCLIP outperforms existing methods for modeling RNA-protein binding. Importantly, we demonstrate that DeepCLIP predictions correlate with the functional outcomes of nucleotide variants in independent wet lab experiments. Furthermore, we show how DeepCLIP binding profiles can be used in the design of therapeutically relevant antisense oligonucleotides, and to uncover possible position-dependent regulation in a tissue-specific manner. DeepCLIP is freely available as a stand-alone application and as a webtool at http://deepclip.compbio.sdu.dk.

RNA-Protein Binding Sites Prediction via Multi Scale Convolutional Gated Recurrent Unit Networks

RNA-Protein binding plays important roles in the field of gene expression. With the development of high throughput sequencing, several conventional methods and deep learning-based methods have been proposed to predict the binding preference of RNA-protein binding. These methods can hardly meet the need of consideration of the dependencies between subsequence and the various motif lengths of different translation factors (TFs). To overcome such limitations, we propose a predictive model that utilizes a combination of multi-scale convolutional layers and bidirectional gated recurrent unit (GRU) layer. Multi-scale convolution layer has the ability to capture the motif features of different lengths, and bidirectional GRU layer is able to capture the dependencies among subsequence. Experimental results show that the proposed method performs better than four state-of-the-art methods in this field. In addition, we investigate the effect of model structure on model performance by performing our proposed method with a different convolution layer and a different number of kernel size. We also demonstrate the effectiveness of bidirectional GRU in improving model performance through comparative experiments.

The Transcriptome-wide Landscape and Modalities of EJC Binding in Adult Drosophila

Exon junction complex (EJC) assembles after splicing at specific positions upstream of exon-exon junctions in mRNAs of all higher eukaryotes, affecting major regulatory events. In mammalian cell cytoplasm, EJC is essential for efficient RNA surveillance, while in Drosophila, EJC is essential for localization of oskar mRNA. Here we developed a method for isolation of protein complexes and associated RNA targets (ipaRt) to explore the EJC RNA-binding landscape in a transcriptome-wide manner in adult Drosophila. We find the EJC at canonical positions, preferably on mRNAs from genes comprising multiple splice sites and long introns. Moreover, EJC occupancy is highest at junctions adjacent to strong splice sites, CG-rich hexamers, and RNA structures. Highly occupied mRNAs tend to be maternally localized and derive from genes involved in differentiation or development. These modalities, which have not been reported in mammals, specify EJC assembly on a biologically coherent set of transcripts in Drosophila.

Plant Individual Nucleotide Resolution Cross-Linking and Immunoprecipitation to Characterize RNA-Protein Complexes

In recent years, it has become increasingly recognized that regulation at the RNA level pervasively shapes the transcriptome in eukaryotic cells. This has fostered an interest in the mode of action of RNA-binding proteins that, via interaction with specific RNA sequence motifs, modulate gene expression. Understanding such posttranscriptional networks controlled by an RNA-binding protein requires a comprehensive identification of its in vivo targets. In metazoans and yeast, methods have been devised to stabilize RNA-protein interactions by UV cross-linking before isolating RNA-protein complexes using antibodies, followed by identification of associated RNAs by next-generation sequencing. These methods are collectively referred to as CLIP-Seq (cross-linking immunoprecipitation-high-throughput sequencing). Here, we present a version of the individual nucleotide resolution cross-linking and immunoprecipitation procedure that is suitable for use in the model plant Arabidopsis thaliana.

Intrinsic Regulatory Role of RNA Structural Arrangement in Alternative Splicing Control

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

A large-scale binding and functional map of human RNA-binding proteins

Many proteins regulate the expression of genes by binding to specific regions encoded in the genome1. Here we introduce a new data set of RNA elements in the human genome that are recognized by RNA-binding proteins (RBPs), generated as part of the Encyclopedia of DNA Elements (ENCODE) project phase III. This class of regulatory elements functions only when transcribed into RNA, as they serve as the binding sites for RBPs that control post-transcriptional processes such as splicing, cleavage and polyadenylation, and the editing, localization, stability and translation of mRNAs. We describe the mapping and characterization of RNA elements recognized by a large collection of human RBPs in K562 and HepG2 cells. Integrative analyses using five assays identify RBP binding sites on RNA and chromatin in vivo, the in vitro binding preferences of RBPs, the function of RBP binding sites and the subcellular localization of RBPs, producing 1,223 replicated data sets for 356 RBPs. We describe the spectrum of RBP binding throughout the transcriptome and the connections between these interactions and various aspects of RNA biology, including RNA stability, splicing regulation and RNA localization. These data expand the catalogue of functional elements encoded in the human genome by the addition of a large set of elements that function at the RNA level by interacting with RBPs.

SRSF7 maintains its homeostasis through the expression of Split-ORFs and nuclear body assembly

SRSF7 is an essential RNA-binding protein whose misexpression promotes cancer. Here, we describe how SRSF7 maintains its protein homeostasis in murine P19 cells using an intricate negative feedback mechanism. SRSF7 binding to its premessenger RNA promotes inclusion of a poison cassette exon and transcript degradation via nonsense-mediated decay (NMD). However, elevated SRSF7 levels inhibit NMD and promote translation of two protein halves, termed Split-ORFs, from the bicistronic SRSF7-PCE transcript. The first half acts as dominant-negative isoform suppressing poison cassette exon inclusion and instead promoting the retention of flanking introns containing repeated SRSF7 binding sites. Massive SRSF7 binding to these sites and its oligomerization promote the assembly of large nuclear bodies, which sequester SRSF7 transcripts at their transcription site, preventing their export and restoring normal SRSF7 protein levels. We further show that hundreds of human and mouse NMD targets, especially RNA-binding proteins, encode potential Split-ORFs, some of which are expressed under specific cellular conditions.

Alternative Splicing Regulatory Networks: Functions, Mechanisms, and Evolution

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

Dysregulation of RNA Splicing in Tauopathies

Pathological aggregation of RNA binding proteins (RBPs) is associated with dysregulation of RNA splicing in PS19 P301S tau transgenic mice and in Alzheimer's disease brain tissues. The dysregulated splicing particularly affects genes involved in synaptic transmission. The effects of neuroprotective TIA1 reduction on PS19 mice are also examined. TIA1 reduction reduces disease-linked alternative splicing events for the major synaptic mRNA transcripts examined, suggesting that normalization of RBP functions is associated with the neuroprotection. Use of the NetDecoder informatics algorithm identifies key upstream biological targets, including MYC and EGFR, underlying the transcriptional and splicing changes in the protected compared to tauopathy mice. Pharmacological inhibition of MYC and EGFR activity in neuronal cultures tau recapitulates the neuroprotective effects of TIA1 reduction. These results demonstrate that dysfunction of RBPs and RNA splicing processes are major elements of the pathophysiology of tauopathies, as well as potential therapeutic targets for tauopathies.

FICC-Seq: a method for enzyme-specified profiling of methyl-5-uridine in cellular RNA

Methyl-5-uridine (m5U) is one the most abundant non-canonical bases present in cellular RNA, and in yeast is found at position U54 of tRNAs where modification is catalysed by the methyltransferase Trm2. Although the mammalian enzymes that catalyse m5U formation are yet to be identified via experimental evidence, based on sequence homology to Trm2, two candidates currently exist, TRMT2A and TRMT2B. Here we developed a genome-wide single-nucleotide resolution mapping method, Fluorouracil-Induced-Catalytic-Crosslinking-Sequencing (FICC-Seq), in order to identify the relevant enzymatic targets. We demonstrate that TRMT2A is responsible for the majority of m5U present in human RNA, and that it commonly targets U54 of cytosolic tRNAs. By comparison to current methods, we show that FICC-Seq is a particularly robust method for accurate and reliable detection of relevant enzymatic target sites. Our associated finding of extensive irreversible TRMT2A-tRNA crosslinking in vivo following 5-Fluorouracil exposure is also intriguing, as it suggests a tangible mechanism for a previously suspected RNA-dependent route of Fluorouracil-mediated cytotoxicity.

Investigation of RNA Editing Sites within Bound Regions of RNA-Binding Proteins

Studies in epitranscriptomics indicate that RNA is modified by a variety of enzymes. Among these RNA modifications, adenosine to inosine (A-to-I) RNA editing occurs frequently in the mammalian transcriptome. These RNA editing sites can be detected directly from RNA sequencing (RNA-seq) data by examining nucleotide changes from adenosine (A) to guanine (G), which substitutes for inosine (I). However, a careful investigation of such nucleotide changes must be conducted to distinguish sequencing errors and genomic mutations from the genuine editing sites. Building upon our recent introduction of an easy-to-use bioinformatics tool, RNA Editor, to detect RNA editing events from RNA-seq data, we examined the extent by which RNA editing events affect the binding of RNA-binding proteins (RBP). Through employing bioinformatic techniques, we uncovered that RNA editing sites occur frequently in RBP-bound regions. Moreover, the presence of RNA editing sites are more frequent when RNA editing islands were examined, which are regions in which RNA editing sites are present in clusters. When the binding of one RBP, human antigen R [HuR; encoded by ELAV-like protein 1 (ELAV1)], was quantified experimentally, its binding was reduced upon silencing of the RNA editing enzyme adenosine deaminases acting on RNA (ADAR) compared to the control-suggesting that the presence of RNA editing islands influence HuR binding to its target regions. These data indicate RNA editing as an important mediator of RBP-RNA interactions-a mechanism which likely constitutes an additional mode of post-transcription gene regulation in biological systems.

The RNA-binding ubiquitin ligase MKRN1 functions in ribosome-associated quality control of poly(A) translation

BACKGROUND

Cells have evolved quality control mechanisms to ensure protein homeostasis by detecting and degrading aberrant mRNAs and proteins. A common source of aberrant mRNAs is premature polyadenylation, which can result in non-functional protein products. Translating ribosomes that encounter poly(A) sequences are terminally stalled, followed by ribosome recycling and decay of the truncated nascent polypeptide via ribosome-associated quality control.

RESULTS

Here, we demonstrate that the conserved RNA-binding E3 ubiquitin ligase Makorin Ring Finger Protein 1 (MKRN1) promotes ribosome stalling at poly(A) sequences during ribosome-associated quality control. We show that MKRN1 directly binds to the cytoplasmic poly(A)-binding protein (PABPC1) and associates with polysomes. MKRN1 is positioned upstream of poly(A) tails in mRNAs in a PABPC1-dependent manner. Ubiquitin remnant profiling and in vitro ubiquitylation assays uncover PABPC1 and ribosomal protein RPS10 as direct ubiquitylation substrates of MKRN1.

CONCLUSIONS

We propose that MKRN1 mediates the recognition of poly(A) tails to prevent the production of erroneous proteins from prematurely polyadenylated transcripts, thereby maintaining proteome integrity.

SMARTIV: combined sequence and structure de-novo motif discovery for in-vivo RNA binding data

Gene expression regulation is highly dependent on binding of RNA-binding proteins (RBPs) to their RNA targets. Growing evidence supports the notion that both RNA primary sequence and its local secondary structure play a role in specific Protein-RNA recognition and binding. Despite the great advance in high-throughput experimental methods for identifying sequence targets of RBPs, predicting the specific sequence and structure binding preferences of RBPs remains a major challenge. We present a novel webserver, SMARTIV, designed for discovering and visualizing combined RNA sequence and structure motifs from high-throughput RNA-binding data, generated from in-vivo experiments. The uniqueness of SMARTIV is that it predicts motifs from enriched k-mers that combine information from ranked RNA sequences and their predicted secondary structure, obtained using various folding methods. Consequently, SMARTIV generates Position Weight Matrices (PWMs) in a combined sequence and structure alphabet with assigned P-values. SMARTIV concisely represents the sequence and structure motif content as a single graphical logo, which is informative and easy for visual perception. SMARTIV was examined extensively on a variety of high-throughput binding experiments for RBPs from different families, generated from different technologies, showing consistent and accurate results. Finally, SMARTIV is a user-friendly webserver, highly efficient in run-time and freely accessible via http://smartiv.technion.ac.il/.

PUF60-activated exons uncover altered 3' splice-site selection by germline missense mutations in a single RRM

PUF60 is a splicing factor that binds uridine (U)-rich tracts and facilitates association of the U2 small nuclear ribonucleoprotein with primary transcripts. PUF60 deficiency (PD) causes a developmental delay coupled with intellectual disability and spinal, cardiac, ocular and renal defects, but PD pathogenesis is not understood. Using RNA-Seq, we identify human PUF60-regulated exons and show that PUF60 preferentially acts as their activator. PUF60-activated internal exons are enriched for Us upstream of their 3′ splice sites (3′ss), are preceded by longer AG dinucleotide exclusion zones and more distant branch sites, with a higher probability of unpaired interactions across a typical branch site location as compared to control exons. In contrast, PUF60-repressed exons show U-depletion with lower estimates of RNA single-strandedness. We also describe PUF60-regulated, alternatively spliced isoforms encoding other U-bound splicing factors, including PUF60 partners, suggesting that they are co-regulated in the cell, and identify PUF60-regulated exons derived from transposed elements. PD-associated amino-acid substitutions, even within a single RNA recognition motif (RRM), altered selection of competing 3′ss and branch points of a PUF60-dependent exon and the 3′ss choice was also influenced by alternative splicing of PUF60. Finally, we propose that differential distribution of RNA processing steps detected in cells lacking PUF60 and the PUF60-paralog RBM39 is due to the RBM39 RS domain interactions. Together, these results provide new insights into regulation of exon usage by the 3′ss organization and reveal that germline mutation heterogeneity in RRMs can enhance phenotypic variability at the level of splice-site and branch-site selection.

How RNA structure dictates the usage of a critical exon of spinal muscular atrophy gene

Role of RNA structure in pre-mRNA splicing has been implicated for several critical exons associated with genetic disorders. However, much of the structural studies linked to pre-mRNA splicing regulation are limited to terminal stem-loop structures (hairpins) sequestering splice sites. In few instances, role of long-distance interactions is implicated as the major determinant of splicing regulation. With the recent surge of reports of circular RNA (circRNAs) generated by backsplicing, role of Alu-associated RNA structures formed by long-range interactions are taking central stage. Humans contain two nearly identical copies of Survival Motor Neuron (SMN) genes, SMN1 and SMN2. Deletion or mutation of SMN1 coupled with the inability of SMN2 to compensate for the loss of SMN1 due to exon 7 skipping causes spinal muscular atrophy (SMA), one of the leading genetic diseases of children. In this review, we describe how structural elements formed by both local and long-distance interactions are being exploited to modulate SMN2 exon 7 splicing as a potential therapy for SMA. We also discuss how Alu-associated secondary structure modulates generation of a vast repertoire of SMN circRNAs. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.