|
Advertisement | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 281, Issue 35, 25134-25142, September 1, 2006
Identification of Cell Cycle Regulatory Genes as Principal Targets of p53-mediated Transcriptional Repression*
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| ABSTRACT |
|---|
|
|
|---|
0.05). Validation of the array data, using reverse transcription-PCR of 20 randomly selected genes, yielded a confirmation rate of >95.5% for the complete data set. Functional over-representation analysis revealed that cell cycle regulatory genes exhibited a highly significant enrichment (p
5 x 1028) within the transrepressed targets. 41% of the repressed targets are cell cycle regulators. A subset of these genes exhibited repression following DNA damage, preceding cell cycle arrest, in LNCaP cells. The use of a p53 small interfering RNA strategy in LNCaP cells and the use of p53-null cell lines demonstrated that this repression is p53-dependent. These findings identify a set of genes not known previously to be down-regulated by p53 and indicate that p53-induced cell cycle arrest is a function of not only the transactivation of cell cycle inhibitors (e.g. p21) but also the repression of targets that regulate proliferation at several distinct phases of the cell cycle. | INTRODUCTION |
|---|
|
|
|---|
Although transcription-independent functions exist, p53 mediates its effects largely by regulating the expression of downstream target genes (6, 7). Most studies investigating p53 function have focused attention on the genes transactivated by p53. However, it is recognized that repression of target genes may be important for p53-induced cell death and cell cycle arrest. For example, trichostatin A, a histone deacetylase inhibitor, inhibits the repression of p53 targets and the ability of p53 to induce apoptosis (8). Additionally, cell cycle regulators such as cdc2 and cyclin B can be repressed by p53 (9, 10). Beyond apoptosis and cell cycle arrest, p53 can affect other cellular processes including DNA repair, senescence, and differentiation (1). Despite documented examples (11, 12), relatively little is known regarding the genes transcriptionally repressed by p53 and what roles they play in mediating p53 function. These repressed genes are of further interest in that they may potentially be overexpressed as a consequence of somatic p53 mutation during cancer progression.
Using high density oligonucleotide microarrays we identified a total of 111 genes that were significantly repressed following adenoviral p53 gene transfer (Ad-p53) in PC3 prostate cancer cells. Functional over-representation analysis was used to objectively identify which cellular processes were influenced by p53-mediated repression. Notably, nearly half of the repressed genes (45 of 111 genes) are involved in cell cycle regulation. This represents a highly significant enrichment of this functional category within the repressed targets compared with the microarray as a whole. Importantly, several cell cycle genes were repressed after genotoxic stress only in cells harboring wild type p53 alleles and before cell cycle arrest is evident. The p53-dependent nature of this repression was further demonstrated using a p53 knockdown strategy. These findings identify cell cycle regulatory genes previously unrecognized as p53-responsive and indicate that p53-induced cell cycle arrest is a function of not only the transactivation of cell cycle inhibitors (e.g. p21), but also the repression of targets that regulate proliferation at several distinct phases of the cell cycle.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
Adenovirus Infection and Microarray HybridizationAll of the adenoviral vector constructs were maintained by the Vector Core Facility of the University of Texas M. D. Anderson Cancer Center (Houston, TX). Construction of the recombinant, replication-deficient, type 5 adenovirus expressing human p53 (Ad-p53) and the experimental parameters and procedures used in this study have been described (13, 15). Briefly, PC3 cells were infected with control virus expressing
-galactosidase (Ad-LacZ) or with Ad-p53 at 1,000 virus particles/cell. RNA was harvested 19 h after infection. Target preparation and Affymetrix (Santa Clara, CA) U133A GeneChip hybridizations have been described (16) and were carried out by the M. D. Anderson Cancer Center Microarray Core Facility (Houston, TX).
Array Data AnalysisThe arrays were preprocessed with dCHIP version 1.3 software (www.dchip.org). Gene expression was obtained after normalization with the Perfect Match only model (17). dCHIP array level quality metrics were examined to identify outliers. The arrays were also assessed for brightness, alignment calibration, and spatial variation following the methods of Gold et al. (16).
Initially, a detection filter was applied to retain only those probe sets called present by dCHIP in all six replicates of at least one biological factor combination. Of the 22,283 probe sets on the U133A microarray, 12,308 met this detection criterion. The log2-transformed expressions of the probe sets were modeled linearly,
![]() |
for the presence of a p53 effect,
for a BCL2 effect, and 
for p53-BCL2 interaction. Here µ accounts for expected log2 expression. Unexplained random variation
was assumed normally distributed with expectation 0 and variance
. We explored possible hybridization day and RNA extraction contributions to variation in gene expression across arrays. These were not considerable, agreeing with hierarchical cluster analysis and thus were omitted from Equations 1 and 2. Genes with significantly altered expression were identified using an F-statistic pass at the 0.05 Bonferroni-corrected significance level (p < 0.05/12,038) to test the hypothesis that the control factors in Equation 1 were simultaneously zero,
![]() |
2g|
1 on the log2 scale. Individual t-statistics for the hypotheses Ho: 1 <
2g < 1 Ha:
2g
1 or
2g
1 were tested at a 0.05 Bonferroni adjusted significance level. While accounting for p53-BCL2 interactions, we found 391 probe sets with significant change coincidental with p53 expression beyond a fold change of 2. These 391 probe sets correspond to 335 unique genes (224 up-regulated and 111 repressed following p53 expression). Hierarchical Cluster AnalysisWe compared relative sources of variation between arrays by using hierarchical cluster analysis. We required probe sets to have expression levels above the median signal on at least six arrays to be included in the analysis (n = 10,625 probe sets). Hierarchical cluster analysis was applied with complete linkage and Euclidean distance on log2 expression. Furthermore, the respective profiles for probe set targets of p53 were examined on the log2 scale.
Estimation of Array Data False Positive RateA random number generator following a discrete uniform distribution was used to select 20 differentially expressed genes for verification. To estimate the fraction
of "true positives" on our list of differentially expressed genes, we modeled the number of successful confirmations as a binomial random variable X =
Binom(N,
). A Bayesian statistical method was used to estimate the probability distribution of
. We looked at various beta prior distributions,
=
Beta(2p, 2(1p)), corresponding to different choices of the parameter p, the prior probability of confirmation. The posterior distribution of
, conditional on confirming k of N genes, is another beta distribution, Beta(2p + k, 2(1p) + Nk). All 20 of the genes tested were successfully confirmed, giving posterior distributions of the form Beta(2p + 20, 2(1p) + 20). We computed the expected value and 95% confidence interval for the confirmation rate
for a wide range of prior distributions.
RT-PCRFor RT-PCRs, 3 µg of total RNA was combined with 500 ng of oligo(dT) primer (Invitrogen) in a final volume of 12 µl, heated at 70 °C for 10 min, and then placed on ice. Reverse transcription reactions (20 µl) were then assembled with final concentrations of 1x reaction buffer, 10 mM dithiothreitol, 0.5 mM dNTPs, and 100 units of Superscript II reverse transcriptase (Invitrogen). Two percent of the RT reaction was used as input for each PCR. PCRs (50 µl) were assembled with final concentrations of 1x buffer, 1.5 mM MgCl2, 0.2 mM dNTPs, 0.2 µM of each primer, and 2.5 units of Taq polymerase (Invitrogen). Primers were obtained from Sigma-Genosys (The Woodlands, TX). The sequence of each primer is available upon request. PCR cycle parameters were as follows: initial step of 95 °C for 5 min and then 2035 cycles of 95 °C for 30 s, annealing temperature of 5961 °C for 30 s, 72 °C for 45 s, and a final step of 72 °C for 5 min. The cycle number that allows for end point analysis during linear amplification was determined empirically.
Functional Over-representation AnalysisThe Expression Analysis Systematic Explorer software package was used to identify biological themes that are significantly enriched in our list of p53-repressed genes (18). Gene categories with an Expression Analysis Systematic Explorer score of
0.05 were considered significantly enriched. Bonferroni multiplicity correction was not employed.
siRNA TransfectionLNCaP cells were transfected with p53-specific or p21waf1-specific siRNA SmartPool duplexes or nontargeting siRNA number 1 (Dharmacon, Lafayette, CO) at a final concentration of 100 nM. Transfections were performed using Lipofectamine 2000 (Invitrogen) following the manufacturer's protocol. A total of 30 and 60 µl of Lipofectamine 2000 were used for each 60-mm and 10-cm plate, respectively. Transfection medium was replaced with fresh complete medium after 7 h.
Western BlotWestern blots were generated and probed as described (13). The primary antibodies used were anti-p53 at 1:1,500 (Bp53-12; Santa Cruz Biotechnology, Santa Cruz, CA), anti-p21 at 1:500 (Ab-1; Calbiochem, San Diego, CA), and antiactin at 1:1,000 (AC-40; Sigma).
BrdUrd IncorporationBrdUrd incorporation was measured in LNCaP cells using the BrdUrd Flow Kit (BD Pharmingen, San Diego, CA) following the manufacturer's protocol with some modifications. BrdUrd was added directly to culture medium at a final concentration of 10 µM. Forty minutes after the addition of BrdUrd, the cells were harvested with trypsin, washed once in 1x phosphate-buffered saline, and fixed in 100 µl of Cytofix/Cytoperm for 30 min at room temperature. Incubation with Cytoperm Plus and a second fixation with Cytofix/Cytoperm were omitted. The cells were read on a FACSCalibur flow cytometer (BD Biosciences, San Jose, CA).
| RESULTS |
|---|
|
|
|---|
Microarray experiments and subsequent statistical analysis were designed to minimize the effect of inherent experimental variation. Six replicate hybridizations (with probes from two independent virus infections and RNA extractions) were performed for each cell line/treatment combination. Hierarchical cluster analysis was used to assess sources of variation across the arrays. The most dominant feature affecting the clustering was p53 expression (not shown). RNA extraction or hybridization day were insignificant sources of variation across arrays.
To identify p53-responsive genes, control virus treatments were compared with Ad-p53 treatments. A total of 111 genes were repressed (fold change
2, p
0.05) in both PC3-BN and PC3-BCL2 cell lines following p53 expression. The repression of these 111 genes was not affected by BCL2 overexpression. No genes were detected that were repressed in one cell line and not the other. The complete list of p53-repressed genes is available as supplementary data (supplemental Table S1).
Objective Assessment of Array Data ReliabilityGiven the nature of microarray experiments and the statistical methods employed, the presence of false positives is inherent in any evaluation of differentially expressed genes. With this in mind, we objectively assessed the false positive rate within our final data set. In addition to p53-repressed genes, we also identified genes transcriptionally activated in response to p53. Twenty genes were chosen at random (10 activated and 10 repressed) and tested for a response to p53 by RT-PCR. The expression patterns of all 20 genes were validated (Fig. 1 and not shown). Given this data, we used a Bayesian statistical method to estimate the false positive rate for the complete data set. Using a neutral, uniform, prior probability for confirmation, the estimate of the confirmation rate is 95.5% with a 95% confidence interval of 86.7100%. Thus,
4.5% of the genes identified as differentially expressed are expected to be false positives.
Even with a very low expected frequency of false positives, our data set does contain at least some false negative results. For example, survivin is strongly down-regulated in PC3 cells following Ad-p53 treatment (13). This gene is represented on the U133A microarray but was not identified as a p53-repressed target in our analysis. Despite the presence of some false negatives, our data set shows considerable overlap with published expression profiling studies looking at p53-responsive gene expression and also identifies novel p53 target genes (23, 24).
Functional Over-representation AnalysisFunctional over-representation analysis was performed to objectively identify biological processes potentially affected by p53-mediated transcriptional repression. Specifically, the percentage of p53-repressed genes with a given gene ontology (GO) annotation was compared with the percentage of genes on the Affymetrix U133A GeneChip with the same annotation. A significant p value (p
0.05) indicates that the observed percentage of p53-repressed genes with a given annotation could not likely occur by chance given the frequency of genes on the microarray with the same annotation. Fig. 2 displays nine GO categories that were significantly enriched (p
0.05) in our list of p53-repressed genes. The GO annotation "cell cycle" shows the most significant enrichment. Specifically, 45 of the 111 p53-repressed genes identified in our study are involved in cell cycle regulation (Table 1). Most of these genes also have more specific GO annotations suggesting that several phases of the cell cycle are affected by p53-mediated transcriptional repression. x indicates that the gene product has a GO annotation for that specific cell cycle phase (Table 1). Genes that harbor consensus p53 DNA-binding sites, based on the data of Hoh et al. (25), are also indicated. Other GO categories of interest significantly enriched within the repressed targets include DNA repair, nucleotide metabolism, and DNA packaging (supplemental Table S2).
|
|
|
p53 dependence was assessed in a more direct manner using siRNA to prevent accumulation of p53 following DNA damage. LNCaP cells were transfected with negative control or p53-specific siRNA or not transfected. 48 h after siRNA transfection, the cells were treated with etoposide for 4 h and then analyzed by Western blot. Pretreatment with p53 siRNA prevented the accumulation of p53 protein following etoposide treatment (Fig. 4B). Next, the genes showing a DNA damage response in LNCaP cells were examined by RT-PCR following p53 siRNA transfection and etoposide (24 h) treatment. Each gene showed a reduction in RNA level in nontargeting siRNA transfected, etoposide-treated cells, compared with control cells (Fig. 4C). This reduction did not occur or was dramatically attenuated in etoposide-treated cells that had been transfected with p53 siRNA. These findings are consistent with the interpretation that down-regulation of these genes in LNCaP cells is specifically dependent on p53 protein accumulation.
Repression of Cell Cycle Genes Requires New Protein Synthesis but Not p21 TransactivationIn some cases the transactivation or repression of p53 target genes is known to require new protein synthesis following p53 expression (28). In these cases, changes in gene expression are thought to occur because of effects of primary p53 targets. In particular, some groups have suggested that p53-mediated repression of cell cycle genes requires the up-regulation of p21waf1 protein (24). We determined whether select cell cycle regulatory genes could be repressed in the presence of the protein synthesis inhibitor CHX. LNCaP cells were treated with etoposide in the presence or absence of CHX. Without CHX, as expected, p53 and p21 protein levels increased 4 h after etoposide treatment (Fig. 5A). The addition of CHX prevented the accumulation of p21 protein but had little effect on p53. CHX did not prevent the increase in p21 mRNA following etoposide treatment, demonstrating a specific inhibition of protein synthesis and not p53-regulated transcription (Fig. 5B). Significantly, CHX prevented the transcriptional repression of each p53-repressed, cell cycle regulatory gene examined (Fig. 5B). These findings suggest that new protein synthesis is necessary for repression and that repression likely occurs through the action of other p53-regulated genes.
We next specifically investigated the requirement of p21 protein up-regulation in the repression of cell cycle genes following DNA damage. LNCaP cells were transfected with negative control siRNA or siRNA specific for p21. The cells were then treated with Me2SO or etoposide after 48 h. Western blot and RT-PCR analyses were carried out 8 h after the addition of etoposide. In the presence of negative control siRNA, p21 was efficiently up-regulated following DNA damage (Fig. 6A). This accumulation of p21 protein was inhibited by p21 siRNA transfection. However, p21 siRNA had little or no effect on the transcriptional repression of select cell cycle genes after etoposide treatment (Fig. 6B). These data suggest that the p53-dependent repression of cell cycle genes following DNA damage does not require the accumulation of high levels of p21 protein.
| DISCUSSION |
|---|
|
|
|---|
|
|
(9, 23, 3234). However, to our knowledge, the majority of the repressed, cell cycle regulatory genes identified in our study have not been previously recognized as p53-responsive. Other processes possibly affected by p53-mediated transcriptional repression include DNA repair, nucleotide metabolism, and DNA packaging.
|
|
|
p53 is thought to repress the expression of some target genes through direct binding to consensus DNA-binding elements (11). p53 can also bind to novel DNA sequences that do not conform to the reported consensus sequence (37, 38). Based on the data of Hoh et al. (25), most of the genes we identified do not exhibit consensus p53 DNA-binding sites, yet many of these genes still undergo p53-dependent repression in response to DNA damage. Even without consensus DNA-binding sites, p53 may act directly, or indirectly via a DNA-binding protein complex, at the promoters of these repressed genes. For example, our analysis, consistent with previous reports, identifies c-Myc, and PRC1 as p53-repressed targets (32, 35). Using a temperature-sensitive p53 expression system, p53 can be detected at the promoters of both of these genes. At the region of p53 enrichment, p53 DNA-binding sequences are not present. The same is true for the p53-repressed genes MAD1L1 and MAP4 (8, 39). The absence of consensus p53 DNA-binding sequences suggests that p53 binds to a nonconsensus sequence or is acting in a complex with other DNA-binding proteins. Experiments are ongoing to determine whether p53 is present at the promoters of the repressed cell cycle genes identified in our study and what DNA sequences mediate this recruitment.
Some progress has been made in determining how p53 mediates repression once it has been recruited to the promoter of an affected gene. p53 can bind directly to the corepressor sin3a, which can recruit histone deacetylase activity to targeted promoters. By immunoprecipitation, p53 was shown to interact with sin3a and HDAC1 in MCF7 cells (8). Similarly, upon p53 expression, p53, sin3a, and HDAC1 can all be detected at the MAP4 promoter. A decrease in histone acetylation, presumably aiding in the repression of transcription, is also found on the promoter at this p53-enriched region. p53 may also influence promoter methylation to influence transcription. A direct interaction between p53 and the DNA methyltransferase DNMT1 has been detected (40). Experiments with DNMT1 null cells suggest that this protein is required for p53-mediated repression of survivin. Further, in response to DNA damage, both p53 and DNMT1 can be detected at the survivin promoter (40).
Repression of some p53-responsive genes requires the transactivation of p21, a cyclin-dependent kinase inhibitor (24). p21 expression can aid the formation of RB·E2F complexes that act to repress E2F target genes. RB family and E2F family proteins are required for the repression of cell cycle regulatory genes and cell cycle arrest following DNA damage or forced p53 expression in some cases (4143). Importantly, a large number of the p53-repressed, cell cycle regulatory genes identified in this study are known E2F targets (44). Experiments with CHX suggest that new protein synthesis is required for the repression observed in LNCaP cells following DNA damage. These data are consistent with a requirement for p21 transactivation. However, when p21 protein accumulation after DNA damage is attenuated with siRNA, we still see the repression of cell cycle genes. This is consistent with reports that p53-dependent repression of c-Myc and CDC25C can occur in the absence of p21 (35, 36). As demonstrated by Lohr et al. (24), several genes identified in our study exhibit p21-dependent repression in doxorubicin-treated HCT116 cells. These particular genes have not yet been tested for a requirement of p21 in our experimental system (etoposide-treated LNCaP cells). Further studies are needed to determine what other proteins are necessary for the repression of particular cell cycle genes in response to p53 expression.
Another possible mechanism of repression involves the NF-Y transcription factor. The promoters of many genes that control the G2/M transition contain multiple CCAAT boxes and are regulated in a cell cycle-dependent manner by NF-Y (45). Several groups have suggested that the negative regulation of G2/M promoters following DNA damage or p53 expression depends on NF-Y (46, 47). Also dependent on NF-Y, p53 can be found at the promoters of cyclin B2, CDC2, and CDC25C, before and after DNA damage, at regions that harbor CCAAT boxes but no consensus p53 DNA-binding sites (48). Several genes identified in our study contain NF-Y-binding sites.
G1 and G2 cell cycle arrest checkpoints in response to DNA damage have been described in detail (49). Importantly, several mechanisms for this arrest have been elucidated, some dependent on p53-mediated transcriptional changes and some p53-independent. Data presented here, and data from other groups, support the idea that a cell cycle target gene can be impacted by transcription-dependent and -independent mechanisms. CDC25A protein, the activity of which is required for S phase entry, is rapidly degraded in response to DNA damage, resulting in G1 arrest (50). To our knowledge, we provide the first evidence that CDC25A is a target for p53-dependent repression at the transcript level. As a critical regulator of cell cycle progression, it is not unreasonable for CDC25A to be the target of two arms of the DNA damage response, rapid degradation via the proteasome, and transcriptional repression by p53.
Our data further suggest that the initiation of DNA synthesis is a prominent target of p53 during the course of p53-induced cell cycle arrest. We identified several MCM (minichromosome maintenance-deficient) proteins, as well as CDC6, as targets for p53-mediated repression. These proteins participate in replication origin licensing and ensure that origins fire only once during a complete cell cycle. CDC6 binds to the origin recognition complex and aids the loading of the MCM27 hexamer (51). Several MCM proteins have been identified as p53-responsive in other expression profiling efforts (23).
Although there are reports that p53 can repress the transcription of a small number of cell cycle regulatory genes, this function of p53 is not widely appreciated. Data presented here substantially expand the known repertoire of cell cycle regulatory genes repressed in response to p53 activation. Importantly, p53-mediated repression of target genes likely impacts cell cycle progression at several distinct points. It can be anticipated that further characterization of new p53 transcriptionally regulated target genes will improve our understanding of the DNA damage response, multistep carcinogenesis, and response to therapy.
| FOOTNOTES |
|---|
The on-line version of this article (available at http://www.jbc.org) contains supplemental tables and a supplemental figure. ![]()
1 Supported in part by the American Legion Auxiliary Fellowship in Cancer Research. ![]()
2 To whom correspondence should be addressed: Dept. of Molecular Pathology, Box 089, University of Texas M. D. Anderson Cancer Center, 1515 Holcombe Blvd., Houston, TX 77030. Tel.: 713-834-6028; Fax: 713-745-6696; E-mail: tmcdonne{at}mdanderson.org.
3 The abbreviations used are: CHX, cycloheximide; RT, reverse transcription; siRNA, small interfering RNA; BrdUrd, bromodeoxyuridine; GO, gene ontology. ![]()
| REFERENCES |
|---|
|
|
|---|