Insights into the design and interpretation of iCLIP experiments

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

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

Serpine1 mRNA confers mesenchymal characteristics to the cell and promotes CD8+ T cells exclusion from colon adenocarcinomas

Serine protease inhibitor clade E member 1 (SERPINE1) inhibits extracellular matrix proteolysis and cell detachment. However, SERPINE1 expression also promotes tumor progression and plays a crucial role in metastasis. Here, we solve this apparent paradox and report that Serpine1 mRNA per se, independent of its protein-coding function, confers mesenchymal properties to the cell, promoting migration, invasiveness, and resistance to anoikis and increasing glycolytic activity by sequestering miRNAs. Expression of Serpine1 mRNA upregulates the expression of the TRA2B splicing factor without affecting its mRNA levels. Through transcriptional profiling, we found that Serpine1 mRNA expression downregulates through TRA2B the expression of genes involved in the immune response. Analysis of human colon tumor samples showed an inverse correlation between SERPINE1 mRNA expression and CD8+ T cell infiltration, unveiling the potential value of SERPINE1 mRNA as a promising therapeutic target for colon tumors.

Comprehensive mapping of exon junction complex binding sites reveals universal EJC deposition in Drosophila

BACKGROUND

The exon junction complex (EJC) is involved in most steps of the mRNA life cycle, ranging from splicing to nonsense-mediated mRNA decay (NMD). It is assembled by the splicing machinery onto mRNA in a sequence-independent manner. A fundamental open question is whether the EJC is deposited onto all exon‒exon junctions or only on a subset of them. Several previous studies have made observations supportive of the latter, yet these have been limited by methodological constraints.

RESULTS

In this study, we sought to overcome these limitations via the integration of two different approaches for transcriptome-wide mapping of EJCs. Our results revealed that nearly all, if not all, internal exons consistently harbor an EJC in Drosophila, demonstrating that EJC presence is an inherent consequence of the splicing reaction. Furthermore, our study underscores the limitations of eCLIP methods in fully elucidating the landscape of RBP binding sites. Our findings highlight how highly specific (low false positive) methodologies can lead to erroneous interpretations due to partial sensitivity (high false negatives).

CONCLUSIONS

This study contributes to our understanding of EJC deposition and its association with pre-mRNA splicing. The universal presence of EJC on internal exons underscores its significance in ensuring proper mRNA processing. Additionally, our observations highlight the need to consider both specificity and sensitivity in RBP mapping methodologies.

Optimized infrared photoactivatable ribonucleoside-enhanced crosslinking and immunoprecipitation (IR-PAR-CLIP) protocol identifies novel IGF2BP3-interacting RNAs in colon cancer cells

The conserved family of RNA-binding proteins (RBPs), IGF2BPs, plays an essential role in posttranscriptional regulation controlling mRNA stability, localization, and translation. Mammalian cells express three isoforms of IGF2BPs: IGF2BP1-3. IGF2BP3 is highly overexpressed in cancer cells, and its expression correlates with a poor prognosis in various tumors. Therefore, revealing its target RNAs with high specificity in healthy tissues and in cancer cells is of crucial importance. Photoactivatable-ribonucleoside-enhanced crosslinking and immunoprecipitation (PAR-CLIP) identifies the binding sites of RBPs on their target RNAs at nucleotide resolution in a transcriptome-wide manner. Here, we optimized the PAR-CLIP protocol to study RNA targets of endogenous IGF2BP3 in a human colorectal carcinoma cell line. To this end, we first established an immunoprecipitation protocol to obtain highly pure endogenous IGF2BP3-RNA complexes. Second, we modified the protocol to use highly sensitive infrared (IR) fluorescent dyes instead of radioactive probes to visualize IGF2BP3-crosslinked RNAs. We named the modified method "IR-PAR-CLIP." Third, we compared RNase cleavage conditions and found that sequence preferences of the RNases impact the number of the identified IGF2BP3 targets and introduce a systematic bias in the identified RNA motifs. Fourth, we adapted the single adapter circular ligation approach to increase the efficiency in library preparation. The optimized IR-PAR-CLIP protocol revealed novel RNA targets of IGF2BP3 in a human colorectal carcinoma cell line. We anticipate that our IR-PAR-CLIP approach provides a framework for studies of other RBPs.

Improved analysis of (e)CLIP data with RCRUNCH yields a compendium of RNA-binding protein binding sites and motifs

We present RCRUNCH, an end-to-end solution to CLIP data analysis for identification of binding sites and sequence specificity of RNA-binding proteins. RCRUNCH can analyze not only reads that map uniquely to the genome but also those that map to multiple genome locations or across splice boundaries and can consider various types of background in the estimation of read enrichment. By applying RCRUNCH to the eCLIP data from the ENCODE project, we have constructed a comprehensive and homogeneous resource of in-vivo-bound RBP sequence motifs. RCRUNCH automates the reproducible analysis of CLIP data, enabling studies of post-transcriptional control of gene expression.

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.

Liver RBFOX2 regulates cholesterol homeostasis via Scarb1 alternative splicing in mice

RNA alternative splicing (AS) expands the regulatory potential of eukaryotic genomes. The mechanisms regulating liver-specific AS profiles and their contribution to liver function are poorly understood. Here, we identify a key role for the splicing factor RNA-binding Fox protein 2 (RBFOX2) in maintaining cholesterol homeostasis in a lipogenic environment in the liver. Using enhanced individual-nucleotide-resolution ultra-violet cross-linking and immunoprecipitation, we identify physiologically relevant targets of RBFOX2 in mouse liver, including the scavenger receptor class B type I (Scarb1). RBFOX2 function is decreased in the liver in diet-induced obesity, causing a Scarb1 isoform switch and alteration of hepatocyte lipid homeostasis. Our findings demonstrate that specific AS programmes actively maintain liver physiology, and underlie the lipotoxic effects of obesogenic diets when dysregulated. Splice-switching oligonucleotides targeting this network alleviate obesity-induced inflammation in the liver and promote an anti-atherogenic lipoprotein profile in the blood, underscoring the potential of isoform-specific RNA therapeutics for treating metabolism-associated diseases.

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

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

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

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

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02755-2.

