Skip to main content


You are viewing the new article page. Let us know what you think. Return to old version

Dissecting microregulation of a master regulatory network



The master regulator p53 tumor-suppressor protein through coordination of several downstream target genes and upstream transcription factors controls many pathways important for tumor suppression. While it has been reported that some of the p53's functions are microRNA-mediated, it is not known as to how many other microRNAs might contribute to the p53-mediated tumorigenesis.


Here, we use bioinformatics-based integrative approach to identify and prioritize putative p53-regulated miRNAs, and unravel the miRNA-based microregulation of the p53 master regulatory network. Specifically, we identify putative microRNA regulators of a) transcription factors that are upstream or downstream to p53 and b) p53 interactants. The putative p53-miRs and their targets are prioritized using current knowledge of cancer biology and literature-reported cancer-miRNAs.


Our predicted p53-miRNA-gene networks strongly suggest that coordinated transcriptional and p53-miR mediated networks could be integral to tumorigenesis and the underlying processes and pathways.


The p53 tumor suppressor gene has been implicated as a master regulator of genomic stability, cell cycle, apoptosis, and DNA repair [1, 2]. p53 is known to act as both transcriptional activator and repressor of expression of specific genes [3]. Studies of p53-mediated repression have shown that both genes that modulate apoptotic responses and genes that promote cell cycle progression can be repressed by p53 [4].

MicroRNAs (miRNAs or miRs) are small non-coding RNAs typically of 21–25 nucleotides in length and reported to be involved in the regulation of a variety of biological processes, including development, cell death, cell proliferation, hematopoiesis and nervous systems patterning [5]. miRNAs enforce posttranscriptional silencing through the RNA interference pathway [6]. Currently more than 450 human miRNAs have been identified and deposited in the miRBase [7]. Computational prediction of miRNA targets reveal that each human miRNA can potentially target several genes [8, 9], underscoring the importance of these tiny noncoding RNAs in gene regulation.

Several previous studies have identified ten positive or negative feedback loops in the p53 pathway (reviewed in [10]). All of these networks or circuits are autoregulatory in that they are either induced by p53 activity at the transcriptional level, transcriptionally repressed by p53 or are regulated by p53-induced proteins [10]. Recently, five independent studies have shown that p53 upregulates miR-34a in response to DNA damage leading to cell cycle arrest and apoptosis [1115]. In light of these reports and based on the growing evidence that microRNAs themselves act to either promote cancer or to stop cancer spread [1619], it is reasonable to conceive that the tumor suppressor p53 might be involved in this microRNA-related network in cancer cells. In support of this, two miRNAs (miR-372 and miR-373) were found to function as oncogenes in human testicular germ cell tumors, probably by numbing the p53 pathway and thereby allowing tumorigenesis [20]. Yang et al recently suggested that analysis of microRNA regulation in the transcriptional network of the p53 gene can reveal relationships between oncogenic microRNAs and p53 [21]. We therefore sought to explore if there are other microRNAs that are components of the p53 regulatory network. In an earlier study (Kaimal and Sinha et al., unpublished), scanning the 474 human miRNAs for putative p53 sites, we identified 143 miRNAs that have at least one p53 site (using p53MH algorithm [22]) and are predicted to target at least one known gene. We termed these 143 miRNAs as p53-miRs (Kaimal and Sinha et al., unpublished). In the current study, apart from correlating these putative p53-miRs with the upstream regulators and downstream targets of p53, we prioritize our predictions using (a) literature reported differentially expressed miRNAs in cancer; (b) miRNAs reported as induced or repressed following p53 activation [11]; and (c) miRNA targets that are functionally enriched with biological processes like apoptosis, cell cycle.

Results and discussion

Prediction of p53-miRsand the target genes

A significant fraction of human intergenic miRNAs are reported to possess TSS sites within 2 kb of the pre-miRNA [23]. Although it is accepted that intronic miRNAs are generally transcribed along with their host genes [24], TSS predictions have shown them to lie predominantly in the region from -2 kb to -6 kb [23]. Hence, we limited our search space to 10 kb flanking regions of the miRNAs. Determination of exact pri-miRNA transcript lengths however is possible only through experimental work.

We used p53MH [22], an algorithm that identifies p53-responsive genes in the human and mouse genome, to identify putative p53-binding site within miRNAs (see Methods). Out of 474 human miRNAs analyzed, 227 had at least one putative p53 site (and a total of 319 p53 sites; some miRNAs have multiple p53 sites) within their flanking regions of 10 kb (Additional File 1). Of the 319 p53 sites, 5 sites were overlapping with exonic region and we did not consider these for further analysis. The 314 p53 sites (Additional File 2) occur within 10 kb flanking sequences of 223 miRNAs. We identified two p53 binding sites in mir-34a; however neither of these corresponds to the experimental validated p53 site [12] in mir-34a. This is because in the current analysis we limit the miRNA flanking sequence search space to 10 kb only while the experimentally validated p53 site in mir-34a occurs at approximately 30 kb upstream to mir-34a [12].

To predict target genes of the 223 miRNAs, we used MAMI server [25], which has compilation of five [8, 9, 2628] target prediction algorithms. A recent study showed that about 82% of unique miRNA-target pairs are predicted by only a single algorithm [19]. Therefore, a combination of predictions from all algorithms might provide a much more comprehensive list of putative miRNA-target pairs than do any single prediction. 154/223 miRNAs are predicted to target at least one gene while 69/223 have no predicted targets. The 154 miRNAs mapped to 143 miRNAs. This is because some of the miRNAs have more than one predicted hairpin precursor sequences. For example mir-135a has two predicted hairpin precursors (mir-135a-1 on chromosome 3 and mir-135a-2 on chromosome 12). The targets are however based on mir-135a only. Thus, our final set consists of 143 miRNAs; each of which has at least one p53 site and is predicted to target at least one known gene. We call these 143 miRNAs as p53-miRs. There were 12497 genes predicted as targets of these 143 p53-miRs! Assuming that most of these predicted targets are purely speculative (false positive rate of at least 4 of these algorithms was estimated to be 20–30% [29]), we sought to integrate multiple types of data resources and explore as to what number of p53-miRs could be targeting cancer-related genes.

