PIPE-CLIP: a comprehensive online tool for CLIP-seq data analysis

Plant immunity suppressor SKRP encodes a novel RNA-binding protein that targets exon 3' end of unspliced RNA

The regulatory roles of RNA splicing in plant immunity are emerging but still largely obscure. We reported previously that Phytophthora pathogen effector Avr3c targets a soybean protein SKRP (serine/lysine/arginine-rich protein) to impair soybean basal immunity by regulating host pre-mRNA alternative splicing, while the biochemical nature of SKRP remains unknown. Here, by using Arabidopsis as a model, we studied the mechanism of SKRP in regulating pre-mRNA splicing and plant immunity. AtSKRP confers impaired plant immunity against Phytophthora capsici and associates with spliceosome component PRP8 and splicing factor SR45, which positively and negatively regulate plant immunity, respectively. Enhanced crosslinking and immunoprecipitation followed by high-throughput sequencing (eCLIP-seq) showed AtSKRP is a novel RNA-binding protein that targets exon 3' end of unspliced RNA. Such position-specific binding of SKRP is associated with its activity in suppressing intron retention, including at positive immune regulatory genes UBP25 and RAR1. In addition, we found AtSKRP self-interact and forms oligomer, and these properties are associated with its function in plant immunity. Overall, our findings reveal that the immune repressor SKRP is a spliceosome-associated protein that targets exon 3' end to regulate pre-mRNA splicing in Arabidopsis.

htseq-clip: a toolset for the preprocessing of eCLIP/iCLIP datasets

SUMMARY

Transcriptome-wide detection of binding sites of RNA-binding proteins is achieved using Individual-nucleotide crosslinking and immunoprecipitation (iCLIP) and its derivative enhanced CLIP (eCLIP) sequencing methods. Here, we introduce htseq-clip, a python package developed for preprocessing, extracting and summarizing crosslink site counts from i/eCLIP experimental data. The package delivers crosslink site count matrices along with other metrics, which can be directly used for filtering and downstream analyses such as the identification of differential binding sites.

AVAILABILITY AND IMPLEMENTATION

The Python package htseq-clip is available via pypi (python package index), bioconda and the Galaxy Tool Shed under the open source MIT License. The code is hosted at https://github.com/EMBL-Hentze-group/htseq-clip and documentation is available under https://htseq-clip.readthedocs.io/en/latest.

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.

Forming cytoplasmic stress granules PURα suppresses mRNA translation initiation of IGFBP3 to promote esophageal squamous cell carcinoma progression

Esophageal squamous cell carcinoma (ESCC) is one of the most fatal malignancies worldwide. Recently, our group identified purine-rich element binding protein alpha (PURα), a single-stranded DNA/RNA-binding protein, to be significantly associated with the progression of ESCC. Additional immunofluorescence staining demonstrated that PURα forms cytoplasmic stress granules to suppress mRNA translation initiation. The expression level of cytoplasmic PURα in ESCC tumor tissues was significantly higher than that in adjacent epithelia and correlated with a worse patient survival rate by immunohistochemistry. Functionally, PURα strongly preferred to bind to UG-/U-rich motifs and mRNA 3´UTR by CLIP-seq analysis. Moreover, PURα knockout significantly increased the protein level of insulin-like growth factor binding protein 3 (IGFBP3). In addition, it was further demonstrated that PURα-interacting proteins are remarkably associated with translation initiation factors and ribosome-related proteins and that PURα regulates protein expression by interacting with translation initiation factors, such as PABPC1, eIF3B and eIF3F, in an RNA-independent manner, while the interaction with ribosome-related proteins is significantly dependent on RNA. Specifically, PURα was shown to interact with the mRNA 3´UTR of IGFBP3 and inhibit its expression by suppressing mRNA translation initiation. Together, this study identifies cytoplasmic PURα as a modulator of IGFBP3, which could be a promising therapeutic target for ESCC treatment.

Multiple Viral microRNAs Regulate Interferon Release and Signaling Early during Infection with Epstein-Barr Virus

Epstein-Barr virus (EBV), a human herpesvirus, encodes 44 microRNAs (miRNAs), which regulate many genes with various functions in EBV-infected cells. Multiple target genes of the EBV miRNAs have been identified, some of which play important roles in adaptive antiviral immune responses. Using EBV mutant derivatives, we identified additional roles of viral miRNAs in governing versatile type I interferon (IFN) responses upon infection of human primary mature B cells. We also found that Epstein-Barr virus-encoded small RNAs (EBERs) and LF2, viral genes with previously reported functions in inducing or regulating IFN-I pathways, had negligible or even contrary effects on secreted IFN-α in our model. Data mining and Ago PAR-CLIP experiments uncovered more than a dozen previously uncharacterized, direct cellular targets of EBV miRNA associated with type I IFN pathways. We also identified indirect targets of EBV miRNAs in B cells, such as TRL7 and TLR9, in the prelatent phase of infection. The presence of epigenetically naive, non-CpG methylated viral DNA was essential to induce IFN-α secretion during EBV infection in a TLR9-dependent manner. In a newly established fusion assay, we verified that EBV virions enter a subset of plasmacytoid dendritic cells (pDCs) and determined that these infected pDCs are the primary producers of IFN-α in EBV-infected peripheral blood mononuclear cells. Our findings document that many EBV-encoded miRNAs regulate type I IFN response in newly EBV infected primary human B cells in the prelatent phase of infection and dampen the acute release of IFN-α in pDCs upon their encounter with EBV.IMPORTANCE Acute antiviral functions of all nucleated cells rely on type I interferon (IFN-I) pathways triggered upon viral infection. Host responses encompass the sensing of incoming viruses, the activation of specific transcription factors that induce the transcription of IFN-I genes, the secretion of different IFN-I types and their recognition by the heterodimeric IFN-α/β receptor, the subsequent activation of JAK/STAT signaling pathways, and, finally, the transcription of many IFN-stimulated genes (ISGs). In sum, these cellular functions establish a so-called antiviral state in infected and neighboring cells. To counteract these cellular defense mechanisms, viruses have evolved diverse strategies and encode gene products that target antiviral responses. Among such immune-evasive factors are viral microRNAs (miRNAs) that can interfere with host gene expression. We discovered that multiple miRNAs of Epstein-Barr virus (EBV) control over a dozen cellular genes that contribute to the antiviral states of immune cells, specifically B cells and plasmacytoid dendritic cells (pDCs). We identified the viral DNA genome as the activator of IFN-α and question the role of abundant EBV EBERs, that, contrary to previous reports, do not have an apparent inducing function in the IFN-I pathway early after infection.