Transcriptome-wide identification of RNA-binding protein binding sites using seCLIP-seq

Discovery of interaction sites between RNA-binding proteins (RBPs) and their RNA targets plays a critical role in enabling our understanding of how these RBPs control RNA processing and regulation. Cross-linking and immunoprecipitation (CLIP) provides a generalizable, transcriptome-wide method by which RBP/RNA complexes are purified and sequenced to identify sites of intermolecular contact. By simplifying technical challenges in prior CLIP methods and incorporating the generation of and quantitative comparison against size-matched input controls, the single-end enhanced CLIP (seCLIP) protocol allows for the profiling of these interactions with high resolution, efficiency and scalability. Here, we present a step-by-step guide to the seCLIP method, detailing critical steps and offering insights regarding troubleshooting and expected results while carrying out the ~4-d protocol. Furthermore, we describe a comprehensive bioinformatics pipeline that offers users the tools necessary to process two replicate datasets and identify reproducible and significant peaks for an RBP of interest in ~2 d.

TDP-43 condensation properties specify its RNA-binding and regulatory repertoire

Mutations causing amyotrophic lateral sclerosis (ALS) often affect the condensation properties of RNA-binding proteins (RBPs). However, the role of RBP condensation in the specificity and function of protein-RNA complexes remains unclear. We created a series of TDP-43 C-terminal domain (CTD) variants that exhibited a gradient of low to high condensation propensity, as observed in vitro and by nuclear mobility and foci formation. Notably, a capacity for condensation was required for efficient TDP-43 assembly on subsets of RNA-binding regions, which contain unusually long clusters of motifs of characteristic types and density. These "binding-region condensates" are promoted by homomeric CTD-driven interactions and required for efficient regulation of a subset of bound transcripts, including autoregulation of TDP-43 mRNA. We establish that RBP condensation can occur in a binding-region-specific manner to selectively modulate transcriptome-wide RNA regulation, which has implications for remodeling RNA networks in the context of signaling, disease, and evolution.

The RNA regulatory programs that govern lymphocyte development and function

Lymphocytes require of constant and dynamic changes in their transcriptome for timely activation and production of effector molecules to combat external pathogens. Synthesis and translation of messenger (m)RNAs into these effector proteins is controlled both quantitatively and qualitatively by RNA binding proteins (RBPs). RBP-dependent regulation of RNA editing, subcellular location, stability, and translation shapes immune cell development and immunity. Extensive evidences have now been gathered from few model RBPs, HuR, PTBP1, ZFP36, and Roquin. However, recently developed methodologies for global characterization of protein:RNA interactions suggest the existence of complex RNA regulatory networks in which RBPs co-ordinately regulate the fate of sets of RNAs controlling cellular pathways and functions. In turn, RNA can also act as scaffolding of functionally related proteins modulating their activation and function. Here we review current knowledge about how RBP-dependent regulation of RNA shapes our immune system and discuss about the existence of a hidden immune cell epitranscriptome. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications.

Recognizing binding sites of poorly characterized RNA-binding proteins on circular RNAs using attention Siamese network

Circular RNAs (circRNAs) interact with RNA-binding proteins (RBPs) to play crucial roles in gene regulation and disease development. Computational approaches have attracted much attention to quickly predict highly potential RBP binding sites on circRNAs using the sequence or structure statistical binding knowledge. Deep learning is one of the popular learning models in this area but usually requires a lot of labeled training data. It would perform unsatisfactorily for the less characterized RBPs with a limited number of known target circRNAs. How to improve the prediction performance for such small-size labeled characterized RBPs is a challenging task for deep learning-based models. In this study, we propose an RBP-specific method iDeepC for predicting RBP binding sites on circRNAs from sequences. It adopts a Siamese neural network consisting of a lightweight attention module and a metric module. We have found that Siamese neural network effectively enhances the network capability of capturing mutual information between circRNAs with pairwise metric learning. To further deal with the small-sample size problem, we have performed the pretraining using available labeled data from other RBPs and also demonstrate the efficacy of this transfer-learning pipeline. We comprehensively evaluated iDeepC on the benchmark datasets of RBP-binding circRNAs, and the results suggest iDeepC achieving promising results on the poorly characterized RBPs. The source code is available at https://github.com/hehew321/iDeepC.

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.

Respiratory syncytial virus M2-1 protein associates non-specifically with viral messenger RNA and with specific cellular messenger RNA transcripts