MicroRNAs associated with p53 upstream and downstream genes

Transcriptional networks commonly contain positive- and negative-feedback loops, which provide robustness and fine-tuning to gene programs [30, 31]. Recent studies have suggested a possible coordinated transcriptional and miRNA-mediated regulation as a recurrent motif to enhance the robustness of gene regulation in mammalian genomes [3235]. It was also reported that miRNAs predominantly target positive regulatory motifs, highly connected scaffolds and most downstream network components such as signaling transcription factors, but less frequently target negative regulatory motifs, common components of basic cellular machines and most upstream network components such as ligands [32]. Using the p53 master regulatory network as a case study, we investigated the potential crosstalk between the miRNAs and the p53 transcriptional network itself involving transcription factors upstream and downstream to p53. At least four possible scenarios can be envisaged for the upstream regulators (Figures 1A–D) and downstream targets (Figures 1E–H) of p53.

Figure 1

Four possible scenarios for the upstream regulators (A-D) and downstream targets (E-H) of p53. For the downstream coherent models, the p53 action on downstream target gene and miRNA is opposite (F and H). In case of downstream incoherent models (E and G), the p53 action on downstream target gene and miRNA is same. There may be a sequential gap in the activation time of the target gene and miRNA facilitating maintenance of steady states.

We downloaded the known upstream transcription factor (TF) regulators and downstream target TFs of p53 from the p53 KnowledgeBase [36]. For the 23 upstream regulators and 48 downstream TFs of p53, we then extracted the putative miRNA regulators using MAMI server [25], which has compilation of five target prediction algorithms [8, 9, 2628]. Since each of these five miRNA target prediction approaches generates an unpredictable number of false positives (an estimate puts it at 20–30% [29]), we intersected the results to identify the genes commonly predicted by at least two of the five algorithms [8, 9, 2628] (see Table 1; for a complete list, see Additional file 3). Although it is important to confirm that a miRNA is not a potential regulatory miRNA of a specific gene by using different bioinformatics algorithms [37], there are several experimentally validated microRNA target genes that are not predicted by more than one of the current prediction algorithms. For instance the experimentally validated mir-21-PTEN miRNA-target pair [38] is not predicted by any of the five algorithms (listed in Methods and used in MAMI compilation). To overcome this, we also checked the TarBase [39] which has experimentally validated microRNA targets irrespective of the computational predictions.

Table 1 p53-miRs (microRNAs which have a conserved p53 binding site in their flanking 10 kb genomic sequence) which putatively target known transcription factors that are either upstream or downstream to p53 network (based on p53 knowledge base [36]) (A: activator; R: Repressor; U: uncharacterized). Only those targets are shown that are predicted as miRNA targets by 2 or more than two algorithms. miRNAs that are in bold represent miRNAs induced more than 1.5 fold (p < 0.05) after p53 activation [11]. The last column indicates whether a miRNA has been reported in literature as differentially expressed in cancer or cancer cell lines.

Based on MAMI predictions, a total of 242 different miRNAs suppress p53 regulators or targets. Of these, 115 (~48%) are p53-miRs (Figure 2A). Some of the p53-miRs (32/143; Figure 2B) are found to exclusively target p53 upstream activators (13/143 p53-miRs; Figure 3A) or repressors (5/143; Figure 3B) or downstream activators (9/143; Figure 3C) or repressors (5/143; Figure 3D). Interestingly, many of these miRNAs have been reported as differentially expressed in various cancerous tissues or cell lines (see last column in Table 1).

Figure 2

Heat map representation of p53-miRs and their target genes relative to the p53 network. 115 p53-miRs are predicted to target transcription factors that are either p53 regulators (upstream) or p53 targets (downstream) (A).

Figure 3

Heat map representation of p53-miRs and their target genes relative to the p53 network. Panels A, B, C and D depict p53-miRs exclusively targeting p53 upstream activators, upstream repressors, downstream activators and downstream repressors respectively along with their predicted target transcription factors.

miRNAs upstream and downstream to p53: coherent and incoherent circuits

If there is a common upstream TF that regulates both the miRNA and its target gene, ideally, the transcription of the miRNAs and their targets should be oppositely regulated by this common upstream TF. For example, p53 which is known to repress the transcription of IER3 [40] may simultaneously activate the transcription of mir-142 (predicted to target IER3) that may further inhibit IER3 posttranscriptionally. This is a coherent model for all downstream targets that are repressed by p53. Previous genome-scale studies have in fact shown that predicted target transcripts of several tissue-specific miRNAs tend to be expressed at a lower level in tissues where the miRNAs are expressed [41, 42].

Alternately, a common TF (p53 in this case) can activate a downstream target gene and also a miRNA which then suppresses the downstream target gene. A well known example is the experimentally confirmed microregulatory network comprising miR-17-5p and E2F1 – both of which are transcriptionally activated by c-Myc in human cells [43]. We report a similar example of mir-122a, a putative p53-miR suppressing CCNG1, a downstream target of p53. CCNG1 is a known transcriptional target (activator) of the p53 tumor suppressor protein [44] and it was recently reported that mir-122a (predicted as p53-miR in our analysis) has an inverse correlation with cyclin G1 expression in primary liver carcinomas [45]. In other words, p53 activates CCNG1 and also mir-122a which in turn suppresses CCNG1. This is an incoherent model. Adding to this complexity further, it is also reported that CCNG1 negatively regulates the stabilization of p53 in a possible negative feedback loop [46] (Figure 4A). Thus, the regulatory network wherein p53 activates a downstream target and a miR (which target the downstream target) simultaneously appears "inefficient". However, feedforward loops have the potential to provide temporal control, because expression of the ultimate target may depend on the accumulation of adequate levels of master regulator and the secondary regulators [47]. Therefore, if we consider a sequential gap in the activation time of the target genes and the miRNA, then this downstream p53-gene and p53-miR regulation appears coherent. Feedforward loops may provide a form of multistep ultrasensitivity [48], as any deviation from the p53's steady state would drive the downstream targets and miRNA (p53-miR) away from their steady states in the same direction. p53-miRs could therefore tune the production rate of the downstream p53 target gene opposite to the direction of p53's fluctuation. Such noise buffering probably helps to maintain target protein homeostasis and ensures more uniform expression [49] of the p53 downstream target genes within a cell population. Additionally, since the level of p53-miR defines the p53 downstream target's translation rate, their coexpression may allow p53-miR to fine-tune the downstream target's steady state. Thus, p53-miRs acting on p53 downstream targets could significantly shorten the response delay, leading to more effective noise buffering, as well as precise definition and maintenance of steady states.