Techniques for Characterizing Cytomegalovirus-Encoded miRNAs

microRNAs (miRNAs) are small noncoding RNAs that regulate gene expression at the posttranscriptional level by binding to sites within the 3' untranslated regions of messenger RNA (mRNA) transcripts. The discovery of this completely new mechanism of gene regulation necessitated the development of a variety of techniques to further characterize miRNAs, their expression, and function. In this chapter, we will discuss techniques currently used in the miRNA field to detect, express and inhibit miRNAs, as well as methods used to identify and validate their targets, specifically with respect to the miRNAs encoded by human cytomegalovirus.

Galaxy CLIP-Explorer: a web server for CLIP-Seq data analysis

BACKGROUND

Post-transcriptional regulation via RNA-binding proteins plays a fundamental role in every organism, but the regulatory mechanisms lack important understanding. Nevertheless, they can be elucidated by cross-linking immunoprecipitation in combination with high-throughput sequencing (CLIP-Seq). CLIP-Seq answers questions about the functional role of an RNA-binding protein and its targets by determining binding sites on a nucleotide level and associated sequence and structural binding patterns. In recent years the amount of CLIP-Seq data skyrocketed, urging the need for an automatic data analysis that can deal with different experimental set-ups. However, noncanonical data, new protocols, and a huge variety of tools, especially for peak calling, made it difficult to define a standard.

FINDINGS

CLIP-Explorer is a flexible and reproducible data analysis pipeline for iCLIP data that supports for the first time eCLIP, FLASH, and uvCLAP data. Individual steps like peak calling can be changed to adapt to different experimental settings. We validate CLIP-Explorer on eCLIP data, finding similar or nearly identical motifs for various proteins in comparison with other databases. In addition, we detect new sequence motifs for PTBP1 and U2AF2. Finally, we optimize the peak calling with 3 different peak callers on RBFOX2 data, discuss the difficulty of the peak-calling step, and give advice for different experimental set-ups.

CONCLUSION

CLIP-Explorer finally fills the demand for a flexible CLIP-Seq data analysis pipeline that is applicable to the up-to-date CLIP protocols. The article further shows the limitations of current peak-calling algorithms and the importance of a robust peak detection.

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.

rMAPS2: an update of the RNA map analysis and plotting server for alternative splicing regulation

The rMAPS2 (RNA Map Analysis and Plotting Server 2) web server, freely available at http://rmaps.cecsresearch.org/, has provided the high-throughput sequencing data research community with curated tools for the identification of RNA binding protein sites. rMAPS2 analyzes differential alternative splicing or CLIP peak data obtained from high-throughput sequencing data analysis tools like MISO, rMATS, Piranha, PIPE-CLIP and PARalyzer, and then, graphically displays enriched RNA-binding protein target sites. The initial release of rMAPS focused only on the most common alternative splicing event, skipped exon or exon skipping. However, there was a high demand for the analysis of other major types of alternative splicing events, especially for retained intron events since this is the most common type of alternative splicing in plants, such as Arabidopsis thaliana. Here, we expanded the implementation of rMAPS2 to facilitate analyses for all five major types of alternative splicing events: skipped exon, mutually exclusive exons, alternative 5′ splice site, alternative 3′ splice site and retained intron. In addition, by employing multi-threading, rMAPS2 has vastly improved the user experience with significant reductions in running time, ∼3.5 min for the analysis of all five major alternative splicing types at once.

Loss of Nuclear TDP-43 Is Associated with Decondensation of LINE Retrotransposons

Loss of the nuclear RNA binding protein TAR DNA binding protein-43 (TDP-43) into cytoplasmic aggregates is the strongest correlate to neurodegeneration in amyotrophic lateral sclerosis and frontotemporal degeneration. The molecular changes associated with the loss of nuclear TDP-43 in human tissues are not entirely known. Using subcellular fractionation andfluorescent-activated cell sorting to enrich for diseased neuronal nuclei without TDP-43 from post-mortem frontotemporal degeneration-amyotro-phic lateral sclerosis (FTD-ALS) human brain, we characterized the effects of TDP-43 loss on the transcriptome and chromatin accessibility. Nuclear TDP-43 loss is associated with gene expression changes that affect RNA processing, nucleocytoplas-mic transport, histone processing, and DNA damage. Loss of nuclear TDP-43 is also associated with chromatin decondensation around long interspersed nuclear elements (LINEs) and increased LINE1 DNA content. Moreover, loss of TDP-43 leads to increased retrotransposition that can be inhibited with antiretro-viral drugs, suggesting that TDP-43 neuropathology is associated with altered chromatin structure including decondensation of LINEs.

Liu et al. fractionated and sorted for diseased neuronal nuclei from post-mortem FTD-ALS human brains and showed that loss of an RNA-binding protein, TDP-43, altered the transcriptome and chromatin accessibility. Their results suggest that loss of nuclear TDP-43 is associated with decondensation of LINE retrotransposons.

Large-scale analysis of the position-dependent binding and regulation of human RNA binding proteins

Background

RNA binding proteins (RBPs) play essential roles in the regulation of RNA metabolism. Recent studies have disclosed that RBPs achieve their functions via binding to their targets in a position-dependent pattern on RNAs. However, few studies have systematically addressed the associations between the RBP's functions and their positional binding preferences.

Methods

Here, we present large-scale analyses on the functional targets of human RBPs by integrating the enhanced cross-linking and immunoprecipitation followed by sequencing (eCLIP-seq) datasets and the shRNA knockdown followed by RNA-seq datasets that are deposited in the integrated ENCyclopedia of DNA Elements in the human genome (ENCODE) data portal.