Respiratory syncytial virus (RSV) is a major cause of respiratory disease in infants and the elderly. RSV is a non-segmented negative strand RNA virus. The viral M2-1 protein plays a key role in viral transcription, serving as an elongation factor to enable synthesis of full-length mRNAs. M2-1 contains an unusual CCCH zinc-finger motif that is conserved in the related human metapneumovirus M2-1 protein and filovirus VP30 proteins. Previous biochemical studies have suggested that RSV M2-1 might bind to specific virus RNA sequences, such as the transcription gene end signals or poly A tails, but there was no clear consensus on what RSV sequences it binds. To determine if M2-1 binds to specific RSV RNA sequences during infection, we mapped points of M2-1:RNA interactions in RSV-infected cells at 8 and 18 hours post infection using crosslinking immunoprecipitation with RNA sequencing (CLIP-Seq). This analysis revealed that M2-1 interacts specifically with positive sense RSV RNA, but not negative sense genome RNA. It also showed that M2-1 makes contacts along the length of each viral mRNA, indicating that M2-1 functions as a component of the transcriptase complex, transiently associating with nascent mRNA being extruded from the polymerase. In addition, we found that M2-1 binds specific cellular mRNAs. In contrast to the situation with RSV mRNA, M2-1 binds discrete sites within cellular mRNAs, with a preference for A/U rich sequences. These results suggest that in addition to its previously described role in transcription elongation, M2-1 might have an additional role involving cellular RNA interactions.

Approaches to Identify and Characterise the Post-Transcriptional Roles of lncRNAs in Cancer

It is becoming increasingly evident that the non-coding genome and transcriptome exert great influence over their coding counterparts through complex molecular interactions. Among non-coding RNAs (ncRNA), long non-coding RNAs (lncRNAs) in particular present increased potential to participate in dysregulation of post-transcriptional processes through both RNA and protein interactions. Since such processes can play key roles in contributing to cancer progression, it is desirable to continue expanding the search for lncRNAs impacting cancer through post-transcriptional mechanisms. The sheer diversity of mechanisms requires diverse resources and methods that have been developed and refined over the past decade. We provide an overview of computational resources as well as proven low-to-high throughput techniques to enable identification and characterisation of lncRNAs in their complex interactive contexts. As more cancer research strategies evolve to explore the non-coding genome and transcriptome, we anticipate this will provide a valuable primer and perspective of how these technologies have matured and will continue to evolve to assist researchers in elucidating post-transcriptional roles of lncRNAs in cancer.

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

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

The kinetic landscape of an RNA-binding protein in cells

Gene expression in higher eukaryotic cells orchestrates interactions between thousands of RNA-binding proteins (RBPs) and tens of thousands of RNAs¹. The kinetics by which RBPs bind to and dissociate from their RNA sites are critical for the coordination of cellular RNA-protein interactions². However, these kinetic parameters have not been experimentally measured in cells. Here we show that time-resolved RNA-protein cross-linking with a pulsed femtosecond ultraviolet laser, followed by immunoprecipitation and high-throughput sequencing, allows the determination of binding and dissociation kinetics of the RBP DAZL for thousands of individual RNA-binding sites in cells. This kinetic cross-linking and immunoprecipitation (KIN-CLIP) approach reveals that DAZL resides at individual binding sites for time periods of only seconds or shorter, whereas the binding sites remain DAZL-free for markedly longer. The data also indicate that DAZL binds to many RNAs in clusters of multiple proximal sites. The effect of DAZL on mRNA levels and ribosome association correlates with the cumulative probability of DAZL binding in these clusters. Integrating kinetic data with mRNA features quantitatively connects DAZL-RNA binding to DAZL function. Our results show how kinetic parameters for RNA-protein interactions can be measured in cells, and how these data link RBP-RNA binding to the cellular function of RBPs.

Emerging roles of N6-methyladenosine (m6A) modification in breast cancer

N6-Methyladenosine (m6A) is the most abundant, dynamic, and reversible epigenetic RNA modification that is found in coding and non-coding RNAs. Emerging studies have shown that m6A and its regulators affect multiple steps in RNA metabolism and play broad roles in various cancers. Worldwide, breast cancer is the most prevalent cancer in female. It is a very heterogeneous disease characterized by genetic and epigenetic variations in tumor cells. Increasing evidence has shown that the dysregulation of m6A-related effectors, as methyltransferases, demethylases, and m6A binding proteins, is pivotal in breast cancer pathogenesis. In this review, we have summarized the most up-to-date research on the biological functions of m6A modification in breast cancer and have discussed the potential clinical applications and future directions of m6A modification as a biomarker as well as a therapeutic target of breast cancer.

Clip for studying protein-RNA interactions that regulate virus replication

Viral and cellular RNA-binding proteins regulate numerous key steps in the replication of diverse virus genera. Viruses efficiently co-opt the host cell machinery for purposes such as transcription, splicing and subcellular localization of viral genomes. Though viral RNAs often need to resemble cellular RNAs to effectively utilize the cellular machinery, they still retain unique sequence and structural features for recognition by viral proteins for processes such as RNA polymerization, RNA export and selective packaging into virus particles. While beneficial for virus replication, distinct features of viral nucleic acids can also be recognized as foreign by several host defense proteins. Development of the crosslinking immunoprecipitation coupled with sequencing (CLIP) approach has allowed the study of viral and cellular RNA binding proteins that regulate critical aspects of viral replication in unprecedented detail. By combining immunoprecipitation of covalently crosslinked protein-RNA complexes with high-throughput sequencing, CLIP provides a global account of RNA sequences bound by RNA-binding proteins of interest in physiological settings and at near-nucleotide resolution. Here, we describe the step-by-step application of the CLIP methodology within the context of two cellular splicing regulatory proteins, hnRNP A1 and hnRNP H1 that regulate HIV-1 splicing. In principle, this versatile protocol can be applied to many other viral and cellular RNA-binding proteins.

Host-derived circular RNAs display proviral activities in Hepatitis C virus-infected cells