Figure 4

Schematic representation of p53-miR-mediated models downstream (A and C) and upstream (B and D) to p53.

We found two entries in TarBase (database of experimentally validated microRNA target genes) [39] which are relevant and could be interesting to p53 network (a) mir-127 targeting BCL6 [50], an upstream repressor of p53 [51] (Figure 4B) and (b) mir-106a reported to suppress RB1 [52]. A previous study has shown that mir-106a is induced more than 1.5 fold (p < 0.05) after p53 activation [11]. p53 is known to suppress RB1 transcription through inhibition of the basal promoter activity [53]. Additionally, posttranscriptionally, p53 might suppress RB1 further by inducing mir-106a, a known suppressor of RB1 (Figure 4C).

As a master transcriptional regulator, p53 is also known to trans-repress genes including those affecting signaling pathways such as cell proliferation, apoptosis and cytoskeleton organization [54, 55]. Although, the mechanism(s) of repression is often unclear and may not be dependent on site-specific DNA-binding activity [5457], the responsiveness of some genes to p53 can be indirect in that target genes (transactivated or transrepressed) might be other transcription factors [55, 58] or microRNAs, which then modulate downstream responses to p53 signaling. We, therefore, hypothesize that in case of p53 downstream targets, the foremost activated positive feedback in a coupled feedback circuit can rapidly induce the "on" state transition of the signaling system (p53 activating downstream target genes), and that then another delayed positive feedback (p53 suppressing miRNAs that target p53 downstream target genes) robustly maintains this "on" state. Finally, the most delayed negative feedback reinstates the system in the original "off " state, preventing any further excessive response specific to the applied stimulus.

Loss of p53 expression in some tumors has also been shown to occur through inhibition of transcription of the p53 gene itself [51, 59, 60]. For example, mir-31, one of the p53-miRs predicted to suppress the upstream activator (YY1) of p53, is reported to be upregulated in colorectal cancer [61]. Shalgi et al recently reported that the extent of miR mediated regulation varies extensively among different genes, some of which, especially those who serve as regulators themselves, are subject to enhanced miRNA-based silencing [34].

p53-miRtarget genes – functional over-representation analysis

Functional over-representation analysis was performed to objectively identify biological processes potentially affected by p53-miR target genes. Specifically, the percentage of p53-miR target genes with a given gene ontology (GO) annotation was compared with the percentage of whole miRNA-target genes genome-wide with the same annotation. A significant p value (p < 0.05) indicates that the observed percentage of p53-miR target genes with a given annotation could not likely occur by chance given the frequency of whole miRNA-target genes genome-wide with the same annotation. MAMI predicted target genes for the 143 p53-miRs were 12497. We used DAVID [62], and our in-house developed ToppGene server [63] for the functional enrichment analysis. Since apoptosis and growth arrest are common known consequences of p53 activation, we also tested whether p53-miRs tend to target apoptotic and cell proliferation related genes. Interestingly, the enrichment analysis revealed that while cell cycle regulation (p = 8.41E-4) was among the top in the p53-miR targets functional enrichment (p < 0.05), apoptosis was not significant. The other significantly enriched GO terms include transcription, vesicle-mediated transport, intracellular signaling cascade and cellular metabolism (Additional File 4). Interestingly, vesicle-mediated transport (p = 3.9 E-6) and endocytosis (p = 5.1 E-4) were specifically enriched in cancer-associated p53-miR target genes. Abnormal expression or mutation of endocytosis proteins has been reported in human cancers [6466]. The possibility of harnessing the endocytic machinery to cancer through miRNAs could offer viable and promising leads in cancer therapy.

We extended p53-miR target gene functional enrichment analysis further to identify specific p53-miRs that significantly target apoptotic and cell proliferation genes (see Tables 2 and 3 for a list of top 25 p53-miRs and their target apoptotic and cell proliferation genes respectively). There was an overlap of 10 p53-miRs (mir-181c, mir-200a, mir-200b, mir-19b, mir-101, mir-106a, mir-142-3p, mir-195, mir-124a and mir-24) between these two groups. Earlier studies have shown that mir-23 is indeed critical for cell proliferation [67]. A recent study also implicated mir-23a, mir-24 and mir-195 in cardiac hypertrophy and reported that these miRNAs were regulated in response to stress signaling in the heart [68]. The Additional File 4 enlists the other functional categories significantly enriched within the putative p53-miR targets. It is to be noted that 18/25 (72%) apoptotic p53-miRs and 22/25 (88%) cell proliferation p53-miRs have been reported in literature as differentially regulated in various cancerous tissues or cancer cell lines (see last columns in Tables 2 and 3).

Table 2 Top 25 most significant p53-miRs that target known apoptotic genes (n = 685). 18/25 apoptotic p53-miRs have been reported in literature as differentially regulated in various cancerous tissues or cancer cell lines (see last column). The p-values represent the probability of a p53-miR targeting an apoptotic gene purely by chance.
Table 3 Top 25 most significant p53-miRs that target known cell proliferation genes (n = 712). 22/25 of these have been reported in literature as differentially regulated in various cancerous tissues or cancer cell lines (see last column). The p-values represent the probability of a p53-miR targeting a cell proliferation gene purely by chance.

