Clustered microRNAs' coordination in regulating protein-protein interaction network
- Xiongying Yuan†1, 2,
- Changning Liu†1,
- Pengcheng Yang1, 2,
- Shunmin He1, 2,
- Qi Liao1, 3,
- Shuli Kang1 and
- Yi Zhao1Email author
© Yuan et al; licensee BioMed Central Ltd. 2009
Received: 12 November 2008
Accepted: 26 June 2009
Published: 26 June 2009
MicroRNAs (miRNAs), a growing class of small RNAs with crucial regulatory roles at the post-transcriptional level, are usually found to be clustered on chromosomes. However, with the exception of a few individual cases, so far little is known about the functional consequence of this conserved clustering of miRNA loci. In animal genomes such clusters often contain non-homologous miRNA genes. One hypothesis to explain this heterogeneity suggests that clustered miRNAs are functionally related by virtue of co-targeting downstream pathways.
Integrating of miRNA cluster information with protein protein interaction (PPI) network data, our research supports the hypothesis of the functional coordination of clustered miRNAs and links it to the topological features of miRNAs' targets in PPI network. Specifically, our results demonstrate that clustered miRNAs jointly regulate proteins in close proximity of the PPI network. The possibility that two proteins yield to this coordinated regulation is negatively correlated with their distance in PPI network. Guided by the knowledge of this preference, we found several network communities enriched with target genes of miRNA clusters. In addition, our results demonstrate that the variance of this propensity can also partly be explained by protein's connectivity and miRNA's conservation.
In summary, this work supports the hypothesis of intra-cluster coordination and investigates the extent of this coordination.
MicroRNAs (miRNAs) are small (22 nt) single-stranded non-coding RNA (ncRNA) molecules. They are processed from hairpin precursors of approximately 70 nt (pre-miRNAs), which in turn are extracted from primary transcripts (pri-miRNA) [1, 2]. miRNAs can repress gene expression post-transcriptionally by binding to the 3' untranslated regions (3' UTRs) of their target mRNAs [3–5]. In animal genomes, miRNAs are found in various genomic locations; while most are located in intergenic regions, some are found to be hosted within the introns of pre-mRNAs or within longer ncRNA genes [6, 7]. Interestingly, known miRNA genes, both hosted and non-hosted, are often observed to be clustered . A cluster usually includes two or three miRNA genes, but larger clusters have also been identified, including a human mir-17 cluster comprising 6 miRNA genes and a human mir-302 cluster comprising 8 miRNA genes. Baskerville and Bartel found that miRNAs within 50 Kb are highly correlated in expression across 24 different human organs, indicating that 50 Kb might be used as preliminary definition for miRNA cluster .
Although miRNAs in a given cluster are usually phylogenetically related in sequence, many human clusters containing miRNAs without apparent sequence homology are also found . A plausible but yet-to-be validated possibility is that the clustered miRNAs are functionally related by virtue of targeting the same gene or different genes in the same pathway [8, 10–12]. Studies of individual cases have to some extent supported this hypothesis. For instance, mir-15a-16 cluster act as tumor suppressor genes in prostate cancer by controlling cell survival, proliferation and invasion . In gastric cancer, mir-106b cluster and mir-222 cluster are upregulated and modulate cell cycle by targeting the Cip/Kip family proteins . The frequently studied mir-17-92 cluster acts not only in tumor formation but also in development of heart, lungs, and immune system [15–18]. Three individual miRNAs of the mir-379-410 cluster are required for activity-dependent development of hippocampal neurons . mir-48 and mir-241 of a let-7 family cluster function together to regulate developmental timing in Caenorhabditis elegans . However, apart from these individual cases, the hypothesis has not yet been systematically investigated.
Here, we conduct a comprehensive analysis of the conservation and functional consequence of miRNA clustering. We show that clustered miRNAs have functional relatedness through co-targeting proteins in close proximity of PPI network.
MiRNA clusters' construction and property
It has been reported that clusters may contain unrelated miRNAs [21, 22]. But the frequency of such heterogeneity has not been comprehensively evaluated. Prior research  has identified the clustering of 326 registered human miRNA genes, using the strategy proposed by Altuvia et al. . According to miRBase , the number of registered miRNA genes has doubled. Still more homologs of the known miRNAs are expected in human genome, especially when more and more miRNA genes derived from TEs (transcriptional elements) have recently been identified [25–27].
In order to analyze the genomic clustering of miRNAs as systematically and exhaustively as possible, we conducted a cross-species Blast search of known pre-miRNA sequences in the entire human, mouse and rat genomes. In rat genome we found 152 novel miRNA loci that are homologues of previously annotated human or mouse miRNAs, whereas the corresponding numbers of loci in the mouse and human genomes are 33 and 19, respectively (see Additional File 1). We clustered miRNA genes within 50 kb of each other according to Bartel et al. . After filtering with RNA polymerase II promoter analysis software and expression profiles (for details, see Material and Methods), a total of 104 miRNA clusters representing 67 unique miRNA combinations have been identified in human genome. Clusters composed of only hsa-mir-566 were not included, since hsa-mir-566 is thought to be derived from repeat sequence Alu, which is pervasive in the human genome . Similarly, 91 clusters representing 64 unique miRNA combinations were found in mouse (mmu-mir-680 excluded), and 56 miRNA clusters, all unique, were identified in rat.
In recent releases of miRBase, similar miRNA precursors are assembled into miRNA families (miFams) based on computational analysis and manual inspection. miRNAs in a family have been shown to be phylogenetically related [29, 30]. We further mapped the miRNAs in our data to their families according to miRBase's assignment. As results, 55 distinct miRNA clusters were identified at family level in human genome, along with 51 in mouse, and 42 in rat (see Additional File 2). After mapping to families, 31 clusters have been found to be conserved among the three species, which take up 56% in human, 61% in mouse and 74% in rat (see Additional File 3).
Heterogeneity analysis of the miRNA clusters.
Functional Coordination on PPI network
Coordination among regulators from the transcriptional and miRNA regulatory layers has recently drawn great interests [32, 33]. To reveal why non-homologous miRNAs are found in the same cluster, we analyzed their potential coordination in the context of PPI networks. Generally, the coordination of miRNAs in a cluster should be analyzed directly through justifying the relative proximity between their targets in the network. However, because of the incompleteness of credible PPI networks (e.g. Human Protein Reference Database, HPRD) and the low reliability of other integrated databases, we primarily chose to reverse the procedure to analyze whether proteins in close proximity incline to be regulated by miRNAs from the same cluster (henceforth referred to as "sc-miRNAs").
We first investigated directly interacting protein pairs. It has recently been reported that interacting protein pairs have a certain propensity to be regulated by the same miRNA . There may exist a similar tendency that two interacting proteins are also regulated by sc-miRNAs. Sc-miRNAs can fulfil the same function through a complementary approach – by assembling two non-related miRNA loci in the same transcriptional unit. Before analysis, we prepare the data as follows:
We used HPRD as the golden standard positive (GSP, interacting protein pairs) dataset, and the dataset generated by Rhodes et al.  according to protein subcellular localization as the golden standard negative (GSN, non-interacting protein pairs) dataset. Comparisons of various miRNA target prediction programs suggest that TargetScan could achieve both high sensitivity and specificity [36, 37]. We therefore employed TargetScan predictions for clustered miRNAs. For targets of a miRNA family, common targets shared by all its miRNA members were used. The number of miRNA sites is 9.7 ± 8.9 per transcript for GSP, and 10.1 ± 9.4 for GSN. To specifically evaluate the coordination among different miRNA families in a cluster, we constructed a family-represented version of hetero-clusters by mapping miRNAs in a cluster to their families, and then having redundancies removed (We will later denote these clusters as family-represented hetero-clusters and clusters before collapsing as original clusters). As a result, the 53 hetero-clusters have been collapsed to 37 family-represented hetero-clusters. To curb possible bias in the target prediction and isolate the effects of clustering, a control was generated by randomly reshuffling miRNAs (or families, when studying family-represented hetero-clusters) among the clusters.
Since the interacting pairs demonstrate stronger propensity to be regulated by sc-miRNAs than non-interacting pairs, we next asked whether there are some correlation between the propensity and the distance between two proteins in the network. The shortest paths were calculated between all the proteins in GSP dataset. We found a significant negative correlation between the number of clusters regulating the two proteins and the length of the shortest path between them in the PPI network (Spearman's rank correlation, correlation coefficient Rs = -0.24, P < 2.2e-16). This result suggests that distance in the PPI network is an important factor with respect to the coordinated regulation from sc-miRNAs, and lend further support to our previous finding.
We have also noticed the great variance in the number of clusters regulating an interacting protein pair. Still 50% of the directly interacting protein pairs are not co-regulated by miRNAs from same cluster, and the numbers of clusters that do regulate them vary from 1 to 14. Despite the noises and data incompleteness, directly interacting pairs are still likely to enjoy different levels of tendency. Therefore, we studied another factor, connectivity, that might influence this propensity.
A positive correlation between the sum of each two proteins' connectivity and the number of clusters regulating them.
To see whether these observations are robust, the same analysis was performed with PicTar prediction set, the overlap of TargetScan and PicTar prediction set as recommended by Sethupathy , as well as miRanda prediction set with alignment score greater than 155 (see Additional File 4). They all gave very consistent results.
Modularity of the targets
Significant network communities that are enriched with target genes of miRNA clusters.
ACTA1, DRP2, SNTA1, DMD, DAG1, SNTB2, UTRN, SNTG1, SNTG2, KCNJ12, SNTB1, PGM5P1, DTNA, DTNB, DGKZ
muscle development and contraction
AR, ESR1, RB1, TP53, CREBBP, EP300, JUN, STAT3, NCOA1, SMAD4, SMAD2, SMAD3, BRCA1, STAT1, TRIP4, RELA, CCND1, SP1
positive regulation of transcription
ARNTL, PER1, TIMELESS, CRY1, CLOCK, PER3, CRY2, PER2, CSNK1D, CSNK1E
HDAC2, RBBP7, RBP1, SAP30, HDAC1, ING1, RBBP4, BRMS1, SIN3A, BRMS1L, C2ORF59
TP53, CUL5, COPS3, GPS1, COPS2, COPS8, COPS6, COPS5, COPS4, COPS7A
cancer suppressor p53 with COP9 signalosome
Although our work mainly focus on protein-protein interaction data, for it is more complete and more suitable for global analysis, we re-ran our analysis with KEGG human pathway data, and presented the results in Additional File 4.
Coordination and conservation
Cluster provides a neat mechanism to transcribe many cooperative miRNAs simultaneously. It could hardly be imaged that this exquisite coordination has been generated in a single event; rather, it should have been consolidated through a long period of evolutionary process. To investigate the relationship between coordination and conservation, we have to calculate the target interactions between miRNAs in a cluster.
Although results in previous section indicate that miRNAs in a cluster may co-regulate proteins in close proximity of PPI network, they provide no information on each miRNA cluster's coordination. Here, those highly coordinated miRNA pairs (-lgP > 1.5), conserved and non-conserved, come from 17 (30.1%) hetero-clusters. We have found literature support for the function of 15 of them, all listed in Additional File 2. Thus it can be estimated that at least 30% of the miRNA hetero-clusters we can find evidence in current HPRD network for their intra-cluster coordination.
Recently, global analysis of miRNA targets in the context of sets of genes or PPI network has drawn great interests [34, 38, 44–47]. It could also provide important insights into how miRNAs in a cluster cooperate with each other. In this work, we adopted a protein-centered perspective to start our analysis. We fixed a protein set and analyzed the relationship of miRNAs that target them. This perspective has advantages in avoiding the incompleteness of protein-protein interaction data, and allowing us to analyze miRNAs' functionality according to the features of their targets.
By using the protein-centered perspective, we observed that directly interacting proteins incline to be regulated by miRNAs in the same cluster. The closer these proteins are in the PPI network, the more likely the targeting miRNAs are located in the same cluster. Proteins usually fulfill certain functions by means of interaction. If their functions are temporarily not needed, then their expressions should be reduced simultaneously. miRNAs in a cluster transcribed as polycistron would provide a competent and efficient mechanism to achieve this goal. Our further analysis reveals that connectivity is another factor that matters only when proteins interact or are at least close in PPI network. Interacting proteins, both with great connectivity, usually connect two functional protein groups. They would be under stronger selective pressure to consent to a miRNA cluster's regulation. If one of them was not regulated in time, the module it involved would likely be still in function, and the coordination of the whole system would be greatly impaired.
The finding that distance is more important than connectivity suggests that sc-miRNAs preferentially co-regulate proteins in close proximity of a PPI network. However, this may also be due to the fact that PPI network are mainly gathered from high-throughput experiments which focus on physical interactions, such as yeast two hybrid system and tandem affinity purification. Therefore, for two proteins with long distance in PPI network, even their concurrence in time and space could not be guaranteed. To overcome this limitation, other information, such as signaling network or transcriptional regulatory network, would be needed [46, 47].
This study supports the putative hypothesis of internal coordination among sc-miRNAs to regulate downstream biological networks. The linkage of sc-miRNAs' functional coordination and their targets' topological features we found highlights the potential to further investigate their subtle relationship.
Datasets of PPI
Gold standard positive dataset: all the information in the Human Protein Reference Database (HPRD) has been manually extracted from the literature by expert biologists who read, interpret and analyze the published data. It contains 38,167 distinct interactions among 9,465 proteins . Our analysis focused exclusively on the giant connected component of 9,134 proteins and 31,909 interactions. The remaining proteins were tiny clusters with sizes between two and eight. According to TargetScan's prediction, 8,264 of the 9,134 proteins are targeted by miRNAs. 6,737 of them are targeted by clustered miRNAs, and 4,835 are targeted by miRNAs that appear in hetero-clusters.
Gold standard negative dataset: it was generated by Rhodes et al. , which includes all possible pair-wise combinations between two sets of proteins that are assigned a subcellular localization of the plasma membrane (1,397 proteins) and the nucleus (2,224 proteins), respectively, by the Gene Ontology (GO) Consortium. There are 2633 overlapping proteins between GSP and GSN.
miRNA target prediction
TargetScan's release 4.2 was downloaded from TargetScan site http://www.targetscan.org/. miRanda's prediction for 418 human miRNAs in microRNA registry release 10.0 were downloaded from microRNA.org Jan 2008 release with align score higher than 155. PicTar prediction was downloaded from UCSC genome browser PicTar miRNA track.
Determination of miRNA clusters
We download all pre-miRNA sequences of Homo sapiens (hsa, 533 miRNAs), Mus musculus (mmu, 442 miRNAs) and Rattus norvegicus (rno, 285 miRNAs) from the miRNA registry release 10.0 http://www.sanger.ac.uk/software/Rfam/mirna/, and blast them against the whole genome of human, mouse and rat. Only the hits with 90% coverage and 90% identity were selected. All miRNA genes whose Blast distance is smaller than 100 bp were treated as one locus. According to Baskerville and Bartel , we considered loci within 50 kb of each other as belonging to the same cluster. miRNAs in a cluster were further screened using following filters:RNA polymerase II filter: Clusters were classified as intronic or intergenic according to their hostages. miRNA genes located at the introns or UTRs of protein-coding genes are transcribed together with their host genes, whereas miRNA genes dispersed in an intergenic region are generally believed to be transcribed by RNA polymerase II . Therefore, we screened the regions surrounding intergenic miRNA clusters using the following RNA polymerase II promoter analysis software: Promoter Scan http://thr.cit.nih.gov/molbio/proscan for predicting promoter regions based on scoring homologies with putative eukaryotic Pol II promoter sequences, Promoter Prediction 2.0 http://www.cbs.dtu.dk/services/Promoter/ for predicting transcription start sites of vertebrate Pol II promoters in DNA sequences, and NNPP (neural network promoter prediction) http://www.fruitfly.org/seq_tools/promoter.html for predicting eukaryotic Pol II promoter sequences. Intronic miRNA clusters were directly reserved for transcription along with their host genes.
Expression profile filter
miRNAs in the clusters were also filtered according to expressional correlation. Three expression profiles [50–53] were used in the screening. Only miRNAs with a Pearson Correlation Coefficient higher than 0.3 were retained in the cluster.
Quantification of the functional coordination
In our research, we have performed two types of randomization to construct control set: one is randomization of miRNA clusters. To obtain random controls for testing the propensity of interacting proteins to be regulated by sc-miRNAs, we generated a cluster set by randomly shuffling miRNAs among clusters, while keeping the size of each cluster unchanged. The other is randomization of miRNA-target protein associations. To obtain random controls for testing the proximity of sc-miRNAs' targets, we generated targets for each miRNA by randomly shuffling the miRNA-target protein associations, while keeping unchanged the number of proteins that a miRNA targets. The empirical P-value for target interactions was calculated against 10,000 independent randomized samples. The negative logarithm of P-value (-lgP) was then used to assess the coordination. The greater -lgP is, the stronger the coordination is.
All the statistical tests and kernel density estimation were done in R. Connectivity and the shortest paths between each protein were calculated using iGraph package http://cneurocvs.rmki.kfki.hu/igraph/.
This work was supported by the National Sciences Foundation of China, under Grant Nos. 30570393 and 30600729. The authors express special thanks to Dr. Dongbo Bu and Dr. Geir Skogerbo for editing and proofreading of this manuscript. The authors also thank anonymous reviewers for their valuable discussions and comments on this manuscript.
- Lee Y, Jeon K, Lee JT, Kim S, Kim VN: MicroRNA maturation: stepwise processing and subcellular localization. Embo J. 2002, 21 (17): 4663-4670. 10.1093/emboj/cdf476PubMed CentralView ArticlePubMedGoogle Scholar
- Cullen BR: Transcription and processing of human microRNA precursors. Mol Cell. 2004, 16 (6): 861-865. 10.1016/j.molcel.2004.12.002View ArticlePubMedGoogle Scholar
- Lim LP, Glasner ME, Yekta S, Burge CB, Bartel DP: Vertebrate microRNA genes. Science. 2003, 299 (5612): 1540- 10.1126/science.1080372View ArticlePubMedGoogle Scholar
- Ambros V: The functions of animal microRNAs. Nature. 2004, 431 (7006): 350-355. 10.1038/nature02871View ArticlePubMedGoogle Scholar
- Bartel DP, Chen CZ: Micromanagers of gene expression: the potentially widespread influence of metazoan microRNAs. Nat Rev Genet. 2004, 5 (5): 396-400. 10.1038/nrg1328View ArticlePubMedGoogle Scholar
- Rodriguez A, Griffiths-Jones S, Ashurst JL, Bradley A: Identification of mammalian microRNA host genes and transcription units. Genome Res. 2004, 14 (10A): 1902-1910. 10.1101/gr.2722704PubMed CentralView ArticlePubMedGoogle Scholar
- Weber MJ: New human and mouse microRNA genes found by homology search. Febs J. 2005, 272 (1): 59-73. 10.1111/j.1432-1033.2004.04389.xView ArticlePubMedGoogle Scholar
- Kim VN, Nam JW: Genomics of microRNA. Trends Genet. 2006, 22 (3): 165-173. 10.1016/j.tig.2006.01.003View ArticlePubMedGoogle Scholar
- Baskerville S, Bartel DP: Microarray profiling of microRNAs reveals frequent coexpression with neighboring miRNAs and host genes. Rna. 2005, 11 (3): 241-247. 10.1261/rna.7240905PubMed CentralView ArticlePubMedGoogle Scholar
- Yu J, Wang F, Yang GH, Wang FL, Ma YN, Du ZW, Zhang JW: Human microRNA clusters: genomic organization and expression profile in leukemia cell lines. Biochem Biophys Res Commun. 2006, 349 (1): 59-68. 10.1016/j.bbrc.2006.07.207View ArticlePubMedGoogle Scholar
- Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297. 10.1016/S0092-8674(04)00045-5View ArticlePubMedGoogle Scholar
- Grun D, Wang YL, Langenberger D, Gunsalus KC, Rajewsky N: microRNA target predictions across seven Drosophila species and comparison to mammalian targets. PLoS Comput Biol. 2005, 1 (1): e13- 10.1371/journal.pcbi.0010013PubMed CentralView ArticlePubMedGoogle Scholar
- Bonci D, Coppola V, Musumeci M, Addario A, Giuffrida R, Memeo L, D'Urso L, Pagliuca A, Biffoni M, Labbaye C, et al.: The miR-15a-miR-16-1 cluster controls prostate cancer by targeting multiple oncogenic activities. Nat Med. 2008, 14 (11): 1271-1277. 10.1038/nm.1880View ArticlePubMedGoogle Scholar
- Kim YK, Yu J, Han TS, Park SY, Namkoong B, Kim DH, Hur K, Yoo MW, Lee HJ, Yang HK, et al.: Functional links between clustered microRNAs: suppression of cell-cycle inhibitors by microRNA clusters in gastric cancer. Nucleic Acids Res. 2009, 37 (5): 1672-1681. 10.1093/nar/gkp002PubMed CentralView ArticlePubMedGoogle Scholar
- O'Donnell KA, Wentzel EA, Zeller KI, Dang CV, Mendell JT: c-Myc-regulated microRNAs modulate E2F1 expression. Nature. 2005, 435 (7043): 839-843. 10.1038/nature03677View ArticlePubMedGoogle Scholar
- Ventura A, Young AG, Winslow MM, Lintault L, Meissner A, Erkeland SJ, Newman J, Bronson RT, Crowley D, Stone JR, et al.: Targeted deletion reveals essential and overlapping functions of the miR-17 through 92 family of miRNA clusters. Cell. 2008, 132 (5): 875-886. 10.1016/j.cell.2008.02.019PubMed CentralView ArticlePubMedGoogle Scholar
- Xiao C, Srinivasan L, Calado DP, Patterson HC, Zhang B, Wang J, Henderson JM, Kutok JL, Rajewsky K: Lymphoproliferative disease and autoimmunity in mice with increased miR-17-92 expression in lymphocytes. Nat Immunol. 2008, 9 (4): 405-414. 10.1038/ni1575PubMed CentralView ArticlePubMedGoogle Scholar
- Lu Y, Thomson JM, Wong HY, Hammond SM, Hogan BL: Transgenic over-expression of the microRNA miR-17-92 cluster promotes proliferation and inhibits differentiation of lung epithelial progenitor cells. Dev Biol. 2007, 310 (2): 442-453. 10.1016/j.ydbio.2007.08.007PubMed CentralView ArticlePubMedGoogle Scholar
- Fiore R, Khudayberdiev S, Christensen M, Siegel G, Flavell SW, Kim TK, Greenberg ME, Schratt G: Mef2-mediated transcription of the miR379-410 cluster regulates activity-dependent dendritogenesis by fine-tuning Pumilio2 protein levels. Embo J. 2009, 28 (6): 697-710. 10.1038/emboj.2009.10PubMed CentralView ArticlePubMedGoogle Scholar
- Abbott AL, Alvarez-Saavedra E, Miska EA, Lau NC, Bartel DP, Horvitz HR, Ambros V: The let-7 MicroRNA family members mir-48, mir-84, and mir-241 function together to regulate developmental timing in Caenorhabditis elegans. Dev Cell. 2005, 9 (3): 403-414. 10.1016/j.devcel.2005.07.009PubMed CentralView ArticlePubMedGoogle Scholar
- Hertel J, Lindemeyer M, Missal K, Fried C, Tanzer A, Flamm C, Hofacker IL, Stadler PF: The expansion of the metazoan microRNA repertoire. BMC Genomics. 2006, 7: 25- 10.1186/1471-2164-7-25PubMed CentralView ArticlePubMedGoogle Scholar
- Hayashita Y, Osada H, Tatematsu Y, Yamada H, Yanagisawa K, Tomida S, Yatabe Y, Kawahara K, Sekido Y, Takahashi T: A polycistronic microRNA cluster, miR-17-92, is overexpressed in human lung cancers and enhances cell proliferation. Cancer Res. 2005, 65 (21): 9628-9632. 10.1158/0008-5472.CAN-05-2352View ArticlePubMedGoogle Scholar
- Altuvia Y, Landgraf P, Lithwick G, Elefant N, Pfeffer S, Aravin A, Brownstein MJ, Tuschl T, Margalit H: Clustering and conservation patterns of human microRNAs. Nucleic Acids Res. 2005, 33 (8): 2697-2706. 10.1093/nar/gki567PubMed CentralView ArticlePubMedGoogle Scholar
- Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ: miRBase: tools for microRNA genomics. Nucleic Acids Res. 2008, D154-158. 36 Database
- Smalheiser NR, Torvik VI: Mammalian microRNAs derived from genomic repeats. Trends Genet. 2005, 21 (6): 322-326. 10.1016/j.tig.2005.04.008View ArticlePubMedGoogle Scholar
- Borchert GM, Lanier W, Davidson BL: RNA polymerase III transcribes human microRNAs. Nat Struct Mol Biol. 2006, 13 (12): 1097-1101. 10.1038/nsmb1167View ArticlePubMedGoogle Scholar
- Piriyapongsa J, Marino-Ramirez L, Jordan IK: Origin and evolution of human microRNAs from transposable elements. Genetics. 2007, 176 (2): 1323-1337. 10.1534/genetics.107.072553PubMed CentralView ArticlePubMedGoogle Scholar
- Piriyapongsa J, Jordan IK: A family of human microRNA genes from miniature inverted-repeat transposable elements. PLoS ONE. 2007, 2 (2): e203- 10.1371/journal.pone.0000203PubMed CentralView ArticlePubMedGoogle Scholar
- Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ: miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006, D140-144. 34 Database
- Huang Y, Gu X: A bootstrap based analysis pipeline for efficient classification of phylogenetically related animal miRNAs. BMC Genomics. 2007, 8: 66- 10.1186/1471-2164-8-66PubMed CentralView ArticlePubMedGoogle Scholar
- Li A, Mao L: Evolution of plant microRNA gene families. Cell Res. 2007, 17 (3): 212-218.PubMedGoogle Scholar
- Shalgi R, Lieber D, Oren M, Pilpel Y: Global and local architecture of the mammalian microRNA-transcription factor regulatory network. PLoS Comput Biol. 2007, 3 (7): e131- 10.1371/journal.pcbi.0030131PubMed CentralView ArticlePubMedGoogle Scholar
- Zhou Y, Ferguson J, Chang JT, Kluger Y: Inter- and intra-combinatorial regulation by transcription factors and microRNAs. BMC Genomics. 2007, 8: 396- 10.1186/1471-2164-8-396PubMed CentralView ArticlePubMedGoogle Scholar
- Liang H, Li WH: MicroRNA regulation of human protein protein interaction network. Rna. 2007, 13 (9): 1402-1408. 10.1261/rna.634607PubMed CentralView ArticlePubMedGoogle Scholar
- Rhodes DR, Tomlins SA, Varambally S, Mahavisno V, Barrette T, Kalyana-Sundaram S, Ghosh D, Pandey A, Chinnaiyan AM: Probabilistic model of the human protein-protein interaction network. Nat Biotechnol. 2005, 23 (8): 951-959. 10.1038/nbt1103View ArticlePubMedGoogle Scholar
- Sethupathy P, Megraw M, Hatzigeorgiou AG: A guide through present computational approaches for the identification of mammalian microRNA targets. Nat Methods. 2006, 3 (11): 881-886. 10.1038/nmeth954View ArticlePubMedGoogle Scholar
- Baek D, Villen J, Shin C, Camargo FD, Gygi SP, Bartel DP: The impact of microRNAs on protein output. Nature. 2008, 455 (7209): 64-71. 10.1038/nature07242PubMed CentralView ArticlePubMedGoogle Scholar
- Xu J, Wong C: A computational screen for mouse signaling pathways targeted by microRNA clusters. Rna. 2008, 14 (7): 1276-1283. 10.1261/rna.997708PubMed CentralView ArticlePubMedGoogle Scholar
- Newman ME, Girvan M: Finding and evaluating community structure in networks. Phys Rev E Stat Nonlin Soft Matter Phys. 2004, 69 (2 Pt 2): 026113-View ArticlePubMedGoogle Scholar
- Adamcsek B, Palla G, Farkas IJ, Derenyi I, Vicsek T: CFinder: locating cliques and overlapping modules in biological networks. Bioinformatics. 2006, 22 (8): 1021-1023. 10.1093/bioinformatics/btl039View ArticlePubMedGoogle Scholar
- Petrocca F, Vecchione A, Croce CM: Emerging role of miR-106b-25/miR-17-92 clusters in the control of transforming growth factor beta signaling. Cancer Res. 2008, 68 (20): 8191-8194. 10.1158/0008-5472.CAN-08-1768View ArticlePubMedGoogle Scholar
- Yu Z, Wang C, Wang M, Li Z, Casimiro MC, Liu M, Wu K, Whittle J, Ju X, Hyslop T, et al.: A cyclin D1/microRNA 17/20 regulatory feedback loop in control of breast cancer cell proliferation. J Cell Biol. 2008, 182 (3): 509-517. 10.1083/jcb.200801079PubMed CentralView ArticlePubMedGoogle Scholar
- Xie W, Song C, Young NL, Sperling AS, Xu F, Sridharan R, Conway AE, Garcia BA, Plath K, Clark AT, et al.: Histone h3 lysine 56 acetylation is linked to the core transcriptional network in human embryonic stem cells. Mol Cell. 2009, 33 (4): 417-427. 10.1016/j.molcel.2009.02.004PubMed CentralView ArticlePubMedGoogle Scholar
- Hsu CW, Juan HF, Huang HC: Characterization of microRNA-regulated protein-protein interaction network. Proteomics. 2008, 8 (10): 1975-1979. 10.1002/pmic.200701004View ArticlePubMedGoogle Scholar
- Gusev Y, Schmittgen TD, Lerner M, Postier R, Brackett D: Computational analysis of biological functions and pathways collectively targeted by co-expressed microRNAs in cancer. BMC Bioinformatics. 2007, 8 (Suppl 7): S16- 10.1186/1471-2105-8-S7-S16PubMed CentralView ArticlePubMedGoogle Scholar
- Cui Q, Yu Z, Purisima EO, Wang E: Principles of microRNA regulation of a human cellular signaling network. Mol Syst Biol. 2006, 2: 46- 10.1038/msb4100089PubMed CentralView ArticlePubMedGoogle Scholar
- Cui Q, Yu Z, Pan Y, Purisima EO, Wang E: MicroRNAs preferentially target the genes with high transcriptional regulation complexity. Biochem Biophys Res Commun. 2007, 352 (3): 733-738. 10.1016/j.bbrc.2006.11.080View ArticlePubMedGoogle Scholar
- Peri S, Navarro JD, Amanchy R, Kristiansen TZ, Jonnalagadda CK, Surendranath V, Niranjan V, Muthusamy B, Gandhi TK, Gronborg M, et al.: Development of human protein reference database as an initial platform for approaching systems biology in humans. Genome Res. 2003, 13 (10): 2363-2371. 10.1101/gr.1680803PubMed CentralView ArticlePubMedGoogle Scholar
- Lee Y, Kim M, Han J, Yeom KH, Lee S, Baek SH, Kim VN: MicroRNA genes are transcribed by RNA polymerase II. Embo J. 2004, 23 (20): 4051-4060. 10.1038/sj.emboj.7600385PubMed CentralView ArticlePubMedGoogle Scholar
- Liang Y, Ridzon D, Wong L, Chen C: Characterization of microRNA expression profiles in normal human tissues. BMC Genomics. 2007, 8: 166- 10.1186/1471-2164-8-166PubMed CentralView ArticlePubMedGoogle Scholar
- Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, Sweet-Cordero A, Ebert BL, Mak RH, Ferrando AA, et al.: MicroRNA expression profiles classify human cancers. Nature. 2005, 435 (7043): 834-838. 10.1038/nature03702View ArticlePubMedGoogle Scholar
- Williams AE, Perry MM, Moschos SA, Lindsay MA: microRNA expression in the aging mouse lung. BMC Genomics. 2007, 8: 172- 10.1186/1471-2164-8-172PubMed CentralView ArticlePubMedGoogle Scholar
- Lee EJ, Baek M, Gusev Y, Brackett DJ, Nuovo GJ, Schmittgen TD: Systematic evaluation of microRNA processing patterns in tissues, cell lines, and tumors. Rna. 2008, 14 (1): 35-42. 10.1261/rna.804508PubMed CentralView ArticlePubMedGoogle Scholar