Improvements to the HITS-CLIP protocol eliminate widespread mispriming artifacts
© Gillen et al. 2016
Received: 24 November 2015
Accepted: 28 April 2016
Published: 5 May 2016
High-throughput sequencing of RNA isolated by crosslinking immunoprecipitation (HITS-CLIP) allows for high resolution, genome-wide mapping of RNA-binding proteins. This methodology is frequently used to validate predicted targets of microRNA binding, as well as direct targets of other RNA-binding proteins. Hence, the accuracy and sensitivity of binding site identification is critical.
We found that substantial mispriming during reverse transcription results in the overrepresentation of sequences complementary to the primer used for reverse transcription. Up to 45 % of peaks in publicly available HITS-CLIP libraries are attributable to this mispriming artifact, and the majority of libraries have detectable levels of mispriming. We also found that standard techniques for validating microRNA-target interactions fail to differentiate between artifactual peaks and physiologically relevant peaks.
Here, we present a modification to the HITS-CLIP protocol that effectively eliminates this artifact and improves the sensitivity and complexity of resulting libraries.
KeywordsmicroRNA HITS-CLIP CLIP-seq iCLIP PAR-CLIP
The advent of HITS-CLIP and its derivative techniques has facilitated the mapping of RNA-protein interactions genome-wide in a variety of tissues and organisms [1, 2]. These techniques have made important contributions to the understanding of the function of a range of RNA-binding proteins, from the microRNA-related argonaute (AGO) to other regulatory factors, including FUS and PTBPs [3–5]. However, while useful, these techniques suffer from substantial mispriming during reverse transcription that can contaminate the resulting libraries with false peaks that may be difficult to distinguish from physiologically relevant peaks.
HITS-CLIP library construction  begins with the cross-linking of RNA-binding proteins to RNA using UV light, partial RNA digestion with RNAse A, and immunoprecipitation of the RNA-binding protein of interest and attached RNAs. The protein-bound RNAs are then dephosphorylated and a 3′ RNA adaptor is ligated with T4 RNA ligase. This step is performed “on-bead”, and the efficiency is highly dependent on the structure of the RNA-adaptor interaction . The RNAs are then treated with PNK, again “on-bead”, and the protein-RNA complexes are resolved on an SDS-page gel to select only complexes of the expected size. The resulting size-selected protein-RNA complexes are then digested with Proteinase K to remove the protein and the remaining RNA is purified by phenol extraction and ethanol precipitation. A 5′ RNA linker is then ligated to the purified RNAs, and the RNA is reverse-transcribed using a primer perfectly complementary to the 3′ adaptor.
The use of a perfectly complementary reverse transcription primer results in substantial mispriming to suboptimal matches of the 5′ end of 3′ adaptor sequence in alternative RNAs during reverse transcription. Such mispriming events produce cDNAs from RNAs that may lack 3′ adaptor sequences. These cDNAs, which will be included in final sequencing libraries if they received a 5′ adaptor, contaminate the library with read pileups anchored to genomic occurrences of the adaptor sequence, rather than bona fide sites of protein binding. Recovery of irrelevant misprimed sequences is exacerbated by the fact that proteinase K digestion leaves one or two amino acids crosslinked to the RNA . These crosslinked bases are frequent sites of termination during reverse transcription, meaning that only a fraction of reverse transcription events on crosslinked RNAs will read through to the 5′ adaptor and produce PCR-competent templates. This frequent early termination on crosslinked templates results in a bias towards RNAs that do not have crosslinking sites. This bias leads to overrepresentation of sequences not bound by the RNA-binding protein in the final library. Together, these problems introduce significant artifacts that make downstream analysis challenging.
HITS-CLIP does not provide direct insight into the microRNA (miRNA) involved in a given AGO-mRNA interaction, and the identification of putative miRNA:mRNA targets relies on in silico modeling of miRNA-mRNA interactions. Thus high levels of artifactual peaks would reduce the sensitivity of these computational approaches, possibly identifying false positive mRNA:miRNA pairings. Indeed, we show that widely accepted “wet” validation techniques that rely on overexpression of candidate miRNAs can confirm artifactual miRNA:mRNA interactions gleaned from HITS-CLIP. Here, we summarize our experiences with HITS-CLIP and propose effective experimental changes and bioinformatic steps to improve the construction and analysis of HITS-CLIP libraries. These improvements yield high confidence peaks that provide better insight into the biological system under study.
Contamination of HITS-CLIP libraries with mispriming artifacts due to inefficient ligation
Next, we attempted to estimate the number of artifactual peaks in each sample, based on the over- or underrepresentation of the 3′ adaptor sequence between positions −25 and +75 in each peak (this region, which shows pronounced overrepresentation of the adaptor sequence in affected samples, is highlighted in grey in Fig. 1a). We found that an average of 14.2 % of peaks in the overrepresentation group (Fig. 1a, blue line) are likely artifacts, with a low of 5.5 % and a high of 45.8 % (Fig. 1b, blue line). This is comparable to our samples (Fig. 1b, emerald line), in which 15.3 % of peaks are likely artifacts In contrast, the underrepresentation group (Fig. 1a, vermilion line) is effectively free of the artifact, with the 3′ adaptor sequence observed at a modest 2 % higher than expected (Fig. 1b). Finally, we attempted to see if any obvious factors could account for the variable presence of the artifact in the public datasets. We assessed the contribution of target species (human vs. mouse), target protein (AGO vs. non-AGO), RT conditions and the number of peaks per library, but none were significantly correlated with overrepresentation of the 3′ adaptor sequence (data not shown). We also considered the sequence of the RT primer. Only two reverse transcription primers were used in more than one publication in StarBase: the primer from the 2005 CLIP protocol  (5′-CCGCTGGAAGTGACTGACAC-3′) and the primer from the Illumina Small RNA Kits (5′-CAAGCAGAAGACGGCATACGA-3′). Each primer was used in 17 samples, and no significant difference in artifact representation was observed between the two primers (5.1 % vs. 5.2 % overrepresentation; Additional file 1: Figure S1). The four remaining primers were used by only one group each (1–3 samples each), making it impossible to separate the sequence of the primers from other factors affecting artifact representation. Finally, we considered the impact of the two major modifications to the HITS-CLIP protocol (iCLIP and PAR-CLIP) on mispriming during RT. We analyzed the 10 iCLIP datasets in StarBase 2.0 (10 samples from two research groups), and found only slight overrepresentation (25 %) of the complement of the reverse transcription primer on average (Additional file 2: Figure S2). While the artifact is much better controlled in the iCLIP libraries than in HITS-CLIP libraries, it is still present. PAR-CLIP, which uses photoreactive ribonucleoside analogs to map crosslinking sites, allows legitimately crosslinked sequences to be identified based on the presence of mutated bases. While mispriming may still occur in these datasets, non-crosslinked sequences (including products of mispriming) can be effectively filtered bioinformatically.
A nested reverse transcription primer, combined with a protected reverse PCR primer, eliminates mispriming artifacts
Removal of artifactual reads eliminates false peaks and improves sensitivity, but does not affect physiologically relevant peaks
We found that the majority of published HITS-CLIP libraries either contain substantial mispriming artifacts or appear to have been processed to remove these artifacts, and that these artifacts can lead to false positive peaks that pass standard validation techniques. The mispriming in the most severely affected libraries may simply be a symptom of generally poor library quality, but without proper controls (libraries prepared in parallel with and without our modifications), it is impossible to differentiate between mispriming resulting from poor library quality and mispriming causing poor library quality. Whatever the cause, this issue is underappreciated, as none of the publications represented in StarBase make any mention of this pervasive artifact, or of any explicit steps taken to mitigate the artifact. It is important to recognize and address this problem, as HITS-CLIP is routinely used to provide evidence of miRNA-mRNA interactions under physiological conditions, in addition to many other RNA-protein interactions. Without such evidence, it is difficult to assess the validity of reported miRNA-mRNA interactions, as most commonly accepted validation techniques (reporter assays, target transcript and protein quantification, etc.) merely provide evidence for binding capacity rather than biologically relevant interaction. As we have demonstrated, validation techniques can yield false positive results when supported by artifactual HITS-CLIP data. We show that an artifactual HITS-CLIP peak (Fig. 4) leads to “validation” of a robust interaction between (overexpressed) miR-888 and endogenous Progesterone Receptor in MCF-7 cells (Fig. 5), including repression of downstream targets of PGR (CK5; Fig. 5c). Overexpressing miR-888 represses PGR expression, but this is not surprising, as the miRNA is highly complementary to the target. While miR-888 is clearly capable of regulating PGR expression, and may do so in other cell types, there is no evidence that it does under normal conditions in MCF-7, where it is not expressed at detectable levels. Here, we present a modification that mitigates the problem while making minimal changes to the original HITS-CLIP protocol. In the case of AGO HITS-CLIP, our modifications to the method and downstream analysis result in high confidence measurements of the miRNA-mRNA interactome under physiologic conditions. These measurements provide direct support for miRNA-target interactions inferred by other methods, and provide biological relevance to downstream functional assays.
There are two existing variants of HITS-CLIP that partially address the problems described here: PAR-CLIP and iCLIP [30, 36]. PAR-CLIP (Photoactivatable-Ribonucleoside-Enhanced Crosslinking and Immunoprecipitation) differs from HITS-CLIP in that it uses the incorporation of photoreactive ribonucleoside analogs, such as 4-thiouridine (4-SU) and 6-thioguanosine (6-SG), to map sites of crosslinking. Cells are grown in the presence of these analogs, which are incorporated into nascent RNA transcripts. The use of these ribonucleoside analogs results in both improved RNA recovery and the introduction of a marker of crosslinking, in the form of mutations at the site of crosslinking, that is detectable in sequencing. While PAR-CLIP is subject to the same potential for mispriming as HITS-CLIP, excluding reads without mutations introduced by crosslinking can be used to effectively retain only legitimately crosslinked sequences. This will eliminate the majority of misprimed sequences, which are not crosslinked. However, while non-crosslinked sequences can be removed from PAR-CLIP libraries, the presence of substantial mispriming artifacts will still reduce effective library size and library quality. Thus, the same nested RT primer approach suggested here for HITS-CLIP can be applied to PAR-CLIP, with similar benefits.
iCLIP (individual-nucleotide resolution Cross-Linking and Immunoprecipitation) is much more effective than either unmodified HITS-CLIP or PAR-CLIP at mitigating the artifact described here because it addresses the issue of reverse transcription termination at sites of crosslinking. This is accomplished by reverse transcribing with a primer complementary to the 3′ adaptor that also includes the 5′ adaptor sequence. The resulting cDNAs are then circularized and cut at a site within the RT primer to produce a linear cDNA with 5′ and 3′ adaptor sequences. Unlike HITS-CLIP and PAR-CLIP, the resulting cDNAs are PCR-competent regardless of where reverse transcription terminated, removing the bias introduced by requiring the reverse transcriptase to read through the crosslinking site to the ligated 5′ adaptor. This increased efficiency, along with the use of a nested RT primer, likely explains the lower rate of mispriming artifacts observed in iCLIP libraries (Additional file 2: Figure S2). Termination at the site of crosslinking also allows for the identification of crosslinking sites indirectly through the mapping of 5′ ends of cDNAs (reverse transcription termination sites). Despite the fact that the mispriming artifact is present at low levels in the iCLIP libraries we analyzed (Additional file 2: Figure S2), there has been at least one attempt to mitigate the inefficiency of RNA ligation in iCLIP, though no explicit rationale is given for the change in protocol . In this case, the authors used a biotinylated 3′ adaptor, which allowed for enrichment of RNA-protein complexes with ligated adaptors using an avidin pull-down. While this should improve the efficiency of reverse transcription on legitimate templates by enriching for adaptor-ligated RNAs, this represents a substantial alteration to the iCLIP protocol, and has the undesirable effect of adding another pull-down to an already complex protocol while still not taking full advantage of the nested iCLIP reverse transcription primer. A much more satisfying solution, either alone or in combination with the biotinlyated 3′ adaptor, is to use a 3′ protected reverse PCR primer, as we have described here for HITS-CLIP. No matter which approach is taken to improve HITS-CLIP, it is important to recognize the limitations of the technique and process the data appropriately.
Here, we describe the use of a nested reverse transcription primer and a protected reverse PCR primer in place of the original oligonucleotide primers to eliminate widespread mispriming artifacts in HITS-CLIP. These simple but effective modifications help to improve HITS-CLIP library quality by producing high-confidence peak calls that accurately represent the RNA binding protein-RNA interactome under physiological conditions. Specifically, removal of artifactual reads both eliminates false peaks (while not affecting physiologically relevant peaks) and improves sensitivity, allowing for the identification of weaker interactions. These modifications have the potential to improve library quality not just for AGO HITS-CLIP, as described here, but for all HITS-CLIP targets and all library construction techniques which rely on the ligation of a 3′ adaptor to RNA.
For each set of peak calls generated in our lab or obtained from StarBase v2.0 , we expanded the peaks to a uniform 200 bp (around the existing center of the peak). We then counted the number of occurrences of the first six bases of the 3′ adaptor (allowing for one mismatch) in the uniform peaks using 20 bp sliding windows (with a 6 bp shift between each window), and plotted this data as the fold-difference relative to the expected frequency of each adaptor-complement (using 1 million randomly sampled genomic sequences of 200 bp). We confirmed the 3′ adapter sequence in each library by downloading the raw data from NCBI SRA (www.ncbi.nlm.nih.gov/sra) and examining the reads directly and with FASTQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/).
HITS-CLIP was performed as in [3 N], with the following modifications: A reverse transcription primer with the sequence 5′-CCGCTGGAAGTGACTGA-3′ was used in place of the “DP3” primer in the original method. Following reverse transcription, the reaction was treated with Exonuclease I (NEB) to remove unincorporated reverse transcription primer. Specifically, 10 U of ExoI was added to each sample and they were incubated at 37C for 30 min followed by heat inactivation of ExoI at 80C for 20 min. A 5′ extension was added to the reverse PCR primer to allow for multiplexing and use on the Illumina Hiseq platform, and phosphorothioate bonds were added between the final four bases to prevent cleavage by exonuclease activity of the polymerase (5′-GTGACTGGAGTTCAGACGTGTGCTCTTCCGATCCCGCTGGAAGTGACTGA*C*A*C-3′, “*” indicate phosphorothioate bonds). The second PCR was performed with index primers compatible with the new first-stage primers (5′-CAAGCAGAAGACGGCATACGAGATCGGT NNNNNN GTGACTGGAGTTCAGACGTGTGCTCTTCCGATC-3′).
Reads were processed with cutadapt (https://code.google.com/p/cutadapt/)  to remove 3′ adaptor sequences and 3′ bases with QUAL < 10. Trimmed reads were aligned to the genome using GSNAP (http://research-pub.gene.com/gmap/) [39, 40], and peaks were called on alignments with MAPQ >10 using piranha (http://smithlabresearch.org/software/piranha/). Overlapping peaks were then merged and all peaks were resized to a uniform 60 bp centered on the summit of each peak using peaktools (https://github.com/jayhesselberth/peaktools). To normalize library sizes for peak calling sensitivity assessment, 4 million aligned reads were randomly sampled without replacement from each aligned library by shuffling the alignment order with GNU shuf (seeded with/dev/random) and retaining the first 4 million alignments.
Cell culture, transfection and treatment assays
The human breast cancer cell lines T47D and MCF7 were originally obtained from Iafa Keydar , and Sam Brooks, the Michigan Cancer Foundation, respectively. Cells were maintained in MEM supplemented with 5 % fetal bovine serum (FBS) and penicillin/streptomycin. BT474 cells were obtained from ATCC (Manassas, VA) and maintained in DMEM medium supplemented with 5 % FBS and penicillin/streptomycin. HEK-293 T cell line was obtained from ATCC (Manassas, VA) and maintained in DMEM with 10 % FBS. Cell line authenticity was confirmed by short tandem repeat analyses in the University of Colorado DNA Sequencing Core Laboratory .
Plasmids, viral constructs and DNA cloning
Luciferase reporter vectors were created by cloning of the miRNA Responsive Element (MRE) sequences into the XhoI-NotI restriction site of psi-Check2 vector (Promega, Madison, WI) as previously reported . Retroviral constructs for stable overexpression of miR-888 were constructed as previously reported .
HEK-293 T cells were plated at approximately 40 % confluence (50,000 cells per well) in 24 well plates. Sixteen hours later the psi-Check2 vectors and miR mimic/control mimic were co-transfected using Lipofectamine 2000 as per manufacturer’s protocol. Cells were harvested by passive lysis and assayed for firefly and renilla luciferase activity using Dual Luciferase Reporter Assay kit on a Glomax luminometer (Promega) following manufacturer’s instructions. The student’s t-test was used to determine significance.
Cell lysates were quantitated for protein content, and approximately 25 mg lysates were resolved by standard SDS-PAGE and transferred to PVDF membrane. Blots were probed by antibodies against PGR (Cell Signaling Technology, Danvers, MA) and an appropriate HRP-conjugated secondary antibody. Blots were then stripped and reprobed for alpha tubulin (Sigma-Aldrich, St. Louis, MO) to ensure equal protein loading. The blots were then visualized using enhanced chemoluminescence.
Cells were fixed onto glass coverslips in ice-cold 70 % acetone/30 % methanol for 5 min and stained by immunocytochemistry (ICC) with antibodies to CK5 (Dako, Glostrup, Denmark). Secondary antibodies were anti-mouse and anti-rabbit Alexa Fluors 488 (green). Nuclei of cells were counterstained with DAPI. A Nikon E600 microscope (Nikon Corporation, Tokyo, Japan) was used for photography. Images were shot in black and white using ImagePro software (Media Cybernetics, Rockville, MD) and merged and pseudo-colored in Adobe Photoshop CS3 (Adobe Systems, San Jose, CA).
Availability of supporting data
All publicly available HITS-CLIP and iCLIP data supporting the results of this article are available at StarBase 2.0 (http://starbase.sysu.edu.cn/).
The original HITS-CLIP datasets supporting the conclusions of this article are available in NCBI Sequence Read Archive (SRA) studies SRP042140 (http://trace.ncbi.nlm.nih.gov/Traces/sra/?study=SRP042140) and SRP070520 (http://trace.ncbi.nlm.nih.gov/Traces/sra/?study=SRP070520).
This work was supported by the American Cancer Society (to JH), Grohne Cancer Research Fund (to PK) and National Institutes of Health (NCI CA164048 to PK). Funding for open access charge: National Institutes of Health.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Ule J, Jensen K, Mele A, Darnell RB. CLIP: a method for identifying protein-RNA interaction sites in living cells. Methods. 2005;37(4):376–86. doi:10.1016/j.ymeth.2005.07.018.View ArticlePubMedGoogle Scholar
- Ule J, Jensen KB, Ruggiu M, Mele A, Ule A, Darnell RB. CLIP identifies Nova-regulated RNA networks in the brain. Science. 2003;302(5648):1212–5. doi:10.1126/science.1090095.View ArticlePubMedGoogle Scholar
- Chi SW, Zang JB, Mele A, Darnell RB. Argonaute HITS-CLIP decodes microRNA-mRNA interaction maps. Nature. 2009;460(7254):479–86. doi:10.1038/nature08170.PubMedPubMed CentralGoogle Scholar
- Xue Y, Ouyang K, Huang J, Zhou Y, Ouyang H, Li H, et al. Direct conversion of fibroblasts to neurons by reprogramming PTB-regulated microRNA circuits. Cell. 2013;152(1-2):82–96. doi:10.1016/j.cell.2012.11.045.View ArticlePubMedPubMed CentralGoogle Scholar
- Nakaya T, Alexiou P, Maragkakis M, Chang A, Mourelatos Z. FUS regulates genes coding for RNA-binding proteins in neurons by binding to their highly conserved introns. RNA. 2013;19(4):498–509. doi:10.1261/rna.037804.112.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhuang F, Fuchs RT, Sun Z, Zheng Y, Robb GB. Structural bias in T4 RNA ligase-mediated 3′-adapter ligation. Nucleic Acids Res. 2012;40(7):e54. doi:10.1093/nar/gkr1263.View ArticlePubMedPubMed CentralGoogle Scholar
- Urlaub H, Hartmuth K, Luhrmann R. A two-tracked approach to analyze RNA-protein crosslinking sites in native, nonlabeled small nuclear ribonucleoprotein particles. Methods. 2002;26(2):170–81. doi:10.1016/S1046-2023(02)00020-8.View ArticlePubMedGoogle 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. 2014;42(Database issue):D92–7. doi:10.1093/nar/gkt1248.View ArticlePubMedPubMed CentralGoogle Scholar
- Keydar I, Chen L, Karby S, Weiss FR, Delarea J, Radu M, et al. Establishment and characterization of a cell line of human breast carcinoma origin. Eur J Cancer. 1979;15(5):659–70.View ArticlePubMedGoogle Scholar
- Chen CZ, Li L, Lodish HF, Bartel DP. MicroRNAs modulate hematopoietic lineage differentiation. Science. 2004;303(5654):83–6. doi:10.1126/science.1091903.View ArticlePubMedGoogle Scholar
- Sanford JR, Wang X, Mort M, Vanduyn N, Cooper DN, Mooney SD, et al. Splicing factor SFRS1 recognizes a functionally diverse landscape of RNA transcripts. Genome Res. 2009;19(3):381–94. doi:10.1101/gr.082503.108.View ArticlePubMedPubMed CentralGoogle Scholar
- Kishore S, Jaskiewicz L, Burger L, Hausser J, Khorshid M, Zavolan M. A quantitative analysis of CLIP methods for identifying binding sites of RNA-binding proteins. Nat Methods. 2011;8(7):559–64. doi:10.1038/nmeth.1608.View ArticlePubMedGoogle Scholar
- Leung AK, Young AG, Bhutkar A, Zheng GX, Bosson AD, Nielsen CB, et al. Genome-wide identification of Ago2 binding sites from mouse embryonic stem cells with and without mature microRNAs. Nat Struct Mol Biol. 2011;18(2):237–44. doi:10.1038/nsmb.1991.View ArticlePubMedPubMed CentralGoogle Scholar
- Polymenidou M, Lagier-Tourenne C, Hutt KR, Huelga SC, Moran J, Liang TY, et al. Long pre-mRNA depletion and RNA missplicing contribute to neuronal vulnerability from loss of TDP-43. Nat Neurosci. 2011;14(4):459–68. doi:10.1038/nn.2779.View ArticlePubMedPubMed CentralGoogle 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. doi:10.1038/nbt.1873.View ArticlePubMedPubMed CentralGoogle Scholar
- Charizanis K, Lee KY, Batra R, Goodwin M, Zhang C, Yuan Y, et al. Muscleblind-like 2-mediated alternative splicing in the developing brain and dysregulation in myotonic dystrophy. Neuron. 2012;75(3):437–50. doi:10.1016/j.neuron.2012.05.029.View ArticlePubMedPubMed CentralGoogle Scholar
- Haecker I, Gay LA, Yang Y, Hu J, Morse AM, McIntyre LM, et al. Ago HITS-CLIP expands understanding of Kaposi's sarcoma-associated herpesvirus miRNA function in primary effusion lymphomas. PLoS Pathog. 2012;8(8):e1002884. doi:10.1371/journal.ppat.1002884.View ArticlePubMedPubMed CentralGoogle Scholar
- Ince-Dunn G, Okano HJ, Jensen KB, Park WY, Zhong R, Ule J, et al. Neuronal Elav-like (Hu) proteins regulate RNA splicing and abundance to control glutamate levels and neuronal excitability. Neuron. 2012;75(6):1067–80. doi:10.1016/j.neuron.2012.07.009.View ArticlePubMedPubMed CentralGoogle Scholar
- Ishigaki S, Masuda A, Fujioka Y, Iguchi Y, Katsuno M, Shibata A, et al. Position-dependent FUS-RNA interactions regulate alternative splicing events and transcriptions. Sci Rep. 2012;2:529. doi:10.1038/srep00529.View ArticlePubMedPubMed CentralGoogle Scholar
- Korch C, Spillman MA, Jackson TA, Jacobsen BM, Murphy SK, Lessey BA, et al. DNA profiling analysis of endometrial and ovarian cell lines reveals misidentification, redundancy and contamination. Gynecol Oncol. 2012;127(1):241–8. doi:10.1016/j.ygyno.2012.06.017.View ArticlePubMedPubMed CentralGoogle Scholar
- Loeb GB, Khan AA, Canner D, Hiatt JB, Shendure J, Darnell RB, et al. Transcriptome-wide miR-155 binding map reveals widespread noncanonical microRNA targeting. Mol Cell. 2012;48(5):760–70. doi:10.1016/j.molcel.2012.10.002.View ArticlePubMedPubMed CentralGoogle Scholar
- Licatalosi DD, Yano M, Fak JJ, Mele A, Grabinski SE, Zhang C, et al. Ptbp2 represses adult-specific splicing to regulate the generation of neuronal precursors in the embryonic brain. Genes Dev. 2012;26(14):1626–42. doi:10.1101/gad.191338.112.View ArticlePubMedPubMed CentralGoogle Scholar
- Macias S, Plass M, Stajuda A, Michlewski G, Eyras E, Caceres JF. DGCR8 HITS-CLIP reveals novel functions for the Microprocessor. Nat Struct Mol Biol. 2012;19(8):760–6. doi:10.1038/nsmb.2344.View ArticlePubMedPubMed CentralGoogle Scholar
- Riley KJ, Rabinowitz GS, Yario TA, Luna JM, Darnell RB, Steitz JA. EBV and human microRNAs co-target oncogenic and apoptotic viral and human genes during latency. EMBO J. 2012;31(9):2207–21. doi:10.1038/emboj.2012.63.View ArticlePubMedPubMed CentralGoogle Scholar
- Sauliere J, Murigneux V, Wang Z, Marquenet E, Barbosa I, Le Tonqueze O, et al. CLIP-seq of eIF4AIII reveals transcriptome-wide mapping of the human exon junction complex. Nat Struct Mol Biol. 2012;19(11):1124–31. doi:10.1038/nsmb.2420.View ArticlePubMedGoogle Scholar
- Wilbert ML, Huelga SC, Kapeli K, Stark TJ, Liang TY, Chen SX, et al. LIN28 binds messenger RNAs at GGAGA motifs and regulates splicing factor abundance. Mol Cell. 2012;48(2):195–206. doi:10.1016/j.molcel.2012.08.004.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang ET, Cody NA, Jog S, Biancolella M, Wang TT, Treacy DJ, et al. Transcriptome-wide regulation of pre-mRNA splicing and mRNA localization by muscleblind proteins. Cell. 2012;150(4):710–24. doi:10.1016/j.cell.2012.06.041.View ArticlePubMedPubMed CentralGoogle Scholar
- Karginov FV, Hannon GJ. Remodeling of Ago2-mRNA interactions upon cellular stress reflects miRNA complementarity and correlates with altered translation rates. Genes Dev. 2013;27(14):1624–32. doi:10.1101/gad.215939.113.View ArticlePubMedPubMed CentralGoogle Scholar
- Pandit S, Zhou Y, Shiue L, Coutinho-Mansfield G, Li H, Qiu J, et al. Genome-wide analysis reveals SR protein cooperation and competition in regulated splicing. Mol Cell. 2013;50(2):223–35. doi:10.1016/j.molcel.2013.03.001.View ArticlePubMedPubMed CentralGoogle Scholar
- Konig J, Zarnack K, Rot G, Curk T, Kayikci M, Zupan B, et al. iCLIP--transcriptome-wide mapping of protein-RNA interactions with individual nucleotide resolution. J Vis Exp. 2011;50. doi:10.3791/2638.Google Scholar
- Pillai MM, Gillen AE, Yamamoto TM, Kline E, Brown J, Flory K, et al. HITS-CLIP reveals key regulators of nuclear receptor signaling in breast cancer. Breast Cancer Res Treat. 2014;146(1):85–97. doi:10.1007/s10549-014-3004-9.View ArticlePubMedPubMed CentralGoogle Scholar
- Leucci E, Cocco M, Onnis A, De Falco G, van Cleef P, Bellan C, et al. MYC translocation-negative classical Burkitt lymphoma cases: an alternative pathogenetic mechanism involving miRNA deregulation. J Pathol. 2008;216(4):440–50. doi:10.1002/path.2410.View ArticlePubMedGoogle Scholar
- Yamakuchi M, Ferlito M, Lowenstein CJ. miR-34a repression of SIRT1 regulates apoptosis. Proc Natl Acad Sci U S A. 2008;105(36):13421–6. doi:10.1073/pnas.0801613105.View ArticlePubMedPubMed CentralGoogle Scholar
- Kong YW, Cannell IG, de Moor CH, Hill K, Garside PG, Hamilton TL, et al. The mechanism of micro-RNA-mediated translation repression is determined by the promoter of the target gene. Proc Natl Acad Sci U S A. 2008;105(26):8866–71. doi:10.1073/pnas.0800650105.View ArticlePubMedPubMed CentralGoogle Scholar
- Cannell IG, Kong YW, Johnston SJ, Chen ML, Collins HM, Dobbyn HC, et al. p38 MAPK/MK2-mediated induction of miR-34c following DNA damage prevents Myc-dependent DNA replication. Proc Natl Acad Sci U S A. 2010;107(12):5375–80. doi:10.1073/pnas.0910015107.View ArticlePubMedPubMed CentralGoogle Scholar
- Hafner M, Landthaler M, Burger L, Khorshid M, Hausser J, Berninger P, et al. PAR-CliP--a method to identify transcriptome-wide the binding sites of RNA binding proteins. J Vis Exp. 2010;41. doi:10.3791/2034.Google Scholar
- Chen T, Xiang JF, Zhu S, Chen S, Yin QF, Zhang XO, et al. ADAR1 is required for differentiation and neural induction by regulating microRNA processing in a catalytically independent manner. Cell Res. 2015;25(4):459–76. doi:10.1038/cr.2015.24.View ArticlePubMedPubMed CentralGoogle Scholar
- Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnetjournal. 2011;17(1):10–2.Google Scholar
- Wu TD, Watanabe CK. GMAP: a genomic mapping and alignment program for mRNA and EST sequences. Bioinformatics. 2005;21(9):1859–75. doi:10.1093/bioinformatics/bti310.View ArticlePubMedGoogle Scholar
- Wu TD, Nacu S. Fast and SNP-tolerant detection of complex variants and splicing in short reads. Bioinformatics. 2010;26(7):873–81. doi:10.1093/bioinformatics/btq057.View ArticlePubMedPubMed CentralGoogle Scholar
- Pillai MM, Yang X, Balakrishnan I, Bemis L, Torok-Storb B. MiR-886-3p down regulates CXCL12 (SDF1) expression in human marrow stromal cells. PLoS One. 2010;5(12):e14304. doi:10.1371/journal.pone.0014304.View ArticlePubMedPubMed CentralGoogle Scholar
- Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, et al. The human genome browser at UCSC. Genome Res. 2002;12(6):996–1006. doi:10.1101/gr.229102. Article published online before print in May 2002.View ArticlePubMedPubMed CentralGoogle Scholar
- Meyer LR, Zweig AS, Hinrichs AS, Karolchik D, Kuhn RM, Wong M, et al. The UCSC Genome Browser database: extensions and updates 2013. Nucleic Acids Res. 2013;41(Database issue):D64–9. doi:10.1093/nar/gks1048.View ArticlePubMedPubMed CentralGoogle Scholar