Another interesting model is SIRT1-mir-34a-p53; SIRT1 is known to inhibit p53-mediated apoptosis by attenuation of the transcriptional activation potential of p53 primarily through deacetylation of p53 [69]. Several recent studies have shown that p53 activates mir-34a. A speculative model (Figure 4D) therefore is, in normal state mir-34a, a known p53-miR, suppresses SIRT1 (which is known to suppress p53) balancing the SIRT1-mediated transcriptional repression of p53. However, during tumorigenesis, increased levels of SIRT1 will suppress p53 triggering further suppression of p53 downstream mir-34a in a feedback loop. Examples of such models may have potential implications for cancer therapy. For example, combining DNA damage drugs (since p53 is strongly activated in response to DNA damage), SIRT1-mediated deacetylase inhibitors, and mir-34a activation (with "agomiR") may have synergistic effects in cancer therapy for maximally activating p53.

p53-miRand p53 protein interaction networks

Since the functional state of a protein-protein interaction network depends on gene expression, a fundamental question is what relationships exist between protein interaction network and gene regulation [70]. Liang and Li in their recent study of evidence for global correlation between microRNA repression and protein-protein interactions reported that interacting proteins tend to share more microRNA target-site types than random pairs [70]. In the current study, using the known p53-interactants, we calculated the probability of p53-miRs regulating the p53 interactants. We mined the BioGRID database [71], and currently there are 141 known interactants of p53. Surprisingly, 114 of these 141 p53-interactants were predicted to be regulated by p53-miRs and were shown to be statistically significant (p = 1.41e-21) (see Additional File 5 for details). Using the known data of miRNAs induced or repressed following p53 activation [11], we further filtered the p53-miRs that could be associated with the p53 interactome. Figures 5 and 6 show the p53 interactome along with the putative p53-miR regulators that are shown to be induced or repressed [11] respectively following the activation of p53.

Figure 5

p53 interactome along with the putative miR regulators that are shown to be induced following the activation of p53 [11]. miRNAs are represented by yellow ellipses, target genes are represented by green boxes. p53 is represented in the center as blue circle. Induction of the miRNAs by p53 is represented by directed red lines. Negative regulation of the target genes by miRNAs are represented by dark grey lines. The protein interactions are represented by undirected light blue lines. Networks (using force-directed layouts) were generated using aiSee [76] network visualization software.

Figure 6

p53 interactome along with the putative miR regulators that are shown to be repressed following the activation of p53 [11]. miRNAs are represented by yellow ellipses, target genes are represented by green boxes. p53 is represented in the center as blue circle. Induction of the miRNAs by p53 is represented by directed red lines. Negative regulation of the target genes by miRNAs are represented by dark grey lines. The protein interactions are represented by undirected light blue lines. Networks (using force-directed layouts) were generated using aiSee [76] network visualization software.


In conclusion, our bioinformatic-based analysis yielded several putative regulatory circuits that involve p53 and microRNAs along with their target genes. Although these predicted networks need to be experimentally validated to ascertain a bona fide regulatory relationship, it is still exciting to speculate that these p53-miRs may be critical to tumorigenesis through their regulation of the p53 master regulatory network. We strongly believe that our results substantially expand the hitherto little known repertoire of miRNAs in p53 regulatory network. It can be anticipated that further characterization of our predicted p53-miRs and their target genes will improve our understanding of the underlying processes of cancer, tumorigenesis, and ultimately impact response to therapy.


Human miRNA genomic analysis

The human miRNAs analyzed in this study (474 in total) were derived from the miRNA registry release 9.2 (May 2007; Genome Assembly NCBI-36 [72]). The genomic sequences, including 10 kb flanking sequences 5' and 3' for each of the miRNAs, were downloaded using Table Browser from the UCSC Golden Path [73].

Out of a total 474 human miRNAs analyzed, 100 microRNAs were part of a cluster (2 or more miRNAs occurring together on the genome). In all, there were 100 miRNAs organized into 37 clusters (see Additional Files 6 and 7 for details). The cluster definitions were based on those defined by Yu et al [74]. The nearest neighboring gene information (gene name, distance from the nearest miRNA, orientation with respect to the host or neighboring gene) for each of the miRNA was extracted from the public database RefFlat provided by the UCSC Golden Path [73].

Screening miRNA genomic sequences for p53 binding sites

The genomic sequences (10 kb flanking region, 5' and 3') of the 474 human miRNAs were scanned for potential p53-binding motifs using the p53MH algorithm [22]. Briefly, p53-binds to a specific sequence comprised of two decamers (RRRCWWGYYY) separated by 0–13 bp [75]. Nearly all established human response element sequences (REs) have at least one mismatch in this degenerate sequence [55]. The computer algorithm p53MH identifies potential p53-binding sites by a scoring system based on the percentage of similarity to a consensus sequence. We set the parameter conditions to include the following restrictions: (i) a p53MH score ≥ 85%, and (ii) the spacer between the 10-bp motifs is ≤ 4 bp. Although p53 is known to bind to sites with a spacer as high as 13 bp [75], for the current analysis we considered 4 bp spacer as the maximum limit because majority of the validated p53 binding sites are reported to have a spacer less than 4 bp (see Additional Files 1 and 2 for the putative p53-miR distribution).

p53 upstream and downstream activators and repressors

The group of TFs that act upstream to p53 either as activators or repressors and the p53 target genes involved in specific p53 downstream effects (activator or repressor) were downloaded from the compiled and curated dataset of the p53 knowledge-base [36]. We also mined the BioGRID database [71] for all known p53 interactants.


We used MAMI [25] to download the predicted targets for the p53-miRs. MAMI has a compilation of target prediction tools, DIANA-microT (version 5/05) [26], miRanda (version 4/06) [8], TargetScanS (version 4/05) [9], miRtarget (version 3/06) [28] and PicTar (version 5/05) [27]. The miRNAs are based on miRs registry version 8.1 (5/06).

Functional annotation and prioritization of p53-miRs

To annotate the p53-miRs and to identify/prioritize p53-miRs that might have potential implications in tumorigenesis we performed functional enrichment analysis and also intersected our predictions with cancer-associated miRNAs. Specifically, we used (a) literature reported differentially expressed miRNAs in cancer; (b) miRNAs reported as induced or repressed following p53 activation; and (c) miRNA targets that were functionally enriched with biological processes like apoptosis, cell cycle, etc.

Known tumor suppressor and oncogenic microRNAs