Viruses subvert macromolecular pathways in infected host cells to aid in viral gene amplification or to counteract innate immune responses. Roles for host-encoded, noncoding RNAs, including microRNAs, have been found to provide pro- and anti-viral functions. Recently, circular RNAs (circRNAs), that are generated by a nuclear back-splicing mechanism of pre-mRNAs, have been implicated to have roles in DNA virus-infected cells. This study examines the circular RNA landscape in uninfected and hepatitis C virus (HCV)-infected liver cells. Results showed that the abundances of distinct classes of circRNAs were up-regulated or down-regulated in infected cells. Identified circRNAs displayed pro-viral effects. One particular up-regulated circRNA, circPSD3, displayed a very pronounced effect on viral RNA abundances in both hepatitis C virus- and Dengue virus-infected cells. Though circPSD3 has been shown to bind factor eIF4A3 that modulates the cellular nonsense-mediated decay (NMD) pathway, circPSD3 regulates RNA amplification in a pro-viral manner at a post-translational step, while eIF4A3 exhibits the anti-viral property of the NMD pathway. Findings from the global analyses of the circular RNA landscape argue that pro-, and likely, anti-viral functions are executed by circRNAs that modulate viral gene expression as well as host pathways. Because of their long half-lives, circRNAs likely play hitherto unknown, important roles in viral pathogenesis.

A Day in the Life of the Exon Junction Complex

The exon junction complex (EJC) is an abundant messenger ribonucleoprotein (mRNP) component that is assembled during splicing and binds to mRNAs upstream of exon-exon junctions. EJCs accompany the mRNA during its entire life in the nucleus and the cytoplasm and communicate the information about the splicing process and the position of introns. Specifically, the EJC’s core components and its associated proteins regulate different steps of gene expression, including pre-mRNA splicing, mRNA export, translation, and nonsense-mediated mRNA decay (NMD). This review summarizes the most important functions and main protagonists in the life of the EJC. It also provides an overview of the latest findings on the assembly, composition and molecular activities of the EJC and presents them in the chronological order, in which they play a role in the EJC’s life cycle.

How Do You Identify m6 A Methylation in Transcriptomes at High Resolution? A Comparison of Recent Datasets

A flurry of methods has been developed in recent years to identify N6-methyladenosine (m6A) sites across transcriptomes at high resolution. This raises the need to understand both the common features and those that are unique to each method. Here, we complement the analyses presented in the original papers by reviewing their various technical aspects and comparing the overlap between m6A-methylated messenger RNAs (mRNAs) identified by each. Specifically, we examine eight different methods that identify m6A sites in human cells with high resolution: two antibody-based crosslinking and immunoprecipitation (CLIP) approaches, two using endoribonuclease MazF, one based on deamination, two using Nanopore direct RNA sequencing, and finally, one based on computational predictions. We contrast the respective datasets and discuss the challenges in interpreting the overlap between them, including a prominent expression bias in detected genes. This overview will help guide researchers in making informed choices about using the available data and assist with the design of future experiments to expand our understanding of m6A and its regulation.

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.

A systems view of spliceosomal assembly and branchpoints with iCLIP

Studies of spliceosomal interactions are challenging due to their dynamic nature. Here we used spliceosome iCLIP, which immunoprecipitates SmB along with small nuclear ribonucleoprotein particles and auxiliary RNA binding proteins, to map spliceosome engagement with pre-messenger RNAs in human cell lines. This revealed seven peaks of spliceosomal crosslinking around branchpoints (BPs) and splice sites. We identified RNA binding proteins that crosslink to each peak, including known and candidate splicing factors. Moreover, we detected the use of over 40,000 BPs with strong sequence consensus and structural accessibility, which align well to nearby crosslinking peaks. We show how the position and strength of BPs affect the crosslinking patterns of spliceosomal factors, which bind more efficiently upstream of strong or proximally located BPs and downstream of weak or distally located BPs. These insights exemplify spliceosome iCLIP as a broadly applicable method for transcriptomic studies of splicing mechanisms.

CLIP and RNA interactome studies to unravel genome-wide RNA-protein interactions in vivo in Arabidopsis thaliana

Post-transcriptional regulation makes an important contribution to adjusting the transcriptome to environmental changes in plants. RNA-binding proteins are key players that interact specifically with mRNAs to co-ordinate their fate. While the regulatory interactions between proteins and RNA are well understood in animals, until recently little information was available on the global binding landscape of RNA-binding proteins in higher plants. This is not least due to technical challenges in plants. In turn, while numerous RNA-binding proteins have been identified through mutant analysis and homology-based searches in plants, only recently a full compendium of proteins with RNA-binding activity has been experimentally determined for the reference plant Arabidopsis thaliana. State-of-the-art techniques to determine RNA-protein interactions genome-wide in animals are based on the covalent fixation of RNA and protein in vivo by UV light. This has only recently been successfully applied to plants. Here, we present practical considerations on the application of UV irradiation based methods to comprehensively determine in vivo RNA-protein interactions in Arabidopsis thaliana, focussing on individual nucleotide resolution crosslinking immunoprecipitation (iCLIP) and mRNA interactome capture.

The landscape of the non-canonical RNA-binding site of Gemin5 unveils a feedback loop counteracting the negative effect on translation

Gemin5 is a predominantly cytoplasmic protein that downregulates translation, beyond controlling snRNPs assembly. The C-terminal region harbors a non-canonical RNA-binding site consisting of two domains, RBS1 and RBS2, which differ in RNA-binding capacity and the ability to modulate translation. Here, we show that these domains recognize distinct RNA targets in living cells. Interestingly, the most abundant and exclusive RNA target of the RBS1 domain was Gemin5 mRNA. Biochemical and functional characterization of this target demonstrated that RBS1 polypeptide physically interacts with a predicted thermodynamically stable stem–loop upregulating mRNA translation, thereby counteracting the negative effect of Gemin5 protein on global protein synthesis. In support of this result, destabilization of the stem–loop impairs the stimulatory effect on translation. Moreover, RBS1 stimulates translation of the endogenous Gemin5 mRNA. Hence, although the RBS1 domain downregulates global translation, it positively enhances translation of RNA targets carrying thermodynamically stable secondary structure motifs. This mechanism allows fine-tuning the availability of Gemin5 to play its multiple roles in gene expression control.