Results

We found that (1) binding to the translation termination site and the 3'untranslated region is important to most human RBPs in the RNA decay regulation; (2) RBPs' binding and regulation follow a cell-type specific pattern.

Conclusions

These analysis results show the strong relationship between the binding position and the functions of RBPs, which provides novel insights into the RBPs' regulation mechanisms.

Neuronal Transcriptome from C9orf72 Repeat Expanded Human Tissue is Associated with Loss of C9orf72 Function

A hexanucleotide G₄C₂ repeat expansion in C9orf72 is the most common genetic cause of familial and sporadic cases of amyotrophic lateral sclerosis (ALS) and frontotemporal degeneration (FTD). The mutation is associated with a reduction of C9orf72 protein and accumulation of toxic RNA and dipeptide repeat aggregates. The accumulation of toxic RNA has been proposed to sequester RNA binding proteins thereby altering RNA processing, consistent with previous transcriptome studies that have shown that the C9orf72 repeat expansion is linked to abundant splicing alterations and transcriptome changes. Here, we used a subcellular fractionation method and FACS to enrich for neuronal nuclei from C9orf72 repeat expanded post-mortem human ALS/FTD brains, and to remove neuronal nuclei with TDP-43 pathology which are observed in nearly all symptomatic C9orf72 repeat expanded cases. We show that the C9orf72 expansion is associated with relatively mild gene expression changes. Dysregulated genes were enriched for vesicle transport pathways, which is consistent with the known functions of C9orf72 protein. Further analysis suggests that the C9orf72 transcriptome is not driven by toxic RNA but is rather shaped by the depletion of pathologic TDP-43 nuclei and the loss of C9orf72 expression. These findings argue against RNA binding protein sequestration in neurons as a major contributor to C9orf72 mediated toxicity.

Computational analysis of ribonomics datasets identifies long non-coding RNA targets of γ-herpesviral miRNAs

Ribonomics experiments involving crosslinking and immuno-precipitation (CLIP) of Ago proteins have expanded the understanding of the miRNA targetome of several organisms. These techniques, collectively referred to as CLIP-seq, have been applied to identifying the mRNA targets of miRNAs expressed by Kaposi’s Sarcoma-associated herpes virus (KSHV) and Epstein–Barr virus (EBV). However, these studies focused on identifying only those RNA targets of KSHV and EBV miRNAs that are known to encode proteins. Recent studies have demonstrated that long non-coding RNAs (lncRNAs) are also targeted by miRNAs. In this study, we performed a systematic re-analysis of published datasets from KSHV- and EBV-driven cancers. We used CLIP-seq data from lymphoma cells or EBV-transformed B cells, and a crosslinking, ligation and sequencing of hybrids dataset from KSHV-infected endothelial cells, to identify novel lncRNA targets of viral miRNAs. Here, we catalog the lncRNA targetome of KSHV and EBV miRNAs, and provide a detailed in silico analysis of lncRNA–miRNA binding interactions. Viral miRNAs target several hundred lncRNAs, including a subset previously shown to be aberrantly expressed in human malignancies. In addition, we identified thousands of lncRNAs to be putative targets of human miRNAs, suggesting that miRNA–lncRNA interactions broadly contribute to the regulation of gene expression.

PRAS: Predicting functional targets of RNA binding proteins based on CLIP-seq peaks

RNA-protein interaction plays important roles in post-transcriptional regulation. Recent advancements in cross-linking and immunoprecipitation followed by sequencing (CLIP-seq) technologies make it possible to detect the binding peaks of a given RNA binding protein (RBP) at transcriptome scale. However, it is still challenging to predict the functional consequences of RBP binding peaks. In this study, we propose the Protein-RNA Association Strength (PRAS), which integrates the intensities and positions of the binding peaks of RBPs for functional mRNA targets prediction. We illustrate the superiority of PRAS over existing approaches on predicting the functional targets of two related but divergent CELF (CUGBP, ELAV-like factor) RBPs in mouse brain and muscle. We also demonstrate the potential of PRAS for wide adoption by applying it to the enhanced CLIP-seq (eCLIP) datasets of 37 RNA decay related RBPs in two human cell lines. PRAS can be utilized to investigate any RBPs with available CLIP-seq peaks. PRAS is freely available at http://ouyanglab.jax.org/pras/.

Recent computational developments on CLIP-seq data analysis and microRNA targeting implications

Cross-Linking Immunoprecipitation associated to high-throughput sequencing (CLIP-seq) is a technique used to identify RNA directly bound to RNA-binding proteins across the entire transcriptome in cell or tissue samples. Recent technological and computational advances permit the analysis of many CLIP-seq samples simultaneously, allowing us to reveal the comprehensive network of RNA–protein interaction and to integrate it to other genome-wide analyses. Therefore, the design and quality management of the CLIP-seq analyses are of critical importance to extract clean and biological meaningful information from CLIP-seq experiments. The application of CLIP-seq technique to Argonaute 2 (Ago2) protein, the main component of the microRNA (miRNA)-induced silencing complex, reveals the direct binding sites of miRNAs, thus providing insightful information about the role played by miRNA(s). In this review, we summarize and discuss the most recent computational methods for CLIP-seq analysis, and discuss their impact on Ago2/miRNA-binding site identification and prediction with a regard toward human pathologies.

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/

BayMAP: a Bayesian hierarchical model for the analysis of PAR-CLIP data

MOTIVATION

Photoactivatable-Ribonucleoside-Enhanced Crosslinking and Immunoprecipitation (PAR-CLIP) is a biochemical method for detecting interaction sites of proteins with mRNA. This method introduces T-to-C substitutions at sequenced cDNA that help to detect binding sites on mRNA. However, T-to-C substitutions can also occur due to other reasons such as mismatches or SNPs. Only few statistical procedures exist for detecting binding sites in PAR-CLIP data. Most of these methods do not account for other types of substitutions than those induced by PAR-CLIP, and therefore, also report positions with high T-to-C substitution rates, e.g. SNPs, as binding sites. Moreover, none of these procedures allow to include additional information, e.g. the type of mRNA region, relevant for the biology of microRNA-binding sites.