We manually compiled the list of microRNAs that were reported in literature as differentially expressed (up- or down-regulated) in various cancers in human (Additional File 8).

miRNAs induced or repressed following p53 activation

We intersected our predicted p53-miRs with the data from Tarasov et al [11] to identify potential p53 induced/repressed miRNAs. Following a genome-wide screen for microRNAs regulated by the transcription factor encoded by the p53 tumor suppressor gene, the authors reported that after p53-activation the abundance of thirty-four miRNAs was significantly increased, whereas sixteen miRNAs were suppressed [11].

p53-miRtarget genes – Functional enrichment

Functional over-representation analysis was performed to objectively identify biological processes potentially affected by p53-miR target genes. Specifically, the percentage of p53-miR target genes with a given gene ontology (GO) or pathway annotation was compared with the percentage of all miRNA-target genes genome-wide with the same annotation. A significant p value (p < 0.05) indicates that the observed percentage of p53-miR target genes with a given annotation could not likely occur by chance given the frequency of miRNA-target genes genome-wide with the same annotation. For functional enrichment analysis, we used DAVID [62] and ToppGene [63].

p53-miRand p53-interactant gene networks

p53-miR and p53-interactant gene networks were generated using aiSee [76] with force-directed layouts. The GDL (Graphics Description Language) files, generated by a locally-implemented JAVA GDL library, were used as input to aiSee to create SVG (Scalable Vector Graphics) network graphs. The p53-related interactions were downloaded from BioGRID [71].

Statistical methods