RNA 5-Methylcytosine Facilitates the Maternal-to-Zygotic Transition by Preventing Maternal mRNA Decay

The maternal-to-zygotic transition (MZT) is a conserved and fundamental process during which the maternal environment is converted to an environment of embryonic-driven development through dramatic reprogramming. However, how maternally supplied transcripts are dynamically regulated during MZT remains largely unknown. Herein, through genome-wide profiling of RNA 5-methylcytosine (m⁵C) modification in zebrafish early embryos, we found that m⁵C-modified maternal mRNAs display higher stability than non-m⁵C-modified mRNAs during MZT. We discovered that Y-box binding protein 1 (Ybx1) preferentially recognizes m⁵C-modified mRNAs through π-π interactions with a key residue, Trp45, in Ybx1's cold shock domain (CSD), which plays essential roles in maternal mRNA stability and early embryogenesis of zebrafish. Together with the mRNA stabilizer Pabpc1a, Ybx1 promotes the stability of its target mRNAs in an m⁵C-dependent manner. Our study demonstrates an unexpected mechanism of RNA m⁵C-regulated maternal mRNA stabilization during zebrafish MZT, highlighting the critical role of m⁵C mRNA modification in early development.

Monitored eCLIP: high accuracy mapping of RNA-protein interactions

CLIP-seq methods provide transcriptome-wide snapshots of RNA-protein interactions in live cells. Reverse transcriptases stopping at cross-linked nucleotides sign for RNA-protein binding sites. Reading through cross-linked positions results in false binding site assignments. In the ‘monitored enhanced CLIP’ (meCLIP) method, a barcoded biotinylated linker is ligated at the 5′ end of cross-linked RNA fragments to purify RNA prior to the reverse transcription. cDNAs keeping the barcode sequence correspond to reverse transcription read-throughs. Read through occurs in unpredictable proportions, representing up to one fourth of total reads. Filtering out those reads strongly improves reliability and precision in protein binding site assignment.

CLIPick: a sensitive peak caller for expression-based deconvolution of HITS-CLIP signals

High-throughput sequencing of RNAs isolated by crosslinking immunoprecipitation (HITS-CLIP, also called CLIP-Seq) has been used to map global RNA–protein interactions. However, a critical caveat of HITS-CLIP results is that they contain non-linear background noise—different extent of non-specific interactions caused by individual transcript abundance—that has been inconsiderately normalized, resulting in sacrifice of sensitivity. To properly deconvolute RNA–protein interactions, we have implemented CLIPick, a flexible peak calling pipeline for analyzing HITS-CLIP data, which statistically determines the signal-to-noise ratio for each transcript based on the expression-dependent background simulation. Comprising of streamlined Python modules with an easy-to-use standalone graphical user interface, CLIPick robustly identifies significant peaks and quantitatively defines footprint regions within which RNA–protein interactions were occurred. CLIPick outperforms other peak callers in accuracy and sensitivity, selecting the largest number of peaks particularly in lowly expressed transcripts where such marginal signals are hard to discriminate. Specifically, the application of CLIPick to Argonaute (Ago) HITS-CLIP data were sensitive enough to uncover extended features of microRNA target sites, and these sites were experimentally validated. CLIPick enables to resolve critical interactions in a wide spectrum of transcript levels and extends the scope of HITS-CLIP analysis. CLIPick is available at: http://clip.korea.ac.kr/clipick/

SpyCLIP: an easy-to-use and high-throughput compatible CLIP platform for the characterization of protein-RNA interactions with high accuracy

UV crosslinking and immunoprecipitation (CLIP) coupled with high-throughput sequencing is the most state-of-the-art technology to characterize protein-RNA interactions, yet only a small portion of RNA-binding proteins (RBPs) have been studied by CLIP due to its complex procedures and high level of false-positive signals. Herein, we report a SpyCLIP method that employs a covalent linkage formed between the RBP-fused SpyTag and SpyCatcher, which can withstand the harshest washing conditions for removing nonspecific interactions. Moreover, SpyCLIP circumvents the radioactive labeling and PAGE-membrane purification steps, and the whole procedure can be performed on beads and is readily amenable to automation. We investigated multiple RBPs by SpyCLIP and generated high-quality RNA binding maps with significantly improved reproductivity and accuracy. Therefore, the small tag size and convenient protocol of SpyCLIP provides a robust method for both routine characterization and high-throughput studies of protein-RNA interactions.

Splicing of long non-coding RNAs primarily depends on polypyrimidine tract and 5' splice-site sequences due to weak interactions with SR proteins

Many nascent long non-coding RNAs (lncRNAs) undergo the same maturation steps as pre-mRNAs of protein-coding genes (PCGs), but they are often poorly spliced. To identify the underlying mechanisms for this phenomenon, we searched for putative splicing inhibitory sequences using the ncRNA-a2 as a model. Genome-wide analyses of intergenic lncRNAs (lincRNAs) revealed that lincRNA splicing efficiency positively correlates with 5'ss strength while no such correlation was identified for PCGs. In addition, efficiently spliced lincRNAs have higher thymidine content in the polypyrimidine tract (PPT) compared to efficiently spliced PCGs. Using model lincRNAs, we provide experimental evidence that strengthening the 5'ss and increasing the T content in PPT significantly enhances lincRNA splicing. We further showed that lincRNA exons contain less putative binding sites for SR proteins. To map binding of SR proteins to lincRNAs, we performed iCLIP with SRSF2, SRSF5 and SRSF6 and analyzed eCLIP data for SRSF1, SRSF7 and SRSF9. All examined SR proteins bind lincRNA exons to a much lower extent than expression-matched PCGs. We propose that lincRNAs lack the cooperative interaction network that enhances splicing, which renders their splicing outcome more dependent on the optimality of splice sites.