RESULTS

We have developed BayMAP, a procedure based on a fully Bayesian hierarchical model that takes other sources of substitutions into account. Furthermore, this model enables the incorporation of additional information into the analysis of PAR-CLIP data. This incorporation does not only permit a better detection of binding sites, but also a better understanding of the data and the biology of binding sites. In applications to simulated PAR-CLIP data, BayMAP distinguishes binding sites from noise better than existing methods. Additionally, it yields good estimates of the influence of the additional information. We here demonstrate BayMAP's usability for real datasets even when noisy data is present.

AVAILABILITY AND IMPLEMENTATION

BayMAP is freely available as an R package at http://stat.math.uni-duesseldorf.de/baymap.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Practical considerations on performing and analyzing CLIP-seq experiments to identify transcriptomic-wide RNA-protein interactions

RNA-binding proteins are important players in post-transcriptional regulation, such as modulating mRNA splicing, translation, and degradation under diverse biological settings. Identifying and characterizing the RNA substrates is a critical step in deciphering the function and molecular mechanisms of the target RNA-binding proteins. High-throughput sequencing of the RNA fragments isolated by crosslinking immunoprecipitation (CLIP-seq) is one of the standard techniques to identify the in vivo transcriptome-wide binding sites of the target RNA-binding protein. This method is widely used in functional and mechanistic characterizations of RNA-binding proteins. In this review, we provide several practical considerations on performing and analyzing CLIP-seq experiments. Particularly, we focus on how to perform CLIP-seq experiments on endogenous RNA-binding proteins. In addition, we provide a practical summary on how to choose and use computational pipelines from an increasing number of computational methods and packages that are available for analyzing the sequencing datasets from the CLIP-seq experiments. We hope these practical considerations will facilitate experimental biologists in performing and analyzing CLIP-seq experiment to obtain biologically relevant mechanistic insights.

Computational Analysis of RNA-Protein Interactions via Deep Sequencing

RNA-binding proteins (RBPs) function in all aspects of RNA processes including stability, structure, export, localization and translation, and control gene expression at the posttranscriptional level. To investigate the roles of RBPs and their direct RNA ligands in vivo, recent global approaches combining RNA immunoprecipitation and deep sequencing (RIP-seq) as well as UV-cross-linking (CLIP-seq) have become instrumental in dissecting RNA-protein interactions. However, the computational analysis of these high-throughput sequencing data is still challenging. Here, we provide a computational pipeline to analyze CLIP-seq and RIP-seq datasets. This generic analytic procedure may help accelerate the identification of direct RNA-protein interactions from high-throughput RBP profiling experiments in a variety of bacterial species.

Influenza Virus NS1 Protein-RNA Interactome Reveals Intron Targeting

The NS1 protein of influenza A virus is a multifunctional virulence factor that inhibits cellular processes to facilitate viral gene expression. While NS1 is known to interact with RNA and proteins to execute these functions, the cellular RNAs that physically interact with NS1 have not been systematically identified. Here we reveal a NS1 protein-RNA interactome and show that NS1 primarily binds intronic sequences. Among this subset of pre-mRNAs is the RIG-I pre-mRNA, which encodes the main cytoplasmic antiviral sensor of influenza virus infection. This suggested that NS1 interferes with the antiviral response at a posttranscriptional level by virtue of its RNA binding properties. Indeed, we show that NS1 is necessary in the context of viral infection and sufficient upon transfection to decrease the rate of RIG-I intron removal. This NS1 function requires a functional RNA binding domain and is independent of the NS1 interaction with the cleavage and polyadenylation specificity factor CPSF30. NS1 has been previously shown to abrogate RIG-I-mediated antiviral immunity by inhibiting its protein function. Our data further suggest that NS1 also posttranscriptionally alters RIG-I pre-mRNA processing by binding to the RIG-I pre-mRNA.IMPORTANCE A key virulence factor of influenza A virus is the NS1 protein, which inhibits various cellular processes to facilitate viral gene expression. The NS1 protein is localized in the nucleus and in the cytoplasm during infection. In the nucleus, NS1 has functions related to inhibition of gene expression that involve protein-protein and protein-RNA interactions. While several studies have elucidated the protein interactome of NS1, we still lack a clear and systematic understanding of the NS1-RNA interactome. Here we reveal a nuclear NS1-RNA interactome and show that NS1 primarily binds intronic sequences within a subset of pre-mRNAs, including the RIG-I pre-mRNA that encodes the main cytoplasmic antiviral sensor of influenza virus infection. Our data here further suggest that NS1 is necessary and sufficient to impair intron processing of the RIG-I pre-mRNA. These findings support a posttranscriptional role for NS1 in the inhibition of RIG-I expression.

microCLIP super learning framework uncovers functional transcriptome-wide miRNA interactions

Argonaute crosslinking and immunoprecipitation (CLIP) experiments are the most widely used high-throughput methodologies for miRNA targetome characterization. The analysis of Photoactivatable Ribonucleoside-Enhanced (PAR) CLIP methodology focuses on sequence clusters containing T-to-C conversions. Here, we demonstrate for the first time that the non-T-to-C clusters, frequently observed in PAR-CLIP experiments, exhibit functional miRNA-binding events and strong RNA accessibility. This discovery is based on the analysis of an extensive compendium of bona fide miRNA-binding events, and is further supported by numerous miRNA perturbation experiments and structural sequencing data. The incorporation of these previously neglected clusters yields an average of 14% increase in miRNA-target interactions per PAR-CLIP library. Our findings are integrated in microCLIP (www.microrna.gr/microCLIP), a cutting-edge framework that combines deep learning classifiers under a super learning scheme. The increased performance of microCLIP in CLIP-Seq-guided detection of miRNA interactions, uncovers previously elusive regulatory events and miRNA-controlled pathways.

