- Open Access
Potential roles of microRNAs in regulating long intergenic noncoding RNAs
© Juan et al.; licensee BioMed Central Ltd. 2013
- Published: 23 January 2013
Over 10,000 long intergenic non-coding RNAs (lincRNAs) have been identified in the human genome. Some have been well characterized and known to participate in various stages of gene regulation. In the post-transcriptional process, another class of well-known small non-coding RNA, or microRNA (miRNA), is very active in inhibiting mRNA. Though similar features between mRNA and lincRNA have been revealed in several recent studies, and a few isolated miRNA-lincRNA relationships have been observed. Despite these advances, the comprehensive miRNA regulation pattern of lincRNA has not been clarified.
In this study, we investigated the possible interaction between the two classes of non-coding RNAs. Instead of using the existing long non-coding database, we employed an ab initio method to annotate lincRNAs expressed in a group of normal breast tissues and breast tumors.
Approximately 90 lincRNAs show strong reverse expression correlation with miRNAs, which have at least one predicted target site presented. These target sites are statistically more conserved than their neighboring genetic regions and other predicted target sites. Several miRNAs that target to these lincRNAs are known to play an essential role in breast cancer.
Similar to inhibiting mRNAs, miRNAs show potential in promoting the degeneration of lincRNAs. Breast-cancer-related miRNAs may influence their target lincRNAs resulting in differential expression in normal and malignant breast tissues. This implies the miRNA regulation of lincRNAs may be involved in the regulatory process in tumor cells.
- miRNA Target Site
- Generalize Linear Model Model
- Putative Exon
- Predict Target Site
- Intergenic Noncoding RNAs
Deep sequencing data from the Encyclopedia of DNA Elements Consortium (ENCODE) suggests that 70-80% of the human genome can be transcribed, and non-protein-coding RNAs (ncRNA) exceed the number of protein-coding genes . The recent discovery of a large number of non-coding RNAs (ncRNAs) significantly enriches the portfolio of potential genetic factors. Rather than being transcriptional noise, many ncRNAs serve as master regulators that affect expression levels of dozens or even hundreds of target genes [2–7]. These regulatory RNAs integrate signals from both genetic and environmental factors, and therefore can play major roles in controlling biological processes . Most notably, a strong association of epigenetic marks with long intergenic non-coding RNAs (lincRNAs >200 nucleotides) in humans and mice was recently described [9–11]. The lincRNAs show evolutionary conservation and spatiotemporally restricted expression patterns, implying that they are functional and regulated. These lincRNAs are reported to regulate dosage compensation, imprinting, and development by establishing chromatin domains in an allele- and cell-type specific manner [12, 13]. According to an analysis on half-lives, lincRNAs are more stable than intron-derived long noncoding RNAs, as they are usually spliced, which also suggest widespread functionality .
In two recent studies, over 10,000 lincRNA regions were identified in human and mouse genomes. These regions were discovered based on epigenetic marks, where the promoter mark H3K4me3 is followed by the transcript mark H3K36me3. Comparative genomics analysis suggests that a significant portion of these regions do not have strong protein coding potentials, and do not have sequence homology with other known proteins . The K4-K36 signature indicates that the transcription of these non-coding RNAs is regulated in a similar fashion as protein-coding genes. Although the molecular mechanisms through which most lincRNAs function were unknown, diverse regulatory mechanisms have been reported for several well-characterized lincRNAs [15, 16]. By targeting the chromatin modification complexes and RNA-binding proteins (XIST, AIR) [17, 18], they can inhibit gene expression (HOTAIR, DBE-T) [12, 19], and control alternative splicing (MALAT-1) . In addition, some lincRNAs are found to be differentially expressed (HOX) or coactivate other proteins (SRA) in tumor cells [10, 21], and are reported to be strongly associated with tumorigenesis [22, 23].
MicroRNA (miRNA) is another type of noncoding RNA whose biological functions have been extensively studied. MiRNAs are ~22nt small non-coding RNAs  that promote mRNA degeneration and/or inhibit their translation by complementarily binding on the 3' un-translated region (3'UTR) of mRNAs  and attracting RNA-induced silencing complex (RISC). In the past decades, over one thousand miRNAs have been identified [25, 26]. They are broadly related with regular cell functions and diseases [27–30].
Despite the extensive investigation on the roles of miRNAs in regulating protein-coding genes, only a few isolated non-coding transcripts (CDR1 antisense and MEG3) are reported to be regulated by miRNA [31, 32]. Though CDR1 antisense is not an intergenic ncRNA, and MEG3 is indirectly regulated by miR-29, these discoveries encouraged us to explore the potential roles of miRNA in regulating non-specific lincRNA expression. Because the primary mechanism of miRNA regulation is by targeting RNA for degradation , they are likely to also broadly regulate lincRNA expression. In fact, 5% of genomic regions interacting with argonaute proteins are located in non-coding RNAs . The argonaute protein family is one of the major components of the RNA-induced silencing complex (RISC) that is responsible for gene silencing due to miRNA expression and RNA interference . In this study, we intend to generally investigate the potential roles of miRNAs in regulating the lincRNA expression.
GENCODE project and NONCODE database have both collected and integrated long noncoding RNA annotations[36, 37]. Based on GENCODE annotation, miRcode provides a map of possible miRNA targets on long noncoding RNAs by using multiple bioinformatics prediction software . However, in tissue-specific lincRNA studies, the ab initio approach is usually employed to reconstruct the annotation of lincRNAs from the RNA-seq data [2, 3, 39]. Considering that lincRNA shows a high tissue or cell-specific expression pattern , this approach is also appropriate for our study.
RNA-seq experiments sequence millions to billions of short RNA fragments in a single experiment, by which it measures the expression levels of the whole transcriptome, including noncoding RNAs, without relying on detailed annotation. To explore the potential miRNA-lincRNA expression relationship, we analyzed RNA-seq data from a group of normal and malignant breast tissues. Random primers were used in the reverse transcription step, and the strands of the transcripts were preserved to aid in identifying non-coding RNAs that are not poly-adenylated. Due to their small size (~22nt), the mature miRNAs were not included in the sequencing libraries. However, our data suggested that we can detect precursor miRNAs (80-120nt) from these RNA-seq data sets. For lincRNAs, the sequence fragments across two exons of the transcripts (junction reads) are able to reveal exon boundaries, especially novel ones. This feature helps us to understand the splicing patterns of unannotated lincRNAs. Several tools have been developed for mapping junction reads , detecting the splicing events, assembling isoforms and estimating expression [41, 42].
In this study, for each putative miRNA-lincRNA pair, (1) we required reverse correlation between the expression levels of pre-miRNA and lincRNA within the 20 examined tissue samples. This is due to the well-established roles of miRNAs in facilitating RNA degradation on protein coding genes; and (2) we also required seed sequence of the miRNA (positions 2-7 on the miRNA) should be present in the exon regions of identified lincRNA, and the surrounding sequences on the lincRNA should also favor miRNA targeting conditions. We further examined the evolutionary conservation of the identified miRNA seeds in lincRNAs.
Based on the RNA-seq data, the discovery of lincRNAs in breast tissue includes the following steps: (1) Epigenetic marks and comparative genetics. As reported previously, thousands of lincRNAs have been discovered by searching for the genomic regions with enhanced promoter marks (H3K4me3, or K4) following the transcript marks (H3K36me3, or K36). Our analysis started from 2,513 intergenic K4-K36 regions in the human genome (hg18) and 1,665 intergenic K4-K36 regions in the mouse genome (mm9) from previous studies [2, 3]. We further mapped the mouse intergenic K4-K36 regions to 1,502 genomic loci in human genome using the liftOver tool in the UCSC Genome Browser . LiftOver can convert genomic coordinates between selected assemblies of the same or different species. After merging the two sources of K4-K36 regions and removing the ones near or overlapping with known genes, 3,096 candidate lincRNA regions remained.
(2) Evidence of transcriptional activity in other breast cancer cell lines. We have previously reported genome-wide binding patterns of RNA polymerase II (RPolII) in MCF7 cell lines , derived from chromatin immunoprecipitation following high throughput sequencing technology (ChIP-seq). Although not directly measure the transcriptional activity for the tissues being measured, the RPolII signals in the MCF7 provide fair estimation in breast cells. We infer the lincRNA transcriptional activity by counting the number of ChIP-seq-derived RPolII tags falling into the lincRNA genomic regions; these counts were further normalized by the length of the lincRNA and the total number of reads for the sample (similar to the RPKM measure in RNA-seq studies [46–48]). Consistent to the results from previous studies , the lincRNAs show intermediate transcriptional activities between protein coding genes and random genomic regions (Figure 2a). We focused our further analysis on the 2,120 lincRNAs (out of 3,096) whose RPolII signal intensities were more than 0.1 RPKM.
(3) Coding potential. We evaluated coding potentials of the identified lincRNA regions using the methods presented in [49, 50]. By modeling the mammalian Codon Substitution Frequency (CSF) of transcript regions and random genomic regions, a CSF score was calculated for each region to represent the codon substitution pattern of the region. A higher CSF score indicates higher potential for being a protein coding gene. The coding potentials of the lincRNA regions are slightly but significantly higher than the ones in random genomic regions (p < 5.4e-07, Wilcoxon test on CSF scores), while the coding potential of protein-coding genes is much higher than lincRNA regions (Figure 2b). For further analysis, we excluded 47 regions with high CSF scores >20; these regions may represent protein-coding genes that are not included in the current annotation database.
We used Scripture  for reconstructing lincRNA exon structures from RNA-seq data. Scripture scans read-enriched regions as putative exons and finds exon boundaries supported by the reads across potential junctions. Scripture identified 525 lincRNAs in the 2,073 candidate regions; this percentage is comparable with the lincRNAs identified by previous studies in other tissues. The genomic coordinates of these lincRNAs are listed in Additional file 1. The average lincRNA length is 1201.7 bases, 75.2% of which are shorter than 1000 bases. The longest lincRNA has 32,178 bases. The lincRNAs are composed by 7.12 exons on average. The mean exon length is 168.8 bases. Almost half (46.7%) of the exons are shorter than 100 bases, with the longest has 5,242 bases. The annotated lincRNAs show a very similar expression patterns comparing with exons and introns with protein-coding genes (Figure 2c). Although lincRNAs are less conserved comparing to protein coding genes, their exons are significantly more conserved than introns (p < 0.0035) and random genome regions (p < 1.35e-14) (Figure 2d). These results provide strong evidence that the lincRNA is functional.
We employed TargetScan  and PITA  to predict putative miRNA target sites on the exonic regions of lincRNAs by 7-mer and 8-mer seed matching (default parameters). We identified 44,887 putative binding sites that belong to 39,384 pairs of miRNA and lincRNA. All 525 lincRNAs and 677 miRNAs are involved in these predicted results. The same seed sequence could be shared in a miRNA family. A miRNA can also bind to multiple sites of a single lincRNA. We also downloaded conserved target prediction of miRNAs and genes for further analysis from the website of TargetScan, including 110,284 predicted pairs of 9,448 genes and 249 conserved human miRNAs.
Expression reverse correlation
We quantified the expression levels of mRNA, lincRNA and precursor miRNAs using RPKM measures (reads per kilo-base exon model per million mappable reads). In total, there are 15,381 genes, 303 lincRNAs and 286 miRNAs expressed more than 0.5 RPKM in more than 15 samples. Our further analysis focuses on these sets of genes whose expression levels are detectable.
A generalized linear model (GLM) was used for modeling the potential effects of miRNA in down-regulating the expression levels of genes and lincRNAs in tumor and normal breast samples. Among the expressed mRNA, lincRNAs and miRNAs, 38,828 pairs are predicted to have regulatory relationships by miRNA target prediction algorithms. Among these potential target pairs, 1,742 and 213 mRNA-miRNA and lincRNA-miRNA pairs also showed significant reverse correlation, as derived from the GLM model (FDR<0.2).
As a class of newly discovered non-coding RNA, genomic loci of lincRNAs are not well annotated. As an ab initio study, we inferred their sequences, strands, exon positions, and regulatory elements from RNA-seq data. At the current stage, the functions and regulatory mechanisms of most lincRNAs are unexplored. Several well-characterized lincRNAs are involved in the establishment of the chromatin state and it is believed that many lincRNAs execute their functions by closely associating with chromatin-modifying proteins . This suggests trans-acting may be the primarily method for the affects of lincRNAs on gene expression .
Several recent studies revealed that thousands of lincRNAs exist in the cell [2, 3, 37, 42]. As a transcript, the nature of lincRNAs are quite similar to mRNA, for instance, at least a fraction of lincRNAs are poly-adenylated and have exonic structures. This as evidence of the possibility that miRNA could bind to lincRNAs and trigger degradation. The repression of lincRNAs could be an unknown part of miRNA regulation. Since some lincRNAs have been observed to be differentially expressed in tumor cells compared to normal tissues , the miRNA-derived dysregulation of lincRNA expression can be another potential mechanism of tumorigenesis.
MiRNAs regulate degradation of its target mRNA through binding on the 3'-un-translated regions. Since the entire noncoding RNA is un-translated, the whole mature transcript can serve as potential miRNA target sites. We therefore predicted miRNA binding sites on the entire exonic regions of the identified lincRNAs.
One advantage of RNA-seq data is that it can simultaneously measure the expression level of both mRNAs/lincRNAs and precursor miRNAs in the same sample. This provides opportunity for revealing the relationship of the mRNA/lincRNA-miRNA in their expression levels. As a precursor stage of mature miRNA, pre-miRNA expression levels can be different from their final functional form. Therefore, major conclusions from this analysis are subject to further experimental validation.
Among 30,069 mRNA-miRNA pairs whose relationships were established based on target prediction alone, 31 pairs were experimentally validated, based on the records documented in the miR2Disease database and TarBase database [54, 55]. Ten of these 31 (32.3%) pairs were among the 1,742 mRNA-miRNA pairs showed significant reverse correlation. Reported in 5 previous studies [56–60], these validated mRNA-miRNA pairs involve 5 miRNAs (hsa-miR-155, hsa-miR-29c, hsa-let-7b, hsa-miR-17, and hsa-miR-222) and 10 genes. In addition, we found many miRNAs showed difference in expression relationship with their potential target lincRNAs in tumor and normal tissues (Figure 4). This suggests that miRNA-lincRNA target can be specific to different cellular states.
RNA-seq data analysis
cDNA libraries from 10 normal breast tissues from the Susan G. Komen Tissue Bank at the IU Simon Cancer Center and 10 triple-negative breast cancer tumors were sequenced on an Applied Biosystems (ABI) SOLiD3 sequencer. We used a customized pipeline for RNA-seq data analysis, which includes three steps, QC filtering, sequence alignment, and gene expression quantification. 1) QC filtering. We first used SOLiD™ Instrument Control Software (ICS) and SOLiD™ Experiment Tracking System (SETS) software for the read quality recalibration. Sequences containing more than two 'N' or wildcards were discarded. Each sequence was scanned for low quality regions, and if a 5 base sliding window had an average quality score less than 20, the read would be truncated at that position. Any read with a length of less than 35 bases was discarded. Our experience suggests that such strategy effectively eliminates low quality reads while retaining high quality regions. 2) Alignment. We used BFAST  as our primary alignment algorithm, because it has high sensitivity for color space data . We used a TopHat-like strategy  to align the sequencing reads that cross splicing junctions. After aligning the sequence reads to a filtering index including repeats, rRNAs (ribosomal RNA), and other sequences not of interest, we conduct sequence alignment on three levels: genomic, known junctions, and novel junctions (based on the enriched regions identified in the genomic alignment). 3) Gene expression level quantification. The expression levels of protein-coding genes and pre-miRNAs were quantified based on the total number of RNA-seq reads falling into their genomic regions. The data were further normalized as RPKM (reads per kilo-base exon model per million mappable reads) based on their length and sequencing depth.
In potential lincRNA regions, we systematically searched regions with short read enrichment as putative lincRNA exons. A putative exon must contain at least 8 reads. If the distance between two exons was less than 10 bases, the two exons would be merged together. To define the exact exon boundaries, we searched for potential splicing sites around the putative exons. Splicing donor/acceptor (GT-AG or CT-AC for RNAs on Watson or Crick strand respectively) sites within 25 bases from the 5'- and 3'- ends of two putative exons would be considered as possible splicing sites. A lincRNA junction library was constructed by connecting all the possible donor and acceptor sites within the 100,000 bases span. The minimum intronic length is set to be 70 bases. After aligning the remaining unmapped reads to the new lincRNA junction library, we used Scripture  for predicting the lincRNA exonic structures with default settings. Similar as for mRNAs and pre-miRNAs, the expression levels of lincRNAs were further quantified by the total number of sequencing reads falling into the lincRNA exonic regions. RPKM values were calculated for each lincRNA.
Conservation, protein coding potential, and miRNA target prediction
We calculated the conservation score across placental mammals using Siphy  on various genomic regions, including exons and introns for protein coding genes, lincRNA exons and introns, random genomic regions, and predicted miRNA binding sites on lincRNAs. We used the as the conservation score, which models the overall mutation rates across placental mammals . The input of Siphy is a multiple alignment file (MAF) and a placental mammals model file. All known gene annotations (RefSeq), genome sequences (hg18 and mm9), liftOver tool, multiple alignment files, and placental mammals' model file were downloaded from the UCSC Genome browser.
We used the Codon Substantial Frequency (CSF) score to measure the coding potential of candidate lincRNA regions [49, 50, 63]. We set a 90 bases (30 Amino Acids) sliding window across the candidate lincRNA regions. The maximum score of any window of all six possible reading frames were used as the CSF score of the whole region . A higher CSF score indicates higher potential for encoding protein-coding genes.
Generalized linear model, Wilcoxon test, and false discovery rate
In which denotes the number of reads that fall into the regions of gene/lincRNA; denotes pre-miRNA expression, represented by the natural logarithm value of reads number normalized by total mapped reads, ; is a categorical variable having two possible values, indicating normal or tumor status; is an interaction factor that represents the influence of the biological conditions to miRNA regulation patterns; is the natural logarithm value of the total reads number as an offset; are intercept and slopes of the predictors and their combinations; denotes the variation that cannot be explained by the independent factors.
To allow over-dispersion, in practice, we employed a quasi-Poisson distribution. For the significance tests on the global gene/lincRNA-miRNA expression correlation, we employed a t-test on the miRNA expression coefficient . If a gene/lincRNA-miRNA pair did not show significance for , which means the relationship between miRNA and its target is not affected by the cellular states (tumor or normal), we conducted GLM model by removing the interaction term.
To correct for multiple hypothesis testing, false discovery rate (FDR) was calculated based on Benjamini-Hochberg correction . All statistical models and tests (GLM, Wilcoxon test, Chi-square test) were implemented in R: http://www.r-project.org/.
This work was supported by grant from the National Natural Science Foundation of China (60901075, 60973078), the Natural Science Foundation of Heilongjiang Province of China (LC2009C35), the Special Research Foundation for technological Innovation Talents of Harbin (2011RFLXG011), and funds from the United States National Cancer Institute grants CA113001.
This article has been published as part of BMC Medical Genomics Volume 6 Supplement 1, 2013: Proceedings of the 2011 International Conference on Bioinformatics and Computational Biology (BIOCOMP'11). The full contents of the supplement are available online at http://www.biomedcentral.com/bmcmedgenomics/supplements/6/S1. Publication of this supplement has been supported by the International Society of Intelligent Biological Medicine.
- The ENCODE (ENCyclopedia Of DNA Elements) Project. Science. 2004, 306 (5696): 636-40.Google Scholar
- Guttman M, et al: Chromatin signature reveals over a thousand highly conserved large non-coding RNAs in mammals. Nature. 2009, 458 (7235): 223-7. 10.1038/nature07672.PubMed CentralView ArticlePubMedGoogle Scholar
- Khalil AM, et al: Many human large intergenic noncoding RNAs associate with chromatin-modifying complexes and affect gene expression. Proc Natl Acad Sci USA. 2009, 106 (28): 11667-72. 10.1073/pnas.0904715106.PubMed CentralView ArticlePubMedGoogle Scholar
- Huarte M, et al: A large intergenic noncoding RNA induced by p53 mediates global gene repression in the p53 response. Cell. 2010, 142 (3): 409-19. 10.1016/j.cell.2010.06.040.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen Y, Dhupelia A, Schoenherr CJ: The Igf2/H19 imprinting control region exhibits sequence-specific and cell-type-dependent DNA methylation-mediated repression. Nucleic Acids Res. 2009, 37 (3): 793-803. 10.1093/nar/gkn985.PubMed CentralView ArticlePubMedGoogle Scholar
- Jia H, et al: Genome-wide computational identification and manual annotation of human long noncoding RNA genes. RNA. 2010, 16 (8): 1478-87. 10.1261/rna.1951310.PubMed CentralView ArticlePubMedGoogle Scholar
- Alexander RP, et al: Annotating non-coding regions of the genome. Nat Rev Genet. 2010, 11 (8): 559-71. 10.1038/nrg2814.View ArticlePubMedGoogle Scholar
- Parker BJ, et al: New families of human regulatory RNA structures identified by comparative analysis of vertebrate genomes. Genome Res. 2011, 21 (11): 1929-43. 10.1101/gr.112516.110.PubMed CentralView ArticlePubMedGoogle Scholar
- Barsotti AM, Prives C: Noncoding RNAs: the missing "linc" in p53-mediated repression. Cell. 2010, 142 (3): 358-60. 10.1016/j.cell.2010.07.029.View ArticlePubMedGoogle Scholar
- Gupta RA, et al: Long non-coding RNA HOTAIR reprograms chromatin state to promote cancer metastasis. Nature. 2010, 464 (7291): 1071-6. 10.1038/nature08975.PubMed CentralView ArticlePubMedGoogle Scholar
- Tsai MC, et al: Long noncoding RNA as modular scaffold of histone modification complexes. Science. 2010, 329 (5992): 689-93. 10.1126/science.1192002.PubMed CentralView ArticlePubMedGoogle Scholar
- Rinn JL, et al: Functional demarcation of active and silent chromatin domains in human HOX loci by noncoding RNAs. Cell. 2007, 129 (7): 1311-23. 10.1016/j.cell.2007.05.022.PubMed CentralView ArticlePubMedGoogle Scholar
- Amaral PP, Mattick JS: Noncoding RNA in development. Mamm Genome. 2008, 19 (7-8): 454-92. 10.1007/s00335-008-9136-7.View ArticlePubMedGoogle Scholar
- Clark MB, et al: Genome-wide analysis of long noncoding RNA stability. Genome Res. 2012, 22 (5): 885-98. 10.1101/gr.131037.111.PubMed CentralView ArticlePubMedGoogle Scholar
- Nagano T, Fraser P: No-nonsense functions for long noncoding RNAs. Cell. 2011, 145 (2): 178-81. 10.1016/j.cell.2011.03.014.View ArticlePubMedGoogle Scholar
- Guttman M, et al: lincRNAs act in the circuitry controlling pluripotency and differentiation. Nature. 2011, 477 (7364): 295-300. 10.1038/nature10398.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhao J, et al: Polycomb proteins targeted by a short repeat RNA to the mouse X chromosome. Science. 2008, 322 (5902): 750-6. 10.1126/science.1163045.PubMed CentralView ArticlePubMedGoogle Scholar
- Nagano T, et al: The Air noncoding RNA epigenetically silences transcription by targeting G9a to chromatin. Science. 2008, 322 (5908): 1717-20. 10.1126/science.1163802.View ArticlePubMedGoogle Scholar
- Cabianca DS, et al: A long ncRNA links copy number variation to a polycomb/trithorax epigenetic switch in FSHD muscular dystrophy. Cell. 2012, 149 (4): 819-31. 10.1016/j.cell.2012.03.035.PubMed CentralView ArticlePubMedGoogle Scholar
- Tripathi V, et al: The nuclear-retained noncoding RNA MALAT1 regulates alternative splicing by modulating SR splicing factor phosphorylation. Mol Cell. 2010, 39 (6): 925-38. 10.1016/j.molcel.2010.08.011.PubMed CentralView ArticlePubMedGoogle Scholar
- Novikova IV, Hennelly SP, Sanbonmatsu KY: Structural architecture of the human long non-coding RNA, steroid receptor RNA activator. Nucleic Acids Res. 2012, 40 (11): 5034-51. 10.1093/nar/gks071.PubMed CentralView ArticlePubMedGoogle Scholar
- Huarte M, Rinn JL: Large non-coding RNAs: missing links in cancer?. Hum Mol Genet. 2010, 19 (R2): R152-61. 10.1093/hmg/ddq353.PubMed CentralView ArticlePubMedGoogle Scholar
- Moran VA, Perera RJ, Khalil AM: Emerging functional and mechanistic paradigms of mammalian long non-coding RNAs. Nucleic Acids Res. 2012, 40 (14): 6391-400. 10.1093/nar/gks296.PubMed CentralView ArticlePubMedGoogle Scholar
- Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-97. 10.1016/S0092-8674(04)00045-5.View ArticlePubMedGoogle Scholar
- O'Connell RM, et al: Sustained expression of microRNA-155 in hematopoietic stem cells causes a myeloproliferative disorder. J Exp Med. 2008, 205 (3): 585-94. 10.1084/jem.20072108.PubMed CentralView ArticlePubMedGoogle Scholar
- Griffiths-Jones S: The microRNA Registry. Nucleic Acids Res. 2004, 32 (Database issue): D109-11.PubMed CentralView ArticlePubMedGoogle Scholar
- Brennecke J, et al: bantam encodes a developmentally regulated microRNA that controls cell proliferation and regulates the proapoptotic gene hid in Drosophila. Cell. 2003, 113 (1): 25-36. 10.1016/S0092-8674(03)00231-9.View ArticlePubMedGoogle Scholar
- Krichevsky AM, et al: A microRNA array reveals extensive regulation of microRNAs during brain development. RNA. 2003, 9 (10): 1274-81. 10.1261/rna.5980303.PubMed CentralView ArticlePubMedGoogle Scholar
- Xie X, et al: Systematic discovery of regulatory motifs in human promoters and 3' UTRs by comparison of several mammals. Nature. 2005, 434 (7031): 338-45. 10.1038/nature03441.PubMed CentralView ArticlePubMedGoogle Scholar
- Wienholds E, et al: MicroRNA expression in zebrafish embryonic development. Science. 2005, 309 (5732): 310-1. 10.1126/science.1114519.View ArticlePubMedGoogle Scholar
- Hansen TB, et al: miRNA-dependent gene silencing involving Ago2-mediated cleavage of a circular antisense RNA. EMBO J. 2011, 30 (21): 4414-22. 10.1038/emboj.2011.359.PubMed CentralView ArticlePubMedGoogle Scholar
- Braconi C, et al: microRNA-29 can regulate expression of the long non-coding RNA gene MEG3 in hepatocellular cancer. Oncogene. 2011, 30 (47): 4750-6. 10.1038/onc.2011.193.PubMed CentralView ArticlePubMedGoogle Scholar
- Guo H, et al: Mammalian microRNAs predominantly act to decrease target mRNA levels. Nature. 2010, 466 (7308): 835-40. 10.1038/nature09267.PubMed CentralView ArticlePubMedGoogle Scholar
- Chi SW, et al: Argonaute HITS-CLIP decodes microRNA-mRNA interaction maps. Nature. 2009, 460 (7254): 479-86.PubMed CentralPubMedGoogle Scholar
- Hendrickson DG, et al: Systematic identification of mRNAs recruited to argonaute 2 by specific microRNAs and corresponding changes in transcript abundance. PLoS One. 2008, 3 (5): e2126-10.1371/journal.pone.0002126.PubMed CentralView ArticlePubMedGoogle Scholar
- Harrow J, et al: GENCODE: producing a reference annotation for ENCODE. Genome Biol. 2006, 7 (Suppl 1): S4 1-9.View ArticleGoogle Scholar
- Bu D, et al: NONCODE v3.0: integrative annotation of long noncoding RNAs. Nucleic Acids Res. 2012, 40 (Database issue): D210-5.PubMed CentralView ArticlePubMedGoogle Scholar
- Jeggari A, Marks DS, Larsson E: miRcode: a map of putative microRNA target sites in the long non-coding transcriptome. Bioinformatics. 2012, 28 (15): 2062-3. 10.1093/bioinformatics/bts344.PubMed CentralView ArticlePubMedGoogle Scholar
- Prensner JR, et al: Transcriptome sequencing across a prostate cancer cohort identifies PCAT-1, an unannotated lincRNA implicated in disease progression. Nat Biotechnol. 2011, 29 (8): 742-9. 10.1038/nbt.1914.PubMed CentralView ArticlePubMedGoogle Scholar
- Trapnell C, Pachter L, Salzberg SL: TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009, 25 (9): 1105-11. 10.1093/bioinformatics/btp120.PubMed CentralView ArticlePubMedGoogle Scholar
- Trapnell C, et al: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010, 28 (5): 511-5. 10.1038/nbt.1621.PubMed CentralView ArticlePubMedGoogle Scholar
- Guttman M, et al: Ab initio reconstruction of cell type-specific transcriptomes in mouse reveals the conserved multi-exonic structure of lincRNAs. Nat Biotechnol. 2010, 28 (5): 503-10. 10.1038/nbt.1633.PubMed CentralView ArticlePubMedGoogle Scholar
- Homer N, Merriman B, Nelson SF: BFAST: an alignment tool for large scale genome resequencing. PLoS One. 2009, 4 (11): e7767-10.1371/journal.pone.0007767.PubMed CentralView ArticlePubMedGoogle Scholar
- Kent WJ, et al: The human genome browser at UCSC. Genome Res. 2002, 12 (6): 996-1006.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang G, et al: RNA polymerase II binding patterns reveal genomic regions involved in microRNA gene regulation. PLoS One. 2010, 5 (11): e13798-10.1371/journal.pone.0013798.PubMed CentralView ArticlePubMedGoogle Scholar
- Hiller D, et al: Identifiability of isoform deconvolution from junction arrays and RNA-Seq. Bioinformatics. 2009, 25 (23): 3056-9. 10.1093/bioinformatics/btp544.PubMed CentralView ArticlePubMedGoogle Scholar
- Jiang H, Wong WH: Statistical inferences for isoform expression in RNA-Seq. Bioinformatics. 2009, 25 (8): 1026-32. 10.1093/bioinformatics/btp113.PubMed CentralView ArticlePubMedGoogle Scholar
- Li B, et al: RNA-Seq gene expression estimation with read mapping uncertainty. Bioinformatics. 2010, 26 (4): 493-500. 10.1093/bioinformatics/btp692.PubMed CentralView ArticlePubMedGoogle Scholar
- Lin MF, et al: Revisiting the protein-coding gene catalog of Drosophila melanogaster using 12 fly genomes. Genome Res. 2007, 17 (12): 1823-36. 10.1101/gr.6679507.PubMed CentralView ArticlePubMedGoogle Scholar
- Lin MF, et al: Performance and scalability of discriminative metrics for comparative gene identification in 12 Drosophila genomes. PLoS Comput Biol. 2008, 4 (4): e1000067-10.1371/journal.pcbi.1000067.PubMed CentralView ArticlePubMedGoogle Scholar
- Friedman RC, et al: Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 2009, 19 (1): 92-105.PubMed CentralView ArticlePubMedGoogle Scholar
- Kertesz M, et al: The role of site accessibility in microRNA target recognition. Nat Genet. 2007, 39 (10): 1278-84. 10.1038/ng2135.View ArticlePubMedGoogle Scholar
- Tsai MC, Spitale RC, Chang HY: Long intergenic noncoding RNAs: new links in cancer progression. Cancer Res. 2011, 71 (1): 3-7. 10.1158/0008-5472.CAN-10-2483.PubMed CentralView ArticlePubMedGoogle Scholar
- Jiang Q, et al: miR2Disease: a manually curated database for microRNA deregulation in human disease. Nucleic Acids Res. 2009, 37 (Database issue): D98-104.PubMed CentralView ArticlePubMedGoogle Scholar
- Papadopoulos GL, et al: The database of experimentally supported targets: a functional update of TarBase. Nucleic Acids Res. 2009, 37 (Database issue): D155-8.PubMed CentralView ArticlePubMedGoogle Scholar
- Medina R, et al: MicroRNAs 221 and 222 bypass quiescence and compromise cell survival. Cancer Res. 2008, 68 (8): 2773-80. 10.1158/0008-5472.CAN-07-6754.PubMed CentralView ArticlePubMedGoogle Scholar
- Selbach M, et al: Widespread changes in protein synthesis induced by microRNAs. Nature. 2008, 455 (7209): 58-63. 10.1038/nature07228.View ArticlePubMedGoogle Scholar
- Sengupta S, et al: MicroRNA 29c is down-regulated in nasopharyngeal carcinomas, up-regulating mRNAs encoding extracellular matrix proteins. Proc Natl Acad Sci USA. 2008, 105 (15): 5874-8. 10.1073/pnas.0801130105.PubMed CentralView ArticlePubMedGoogle Scholar
- Bolisetty MT, et al: Reticuloendotheliosis virus strain T induces miR-155, which targets JARID2 and promotes cell survival. J Virol. 2009, 83 (23): 12009-17. 10.1128/JVI.01182-09.PubMed CentralView ArticlePubMedGoogle Scholar
- Inomata M, et al: MicroRNA-17-92 down-regulates expression of distinct targets in different B-cell lymphoma subtypes. Blood. 2009, 113 (2): 396-402.View ArticlePubMedGoogle Scholar
- Fu Y, et al: SOLiD™ Sequencing and 2-Base Encoding. The Biology of Genomes Meeting. 2008Google Scholar
- Garber M, et al: Identifying novel constrained elements by exploiting biased substitution patterns. Bioinformatics. 2009, 25 (12): i54-62. 10.1093/bioinformatics/btp190.PubMed CentralView ArticlePubMedGoogle Scholar
- Clamp M, et al: Distinguishing protein-coding and noncoding genes in the human genome. Proc Natl Acad Sci USA. 2007, 104 (49): 19428-33. 10.1073/pnas.0709013104.PubMed CentralView ArticlePubMedGoogle Scholar
- Kozomara A, Griffiths-Jones S: miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res. 2011, 39 (Database issue): D152-7.PubMed CentralView ArticlePubMedGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society. Series B (Methodological). 1995, 289-300.Google Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.