The Human RNA-Binding Proteome and Its Dynamics during Translational Arrest

Proteins and RNA functionally and physically intersect in multiple biological processes, however, currently no universal method is available to purify protein-RNA complexes. Here, we introduce XRNAX, a method for the generic purification of protein-crosslinked RNA, and demonstrate its versatility to study the composition and dynamics of protein-RNA interactions by various transcriptomic and proteomic approaches. We show that XRNAX captures all RNA biotypes and use this to characterize the sub-proteomes that interact with coding and non-coding RNAs (ncRNAs) and to identify hundreds of protein-RNA interfaces. Exploiting the quantitative nature of XRNAX, we observe drastic remodeling of the RNA-bound proteome during arsenite-induced stress, distinct from autophagy-related changes in the total proteome. In addition, we combine XRNAX with crosslinking immunoprecipitation sequencing (CLIP-seq) to validate the interaction of ncRNA with lamin B1 and EXOSC2. Thus, XRNAX is a resourceful approach to study structural and compositional aspects of protein-RNA interactions to address fundamental questions in RNA-biology.

Exon Junction Complex Shapes the Transcriptome by Repressing Recursive Splicing

Recursive splicing (RS) starts by defining an "RS-exon," which is then spliced to the preceding exon, thus creating a recursive 5' splice site (RS-5ss). Previous studies focused on cryptic RS-exons, and now we find that the exon junction complex (EJC) represses RS of hundreds of annotated, mainly constitutive RS-exons. The core EJC factors, and the peripheral factors PNN and RNPS1, maintain RS-exon inclusion by repressing spliceosomal assembly on RS-5ss. The EJC also blocks 5ss located near exon-exon junctions, thus repressing inclusion of cryptic microexons. The prevalence of annotated RS-exons is high in deuterostomes, while the cryptic RS-exons are more prevalent in Drosophila, where EJC appears less capable of repressing RS. Notably, incomplete repression of RS also contributes to physiological alternative splicing of several human RS-exons. Finally, haploinsufficiency of the EJC factor Magoh in mice is associated with skipping of RS-exons in the brain, with relevance to the microcephaly phenotype and human diseases.

A Short Tandem Repeat-Enriched RNA Assembles a Nuclear Compartment to Control Alternative Splicing and Promote Cell Survival

Functions of many long noncoding RNAs (lncRNAs) depend on their ability to interact with multiple copies of specific RNA-binding proteins (RBPs). Here, we devised a workflow combining bioinformatics and experimental validation steps to systematically identify RNAs capable of multivalent RBP recruitment. This uncovered a number of previously unknown transcripts encoding high-density RBP recognition arrays within genetically normal short tandem repeats. We show that a top-scoring hit in this screen, lncRNA PNCTR, contains hundreds of pyrimidine tract-binding protein (PTBP1)-specific motifs allowing it to sequester a substantial fraction of PTBP1 in a nuclear body called perinucleolar compartment. Importantly, PNCTR is markedly overexpressed in a variety of cancer cells and its downregulation is sufficient to induce programmed cell death at least in part by stimulating PTBP1 splicing regulation activity. This work expands our understanding of the repeat-containing fraction of the human genome and illuminates a novel mechanism driving malignant transformation of cancer cells.

•

Human genome encodes many transcripts enriched in short tandem repeats (strRNAs)

•

strRNA PNCTR recruits RNA-binding protein PTBP1 to a nuclear body called PNC

•

PNCTR antagonizes splicing regulation function of PTBP1 and promotes cell survival

•

PNCTR is dramatically upregulated in a wide range of cancer cells

RNA-protein interactions in an unstructured context

Despite their importance, our understanding of noncovalent RNA-protein interactions is incomplete. This especially concerns the binding between RNA and unstructured protein regions, a widespread class of such interactions. Here, we review the recent experimental and computational work on RNA-protein interactions in an unstructured context with a particular focus on how such interactions may be shaped by the intrinsic interaction affinities between individual nucleobases and protein side chains. Specifically, we articulate the claim that the universal genetic code reflects the binding specificity between nucleobases and protein side chains and that, in turn, the code may be seen as the Rosetta stone for understanding RNA-protein interactions in general.

Heteromeric RNP Assembly at LINEs Controls Lineage-Specific RNA Processing

Long mammalian introns make it challenging for the RNA processing machinery to identify exons accurately. We find that LINE-derived sequences (LINEs) contribute to this selection by recruiting dozens of RNA-binding proteins (RBPs) to introns. This includes MATR3, which promotes binding of PTBP1 to multivalent binding sites within LINEs. Both RBPs repress splicing and 3' end processing within and around LINEs. Notably, repressive RBPs preferentially bind to evolutionarily young LINEs, which are located far from exons. These RBPs insulate the LINEs and the surrounding intronic regions from RNA processing. Upon evolutionary divergence, changes in RNA motifs within LINEs lead to gradual loss of their insulation. Hence, older LINEs are located closer to exons, are a common source of tissue-specific exons, and increasingly bind to RBPs that enhance RNA processing. Thus, LINEs are hubs for the assembly of repressive RBPs and also contribute to the evolution of new, lineage-specific transcripts in mammals. VIDEO ABSTRACT.