Fisher's exact test was used to calculate a p-value estimating the probability that a particular category of miRNAs is associated with a particular pattern or cluster of miRs more than would be expected by chance. Fisher's exact test was performed online at the MATFORSK [77]. For testing the functional enrichment of p53-miR target genes, DAVID [62] and ToppGene server [63] were used. For statistical analysis, the programming language and statistical environment R [78] were used.


  1. 1.

    Heinrichs S, Deppert W: Apoptosis or growth arrest: modulation of the cellular response to p53 by proliferative signals. Oncogene. 2003, 22 (4): 555-571. 10.1038/sj.onc.1206138.

  2. 2.

    Sionov RV, Haupt Y: The cellular response to p53: the decision between life and death. Oncogene. 1999, 18 (45): 6145-6157. 10.1038/sj.onc.1203130.

  3. 3.

    Dickins RA, Hemann MT, Zilfou JT, Simpson DR, Ibarra I, Hannon GJ, Lowe SW: Probing tumor phenotypes using stable and regulated synthetic microRNA precursors. Nat Genet. 2005, 37 (11): 1289-1295.

  4. 4.

    Spurgers KB, Gold DL, Coombes KR, Bohnenstiehl NL, Mullins B, Meyn RE, Logothetis CJ, McDonnell TJ: Identification of cell cycle regulatory genes as principal targets of p53-mediated transcriptional repression. J Biol Chem. 2006, 281 (35): 25134-25142. 10.1074/jbc.M513901200.

  5. 5.

    Ambros V: MicroRNA pathways in flies and worms: growth, death, fat, stress, and timing. Cell. 2003, 113 (6): 673-676. 10.1016/S0092-8674(03)00428-8.

  6. 6.

    Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297. 10.1016/S0092-8674(04)00045-5.

  7. 7.

    Griffiths-Jones S: miRBase: the microRNA sequence database. Methods Mol Biol. 2006, 342: 129-138.

  8. 8.

    John B, Enright AJ, Aravin A, Tuschl T, Sander C, Marks DS: Human MicroRNA targets. PLoS Biol. 2004, 2 (11): e363-10.1371/journal.pbio.0020363.

  9. 9.

    Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB: Prediction of mammalian microRNA targets. Cell. 2003, 115 (7): 787-798. 10.1016/S0092-8674(03)01018-3.

  10. 10.

    Harris SL, Levine AJ: The p53 pathway: positive and negative feedback loops. Oncogene. 2005, 24 (17): 2899-2908. 10.1038/sj.onc.1208615.

  11. 11.

    Tarasov V, Jung P, Verdoodt B, Lodygin D, Epanchintsev A, Menssen A, Meister G, Hermeking H: Differential Regulation of microRNAs by p53 Revealed by Massively Parallel Sequencing: miR-34a is a p53 Target That Induces Apoptosis and G(1)-arrest. Cell Cycle. 2007, 6 (13):

  12. 12.

    Chang TC, Wentzel EA, Kent OA, Ramachandran K, Mullendore M, Lee KH, Feldmann G, Yamakuchi M, Ferlito M, Lowenstein CJ, Arking DE, Beer MA, Maitra A, Mendell JT: Transactivation of miR-34a by p53 Broadly Influences Gene Expression and Promotes Apoptosis. Mol Cell. 2007, 26 (5): 745-752. 10.1016/j.molcel.2007.05.010.

  13. 13.

    Raver-Shapira N, Marciano E, Meiri E, Spector Y, Rosenfeld N, Moskovits N, Bentwich Z, Oren M: Transcriptional Activation of miR-34a Contributes to p53-Mediated Apoptosis. Mol Cell. 2007, 26 (5): 731-743. 10.1016/j.molcel.2007.05.017.

  14. 14.

    He L, He X, Lim LP, de Stanchina E, Xuan Z, Liang Y, Xue W, Zender L, Magnus J, Ridzon D, Jackson AL, Linsley PS, Chen C, Lowe SW, Cleary MA, Hannon GJ: A microRNA component of the p53 tumour suppressor network. Nature. 2007

  15. 15.

    Bommer GT, Gerin I, Feng Y, Kaczorowski AJ, Kuick R, Love RE, Zhai Y, Giordano TJ, Qin ZS, Moore BB, Macdougald OA, Cho KR, Fearon ER: p53-Mediated Activation of miRNA34 Candidate Tumor-Suppressor Genes. Curr Biol. 2007

  16. 16.

    He L, Thomson JM, Hemann MT, Hernando-Monge E, Mu D, Goodson S, Powers S, Cordon-Cardo C, Lowe SW, Hannon GJ, Hammond SM: A microRNA polycistron as a potential human oncogene. Nature. 2005, 435 (7043): 828-833. 10.1038/nature03552.

  17. 17.

    Hammond SM: MicroRNAs as tumor suppressors. Nat Genet. 2007, 39 (5): 582-583. 10.1038/ng0507-582.

  18. 18.

    Zhang B, Pan X, Cobb GP, Anderson TA: microRNAs as oncogenes and tumor suppressors. Dev Biol. 2007, 302 (1): 1-12. 10.1016/j.ydbio.2006.08.028.

  19. 19.

    Wu W, Sun M, Zou GM, Chen J: MicroRNA and cancer: Current status and prospective. Int J Cancer. 2007, 120 (5): 953-960. 10.1002/ijc.22454.

  20. 20.

    Voorhoeve PM, le Sage C, Schrier M, Gillis AJ, Stoop H, Nagel R, Liu YP, van Duijse J, Drost J, Griekspoor A, Zlotorynski E, Yabuta N, De Vita G, Nojima H, Looijenga LH, Agami R: A genetic screen implicates miRNA-372 and miRNA-373 as oncogenes in testicular germ cell tumors. Cell. 2006, 124 (6): 1169-1181. 10.1016/j.cell.2006.02.037.

  21. 21.

    Yang Y, Tantoso E, Chua GH, Yeo ZX, Ng FS, Wong ST, Chung CW, Li KB: In silico analysis of p53 using the P53 Knowledgebase: mutations, polymorphisms, microRNAs and pathways. In Silico Biol. 2007, 7 (1): 61-75.

  22. 22.

    Hoh J, Jin S, Parrado T, Edington J, Levine AJ, Ott J: The p53MH algorithm and its application in detecting p53-responsive genes. Proc Natl Acad Sci U S A. 2002, 99 (13): 8467-8472. 10.1073/pnas.132268899.

  23. 23.

    Saini HK, Griffiths-Jones S, Enright AJ: Genomic analysis of human microRNA transcripts. Proc Natl Acad Sci U S A. 2007, 104 (45): 17719-17724. 10.1073/pnas.0703890104.

  24. 24.

    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.2722704.

  25. 25.

    MAMI. []

  26. 26.

    Kiriakidou M, Nelson PT, Kouranov A, Fitziev P, Bouyioukos C, Mourelatos Z, Hatzigeorgiou A: A combined computational-experimental approach predicts human microRNA targets. Genes Dev. 2004, 18 (10): 1165-1178. 10.1101/gad.1184704.

  27. 27.

    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 (5): 495-500. 10.1038/ng1536.

  28. 28.

    Wang X, Wang X: Systematic identification of microRNA functions by combining target prediction and expression profiling. Nucleic Acids Res. 2006, 34 (5): 1646-1652. 10.1093/nar/gkl068.

  29. 29.

    Bentwich I: Prediction and validation of microRNAs and their targets. FEBS Lett. 2005, 579 (26): 5904-5910. 10.1016/j.febslet.2005.09.040.

  30. 30.

    Arias AM, Hayward P: Filtering transcriptional noise during development: concepts and mechanisms. Nat Rev Genet. 2006, 7 (1): 34-44. 10.1038/nrg1750.

  31. 31.

    Milo R, Shen-Orr S, Itzkovitz S, Kashtan N, Chklovskii D, Alon U: Network motifs: simple building blocks of complex networks. Science. 2002, 298 (5594): 824-827. 10.1126/science.298.5594.824.

  32. 32.

    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/msb4100089.

  33. 33.

    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.080.

  34. 34.

    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.0030131.

  35. 35.

    Tsang J, Zhu J, van Oudenaarden A: MicroRNA-mediated feedback and feedforward loops are recurrent network motifs in mammals. Mol Cell. 2007, 26 (5): 753-767. 10.1016/j.molcel.2007.05.018.

  36. 36.

    Lim YP, Lim TT, Chan YL, Song AC, Yeo BH, Vojtesek B, Coomber D, Rajagopal G, Lane D: The p53 knowledgebase: an integrated information resource for p53 research. Oncogene. 2007, 26 (11): 1517-1521. 10.1038/sj.onc.1209952.

  37. 37.

    Hua Z, Lv Q, Ye W, Wong CK, Cai G, Gu D, Ji Y, Zhao C, Wang J, Yang BB, Zhang Y: MiRNA-Directed Regulation of VEGF and Other Angiogenic Factors under Hypoxia. PLoS ONE. 2006, 1: e116-10.1371/journal.pone.0000116.

  38. 38.

    Meng F, Henson R, Wehbe-Janek H, Ghoshal K, Jacob ST, Patel T: MicroRNA-21 regulates expression of the PTEN tumor suppressor gene in human hepatocellular cancer. Gastroenterology. 2007, 133 (2): 647-658. 10.1053/j.gastro.2007.05.022.

  39. 39.

    Sethupathy P, Corda B, Hatzigeorgiou AG: TarBase: A comprehensive database of experimentally supported animal microRNA targets. Rna. 2006, 12 (2): 192-197. 10.1261/rna.2239606.

  40. 40.

    Im HJ, Pittelkow MR, Kumar R: Divergent regulation of the growth-promoting gene IEX-1 by the p53 tumor suppressor and Sp1. J Biol Chem. 2002, 277 (17): 14612-14621.

  41. 41.

    Farh KK, Grimson A, Jan C, Lewis BP, Johnston WK, Lim LP, Burge CB, Bartel DP: The widespread impact of mammalian MicroRNAs on mRNA repression and evolution. Science. 2005, 310 (5755): 1817-1821. 10.1126/science.1121158.

  42. 42.

    Sood P, Krek A, Zavolan M, Macino G, Rajewsky N: Cell-type-specific signatures of microRNAs on target mRNA expression. Proc Natl Acad Sci U S A. 2006, 103 (8): 2746-2751. 10.1073/pnas.0511045103.

  43. 43.

    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/nature03677.

  44. 44.

    Okamoto K, Beach D: Cyclin G is a transcriptional target of the p53 tumor suppressor protein. Embo J. 1994, 13 (20): 4816-4822.

  45. 45.

    Gramantieri L, Ferracin M, Fornari F, Veronese A, Sabbioni S, Liu CG, Calin GA, Giovannini C, Ferrazzi E, Grazi GL, Croce CM, Bolondi L, Negrini M: Cyclin G1 is a target of miR-122a, a microRNA frequently down-regulated in human hepatocellular carcinoma. Cancer Res. 2007, 67 (13): 6092-6099. 10.1158/0008-5472.CAN-06-4607.

  46. 46.

    Ohtsuka T, Ryu H, Minamishima YA, Ryo A, Lee SW: Modulation of p53 and p73 levels by cyclin G: implication of a negative feedback regulation. Oncogene. 2003, 22 (11): 1678-1687. 10.1038/sj.onc.1206306.

  47. 47.

    Lee TI, Rinaldi NJ, Robert F, Odom DT, Bar-Joseph Z, Gerber GK, Hannett NM, Harbison CT, Thompson CM, Simon I, Zeitlinger J, Jennings EG, Murray HL, Gordon DB, Ren B, Wyrick JJ, Tagne JB, Volkert TL, Fraenkel E, Gifford DK, Young RA: Transcriptional regulatory networks in Saccharomyces cerevisiae. Science. 2002, 298 (5594): 799-804. 10.1126/science.1075090.

  48. 48.

    Goldbeter A, Koshland DE: Ultrasensitivity in biochemical systems controlled by covalent modification. Interplay between zero-order and multistep effects. J Biol Chem. 1984, 259 (23): 14441-14447.

  49. 49.

    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/nrg1328.

  50. 50.

    Saito Y, Liang G, Egger G, Friedman JM, Chuang JC, Coetzee GA, Jones PA: Specific activation of microRNA-127 with downregulation of the proto-oncogene BCL6 by chromatin-modifying drugs in human cancer cells. Cancer Cell. 2006, 9 (6): 435-443. 10.1016/j.ccr.2006.04.020.

  51. 51.

    Phan RT, Dalla-Favera R: The BCL6 proto-oncogene suppresses p53 expression in germinal-centre B cells. Nature. 2004, 432 (7017): 635-639. 10.1038/nature03147.

  52. 52.

    Volinia S, Calin GA, Liu CG, Ambs S, Cimmino A, Petrocca F, Visone R, Iorio M, Roldo C, Ferracin M, Prueitt RL, Yanaihara N, Lanza G, Scarpa A, Vecchione A, Negrini M, Harris CC, Croce CM: A microRNA expression signature of human solid tumors defines cancer gene targets. Proc Natl Acad Sci U S A. 2006, 103 (7): 2257-2261. 10.1073/pnas.0510565103.

  53. 53.

    Shiio Y, Yamamoto T, Yamaguchi N: Negative regulation of Rb expression by the p53 gene product. Proc Natl Acad Sci U S A. 1992, 89 (12): 5206-5210. 10.1073/pnas.89.12.5206.

  54. 54.

    Ho J, Benchimol S: Transcriptional repression mediated by the p53 tumour suppressor. Cell Death Differ. 2003, 10 (4): 404-408. 10.1038/sj.cdd.4401191.

  55. 55.

    Menendez D, Inga A, Jordan JJ, Resnick MA: Changing the p53 master regulatory network: ELEMENTary, my dear Mr Watson. Oncogene. 2007, 26 (15): 2191-2201. 10.1038/sj.onc.1210277.

  56. 56.

    Ho JS, Ma W, Mao DY, Benchimol S: p53-Dependent transcriptional repression of c-myc is required for G1 cell cycle arrest. Mol Cell Biol. 2005, 25 (17): 7423-7431. 10.1128/MCB.25.17.7423-7431.2005.

  57. 57.

    Zhai W, Comai L: Repression of RNA polymerase I transcription by the tumor suppressor p53. Mol Cell Biol. 2000, 20 (16): 5930-5938. 10.1128/MCB.20.16.5930-5938.2000.

  58. 58.

    Xu H, El-Gewely MR: Differentially expressed downstream genes in cells with normal or mutated p53. Oncol Res. 2003, 13 (6-10): 429-436.

  59. 59.

    Kang JH, Kim SJ, Noh DY, Park IA, Choe KJ, Yoo OJ, Kang HS: Methylation in the p53 promoter is a supplementary route to breast carcinogenesis: correlation between CpG methylation in the p53 promoter and the mutation of the p53 gene in the progression from ductal carcinoma in situ to invasive ductal carcinoma. Lab Invest. 2001, 81 (4): 573-579.

  60. 60.

    Raman V, Martensen SA, Reisman D, Evron E, Odenwald WF, Jaffee E, Marks J, Sukumar S: Compromised HOXA5 function can limit p53 expression in human breast tumours. Nature. 2000, 405 (6789): 974-978. 10.1038/35016125.

  61. 61.

    Bandres E, Cubedo E, Agirre X, Malumbres R, Zarate R, Ramirez N, Abajo A, Navarro A, Moreno I, Monzo M, Garcia-Foncillas J: Identification by Real-time PCR of 13 mature microRNAs differentially expressed in colorectal cancer and non-tumoral tissues. Mol Cancer. 2006, 5: 29-10.1186/1476-4598-5-29.

  62. 62.

    Sherman BT, Huang DW, Tan Q, Guo Y, Bour S, Liu D, Stephens R, Baseler MW, Lane HC, Lempicki RA: DAVID Knowledgebase: a gene-centered database integrating heterogeneous gene annotation resources to facilitate high-throughput gene functional analysis. BMC Bioinformatics. 2007, 8 (1): 426-10.1186/1471-2105-8-426.

  63. 63.

    Chen J, Xu H, Aronow BJ, Jegga AG: Improved human disease candidate gene prioritization using mouse phenotype. BMC Bioinformatics. 2007, 8 (1): 392-10.1186/1471-2105-8-392.

  64. 64.

    Bonaccorsi L, Nosi D, Muratori M, Formigli L, Forti G, Baldi E: Altered endocytosis of epidermal growth factor receptor in androgen receptor positive prostate cancer cell lines. J Mol Endocrinol. 2007, 38 (1-2): 51-66. 10.1677/jme.1.02155.

  65. 65.

    Floyd S, De Camilli P: Endocytosis proteins and cancer: a potential link?. Trends Cell Biol. 1998, 8 (8): 299-301. 10.1016/S0962-8924(98)01316-6.

  66. 66.

    Polo S, Pece S, Di Fiore PP: Endocytosis and cancer. Curr Opin Cell Biol. 2004, 16 (2): 156-161. 10.1016/

  67. 67.

    Kawasaki H, Taira K: Functional analysis of microRNAs during the retinoic acid-induced neuronal differentiation of human NT2 cells. Nucleic Acids Res Suppl. 2003, 243-244.

  68. 68.

    van Rooij E, Sutherland LB, Liu N, Williams AH, McAnally J, Gerard RD, Richardson JA, Olson EN: A signature pattern of stress-responsive microRNAs that can evoke cardiac hypertrophy and heart failure. Proc Natl Acad Sci U S A. 2006, 103 (48): 18255-18260. 10.1073/pnas.0608791103.

  69. 69.

    Luo J, Nikolaev AY, Imai S, Chen D, Su F, Shiloh A, Guarente L, Gu W: Negative control of p53 by Sir2alpha promotes cell survival under stress. Cell. 2001, 107 (2): 137-148. 10.1016/S0092-8674(01)00524-4.

  70. 70.

    Liang H, Li WH: MicroRNA regulation of human protein protein interaction network. Rna. 2007, 13 (9): 1402-1408. 10.1261/rna.634607.

  71. 71.

    Stark C, Breitkreutz BJ, Reguly T, Boucher L, Breitkreutz A, Tyers M: BioGRID: a general repository for interaction datasets. Nucleic Acids Res. 2006, 34 (Database issue): D535-9. 10.1093/nar/gkj109.

  72. 72.

    Human miRNA registry. []

  73. 73.

    Kuhn RM, Karolchik D, Zweig AS, Trumbower H, Thomas DJ, Thakkapallayil A, Sugnet CW, Stanke M, Smith KE, Siepel A, Rosenbloom KR, Rhead B, Raney BJ, Pohl A, Pedersen JS, Hsu F, Hinrichs AS, Harte RA, Diekhans M, Clawson H, Bejerano G, Barber GP, Baertsch R, Haussler D, Kent WJ: The UCSC genome browser database: update 2007. Nucleic Acids Res. 2007, 35 (Database issue): D668-73. 10.1093/nar/gkl928.

  74. 74.

    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.207.

  75. 75.

    el-Deiry WS, Kern SE, Pietenpol JA, Kinzler KW, Vogelstein B: Definition of a consensus binding site for p53. Nat Genet. 1992, 1 (1): 45-49. 10.1038/ng0492-45.

  76. 76.

    aiSee. []

  77. 77.

    MATFORSK. []

  78. 78.

    R programming language and statistical environment. []

