CLIPdb: a CLIP-seq database for protein-RNA interactions
- Yu-Cheng T Yang†1,
- Chao Di†1,
- Boqin Hu†1,
- Meifeng Zhou1,
- Yifang Liu1,
- Nanxi Song1,
- Yang Li1,
- Jumpei Umetsu1, 2 and
- Zhi John Lu1Email author
© Yang et al.; licensee BioMed Central. 2015
Received: 25 October 2014
Accepted: 22 January 2015
Published: 5 February 2015
RNA-binding proteins (RBPs) play essential roles in gene expression regulation through their interactions with RNA transcripts, including coding, canonical non-coding and long non-coding RNAs. Large amounts of crosslinking immunoprecipitation (CLIP)-seq data (including HITS-CLIP, PAR-CLIP, and iCLIP) have been recently produced to reveal transcriptome-wide binding sites of RBPs at the single-nucleotide level.
Here, we constructed a database, CLIPdb, to describe RBP-RNA interactions based on 395 publicly available CLIP-seq data sets for 111 RBPs from four organisms: human, mouse, worm and yeast. We consistently annotated the CLIP-seq data sets and RBPs, and developed a user-friendly interface for rapid navigation of the CLIP-seq data. We applied a unified computational method to identify transcriptome-wide binding sites, making the binding sites directly comparable and the data available for integration across different CLIP-seq studies. The high-resolution binding sites of the RBPs can be visualized on the whole-genome scale using a browser. In addition, users can browse and download the identified binding sites of all profiled RBPs by querying genes of interest, including both protein coding genes and non-coding RNAs.
Manually curated metadata and uniformly identified binding sites of publicly available CLIP-seq data sets will be a foundation for further integrative and comparative analyses. With maintained up-to-date data sets and improved functionality, CLIPdb (http://clipdb.ncrnalab.org) will be a valuable resource for improving the understanding of post-transcriptional regulatory networks.
KeywordsCLIP-seq RNA-binding protein RNA Regulatory networks Data integration
RNA binding proteins (RBPs) play essential roles in the co-transcriptional and post-transcriptional regulation of gene expression [1-3]. In recent years, much progress has been made in the field of ribonomics, which uses high-throughput technologies to investigate the interactions between RBPs and their target RNAs, including coding and non-coding RNAs, in a quantitative and high-resolution manner. With the development of next-generation sequencing technologies, crosslinking immunoprecipitation (CLIP)-seq technology , which includes high-throughput sequencing (HITS)-CLIP , photoactivatable ribonucleoside-enhanced (PAR)-CLIP  and individual-nucleotide resolution crosslinking immunoprecipitation (iCLIP) , has become a powerful tool to study the transcriptome-wide in vivo binding sites of RBPs at the single-nucleotide level. CLIP-seq provides higher resolution than previous technologies (e.g., RNA Immunoprecipitation) for identifying protein binding sites on RNAs.
Approximately 400 CLIP-seq samples from about 80 publications are publicly available from various model organisms. However, the many inconsistencies existing in the metadata annotations for these samples, such as RBP names, tissue types, cell types and disease states, create challenges for public data sharing and reuse. Notably, the analyses of most CLIP-seq studies emphasize only the binding signatures of the RBPs on mRNAs, whereas the binding information is often ignored for long non-coding RNAs (lncRNAs) in raw CLIP-seq data. With the emerging important roles of lncRNAs [8-10], a database without metadata annotation inconsistencies will help researchers to examine the post-transcriptional regulatory mechanisms of lncRNAs.
Several resources, such as CLIPZ , doRiNA , starBase v2.0  and AURA 2 , have provided annotations for some CLIP-seq studies and locations of RBP binding sites. However, most of these resources focus only on a small number of the existing CLIP-seq experiments. For instance, doRiNA contains data for 27 human RBPs, AURA 2 provides data for 47 RBPs, and starBase v2.0 provides data for 49 RBPs. Other databases, such as NPInter v2.0, which focuses on ncRNA-centric interactions , integrate various resources (e.g., starBase). Therefore, although these integrative databases cover more RBPs and species and have a broader perspective on various interactions, they usually include fewer CLIP-seq data sets and provide few high-resolution binding sites on RNAs. In addition, the heterogeneous results obtained using different resources hinder further comparison and integration of the data sets across different cell lines, tissues and developmental stages. Thus, a more complete and better-curated database of CLIP-seq studies is needed.
Construction and content
Collection and annotation of published CLIP-seq data sets
We curated published CLIP-seq data sets for four model organisms, Homo sapiens, Mus musculus, Caenorhabditis elegans and Saccharomyces cerevisiae, from public data repositories, including Gene Expression Omnibus (GEO), European Nucleotide Archive (ENA) and authors’ websites. In total, CLIPdb collected 395 CLIP-seq samples. We consistently annotated the samples with the following categories: RBP factor, species, tissue type, material (i.e., cell line, tissue and cell), cell type, cell line ID, disease state, treatment, assay method, data accession and reference.
Pre-processing of CLIP-seq data sets
Different CLIP libraries may have very different complexities and signal-to-noise ratios depending on many experimental factors. We used multi-step and uniform procedures to pre-process the collected raw CLIP-seq data. First, we combined sequencing runs from the same sample. Then, we trimmed the adaptor and barcode sequences from the raw reads using the FASTX-Toolkit package . We set a stringent read quality for each data set, retaining only those reads with a quality score above 20 in 80% of their nucleotides, and we also restricted the read length to 13 nucleotides (nt) after adapter trimming. Finally, we collapsed identical reads to minimize PCR duplicates. Detailed information for the CLIP-seq data sets we used and processed is shown in Additional file 1.
Identification of transcriptome-wide binding sites
After pre-processing, retained reads from all samples were aligned to their respective genomes (human, hg19; mouse, mm10; yeast, R64-1-1; worm, ws220) using Bowtie-1.0.0 . We retained only those reads with unique mapping locations by setting the Bowtie parameter as “-m 1 --best --strata.” To identify RBP binding sites, we applied Piranha (v1.2.0)  to all the samples, using the following parameters: −b 20 -d ZeroTruncatedNegativeBinomial -p 0.01 or 0.001. Piranha is a statistically robust and flexible computational method for identifying binding sites from CLIP-seq data, and is applicable to all variations of the CLIP-seq technologies. Piranha allows for the direct comparison of binding sites across different cell lines or tissues . Notably, the binding sites obtained from Piranha software are not strand-specific because the current Piranha version does not provide strandness information. In addition, for some CLIP-seq samples (14 data sets), the raw data could not be easily accessed or were of such low quality, according to our results, that we did not provide their binding sites (see Additional file 1).
In addition to the peaks indiscriminately called by Piranha, we also provided the binding peaks called by specialized tools for different CLIP-seq technologies. For PAR-CLIP, we used PARalyzer to identify the T-C transitions . For HITS-CLIP, we used CIMS to identify the mutation sites generated by UV crosslinking [20,21]. For iCLIP, we used CIMS/CITS (included in CIMS v1.0.3) to discriminate the truncation sties . The default parameters were used when applying these specialized tools to identify binding sites.
Annotation for binding sites of RBPs
We first used Gencode human (version 19), Gencode mouse (version M2), SGD (version R64-1-1) and WormBase (version ws220) to annotate the binding sites in the human, mouse, yeast and worm, respectively. In addition, for miRNAs, we used the latest version of annotations from miRBase . For transposable elements (TEs) in the human and mouse, we downloaded the annotation from UCSC Genome Browser. For lncRNAs of the worm, we used the annotations from Nam et al., 2012 . Binding sites with a significance cutoff of 0.01 by Piranha were overlapped with the genomic annotations of the coding sequences (CDS), 5′ untranslated regions (UTRs), 3′ UTRs, TEs, introns, lncRNAs, microRNA primary transcripts, canonical ncRNAs (snRNA, snoRNA, rRNA, tRNA, 7SK_RNA, Y_RNA), intergenic regions and other locations. Here an “intergenic region” denoted a constant distance from any genic region (coding genes, ncRNAs, TEs and pseudogenes), fixed at 2,000 nt for the human and mouse, and 500 nt for the worm and yeast. When 50% of the nucleotides in the binding sites overlapped with the genomic annotations, they were correspondingly annotated. We annotated the binding sites for the RBPs based on the following priority scheme: 3′ UTR, CDS, 5′ UTR, intron, microRNA primary transcript, canonical ncRNA, lncRNA, TE, pseudogene, intergenic region and others.
Annotation of RBPs
The RBPs corresponding to these CLIP-seq data were further annotated with detailed information that was manually curated from other databases, such as NCBI , CISBP-RNA  and RBPDB . Each RBP contained the following terms: the gene names for the RBP factors, species, Ensembl ID, synonyms and RNA binding domain. The RNA binding motifs of some RBPs were depicted in greater detail, including the RNA recognition motif sequence, position weight matrix of the motif, assay method and source literature.
In addition, all RBPs were classified into three major functional groups according to their molecular functional annotations provided in NCBI: (i) splicing factors, (ii) 3′ UTR/poly(A) binding factors and (iii) microRNA binding factors. The remaining RBPs were not grouped. We annotated the functions of the RBPs based on the following priority scheme: splicing factors, 3′ UTR/poly(A) binding factors, microRNA binding factors, and others (Additional file 2).
Utility and discussion
Basic characteristics of CLIPdb
From statistical analyses, we observed that the total number of CLIP-seq studies published increased dramatically after 2008 and exhibited a gradual upward trend thereafter (Figure 2C). These results suggest that the investigation of protein-RNA interactions and CLIP-seq technology is becoming increasingly popular and important within the scientific community. We anticipate that more CLIP-seq studies in various cell types and tissues will be published in the future.
We analyzed the distribution of the genomic elements for the RBP binding sites from each sample in the four species (summarized in Figure 2D and Additional file 3: Figures S1 and S2). After pooling the binding sites in all RBPs for each species, we found that the human and mouse exhibited similar genomic elements, suggesting that functional binding patterns are conserved between mammals. The 3′ UTRs have more binding sites than the 5′ UTRs across all four species. Interestingly, this trend was most obvious for humans and mice, suggesting that 3′ UTR binding may have important regulatory functions in higher species. In addition to protein-coding regions, approximately 2-6% of all binding sites were located in mammalian lncRNAs, indicating that RBP binding may regulate the cellular functions of lncRNAs.
The “binding sites navigation” module of CLIPdb contains four view tabs. The first tab is a “factor” view. Users select the RBPs of interest and obtain detailed information from the studies displayed in the data matrix. Users are also provided links to data sources and the appropriate literature. The RBP annotation matrices are available only under this factor view tab. The second tab is a “species” view in which users search CLIP-seq data sets in specific species. The third tab is a “cell line” view for searching the RBPs that have been profiled in specific cell lines, such as HeLa and HEK293 cells. The fourth tab is a “technology” view for users to search CLIP-seq studies according to the CLIP technologies used, including HITS-CLIP, iCLIP, PAR-CLIP and gPAR-CLIP.
Downloading and visualizing binding sites
The user may visualize the transcriptome-wide binding sites with different p-value cutoffs (0.01 or 0.001, assigned by Piranha software ) through the “browser” module. The smaller the p-value, the less likely the binding site represents background noise. For example, three binding sites in the UTR region of the Fbxo3 gene are shown using Jbrowse  (Figure 3B). The user may simultaneously select multiple or all CLIP-seq data sets to investigate RBP co-binding patterns in the Fbxo3 gene.
Searching binding sites for a given gene
Users often need to know the RBPs that bound the genes of interest profiled by public CLIP-seq data. Alternatively, users may be interested in knowing whether two proteins frequently compete or co-bind for the same site, in which case each of the proteins may be expected to occupy the binding site under different conditions. The binding target search tool enables users to explore the answers to these questions.
As another example, if a user is interested in the mRNA binding sites for human phosphatase and tensin homolog (PTEN), the user enters or selects “PTEN” in the search box for humans in the “binding target search” module. The database returns 33 RBPs that have binding sites in the PTEN transcript (Figure 3C). The user can obtain detailed information about these binding sites, including the chromosome location, binding strength, p-value, and the CLIP-seq sample and study from which the results were acquired. Importantly, the interactions between RBPs and lncRNAs are not well analyzed in the currently available CLIP-seq studies. With increasing evidence indicating that lncRNAs interact with various RBPs to form complexes , we anticipate that the integrated information provided by CLIPdb may lead to novel hypotheses and interesting discoveries.
As the number of publicly available CLIP-seq data sets continues to increase, we will continue to update the database for newly published CLIP-seq studies. As CLIP-seq assays are applied to more RBPs obtained from a broader set of species, cell lines and tissues, data reuse and integrative analysis may become a greater challenge. Thus, we will provide user-friendly CLIP-seq data analysis pipelines for users, including pre-processing, mapping, peak calling and motif identification. In addition to CLIP-seq data, we plan to include ribosome profiling data in the future. Integrating CLIP-seq data and ribosome profiling data will greatly improve our understanding of post-transcriptional regulatory mechanisms. In addition to CLIP-seq experiments, many other low-resolution data are included in other databases (e.g., NPInter ). Because these resources provide valuable and more complete information on RBP-RNA interactions, they may also be included in future versions of CLIPdb.
CLIPdb is a more complete and better-curated database focused on CLIP-seq data sets than those currently available, providing genome-wide, high-resolution binding sites. Characterizing the post-transcriptional regulatory networks mediated by RNA-binding proteins is a problem of great interest. CLIP-seq technologies (i.e., HITS-CLIP, PAR-CLIP and iCLIP) have proven their power in deciphering complex codes at the single nucleotide level. Here, we presented a useful resource for the reuse and mining of publicly available CLIP-seq data. To our knowledge, CLIPdb is the most comprehensive database available for use in examining published CLIP-seq studies. It provides high-resolution RBP binding sites across whole genomes, including both mRNA (CDSs plus UTRs) and non-coding RNA.
To make the binding sites comparable and the data integrative across multiple samples from different batches, we used uniform pre-processing and peak-calling procedures to identify RBP binding sites from raw CLIP-seq data. CLIPdb embraces all variations of CLIP technologies, including PAR-CLIP with its T-to-C conversions, HITS-CLIP with its deletions at the crosslinked sites, and iCLIP with its truncations at the crosslinked sites. We systematically identified binding sites using Piranha, a tool that has demonstrated success across three CLIP-seq variants . Thus, CLIPdb uniformly scores or ranks binding sites from various CLIP-seq experiments, which is important for the integrative analysis of these published CLIP-seq data. However, Piranha was not optimized for specific CLIP-seq technology, such as PAR-CLIP for which T-to-C conversions indicate binding events. Thus, we also provided binding sites identified by specialized peak-calling tools for specific CLIP-seq technology. In addition, the various sequencing depths of the collected CLIP-seq data sets may significantly affect peak calling. We found very few (<10) binding sites generated from some CLIP-seq data sets because of low sequencing quality and depth (usually less than 1 M unique mapped reads). We do not include these data sets in CLIPdb because the mapping size would have significant impact on the binding sites being called. Users should exercise caution while optimizing the peak-calling procedure when reusing these published CLIP-seq data.
Re-analysis of publicly available CLIP-seq data may generate novel hypotheses. For example, we found that more than 10 orthologous RBPs have been studied using CLIP-seq technologies. Our data will be useful for comparative analyses of their binding sites and may uncover evolutionary signatures for RBPs binding. In addition, CLIPdb provides a unified data set annotation vocabulary as well as several queries and views of the data as a convenience for users. Therefore, CLIPdb should provide additional insights into RNA-protein interaction networks, such as lncRNA-protein interactions and their cellular functions.
Availability and requirements
We thank Chaolin Zhang for help with CIMS usage. We also thank Peipei Yin and Yuchuan Wang for their help in data collection. Finally, we thank the two anonymous reviewers for their extensive comments that have helped us to greatly improve this manuscript and our database. This work is supported by grants from the National Key Basic Research Program (2012CB316503), National High-tech Research and Development Program (2014AA021103) and National Natural Science Foundation of China (31271402, 31100601).
- Keene JD. RNA regulons: coordination of post-transcriptional events. Nat Rev Genet. 2007;8(7):533–43.View ArticlePubMedGoogle Scholar
- Mitchell SF, Parker R. Principles and properties of eukaryotic mRNPs. Mol Cell. 2014;54(4):547–58.View ArticlePubMedGoogle Scholar
- Lunde BM, Moore C, Varani G. RNA-binding proteins: modular design for efficient function. Nat Rev Mol Cell Biol. 2007;8(6):479–90.View ArticlePubMedGoogle Scholar
- Konig J, Zarnack K, Luscombe NM, Ule J. Protein-RNA interactions: new genomic technologies and perspectives. Nat Rev Genet. 2011;13(2):77–83.View ArticleGoogle Scholar
- Licatalosi DD, Mele A, Fak JJ, Ule J, Kayikci M, Chi SW, et al. HITS-CLIP yields genome-wide insights into brain alternative RNA processing. Nature. 2008;456(7221):464–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Hafner M, Landthaler M, Burger L, Khorshid M, Hausser J, Berninger P, et al. Transcriptome-wide identification of RNA-binding protein and microRNA target sites by PAR-CLIP. Cell. 2010;141(1):129–141.View ArticlePubMed CentralPubMedGoogle Scholar
- Konig J, Zarnack K, Rot G, Curk T, Kayikci M, Zupan B, et al. iCLIP reveals the function of hnRNP particles in splicing at individual nucleotide resolution. Nat Struct Mol Biol. 2010;17(7):909–15.View ArticlePubMed CentralPubMedGoogle Scholar
- Ulitsky I, Bartel DP. lincRNAs: genomics, evolution, and mechanisms. Cell. 2013;154(1):26–46.View ArticlePubMed CentralPubMedGoogle Scholar
- Fatica A, Bozzoni I. Long non-coding RNAs: new players in cell differentiation and development. Nat Rev Genet. 2014;15(1):7–21.View ArticlePubMedGoogle Scholar
- Cech TR, Steitz JA. The noncoding RNA revolution-trashing old rules to forge new ones. Cell. 2014;157(1):77–94.View ArticlePubMedGoogle Scholar
- Khorshid M, Rodak C, Zavolan M. CLIPZ: a database and analysis environment for experimentally determined binding sites of RNA-binding proteins. Nucleic Acids Res. 2011;39(Database issue):D245–252.View ArticlePubMed CentralPubMedGoogle Scholar
- Anders G, Mackowiak SD, Jens M, Maaskola J, Kuntzagk A, Rajewsky N, et al. doRiNA: a database of RNA interactions in post-transcriptional regulation. Nucleic Acids Res. 2012;40(Database issue):D180–186.View ArticlePubMed CentralPubMedGoogle Scholar
- Li JH, Liu S, Zhou H, Qu LH, Yang JH. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 2013;42(Database issue):D92–97.PubMed CentralPubMedGoogle Scholar
- Dassi E, Re A, Leo S, Tebaldi T, Pasini L, Peroni D, et al. AURA 2: empowering discovery of post-transcriptional networks. Translation. 2014;2(1):e27738.View ArticleGoogle Scholar
- Yuan J, Wu W, Xie C, Zhao G, Zhao Y, Chen R. NPInter v2.0: an updated database of ncRNA interactions. Nucleic Acids Res. 2013;42(Database issue):D104–108.PubMed CentralPubMedGoogle Scholar
- FASTX-Toolkit package [http://hannonlab.cshl.edu/fastx_toolkit/index.html].
- Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3):R25.View ArticlePubMed CentralPubMedGoogle Scholar
- Uren PJ, Bahrami-Samani E, Burns SC, Qiao M, Karginov FV, Hodges E, et al. Site identification in high-throughput RNA-protein interaction data. Bioinformatics. 2012;28(23):3013–20.View ArticlePubMed CentralPubMedGoogle Scholar
- Corcoran DL, Georgiev S, Mukherjee N, Gottwein E, Skalsky RL, Keene JD, et al. PARalyzer: definition of RNA binding sites from PAR-CLIP short-read sequence data. Genome Biol. 2011;12(8):R79.View ArticlePubMed CentralPubMedGoogle Scholar
- Moore MJ, Zhang C, Gantman EC, Mele A, Darnell JC, Darnell RB. Mapping Argonaute and conventional RNA-binding protein interactions with RNA at single-nucleotide resolution using HITS-CLIP and CIMS analysis. Nat Protoc. 2014;9(2):263–93.View ArticlePubMed CentralPubMedGoogle Scholar
- Zhang C, Darnell RB. Mapping in vivo protein-RNA interactions at single-nucleotide resolution from HITS-CLIP data. Nat Biotechnol. 2011;29(7):607–14.View ArticlePubMed CentralPubMedGoogle Scholar
- Weyn-Vanhentenryck SM, Mele A, Yan Q, Sun S, Farny N, Zhang Z, et al. HITS-CLIP and integrative modeling define the Rbfox splicing-regulatory network linked to brain development and autism. Cell Rep. 2014;6(6):1139–52.View ArticlePubMedGoogle Scholar
- Kozomara A, Griffiths-Jones S. miRBase: annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014;42(Database issue):D68–73.View ArticlePubMed CentralPubMedGoogle Scholar
- Nam JW, Bartel DP DP. Long noncoding RNAs in C. elegans. Genome Res. 2012;22(12):2529–40.View ArticlePubMed CentralPubMedGoogle Scholar
- Coordinators NR. Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2014;42(Database issue):D7–17.View ArticleGoogle Scholar
- Ray D, Kazan H, Cook KB, Weirauch MT, Najafabadi HS, Li X, et al. A compendium of RNA-binding motifs for decoding gene regulation. Nature. 2013;499(7457):172–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Cook KB, Kazan H, Zuberi K, Morris Q, Hughes TR. RBPDB: a database of RNA-binding specificities. Nucleic Acids Res. 2011;39(Database issue):D301–308.View ArticlePubMed CentralPubMedGoogle Scholar
- Baltz AG, Munschauer M, Schwanhausser B, Vasile A, Murakawa Y, Schueler M, et al. The mRNA-bound proteome and its global occupancy profile on protein-coding transcripts. Mol Cell. 2012;46(5):674–90.View ArticlePubMedGoogle Scholar
- Freeberg MA, Han T, Moresco JJ, Kong A, Yang YC, Lu ZJ, et al. Pervasive and dynamic protein binding sites of the mRNA transcriptome in Saccharomyces cerevisiae. Genome Biol. 2013;14(2):R13.View ArticlePubMed CentralPubMedGoogle Scholar
- Skinner ME, Uzilov AV, Stein LD, Mungall CJ, Holmes IH. JBrowse: a next-generation genome browser. Genome Res. 2009;19(9):1630–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Guttman M, Rinn JL. Modular regulatory principles of large non-coding RNAs. Nature. 2012;482(7385):339–46.View ArticlePubMed CentralPubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.