The Future of Cross-Linking and Immunoprecipitation (CLIP)

To understand the assembly and functional outcomes of protein-RNA regulation, it is crucial to precisely identify the positions of such interactions. Cross-linking and immunoprecipitation (CLIP) serves this purpose by exploiting covalent protein-RNA cross-linking and RNA fragmentation, along with a series of stringent purification and quality control steps to prepare complementary DNA (cDNA) libraries for sequencing. Here we describe the core steps of CLIP, its primary variations, and the approaches to data analysis. We present the application of CLIP to studies of specific cell types in genetically engineered mice and discuss the mechanistic and physiologic insights that have already been gained from studies using CLIP. We conclude by discussing the future opportunities for CLIP, including studies of human postmortem tissues from disease patients and controls, RNA epigenetic modifications, and RNA structure. These and other applications of CLIP will continue to unravel fundamental gene regulatory mechanisms while providing important biologic and clinically relevant insights.

Data Science Issues in Studying Protein–RNA Interactions with CLIP Technologies

An interplay of experimental and computational methods is required to achieve a comprehensive understanding of protein–RNA interactions. UV crosslinking and immunoprecipitation (CLIP) identifies endogenous interactions by sequencing RNA fragments that copurify with a selected RNA-binding protein under stringent conditions. Here we focus on approaches for the analysis of the resulting data and appraise the methods for peak calling, visualization, analysis, and computational modeling of protein–RNA binding sites. We advocate that the sensitivity and specificity of data be assessed in combination for computational quality control. Moreover, we demonstrate the value of analyzing sequence motif enrichment in peaks assigned from CLIP data and of visualizing RNA maps, which examine the positional distribution of peaks around regulated landmarks in transcripts. We use these to assess how variations in CLIP data quality and in different peak calling methods affect the insights into regulatory mechanisms. We conclude by discussing future opportunities for the computational analysis of protein–RNA interaction experiments.

Individual Nucleotide Resolution UV Cross-Linking and Immunoprecipitation (iCLIP) to Determine Protein-RNA Interactions

RNA-binding proteins (RBPs) interact with and determine the fate of many cellular RNA transcripts. In doing so they help direct many essential roles in cellular physiology, while their perturbed activity can contribute to disease etiology. In this chapter we detail a functional genomics approach, termed individual nucleotide resolution UV cross-linking and immunoprecipitation (iCLIP), that can determine the interactions of RBPs with their RNA targets in high throughput and at nucleotide resolution. iCLIP achieves this by exploiting UV-induced covalent cross-links formed between RBPs and their target RNAs to both purify the RBP-RNA complexes under stringent conditions, and to cause reverse transcription stalling that then identifies the direct cross-link sites in the high throughput sequenced cDNA libraries.

'Read-through marking' reveals differential nucleotide composition of read-through and truncated cDNAs in iCLIP

We established a modified iCLIP protocol, called 'read-through marking', which facilitates the detection of cDNAs that have not been truncated upon encountering the RNA-peptide complex during reverse transcription (read-through cDNAs). A large proportion of these cDNAs would be undesirable in an iCLIP library, as it could affect the resolution of the method. To this end, we added an oligonucleotide to the 5'-end of RNA fragments-a 5'-marker-to mark the read-through cDNAs. By applying this modified iCLIP protocol to PTBP1 and eIF4A3, we found that the start sites of read-through cDNAs are enriched in adenosines, while the remaining cDNAs have a markedly different sequence content at their starts, preferentially containing thymidines. This finding in turn indicates that most of the reads in our iCLIP libraries are a product of truncation with valuable information regarding the proteins' RNA-binding sites. Thus, cDNA start sites confidently identify a protein's RNA-crosslink sites and we can account for the impact of read-through cDNAs by commonly adding a 5'-marker.

CLIP-related methodologies and their application to retrovirology

Virtually every step of HIV-1 replication and numerous cellular antiviral defense mechanisms are regulated by the binding of a viral or cellular RNA-binding protein (RBP) to distinct sequence or structural elements on HIV-1 RNAs. Until recently, these protein-RNA interactions were studied largely by in vitro binding assays complemented with genetics approaches. However, these methods are highly limited in the identification of the relevant targets of RBPs in physiologically relevant settings. Development of crosslinking-immunoprecipitation sequencing (CLIP) methodology has revolutionized the analysis of protein-nucleic acid complexes. CLIP combines immunoprecipitation of covalently crosslinked protein-RNA complexes with high-throughput sequencing, providing a global account of RNA sequences bound by a RBP of interest in cells (or virions) at near-nucleotide resolution. Numerous variants of the CLIP protocol have recently been developed, some with major improvements over the original. Herein, we briefly review these methodologies and give examples of how CLIP has been successfully applied to retrovirology research.

GoldCLIP: Gel-omitted Ligation-dependent CLIP

Protein-RNA interaction networks are essential to understand gene regulation control. Identifying binding sites of RNA-binding proteins (RBPs) by the UV-crosslinking and immunoprecipitation (CLIP) represents one of the most powerful methods to map protein-RNA interactions in vivo. However, the traditional CLIP protocol is technically challenging, which requires radioactive labeling and suffers from material loss during PAGE-membrane transfer procedures. Here we introduce a super-efficient CLIP method (GoldCLIP) that omits all gel purification steps. This nonisotopic method allows us to perform highly reproducible CLIP experiments with polypyrimidine tract-binding protein (PTB), a classical RBP in human cell lines. In principle, our method guarantees sequencing library constructions, providing the protein of interest can be successfully crosslinked to RNAs in living cells. GoldCLIP is readily applicable to diverse proteins to uncover their endogenous RNA targets.