Download references


This study was supported by the State of Ohio Computational Medicine Center (ODD TECH 04–042) and NCI UO1 CA84291-07 (Mouse Models of Human Cancer Consortium). The authors also acknowledge the support of Drs. Bruce Aronow (Division of Biomedical Informatics, Cincinnati Children's Hospital Medical Center) and Raj Bhatnagar (Department of Computer Science, University of Cincinnati).

Author information

Correspondence to Anil G Jegga.

Additional information

Authors' contributions

AS, VK and AJ conceived the study design which was coordinated by AJ. AS and VK designed and implemented the p53-miR discovery and prioritization approaches and along with AJ participated in the analysis and interpretation of results. JC carried out the p53 interactome data analysis. AJ drafted the manuscript. All the authors have read and approved the final manuscript

Amit U Sinha, Vivek Kaimal contributed equally to this work.

Electronic supplementary material

Additional File 1: Systematic break-down of the putative p53-regulated miRNAs (p53-miRs). Schematic classification of putative p53 sites and p53-miRs based on the miRNA types (intragenic or intergenic). (PDF 185 KB)

Additional File 2: Genomic coordinates of putative p53 sites within 10 kb flanking regions of human miRNAs. Chromosome cooridnates (based on Human Mar. 2006 – hg18 – assembly) of putative p53 sites in the flanking 10 kb regions of human microRNAs. (PDF 369 KB)

Additional File 3: Putative miRNAs regulating known 23 upstream regulators and 48 downstream TFs of p53. For the known 23 upstream regulators and 48 downstream TFs of p53 (based on p53 Knowledgebase;, putative miRs regulating them were extracted using MAMI server and the database (PDF 260 KB)

Additional File 4: Functional enrichment (based on GO and pathways) in target genes of 143 p53-miRs. Results of functional over-representation analysis, performed to objectively identify biological processes potentially affected by p53-miR target genes. (PDF 46 KB)

Additional File 5: Details of Fisher's exact test used to test whether p53-miRs tend to target p53-interactants. Fisher's exact test was used to calculate the probability of p53-miRs regulating the p53 interactants (downloaded from the BioGRID database). (PDF 37 KB)

Additional File 6: Systematic break-down of the 474 human miRNAs based on their genomic location. Classification of the known 474 human miRNAs based on their genomic location (intergenic or intragenic/intronic). (PDF 130 KB)

Additional File 7: Details of the 474 human miRNAs analyzed. Additional details of genomic locations of the 474 human miRNAs analyzed (whether part of a cluster, i.e., 2 or more miRNAs occurring together on the genome; or located in the intergenic or intragenic/intronic regions of known genes). (PDF 154 KB)

Additional File 8: Manually compiled list of miRNAs that are reported in the literature as either up or down regulated in various human cancers or cancer cell lines. This file has a list of miRNAs that are reported in the literature as either up or down regulated in various human cancers or cancer cell lines. These "cancer-associated" miRNAs were used to prioritize putative p53-miR targets. (PDF 136 KB)

Authors’ original submitted files for images

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Authors’ original file for figure 4

Authors’ original file for figure 5

Authors’ original file for figure 6

Rights and permissions

Reprints and Permissions

About this article


  • Downstream Target Gene
  • Functional Enrichment Analysis
  • Human miRNAs
  • Target Prediction Algorithm
  • Cell Proliferation Gene