AGO-PAR-CLIP is widely used for high-throughput miRNA target characterization. Here, the authors show that the previously neglected non-T-to-C clusters denote functional miRNA binding events, and develop microCLIP, a super learning framework that accurately detects miRNA interactions.

Biological functions, regulatory mechanisms, and disease relevance of RNA localization pathways

The asymmetric subcellular distribution of RNA molecules from their sites of transcription to specific compartments of the cell is an important aspect of post-transcriptional gene regulation. This involves the interplay of intrinsic cis-regulatory elements within the RNA molecules with trans-acting RNA-binding proteins and associated factors. Together, these interactions dictate the intracellular localization route of RNAs, whose downstream impacts have wide-ranging implications in cellular physiology. In this review, we examine the mechanisms underlying RNA localization and discuss their biological significance. We also review the growing body of evidence pointing to aberrant RNA localization pathways in the development and progression of diseases.

Mapping transcriptome-wide protein-RNA interactions to elucidate RNA regulatory programs

BACKGROUND

Our understanding of post-transcriptional gene regulation has increased exponentially with the development of robust methods to define protein-RNA interactions across the transcriptome. In this review, we highlight the evolution and successful applications of crosslinking and immunoprecipitation (CLIP) methods to interrogate protein-RNA interactions in a transcriptome-wide manner.

RESULTS

Here, we survey the vast array of in vitro and in vivo approaches used to identify protein-RNA interactions, including but not limited to electrophoretic mobility shift assays, systematic evolution of ligands by exponential enrichment (SELEX), and RIP-seq. We particularly emphasize the advancement of CLIP technologies, and detail protocol improvements and computational tools used to analyze the output data. Importantly, we discuss how profiling protein-RNA interactions can delineate biological functions including splicing regulation, alternative polyadenylation, cytoplasmic decay substrates, and miRNA targets.

CONCLUSIONS

In summary, this review summarizes the benefits of characterizing RNA-protein networks to further understand the regulation of gene expression and disease pathogenesis. Our review comments on how future CLIP technologies can be adapted to address outstanding questions related to many aspects of RNA metabolism and further advance our understanding of RNA biology.

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.

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.

Multiple nuclear-replicating viruses require the stress-induced protein ZC3H11A for efficient growth

There is a strong need for the development of new antiviral therapies, and this study sheds light on a host–virus interaction that is significant for a number of medically important human viruses. The study also suggests that the RNA-binding zinc finger CCCH-type containing 11A (ZC3H11A) protein takes part in a mechanism that facilitates nuclear export of mRNA, particularly under cellular stress, a mechanism that has been “hijacked” by several nuclear-replicating viruses to promote their replication. ZC3H11A is therefore a potential target for development of an antiviral therapy.

The zinc finger CCCH-type containing 11A (ZC3H11A) gene encodes a well-conserved zinc finger protein that may function in mRNA export as it has been shown to associate with the transcription export (TREX) complex in proteomic screens. Here, we report that ZC3H11A is a stress-induced nuclear protein with RNA-binding capacity that localizes to nuclear splicing speckles. During an adenovirus infection, the ZC3H11A protein and splicing factor SRSF2 relocalize to nuclear regions where viral DNA replication and transcription take place. Knockout (KO) of ZC3H11A in HeLa cells demonstrated that several nuclear-replicating viruses are dependent on ZC3H11A for efficient growth (HIV, influenza virus, herpes simplex virus, and adenovirus), whereas cytoplasmic replicating viruses are not (vaccinia virus and Semliki Forest virus). High-throughput sequencing of ZC3H11A–cross-linked RNA showed that ZC3H11A binds to short purine-rich ribonucleotide stretches in cellular and adenoviral transcripts. We show that the RNA-binding property of ZC3H11A is crucial for its function and localization. In ZC3H11A KO cells, the adenovirus fiber mRNA accumulates in the cell nucleus. Our results suggest that ZC3H11A is important for maintaining nuclear export of mRNAs during stress and that several nuclear-replicating viruses take advantage of this mechanism to facilitate their replication.

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.

PAR-CLIP for Discovering Target Sites of RNA-Binding Proteins

RNA-binding proteins (RBPs) establish posttranscriptional gene regulation (PTGR) by coordinating the maturation, editing, transport, stability, and translation of cellular RNAs. A variety of experimental approaches have been developed to characterize the RNAs associated with RBPs in vitro as well as in vivo. Our laboratory developed Photoactivatable-Ribonucleoside-Enhanced Cross-Linking and Immunoprecipitation (PAR-CLIP), which in combination with next-generation sequencing enables the identification of RNA targets of RBPs at a nucleotide-level resolution. Here we present an updated and condensed step-by-step PAR-CLIP protocol followed by the description of our RNA-seq data analysis pipeline.

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.

Global role of the bacterial post-transcriptional regulator CsrA revealed by integrated transcriptomics

CsrA is a post-transcriptional regulatory protein that is widely distributed among bacteria. This protein influences bacterial lifestyle decisions by binding to the 5′ untranslated and/or early coding regions of mRNA targets, causing changes in translation initiation, RNA stability, and/or transcription elongation. Here, we assess the contribution of CsrA to gene expression in Escherichia coli on a global scale. UV crosslinking immunoprecipitation and sequencing (CLIP-seq) identify RNAs that interact directly with CsrA in vivo, while ribosome profiling and RNA-seq uncover the impact of CsrA on translation, RNA abundance, and RNA stability. This combination of approaches reveals unprecedented detail about the regulatory role of CsrA, including novel binding targets and physiological roles, such as in envelope function and iron homeostasis. Our findings highlight the integration of CsrA throughout the E. coli regulatory network, where it orchestrates vast effects on gene expression.

The RNA-binding protein CsrA regulates the expression of hundreds of bacterial genes. Here, Potts et al. use several approaches to assess the contribution of CsrA to global gene expression in E. coli, revealing new binding targets and physiological roles such as in envelope function and iron homeostasis.

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 Tool Kit (CTK): a flexible and robust pipeline to analyze CLIP sequencing data

Summary