uvCLAP is a fast and non-radioactive method to identify in vivo targets of RNA-binding proteins

RNA-binding proteins (RBPs) play important and essential roles in eukaryotic gene expression regulating splicing, localization, translation, and stability of mRNAs. We describe ultraviolet crosslinking and affinity purification (uvCLAP), an easy-to-use, robust, reproducible, and high-throughput method to determine in vivo targets of RBPs. uvCLAP is fast and does not rely on radioactive labeling of RNA. We investigate binding of 15 RBPs from fly, mouse, and human cells to test the method's performance and applicability. Multiplexing of signal and control libraries enables straightforward comparison of samples. Experiments for most proteins achieve high enrichment of signal over background. A point mutation and a natural splice isoform that change the RBP subcellular localization dramatically alter target selection without changing the targeted RNA motif, showing that compartmentalization of RBPs can be used as an elegant means to generate RNA target specificity.

A brave new world of RNA-binding proteins

RNA-binding proteins (RBPs) are typically thought of as proteins that bind RNA through one or multiple globular RNA-binding domains (RBDs) and change the fate or function of the bound RNAs. Several hundred such RBPs have been discovered and investigated over the years. Recent proteome-wide studies have more than doubled the number of proteins implicated in RNA binding and uncovered hundreds of additional RBPs lacking conventional RBDs. In this Review, we discuss these new RBPs and the emerging understanding of their unexpected modes of RNA binding, which can be mediated by intrinsically disordered regions, protein-protein interaction interfaces and enzymatic cores, among others. We also discuss the RNA targets and molecular and cellular functions of the new RBPs, as well as the possibility that some RBPs may be regulated by RNA rather than regulate RNA.

PureCLIP: capturing target-specific protein-RNA interaction footprints from single-nucleotide CLIP-seq data

The iCLIP and eCLIP techniques facilitate the detection of protein–RNA interaction sites at high resolution, based on diagnostic events at crosslink sites. However, previous methods do not explicitly model the specifics of iCLIP and eCLIP truncation patterns and possible biases. We developed PureCLIP (https://github.com/skrakau/PureCLIP), a hidden Markov model based approach, which simultaneously performs peak-calling and individual crosslink site detection. It explicitly incorporates a non-specific background signal and, for the first time, non-specific sequence biases. On both simulated and real data, PureCLIP is more accurate in calling crosslink sites than other state-of-the-art methods and has a higher agreement across replicates.

Adaptation of iCLIP to plants determines the binding landscape of the clock-regulated RNA-binding protein AtGRP7

Background

Functions for RNA-binding proteins in orchestrating plant development and environmental responses are well established. However, the lack of a genome-wide view of their in vivo binding targets and binding landscapes represents a gap in understanding the mode of action of plant RNA-binding proteins. Here, we adapt individual nucleotide resolution crosslinking and immunoprecipitation (iCLIP) genome-wide to determine the binding repertoire of the circadian clock-regulated Arabidopsis thaliana glycine-rich RNA-binding protein AtGRP7.

Results

iCLIP identifies 858 transcripts with significantly enriched crosslink sites in plants expressing AtGRP7-GFP that are absent in plants expressing an RNA-binding-dead AtGRP7 variant or GFP alone. To independently validate the targets, we performed RNA immunoprecipitation (RIP)-sequencing of AtGRP7-GFP plants subjected to formaldehyde fixation. Of the iCLIP targets, 452 were also identified by RIP-seq and represent a set of high-confidence binders. AtGRP7 can bind to all transcript regions, with a preference for 3′ untranslated regions. In the vicinity of crosslink sites, U/C-rich motifs are overrepresented. Cross-referencing the targets against transcriptome changes in AtGRP7 loss-of-function mutants or AtGRP7-overexpressing plants reveals a predominantly negative effect of AtGRP7 on its targets. In particular, elevated AtGRP7 levels lead to damping of circadian oscillations of transcripts, including DORMANCY/AUXIN ASSOCIATED FAMILY PROTEIN2 and CCR-LIKE. Furthermore, several targets show changes in alternative splicing or polyadenylation in response to altered AtGRP7 levels.

Conclusions

We have established iCLIP for plants to identify target transcripts of the RNA-binding protein AtGRP7. This paves the way to investigate the dynamics of posttranscriptional networks in response to exogenous and endogenous cues.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1332-x) contains supplementary material, which is available to authorized users.

CLIP-seq analysis of multi-mapped reads discovers novel functional RNA regulatory sites in the human transcriptome

Crosslinking or RNA immunoprecipitation followed by sequencing (CLIP-seq or RIP-seq) allows transcriptome-wide discovery of RNA regulatory sites. As CLIP-seq/RIP-seq reads are short, existing computational tools focus on uniquely mapped reads, while reads mapped to multiple loci are discarded. We present CLAM (CLIP-seq Analysis of Multi-mapped reads). CLAM uses an expectation–maximization algorithm to assign multi-mapped reads and calls peaks combining uniquely and multi-mapped reads. To demonstrate the utility of CLAM, we applied it to a wide range of public CLIP-seq/RIP-seq datasets involving numerous splicing factors, microRNAs and m⁶A RNA methylation. CLAM recovered a large number of novel RNA regulatory sites inaccessible by uniquely mapped reads. The functional significance of these sites was demonstrated by consensus motif patterns and association with alternative splicing (splicing factors), transcript abundance (AGO2) and mRNA half-life (m⁶A). CLAM provides a useful tool to discover novel protein–RNA interactions and RNA modification sites from CLIP-seq and RIP-seq data, and reveals the significant contribution of repetitive elements to the RNA regulatory landscape of the human transcriptome.