Regulatory coordination of clustered microRNAs based on microRNA-transcription factor regulatory network
© Wang et al; licensee BioMed Central Ltd. 2011
Received: 12 April 2011
Accepted: 16 December 2011
Published: 16 December 2011
MicroRNA (miRNA) is a class of small RNAs of ~22nt which play essential roles in many crucial biological processes and numerous human diseases at post-transcriptional level of gene expression. It has been revealed that miRNA genes tend to be clustered, and the miRNAs organized into one cluster are usually transcribed coordinately. This implies a coordinated regulation mode exerted by clustered miRNAs. However, how the clustered miRNAs coordinate their regulations on large scale gene expression is still unclear.
We constructed the miRNA-transcription factor regulatory network that contains the interactions between transcription factors (TFs), miRNAs and non-TF protein-coding genes, and made a genome-wide study on the regulatory coordination of clustered miRNAs. We found that there are two types of miRNA clusters, i.e. homo-clusters that contain miRNAs of the same family and hetero-clusters that contain miRNAs of various families. In general, the homo-clustered as well as the hetero-clustered miRNAs both exhibit coordinated regulation since the miRNAs belonging to one cluster tend to be involved in the same network module, which performs a relatively isolated biological function. However, the homo-clustered miRNAs show a direct regulatory coordination that is realized by one-step regulation (i.e. the direct regulation of the coordinated targets), whereas the hetero-clustered miRNAs show an indirect regulatory coordination that is realized by a regulation comprising at least three steps (e.g. the regulation on the coordinated targets by a miRNA through a sequential action of two TFs). The direct and indirect regulation target different categories of genes, the former predominantly regulating genes involved in emergent responses, the latter targeting genes that imply long-term effects.
The genomic clustering of miRNAs is closely related to the coordinated regulation in the gene regulatory network. The pattern of regulatory coordination is dependent on the composition of the miRNA cluster. The homo-clustered miRNAs mainly coordinate their regulation rapidly, while the hetero-clustered miRNAs exert control with a delay. The diverse pattern of regulatory coordination suggests distinct roles of the homo-clustered and the hetero-clustered miRNAs in biological processes.
MicroRNA (miRNA) is a class of ~22 nt small non-coding RNAs (ncRNAs) which inhibit gene expression at post-transcriptional stage by binding to the 3'UTR of mRNAs. They play essential roles in many crucial biological processes, including development, differentiation, apoptosis and cell proliferation [1–4], as well as numerous human diseases, such as chronic lymphocytic leukemia, fragile X syndrome, and various types of cancers [5–8].
The studies on the biogenesis of miRNAs [9, 10] show that miRNA is firstly transcribed as pri-miRNA (i.e. primary miRNA) in the nucleus, then exported to the cytoplasm after being cleaved into pre-miRNA (precursor miRNA). In cytoplasm, pre-miRNA is processed into mature miRNA and incorporated into the RNA-induced silencing complex (RISC) which subsequently binds to the 3'UTR of mRNAs. The conversion of pre-miRNA to mature miRNA is generally one-to-one (i.e. one pre-miRNA once generates one mature miRNA), although both strands of pre-miRNAs could potentially become mature miRNAs . However, the splicing of one pri-miRNA could result in multiple pre-miRNAs. This is because miRNA genes tend to be clustered in the genome  and one pri-miRNA could be a cluster of several miRNAs .
A set of miRNAs that are closely distributed in genome is termed as the miRNA cluster. The clustering propensity of miRNAs was first been discovered by large-scale surveys of small ncRNAs [14, 15]. At present, it has been confirmed that miRNA clusters are widely distributed in animal genomes [16, 17]. The conservation of miRNA clusters across species  indicates that miRNA clusters adapt special regulatory functions in biological processes. In addition, it has been shown by expression studies that the clustered miRNAs are often co-expressed [19–21], suggesting that they are jointly transcribed as a polycistron. Thus, the hypothesis comes up that the genomic coordination of clustered miRNA genes, which further leads to their coordinated transcription, will consequently result in a functional coordination. However, it is still unclear how the clustered miRNAs function coordinately.
Recently, Yuan X-Y et al. studied the functional coordination of clustered miRNAs based on the protein-protein interaction (PPI) network . They found that the clustered miRNAs tend to target the mRNA of proteins that are located in the same functional module. While this kind of correlation supports the view of functional coordination of clustered miRNAs, little is known about the underlying mechanisms. In addition, the PPI network which is composed of direct protein-protein interactions cannot provide the successive regulation details of miRNAs. For example, the successive regulation on a protein by a miRNA through one or more TFs is not included in the PPI network.
Here, we studied the regulatory coordination of clustered miRNAs based on the miRNA-Transcription factor (miRNA-TF) regulatory network that comprises the interactions between transcription factors (TF), miRNAs and non-TF protein-coding genes. We found that there are two types of miRNA clusters, i.e. homo-clusters (miRNA clusters composed of miRNAs from a single miRNA family) and hetero-clusters (miRNA clusters composed of miRNAs from multiple miRNA families). In general, both the homo- and the hetero-clustered miRNAs show the behavior of regulatory coordination. However, the ways of regulatory coordination of both types of clustered miRNAs are different. The homo-clustered miRNAs show a direct regulatory coordination which is realized by a single regulation step (i.e. direct regulation), and tend to be involved in emergency processes, whereas the hetero-clustered miRNAs show an indirect regulatory coordination which is accomplished by 3 or more steps, and tend to participate in more complex processes.
Homologous and heterologous miRNA clusters
Usually, the miRNAs of the same family target a similar set of mRNAs, since they have the same seed regions  which predominantly determine the targets of miRNA . Thus, it is suggested that the homo-clustered and the hetero-clustered miRNAs may have quite different regulatory features.
Clustered miRNAs in network modules
In the miRNA-TF regulatory network containing the interactions between TF, miRNA and non-TF protein-coding genes (see Methods), 39 modules that have dense interactions were found using MCODE plus  of the network tool Cytoscape . These modules altogether contain 47 homo-clustered miRNAs, 132 hetero-clustered miRNAs and 232 isolated miRNAs that are far from other miRNAs and do not belong to any clusters. The statistical analysis (see Methods) shows that the homo-clustered miRNAs and the hetero-clustered miRNAs are both significantly enriched (P < 0.01) in these 39 modules. Moreover, more than 50% of the modules (i.e. 22 out of 39) contain at least one miRNA pair that comes from the same miRNA cluster. The average cluster entropy (see Methods) of the 39 modules which describes the diversity of the miRNA clusters within a module is significantly lower (P < 0.01) than that in the random case in which miRNAs were randomly assigned to groups of the same number as human miRNA clusters (see Methods). These findings suggest that the miRNAs located within one genomic cluster tend to be involved preferentially in the same module. Such preference means that, in general, the clustered miRNAs exert a coordinated regulation since a module in a biological network usually represents a special or relatively independent biological function .
Distinctive regulatory coordination of homo- and hetero-clusters
where l i is the number of targets that are regulated by miRNA i, N is the total number of targets, and s ik is the element of the connecting matrix. s ik equals to 1 when miRNA i regulates target k directly (or indirectly), otherwise it equals to 0. Clearly, TO ij varies in the range of 0 to 1. The closer TO ij is to 1, the more targets miRNA i and miRNA j share, and the stronger is the coordinated regulation by miRNA i and miRNA j.
Functional analysis of targets of homo- and hetero-clusters
The function of miRNA targets was analyzed using the tool DAVID (http://david.abcc.ncifcrf.gov/, see Methods). We found that the targets of homo-clustered miRNAs are significantly involved in emergency processes that need to be preceded rapidly, such as response to stimuli and the intrinsic apoptotic pathway that involves mitochondria. The latter represents an emergency process since the intrinsic apoptotic pathway is usually activated to induce a rapid cell death . In contrast, the hetero-clustered miRNAs are involved in complex biological processes including the metabolism and the extrinsic apoptotic pathway that happens primarily in the cytoplasm. These processes are generally of less urgency, but they are more complex than those of homo-clusters. It has been revealed that thousands of genes and reactions are involved in metabolic processes , and the extrinsic apoptotic pathway is composed of several complex caspase processes .
Moreover, the extrinsic apoptotic process is also more complex than the internal apoptotic process since it comprises significantly more reactions . Besides, three functions (i.e. signal transduction, development and transport) are shared targets of homo- and hetero-clustered miRNAs.
MicroRNA clusters, which are groups of tandem miRNA genes that are closely located in the genome, are abundantly and widely distributed in animal genomes. It has been revealed that about 50% of the miRNA genes in Drosophila  and over 30% of the miRNA genes in human, mouse, rat and chicken are located in clusters . Co-expression experiments of clustered miRNA genes show that one miRNA cluster is usually transcribed as a single transcriptional unit . This suggests the existence of some kind of regulatory coordination between the clustered miRNAs. However, it is still unclear how the clustered miRNAs coordinate their regulation. Here, we describe a genome-wide study on the regulatory coordination of clustered miRNAs based on the miRNA-TF regulatory network.
The miRNA clusters are classified into homo-clusters that contain miRNAs of the same family and hetero-clusters that contain miRNAs of multiple families. Most of the miRNA clusters are either homo-clusters or hetero-clusters of miRNAs with completely different families. Such polarized behavior indicates that the family composition of clustered miRNAs may be an important characteristic that is closely related to the regulatory features of miRNA clusters.
In this study, we have used a miRNA-TF regulatory network that represents the regulation exerted by TFs and miRNAs on gene expression of target genes. This regulatory network presents a comprehensive view on the regulations of miRNAs since it involves the transcriptional regulation of miRNAs genes by TFs as well as the direct regulation of the targeted mRNAs by miRNAs. Such regulatory network has previously been used to study the combinatory regulation of miRNAs and TFs on gene expression. For example, different types of miRNA-TF co-regulations have been revealed based on the miRNA-TF regulatory network [31, 32]. In addition, it is reported by Kang et al. that there are two-layer regulations on the gene expression, where TFs function as important mediators of miRNA-initiated regulatory effects . These studies suggest that the miRNA-TF regulatory network is a good substrate for studying the complex regulatory features of miRNAs.
The result that the clustered miRNAs, whether they are from homo- or hetero-clusters, preferably exert their effect in one module suggests a general regulatory coordination of clustered miRNAs. Intuitively, the regulatory coordination of homo-clustered miRNAs is ascribed to the high sequence similarity of homologous miRNAs. However, not all the sequences of homologous miRNAs are similar enough to bind to the same targets. More than half of the homologous miRNA pairs share less than 50% of the targets (see Additional File 2). An example is miR-329 family, in which any pairs of the three members (i.e. hsa-mir-543, hsa-mir-329, hsa-mir-495) share less than 20% of the targets. The homo-clustered miRNAs are the homologous ones that share large amount of targets. This suggests that the homo-clustered miRNAs are not the arbitrary homologous miRNAs, but the ones finely designed for the regulatory coordination. The target overlaps of hetero-clustered miRNAs are much smaller than the homo-clustered miRNAs, but they similarly appear in the same modules. This indicates that the homo-clustered miRNAs and the hetero-clustered miRNAs have distinctive ways of coordinated regulation. In addition, there are some cases that not all the members of miRNAs in homo-/hetero-clusters are found in the same modules. One possible reason is that the size of modules depends on the parameter that scales the density of interactions. If the parameter is strict, the resulting modules, which are generally of small size, will include few miRNA cluster members. Another reason may be that 10 kb is not an accurate cutoff for the definition of miRNA clusters. The miRNA cluster members that are not found in the same modules may not be included anymore in the cluster when there is a little deviation of the cutoff.
It is clear that the regulatory coordination of homo-clustered miRNAs is achieved by one regulation step (i.e. at the level of direct targets of miRNAs), since the homologous miRNAs in homo-clusters have almost the same targets. However, the regulatory coordination of hetero-clustered miRNAs is realized by at least 3 steps. Thus, the regulatory coordination of homo-clustered miRNAs is direct, while that of hetero-clustered miRNAs is indirect. These two types of coordinated regulation both have their own advantage. Since there is no intermediate regulator between the miRNAs and their targets, the direct coordinated regulation has the advantage of accuracy and quickness. Whereas, the indirect coordinated regulation has the advantage of variety since diverse types of coordinated regulation can be realized by inducing additional cross-regulations between intermediated regulators (see Additional File 3).
Our study is focused on the internal coordination of genomic-clustered miRNAs. The results suggest that there are two types of miRNA clusters, i.e. homo-clusters that contain miRNAs of the same family and hetero-clusters that contain miRNAs of multiple families. These two types of miRNA clusters show distinct behaviors of regulatory coordination in the gene regulatory network that represents the direct interactions between miRNA, TF and non-TF protein-coding genes. The homo-clusters show a direct regulatory coordination and tend to be involved in biological processes of emergency situations, while the hetero-clusters show an indirect regulatory coordination and tend to take part in more complex biological processes.
Our study shows the diversity of miRNA regulations responding to the complex requirements of biological functions and contributes to understand the complex function and regulatory mechanism of miRNAs at a network scale.
All the human miRNAs along with their genomic information were retrieved from miRBase 13.0 . Two miRNAs that are consecutively located within 10 kb of each other were considered to belong to one miRNA cluster. This definition about miRNA clusters is based on the study of miRNA genomic distribution. In previous, it has been revealed that the distances between miRNA pairs located consecutively in genome are following a biomodal distribution. The valley between the two peaks is located at around 10 kb, suggesting that 10 kb may be the reasonable cutoff to define miRNA clusters . In total, there are 718 human precursor miRNAs (pre-miRNAs) of distinct genomic locations. Among these 718 pre-miRNAs, about 36% form 66 miRNA clusters. The number of miRNA clusters is stable as the cutoff deviates from 10 kb. It varies at the range of 64-67 as the cutoff increases from 5 kb to 15 kb. The miRNAs in the same miRNA cluster composed of intronic/extronic miRNAs are of the same host genes (see Additional File 4), suggesting that the miRNA clusters tend to be transcribed into one transcript. For the intergenic miRNA clusters, it is suggested that the intergenic miRNA cluster may have short transcript less than 4 kbp . This means our definition about intergenic miRNA clusters has the risk to draw an inaccurate conclusions since some of our intergenic miRNA clusters (about 30%) are much longer than 4 kbp. To detect whether our definition about intergenic miRNA clusters take a bias impact on the conclusion, we compared the coordinated behavior of our intergenic miRNA clusters with that of a newly-defined clusters which follow the constraint that the length of every intergenic miRNA cluster should be less than 4 kbp (i.e. the 4 kbp-constraint intergenic miRNA cluster). The results show that the direct and indirect coordinated behaviors of these two types of intergenic clusters are both very similar (see Additional File 5 and 6 for details). This suggests our definition on intergenic miRNA clusters is good enough for our study, and it does not take any bias impact on the final conclusion.
The information on miRNA families has been derived from Rfam database, in which various kinds of RNAs are classified into families based on their sequence and structural alignment . A miRNA family generally means a collection of miRNAs that are derived from a common ancestor.
Family entropy of a miRNA cluster
Given that a cluster comprises miRNAs of N fam families, p i is the probability that any miRNA in the cluster belongs to family i. E fam varies in the range of 0-1. E fam = 0 means that all the miRNAs in the cluster come from the same family, and E fam = 1 means that all clustered miRNAs come from different families.
Construction of the regulatory network
To study the regulatory feature of the clustered miRNAs, we constructed the miRNA-TF regulatory network, which contains the regulatory relationships between TFs, miRNAs and non-TF protein-coding genes for human (see Additional File 7). Generally, there are two kinds of interactions in the network which respectively start from TFs and miRNAs. An interaction starting from a TF means that the TF regulates the transcription of the target, while an interaction starting from a miRNA means that the miRNA represses the translation of the target. We predicted the interactions starting from a TF by searching the conserved TF binding sites (TFBSs) within a putative promoter area 1 kb upstream the transcriptional start site of the target. Firstly, all the potential human TFBSs are collected from TRANSFAC (version 2009.4)  based on the position weighted matrix (PWM). Secondly, the conserved TFBSs are derived from the conserved promoter area across the 5 species of human, mouse, dog, cow and opossum based on RefSeq annotation of UCSC hg18 http://genome.ucsc.edu/index.html. Finally, the conserved relationships between TFBSs and TFs are predicted using the Match-algorithm provided by TRANSFAC. Note that the prediction of the TFBSs is done based on phylogenetic footprinting from five mammalian genomes in order to limit the false positive of the TFBS prediction. Since it has been discussed that the miRNAs within the same cluster are often transcribed simultaneously, we took a whole cluster of miRNAs as a single unit, and searched the TFBSs in the 10 kb upstream to the start point of the first miRNA in the cluster. In addition, we predicted the interactions starting from miRNAs by the three tools of Targetscan , Pictar  and Tarbase . The union results predicted by these three tools are taken to give a comprehensive regulatory network so that many potential instants of regulatory coordination will be dug out to help further studies in experiment. Replacing the union results with the results of Tarbase which are verified by experiments does not affect our conclusion.
Statistic analysis on the clustered miRNAs in modules
Enrichment of homo-/hetero- clusters
The enrichment of homo-/hetero-clusters in the 39 modules is evaluated based on the random distribution of the average homo-/hetero-clustered miRNA numbers in modules. Firstly, 1000 random cases are generated by randomly re-participating miRNAs into the 39 modules with the total miRNA number in each module kept. Then the random distribution is obtained by calculating the average homo-/hetero-clustered miRNA numbers in modules. Finally, the significance of the enrichment of homo-/hetero-clustered miRNAs in modules is evaluated based on the random distribution.
miRNA cluster entropy of a module
Note that p c i is the probability of miRNA cluster i occurring in a module and N c is the total number of miRNA clusters in the module. Similar as E fam , E c varies at the range of 0-1. E c = 0 means that all the miRNAs in the module are in a same cluster, and E c = 1 means that the miRNAs are from completely different clusters.
Random clusters of miRNAs
There are totally 458 isolated miRNAs and 66 miRNA clusters in genome. If the isolated miRNAs are considered as the pseudo-cluster, there are 524 (i.e. 458 + 66) clusters of miRNAs. Thus, a set of random clusters are generated by re-assigning miRNAs into these 524 clusters with the number of miRNAs in each cluster maintained. All the distributions for the random miRNA clusters are the ones averaged over 1000 random sets.
Functional analysis of miRNA targets
Enriched functions of miRNA targets.
Response to stimuli
Internal apoptotic pathway
External apoptotic pathway
This work was supported by the National Natural Science Foundation of China (Grant no. 30890044), the National Basic Research Program of China (Grant no. 2007CB814806), and the German-China joint project (Grant no. CHN08/031).
- He L, Hannon GJ: MicroRNAs: small RNAs with a big role in gene regulation. Nat Rev Genet. 2004, 5: 522-531. 10.1038/nrg1379.View ArticlePubMedGoogle Scholar
- Harfe BD: MicroRNAs in vertebrate development. Curr Opin Genet Dev. 2005, 15: 410-415. 10.1016/j.gde.2005.06.012.View ArticlePubMedGoogle Scholar
- Bushati N, Cohen SM: microRNA functions. Annu Rev Cell Dev Biol. 2007, 23: 175-205. 10.1146/annurev.cellbio.23.090506.123406.View ArticlePubMedGoogle Scholar
- Krol J, Loedige I, Filipowicz W: The widespread regulation of microRNA biogenesis, function and decay. Nat Rev Genet. 2010, 11: 597-610.PubMedGoogle Scholar
- Calin GA, Ferracin M, Cimmino A, Di Leva G, Shimizu M, Wojcik SE, Iorio MV, Visone R, Sever NI, Fabbri M, Luliano R, Palumbo T, Pichiorri F, Roldo C, Garzon R, Sevignani C, Rassenti L, Alder H, Volinia S, Liu CG, Kipps TJ, Negrini M, Croce CM: A MicroRNA signature associated with prognosis and progression in chronic lymphocytic leukemia. N Engl J Med. 2005, 353: 1793-1801. 10.1056/NEJMoa050995.View ArticlePubMedGoogle Scholar
- Alvarez-Garcia I, Miska EA: MicroRNA functions in animal development and human disease. Development. 2005, 132: 4653-4662. 10.1242/dev.02073.View ArticlePubMedGoogle Scholar
- Mitchell PS, Parkin RK, Kroh EM, Fritz BR, Wyman SK, Pogosova-Agadjanyan EL, Peterson A, Noteboom J, O'Briant KC, Allen A, Lin DW, Urban N, Drescher CW, Knudsen BS, Stirewalt DL, Gentleman R, Vessella RL, Nelson PS, Martin DB, Tewari M: Circulating microRNAs as stable blood-based markers for cancer detection. Proc Natl Acad Sci USA. 2008, 105: 10513-10518. 10.1073/pnas.0804549105.PubMed CentralView ArticlePubMedGoogle Scholar
- Beezhold KJ, Castranova V, Chen F: Microprocessor of microRNAs: regulation and potential for therapeutic intervention. Mol Cancer. 2010, 9: 134-PubMed CentralView ArticlePubMedGoogle Scholar
- Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116: 281-297. 10.1016/S0092-8674(04)00045-5.View ArticlePubMedGoogle Scholar
- Kim VN: MicroRNA biogenesis: coordinated cropping and dicing. Nat Rev Mol Cell Biol. 2005, 6: 376-385.View ArticlePubMedGoogle Scholar
- Ro S, Park C, Young D, Sanders KM, Yan W: Tissue-dependent paired expression of miRNAs. Nucleic Acids Res. 2007, 35: 5944-5953. 10.1093/nar/gkm641.PubMed CentralView ArticlePubMedGoogle Scholar
- Kim VN, Nam JW: Genomics of microRNA. Trends Genet. 2006, 22: 165-173. 10.1016/j.tig.2006.01.003.View ArticlePubMedGoogle Scholar
- Rodriguez A, Griffiths-Jones S, Ashurst JL, Bradley A: Identification of mammalian microRNA host genes and transcription units. Genome Res. 2004, 14: 1902-1910. 10.1101/gr.2722704.PubMed CentralView ArticlePubMedGoogle Scholar
- Lagos-Quintana M, Rauhut R, Meyer J, Borkhardt A, Tuschl T: New microRNAs from mouse and human. RNA. 2003, 9: 175-179. 10.1261/rna.2146903.PubMed CentralView ArticlePubMedGoogle Scholar
- Lai EC, Tomancak P, Williams RW, Rubin GM: Computational identification of Drosophila microRNA genes. Genome Biol. 2003, 4: R42-10.1186/gb-2003-4-7-r42.PubMed CentralView ArticlePubMedGoogle Scholar
- Thatcher EJ, Bond J, Paydar I, Patton JG: Genomic organization of zebrafish microRNAs. BMC Genomics. 2008, 9: 253-10.1186/1471-2164-9-253.PubMed CentralView ArticlePubMedGoogle Scholar
- Olena AF, Patton JG: Genomic organization of microRNAs. J Cell Physiol. 2010, 222: 540-545.PubMed CentralPubMedGoogle 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: 2697-2706. 10.1093/nar/gki567.PubMed CentralView ArticlePubMedGoogle Scholar
- Ambros V: The functions of animal microRNAs. Nature. 2004, 431: 350-355. 10.1038/nature02871.View ArticlePubMedGoogle Scholar
- Cullen BR: Transcription and processing of human microRNA precursors. Mol Cell. 2004, 16: 861-865.View ArticlePubMedGoogle Scholar
- Baskerville S, Bartel DP: Microarray profiling of microRNAs reveals frequent coexpression with neighboring miRNAs and host genes. RNA. 2005, 11: 241-247. 10.1261/rna.7240905.PubMed CentralView ArticlePubMedGoogle Scholar
- Yuan X, Liu C, Yang P, He S, Liao Q, Kang S, Zhao Y: Clustered microRNAs' coordination in regulating protein-protein interaction network. BMC Syst Biol. 2009, 3: 65-10.1186/1752-0509-3-65.PubMed 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, 34: D140-144. 10.1093/nar/gkj112.PubMed CentralView ArticlePubMedGoogle Scholar
- Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB: Prediction of mammalian microRNA targets. Cell. 2003, 115: 787-798. 10.1016/S0092-8674(03)01018-3.View ArticlePubMedGoogle Scholar
- Bader GD, Hogue CW: An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinformatics. 2003, 4: 2-10.1186/1471-2105-4-2.PubMed CentralView ArticlePubMedGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13: 2498-2504. 10.1101/gr.1239303.PubMed CentralView ArticlePubMedGoogle Scholar
- Rives AW, Galitski T: Modular organization of cellular networks. Proc Natl Acad Sci USA. 2003, 100: 1128-1133. 10.1073/pnas.0237338100.PubMed CentralView ArticlePubMedGoogle Scholar
- Jin Z, El-Deiry WS: Overview of cell death signaling pathways. Cancer Biol Ther. 2005, 4: 139-163. 10.4161/cbt.4.2.1508.View ArticlePubMedGoogle Scholar
- Ravasz E, Somera AL, Mongru DA, Oltvai ZN, Barabasi AL: Hierarchical organization of modularity in metabolic networks. Science. 2002, 297: 1551-1555. 10.1126/science.1073374.View ArticlePubMedGoogle Scholar
- Megraw M, Sethupathy P, Corda B, Hatzigeorgiou AG: miRGen: a database for the study of animal microRNA genomic organization and function. Nucleic Acids Res. 2007, 35: D149-155. 10.1093/nar/gkl904.PubMed CentralView ArticlePubMedGoogle 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: e131-10.1371/journal.pcbi.0030131.PubMed 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-396.PubMed CentralView ArticlePubMedGoogle Scholar
- Tu K, Yu H, Hua YJ, Li YY, Liu L, Xie L, Li YX: Combinatorial network of primary and secondary microRNA-driven regulatory mechanisms. Nucleic Acids Res. 2009, 37: 5969-5980. 10.1093/nar/gkp638.PubMed CentralView ArticlePubMedGoogle Scholar
- Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ: miRBase: tools for microRNA genomics. Nucleic Acids Res. 2008, 36: D154-158. 10.1093/nar/gkn221.PubMed CentralView ArticlePubMedGoogle Scholar
- Saini HK, Griffiths-Jones S, Enright AJ: Genomic analysis of human microRNA transcripts. Proc Natl Acad Sci USA. 2007, 104: 17719-17724. 10.1073/pnas.0703890104.PubMed CentralView ArticlePubMedGoogle Scholar
- Griffiths-Jones S, Bateman A, Marshall M, Khanna A, Eddy SR: Rfam: an RNA family database. Nucleic Acids Res. 2003, 31: 439-441. 10.1093/nar/gkg006.PubMed CentralView ArticlePubMedGoogle Scholar
- Matys V, Fricke E, Geffers R, Gossling E, Haubrock M, Hehl R, Hornischer K, Karas D, Kel AE, Kel-Margoulis OV, Kloos DU, Land S, Lewicki-Potapov B, Michael H, Munch R, Reuter I, Rotert S, Saxel H, Scheer M, Thiele S, Wingender E: TRANSFAC: transcriptional regulation, from patterns to profiles. Nucleic Acids Res. 2003, 31: 374-378. 10.1093/nar/gkg108.PubMed CentralView ArticlePubMedGoogle Scholar
- Krek A, Grun D, Poy MN, Wolf R, Rosenberg L, Epstein EJ, MacMenamin P, da Piedade I, Gunsalus KC, Stoffel M, Rajewsky N: Combinatorial microRNA target predictions. Nat Genet. 2005, 37: 495-500. 10.1038/ng1536.View ArticlePubMedGoogle Scholar
- Sethupathy P, Corda B, Hatzigeorgiou AG: TarBase: A comprehensive database of experimentally supported animal microRNA targets. RNA. 2006, 12: 192-197.PubMed CentralView ArticlePubMedGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29. 10.1038/75556.PubMed CentralView ArticlePubMedGoogle Scholar
- Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4: P3-10.1186/gb-2003-4-5-p3.View ArticlePubMedGoogle Scholar