UV cross-linking and immunoprecipitation (CLIP), followed by high-throughput sequencing, is a powerful biochemical assay that maps in vivo protein-RNA interactions on a genome-wide scale. The CLIP Tool Kit (CTK) aims at providing a set of tools for flexible, streamlined and comprehensive CLIP data analysis. This software package extends the scope of our original CIMS package.

Availability and Implementation

The software is implemented in Perl. The source code and detailed documentation are available at http://zhanglab.c2b2.columbia.edu/index.php/CTK .

Contact

cz2294@columbia.edu.

Dysregulation of spliceosome gene expression in advanced prostate cancer by RNA-binding protein PSF

Developing therapeutic approaches are necessary for treating hormone-refractory prostate cancer. Activation of androgen receptor (AR) and its variants' expression along with the downstream signals are mostly important for disease progression. However, the mechanism for marked increases of AR signals and its expression is still unclear. Here, we revealed that various spliceosome genes are aberrantly induced by RNA-binding protein PSF, leading to enhancement of the splicing activities for AR expression. Our high-speed sequence analyses identified global PSF-binding transcripts. PSF was shown to stabilize and activate key long noncoding RNAs and AR-regulated gene expressions in prostate cancer cells. Interestingly, mRNAs of spliceosome-related genes are putative primary targets of PSF. Their gene expressions are up-regulated by PSF in hormone-refractory prostate cancer. Moreover, PSF coordinated these spliceosome proteins to form a complex to promote AR splicing and expression. Thus, targeting PSF and its related pathways implicates the therapeutic possibility for hormone-refractory prostate cancer.

From benchmarking HITS-CLIP peak detection programs to a new method for identification of miRNA-binding sites from Ago2-CLIP data

Experimental evidence indicates that about 60% of miRNA-binding activity does not follow the canonical rule about the seed matching between miRNA and target mRNAs, but rather a non-canonical miRNA targeting activity outside the seed or with a seed-like motifs. Here, we propose a new unbiased method to identify canonical and non-canonical miRNA-binding sites from peaks identified by Ago2 Cross-Linked ImmunoPrecipitation associated to high-throughput sequencing (CLIP-seq). Since the quality of peaks is of pivotal importance for the final output of the proposed method, we provide a comprehensive benchmarking of four peak detection programs, namely CIMS, PIPE-CLIP, Piranha and Pyicoclip, on four publicly available Ago2-HITS-CLIP datasets and one unpublished in-house Ago2-dataset in stem cells. We measured the sensitivity, the specificity and the position accuracy toward miRNA binding sites identification, and the agreement with TargetScan. Secondly, we developed a new pipeline, called miRBShunter, to identify canonical and non-canonical miRNA-binding sites based on de novo motif identification from Ago2 peaks and prediction of miRNA::RNA heteroduplexes. miRBShunter was tested and experimentally validated on the in-house Ago2-dataset and on an Ago2-PAR-CLIP dataset in human stem cells. Overall, we provide guidelines to choose a suitable peak detection program and a new method for miRNA-target identification.

Seten: a tool for systematic identification and comparison of processes, phenotypes, and diseases associated with RNA-binding proteins from condition-specific CLIP-seq profiles

RNA-binding proteins (RBPs) control the regulation of gene expression in eukaryotic genomes at post-transcriptional level by binding to their cognate RNAs. Although several variants of CLIP (crosslinking and immunoprecipitation) protocols are currently available to study the global protein-RNA interaction landscape at single-nucleotide resolution in a cell, currently there are very few tools that can facilitate understanding and dissecting the functional associations of RBPs from the resulting binding maps. Here, we present Seten, a web-based and command line tool, which can identify and compare processes, phenotypes, and diseases associated with RBPs from condition-specific CLIP-seq profiles. Seten uses BED files resulting from most peak calling algorithms, which include scores reflecting the extent of binding of an RBP on the target transcript, to provide both traditional functional enrichment as well as gene set enrichment results for a number of gene set collections including BioCarta, KEGG, Reactome, Gene Ontology (GO), Human Phenotype Ontology (HPO), and MalaCards Disease Ontology for several organisms including fruit fly, human, mouse, rat, worm, and yeast. It also provides an option to dynamically compare the associated gene sets across data sets as bubble charts, to facilitate comparative analysis. Benchmarking of Seten using eCLIP data for IGF2BP1, SRSF7, and PTBP1 against their corresponding CRISPR RNA-seq in K562 cells as well as randomized negative controls, demonstrated that its gene set enrichment method outperforms functional enrichment, with scores significantly contributing to the discovery of true annotations. Comparative performance analysis using these CRISPR control data sets revealed significantly higher precision and comparable recall to that observed using ChIP-Enrich. Seten's web interface currently provides precomputed results for about 200 CLIP-seq data sets and both command line as well as web interfaces can be used to analyze CLIP-seq data sets. We highlight several examples to show the utility of Seten for rapid profiling of various CLIP-seq data sets. Seten is available on http://www.iupui.edu/∼sysbio/seten/.

Optimization of PAR-CLIP for transcriptome-wide identification of binding sites of RNA-binding proteins

Photoactivatable-Ribonucleoside-Enhanced Crosslinking and Immunoprecipitation (PAR-CLIP) in combination with next-generation sequencing is a powerful method for identifying endogenous targets of RNA-binding proteins (RBPs). Depending on the characteristics of each RBP, key steps in the PAR-CLIP procedure must be optimized. Here we present a comprehensive step-by-step PAR-CLIP protocol with detailed explanations of the critical steps. Furthermore, we report the application of a new PAR-CLIP data analysis pipeline to three distinct RBPs targeting different annotation categories of cellular RNAs.

Six GU-rich (6GUR) FUS-binding motifs detected by normalization of CLIP-seq by Nascent-seq

FUS, an RNA-binding protein (RBP), is mutated or abnormally regulated in neurodegenerative disorders. FUS regulates various aspects of RNA metabolisms. FUS-binding sites are rich in GU contents and are highly degenerative. FUS-binding motifs of GGU, GGUG, GUGGU and CGCGC have been previously reported. These motifs, however, are applicable to a small fraction of FUS-binding sites. As CLIP-seq tags are enriched in genes that are highly expressed, we normalized CLIP-seq tags by Nascent-seq tags or RNA-seq tags of mouse N2a cells. Nascent-seq identifies nascent transcripts before being processed for splicing and polyadenylation. We extracted frequently observed 4-nt motifs from Nascent-seq-normalized CLIP regions, RNA-seq-normalized CLIP regions, and native CLIP regions. Specific GU-rich motifs were best detected in Nascent-seq-normalized CLIP regions. Analysis of structural motifs using Nascent-seq-normalized CLIP regions also predicted GU-rich sequence forming a stem structure. Sensitivity and specificity were calculated by examining whether the extracted motifs were present at the cross-linking-induced mutation sites (CIMS), where FUS was directly bound. We found that a combination of six motifs (UGUG, CUGG, UGGU, GCUG, GUGG, and UUGG), which were extracted from Nascent-seq-normalized CLIP-regions, had a better discriminative power than (i) motifs extracted from RNA-seq-normalized CLIP regions, (ii) motifs extracted from native CLIP regions, (iii) previously reported individual motifs, or (iv) 15 motifs in SpliceAid 2. Validation of the 6 GU-rich (6GU_R) motifs using CLIP-seq of the cerebrum and the whole brain showed that the 6GU_R motifs were specifically enriched in CIMS. The number of the 6GU_R motifs in an uninterrupted region was counted and multiplied by four to calculate the area, which was defined as the 6GU_R-Score. The 6GU_R-Score of 8 or more best discriminated CIMS from CIMS-flanking regions. We propose that the 6GU_R motifs predict FUS-binding sites more efficiently than previously reported individual motifs or 15 motifs in SpliceAid 2.

Computational and Experimental Identification of Tissue-Specific MicroRNA Targets

In this chapter we discuss computational methods for the prediction of microRNA (miRNA) targets. More specifically, we consider machine learning-based approaches and explain why these methods have been relatively unsuccessful in reducing the number of false positive predictions. Further we suggest approaches designed to improve their performance by considering tissue-specific target regulation. We argue that the miRNA targetome differs depending on the tissue type and introduce a novel algorithm that predicts miRNA targets specifically for colorectal cancer. We discuss features of miRNAs and target sites that affect target recognition, and how next-generation sequencing data can support the identification of novel miRNAs, differentially expressed miRNAs and their tissue-specific mRNA targets. In addition, we introduce some experimental approaches for the validation of miRNA targets as well as web-based resources sharing predicted and validated miRNA target interactions.

Illustrating the Epitranscriptome at Nucleotide Resolution Using Methylation-iCLIP (miCLIP)

Next-generation sequencing technologies have enabled the transcriptome to be profiled at a previously unprecedented speed and depth. This yielded insights into fundamental transcriptomic processes such as gene transcription, RNA processing, and mRNA splicing. Immunoprecipitation-based transcriptomic methods such as individual nucleotide resolution crosslinking immunoprecipitation (iCLIP) have also allowed high-resolution analysis of the RNA interactions of a protein of interest, thus revealing new regulatory mechanisms. We and others have recently modified this method to profile RNA methylation, and we refer to this customized technique as methylation-iCLIP (miCLIP). Variants of miCLIP have been used to map the methyl-5-cytosine (m5C) or methyl-6-adenosine (m6A) modification at nucleotide resolution in the human transcriptome. Here we describe the m5C-miCLIP protocol, discuss how it yields the nucleotide-resolution RNA modification maps, and comment on how these have contributed to the new field of molecular genetics research coined "epitranscriptomics."

Bioinformatic tools for analysis of CLIP ribonucleoprotein data

Investigating the interactions of RNA-binding proteins (RBPs) with RNAs is a complex task for molecular and computational biologists. The molecular biology techniques and the computational approaches to understand RBP-RNA (or ribonucleoprotein, RNP) interactions have advanced considerably over the past few years and numerous and diverse software tools have been developed to analyze these data. Accordingly, laboratories interested in RNP biology face the challenge of choosing adequately among the available software tools those that best address the biological problem they are studying. Here, we focus on state-of-the-art molecular biology techniques that employ crosslinking and immunoprecipitation (CLIP) of an RBP to study and map RNP interactions. We review the different software tools and databases available to analyze the most widely used CLIP methods, HITS-CLIP, PAR-CLIP, and iCLIP. WIREs RNA 2017, 8:e1404. doi: 10.1002/wrna.1404 For further resources related to this article, please visit the WIREs website.

Local synthesis of dynein cofactors matches retrograde transport to acutely changing demands

Cytoplasmic dynein mediates retrograde transport in axons, but it is unknown how its transport characteristics are regulated to meet acutely changing demands. We find that stimulus-induced retrograde transport of different cargos requires the local synthesis of different dynein cofactors. Nerve growth factor (NGF)-induced transport of large vesicles requires local synthesis of Lis1, while smaller signalling endosomes require both Lis1 and p150^Glued. Lis1 synthesis is also triggered by NGF withdrawal and required for the transport of a death signal. Association of Lis1 transcripts with the microtubule plus-end tracking protein APC is required for their translation in response to NGF stimulation but not for their axonal recruitment and translation upon NGF withdrawal. These studies reveal a critical role for local synthesis of dynein cofactors for the transport of specific cargos and identify association with RNA-binding proteins as a mechanism to establish functionally distinct pools of a single transcript species in axons.

The molecular mechanisms underlying retrograde transport in axons are only partially understood. Villarin et al. show that in cultured DRG neurons, extracellular trophic cues such as NGF dynamically regulate local protein synthesis of dynein cofactors, thus controlling retrograde trafficking in neurons.

Integrating Epigenomics into the Understanding of Biomedical Insight

Epigenetics is one of the most rapidly expanding fields in biomedical research, and the popularity of the high-throughput next-generation sequencing (NGS) highlights the accelerating speed of epigenomics discovery over the past decade. Epigenetics studies the heritable phenotypes resulting from chromatin changes but without alteration on DNA sequence. Epigenetic factors and their interactive network regulate almost all of the fundamental biological procedures, and incorrect epigenetic information may lead to complex diseases. A comprehensive understanding of epigenetic mechanisms, their interactions, and alterations in health and diseases genome widely has become a priority in biological research. Bioinformatics is expected to make a remarkable contribution for this purpose, especially in processing and interpreting the large-scale NGS datasets. In this review, we introduce the epigenetics pioneering achievements in health status and complex diseases; next, we give a systematic review of the epigenomics data generation, summarize public resources and integrative analysis approaches, and finally outline the challenges and future directions in computational epigenomics.

Antisense transcription licenses nascent transcripts to mediate transcriptional gene silencing

In this study, Dang et al. use Neurospora to demonstrate a critical role for transcription kinetics in long noncoding RNA-mediated epigenetic modifications and identify ERI-1 as an important regulator of cotranscriptional gene silencing and post-transcriptional RNA metabolism.

In eukaryotes, antisense transcription can regulate sense transcription by induction of epigenetic modifications. We showed previously that antisense transcription triggers Dicer-independent siRNA (disiRNA) production and disiRNA locus DNA methylation (DLDM) in Neurospora crassa. Here we show that the conserved exonuclease ERI-1 (enhanced RNAi-1) is a critical component in this process. Antisense transcription and ERI-1 binding to target RNAs are necessary and sufficient to trigger DLDM. Convergent transcription causes stalling of RNA polymerase II during transcription, which permits ERI-1 to bind nascent RNAs in the nucleus and recruit a histone methyltransferase complex that catalyzes chromatin modifications. Furthermore, we show that, in the cytoplasm, ERI-1 targets hundreds of transcripts from loci without antisense transcription to regulate RNA stability. Together, our results demonstrate a critical role for transcription kinetics in long noncoding RNA-mediated epigenetic modifications and identify ERI-1 as an important regulator of cotranscriptional gene silencing and post-transcriptional RNA metabolism.

AnnoLnc: a web server for systematically annotating novel human lncRNAs

BACKGROUND

Long noncoding RNAs (lncRNAs) have been shown to play essential roles in almost every important biological process through multiple mechanisms. Although the repertoire of human lncRNAs has rapidly expanded, their biological function and regulation remain largely elusive, calling for a systematic and integrative annotation tool.

RESULTS

Here we present AnnoLnc ( http://annolnc.cbi.pku.edu.cn ), a one-stop portal for systematically annotating novel human lncRNAs. Based on more than 700 data sources and various tool chains, AnnoLnc enables a systematic annotation covering genomic location, secondary structure, expression patterns, transcriptional regulation, miRNA interaction, protein interaction, genetic association and evolution. An intuitive web interface is available for interactive analysis through both desktops and mobile devices, and programmers can further integrate AnnoLnc into their pipeline through standard JSON-based Web Service APIs.

CONCLUSIONS

To the best of our knowledge, AnnoLnc is the only web server to provide on-the-fly and systematic annotation for newly identified human lncRNAs. Compared with similar tools, the annotation generated by AnnoLnc covers a much wider spectrum with intuitive visualization. Case studies demonstrate the power of AnnoLnc in not only rediscovering known functions of human lncRNAs but also inspiring novel hypotheses.

rMAPS: RNA map analysis and plotting server for alternative exon regulation

RNA-binding proteins (RBPs) play a critical role in the regulation of alternative splicing (AS), a prevalent mechanism for generating transcriptomic and proteomic diversity in eukaryotic cells. Studies have shown that AS can be regulated by RBPs in a binding-site-position dependent manner. Depending on where RBPs bind, splicing of an alternative exon can be enhanced or suppressed. Therefore, spatial analyses of RBP motifs and binding sites around alternative exons will help elucidate splicing regulation by RBPs. The development of high-throughput sequencing technologies has allowed transcriptome-wide analyses of AS and RBP–RNA interactions. Given a set of differentially regulated alternative exons obtained from RNA sequencing (RNA-seq) experiments, the rMAPS web server (http://rmaps.cecsresearch.org) performs motif analyses of RBPs in the vicinity of alternatively spliced exons and creates RNA maps that depict the spatial patterns of RBP motifs. Similarly, rMAPS can also perform spatial analyses of RBP–RNA binding sites identified by cross-linking immunoprecipitation sequencing (CLIP-seq) experiments. We anticipate rMAPS will be a useful tool for elucidating RBP regulation of alternative exon splicing using high-throughput sequencing data.

Tiny giants of gene regulation: experimental strategies for microRNA functional studies

The discovery over two decades ago of short regulatory microRNAs (miRNAs) has led to the inception of a vast biomedical research field dedicated to understanding these powerful orchestrators of gene expression. Here we aim to provide a comprehensive overview of the methods and techniques underpinning the experimental pipeline employed for exploratory miRNA studies in animals. Some of the greatest challenges in this field have been uncovering the identity of miRNA-target interactions and deciphering their significance with regard to particular physiological or pathological processes. These endeavors relied almost exclusively on the development of powerful research tools encompassing novel bioinformatics pipelines, high-throughput target identification platforms, and functional target validation methodologies. Thus, in an unparalleled manner, the biomedical technology revolution unceasingly enhanced and refined our ability to dissect miRNA regulatory networks and understand their roles in vivo in the context of cells and organisms. Recurring motifs of target recognition have led to the creation of a large number of multifactorial bioinformatics analysis platforms, which have proved instrumental in guiding experimental miRNA studies. Subsequently, the need for discovery of miRNA-target binding events in vivo drove the emergence of a slew of high-throughput multiplex strategies, which now provide a viable prospect for elucidating genome-wide miRNA-target binding maps in a variety of cell types and tissues. Finally, deciphering the functional relevance of miRNA post-transcriptional gene silencing under physiological conditions, prompted the evolution of a host of technologies enabling systemic manipulation of miRNA homeostasis as well as high-precision interference with their direct, endogenous targets. For further resources related to this article, please visit the WIREs website.