Permalink
Cannot retrieve contributors at this time
Fetching contributors…
Cannot retrieve contributors at this time
| <?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.1d1 20130915//EN" "JATS-archivearticle1.dtd"><article article-type="research-article" dtd-version="1.1d1" xmlns:xlink="http://www.w3.org/1999/xlink"><front><journal-meta><journal-id journal-id-type="nlm-ta">elife</journal-id><journal-id journal-id-type="hwp">eLife</journal-id><journal-id journal-id-type="publisher-id">eLife</journal-id><journal-title-group><journal-title>eLife</journal-title></journal-title-group><issn publication-format="electronic">2050-084X</issn><publisher><publisher-name>eLife Sciences Publications, Ltd</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="publisher-id">02200</article-id><article-id pub-id-type="doi">10.7554/eLife.02200</article-id><article-categories><subj-group subj-group-type="display-channel"><subject>Research article</subject></subj-group><subj-group subj-group-type="heading"><subject>Genes and chromosomes</subject></subj-group><subj-group subj-group-type="heading"><subject>Human biology and medicine</subject></subj-group></article-categories><title-group><article-title>Global analysis of p53-regulated transcription identifies its direct targets and unexpected regulatory mechanisms</article-title></title-group><contrib-group><contrib contrib-type="author" id="author-10268"><name><surname>Allen</surname><given-names>Mary Ann</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="aff" rid="aff2"/><xref ref-type="aff" rid="aff3"/><xref ref-type="aff" rid="aff4"/><xref ref-type="other" rid="par-7"/><xref ref-type="fn" rid="con1"/><xref ref-type="fn" rid="conf2"/><xref ref-type="other" rid="dataro1"/></contrib><contrib contrib-type="author" equal-contrib="yes" id="author-10271"><name><surname>Andrysik</surname><given-names>Zdenek</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="aff" rid="aff2"/><xref ref-type="fn" rid="equal-contrib">†</xref><xref ref-type="fn" rid="con4"/><xref ref-type="fn" rid="conf2"/><xref ref-type="other" rid="dataro1"/></contrib><contrib contrib-type="author" equal-contrib="yes" id="author-10270"><name><surname>Dengler</surname><given-names>Veronica L</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="aff" rid="aff2"/><xref ref-type="fn" rid="equal-contrib">†</xref><xref ref-type="fn" rid="con5"/><xref ref-type="fn" rid="conf2"/><xref ref-type="other" rid="dataro1"/></contrib><contrib contrib-type="author" equal-contrib="yes" id="author-10269"><name><surname>Mellert</surname><given-names>Hestia S</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="aff" rid="aff2"/><xref ref-type="fn" rid="equal-contrib">†</xref><xref ref-type="other" rid="par-6"/><xref ref-type="fn" rid="con2"/><xref ref-type="fn" rid="conf2"/><xref ref-type="other" rid="dataro1"/></contrib><contrib contrib-type="author" id="author-10272"><name><surname>Guarnieri</surname><given-names>Anna</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="aff" rid="aff2"/><xref ref-type="fn" rid="con6"/><xref ref-type="fn" rid="conf2"/><xref ref-type="other" rid="dataro1"/></contrib><contrib contrib-type="author" id="author-10273"><name><surname>Freeman</surname><given-names>Justin A</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="aff" rid="aff2"/><xref ref-type="fn" rid="con7"/><xref ref-type="fn" rid="conf2"/><xref ref-type="other" rid="dataro1"/></contrib><contrib contrib-type="author" id="author-12992"><name><surname>Sullivan</surname><given-names>Kelly D</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="aff" rid="aff2"/><xref ref-type="fn" rid="con10"/><xref ref-type="fn" rid="conf2"/></contrib><contrib contrib-type="author" id="author-12993"><name><surname>Galbraith</surname><given-names>Matthew D</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="aff" rid="aff2"/><xref ref-type="fn" rid="con8"/><xref ref-type="fn" rid="conf2"/></contrib><contrib contrib-type="author" id="author-10274"><name><surname>Luo</surname><given-names>Xin</given-names></name><xref ref-type="aff" rid="aff5"/><xref ref-type="fn" rid="con11"/><xref ref-type="fn" rid="conf2"/><xref ref-type="other" rid="dataro1"/></contrib><contrib contrib-type="author" id="author-1420"><name><surname>Kraus</surname><given-names>W Lee</given-names></name><xref ref-type="aff" rid="aff5"/><xref ref-type="fn" rid="con12"/><xref ref-type="fn" rid="conf2"/><xref ref-type="other" rid="dataro1"/></contrib><contrib contrib-type="author" corresp="yes" id="author-10276"><name><surname>Dowell</surname><given-names>Robin D</given-names></name><xref ref-type="aff" rid="aff2"/><xref ref-type="aff" rid="aff3"/><xref ref-type="corresp" rid="cor1">*</xref><xref ref-type="other" rid="par-3"/><xref ref-type="other" rid="par-4"/><xref ref-type="fn" rid="con9"/><xref ref-type="fn" rid="conf2"/><xref ref-type="other" rid="dataro1"/></contrib><contrib contrib-type="author" corresp="yes" id="author-1082"><name><surname>Espinosa</surname><given-names>Joaquin M</given-names></name><xref ref-type="aff" rid="aff1"/><xref ref-type="aff" rid="aff2"/><xref ref-type="corresp" rid="cor2">*</xref><xref ref-type="other" rid="par-1"/><xref ref-type="other" rid="par-2"/><xref ref-type="other" rid="par-3"/><xref ref-type="other" rid="par-5"/><xref ref-type="fn" rid="con3"/><xref ref-type="fn" rid="conf1"/><xref ref-type="other" rid="dataro1"/></contrib><aff id="aff1"><institution content-type="dept">Howard Hughes Medical Institute</institution>, <institution>University of Colorado, Boulder</institution>, <addr-line><named-content content-type="city">Boulder</named-content></addr-line>, <country>United States</country></aff><aff id="aff2"><institution content-type="dept">Department of Molecular, Cellular and Developmental Biology</institution>, <institution>University of Colorado, Boulder</institution>, <addr-line><named-content content-type="city">Boulder</named-content></addr-line>, <country>United States</country></aff><aff id="aff3"><institution>BioFrontiers Institute</institution>, <addr-line><named-content content-type="city">Boulder</named-content></addr-line>, <country>United States</country></aff><aff id="aff4"><institution content-type="dept">Computational Biosciences Program</institution>, <institution>University of Colorado, Denver–Anschutz Medical Campus</institution>, <addr-line><named-content content-type="city">Aurora</named-content></addr-line>, <country>United States</country></aff><aff id="aff5"><institution>Signalling and Gene Regulation Laboratory, Cecil H and Ida Green Center for Reproductive Biology Sciences, University of Texas Southwestern Medical Center</institution>, <addr-line><named-content content-type="city">Dallas</named-content></addr-line>, <country>United States</country></aff></contrib-group><contrib-group content-type="section"><contrib contrib-type="editor"><name><surname>Prives</surname><given-names>Carol</given-names></name><role>Reviewing editor</role><aff><institution>Columbia University</institution>, <country>United States</country></aff></contrib></contrib-group><author-notes><corresp id="cor1"><label>*</label>For correspondence: <email>robin.dowell@colorado.edu</email> (RDD);</corresp><corresp id="cor2"><label>*</label>For correspondence: <email>joaquin.espinosa@colorado.edu</email> (JME)</corresp><fn fn-type="con" id="equal-contrib"><label>†</label><p>These authors contributed equally to this work</p></fn></author-notes><pub-date date-type="pub" publication-format="electronic"><day>27</day><month>05</month><year>2014</year></pub-date><pub-date pub-type="collection"><year>2014</year></pub-date><volume>3</volume><elocation-id>e02200</elocation-id><history><date date-type="received"><day>03</day><month>01</month><year>2014</year></date><date date-type="accepted"><day>25</day><month>04</month><year>2014</year></date></history><permissions><copyright-statement>© 2014, Allen et al</copyright-statement><copyright-year>2014</copyright-year><copyright-holder>Allen et al</copyright-holder><license xlink:href="http://creativecommons.org/licenses/by/3.0/"><license-p>This article is distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/3.0/">Creative Commons Attribution License</ext-link>, which permits unrestricted use and redistribution provided that the original author and source are credited.</license-p></license></permissions><self-uri content-type="pdf" xlink:href="elife02200.pdf"/><abstract><object-id pub-id-type="doi">10.7554/eLife.02200.001</object-id><p>The p53 transcription factor is a potent suppressor of tumor growth. We report here an analysis of its direct transcriptional program using Global Run-On sequencing (GRO-seq). Shortly after MDM2 inhibition by Nutlin-3, low levels of p53 rapidly activate ∼200 genes, most of them not previously established as direct targets. This immediate response involves all canonical p53 effector pathways, including apoptosis. Comparative global analysis of RNA synthesis vs steady state levels revealed that microarray profiling fails to identify low abundance transcripts directly activated by p53. Interestingly, p53 represses a subset of its activation targets before MDM2 inhibition. GRO-seq uncovered a plethora of gene-specific regulatory features affecting key survival and apoptotic genes within the p53 network. p53 regulates hundreds of enhancer-derived RNAs. Strikingly, direct p53 targets harbor pre-activated enhancers highly transcribed in p53 null cells. Altogether, these results enable the study of many uncharacterized p53 target genes and unexpected regulatory mechanisms.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.02200.001">http://dx.doi.org/10.7554/eLife.02200.001</ext-link></p></abstract><abstract abstract-type="executive-summary"><object-id pub-id-type="doi">10.7554/eLife.02200.002</object-id><title>eLife digest</title><p>The growth, division and eventual death of the cells in the body are processes that are tightly controlled by hundreds of genes working together. If any of these genes are switched on (or off) in the wrong cell or at the wrong time, it can lead to cancer.</p><p>It has been known for many years that the protein encoded by one gene in particular—called <italic>p53</italic>—is nearly always switched off in cancer cells. The p53 protein normally acts like a ‘brake’ to slow the uncontrolled division of cells, and some researchers are working to find ways to switch on this protein in cancer cells. However, this approach appears to only work in specific cases of this disease. For better results, we need to understand how p53 is normally switched on, and what other genes this protein controls once it is activated.</p><p>Allen et al. have now identified the genes that are directly switched on when cancer cells are treated with a drug that artificially activates the p53 protein. Nearly 200 genes were switched on, and almost three quarters of these genes had not previously been identified as direct targets of p53.</p><p>Although p53 tends to act as a brake to slow cell division, it is not clear how it distinguishes between its target genes—some of which promote cell survival, while others promote cell death. Allen et al. found that survival genes are switched on more strongly than cell death genes via a range of different mechanisms; this may explain why most cancers can survive drug treatments that reactivate p53. Also, Allen et al. revealed that some p53 target genes are primed to be switched on, even before the p53 protein is activated, by proteins (and other molecules) acting in regions of the DNA outside of the genes.</p><p>By uncovering many new gene targets for the p53 protein, the findings of Allen et al. could help researchers developing new drugs or treatments for cancer.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.02200.002">http://dx.doi.org/10.7554/eLife.02200.002</ext-link></p></abstract><kwd-group kwd-group-type="author-keywords"><title>Author keywords</title><kwd>tumor supressor</kwd><kwd>genomics</kwd><kwd>p21</kwd><kwd>PUMA</kwd><kwd>PIG3</kwd><kwd>eRNA</kwd></kwd-group><kwd-group kwd-group-type="research-organism"><title>Research organism</title><kwd>human</kwd></kwd-group><funding-group><award-group id="par-1"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000011</institution-id><institution>Howard Hughes Medical Institute</institution></institution-wrap></funding-source><award-id>Early Career Award</award-id><principal-award-recipient><name><surname>Espinosa</surname><given-names>Joaquin M</given-names></name></principal-award-recipient></award-group><award-group id="par-2"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000002</institution-id><institution>National Institutes of Health</institution></institution-wrap></funding-source><award-id>RO1 CA117907-07</award-id><principal-award-recipient><name><surname>Espinosa</surname><given-names>Joaquin M</given-names></name></principal-award-recipient></award-group><award-group id="par-3"><funding-source><institution-wrap><institution>Butcher Foundation</institution></institution-wrap></funding-source><principal-award-recipient><name><surname>Dowell</surname><given-names>Robin D</given-names></name><name><surname>Espinosa</surname><given-names>Joaquin M</given-names></name></principal-award-recipient></award-group><award-group id="par-4"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100005508</institution-id><institution>Boettcher Foundation</institution></institution-wrap></funding-source><principal-award-recipient><name><surname>Dowell</surname><given-names>Robin D</given-names></name></principal-award-recipient></award-group><award-group id="par-5"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000001</institution-id><institution>National Science Foundation</institution></institution-wrap></funding-source><award-id>MCB1243522</award-id><principal-award-recipient><name><surname>Espinosa</surname><given-names>Joaquin M</given-names></name></principal-award-recipient></award-group><award-group id="par-6"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000048</institution-id><institution>American Cancer Society</institution></institution-wrap></funding-source><principal-award-recipient><name><surname>Mellert</surname><given-names>Hestia S</given-names></name></principal-award-recipient></award-group><award-group id="par-7"><funding-source><institution-wrap><institution-id institution-id-type="FundRef">http://dx.doi.org/10.13039/100000002</institution-id><institution>National Institutes of Health</institution></institution-wrap></funding-source><award-id>2T15 LM009451</award-id><principal-award-recipient><name><surname>Allen</surname><given-names>Mary Ann</given-names></name></principal-award-recipient></award-group><funding-statement>The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.</funding-statement></funding-group><custom-meta-group><custom-meta><meta-name>elife-xml-version</meta-name><meta-value>2</meta-value></custom-meta><custom-meta specific-use="meta-only"><meta-name>Author impact statement</meta-name><meta-value>The first global analysis of the effects of the transcription factor p53 on RNA synthesis illuminates several novels aspects of the molecular mechanism of action of this potent tumor suppressor.</meta-value></custom-meta></custom-meta-group></article-meta></front><body><sec id="s1" sec-type="intro"><title>Introduction</title><p>The p53 transcription factor is activated by potentially oncogenic stimuli such as ribosomal stress, DNA damage, telomere erosion, nutrient deprivation and oncogene hyperactivation (<xref ref-type="bibr" rid="bib45">Vousden and Prives, 2009</xref>). In the absence of activating signals, p53 is repressed by the oncoproteins MDM2 and MDM4. MDM2 masks the transactivation domain of p53 and is also an E3 ligase that targets p53 for degradation (<xref ref-type="bibr" rid="bib27">Momand et al., 1992</xref>; <xref ref-type="bibr" rid="bib31">Oliner et al., 1993</xref>; <xref ref-type="bibr" rid="bib20">Kubbutat et al., 1997</xref>). MDM4 lacks E3 ligase activity, but represses p53 transactivation potential (<xref ref-type="bibr" rid="bib35">Riemenschneider et al., 1999</xref>). Diverse signaling pathways converge on the p53/MDM2/MDM4 complex to release p53 from its repressors and enable it to regulate transcription of downstream target genes involved in cellular responses such as cell cycle arrest, apoptosis, senescence, autophagy, DNA repair and central metabolism (<xref ref-type="bibr" rid="bib45">Vousden and Prives, 2009</xref>). p53 is inactivated in virtually all human cancers, either by mutations in its DNA binding domain or MDM2/MDM4 overexpression. Significant advances have been made to develop p53-based targeted therapies (<xref ref-type="bibr" rid="bib4">Brown et al., 2009</xref>). One class of small molecules targets the interaction between p53 and its repressors, thus bypassing the need of stress signaling to trigger p53 activation. For example, Nutlin-3, the first-in-class compound, binds to the hydrophobic pocket in MDM2 required for binding to p53, thus acting as a competitive inhibitor (<xref ref-type="bibr" rid="bib44">Vassilev et al., 2004</xref>). A second class of molecules binds to mutant p53 and partially restores its wild type function (<xref ref-type="bibr" rid="bib4">Brown et al., 2009</xref>). As these compounds enter clinical trials, their efficacy is limited by the fact that p53 activation leads to cancer cell death only in specific scenarios. Thus, there is a clear need to understand how these molecules modulate p53 function and how cell fate choice upon p53 activation is defined. A missing piece in this effort is a definitive elucidation of the direct p53 transcriptome.</p><p>Despite its unequivocal importance in cancer biology, our understanding of p53 function as a transcription factor is limited. The protein domains required for DNA binding and transactivation are well characterized, as well as its DNA response elements (p53REs) (<xref ref-type="bibr" rid="bib22">Laptenko and Prives, 2006</xref>). A recent comprehensive survey of the literature identified ∼120 genes for which direct regulation has been established (<xref ref-type="bibr" rid="bib36">Riley et al., 2008</xref>), but a comprehensive analysis of p53-regulated RNAs is still missing. Up to this point, the global p53 transcriptional response has been investigated with techniques that measure steady state RNA levels, mostly microarray profiling. These methods require long time points to observe a significant change in the expression of p53-regulated RNAs, which confounds direct vs indirect effects, and additional experiments are required to ascertain direct transcriptional regulation. A popular approach has been to cross-reference microarray data with p53 binding data derived from ChIP-seq assays. Meta-analysis of four recent papers using this strategy indicates that p53 may directly activate >1200 genes, yet only 26 of these genes were commonly activated in all four studies (<xref ref-type="bibr" rid="bib29">Nikulenkov et al., 2012</xref>; <xref ref-type="bibr" rid="bib26">Menendez et al., 2013</xref>; <xref ref-type="bibr" rid="bib37">Schlereth et al., 2013</xref>; <xref ref-type="bibr" rid="bib46">Wang et al., 2013</xref>) (see later, <xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1</xref>). It is unclear to what extent this lack of overlap is due to methodological differences and/or cell type-specific differences in direct p53 action vs post-transcriptional regulation.</p><p>We report here the first genome-wide analysis of p53-regulated RNA synthesis. Using Global Run-On sequencing (GRO-seq) (<xref ref-type="bibr" rid="bib7">Core et al., 2008</xref>), we ascertained direct regulation by using a short time point of Nutlin-3 treatment in isogenic cell lines with or without p53. Strikingly, Nutlin leads to p53-dependent transcriptional activation of hundreds of genomic loci prior to any significant increase in total p53 levels, thus highlighting the critical role of MDM2 in masking the p53 transactivation domain. Comparative global analysis of RNA synthesis by GRO-seq vs RNA steady state levels by microarray revealed that many p53 target genes transcribed at low levels are missed by microarray experiments. Strikingly, p53 represses the basal expression of a subset of its target genes before MDM2 inhibition. GRO-seq uncovered many gene-specific transcriptional events affecting key survival and apoptotic genes within the network, including the occurrence of bidirectional promoters (both overlapping and non-overlapping), clustered transactivation and various forms of antisense transcription. GRO-seq revealed widespread activation of enhancer-derived RNAs (eRNAs) arising from p53REs. Interestingly, direct p53 target genes harbor ‘pre-activated’ p53REs, as defined by the strong production of eRNAs in the isogenic p53 null cells. These results elucidate novel p53-regulated RNAs as well as gene-specific regulatory events within the p53 network and pave the road for a myriad of future mechanistic studies.</p></sec><sec id="s2" sec-type="results"><title>Results</title><sec id="s2-1"><title>p53 rapidly activates a multifunctional transcriptional program upon MDM2 inhibition</title><p>In order to study the direct transcriptional response upon p53 activation, we performed GRO-seq in isogenic HCT116 p53 +/+ and p53 −/− cell lines. After several hours of treatment with 10 μM Nutlin-3a (referred hereto as Nutlin), only p53 +/+ cells undergo cell cycle arrest and display induction of many known p53 target genes (<xref ref-type="bibr" rid="bib43">Tovar et al., 2006</xref>; <xref ref-type="bibr" rid="bib16">Henry et al., 2012</xref>). Using GRO-seq, we investigated the effects of Nutlin vs vehicle (DMSO) in both cell lines under conditions of exponential cell proliferation, when p53 levels are low and cell cultures display virtually no signs of cell cycle arrest or apoptosis. To minimize the possibility of measuring indirect effects, we chose a 1 hr time point of Nutlin treatment. Typical GRO-seq results are shown for two well-characterized p53 target genes, <italic>CDKN1A</italic> (<italic>p21</italic>) and <italic>TP53I3</italic> (<italic>PIG3</italic>) (<xref ref-type="bibr" rid="bib9">el-Deiry et al., 1993</xref>; <xref ref-type="bibr" rid="bib33">Polyak et al., 1997</xref>), which display significant increase in GRO-seq signals upon Nutlin treatment only in HCT116 p53 +/+ cells (<xref ref-type="fig" rid="fig1">Figure 1A</xref>, <xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1A</xref>, respectively). Not surprisingly, Q-RT-PCR demonstrates that the steady state expression of the CDKN1A and TP53I3 mRNAs does not increase at the 1 hr time point as used for GRO-seq (<xref ref-type="fig" rid="fig1">Figure 1B</xref>). In fact, a significant increase in the steady state levels of both mRNAs requires several hours of p53 activation. Furthermore, Western blot analysis shows that 1 hr of Nutlin treatment does not increase total p53 or p21 protein levels to a significant degree (<xref ref-type="fig" rid="fig1">Figure 1C</xref>). Analysis of cell cycle progression using BrdU incorporation assays shows no signs of cell cycle arrest at the 1 hr time point, but a clear reduction in S phase cells is evident at 12 hr only in the p53 +/+ cells (<xref ref-type="fig" rid="fig1">Figure 1C</xref>, <xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1B</xref>). Overall, these observations indicate that our GRO-seq analysis using a 1 hr time point of p53 activation would largely preclude secondary effects driven by well-established direct p53 target genes, such as <italic>CDKN1A</italic> and <italic>TP53I3</italic>, thus enabling the identification of the direct p53 transcriptome.<fig-group><fig id="fig1" position="float"><object-id pub-id-type="doi">10.7554/eLife.02200.003</object-id><label>Figure 1.</label><caption><title>GRO-seq analysis of the p53 transcriptional program.</title><p>(<bold>A</bold>) GRO-seq results for the p53 target locus <italic>CDKN1A</italic> (<italic>p21</italic>). Isogenic p53 −/− and p53 +/+ HCT116 cells were treated for 1 hr with either 10 μM Nutlin-3a (Nutlin) or vehicle (DMSO, Control). Fragments per kilobase per million reads (fpkm) are shown for the intragenic region. The first kilobase downstream of the transcription start site (TSS) was excluded from the fpkm calculation to minimize effects of RNAPII pausing. The total genomic region displayed is indicated in the top left corner. Blue signals are reads mapping to the sense strand, red signals are reads mapping to the antisense strand. See <xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1A</xref> for results of the <italic>TP53I3</italic> locus. (<bold>B</bold>) GRO-seq detects transactivation of the canonical p53 target genes <italic>CDKN1A</italic> and <italic>TP53I3</italic> at 1 hr of Nutlin treatment, prior to any detectable increase in steady state mRNA levels as measured by Q-RT-PCR. (<bold>C</bold>) A 1 hr time point of Nutlin treatment does not produce significant p53 accumulation, p21 protein induction or a decrease in number of S phase cells as measured by BrdU incorporation assays. * indicates p<0.05. See also <xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1B</xref> for quantification data of BrdU assays. (<bold>D</bold>) Genome-wide analysis using the DESeq algorithm identifies 198 annotated gene loci transactivated upon Nutlin treatment only in HCT116 p53 +/+ cells. See <xref ref-type="supplementary-material" rid="SD1-data">Supplementary file 1</xref> for a detailed annotation of these genes. (<bold>E</bold>) Q-RT-PCR validates induction of novel and predicted direct p53 target genes upon 12 hr of Nutlin treatment. mRNA expression was normalized to 18s rRNA values and expressed as fold change Nutlin/DMSO. Data shown are the average of three biological replicates with standard errors from the mean. (<bold>F</bold>) Flow cytometry analysis using the DO-1 antibody recognizing the MDM2-binding surface in the p53 transcactivation domain 1 (TAD1) reveals increased reactivity as early as 1 hr of Nutlin treatment, indicative of unmasking of the TAD1 at this early time point. (<bold>G</bold>) p53 directly activates a multifunctional transcriptional program at 1 hour of Nutlin treatment, including many canonical apoptotic genes. See <xref ref-type="supplementary-material" rid="SD1-data">Supplementary file 1</xref> for a complete list and annotation.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.02200.003">http://dx.doi.org/10.7554/eLife.02200.003</ext-link></p></caption><graphic xlink:href="elife02200f001"/></fig><fig id="fig1s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.02200.004</object-id><label>Figure 1—figure supplement 1.</label><caption><title>GRO-seq reveals the immediate direct p53 transcriptional response.</title><p>(<bold>A</bold>) GRO-seq results for the well characterized pro-apoptotic p53 target gene <italic>TP53I3</italic> (<italic>PIG3</italic>). (<bold>B</bold>). BrdU incorporation assays for HCT116 p53 −/− and p53 +/+ cells upon Nutlin treatment for the 1 hr time point employed for GRO-seq and the 12 hr time point used for microarray experiments. One hour Nutlin does not induce significant cell cycle arrest. (<bold>C</bold>) GRO-seq shows that 176 of the 198 genes transactivated upon 1 hr of Nutlin are already increased above 1x as early as 30 min after addition of the drug to cell cultures. (<bold>D</bold>) Of the 198 genes transactivated by p53 at 1 hr of Nutlin treatment, 55 were known direct targets, 66 were predicted by published ChIP-seq/microarray studies and 77 are novel. See also <xref ref-type="supplementary-material" rid="SD1-data">Supplementary file 1</xref>. (<bold>E</bold>) Q-RT-PCR shows that most novel p53 targets identified by GRO-seq are also induced by doxorubicin, a genotoxic p53-activating agent. mRNA expression was normalized to 18s rRNA values and is expressed as fold change doxorubicin/DMSO. Data shown are the average of three biological replicates with standard errors from the mean. (<bold>F</bold>) Western blot assays show induction of the novel p53 target genes GRIN2C, PTCDH4 and RINL at the protein level upon Nutlin treatment of HCT116 p53 +/+ cells. Nucleolin and α-tubulin serve as loading controls. Asterisk denotes a non-specific band in the RINL blot. (<bold>G</bold>–<bold>H</bold>) Immunofluroescence staining of p53 shows a similar pattern for control cells (0 hr) and cells treated with Nutlin for 1 hr. Upon long exposure to Nutlin (12 hr), almost the entire cellular population exhibits strong nuclear signal for p53 staining. In H, 'high p53 staining' was defined as the average signal observed in the 12 hr Nutlin cell culture.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.02200.004">http://dx.doi.org/10.7554/eLife.02200.004</ext-link></p></caption><graphic xlink:href="elife02200fs001"/></fig></fig-group></p><p>Next, we used the DESeq algorithm (<xref ref-type="bibr" rid="bib1">Anders and Huber, 2010</xref>) to identify annotated gene loci whose transcription is significantly changed upon Nutlin treatment (see ‘Materials and methods’ for details). Using a cut-off of adjusted <italic>p(a)</italic> <0.1, we identified 198 gene loci whose transcription is significantly induced upon Nutlin treatment in p53 +/+ cells (<xref ref-type="fig" rid="fig1">Figure 1D</xref>; <xref ref-type="supplementary-material" rid="SD1-data">Supplementary file 1</xref>). This analysis identified only four gene loci whose transcription was diminished in the p53 +/+ cells (<italic>FLVCR2</italic>, <italic>NR4A3</italic>, <italic>RELB</italic> and <italic>EGR1</italic>); however, none of these genes showed reductions in steady state mRNA levels upon prolonged p53 activation (see later, <xref ref-type="fig" rid="fig2">Figure 2</xref>). The specificity of Nutlin is demonstrated by the negligible changes observed in p53 −/− cells, where our analysis identified 5 induced and 2 repressed genes, all of which have less than 1.5-fold changes and none of which was among those differentially transcribed in p53 +/+ cells (<xref ref-type="fig" rid="fig1">Figure 1D</xref>). From this point forth, we focused on the 198 genes activated in the p53 +/+ cells, which we considered to be the direct p53 transcriptional program in this cell type. The notion that these genes are indeed direct p53 targets is reinforced by the observation that most of them (176 out of 198) show an increase in transcription as early as 30 min after Nutlin addition to the cell culture (<xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1C</xref>). Of these 198 genes, 55 were known validated direct p53 targets, 66 were targets predicted by one or more published microarray/ChIP-seq studies, and 77 are putative novel direct p53 targets (<xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1D</xref>, a comprehensive annotation of these genes is provided in <xref ref-type="supplementary-material" rid="SD1-data">Supplementary file 1</xref>). Q-RT-PCR validation showed that novel genes are induced at a 12 hr time point of Nutlin treatment at the mRNA steady state level to a degree comparable to those genes predicted by published microarray/ChIP-seq studies (<xref ref-type="fig" rid="fig1">Figure 1E</xref>). Furthermore, 12 out of the 14 novel p53 target genes tested are also induced at the mRNA steady state level when using doxorubicin, a DNA-damaging agent that activates p53 via stress signaling cascades (<xref ref-type="bibr" rid="bib24">Lowe et al., 1994</xref>), thus revealing that transactivation of most novel genes is not unique to pharmacological inhibition of MDM2 (<xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1E</xref>). Finally, we investigated whether activation of novel p53 targets can also be observed at the protein level. Indeed, Western blot analysis demonstrates protein induction for the novel genes GRIN2C, PTCDH4 and RINL (<xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1F</xref>). Thus, our GRO-seq experiment clearly expands the universe of direct p53 target genes, paving the road for mechanistic studies investigating the function of these genes in the p53 network.<fig-group><fig id="fig2" position="float"><object-id pub-id-type="doi">10.7554/eLife.02200.005</object-id><label>Figure 2.</label><caption><title>Global analysis of p53 effects on RNA synthesis vs steady state levels.</title><p>(<bold>A</bold>) MAplots for GRO-seq and microarray gene profiling experiments in HCT116 p53 +/+ cells after 1 hr and 12 hr of Nutlin treatment, respectively. Colors indicate whether genes scored as statistically different in both platforms (purple), in the GRO-seq only (red) or the microarray experiment only (blue). (<bold>B</bold>) Few genes downregulated in the microarray experiment show p53 binding within 25 kb of the gene, suggestive of indirect regulation. (<bold>C</bold>) Bubble plots displaying relative signals derived from the GRO-seq and microarray experiments illustrate how genes with very high basal expression or very low transcription are not significantly affected at the steady state level as measured by microarray. For the <italic>CDC42BPG</italic>, <italic>KLHDC7A</italic>, <italic>ADAMTS7</italic>, <italic>LRP1</italic> and <italic>ASTN2</italic> loci, the signals were replotted at 25-fold magnification. (<bold>D</bold>) Scatter plot showing comparative fold induction for p53 target genes transactivated at 1 hr Nutlin treatment between the GRO-seq and microarray experiments. (<bold>E</bold>) Q-RT-PCR indicates that many low abundance transcripts upregulated by GRO-seq are indeed induced at the steady state level. (<bold>F</bold>) Box and whisker plots showing the expression of various gene sets as detected by microarray.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.02200.005">http://dx.doi.org/10.7554/eLife.02200.005</ext-link></p></caption><graphic xlink:href="elife02200f002"/></fig><fig id="fig2s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.02200.006</object-id><label>Figure 2—figure supplement 1.</label><caption><title>Mechanisms of indirect gene repression by p53.</title><p>(<bold>A</bold> and <bold>B</bold>) Meta-analysis of four recent investigations of the p53 transcriptional program using microarray analysis of RNA steady state levels and ChIP-seq measurements of p53 binding reveals little overlap between the experiments. See <xref ref-type="supplementary-material" rid="SD2-data">Supplementary file 2</xref> for the various gene lists used to generate the Venn diagrams using Venny (<ext-link ext-link-type="uri" xlink:href="http://bioinfogp.cnb.csic.es/tools/venny/">http://bioinfogp.cnb.csic.es/tools/venny/</ext-link>). *This study employed SAOS2 cells, which are p53 null, with overexpression of various p53 isoforms. ** This study employed HCT116 p53 −/− cells with overexpression of natural p53 polymorphic variants. (<bold>C</bold>) 72% of genes downregulated upon 12 hr of Nutlin treatment in HCT116 cells were found to be repressed in this same cell type by overexpression of miR-34a, a p53 inducible miRNA. (<bold>D</bold>) Ingenuity Pathway Analysis of genes downregulated upon 12 hr Nutlin treatment indicates that the top three regulators of this gene set are E2F4, CDKN1A and RB. (<bold>E</bold>) Model of indirect gene repression by p53 via upregulation of <italic>CDKN1A</italic> and miR-34a, both of which inhibited G1-S CDK complexes.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.02200.006">http://dx.doi.org/10.7554/eLife.02200.006</ext-link></p></caption><graphic xlink:href="elife02200fs002"/></fig><fig id="fig2s2" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.02200.007</object-id><label>Figure 2—figure supplement 2.</label><caption><title>ChIP analysis of novel p53 target genes.</title><p>Analysis of p53 ChIP-seq datasets indicated high confidence p53 binding events in the vicinity of novel p53 target genes identified by GRO-seq at the genomic positions labeled with a red asterisk. ChIP-Q-PCR assays in HCT116 p53 +/+ cells confirm p53 binding at these locations above background levels (IgG control).</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.02200.007">http://dx.doi.org/10.7554/eLife.02200.007</ext-link></p></caption><graphic xlink:href="elife02200fs003"/></fig></fig-group></p><p>Although it is known that MDM2 represses p53 by both masking its transactivation domain and also targeting it for degradation (<xref ref-type="bibr" rid="bib27">Momand et al., 1992</xref>; <xref ref-type="bibr" rid="bib31">Oliner et al., 1993</xref>; <xref ref-type="bibr" rid="bib20">Kubbutat et al., 1997</xref>), it has been difficult to dissect to what extent each mechanism contributes to repression of p53 target genes in diverse functional categories. Studies employing steady state mRNA measurements concluded that prolonged p53 activation and/or higher levels of cellular p53 were required for activation of apoptotic genes, some of which display delayed kinetics of induction at the mRNA steady state level as compared to cell cycle arrest genes (<xref ref-type="bibr" rid="bib5">Chen et al., 1996</xref>; <xref ref-type="bibr" rid="bib50">Zhao et al., 2000</xref>; <xref ref-type="bibr" rid="bib40">Szak et al., 2001</xref>; <xref ref-type="bibr" rid="bib11">Espinosa et al., 2003</xref>; <xref ref-type="bibr" rid="bib8">Das et al., 2007</xref>). However, GRO-seq demonstrates that a 1 hr time point of Nutlin treatment induces transcription of genes in every major pathway downstream of p53 (<xref ref-type="supplementary-material" rid="SD1-data">Supplementary file 1</xref>). The observation that key survival and apoptotic genes (e.g., <italic>CDKN1A</italic>, <italic>TP53I3</italic>) show greater than sixfold increase in transcription at a time point preceding a proportional increase in total p53 levels (<xref ref-type="fig" rid="fig1">Figure 1A,C</xref>, <xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1A</xref>), suggests that the mere unmasking of the p53 transactivation domain suffices to activate a multifaceted transcriptional program. To further test this notion, we performed flow cytometry analyses using a monoclonal antibody (DO-1) that recognizes an epitope in the p53 N-terminal transactivation domain 1 (TAD1) that overlaps with the MDM2-binding surface competed by Nutlin (<xref ref-type="bibr" rid="bib32">Picksley et al., 1994</xref>). In fact, the DO-1 antibody competes the p53-MDM2 interaction in vitro in analogous fashion to Nutlin (<xref ref-type="bibr" rid="bib6">Cohen et al., 1998</xref>). Under the denaturing conditions of a Western Blot assay, where p53-MDM2 complexes are fully disrupted, this antibody shows no significant increase in total p53 levels at the 1 hr time point of Nutlin treatment (<xref ref-type="fig" rid="fig1">Figure 1C</xref>). However, we posited that DO-1 reactivity should increase significantly upon Nutlin treatment under the fixed conditions employed in flow cytometry. Expectedly, flow cytometry quantitation shows that, even before Nutlin treatment, p53 +/+ cells have significantly more DO-1 reactivity than p53 −/− cells (<xref ref-type="fig" rid="fig1">Figure 1F</xref>). The functional importance of this ‘basal p53 activity’ will be investigated later in this report (<xref ref-type="fig" rid="fig3">Figure 3</xref>). Interestingly, the p53 +/+ cell population shifts to significantly higher DO-1 reactivity at the 1 hr time point, as predicted by epitope unmasking. A further increase is observed at 12 hr of Nutlin treatment, when total p53 levels have risen considerably as measured by Western blots (<xref ref-type="fig" rid="fig1">Figure 1C,F</xref>). Finally, since GRO-seq is a population average experiment, we performed immunofluorescence assays to test if our GRO-seq results could be explained by massive p53 accumulation in just a few outlier cells within the population at the 1 hr time point. However, these experiments discarded the notion of outlier cells: although ∼3% cells show high p53 staining at the 1 hr time point, this number is not significantly different than observed in control p53 +/+ cells (<xref ref-type="fig" rid="fig1s1">Figure 1—figure supplement 1G,H</xref>).<fig-group><fig id="fig3" position="float"><object-id pub-id-type="doi">10.7554/eLife.02200.008</object-id><label>Figure 3.</label><caption><title>p53 exerts varying activating and repressing effects on its target genes prior to MDM2 inhibition.</title><p>(<bold>A</bold>) 198 genes activated upon 1 hr Nutlin treatment in HCT116 p53 +/+ cells are ranked from left to right based on their basal transcription in p53 +/+ cells over p53 −/− cells. Green indicates genes whose basal transcription is greater than twofold in p53 +/+ cells, red indicates lesser than twofold. Grey dots display the transcription of the same genes in Nutlin-treated p53 +/+ cells. (<bold>B</bold>) Heatmap displaying relative transcriptional activity of direct p53 target genes identified by GRO-seq relative to control p53 −/− cells. Genes are sorted based on their transcription in control p53 +/+ cells. (<bold>C</bold>) Genome browser views of representative genes whose basal transcription is higher (<italic>GDF15</italic>) or lower (<italic>PTP4A1</italic>) in the presence of MDM2-bound p53. See <xref ref-type="fig" rid="fig3s1">Figure 3—figure supplement 1A</xref> for matching RNAPII ChIP data. (<bold>D</bold>) Q-RT-PCR measurements of genes whose basal transcription was found to be 2x higher (green) or lower (red) in the presence of MDM2-bound p53. (<bold>E</bold>) ChIP assays show binding of p53 and MDM2 to the p53REs in the <italic>CDKN1A</italic> and <italic>PTP4A1</italic> gene loci (−2283 bp and +1789 relative to TSS, respectively), prior to inhibition of the p53-MDM2 interaction by Nutlin. Nutlin treatment leads to increased p53 signals with the DO-1 antibody recognizing the p53 TAD1, concurrently with a decrease in MDM2 signals. MDM2 ChIP was performed in SJSA cells carrying a MDM2 gene amplification <bold>F</bold>. Oncomine gene expression analysis of 598 cancer cell lines of varied p53 status shows that <italic>CDKN1A</italic>, <italic>DDB2</italic> and <italic>GDF15</italic> are more highly expressed in wild type p53 cell lines, whereas <italic>GJB5</italic> is more highly expressed in mutant p53 cell lines. The ranking position of these genes is also indicated.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.02200.008">http://dx.doi.org/10.7554/eLife.02200.008</ext-link></p></caption><graphic xlink:href="elife02200f003"/></fig><fig id="fig3s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.02200.009</object-id><label>Figure 3—figure supplement 1.</label><caption><title>Differential effects of p53 on the basal transcription of its target genes.</title><p>(<bold>A</bold>) ChIP analysis of the ‘basally activated’ <italic>GDF15</italic> and ‘basally repressed’ <italic>PTP4A1</italic> gene loci using antibodies against the Serine 5- and Serine 2-phosphorylated forms of the C-terminal domain repeats of RBP1, the largest subunit of RNAPII. (<bold>B</bold>) ChIP assays for p53 and MDM2 at the p53RE in the <italic>HES2</italic> locus (+5160 relative to TSS) performed as in <xref ref-type="fig" rid="fig3">Figure 3E</xref>. (<bold>C</bold>) Western blots for MDM2 using the protein extracts from SJSA cells employed for MDM2 ChIP assays. Long exposures show detectable amounts of MDM2 in both the input and MDM2-ChIP samples. A 12 hr or Nutlin treatment was included as positive control.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.02200.009">http://dx.doi.org/10.7554/eLife.02200.009</ext-link></p></caption><graphic xlink:href="elife02200fs004"/></fig><fig id="fig3s2" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.02200.010</object-id><label>Figure 3—figure supplement 2.</label><caption><title>p53 mutational status affects the basal expression of its target genes.</title><p>(<bold>A</bold>) Gene set enrichment analysis (GSEA) of direct p53 target genes identified by GRO-seq (left) or genes upregulated upon Nutlin treatment in the microarray experiment only (right) relative to a ranked list of 12,624 genes analyzed via Oncomine for their relative expression in WT p53 vs mutant p53 cell lines. (<bold>B</bold>) Scatter plot comparing relative transcription as measured by GRO-seq in HCT116 p53 +/+ over p53 −/− cells vs relative mRNA expression in p53 WT vs p53 mutant cell lines. Red dots are genes whose transcription was lesser than twofold in p53 +/+ cells, green dots are genes whose transcription was greater than twofold in p53 +/+ cells, all other direct p53 targets in HCT116 cells are highlighted in black. (<bold>C</bold>) Oncomine gene expression analysis of 247 breast carcinomas of varied p53 status shows that <italic>CDKN1A</italic>, <italic>DDB2</italic> and <italic>GDF15</italic> are more highly expressed in wild type p53 tumors, whereas <italic>GJB5</italic> is more highly expressed in mutant p53 tumors. (<bold>D</bold>) Model depicting the differential effects of MDM2-p53 complexes on the basal expression of p53 target genes.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.02200.010">http://dx.doi.org/10.7554/eLife.02200.010</ext-link></p></caption><graphic xlink:href="elife02200fs005"/></fig></fig-group></p><p>Altogether, these results indicate that the low levels of p53 present in proliferating cancer cells suffice to directly activate a multifunctional transcriptional program, including many canonical apoptotic genes, upon unmasking of the p53 transactivation domain by Nutlin. However, as discussed later in the paper (<xref ref-type="fig" rid="fig4">Figure 4</xref>), this conclusion does not necessarily conflict with previous reports showing differential timing of mRNA accumulation between arrest and apoptotic genes as seen by steady state RNA measurements.<fig-group><fig id="fig4" position="float"><object-id pub-id-type="doi">10.7554/eLife.02200.011</object-id><label>Figure 4.</label><caption><title>GRO-seq analysis of key survival and death genes within the p53 network.</title><p>(<bold>A</bold>) The 10 most transcribed pro-survival and pro-apoptotic genes identified by GRO-seq ranked by decreasing transcriptional output in Nutlin-treated p53 +/+ cells. The surface of the bubbles represents the GRO-seq signal output relative to the <italic>CDKN1A</italic> locus. (<bold>B</bold>) Transcriptional output of same genes shown in <bold>A</bold> in p53 −/− cells. (<bold>C</bold>) Fold change analysis showing the overall effect of p53 on the transcription of its survival and apoptotic targets. (<bold>D</bold>) Survival genes within the p53 network tend to carry more proximally bound, transcriptionally engaged RNAPII over their promoter regions than apoptotic genes.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.02200.011">http://dx.doi.org/10.7554/eLife.02200.011</ext-link></p></caption><graphic xlink:href="elife02200f004"/></fig><fig id="fig4s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.02200.012</object-id><label>Figure 4—figure supplement 1.</label><caption><title>p53 target genes display a wide range of RNAPII pausing and promoter divergence.</title><p>(<bold>A</bold>) RNAPII pausing is not a pre-requisite for rapid activation among p53 target genes. A ranking of pausing indices show no correlation with fold induction among p53 targets. <italic>POLH</italic> and <italic>PHLDA3</italic> are representative examples of genes with high and low pausing indices. (<bold>B</bold>) p53 target genes show a great variation in the amount of promoter divergence, but a ranking of divergence indices shows no correlation with fold induction. <italic>DRAM1</italic> and <italic>DDB2</italic> are examples of p53 targets with high and low divergence, respectively.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.02200.012">http://dx.doi.org/10.7554/eLife.02200.012</ext-link></p></caption><graphic xlink:href="elife02200fs006"/></fig><fig id="fig4s2" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.02200.013</object-id><label>Figure 4—figure supplement 2.</label><caption><title>Examples of gene-specific features affecting key pro-apoptotic and survival p53 target genes.</title><p>See main text for details.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.02200.013">http://dx.doi.org/10.7554/eLife.02200.013</ext-link></p></caption><graphic xlink:href="elife02200fs007"/></fig></fig-group></p></sec><sec id="s2-2"><title>Global analysis of p53 effects on RNA synthesis vs RNA steady state levels</title><p>The global p53 transcriptional response has been previously investigated using measurements of RNA steady state levels (i.e., microarray profiling) and p53 chromatin binding (e.g., ChIP-seq). Meta-analysis of four recent reports using this approach indicates that >1200 genes are putative direct targets of p53 transactivation, yet only 26 are common between the four studies (<xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1A,B</xref>; <xref ref-type="supplementary-material" rid="SD2-data">Supplementary file 2</xref>) (<xref ref-type="bibr" rid="bib29">Nikulenkov et al., 2012</xref>; <xref ref-type="bibr" rid="bib26">Menendez et al., 2013</xref>; <xref ref-type="bibr" rid="bib37">Schlereth et al., 2013</xref>; <xref ref-type="bibr" rid="bib46">Wang et al., 2013</xref>). Furthermore, these studies suggest 80 genes that could be directly repressed by p53, yet none are shared between any two studies (<xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1A,B</xref>; <xref ref-type="supplementary-material" rid="SD2-data">Supplementary file 2</xref>). In order to investigate how GRO-seq analysis of the immediate p53 transcriptional response would compare to a global analysis of RNA steady state levels, we performed a microarray analysis of HCT116 p53 +/+ cells after 12 hr of Nutlin treatment, a time point similar to that used in the previous studies. Several important observations arise from this comparison.</p><p>First, there is a clear lack of overlap between the two analyses (<xref ref-type="fig" rid="fig2">Figure 2A</xref>). Among the induced genes identified by the two experimental platforms, only 102 are common. 291 genes are called as induced by the microarray experiment only. This group would include genes whose transcription may be stimulated at later time points via indirect mechanisms, but may also include true direct p53 target genes that require higher levels of p53 to be activated. For example, we noted that the canonical p53 target gene <italic>GADD45A</italic> fell in this group<italic>,</italic> as its transcription was mildly induced at 1 hr and thus fell below our statistical cut-off. Interestingly, 72 genes were identified as induced by GRO-seq only, despite the fact that the microarrays utilized harbored multiple probes against these mRNAs. The possible explanations for this finding are discussed below.</p><p>Second, microarrays detect 324 genes repressed upon 12 hr of Nutlin treatment, none of which were called as repressed by GRO-seq. The mechanism of p53-mediated gene repression remains debated in the field. Multiple independent ChIP-seq studies concur in that p53 binds weakly and very distally to those gene loci whose mRNAs are downregulated at the steady state level, and that the p53REs found at these sites match poorly to the consensus DNA sequence (<xref ref-type="bibr" rid="bib29">Nikulenkov et al., 2012</xref>; <xref ref-type="bibr" rid="bib26">Menendez et al., 2013</xref>; <xref ref-type="bibr" rid="bib37">Schlereth et al., 2013</xref>; <xref ref-type="bibr" rid="bib46">Wang et al., 2013</xref>). Using seven different available global ChIP datasets derived from HCT116 and two other cell lines, we created a collection of high confidence p53 binding events to analyze p53 binding in the vicinity of the various gene groups (‘Materials and methods’). Nearly 40% of the 198 genes induced by GRO-seq harbor a p53 binding event within 25 kb, significantly more than expected from random occurrence (p<italic>=</italic>1e−48, Hypergeometric test) (<xref ref-type="fig" rid="fig2">Figure 2B</xref>). Among the genes induced by microarray only, nearly 15% harbored p53 binding within 25 kb, still significantly more than expected by chance (p=8e−11), which suggests that some of these genes may be true direct targets activated at later time points. Most importantly, genes considered as repressed by the microarray profiling show little p53 binding within 25 kb, barely above what is expected by chance (p=3e−2), suggesting that the repression observed is largely indirect.</p><p>Indirect gene repression downstream of p53 activation could be mediated at the post-transcriptional level by p53-inducible miRNAs, and/or at the transcriptional level by the action of direct p53 targets known to repress transcription. Of note, GRO-seq identified 5 miRNAs directly transactivated by p53 (miR-1204, miR-3189, miR-34a, miR4679-1 and miR-4692, see <xref ref-type="supplementary-material" rid="SD1-data">Supplementary file 1</xref>). Most prominent among these is miR-34a, a well characterized p53-inducible miRNA known to mediate indirect repression by p53 at late time points. In fact, we found that nearly 72% of genes repressed in our microarray by Nutlin were previously shown by others (<xref ref-type="bibr" rid="bib21">Lal et al., 2011</xref>) to be downregulated upon overexpression of miR-34a in HCT116 cells (p<2.2e−16, Hypergeometric test, <xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1C</xref>). A recent report demonstrated that p21 and E2F4, a transcriptional repressor of S-phase genes acting coordinately with co-repressors of the RB family, are required for the downregulation of many genes previously characterized as ‘direct’ targets of p53 repression (<xref ref-type="bibr" rid="bib2">Benson et al., 2013</xref>). In agreement with these published findings, Ingenuity Pathway Analysis (IPA) of the genes repressed in our microarray experiment revealed that the top three regulators affecting these genes are indeed E2F4 (p=1.02e−81, Fisher's Exact Test), CDKN1A (p=8.21e−62) and RB (p=8.12e−60) (<xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1D</xref>). Altogether, these data indicate that most gene repression observed in our system is likely to be indirect, either via miRNAs, such as miR-34a, and/or the p21>RB>E2F4 axis (<xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1E</xref>). Of note, GRO-seq also identified the transcriptional repressor PRDM1 (BLIMP1) as a direct target of p53 (<xref ref-type="bibr" rid="bib49">Yan et al., 2007</xref>), revealing yet another possible mechanism for indirect gene repression downstream of p53. However, IPA did not identify PRDM1 as a top regulator of genes repressed in the microarray study (not shown).</p><p>Given that 72 genes were identified as activated only by GRO-seq, we further investigated the possible reasons for this result. Analysis of absolute signals generated by the GRO-seq (fpkm) vs microarray (mean fluorescence intensity, MFI) experiments generated several important insights (<xref ref-type="fig" rid="fig2">Figure 2C</xref>). First of all, robust p53 target genes such as <italic>CDKN1A</italic> and <italic>GDF15</italic> show strong increases in both platforms. When these results are expressed as fold induction, a strong correlation between transcriptional output and steady state levels is evident for these genes (<xref ref-type="fig" rid="fig2">Figure 2D</xref>). However, transcriptional output often does not correlate with steady state mRNA levels. For example, while the <italic>BTG2</italic> locus has a greater transcriptional output upon p53 activation than the <italic>MDM2</italic>, <italic>PTP4A1</italic>, <italic>AMZ2</italic> and <italic>CCNG1</italic> loci, its steady state mRNA levels are much lower, likely due to post-transcriptional repression of <italic>BTG2</italic> (<xref ref-type="fig" rid="fig2">Figure 2C</xref>). This would explain why a small group of ‘GRO-seq only' genes, including known p53 targets such as <italic>PTP4A1</italic> and <italic>CCNG1,</italic> are not called by the microarray experiments: they display very high basal steady state levels, which do not increase significantly despite clear transcriptional induction (<xref ref-type="fig" rid="fig2">Figure 2C</xref>). However, most ‘GRO-seq only’ genes belong to a different category, as marked by their very low levels of transcription and steady state mRNA expression. Genes like <italic>CDC42BPG</italic>, <italic>ADAMTS7</italic> and <italic>LRP1</italic> are clearly induced at the transcriptional level but show no apparent increase in the microarray signals (<xref ref-type="fig" rid="fig2">Figure 2C,D</xref>). Remarkably, the induction of these mRNAs upon p53 activation is evident by Q-RT-PCR at the time point of the microarray experiment (<xref ref-type="fig" rid="fig2">Figure 2E</xref>). In fact, when microarray-derived signals are displayed for the various gene sets, the ‘GRO-seq only’ group shows very significant lower expression as compared to all other groups (<xref ref-type="fig" rid="fig2">Figure 2F</xref>). Altogether, we conclude that microarray profiling is not sensitive enough to detect these low abundance transcripts, which could explain why several published ChIP-seq/microarray studies failed to identify these genes as direct p53 targets. Alternatively, it is possible that p53 binds to these genes from very distal sites outside of the arbitrary window defined during bioinformatics analysis of ChIP-seq data. To discern among these possibilities, we analyzed ChIP-seq data in search of high confidence p53 binding events in the vicinity of several novel genes identified by GRO-seq, and evaluated p53 binding using standard ChIP assays. Indeed, we detected clear p53 binding to all p53REs tested at these novel p53 targets (<xref ref-type="fig" rid="fig2s2">Figure 2—figure supplement 2</xref>). Of note, p53 binds to proximal regions at the <italic>CDC42BPG</italic> and <italic>LRP1</italic> loci (+1373 bp and −694 bp relative to transcription start site [TSS], respectively), indicating that these genes could have been missed in previous studies due to the low abundance of their transcripts. In contrast, p53 binds to very distal sites (i.e., >30 kb from the TSS) at the <italic>ADAMTS7</italic>, <italic>TOB1</italic>, <italic>ASS1</italic> and <italic>CEP85L</italic> loci (<xref ref-type="fig" rid="fig2s2">Figure 2—figure supplement 2</xref>), suggesting that these genes would have been missed as direct targets when setting an arbitrary <30 kb window during ChIP-seq analysis. In summary, GRO-seq enables the identification of novel direct p53 target genes due both to its increased sensitivity and the fact that it does not require proximal p53 binding to ascertain direct regulation.</p></sec><sec id="s2-3"><title>p53 represses a subset of its direct target genes prior to MDM2 inhibition</title><p>Others and we have observed that in proliferating cells with minimal p53 activity, p53 increases the basal expression of some of its target genes (<xref ref-type="bibr" rid="bib41">Tang et al., 1998</xref>; <xref ref-type="bibr" rid="bib11">Espinosa et al., 2003</xref>). This was first recorded for <italic>CDKN1A</italic> (<xref ref-type="bibr" rid="bib41">Tang et al., 1998</xref>), and it is confirmed by our GRO-seq analysis (<xref ref-type="fig" rid="fig1">Figure 1A</xref>, compare 2.6 to 5.7 fpkm in the Control tracks). To investigate whether this is a general phenomenon we analyzed the basal transcription of all p53-activated genes in control p53 +/+ vs p53 −/− cells (<xref ref-type="fig" rid="fig3">Figure 3A,B</xref>). Interestingly, p53 status exerts differential effects among its target genes prior to MDM2 inhibition with Nutlin. While many genes show the same behavior as <italic>CDKN1A</italic> (e.g., <italic>GDF15</italic>, <italic>DDB2</italic>, labeled green throughout <xref ref-type="fig" rid="fig3">Figure 3</xref>), another group shows decreased transcription in the presence of MDM2-bound p53 (e.g., <italic>PTP4A1</italic>, <italic>HES2</italic>, <italic>GJB5,</italic> labeled red throughout <xref ref-type="fig" rid="fig3">Figure 3</xref>). Genome browser views illustrating this phenomena are provided for <italic>GDF15</italic> and <italic>PTP4A1</italic> in <xref ref-type="fig" rid="fig3">Figure 3C</xref>. The differential behavior of RNAPII at these gene loci is also observed in ChIP assays using antibodies against the Serine 5- and Serine 2-phosphorylated forms of the RBP1 C-terminal domain repeats, which mark initiating and elongating RNAPII complexes, respectively (S5P- and S2P-RNAPII, <xref ref-type="fig" rid="fig3s1">Figure 3—figure supplement 1A</xref>). Whereas the ‘basally activated’ <italic>GDF15</italic> locus displays higher GRO-seq and RNAPII ChIP signals in untreated p53 +/+ cells, the ‘basally repressed’ <italic>PTP4A1</italic> locus shows lower signals in the presence of MDM2-p53 complexes. Such differential effects among p53 target genes have a clear impact on their absolute level of transcription upon MDM2 inhibition: whereas those whose basal activity was increased ranked among the most differentially transcribed between Nutlin-treated p53 +/+ and p53 −/− cells (top of the heatmap in <xref ref-type="fig" rid="fig3">Figure 3B</xref>); those basally repressed are virtually no different in expression between Nutlin-treated p53 +/+ and p53 −/− cells (bottom of the heatmap in <xref ref-type="fig" rid="fig3">Figure 3B</xref>). Importantly, Q-RT-PCR shows that the differential effects of p53 on the basal transcription of its targets are generally translated into differences in mRNA steady level (<xref ref-type="fig" rid="fig3">Figure 3D</xref>).</p><p>Overall, these results indicate that p53 acts as a repressor at a subset of its targets in a manner that is relieved by Nutlin, suggesting that MDM2 recruitment by basal levels of p53 may repress transcription at specific loci. To test this hypothesis, we performed ChIP experiments for p53 and MDM2 under conditions matching the GRO-seq experiment. For the p53 ChIP, we employed the monoclonal antibody DO-1 that recognizes the p53 TAD1 and whose reactivity should increase upon displacement of MDM2 by Nutlin. Importantly, ChIP assays show a significant amount of chromatin-bound p53 above background levels at the p53REs in the <italic>CDKN1A</italic> locus (basally activated), and the <italic>PTP4A1</italic> and <italic>HES2</italic> loci (basally repressed), even before Nutlin treatment (<xref ref-type="fig" rid="fig3">Figure 3E</xref>, <xref ref-type="fig" rid="fig3s1">Figure 3—figure supplement 1B</xref>). Of note, the DO-1 ChIP signals increase upon Nutlin treatment, as expected from epitope unmasking. ChIP assays also detect MDM2 chromatin binding above background levels at these three p53REs, with signals decreasing upon Nutlin treatment, as expected by the competitive action of this molecule (<xref ref-type="fig" rid="fig3">Figure 3E</xref>, <xref ref-type="fig" rid="fig3s1">Figure 3—figure supplement 1B</xref>). Of note, although Nutlin disrupts the interaction between the p53 N-terminus and the hydrophobic pocket in the N-terminal domain of MDM2, a second molecular interaction occurs between the p53 C-terminus and the MDM2 N-terminus that is not competed by Nutlin in vitro (<xref ref-type="bibr" rid="bib34">Poyurovsky et al., 2010</xref>), which may explain why the MDM2 signal is not completely abrogated upon a short time point of Nutlin treatment. Western blots demonstrating specific MDM2 immunoprecipitation under the ChIP conditions utilized are shown in <xref ref-type="fig" rid="fig3s1">Figure 3—figure supplement 1C</xref>.</p><p>Next, we investigated whether the differential effects of basal p53 levels on expression of its direct targets could be revealed in an analysis of hundreds of cell lines expressing wild type vs mutant p53. More specifically, we hypothesized that genes that are basally transactivated by p53 would be more highly expressed in p53 WT cells than transrepressed genes. Oncomine analysis of 598 cancer cell lines not treated with p53-activating agents revealed that many genes that are ‘basally activated’ in HCT116 cells such as <italic>CDKN1A</italic>, <italic>DDB2</italic> and <italic>GDF15</italic> indeed show significantly higher mRNA expression in WT p53 cell lines (<xref ref-type="fig" rid="fig3">Figure 3F</xref>). In contrast, the ‘basally repressed’ gene <italic>GJB5</italic> shows significantly higher expression in mutant p53 cell lines. When all 12,624 genes in the Oncomine analysis are ranked according to their relative expression in WT p53 over mutant p53 cell lines, many genes whose basal transcription is upregulated by p53 in HCT116 cells appear at the top of this ranking (e.g., <italic>CDKN1A</italic>, <italic>DDB2</italic> and <italic>GDF15</italic>, ranked 2, 4 and 62, respectively) (<xref ref-type="fig" rid="fig3s2">Figure 3—figure supplement 2A</xref>). However, some direct targets ‘basally repressed’ by p53, such as <italic>GJB5</italic>, show an inverse correlation with WT p53 status. Collectivelly, the direct p53 targets identified by GRO-seq are enriched toward the top of the ranking, which is revealed in a Gene set enrichment analysis (GSEA) (<xref ref-type="fig" rid="fig3s2">Figure 3—figure supplement 2A</xref>). In contrast, genes induced only in the microarray platform (i.e., mostly indirect targets) do not show strong enrichment in a GSEA analysis. When the relative basal transcription between HCT116 p53 +/+ and p53 −/− cells is plotted against the relative mRNA expression in p53 WT vs p53 mutant cell lines, it is apparent that many ‘basally activated’ genes are more highly expressed in p53 WT cells (green dots in the upper right quadrant in <xref ref-type="fig" rid="fig3s2">Figure 3—figure supplement 2B</xref>). Finally, the differential pattern of basal expression among p53 targets is also observed in an analysis of 256 breast tumors for which p53 status was determined, where <italic>CDKN1A</italic>, <italic>DDB2</italic> and <italic>GDF15</italic> (but not <italic>GJB5</italic>) show higher expression in the p53 WT tumors (<xref ref-type="fig" rid="fig3s2">Figure 3—figure supplement 2C</xref>).</p><p>Altogether, these results reveal a qualitative difference among p53 target genes in terms of their sensitivity to basal p53-MDM2 complexes as depicted in <xref ref-type="fig" rid="fig2s2">Figure 3—figure supplement 2D</xref>. Although indirect effects can not be fully ruled out, the fact that we can detect p53 and MDM2 binding to the p53REs near these gene loci suggest direct action. Of note, early in vitro transcription studies demonstrated that MDM2 represses transcription when tethered to DNA independently of p53, which may provide the molecular mechanism behind our observations (<xref ref-type="bibr" rid="bib42">Thut et al., 1997</xref>) (‘Discussion’).</p></sec><sec id="s2-4"><title>GRO-seq reveals gene-specific regulatory mechanisms affecting key survival and apoptotic genes</title><p>Many research efforts have been devoted to the characterization of molecular mechanisms conferring gene-specific regulation within the p53 network, as these mechanisms could be exploited to manipulate the cellular response to p53 activation. Most research has focused on factors that differentially modulate p53 binding or transactivation of survival vs apoptotic genes (<xref ref-type="bibr" rid="bib45">Vousden and Prives, 2009</xref>). GRO-seq identified a plethora of gene-specific regulatory features affecting p53 targets, but our analysis failed to reveal a universal discriminator between survival and death genes within the network.</p><p>When direct p53 target genes with well-established pro-survival (i.e., cell cycle arrest, survival, DNA repair and negative regulation of p53) and pro-death (i.e., extrinsic and intrinsic apoptotic signaling) functions are ranked based on their transcriptional output in Nutlin-treated p53 +/+ cells, it is evident that key pro-survival genes such as <italic>CDKN1A</italic>, <italic>GDF15</italic>, <italic>BTG2</italic> and <italic>MDM2</italic> are transcribed at much higher rates than any apoptotic gene (<xref ref-type="fig" rid="fig4">Figure 4A</xref>). For example, ∼20-fold more RNA is produced from the <italic>CDKN1A</italic> locus than from the <italic>BBC3</italic> locus encoding the BH3-only protein PUMA. The most potently transcribed apoptotic gene is <italic>TP53I3</italic> (<italic>PIG3</italic>), yet its transcriptional output is still 3.4-fold lower than <italic>CDKN1A</italic>. Based on measurements of steady state RNA levels, it was observed that apoptotic genes such as <italic>TP53I3</italic> and <italic>FAS</italic> are induced with a slower kinetics than <italic>CDKN1A</italic> (<xref ref-type="bibr" rid="bib40">Szak et al., 2001</xref>; <xref ref-type="bibr" rid="bib11">Espinosa et al., 2003</xref>). However, GRO-seq suggests that this is not due to a slower kinetics of RNAPII transactivation, but rather to a lower transcriptional output from the apoptotic loci. Although both <italic>CDKN1A</italic> and <italic>TP53I3</italic> display similar transcriptional induction within 1 hr of Nutlin treatment (7.4 and 6.07 fold induction, <xref ref-type="supplementary-material" rid="SD1-data">Supplementary file 1</xref>), the <italic>TP53I3</italic> mRNA takes longer to display increased accumulation over basal levels (see Q-RT-PCR in <xref ref-type="fig" rid="fig1">Figure 1B</xref>) (<xref ref-type="bibr" rid="bib40">Szak et al., 2001</xref>). Thus, differences in mRNA induction kinetics between gene classes could be explained from differential amounts of RNA synthesis. To investigate how much of the differential transcriptional output is due to p53 action vs p53-autonomous mechanisms, we analyzed the activity of these gene loci in p53 −/− cells (<xref ref-type="fig" rid="fig4">Figure 4B</xref>). Several important observations arise from this analysis. First, the basal activity of the <italic>CDKN1A</italic> locus is higher than most pro-apoptotic genes even in p53 −/− cells, which agrees with previous studies revealing the action of strong core promoter elements at this locus (<xref ref-type="bibr" rid="bib10">Espinosa and Emerson, 2001</xref>; <xref ref-type="bibr" rid="bib28">Morachis et al., 2010</xref>). However, the pro-apoptotic genes <italic>TRAF4</italic> and <italic>AEN</italic> are exceptions to this trend, as they display higher basal activity in p53 −/− than <italic>CDKN1A</italic> and most pro-survival genes. Second, p53 action reinforces the distinction between the two classes by leading to ‘super-activation’ (i.e., greater than fivefold) of the survival genes <italic>CDKN1A</italic>, <italic>GDF15</italic>, <italic>BTG2</italic> and <italic>MDM2</italic> (<xref ref-type="fig" rid="fig4">Figure 4C</xref>). Although select apoptotic genes such as <italic>TP53I3</italic>, <italic>PHLDA3</italic> and <italic>FAS</italic> also undergo super-activation, this does not suffice to override their overall lower transcriptional output. In sum, as a group, survival genes tend to be transcribed at a higher extent than apoptotic genes, which is due to a combination of p53-dependent and -independent mechanisms.</p><p>Next, we investigated whether RNAPII pausing exerted gene-specific effects within the p53 transcriptional program. Using ChIP assays, we previously reported that the <italic>CDKN1A</italic> promoter carries significantly higher levels of promoter-bound RNAPII than the apoptotic genes <italic>TP53I3</italic>, <italic>TNFRSF10B</italic>, <italic>BBC3</italic> and <italic>FAS</italic> prior to p53 activation (<xref ref-type="bibr" rid="bib11">Espinosa et al., 2003</xref>). GRO-seq confirms that the promoters of survival genes, including <italic>CDKN1A</italic>, indeed carry more transcriptionally engaged RNAPII than the promoters of apoptotic genes (<xref ref-type="fig" rid="fig4">Figure 4D</xref>); however, there was no obvious correlation between the amount of active RNAPII over the promoter and the degree of transcriptional output or induction. RNAPII pausing was proposed to modulate the timing of signal-induced gene expression, such as during the cellular response to LPS, where primary response genes were found to carry more paused RNAPII than secondary response genes (<xref ref-type="bibr" rid="bib15">Hargreaves et al., 2009</xref>). However, others found that RNAPII pausing is not necessary for rapid gene induction (<xref ref-type="bibr" rid="bib14">Hah et al., 2011</xref>). To investigate this issue more thoroughly within the p53 network, we performed an analysis of pausing indices as previously described (<xref ref-type="bibr" rid="bib7">Core et al., 2008</xref>). This exercise revealed that although p53 target genes vary widely in their pausing indexes, there is no obvious correlation between pausing and the degree of transcriptional induction at the 1 hr time point of Nutlin treatment (<xref ref-type="fig" rid="fig4s1">Figure 4—figure supplement 1A</xref>). For example, the <italic>POLH</italic> gene, which shows a very high pausing index, displays a lower fold induction and overall lower transcriptional output than the <italic>PHLDA3</italic> gene, which shows little signs of RNAPII pausing (<xref ref-type="fig" rid="fig4s1">Figure 4—figure supplement 1A</xref>). Thus, RNAPII pausing is not a pre-requisite for rapid induction within the p53 transcriptional program. The first GRO-seq analysis in human cells revealed widespread divergent transcription at most active promoters (<xref ref-type="bibr" rid="bib7">Core et al., 2008</xref>). We observed that the degree of divergence is highly variable across promoters of p53 target genes. For example, the <italic>DRAM1</italic> promoter displays more divergent transcription than sense transcription (<xref ref-type="fig" rid="fig4s1">Figure 4—figure supplement 1B</xref>), yet divergent transcription is minor at gene loci such as <italic>CDKN1A</italic> or <italic>DDB2</italic> (<xref ref-type="fig" rid="fig1">Figure 1A</xref>, <xref ref-type="fig" rid="fig4s1">Figure 4—figure supplement 1B</xref>, respectively). Using a ‘divergence index’, which evaluates the ratio of productive (sense) vs divergent transcription, we ranked p53 target genes from high to low promoter divergence, which revealed that there is no correlation between divergence and fold induction (<xref ref-type="fig" rid="fig4s1">Figure 4—figure supplement 1B</xref>). Thus, divergence of RNAPII in the unproductive direction does not obviously affect p53 transactivation of its target genes.</p><p>Overall, we did not find evidence of an overarching regulatory feature discriminating between survival vs apoptotic genes. Instead, GRO-seq analysis revealed many instances of potential for gene-specific regulation at key loci. For example, the <italic>APAF1</italic> gene<italic>,</italic> encoding a core component of the apoptosome, is activated from a non-overlapping bidirectional promoter with the <italic>IKBIP</italic> gene, where the two start sites are less than 280 bp apart (<xref ref-type="fig" rid="fig4s2">Figure 4—figure supplement 2A</xref>). The <italic>FAS</italic> gene, which encodes a death receptor involved in p53-dependent apoptosis<italic>,</italic> is transcribed from a bidirectional overlapping promoter with the <italic>ACTA2</italic> gene, where the start site of one gene resides within the first exon of the other gene, with the two transcription units overlapping for ∼850 bp (<xref ref-type="fig" rid="fig4s2">Figure 4—figure supplement 2B</xref>). Several p53 target genes are organized in clusters, where a group of genes within a chromatin domain are concurrently upregulated. One example is the <italic>TNFRS10</italic> cluster encoding the pro-apoptotic TRAIL receptors DR4 (<italic>TNFRSF10A)</italic> and DR5/Killer (<italic>TNFRSF10B</italic>), as well as the counteracting decoy receptors DCR1 (<italic>TNFRSF10C)</italic> and DCR2 (<italic>TNFRSF10D</italic>) (<xref ref-type="fig" rid="fig4s2">Figure 4—figure supplement 2C</xref>). Another interesting feature revealed by GRO-seq is the presence of significant intragenic antisense transcription at many p53 target genes. In many cases such antisense transcription can be attributed to eRNAs produced from intronic p53 enhancers (see later, <xref ref-type="fig" rid="fig5">Figure 5</xref>); however, antisense transcription often originates from an overlapping transcriptional unit that is neither derived from a p53 binding event nor induced by p53. One example of this scenario is the pro-survival p53 target gene <italic>GPR87</italic>, which is embedded within an intron of <italic>MED12L</italic>, which is not a p53 target (<xref ref-type="fig" rid="fig4s2">Figure 4—figure supplement 2D</xref>). Thus, our GRO-seq dataset paves the way for many future mechanistic studies aimed at deciphering the impact of these gene-specific events in the regulation of p53 target genes and the cellular outcome upon p53 activation.<fig-group><fig id="fig5" position="float"><object-id pub-id-type="doi">10.7554/eLife.02200.014</object-id><label>Figure 5.</label><caption><title>Direct p53 target genes harbor pre-activated enhancers.</title><p>(<bold>A</bold>) GRO-seq results for the <italic>DDIT4</italic> locus, representative of p53 target genes that display bidirectional eRNA transcription (arrow) arising near sites of p53 binding (indicated by a purple dot). (<bold>B</bold>) Analysis of nearest p53 binding events relative to the transcription start site (TSS) of direct p53 target genes detected by GRO-seq (grey bars) vs all RefSeq genes (pink). (<bold>C</bold>) GRO-seq results for the <italic>SYTL1</italic> locus, representative of p53 target genes that display intronic antisense eRNA transcription arising near sites of p53 binding. (<bold>D</bold>) Analysis of eRNA transcription at distal p53 binding sites (>25 kb of any gene), proximal sites associate with a gene not activated by p53 (<25 kb of non-target), and those proximal to a p53 target gene identified by GRO-seq. (<bold>E</bold>) p53 binding sites near target genes have higher transcription levels than sites near other genes even in p53 null cells, indicating the likely action of pioneer factors.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.02200.014">http://dx.doi.org/10.7554/eLife.02200.014</ext-link></p></caption><graphic xlink:href="elife02200f005"/></fig><fig id="fig5s1" position="float" specific-use="child-fig"><object-id pub-id-type="doi">10.7554/eLife.02200.015</object-id><label>Figure 5—figure supplement 1.</label><caption><title>p53 stimulates eRNA production at extragenic and intragenic locations.</title><p>(<bold>A</bold> and <bold>B)</bold> GRO-seq results for the region upstream of <italic>CDKN1A</italic> showing eRNAs derived from the p53REs at −1.3 and −2.4 kb (<bold>A</bold>) and from the region encoding <italic>lincRNA-p21</italic> (<bold>B</bold>). (<bold>C</bold>) GRO-seq results for the <italic>DRAM1</italic> locus, a rare example of a p53 target gene whose enhancer is not obviously transcribed in p53 null cells. (<bold>D</bold>) GRO-seq results for the <italic>BTG2</italic> locus, representative of p53 target genes that display intronic antisense eRNA transcription arising near sites of p53 binding. (<bold>E</bold>) Percentage of p53 binding sites that are transcribed as defined by <xref ref-type="bibr" rid="bib7">Core et al. (2008)</xref>. The distal sites are transcribed less often than the sites proximal to target genes. A higher percentage of sites near p53 target genes are transcribed than sites that are distal or near non-targets. (<bold>F</bold>) The fold change in transcription after addition of Nutlin is similar across p53 binding sites in various locations. Expectedly, in p53 −/− cells, there is no change in transcription after Nutlin treatment. In p53 +/+ cells eRNAs are increased in transcription is approximately fourfold and this is true for distal sites, sites proximal to non-targets and sites proximal to targets.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.02200.015">http://dx.doi.org/10.7554/eLife.02200.015</ext-link></p></caption><graphic xlink:href="elife02200fs008"/></fig></fig-group></p></sec><sec id="s2-5"><title>Direct p53 target genes harbor highly transcribed p53 binding sites</title><p>Recently, a novel class of non-coding RNAs has been identified whose transcription originates from active enhancers. These enhancer-derived RNAs (eRNAs) have been shown to contribute to gene activation in a variety of systems (<xref ref-type="bibr" rid="bib18">Jiao and Slack, 2013</xref>). A recent report characterized eRNAs derived from three distal p53 enhancers and showed that they are required for efficient p53 transactivation of neighboring genes (<xref ref-type="bibr" rid="bib25">Melo et al., 2013</xref>). In order to investigate the prevalence of transcriptionally active enhancers within the p53 transcriptional program, we examined our GRO-seq data with respect to hundreds of p53 binding events as defined by ChIP-seq. Of note, we have not employed here data on histone marks or p300 occupancy to define how many of these p53 binding events reside within regions harboring the accepted hallmarks of enhancers, and thus some of these p53 binding sites should be considered as putative enhancers.</p><p>GRO-seq readily detects RNAs originating from most p53 binding events, which we refer hereto as eRNAs. A typical example is shown for the <italic>DDIT4</italic> locus in <xref ref-type="fig" rid="fig5">Figure 5A</xref>, where a distal p53 binding site located downstream of the gene is clearly transcribed in both the sense and antisense directions, with increased signals upon p53 activation. Interestingly, this p53RE is also transcribed in p53 −/− cells (<xref ref-type="fig" rid="fig5">Figure 5A</xref>, top track, arrow). Analysis of the <italic>CDKN1A</italic> locus shows transcription from the well characterized p53REs at −1.3 and −2.4 kb (<xref ref-type="fig" rid="fig5s1">Figure 5—figure supplement 1A</xref>). Analysis of the distal upstream region in this locus encoding the long intragenic ncRNA known as lincRNA-p21 shows transcription in both strands originating from a p53 binding site, with the antisense strand corresponding to the reported lncRNA-p21 sequence (<xref ref-type="fig" rid="fig5s1">Figure 5—figure supplement 1B</xref>). This suggests that lncRNA-p21 could be classified as an eRNA, as it originates from the vicinity of a p53RE associated to a canonical p53 target gene. Once again, transcripts derived from the lincRNA-p21 region can also be detected in p53 −/− cells (<xref ref-type="fig" rid="fig5s1">Figure 5—figure supplement 1B</xref>, top track). A rare example of a p53RE near a target gene not transcribed in p53 −/− cells is that of the <italic>DRAM1</italic> locus, which displays transcription of bidirectional eRNAs in p53 +/+ cells before p53 activation, with signals increasing upon Nutlin treatment (<xref ref-type="fig" rid="fig5s1">Figure 5—figure supplement 1C</xref>).</p><p>Analysis of the spatial distribution of p53 binding events relative to transcription start sites (TSSs) shows that direct p53 target genes display an enrichment in p53 binding close to promoters, but also within genes (<xref ref-type="fig" rid="fig5">Figure 5B</xref>). In fact, it has been estimated that ∼40% of p53 enhancers are intragenic (<xref ref-type="bibr" rid="bib29">Nikulenkov et al., 2012</xref>; <xref ref-type="bibr" rid="bib26">Menendez et al., 2013</xref>; <xref ref-type="bibr" rid="bib37">Schlereth et al., 2013</xref>; <xref ref-type="bibr" rid="bib46">Wang et al., 2013</xref>). Although eRNAs derived from the sense strands can not be distinguished from the protein coding pre-mRNAs at these locations, the eRNAs arising from the antisense strands are clearly distinguishable, as illustrated for the <italic>SYTL</italic> and <italic>BTG2</italic> loci (<xref ref-type="fig" rid="fig5">Figure 5C</xref>, <xref ref-type="fig" rid="fig5s1">Figure 5—figure supplement 1D</xref>, respectively). Thus, p53 activation leads to antisense transcription within a large fraction of its direct target genes concurrently with activation of the protein-coding RNAs, a phenomenon with potential regulatory consequences.</p><p>Next, we analyzed the production of eRNAs at three different sets of p53 binding events: (a) distal binding sites (>25 kb of any gene), (b) proximal binding sites associated with a gene not activated by p53 (<25 kb of non GRO-seq target gene), and (c) proximal binding sites associated with a p53 target gene (<25 kb of a p53 target defined by GRO-seq). First, we asked what fraction of sites in each category is transcriptionally active in the four experimental conditions. This revealed that ∼40% of distal sites and those not associated with p53 targets are transcriptionally active in p53 −/− cells (<xref ref-type="fig" rid="fig5s1">Figure 5—figure supplement 1E</xref>). The number increases to ∼60% in p53 +/+ cells, likely revealing the action of basal p53 levels, as observed for the <italic>DRAM1</italic> locus. Expectedly, Nutlin leads to further increase in the fraction of active enhancers only in p53 +/+ cells. When the analysis is restricted to those p53REs within 25 kb of a direct p53 target, the fraction that are transcribed in p53 −/− cells increases to ∼70%, and nearly all of them are transcriptionally active in Nutlin-treated p53 +/+ cells (<xref ref-type="fig" rid="fig5s1">Figure 5—figure supplement 1E</xref>). Thus, eRNAs are a hallmark of active p53 binding sites. We next investigated whether there are differences across p53 binding events found in the various locations in terms of eRNA activation and overall expression. A fold induction analysis shows that p53 stimulates eRNA transcription to a similar degree regardless of location (<xref ref-type="fig" rid="fig5s1">Figure 5—figure supplement 1F</xref>). This result counters the notion that p53 acts as transcriptional repressor from distal sites by a mechanism described as ‘enhancer interference’ (<xref ref-type="bibr" rid="bib23">Li et al., 2012</xref>). Interestingly, the absolute amount of eRNA produced varies greatly with location (<xref ref-type="fig" rid="fig5">Figure 5D</xref>). First, eRNA transcription is much weaker from distal sites relative to proximal sites (p=1.28e−12 for Nutlin-treated p53 +/+ cells). This could be explained by more efficient communication between enhancers and promoters when in closer proximity. Second, eRNA transcription is significantly higher from proximal sites associated with direct p53 target genes. Strikingly, this difference is already evident in p53 −/− cells (p=1e−10). Thus, a distinctive feature between direct p53 target genes and those genes proximal to bound p53REs but not activated is the strength of eRNA production in the absence of p53. In other words, p53 seems to activate gene expression at genomic locations carrying ‘primed’, pre-activated enhancers, likely revealing the action of pioneer factors (<xref ref-type="fig" rid="fig5">Figure 5E</xref>).</p></sec></sec><sec id="s3" sec-type="discussion"><title>Discussion</title><p>The importance of p53 in cancer biology is undisputed, yet the mechanisms by which this transcription factor suppresses tumor growth remain to be fully elucidated. In particular, it is unclear which p53 target genes contribute to tumor suppression in various contexts. A thorough analysis of the literature up to 2008 revealed ∼120 direct p53 target genes (<xref ref-type="bibr" rid="bib36">Riley et al., 2008</xref>). Since then, genomics experiments using microarrays and ChIP-seq suggest thousands of p53 targets, but very few genes were commonly identified by multiple studies (<xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1A,B</xref>) (<xref ref-type="bibr" rid="bib29">Nikulenkov et al., 2012</xref>; <xref ref-type="bibr" rid="bib26">Menendez et al., 2013</xref>; <xref ref-type="bibr" rid="bib37">Schlereth et al., 2013</xref>; <xref ref-type="bibr" rid="bib46">Wang et al., 2013</xref>). The lack of overlap between these reports could be due to the fact that they employed different cell types and somewhat different experimental platforms. However, measurements of RNA steady state levels may produce a misleading view of direct p53 action, as they confound direct vs indirect effects. Thus, it is possible that cell type-specific secondary effects and post-transcriptional regulatory mechanisms strongly qualify the direct transcriptional response in different cell types. Ascertaining direct p53 action by the mere presence of a p53 binding event within an arbitrary distance to a putative target gene is an imprecise exercise, as p53 may act directly from very distal sites outside of this arbitrary cutoff (leading to false negatives) and because many proximal p53 binding events may be unproductive (leading to false positives). Because of these caveats, we investigated direct transcriptional regulation by p53 using GRO-seq.</p><p>A surprising result from our GRO-seq analysis is that a short time point of MDM2 inhibition suffices for p53 to activate hundreds of genomic loci, even prior to any detectable increase in total p53 levels. Because MDM2 functions as an E3 ligase targeting p53 for degradation (<xref ref-type="bibr" rid="bib20">Kubbutat et al., 1997</xref>), there was no guarantee that the low basal levels of p53 present in a proliferating cell culture would suffice to induce transcription of its target genes. Importantly, ChIP assays demonstrate that p53 and MDM2 occupy p53REs in proliferating cells and that MDM2 binding is decreased upon Nutlin treatment (<xref ref-type="fig" rid="fig3">Figure 3E</xref>). These observations emphasize the role of MDM2 as a transcriptional repressor by masking of the p53 transactivation domains (<xref ref-type="bibr" rid="bib31">Oliner et al., 1993</xref>), but do not negate the importance of p53 degradation as a repressive mechanism, as it is possible that increased p53 levels are required for activation of target genes at later time points. Our results contrast the notion that apoptotic genes require higher levels of p53 for transactivation or that they are transcriptionally induced at later time points, highlighting instead the ‘primed’ nature of a multifunctional p53 transcriptional response. Furthermore, this confirms that the failure of many cell types to undergo apoptosis upon Nutlin treatment is not due to a defect in transactivation of key apoptotic genes (<xref ref-type="bibr" rid="bib16">Henry et al., 2012</xref>; <xref ref-type="bibr" rid="bib39">Sullivan et al., 2012</xref>).</p><p>Although p53 action leads to massive gene repression at a global scale over time, it is unclear how much of these effects are direct vs indirect. Previous genomics experiments identified dozens of genes that are both bound by p53 within a certain arbitrary distance and whose steady state RNA levels decrease at late time points (<xref ref-type="bibr" rid="bib29">Nikulenkov et al., 2012</xref>; <xref ref-type="bibr" rid="bib26">Menendez et al., 2013</xref>; <xref ref-type="bibr" rid="bib37">Schlereth et al., 2013</xref>; <xref ref-type="bibr" rid="bib46">Wang et al., 2013</xref>). However, meta-analysis of these reports failed to identify a single gene commonly repressed in more than one study (<xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1A,B</xref>). Recent work showed that p21 is both necessary and sufficient to downregulate many genes commonly described as direct targets of p53 repression, mostly acting via E2F4 (<xref ref-type="bibr" rid="bib2">Benson et al., 2013</xref>). Other cell cycle inhibitory pathways may also converge on E2F4 repressive complexes, such as the p53-inducible miRNA miR-34a, which targets the mRNAs encoding G1-S cyclins (<xref ref-type="bibr" rid="bib21">Lal et al., 2011</xref>). Our data supports the notion that most repression downstream of p53 activation is indirect. First, MDM2 inhibition by 1 hr Nutlin treatment identified only four repressed genes, none of which showed repression at the steady state levels. In contrast, a microarray experiment at 12 hr showed hundreds of downregulated genes. Analysis of this gene set strongly supports the notion that E2F4, p21, RB and miR-34a largely mediate their repression (<xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1C–E</xref>).</p><p>Interestingly, GRO-seq analysis of p53 null cells revealed that p53-MDM2 complexes might directly repress transcription at a subset of p53 targets. These genes are downregulated in the presence of MDM2-bound p53 but then activated by Nutlin. These results reveal that basal amounts of p53 found in proliferating cells create an uneven landscape among its transactivation targets, pre-activating some and repressing others. Mechanistically, p53-MDM2 complexes may directly repress transcription due to the inhibitory effects of MDM2 on components of the Pre-Initiation Complex (PIC). Early work by Tjian et al. using in vitro transcription assays demonstrated a dual mechanism of transcription inhibition by MDM2 (<xref ref-type="bibr" rid="bib42">Thut et al., 1997</xref>). Their biochemical assays demonstrated that MDM2 not only masks the p53 transactivation domain, but that it also represses transcription when tethered to DNA by a GAL4 DNA binding domain. They identified an inhibitory domain in MDM2 that binds to the PIC components TBP and TFIIE, and hypothesized that MDM2 could repress transcription by targeting the basal transcription machinery. Our GRO-seq results identify specific p53 targets where this mechanism could be taking place and ChIP experiments using p53 and MDM2 antibodies confirm binding of both proteins to the p53REs at these loci. In agreement with these results, others have previously demonstrated that in proliferating cells MDM2 binds to p53REs in a p53-dependent manner, and that MDM2 recruitment to chromatin can be disrupted by Nutlin or DNA damaging agents (<xref ref-type="bibr" rid="bib48">White et al., 2006</xref>). Also, excess MDM2 was shown to exert uneven repressive effects on the expression of p53 target genes, independently of effects on p53 levels or chromatin binding (<xref ref-type="bibr" rid="bib30">Ohkubo et al., 2006</xref>). Altogether, these data support the arising notion that MDM2 works as a gene-specific co-regulator of p53 target genes by mechanisms other than mere p53 inhibition (<xref ref-type="bibr" rid="bib3">Biderman et al., 2012</xref>).</p><p>Many research efforts in the p53 field have been devoted to the characterization of regulatory mechanisms discriminating between survival and apoptotic genes. Our GRO-seq analysis reinforced the notion that <italic>CDKN1A,</italic> a key mediator of arrest, differs from key apoptotic genes in several aspects. <italic>CDKN1A</italic> has outstanding transcriptional output among p53 target genes, which is partly due to the fact that its promoter drives substantial p53-independent transcription, but also due to potent p53-dependent transactivation. In vitro transcription assays demonstrated the <italic>CDKN1A</italic> core promoter initiates transcription more rapidly and effectively than the <italic>FAS</italic> core promoter (<xref ref-type="bibr" rid="bib28">Morachis et al., 2010</xref>), and GRO-seq confirms that <italic>FAS</italic> has weaker transcriptional output than <italic>CDKN1A</italic>. However, our GRO-seq analysis failed to identify a uniform criterion discriminating between the most well studied survival and apoptotic genes. To the contrary, GRO-seq revealed that each individual p53 target gene is subject to various layers of gene-specific regulatory mechanisms, including but not restricted to differential levels of p53-independent transcription, p53 transactivation potential, RNAPII pausing, promoter divergence, extragenic vs intragenic eRNAs, overlapping promoters, clustered activation and antisense transcription.</p><p>A key observation arising from our GRO-seq analysis is that p53 target genes often have ‘primed’ p53REs, as denoted by significantly higher levels of eRNA production in p53 null cells. We interpret this result as the action of unknown pioneering factors acting at these putative enhancers prior to p53 signaling, which would establish enhancer-promoter communication and ready these genes for further transactivation by p53 or other stimulus-induced transcription factors. This notion is supported by a recent analysis of eRNAs at three distal p53 binding sites, which were shown to be involved in long range chromatin loops independently of p53 (<xref ref-type="bibr" rid="bib25">Melo et al., 2013</xref>). This model also agrees with a recent report showing that TNF-responsive enhancers are in physical contact with their target promoters prior to TNF signaling (<xref ref-type="bibr" rid="bib19">Jin et al., 2013</xref>). Thus, it is likely that the p53 transcriptional program is qualified by the action of lineage-specific factors that prepare a subset of p53 enhancers in a cell type-specific manner.</p><p>Altogether, the results presented here provide a significant advance in our understanding of the p53 transcriptional program and pave the way for functional studies of novel p53 target genes and elucidation of unique regulatory mechanisms within this tumor suppressive gene network.</p></sec><sec id="s4" sec-type="materials|methods"><title>Materials and methods</title><sec id="s4-1"><title>Global run-on deep-sequencing</title><p>Global run-on and library preparation for sequencing were basically done as described in <xref ref-type="bibr" rid="bib14">Hah et al. (2011)</xref>. GRO-seq and microarray datasets are available at Gene Expression Omnibus, data series GSE53966.</p><sec id="s4-1-1"><title>Cell culture</title><p>HCT116 cells were grown in McCoys 5A media and passaged 2 days in a row prior to treatment. We found passaging HCT116 cells twice before the experiment resulted in less clumping of the cells and therefore better nuclei isolation. Cells were plated at a concentration of 10 × 10<sup>6</sup> on 15 cm plates and treated 24 hr later with media containing either Nutlin-3a (10 μM) or the equivalent amount of vehicle (DMSO) for 30 min or 1 hr.</p></sec><sec id="s4-1-2"><title>Nuclei preparation</title><p>Cells were washed 3x with ice cold PBS and then treated with 10 ml per 15 cm plate of ice-cold Lysis Buffer (10 mM Tris–HCl pH 7.4, 2 mM MgCl<sub>2</sub>, 3 mM CaCl<sub>2</sub>, 0.5% NP-40, 10% glycerol, 1 mM DTT, 1x Protease Inhibitor Cocktail Tablets (Roche 11,836 153 001 Germany), 4U/ml SUPERase-In) and scrapped from the plates. Cells were centrifuged 1000×<italic>g</italic> for 7 min at 4°C. Supernatant was removed and pellet was resuspended in 1.5 ml of Lysis Buffer to a homogenous mixture by pipetting 20-30X before adding another 8.5 ml of Lysis Buffer. Suspension was centrifuged with a fixed angle rotor at 1000×<italic>g</italic> for 7 min at 4°C. Supernatant was removed and pellet was resuspended in 1 ml of Lysis Buffer and transfered to a 1.7 ml Eppendorf tube. Suspensions were then pelleted in a microcentrifuge at 1000×<italic>g</italic> for 3 min at 4°C. Next, supernatant was removed and pellets were resuspended in 500 µl of Freezing Buffer (50 mM Tris pH 8.3, 40% glycerol, 5 mM MgCl<sub>2</sub>, 0.1 mM EDTA, 4U/ml SUPERase-In). Nuclei were centrifuged 2000×<italic>g</italic> for 2 min at 4°C. Pellets were resuspended in 100 µl Freezing Buffer. To determine concentration, nuclei were counted from 1 µl of suspension and Freezing Buffer was added to create as many 100 µl aliquots of 5 × 10<sup>6</sup> nuclei as possible. Aliquots were quick frozen in liquid nitrogen and stored at −80°C.</p></sec><sec id="s4-1-3"><title>Nuclear run-on and RNA preparation</title><p>After thawing, each 100 µl aliquot of nuclei was added to 100 µl of Reaction Buffer (10 mM Tris pH 8.0, 5 mM MgCl2, 1 mM DTT, 300 mM KCl, 20 units of SUPERase-In, 1% Sarkosyl, 500 μM ATP, GTP, CTP and Br-UTP) and incubated for 5 min at 30°C. To isolate RNA, 1 ml of Trizol was added to the reaction and vortexed to homogeneity. Samples were split in half and another 500 µl of Trizol added to each half. To isolate RNA, 220 µl chloroform was added to each half sample and samples were centrifuged at max speed for 15 min. Aqueous phase was moved into a new tube and 22.5 µl of 5M NaCl was added. Samples were Acid Phenol-Chloroform extracted twice, then Chloroform extracted once. RNA was then precipitated by adding 1 µl glyco-blue and 3 volumes ice cold ethanol to each sample before storing at −20°C for 20 min or more.</p><sec id="s4-1-3-1"><title>Note on phenol and chloroform extractions</title><p>The current volume of the sample is measured and then an equal volume of Phenol-Chloroform, Chloroform or Acid Phenol-Chloroform is added. Then the mixture is vortexed and centrifuged at 12000×<italic>g</italic> for 15 min (Phenol-Chloroform, Acid Phenol-Chloroform) or 10 min (Chloroform) and the top aqueous layer is kept, the lower organic layer and interphase discarded. Acid Phenol-Chloroform was stored at 4°C but was brought to room temperature before use (∼30 min).</p></sec></sec><sec id="s4-1-4"><title>DNAse treatment and removal of 5′ phosphate groups</title><p>Samples were centrifuged at 12,000×<italic>g</italic> for 10 min washed with 70% ethanol, and then centrifuged at 12,000×<italic>g</italic> for 5 min again. Pellets were air dried for 2 min and resuspended in 20 µl DEPC-treated water. Samples were base-hydrolyzed with 5 µl 1M NaOH on ice for 30 min (creating an average fragment size of 150 nt). Samples were neutralized with 25 µl 1M Tris-Cl pH6.8 and then run through a BioRad P-30 column per manufacturer's protocol. Samples were DNAse-treated in 1x RQ1 DNase buffer and 3 µl DNase I (1unit/μl, # M6101; Promega, Madison, WI) at 37°C for 10 min and then run through a BioRad P-30 column per manufacturer's protocol. To each RNA sample 8.5 μl 10 × antarctic phosphatase buffer, 1 μl of SUPERase-In and 5 μl of antarctic phosphatase was added for 1 hr at 37°C, and then run through a BioRad P-30 column per manufacturer's protocol. Final volume of RNA solution was brought up to 100 µl with DEPC-treated water and 1 µl 500 mM EDTA was added.</p></sec><sec id="s4-1-5"><title>Anti-BrU bead preparation</title><p>To prepare beads, 60 μl Anti-BrU agarose beads (Santa Cruz Biotech, Santa Cruz, CA) were washed twice for 5 min in 500 μl of Binding Buffer (0.5 × SSPE, 1 mM EDTA, 0.05% Tween-20). After each wash buffer was removed after centrifugation at 1000×<italic>g</italic> for 2 min. Beads were then blocked in 500 μl Blocking Buffer (0.5 × SSPE, 1 mM EDTA, 0.05% Tween-20, 0.1% PVP, and 1 mg/ml Ultrapure BSA) for 1 hr. Beads were then washed twice for 5 min each in Binding Buffer. Beads were finally resuspended in 400 µl Binding Buffer.</p></sec><sec id="s4-1-6"><title>Nascent RNA isolation</title><p>All washes and incubations in this section were done with rotation of the tubes. RNA (100 μl) was heated to 65°C for 5 min and kept on ice and added to prepared Anti-BrU beads in 400 µl Binding Buffer for 1 hr at room temperature. BrU-labeled nascent RNA will therefore be attached to the beads at this step. Beads were then washed with multiple wash solutions for 3 min each at room temperature then centrifuged for 2 min at 12,000×<italic>g</italic> and resuspended in the next wash. Beads were washed in 1X Binding Buffer, 1X Low Salt buffer (0.2 × SSPE, 1 mM EDTA, 0.05% Tween-20), 1X High Salt Buffer (0.5% SSPE, 1 mM EDTA, 0.05% Tween-20, 150 mM NaCl) and 2X TET buffer (TE pH 7.4, 0.05% Tween-20). BrU-labeled nascent RNA was eluted at 42°C with 4 × 125 μl of Elution Buffer (5 mM Tris pH 7.5, 300 mM NaCl, 20 mM DTT, 1 mM EDTA and 0.1% SDS). RNA was then Phenol/Chloroform extracted, Chloroform extracted and precipitated with 1.0 µl glyco-blue, 15 μl of 5M NaCl, 3 volumes 100% ethanol at −20°C for more than 20 min.</p></sec><sec id="s4-1-7"><title>PNK treatment and second bead-binding</title><p>Samples were centrifuged for 20 min at 12,000×<italic>g</italic> then washed with 70% ethanol and then pellets were resuspended in 50 μl PNK Reaction Buffer (45 μl of DEPC water, 5.2 μl of T4 PNK buffer, 1 μl of SUPERase_In and 1 μl of T4 PNK [New England BiolabsIpswich, MA]) and incubated at 37 C for 1 hr. To this solution 225 µl water, 5 µl 500 mM EDTA and 18 µl 5M NaCl RNA were added and then the sample was Phenol/Chloroform extracted with 300 µl twice, Chloroform extracted once and precipitated with 3 volumes 100% ethanol at 20°C for more than 20 min. Entire bead binding step was then repeated again to precipitation.</p></sec><sec id="s4-1-8"><title>Reverse transcription</title><p>Reverse transcription was performed as follows: RNA was resuspended in 8.0 μl water and the following was added: 1 μl dNTP mix (10 mM), 2.5 μl oNTI223HIseq primer (12.5 μM) (Sequence: 5′-pGATCGTCGGACTGTAGAACTCT/idSp/CCTTGGCACCCGAGAATTCCATTTTTTTTTTTTTTTTTTTTVN; where p indicates 5′ phosphorylation,/idSp/indicates the 1′,2′-Dideoxyribose modification used to introduce a stable abasic site and VN indicates degenerate nucleotides). This mix was then heated for 3 min at 75°C and chilled briefly on ice. Then 0.5 μl SuperRnaseIn, 3.75 μl 0.1M DTT, 2.5 μl 25 mM MgCl2, 5 μl 5X Reverse Transcription Buffer, and 2 μl Superscript III Reverse Transcriptase were added and the reaction was incubated at 48°C for 30 min. To eliminate excess oNTI223HIseq primer, 4 μl Exonuclease I and 3.2 μl 10X Exonuclease I Buffer were added and the reaction was incubated at 37°C for 1 hr . Finally, RNA was eliminated by adding 1.8 μl 1N NaOH and incubated for 20 min at 98°C. The reaction was then neutralized with 2 μl of 1N HCl. Next, the cDNA was Phenol:Chloroform extracted twice, chloroform extracted once and then precipitated with 300 mM NaCl and 3 volumes of ethanol.</p></sec><sec id="s4-1-9"><title>Size selection</title><p>cDNA was resuspended in 8 μl of water and added to 20 μl FLB (80% Formamide, 10 mM EDTA, 1 mg/ml Xylene Cyanol, 1 mg/ml Bromophenol Blue) before loading on an 8% Urea gel. RNAs between 200–650 nt were selected and gel fragments were shattered, eluted from the gel via rotating overnight in 150 mM NaCl, 1x TE and 0.1% Tween. Entire solution was than ran through Spin X column (CLS8163; Sigma-Corning, Pittston, PA) at 10,000 RPM for 2 min at room temperature. Samples were snap frozen in liquid nitrogen and volume was reduced by spinning in speed-vac until under 500 μl. The cDNA was then Phenol:Chloroform extracted twice, Chloroform extracted once and then precipitated with 300 mM NaCl and 3 volumes of ethanol.</p></sec><sec id="s4-1-10"><title>Circularization and relinearization</title><p>cDNA was resuspended in 8 μl and added to 1 μl Circuligase Buffer, 0.5 μl mM ATP and 0.5 μl of 50 mM MnCl2 and 0.5 μl CicLigase (CL4111K; Epicentre, Madison, WI) and incubated at 60°C for 1 hr. Then reaction was heat-inactivated at 80°C for 20 min. To the above reaction 3.8 μl of 4X relinearization supplement (100 mM KCl, 2 mM DTT) and 1.5 μl ApeI was added and the reaction was incubated at 37°C for 1 hr. The cDNA was then Phenol:Chloroform extracted twice, Chloroform extracted once and then precipitated with 300 mM NaCl and 3 volumes of ethanol.</p></sec><sec id="s4-1-11"><title>PCR and sequencing</title><p>cDNA was resuspended in 15 μl of water and added to 10 μl 5x Phusion HF buffer, 0.5 μl 25 μM HighSeq fwd primer, 0.5 μl 25 μM HighSeq rev primer, 2.5 μl 10 mM dNTPs, 1 μl Phusion DNA polymerase and 20 μl water. PCR amplification was done as follows: Step 1: 98°C for 3 min, Step 2: 98°C for 80 s Step 3: 60°C for 90 s, Step 4: 72°C for 30 s, Step 5: repeat Steps 2–4 for 15–20 cycles, dependent on concentration of the final cDNA. Step 6. 65°C for 10 min, Step 7: 4°C. The entire PCR reaction was then Phenol-Chloroform extracted twice, Chloroform extracted once and precipitated with 300 mM NaCl and 3 volumes of ethanol. The DNA was then run on a 6% native gel (1x TBE, 6% acrylamide [19:1]). Products between 250–650 nt were gel extracted and precipitated as described above. Samples were then sequenced using an Illumina HiSeq2000 instrument with 50 bp reads.</p></sec></sec><sec id="s4-2"><title>Computational analysis</title><p>Unless otherwise noted all processing was done with python version 2.6. Graphs were created with either Microsoft Excel or the python package matplotlib version 1.0.1.</p><sec id="s4-2-1"><title>Mapping</title><p>Two biological replicates of Nutlin-treated and DMSO-treated (control) samples were combined. Single samples were created for p53 null cells and for the 30 min Nutlin treatment. Total reads were as follows p53 +/+ control (DMSO 1 hr) = 214019557 (69399243 + 144620314), p53 +/+ Nutlin 1 hr= 325026803 (166083763 + 158943040), p53 −/− control (DMSO 1 hr) = 184022999, p53 −/− Nutlin 1 hr = 162410974, p53 +/+ Nutlin 30 min = 205211352, p53 +/+ Control (DMSO 30 min) = 175418701. With regards to the quality of the bases in each sample, the median Illumina quality score for all bases was always greater than 32, most of them above 39, that is the inferred base call accuracy was greater than 99.9% for all bases in all samples (Fastx v. 0.0.13.2). Reads were mapped to the human genome (hg19 downloaded from bowtie website) using bowtie (version 0.12.7) with the command bowtie -S -t -v 2 –-best. Using this strategy approximately 70% of the reads mapped. Mapping created a sam file that was then processed by samtools view -bS–o and samtools sort to create a sorted bam file (version 0.1.16).</p></sec><sec id="s4-2-2"><title>Creation of igv figures</title><p>The sorted bam files mentioned above were used to make a Bedgraph by running genomeCoverageBed–bg (v2.12.0) on each strand separately. The negative strand values were changed to negative values in the BedGraph. The BedGraph values were then divided by the number of million of mappable reads. The two files (strands) were concatenated back together and igvtools sort and igvtools tile was used to create a tdf file that was loaded into igv for creation of snapshots of genes (IGVtools 1.5.10, IGV version 2.0.34).</p></sec><sec id="s4-2-3"><title>Calculation of activities and pausing indexes</title><p>Calculations were done exactly as in <xref ref-type="bibr" rid="bib7">Core et al. (2008)</xref> unless otherwise noted. Gene annotations (hg19) were downloaded from: <ext-link ext-link-type="uri" xlink:href="http://hgdownload.cse.ucsc.edu/goldenPath/hg19/database/refGene.txt.gz">http://hgdownload.cse.ucsc.edu/goldenPath/hg19/database/refGene.txt.gz</ext-link>.</p><p>Number of reads in the gene body (1 kb from transcription start site [TSS] to the end of the annotation) and number of reads around the promoter (−100 to +400 bp from annotated TSS) were counted by the program coverageBed v2.12.0. A program to calculate fpkm, pausing indexes, gene activity, and promoter activity was written and run on python 2.6. Fisher's exact test was done using the python module fisher 0.1.4 downloaded from <ext-link ext-link-type="uri" xlink:href="https://pypi.python.org/pypi/fisher/">https://pypi.python.org/pypi/fisher/</ext-link>. RefSeq genes shorter than 1 kb were not used. Genes that are differentially expressed were determined in R version 2.13.0 using DEseq v1.4.1 (<xref ref-type="bibr" rid="bib1">Anders and Huber, 2010</xref>). Settings for DEseq were cds <–estimateSizeFactors(cds), method = 'blind', sharingMode = 'fit-only'. Genes were called as differentially transcribed if they had an adjusted p-value less than or equal to 0.1. Manual curation was used to choose the most parsimonious isoform for the Nutlin vs control (DMSO) comparisons. For genes only differentially expressed across cell lines, we utilized the isoform with the highest fold change (p53+/+ control vs p53 −/− controls). For all other genes we used the isoform identifier with the highest fold change between p53+/+ control and p53+/+ Nutlin.</p></sec><sec id="s4-2-4"><title>Microarray analysis</title><p>HCT116 cells were grown in McCoy's 5A and passaged the day prior to treatment. Cells were plated at a concentration of 300,000 cells per well of six well plate and treated 24 hr later with either Nutlin-3 (10 μM) or the equivalent amount of vehicle (DMSO) for 12 hr. Total RNA from HCT116 cells was harvested with an RNeasy kit (Qiagen, Germantown, MD) and analyzed on Affymetrix HuGene 1.0 ST arrays following the manufacturer's instructions. Microarray data were processed using Partek Genomics Suite 6.6. Anova was used to call differentially expressed genes for which any isoform showed a fold change > |+/−1.5| with FDR <0.05. There were 362 genes called as upregulated and 367 genes as downregulated.</p></sec><sec id="s4-2-5"><title>Comparative analysis of GRO-seq vs microarray data</title><p>The microarray analysis provided a list of gene names and their fold change on the microarray. Since many of the genes had multiple isoforms we simplified by keeping only the isoform with the greatest fold change between Control and Nutlin. For comparisons of microarray and GRO-seq, a list of genes common to both analyses was used. If a gene was found in only one analysis (GRO-seq or microarray) it was not used. In the microarray graphs, expression values from the three biological replicates were averaged. Graphs (MAplot, scatter plot, box and wiskers) were created in python by using matplotlib.</p></sec><sec id="s4-2-6"><title>Meta-analysis of published p53 ChIP-seq data</title><p>To create a list of high confidence p53 binding sites, we combined the data from of 7 ChIP assays for p53 (<xref ref-type="bibr" rid="bib47">Wei et al., 2006</xref>; <xref ref-type="bibr" rid="bib37a">Smeenk et al., 2008</xref>; <xref ref-type="bibr" rid="bib38">Smeenk et al., 2011</xref>; <xref ref-type="bibr" rid="bib29">Nikulenkov et al., 2012</xref>) and kept only sites that were found in at least five of the seven assays. The assays covered three cell lines (HCT116, U20S, MCF7) and 6 different conditions (untreated, Nutlin, 5FU, RITA, ActD, Etoposide). First, all samples were lifted over to hg19 coordinates. Second, a file was created that contained all sites from each experiment and the sites were then sorted and merged using BedTools (2.16.1). This created a list of 106,600 total sites, only 1481 of which were found in 5/7 datasets.</p></sec><sec id="s4-2-7"><title>Meta-analysis of published investigations of the p53 transcriptional program using a combination of microarray and ChIP-seq data</title><p>To generate the Venn diagrams in <xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1B</xref>, we analyzed the information provided in the following papers: (<xref ref-type="bibr" rid="bib29">Nikulenkov et al., 2012</xref>; <xref ref-type="bibr" rid="bib26">Menendez et al., 2013</xref>; <xref ref-type="bibr" rid="bib37">Schlereth et al., 2013</xref>; <xref ref-type="bibr" rid="bib46">Wang et al., 2013</xref>). All gene lists are provided in <xref ref-type="supplementary-material" rid="SD2-data">Supplementary file 2</xref>. Genes bound by p53 and upregulated (339) or downregulated (8) upon Nutlin treatment in MCF7 cells were obtained from Tables S7 and S10 of the Nikulenkov et al paper. Genes bound by p53 and upregulated or downregulated by either Nutlin or Doxorubicin treatment in U2OS cells were obtained from Supplemental Dataset 3, (file: nar-00620-v-2013-File011), of the Menendez et al paper. We removed from the analysis of ‘p53 activated genes’ those that were repressed by Nutlin or Doxorubicin, keeping only those that were either activated by both drugs or activated by one drug and unchanged with the other. Genes bound by p53 and upregulated or downregulated upon p53 overexpression in SAOS2 cells were obtained from Table S1 of the Schlereth et al paper. We removed from this list duplicate gene entries. Genes bound by p53 and upregulated or downregulated upon p53 overexpression in HCT116 p53 −/− cells were obtained from Table S3 and Table S5 of the Wang et al paper. We manually curated all gene lists to ensure that gene nomenclature was the same before analyzing the overlap among the various lists using Venny: <ext-link ext-link-type="uri" xlink:href="http://bioinfogp.cnb.csic.es/tools/venny/index.html">http://bioinfogp.cnb.csic.es/tools/venny/index.html</ext-link>. Note: because some these published studies used somewhat different microarray platforms, it is possible that some genes were not actually measured in all four studies.</p></sec><sec id="s4-2-8"><title>Analysis of eRNAs derived from p53 binding sites</title><p>For each merged ChIP site mentioned above a bed file was created that contain two copies of each site, one assigned to each strand (213,200 ‘stranded-sites’). Then sites were removed if overlapped with an annotated refseq gene on the same strand (160,604 stranded-sites left after removal). The number of reads covering each remaining site were counted using coveragebed and fpkm calculate for each site. To create a list of high confidence p53 binding sites we could use for analyses of eRNAs we used only sites that were bound in five of the seven ChIP experiments used above (2283 stranded-sites, 1454 sites). We then divided these sites into three categories. Distal sites (more than 25 kb from any annotated refseq gene) (992 stranded-sites, 496 sites); proximal to non-target gene sites (less than or equal to 25 kb from a gene i.e. not a GRO-seq target, 1187 stranded-sites, 555 sites with only one strand, 316 sites with both strands); and sites near p53 target genes (less than or equal to 25 kb from a gene i.e. a GRO-seq target, 104 stranded-sites, 70 sites with only one strand, 17 sites with both strands). We then determined a percent of the sites that are active by the same method used in <xref ref-type="bibr" rid="bib7">Core et al. (2008)</xref>. Briefly, background was set at 3% of the reads. We estimate the probability of observing at least N reads in an interval of length L using a Poisson distribution. If the probability was above 0.01 the site was called transcriptionally active. For each site, if either of the two strands was transcriptionally active (logical or) the site was counted as active therefore for the <xref ref-type="fig" rid="fig5s1">Figure 5—figure supplement 1E</xref>, n = 496, 496, 496, 496, 871, 871, 871, 871, 87, 87, 87, 87. For <xref ref-type="fig" rid="fig5">Figure 5D</xref> we wanted to include information from both strands when available so stranded-sites were used to determine the fpkm of each site, therefore n = 992, 992, 992, 992, 1187, 1187, 1187, 1187, 104, 104, 104, 104. Similarly, for <xref ref-type="fig" rid="fig4s1">Figure 4—figure supplement 1F</xref> we included all stranded-sites but for each comparison we had to remove any sites in a given sample that had 0 reads, therefore n = 323, 602, 566, 897, 83, 100.</p><sec id="s4-2-8-1"><title>Note on distance to the nearest binding site</title><p>To determine the distance to the nearest p53 binding site, for all genes the program pybedtools and the script closest was used. The sites were the 1481 sites that were in five of seven ChIP experiments as mentioned above. The target genes were the 202 up and down-regulated genes by GRO-seq. The distances were then binned to create the histogram shown in <xref ref-type="fig" rid="fig5">Figure 5B</xref>. The 10 most outer bins on the left and right are bins of 5 kb; the inner bins are bins of 1 kb.</p></sec></sec><sec id="s4-2-9"><title>Overlap of genes downregulated in microarray and miR-34a targets</title><p>The published genes that were downregualted upon miR-34a overexpression in HCT116 cells (<xref ref-type="bibr" rid="bib21">Lal et al., 2011</xref>) (2091 total, 1765 also found in our microarray experiment) were compared to the genes that were downregulated upon Nutlin treatment in our microarray experiment (367 total, 342 also found in the published microarray by <xref ref-type="bibr" rid="bib21">Lal et al. (2011)</xref>). Of the 342, 245 (72%) were downregulated in the miR-34a overexpression experiment. All genes which overlap (16,553) between the two microarrays (miR-34a overexpression n = 21,050 and Nutlin treatment, n = 19,901) were determined assuming Lal et al used the annotations from version 32 of Affymetrix U133 plus 2.0 mRNA microarray. Hypergeometric was used to calculate a p-value.</p></sec><sec id="s4-2-10"><title>Ingenuity pathway analysis (IPA) of genes downregulated upon Nutlin treatment</title><p>The 367 genes shown to be downregulated upon Nutlin treatment in our microarray experiment (‘Microarray analysis’) were subjected to IPA Upstream Regulator Analysis, which identifies upstream transcriptional regulators that can explain the observed gene expression changes in a user's dataset. The top three upstream transcriptional regulators identified in our dataset were E2F4, CDKN1A and RB, all three identified as ‘transcriptional repressors’ by this analysis. Statistical significance and p-values were determined by IPA using a Fisher's Exact Test. Detailed explanation of this analysis is provided by IPA at: <ext-link ext-link-type="uri" xlink:href="http://pages.ingenuity.com/IngenuityUpstreamRegulatorAnalysisWhitepaper.html">http://pages.ingenuity.com/IngenuityUpstreamRegulatorAnalysisWhitepaper.html</ext-link>.</p></sec><sec id="s4-2-11"><title>Oncomine analysis of p53 wild type and p53 mutant cell lines and tumors</title><p>Oncomine (Compendia Bioscience, Ann Arbor, MI) was used for analysis and visualization of expression data from the Garnett Cell Line dataset (<xref ref-type="bibr" rid="bib12">Garnett et al., 2012</xref>) containing gene expression data for hundreds of p53 wild type and p53 mutant cancer cell lines or the Ivshina Breast Carcinoma dataset (<xref ref-type="bibr" rid="bib17">Ivshina et al., 2006</xref>). Filters in the Oncomine database were set to select Garnett Cell Line dataset (or the Ivshina Breast Carcinoma dataset), and <italic>TP53</italic> mutation status. Genes analyzed were individually filtered with the above filters in place, and analysis was adjusted to view under-expression with p53mt vs p53wt (for basal up with p53 genes) and adjusted to view over-expression with p53mt vs p53wt (for basal down with p53 genes). Graphs represent the log2 median centered intensity (MCI) of gene expression from selected genes on the array. The p-value is derived from a Student's <italic>t</italic> test between p53wt and p53mt expression, and the gene rank reveals the rank of that gene in the dataset among all other genes in the analysis according to the p-value for the differential expression analysis.</p></sec><sec id="s4-2-12"><title>Gene set enrichment analysis (GSEA) of GRO-seq p53 target genes analyzed by Oncomine</title><p>To generate the GSEA plots shown in <xref ref-type="fig" rid="fig3s2">Figure 3—figure supplement 2A</xref>, we analyzed the position of 183 GRO-seq p53 target genes (out of 198) which were also present in the gene rank generated by Oncomine analysis of the Garnett Cell Line dataset (see above) using GSEA software from the Broad Institute: (<ext-link ext-link-type="uri" xlink:href="http://www.broadinstitute.org/gsea/index.jsp">http://www.broadinstitute.org/gsea/index.jsp</ext-link>). In this analysis, the ‘molecular profile data’ was the Oncomine gene rank and the ‘gene set database’ was the list the GRO-seq p53 target genes. The analysis was repeated using the ‘microarray only’ genes (i.e., genes upregulated by Nutlin as seen only in our microarray experiment). For the scatter plot log2 fold changes in the Oncomine dataset were plotted against fold changes in GRO-seq data using R.</p></sec></sec><sec id="s4-3"><title>RT-PCR</title><p>Total RNA was isolated using Trizol (Life Technologies, Frederick, MD) following the manufacturer's instructions. cDNA was generated using qScript kit (Quanta Biosciences, Gaithersburg, MD) with random priming. cDNA was subjected to standard or quantitative PCR (Q-PCR) using SYBR green or SYBR select master mix on a Viia 7 instrument (Life Technologies) with the primer pairs listed below. The 18S rRNA was used for normalization. Experiments were done in biological triplicate and error bars represent the SEM.</p></sec><sec id="s4-4"><title>Primers used for Q-RT-PCR</title><p><list list-type="simple"><list-item><p>Gene, Forward Primer, Reverse Primer</p></list-item><list-item><p>Gene, Forward Primer, Reverse Primer</p></list-item><list-item><p><bold>18S rRNA,</bold> GCCGCTAGAGGTGAAATTCTTG, CTTTCGCTCTGGTCCGTCTT</p></list-item><list-item><p><bold>ADAMTS7</bold>, GTCGAGCCTCCCCGCTGTG, CAGCGTCTCGCAGAACCCGA</p></list-item><list-item><p><bold>ALOX5</bold>, AAACGAGCTGTTCCTGGGCATGT, GGCCTCGAGGTTCTTGCGGA</p></list-item><list-item><p><bold>APOBEC3C</bold>, AGTATCCATGTTACCAGGAGGGGCT, TTTAAAATCTTCATAGTCCATGATCTCCACAGCGA</p></list-item><list-item><p><bold>ASCC3</bold>, GCTTGCAGAGAGTCCACTTTGGGT, ACAGCTCTCTGGCTTTCCCTTGTG</p></list-item><list-item><p><bold>ASS1</bold>, TGAGGAGCTGGTGAGCATGAACG, ACCCGGTGGCATCAGTTGGC</p></list-item><list-item><p><bold>ASTN2</bold>, TTTTCTGCCGCAGCGAGGAGG, AGAGTCAAATAATATACGTGATTTTGGTGTCCTTGA</p></list-item><list-item><p><bold>BLOC1S2</bold>, CCGAGGGCGTACTGGCGAC, GGCATCGTCTCGGGCGGG</p></list-item><list-item><p><bold>C19ORF82 (LOC284385),</bold> ATTTAGAGGAGGGACCCGGC, AAGCTTTGGAGGAGCATCCC</p></list-item><list-item><p><bold>CDC42BPG</bold>, TGTCCTGCCCCCAGGGATCG, GGGCCGTTCTGAGACCTGCATTAG</p></list-item><list-item><p><bold>CDKN1A</bold>, CTGGAGACTCTCAGGGTCGAAA, GATTAGGGCTTCCTCTTGGAGAA</p></list-item><list-item><p><bold>CEP85L</bold>, AGCTCCTTGGCAACAGCAGCAA, TGATCCAGTCAAAACCTGCATTGGTGT</p></list-item><list-item><p><bold>CHAC1</bold>, TGAAGATCATGAGGGCTGCACTTGG, CAGGGCCTTGCTTACCTGCTCC</p></list-item><list-item><p><bold>DDB2</bold>, TCTACTCGCTGCCGCACAGG, TCGGGACTGAAACAAGCTGCGT</p></list-item><list-item><p><bold>EFNB1,</bold> TGGGCAAGATCCCAATGCTGTG, TGCTTGCCATCAGAGTCACCCA</p></list-item><list-item><p><bold>FAM210B</bold>, TGGCACTGTTGGCGTGTCAT, CAGGCATGTCCACACCACTTGA</p></list-item><list-item><p><bold>FAM212B</bold>, AATGGGCGCATGACGATGGAGA, TGCAGTGCACCCATCATGCAGT</p></list-item><list-item><p><bold>GDF15</bold>, TAACCAGGCTGCGGGCCAAC, CAGCCGCACTTCTGGCGTGA</p></list-item><list-item><p><bold>GJB5</bold>, CTGCTGGGAGCCAGGAGAGC, CGCGTCAGCTGCTACTGAGTGA</p></list-item><list-item><p><bold>GPR87</bold>, AACCTATGCTGAACCCACGCCT, GCCGTGCAGCTCGTTATTTGGT</p></list-item><list-item><p><bold>GRIN2C</bold>, ACCGTGACATGCACACCCACAT, TGAAGGCATCCAGCTTCCCCAT</p></list-item><list-item><p><bold>HES2</bold>, CTGGGCCGGGAGAACTCCAA, GCAGGAAGCGCACGGTCATT</p></list-item><list-item><p><bold>KCNN4</bold>, CCGCCATCAACGCGTTCCG, CCCGGAGCTTCCGGTGTTTCA</p></list-item><list-item><p><bold>LRP1</bold>, ACCCACTGAGGGGGACCATGT, TCCCATCCATCGCTGCCGTC</p></list-item><list-item><p><bold>PHLDA3</bold>, TAAACCACCGGGCGCACCAT, CAGAGGGAACAACGAAGCTGCC</p></list-item><list-item><p><bold>PTP4A1</bold>, TAAGACAAAAGCGGCGTGGAGC, ACGCAGCCGCATTTTAGGACG</p></list-item><list-item><p><bold>PTPRE</bold>, CTTTTCCCGGCTCACCTGGTTCAG, GGGCATCTTCTTGTCGCTGGTG</p></list-item><list-item><p><bold>RINL</bold>, TGGCTTCCTGCAACCTGACGAT, TGCTTTGTTCACCTGTGGTGGG</p></list-item><list-item><p><bold>SCL30A1</bold>, TTGGACCCCGCAGACCCAGA, TGGTCAGGTTCTCTGACAAGATTTCCATT</p></list-item><list-item><p><bold>SULF2</bold>, CCCCCGGACTCGAAACATGGA, ACTTTCGACGCTGAAACTGCCTGT</p></list-item><list-item><p><bold>TM7SF3</bold>, CTTTCACAAATGTGCCTTTTCAAACTAATGACTTC, AGCATGCCCCATACTGCCAGG</p></list-item><list-item><p><bold>TOB1</bold>, GAAGCAGCCCCTTCCCAGGAG, GCTTTTCAGGATACCAGTGCCCTTCAT</p></list-item><list-item><p><bold>TP53INP1</bold>, AATGAGAAAGAAGATGATGAATGGATT, TGCTGAGAAACCAGTGCAAGTATC</p></list-item><list-item><p><bold>WRAP73</bold>, TGCCTGGGGAAGGCGACTTTG, GAGGGCCATCGAGTCTCCGC</p></list-item></list></p></sec><sec id="s4-5"><title>Western blotting</title><p>Protein extracts were separated by SDS-PAGE and transferred to PVDF membranes. Blots were probed with the primary antibodies listed below, followed by peroxidase-conjugated secondary antibodies (Santa Cruz Biotechnologies). Detection was by enhanced chemiluminescence or Luminata Crescendo (Millipore, Tauton, MA) and images were captured using an ImageQuant LAS4000 digital camera system (GE Healthcare, Sweden). Antibodies: p21 (sc-817; Santa Cruz); p53 (DO-1, OP43; Calbiochem, Tauton, MA); MDM2 (SMP14, sc-965; Santa cruz, used in combination with True Blot anti-mouse secondary antibody (Rockland, Boyertown, PA) to avoid IgG signal in immunoprecipitated samples); GRIN2C (LS-C157785; LifeSpan BioSciences, Seattle, WA); PTCDH4 (sc-139478; Santa Cruz (C-15)); RINL (OAAB08343; Aviva Systems Biology, San Diego, CA); Nucleolin (sc-8031; Santa Cruz); α-tubulin (T9026; Sigma, Israel).</p></sec><sec id="s4-6"><title>BrdU incorporation assay</title><p>Cell proliferation was assessed by quantification of cells actively synthesizing DNA. HCT116 p53+/+ and HCT116 p53−/− cells treated with 10 μM Nutlin-3 as indicated. 60 min prior to harvesting cells, 1 mg/ml of BrdU was added to cultivation medium. After trypsinization, cells were fixed with 70% ethanol. DNA denaturation was achieved by incubation with 2M HCl with 0.5% Triton X-100 for 10 min at 37°C. Before staining with anti-BrdU antibody (sc-51514; Santa Cruz Biotechnology), samples were neutralized with 0.1M Na<sub>2</sub>B<sub>4</sub>O<sub>7</sub>. Following incubation with secondary antibody (A21202, Life Technologies). BrdU positivity was analyzed by flow cytometry. Data represent average numbers of BrdU positive cells from three independent experiments (each run in duplicate). <italic>T</italic> test was used for statistical analysis; an asterisk denotes significance at p<0.001.</p></sec><sec id="s4-7"><title>ChIP assays</title><p>Chromatin immunoprecipitation (ChIP) was carried out as previously described (<xref ref-type="bibr" rid="bib13">Gomes et al., 2006</xref>). Briefly, for quantitative ChIP analysis of the <italic>GDF15 and PTP4A1</italic> locus, cells were cross-linked with 1% formaldehyde and whole-cell lysates were prepared in RIPA buffer. 1 mg of protein extract was immunoprecipitated with specific antibodies. Real-time PCR was carried out on ChIP-enriched DNA against a standard curve of genomic DNA, with amplicons tiling across the locus. Enrichment values for each amplicon were calculated as percentage of the amplicon with maximum signal for each antibody or expressed as equivalents of gDNA in nanograms. For MDM2 ChIPs a modified RIPA pH 7.4 buffer was used (150 mM NaCl, 1% Nonidet P-40, 0.5% sodium deoxycholate, 0.025% SDS, 50 mM Tris pH 7.6, 5 mM EDTA, protease/phosphatase inhibitors).</p><sec id="s4-7-1"><title>Antibodies used for ChIP</title><p>RNAPII S5P (H14): Covance (Dedham, MA) MMS134R; RNAPII S2P (H5): Covance MMS129R; p53: Calbiochem OP43; MDM2: Santa Cruz sc-13161 (D-7); normal mouse IgG: Santa Cruz sc-2025.</p></sec></sec><sec id="s4-8"><title>Primers used for ChIP</title><p>Location of amplicon center relative to TSS, forward primer, reverse primer (5′>3′):</p><sec id="s4-8-1"><title>GDF15 tiling:</title><p><list list-type="simple"><list-item><p>GDF15 -692, AAGGTGAGCGCTGAGCCAGA, AGCAGGGTGCTCCCAAACCC</p></list-item><list-item><p>GDF15 +54, CCAGCTCAGAGCCGCAACCT, CAGCCACGAGAGCACCAGCA</p></list-item><list-item><p>GDF15 +851, GAGGCTCCAGGGGTGCAGG, GTCCCACCACAGCCCAGACC</p></list-item><list-item><p>GDF15 +1596, GGGTGCTTCCTATTTTAGGCAGGCT, ATGGCTGGCTCCTCCTCTTCCT</p></list-item><list-item><p>GDF15 +2236, CGCCTTCACCGGGCTCTGTT, CGTGTCACGTCCCACGACCT</p></list-item><list-item><p>GDF15 +2991, TCGGGGGCTGGTCTGATGGA, CCTGGACACCACAGGGAACAGT</p></list-item><list-item><p>GDF15 +4107, TTATGCCAGGCCCCGGACAC, CCGCAAAGTGGGGTGGTCCT</p></list-item></list></p></sec><sec id="s4-8-2"><title>PTP4A1 tiling:</title><p><list list-type="simple"><list-item><p>PTP4A1 -1524, AGGCTAAGACCCCCAGGCCA, ACCACAAGCAGTGCTGACACGA</p></list-item><list-item><p>PTP4A1 +16, GGGCCCAGGTGCCAAGGTAA, GGCTGGAACCGCGTCTCAGT</p></list-item><list-item><p>PTP4A1 +2484, TGGAGGAGTGTCTTGGGCGG, TGGAGCCACCTAGTGGATGGGG</p></list-item><list-item><p>PTP4A1 +4564, TGCTCCACCAAGAAGCCCCC, ACACTTCTGTGTCCAGGATAACCACTC</p></list-item><list-item><p>PTP4A1 +7208, ACTTCTGAGTTGGGGAATGTTTGGAGAA, AGTACTGGAGCTCTGGTGAGGTAAGAA</p></list-item><list-item><p>PTP4A1 +9477, ACATGCATTTATCATCATTGAGATTGGTTTGCT, AGATGAGGACACACATTAAGTGAAGTGCAA</p></list-item><list-item><p>PTP4A1 +11301, TGCCATACATGTTAATATTCTACATTCTTGCTTCCT, TCAATTCAGGTGTGAGGCACATATATACACA</p></list-item><list-item><p>PTP4A1 +12,375, AGCCACTTGGGACAAGTCAATGCT, CGTTTTGTGCTGTGTAGGAAATACCGA</p></list-item></list></p><p>p53REs at various gene loci:<list list-type="simple"><list-item><p>CDK1NA (p21) −2283, AGCAGGCTGTGGCTCTGATT, CAAAATAGCCACCAGCCTCTTCT</p></list-item><list-item><p>PTP4A1 +1789, AAATGGCCTGGTTCGGAGCGT, ATGTATGCTGGTGCCCGGACG</p></list-item><list-item><p>HES2 +5170, TCTCTGGGGGCTGAGTGGGAGT, CGGGCAGCTCTCTGAAGGGCA</p></list-item><list-item><p>CDC42BPG +1373, CGGCATGGGGGCCCTCTGTT, CTCTGGCCTCCCACAGGGCAT</p></list-item><list-item><p>ADAMTS7 -48815, GTAAGGCGGGCCAAGGAAGGC, TTCAAGGGGTGGGCGCTCTG</p></list-item><list-item><p>LRP1 -694, GGGCGGGCGGGTACTAAGGT, CAGGGACACCCGAGGGGACA</p></list-item><list-item><p>TOB1 -38453, GCTGGCTACGGGCTGGTCTC, CAGGTGCTGGCTTGCTGCCC</p></list-item><list-item><p>ASS1 +116325, GTGGCTAGTGGGAGGGGGCAT, TGACAAGGGCCAGTGCTGAGGA</p></list-item><list-item><p>CEP85L +77725, AGGAGGGGAGGAGGAGATTGTTGC, GCATGTGTTCCTCAGCGTAATCTGT</p></list-item></list></p></sec></sec><sec id="s4-9"><title>p53 N-terminus detection by flow cytometry</title><p>HCT116 cells were fixed and stained as described previously (<xref ref-type="bibr" rid="bib1a">Andrysik, 2013</xref>). Briefly, trypsinized cells were fixed with methanol and incubated with blocking solution (1% BSA-V in PBS/0.1% Triton X-100, 0.01% sodium azide, pH 7.2). p53 level was detected using DO-1 antibody (Calbiochem OP43) and Alexa Fluor A488 secondary antibody. At least 10,000 single cells were measured by Accuri C6 flow cytometer to asses signal distribution in cellular population.</p></sec><sec id="s4-10"><title>Immunofluorescence</title><p>Following the Nutlin treatment period, HCT116 cells cultivated on glass cover slips were fixed with cold methanol and washed briefly with wash buffer (PBS with 0.1% Triton X-100). After incubation with blocking solution (1% BSA-V in wash buffer, 0.01% sodium azide, pH 7.2) cells were labeled with p53 antibody (sc-6343X; 1:1000; Santa Cruz FL393) or isotype control. Three times washed cells were incubated with secondary antibody (1:2000; Alexa Fluor A21202; Invitrogen, Frederick, MD), nuclei were counterstained with DAPI (Invitrogen, 0.2 μg/ml) and cover slips were mounted using anti-fade medium (Fluoromount, Diagnostic Biosystems, Pleasanton, CA). Images were taken using Zeiss 510 Confocal Laser Scanning Microscope.</p></sec></sec></body><back><ack id="ack"><title>Acknowledgements</title><p>This work was supported primarily by NIH RO1CA117907 grant (JME), NSF grant MCB1243522 (JME), a Butcher Foundation Seed Grant (JME, RDD), the Boettcher Foundation's Webb-Waring Biomedical Research program (RDD), a NIH training grant N 2T15 LM009451 (MAA), a post-doctoral fellowship from the American Cancer Society (HM). JME is an HHMI Early Career Scientist. We thank Dr Honnold and Dr Caldwell for their inspiring work.</p></ack><sec sec-type="additional-information"><title>Additional information</title><fn-group content-type="competing-interest"><title>Competing interests</title><fn fn-type="conflict" id="conf1"><p>JME: Reviewing editor, <italic>eLife</italic>.</p></fn><fn fn-type="conflict" id="conf2"><p>The other authors declare that no competing interests exist.</p></fn></fn-group><fn-group content-type="author-contribution"><title>Author contributions</title><fn fn-type="con" id="con1"><p>MAA, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con2"><p>HSM, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con3"><p>JME, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con4"><p>ZA, Acquisition of data, Analysis and interpretation of data</p></fn><fn fn-type="con" id="con5"><p>VLD, Acquisition of data, Analysis and interpretation of data</p></fn><fn fn-type="con" id="con6"><p>AG, Acquisition of data, Analysis and interpretation of data</p></fn><fn fn-type="con" id="con7"><p>JAF, Conception and design, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con8"><p>MDG, Conception and design, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con9"><p>RDD, Conception and design, Analysis and interpretation of data, Drafting or revising the article</p></fn><fn fn-type="con" id="con10"><p>KDS, Conception and design, Acquisition of data, Analysis and interpretation of data</p></fn><fn fn-type="con" id="con11"><p>XL, Conception and design, Drafting or revising the article, Contributed unpublished essential data or reagents</p></fn><fn fn-type="con" id="con12"><p>WLK, Conception and design, Drafting or revising the article, Contributed unpublished essential data or reagents</p></fn></fn-group></sec><sec sec-type="supplementary-material"><title>Additional files</title><supplementary-material id="SD1-data"><object-id pub-id-type="doi">10.7554/eLife.02200.016</object-id><label>Supplementary file 1.</label><caption><p>Genes upregulated at the transcriptional level in HCT116 p53 +/+ cells treated with 10 μM Nutlin-3a for 1 hr as detected by GRO-seq (198 genes). DeSeq algorithm was used to detect annotated gene loci whose GRO-seq signal was statistically different between DMSO- and Nutlin-treated cells (adjusted p<0.1). Columns in this table indicate: (a) Gene name, (b) Whether the gene was previously identified as a direct p53 target gene in the literature, (c–f) Whether the gene was predicted to be a direct p53 target gene by one or more recent studies employing ChIP-seq and microarrays (<xref ref-type="fig" rid="fig2s1 fig2s2">Figure 2—figure supplements 1 and 2</xref>), (g) fpkm in p53 +/+ control, (h) fpkm in p53 +/+ Nutlin-3, (i) Fold induction, (j) Protein Function, (k) Putative downstream pathway within the p53 network, (l) References describing the gene as a direct target, putative target or establishing gene function.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.02200.016">http://dx.doi.org/10.7554/eLife.02200.016</ext-link></p></caption><media mime-subtype="xlsx" mimetype="application" xlink:href="elife02200s001.xlsx"/></supplementary-material><supplementary-material id="SD2-data"><object-id pub-id-type="doi">10.7554/eLife.02200.017</object-id><label>Supplementary file 2.</label><caption><p>Lists of genes bound by p53 as defined by ChIP-seq and concurrently upregulated or downregulated as defined by microarray measurements of RNA steady state levels. Related to <xref ref-type="fig" rid="fig2s1">Figure 2—figure supplement 1A,B</xref>. See ‘Materials and methods’, ‘Computational Analysis–Meta-analysis of published investigations of the p53 transcriptional program using a combination of microarray and ChIP-seq data’ for details.</p><p><bold>DOI:</bold> <ext-link ext-link-type="doi" xlink:href="10.7554/eLife.02200.017">http://dx.doi.org/10.7554/eLife.02200.017</ext-link></p></caption><media mime-subtype="xlsx" mimetype="application" xlink:href="elife02200s002.xlsx"/></supplementary-material><sec sec-type="datasets"><title>Major datasets</title><p>The following dataset was generated:</p><p><related-object content-type="generated-dataset" document-id="Dataset ID and/or url" document-id-type="dataset" document-type="data" id="dataro1"><name><surname>Allen</surname><given-names>Mary Ann</given-names></name>, <name><surname>Mellert</surname><given-names>Hestia</given-names></name>, <name><surname>Dengler</surname><given-names>Veronica</given-names></name>, <name><surname>Andrysik</surname><given-names>Zdenek</given-names></name>, <name><surname>Guarnieri</surname><given-names>Anna</given-names></name>, <name><surname>Freeman</surname><given-names>Justin A</given-names></name>, <name><surname>Luo</surname><given-names>Xin</given-names></name>, <name><surname>Kraus</surname><given-names>W Lee</given-names></name>, <name><surname>Dowell</surname><given-names>Robin D</given-names></name> and <name><surname>Espinosa</surname><given-names>Joaquin M</given-names></name>, <year>2014</year><x>, </x><source>Global analysis of p53-regulated transcription reveals its direct targets and unexpected regulatory mechanisms</source><x>, </x><ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE53966">http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE53966</ext-link><x>, </x><comment>Publicly available at NCBI Gene Expression Omnibus.</comment></related-object></p><p>The following previously published datasets were used:</p><p><related-object content-type="generated-dataset" document-id="Dataset ID and/or url" document-id-type="dataset" document-type="data" id="dataro2"><name><surname>Nikulenkov</surname><given-names>F</given-names></name>, <name><surname>Spinnler</surname><given-names>C</given-names></name>, <name><surname>Li</surname><given-names>H</given-names></name>, <name><surname>Tonelli</surname><given-names>C</given-names></name>, <name><surname>Shi</surname><given-names>Y</given-names></name>, <name><surname>Turunen</surname><given-names>M</given-names></name>, <name><surname>Kivioja</surname><given-names>T</given-names></name>, <name><surname>Ignatiev</surname><given-names>I</given-names></name>, <name><surname>Kel</surname><given-names>A</given-names></name>, <name><surname>Taipale</surname><given-names>J</given-names></name>, <name><surname>Selivanova</surname><given-names>G</given-names></name>, <year>2012</year><x>, </x><source>Microarray and ChIP-seq data from Insights into p53 transcriptional function via genome-wide chromatin occupancy and gene expression analysis</source><x>, </x><object-id pub-id-type="art-access-id">SRP007261</object-id><x>; </x><ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/sra/SRP007261">http://www.ncbi.nlm.nih.gov/sra/SRP007261</ext-link><x>, </x><comment>Publicly available at the NCBI Sequence Read Archive (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/sra">http://www.ncbi.nlm.nih.gov/sra</ext-link>).</comment></related-object></p><p><related-object content-type="generated-dataset" document-id="Dataset ID and/or url" document-id-type="dataset" document-type="data" id="dataro3"><name><surname>Garnett</surname><given-names>MJ</given-names></name>, <name><surname>Edelman</surname><given-names>EJ</given-names></name>, <name><surname>Heidorn</surname><given-names>SJ</given-names></name>, <name><surname>Greenman</surname><given-names>CD</given-names></name>, <name><surname>Dastur</surname><given-names>A</given-names></name>, <name><surname>Lau</surname><given-names>KW</given-names></name>, <name><surname>Greninger</surname><given-names>P</given-names></name>, <name><surname>Thompson</surname><given-names>IR</given-names></name>, <name><surname>Luo</surname><given-names>X</given-names></name>, <name><surname>Soares</surname><given-names>J</given-names></name>, <name><surname>Liu</surname><given-names>Q</given-names></name>, <name><surname>Iorio</surname><given-names>F</given-names></name>, <name><surname>Surdez</surname><given-names>D</given-names></name>, <name><surname>Chen</surname><given-names>L</given-names></name>, <name><surname>Milano</surname><given-names>RJ</given-names></name>, <name><surname>Bignell</surname><given-names>GR</given-names></name>, <name><surname>Tam</surname><given-names>AT</given-names></name>, <name><surname>Davies</surname><given-names>H</given-names></name>, <name><surname>Stevenson</surname><given-names>JA</given-names></name>, <name><surname>Barthorpe</surname><given-names>S</given-names></name>, <name><surname>Lutz</surname><given-names>SR</given-names></name>, <name><surname>Kogera</surname><given-names>F</given-names></name>, <name><surname>Lawrence</surname><given-names>K</given-names></name>, <name><surname>McLaren-Douglas</surname><given-names>A</given-names></name>, <name><surname>Mitropoulos</surname><given-names>X</given-names></name>, <name><surname>Mironenko</surname><given-names>T</given-names></name>, <name><surname>Thi</surname><given-names>H</given-names></name>, <name><surname>Richardson</surname><given-names>L</given-names></name>, <name><surname>Zhou</surname><given-names>W</given-names></name>, <name><surname>Jewitt</surname><given-names>F</given-names></name>, <name><surname>Zhang</surname><given-names>T</given-names></name>, <name><surname>O'Brien</surname><given-names>P</given-names></name>, <name><surname>Boisvert</surname><given-names>JL</given-names></name>, <name><surname>Price</surname><given-names>S</given-names></name>, <name><surname>Hur</surname><given-names>W</given-names></name>, <name><surname>Yang</surname><given-names>W</given-names></name>, <name><surname>Deng</surname><given-names>X</given-names></name>, <name><surname>Butler</surname><given-names>A</given-names></name>, <name><surname>Choi</surname><given-names>HG</given-names></name>, <name><surname>Chang</surname><given-names>JW</given-names></name>, <name><surname>Baselga</surname><given-names>J</given-names></name>, <name><surname>Stamenkovic</surname><given-names>I</given-names></name>, <name><surname>Engelman</surname><given-names>JA</given-names></name>, <name><surname>Sharma</surname><given-names>SV</given-names></name>, <name><surname>Delattre</surname><given-names>O</given-names></name>, <name><surname>Saez-Rodriguez</surname><given-names>J</given-names></name>, <name><surname>Gray</surname><given-names>NS</given-names></name>, <name><surname>Settleman</surname><given-names>J</given-names></name>, <name><surname>Futreal</surname><given-names>PA</given-names></name>, <name><surname>Haber</surname><given-names>DA</given-names></name>, <name><surname>Stratton</surname><given-names>MR</given-names></name>, <name><surname>Ramaswamy</surname><given-names>S</given-names></name>, <name><surname>McDermott</surname><given-names>U</given-names></name>, <name><surname>Benes</surname><given-names>CH</given-names></name>, <year>2011</year><x>, </x><source>Gene expression analysis of 789 cancer cell lines using the Affymetrix HT- HG-U133A v2 platform</source><x>, </x><object-id pub-id-type="art-access-id">E-MTAB-783</object-id><x>; </x><ext-link ext-link-type="uri" xlink:href="http://www.ebi.ac.uk/arrayexpress/experiments/E-MTAB-783/">http://www.ebi.ac.uk/arrayexpress/experiments/E-MTAB-783/</ext-link><x>, </x><comment>Publicly available at ArrayExpress (<ext-link ext-link-type="uri" xlink:href="http://www.ebi.ac.uk/arrayexpress">http://www.ebi.ac.uk/arrayexpress</ext-link>).</comment></related-object></p><p><related-object content-type="generated-dataset" document-id="Dataset ID and/or url" document-id-type="dataset" document-type="data" id="dataro4"><name><surname>Smeenk</surname><given-names>L</given-names></name>, <name><surname>van Heeringen</surname><given-names>SJ</given-names></name>, <name><surname>Koeppel</surname><given-names>M</given-names></name>, <name><surname>van Driel</surname><given-names>MA</given-names></name>, <name><surname>Bartels</surname><given-names>SJ</given-names></name>, <name><surname>Akkers</surname><given-names>RC</given-names></name>, <name><surname>Denissov</surname><given-names>S</given-names></name>, <name><surname>Stunnenberg</surname><given-names>HG</given-names></name>, <name><surname>Lohrum</surname><given-names>M</given-names></name>, <year>2008</year><x>, </x><source>Chromatin immunoprecipitation of p53 in human osteocarcoma cells</source><x>, </x><object-id pub-id-type="art-access-id">E-TABM-442</object-id><x>; </x><ext-link ext-link-type="uri" xlink:href="http://www.ebi.ac.uk/arrayexpress/experiments/E-TABM-442/">http://www.ebi.ac.uk/arrayexpress/experiments/E-TABM-442/</ext-link><x>, </x><comment>Publicly available at ArrayExpress (<ext-link ext-link-type="uri" xlink:href="http://www.ebi.ac.uk/arrayexpress">http://www.ebi.ac.uk/arrayexpress</ext-link>).</comment></related-object></p><p><related-object content-type="generated-dataset" document-id="Dataset ID and/or url" document-id-type="dataset" document-type="data" id="dataro5"><name><surname>Wei</surname><given-names>CL</given-names></name>, <name><surname>Wu</surname><given-names>Q</given-names></name>, <name><surname>Vega</surname><given-names>VB</given-names></name>, <name><surname>Chiu</surname><given-names>KP</given-names></name>, <name><surname>Ng</surname><given-names>P</given-names></name>, <name><surname>Zhang</surname><given-names>T</given-names></name>, <name><surname>Shahab</surname><given-names>A</given-names></name>, <name><surname>Yong</surname><given-names>HC</given-names></name>, <name><surname>Fu</surname><given-names>Y</given-names></name>, <name><surname>Weng</surname><given-names>Z</given-names></name>, <name><surname>Liu</surname><given-names>J</given-names></name>, <name><surname>Zhao</surname><given-names>XD</given-names></name>, <name><surname>Chew</surname><given-names>JL</given-names></name>, <name><surname>Lee</surname><given-names>YL</given-names></name>, <name><surname>Kuznetsov</surname><given-names>VA</given-names></name>, <name><surname>Sung</surname><given-names>WK</given-names></name>, <name><surname>Miller</surname><given-names>LD</given-names></name>, <name><surname>Lim</surname><given-names>B</given-names></name>, <name><surname>Liu</surname><given-names>ET</given-names></name>, <name><surname>Yu</surname><given-names>Q</given-names></name>, <name><surname>Ng</surname><given-names>HH</given-names></name>, <name><surname>Ruan</surname><given-names>Y</given-names></name>, <year>2006</year><x>, </x><source>p53 ChIP data from A global map of p53 transcription-factor binding sites in the human genome</source><x>, </x><ext-link ext-link-type="uri" xlink:href="http://hgdownload.cse.ucsc.edu/goldenPath/hg17/encode/database/encodeGisChipPet.txt.gz">http://hgdownload.cse.ucsc.edu/goldenPath/hg17/encode/database/encodeGisChipPet.txt.gz</ext-link><x>, </x><comment>Available at <ext-link ext-link-type="uri" xlink:href="http://hgdownload.cse.ucsc.edu/downloads.html">http://hgdownload.cse.ucsc.edu/downloads.html</ext-link>.</comment></related-object></p></sec></sec><ref-list><title>References</title><ref id="bib1"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Anders</surname><given-names>S</given-names></name><name><surname>Huber</surname><given-names>W</given-names></name></person-group><year>2010</year><article-title>Differential expression analysis for sequence count data</article-title><source>Genome Biology</source><volume>11</volume><fpage>R106</fpage><pub-id pub-id-type="doi">10.1186/gb-2010-11-10-r106</pub-id></element-citation></ref><ref id="bib1a"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Andrysik</surname><given-names>A</given-names></name><name><surname>Kim</surname><given-names>J</given-names></name><name><surname>Tan</surname><given-names>AC</given-names></name><name><surname>Espinosa</surname><given-names>JM</given-names></name></person-group><year>2013</year><article-title>A genetic screen identifies TCF3/E2A and TRIAP1 as pathway-specific regulators of the cellular response to p53 activation</article-title><source>Cell Reports</source><volume>3</volume><fpage>1</fpage><lpage>9</lpage><pub-id pub-id-type="doi">10.1016/j.celrep.2013.04.014</pub-id></element-citation></ref><ref id="bib2"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Benson</surname><given-names>EK</given-names></name><name><surname>Mungamuri</surname><given-names>SK</given-names></name><name><surname>Attie</surname><given-names>O</given-names></name><name><surname>Kracikova</surname><given-names>M</given-names></name><name><surname>Sachidanandam</surname><given-names>R</given-names></name><name><surname>Manfredi</surname><given-names>JJ</given-names></name><name><surname>Aaronson</surname><given-names>SA</given-names></name></person-group><year>2013</year><article-title>p53-dependent gene repression through p21 is mediated by recruitment of E2F4 repression complexes</article-title><source>Oncogene</source><pub-id pub-id-type="doi">10.1038/onc.2013.378</pub-id></element-citation></ref><ref id="bib3"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Biderman</surname><given-names>L</given-names></name><name><surname>Manley</surname><given-names>JL</given-names></name><name><surname>Prives</surname><given-names>C</given-names></name></person-group><year>2012</year><article-title>Mdm2 and MdmX as regulators of gene expression</article-title><source>Genes & cancer</source><volume>3</volume><fpage>264</fpage><lpage>273</lpage><pub-id pub-id-type="doi">10.1177/1947601912455331</pub-id></element-citation></ref><ref id="bib4"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Brown</surname><given-names>CJ</given-names></name><name><surname>Lain</surname><given-names>S</given-names></name><name><surname>Verma</surname><given-names>CS</given-names></name><name><surname>Fersht</surname><given-names>AR</given-names></name><name><surname>Lane</surname><given-names>DP</given-names></name></person-group><year>2009</year><article-title>Awakening guardian angels: drugging the p53 pathway</article-title><source>Nature Reviews Cancer</source><volume>9</volume><fpage>862</fpage><lpage>873</lpage><pub-id pub-id-type="doi">10.1038/nrc2763</pub-id></element-citation></ref><ref id="bib5"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname><given-names>X</given-names></name><name><surname>Ko</surname><given-names>LJ</given-names></name><name><surname>Jayaraman</surname><given-names>L</given-names></name><name><surname>Prives</surname><given-names>C</given-names></name></person-group><year>1996</year><article-title>p53 levels, functional domains, and DNA damage determine the extent of the apoptotic response of tumor cells</article-title><source>Genes & Development</source><volume>10</volume><fpage>2438</fpage><lpage>2451</lpage><pub-id pub-id-type="doi">10.1101/gad.10.19.2438</pub-id></element-citation></ref><ref id="bib6"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Cohen</surname><given-names>PA</given-names></name><name><surname>Mani</surname><given-names>JC</given-names></name><name><surname>Lane</surname><given-names>DP</given-names></name></person-group><year>1998</year><article-title>Characterization of a new intrabody directed against the N-terminal region of human p53</article-title><source>Oncogene</source><volume>17</volume><fpage>2445</fpage><lpage>2456</lpage><pub-id pub-id-type="doi">10.1038/sj.onc.1202190</pub-id></element-citation></ref><ref id="bib7"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Core</surname><given-names>LJ</given-names></name><name><surname>Waterfall</surname><given-names>JJ</given-names></name><name><surname>Lis</surname><given-names>JT</given-names></name></person-group><year>2008</year><article-title>Nascent RNA sequencing reveals widespread pausing and divergent initiation at human promoters</article-title><source>Science</source><volume>322</volume><fpage>1845</fpage><lpage>1848</lpage><pub-id pub-id-type="doi">10.1126/science.1162228</pub-id></element-citation></ref><ref id="bib8"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Das</surname><given-names>S</given-names></name><name><surname>Raj</surname><given-names>L</given-names></name><name><surname>Zhao</surname><given-names>B</given-names></name><name><surname>Kimura</surname><given-names>Y</given-names></name><name><surname>Bernstein</surname><given-names>A</given-names></name><name><surname>Aaronson</surname><given-names>SA</given-names></name><name><surname>Lee</surname><given-names>SW</given-names></name></person-group><year>2007</year><article-title>Hzf Determines cell survival upon genotoxic stress by modulating p53 transactivation</article-title><source>Cell</source><volume>130</volume><fpage>624</fpage><lpage>637</lpage><pub-id pub-id-type="doi">10.1016/j.cell.2007.06.013</pub-id></element-citation></ref><ref id="bib9"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>el-Deiry</surname><given-names>WS</given-names></name><name><surname>Tokino</surname><given-names>T</given-names></name><name><surname>Velculescu</surname><given-names>VE</given-names></name><name><surname>Levy</surname><given-names>DB</given-names></name><name><surname>Parsons</surname><given-names>R</given-names></name><name><surname>Trent</surname><given-names>JM</given-names></name><name><surname>Lin</surname><given-names>D</given-names></name><name><surname>Mercer</surname><given-names>WE</given-names></name><name><surname>Kinzler</surname><given-names>KW</given-names></name><name><surname>Vogelstein</surname><given-names>B</given-names></name></person-group><year>1993</year><article-title>WAF1, a potential mediator of p53 tumor suppression</article-title><source>Cell</source><volume>75</volume><fpage>817</fpage><lpage>825</lpage><pub-id pub-id-type="doi">10.1016/0092-8674(93)90500-P</pub-id></element-citation></ref><ref id="bib10"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Espinosa</surname><given-names>JM</given-names></name><name><surname>Emerson</surname><given-names>BM</given-names></name></person-group><year>2001</year><article-title>Transcriptional regulation by p53 through intrinsic DNA/chromatin binding and site-directed cofactor recruitment</article-title><source>Molecular Cell</source><volume>8</volume><fpage>57</fpage><lpage>69</lpage><pub-id pub-id-type="doi">10.1016/S1097-2765(01)00283-0</pub-id></element-citation></ref><ref id="bib11"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Espinosa</surname><given-names>JM</given-names></name><name><surname>Verdun</surname><given-names>RE</given-names></name><name><surname>Emerson</surname><given-names>BM</given-names></name></person-group><year>2003</year><article-title>p53 functions through stress- and promoter-specific recruitment of transcription initiation components before and after DNA damage</article-title><source>Molecular Cell</source><volume>12</volume><fpage>1015</fpage><lpage>1027</lpage><pub-id pub-id-type="doi">10.1016/S1097-2765(03)00359-9</pub-id></element-citation></ref><ref id="bib12"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Garnett</surname><given-names>MJ</given-names></name><name><surname>Edelman</surname><given-names>EJ</given-names></name><name><surname>Heidorn</surname><given-names>SJ</given-names></name><name><surname>Greenman</surname><given-names>CD</given-names></name><name><surname>Dastur</surname><given-names>A</given-names></name><name><surname>Lau</surname><given-names>KW</given-names></name><name><surname>Greninger</surname><given-names>P</given-names></name><name><surname>Thompson</surname><given-names>IR</given-names></name><name><surname>Luo</surname><given-names>X</given-names></name><name><surname>Soares</surname><given-names>J</given-names></name><name><surname>Liu</surname><given-names>Q</given-names></name><name><surname>Iorio</surname><given-names>F</given-names></name><name><surname>Surdez</surname><given-names>D</given-names></name><name><surname>Chen</surname><given-names>L</given-names></name><name><surname>Milano</surname><given-names>RJ</given-names></name><name><surname>Bignell</surname><given-names>GR</given-names></name><name><surname>Tam</surname><given-names>AT</given-names></name><name><surname>Davies</surname><given-names>H</given-names></name><name><surname>Stevenson</surname><given-names>JA</given-names></name><name><surname>Barthorpe</surname><given-names>S</given-names></name><name><surname>Lutz</surname><given-names>SR</given-names></name><name><surname>Kogera</surname><given-names>F</given-names></name><name><surname>Lawrence</surname><given-names>K</given-names></name><name><surname>McLaren-Douglas</surname><given-names>A</given-names></name><name><surname>Mitropoulos</surname><given-names>X</given-names></name><name><surname>Mironenko</surname><given-names>T</given-names></name><name><surname>Thi</surname><given-names>H</given-names></name><name><surname>Richardson</surname><given-names>L</given-names></name><name><surname>Zhou</surname><given-names>W</given-names></name><name><surname>Jewitt</surname><given-names>F</given-names></name><name><surname>Zhang</surname><given-names>T</given-names></name><name><surname>O'Brien</surname><given-names>P</given-names></name><name><surname>Boisvert</surname><given-names>JL</given-names></name><name><surname>Price</surname><given-names>S</given-names></name><name><surname>Hur</surname><given-names>W</given-names></name><name><surname>Yang</surname><given-names>W</given-names></name><name><surname>Deng</surname><given-names>X</given-names></name><name><surname>Butler</surname><given-names>A</given-names></name><name><surname>Choi</surname><given-names>HG</given-names></name><name><surname>Chang</surname><given-names>JW</given-names></name><name><surname>Baselga</surname><given-names>J</given-names></name><name><surname>Stamenkovic</surname><given-names>I</given-names></name><name><surname>Engelman</surname><given-names>JA</given-names></name><name><surname>Sharma</surname><given-names>SV</given-names></name><name><surname>Delattre</surname><given-names>O</given-names></name><name><surname>Saez-Rodriguez</surname><given-names>J</given-names></name><name><surname>Gray</surname><given-names>NS</given-names></name><name><surname>Settleman</surname><given-names>J</given-names></name><name><surname>Futreal</surname><given-names>PA</given-names></name><name><surname>Haber</surname><given-names>DA</given-names></name><name><surname>Stratton</surname><given-names>MR</given-names></name><name><surname>Ramaswamy</surname><given-names>S</given-names></name><name><surname>McDermott</surname><given-names>U</given-names></name><name><surname>Benes</surname><given-names>CH</given-names></name></person-group><year>2012</year><article-title>Systematic identification of genomic markers of drug sensitivity in cancer cells</article-title><source>Nature</source><volume>483</volume><fpage>570</fpage><lpage>575</lpage><pub-id pub-id-type="doi">10.1038/nature11005</pub-id></element-citation></ref><ref id="bib13"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Gomes</surname><given-names>NP</given-names></name><name><surname>Bjerke</surname><given-names>G</given-names></name><name><surname>Llorente</surname><given-names>B</given-names></name><name><surname>Szostek</surname><given-names>SA</given-names></name><name><surname>Emerson</surname><given-names>BM</given-names></name><name><surname>Espinosa</surname><given-names>JM</given-names></name></person-group><year>2006</year><article-title>Gene-specific requirement for P-TEFb activity and RNA polymerase II phosphorylation within the p53 transcriptional program</article-title><source>Genes & Development</source><volume>20</volume><fpage>601</fpage><lpage>612</lpage><pub-id pub-id-type="doi">10.1101/gad.1398206</pub-id></element-citation></ref><ref id="bib14"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hah</surname><given-names>N</given-names></name><name><surname>Danko</surname><given-names>CG</given-names></name><name><surname>Core</surname><given-names>L</given-names></name><name><surname>Waterfall</surname><given-names>JJ</given-names></name><name><surname>Siepel</surname><given-names>A</given-names></name><name><surname>Lis</surname><given-names>JT</given-names></name><name><surname>Kraus</surname><given-names>WL</given-names></name></person-group><year>2011</year><article-title>A rapid, extensive, and transient transcriptional response to estrogen signaling in breast cancer cells</article-title><source>Cell</source><volume>145</volume><fpage>622</fpage><lpage>634</lpage><pub-id pub-id-type="doi">10.1016/j.cell.2011.03.042</pub-id></element-citation></ref><ref id="bib15"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Hargreaves</surname><given-names>DC</given-names></name><name><surname>Horng</surname><given-names>T</given-names></name><name><surname>Medzhitov</surname><given-names>R</given-names></name></person-group><year>2009</year><article-title>Control of inducible gene expression by signal-dependent transcriptional elongation</article-title><source>Cell</source><volume>138</volume><fpage>129</fpage><lpage>145</lpage><pub-id pub-id-type="doi">10.1016/j.cell.2009.05.047</pub-id></element-citation></ref><ref id="bib16"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Henry</surname><given-names>RE</given-names></name><name><surname>Andrysik</surname><given-names>Z</given-names></name><name><surname>Paris</surname><given-names>R</given-names></name><name><surname>Galbraith</surname><given-names>MD</given-names></name><name><surname>Espinosa</surname><given-names>JM</given-names></name></person-group><year>2012</year><article-title>A DR4:tBID axis drives the p53 apoptotic response by promoting oligomerization of poised BAX</article-title><source>The EMBO Journal</source><volume>31</volume><fpage>1266</fpage><lpage>1278</lpage><pub-id pub-id-type="doi">10.1038/emboj.2011.498</pub-id></element-citation></ref><ref id="bib17"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ivshina</surname><given-names>AV</given-names></name><name><surname>George</surname><given-names>J</given-names></name><name><surname>Senko</surname><given-names>O</given-names></name><name><surname>Mow</surname><given-names>B</given-names></name><name><surname>Putti</surname><given-names>TC</given-names></name><name><surname>Smeds</surname><given-names>J</given-names></name><name><surname>Lindahl</surname><given-names>T</given-names></name><name><surname>Pawitan</surname><given-names>Y</given-names></name><name><surname>Hall</surname><given-names>P</given-names></name><name><surname>Nordgren</surname><given-names>H</given-names></name><name><surname>Wong</surname><given-names>JE</given-names></name><name><surname>Liu</surname><given-names>ET</given-names></name><name><surname>Bergh</surname><given-names>J</given-names></name><name><surname>Kuznetsov</surname><given-names>VA</given-names></name><name><surname>Miller</surname><given-names>LD</given-names></name></person-group><year>2006</year><article-title>Genetic reclassification of histologic grade delineates new clinical subtypes of breast cancer</article-title><source>Cancer Research</source><volume>66</volume><fpage>10292</fpage><lpage>10301</lpage><pub-id pub-id-type="doi">10.1158/0008-5472.CAN-05-4414</pub-id></element-citation></ref><ref id="bib18"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jiao</surname><given-names>AL</given-names></name><name><surname>Slack</surname><given-names>FJ</given-names></name></person-group><year>2013</year><article-title>RNA-mediated gene activation</article-title><source>Epigenetics</source><volume>9</volume><fpage>27</fpage><lpage>36</lpage><pub-id pub-id-type="doi">10.4161/epi.26942</pub-id></element-citation></ref><ref id="bib19"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jin</surname><given-names>F</given-names></name><name><surname>Li</surname><given-names>Y</given-names></name><name><surname>Dixon</surname><given-names>JR</given-names></name><name><surname>Selvaraj</surname><given-names>S</given-names></name><name><surname>Ye</surname><given-names>Z</given-names></name><name><surname>Lee</surname><given-names>AY</given-names></name><name><surname>Yen</surname><given-names>CA</given-names></name><name><surname>Schmitt</surname><given-names>AD</given-names></name><name><surname>Espinoza</surname><given-names>CA</given-names></name><name><surname>Ren</surname><given-names>B</given-names></name></person-group><year>2013</year><article-title>A high-resolution map of the three-dimensional chromatin interactome in human cells</article-title><source>Nature</source><volume>503</volume><fpage>290</fpage><lpage>294</lpage><pub-id pub-id-type="doi">10.1038/nature12644</pub-id></element-citation></ref><ref id="bib20"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Kubbutat</surname><given-names>MH</given-names></name><name><surname>Jones</surname><given-names>SN</given-names></name><name><surname>Vousden</surname><given-names>KH</given-names></name></person-group><year>1997</year><article-title>Regulation of p53 stability by Mdm2</article-title><source>Nature</source><volume>387</volume><fpage>299</fpage><lpage>303</lpage><pub-id pub-id-type="doi">10.1038/387299a0</pub-id></element-citation></ref><ref id="bib21"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lal</surname><given-names>A</given-names></name><name><surname>Thomas</surname><given-names>MP</given-names></name><name><surname>Altschuler</surname><given-names>G</given-names></name><name><surname>Navarro</surname><given-names>F</given-names></name><name><surname>O'Day</surname><given-names>E</given-names></name><name><surname>Li</surname><given-names>XL</given-names></name><name><surname>Concepcion</surname><given-names>C</given-names></name><name><surname>Han</surname><given-names>YC</given-names></name><name><surname>Thiery</surname><given-names>J</given-names></name><name><surname>Rajani</surname><given-names>DK</given-names></name><name><surname>Deutsch</surname><given-names>A</given-names></name><name><surname>Hofmann</surname><given-names>O</given-names></name><name><surname>Ventura</surname><given-names>A</given-names></name><name><surname>Hide</surname><given-names>W</given-names></name><name><surname>Lieberman</surname><given-names>J</given-names></name></person-group><year>2011</year><article-title>Capture of microRNA-bound mRNAs identifies the tumor suppressor miR-34a as a regulator of growth factor signaling</article-title><source>PLOS Genetics</source><volume>7</volume><fpage>e1002363</fpage><pub-id pub-id-type="doi">10.1371/journal.pgen.1002363</pub-id></element-citation></ref><ref id="bib22"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Laptenko</surname><given-names>O</given-names></name><name><surname>Prives</surname><given-names>C</given-names></name></person-group><year>2006</year><article-title>Transcriptional regulation by p53: one protein, many possibilities</article-title><source>Cell Death and Differentiation</source><volume>13</volume><fpage>951</fpage><lpage>961</lpage><pub-id pub-id-type="doi">10.1038/sj.cdd.4401916</pub-id></element-citation></ref><ref id="bib23"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Li</surname><given-names>M</given-names></name><name><surname>He</surname><given-names>Y</given-names></name><name><surname>Dubois</surname><given-names>W</given-names></name><name><surname>Wu</surname><given-names>X</given-names></name><name><surname>Shi</surname><given-names>J</given-names></name><name><surname>Huang</surname><given-names>J</given-names></name></person-group><year>2012</year><article-title>Distinct regulatory mechanisms and functions for p53-activated and p53-repressed DNA damage response genes in embryonic stem cells</article-title><source>Molecular Cell</source><volume>46</volume><fpage>30</fpage><lpage>42</lpage><pub-id pub-id-type="doi">10.1016/j.molcel.2012.01.020</pub-id></element-citation></ref><ref id="bib24"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lowe</surname><given-names>SW</given-names></name><name><surname>Bodis</surname><given-names>S</given-names></name><name><surname>McClatchey</surname><given-names>A</given-names></name><name><surname>Remington</surname><given-names>L</given-names></name><name><surname>Ruley</surname><given-names>HE</given-names></name><name><surname>Fisher</surname><given-names>DE</given-names></name><name><surname>Housman</surname><given-names>DE</given-names></name><name><surname>Jacks</surname><given-names>T</given-names></name></person-group><year>1994</year><article-title>p53 status and the efficacy of cancer therapy in vivo</article-title><source>Science</source><volume>266</volume><fpage>807</fpage><lpage>810</lpage><pub-id pub-id-type="doi">10.1126/science.7973635</pub-id></element-citation></ref><ref id="bib25"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Melo</surname><given-names>CA</given-names></name><name><surname>Drost</surname><given-names>J</given-names></name><name><surname>Wijchers</surname><given-names>PJ</given-names></name><name><surname>van de Werken</surname><given-names>H</given-names></name><name><surname>de Wit</surname><given-names>E</given-names></name><name><surname>Oude Vrielink</surname><given-names>JA</given-names></name><name><surname>Elkon</surname><given-names>R</given-names></name><name><surname>Melo</surname><given-names>SA</given-names></name><name><surname>Leveille</surname><given-names>N</given-names></name><name><surname>Kalluri</surname><given-names>R</given-names></name><name><surname>de Laat</surname><given-names>W</given-names></name><name><surname>Agami</surname><given-names>R</given-names></name></person-group><year>2013</year><article-title>eRNAs are required for p53-dependent enhancer activity and gene transcription</article-title><source>Molecular Cell</source><volume>49</volume><fpage>524</fpage><lpage>535</lpage><pub-id pub-id-type="doi">10.1016/j.molcel.2012.11.021</pub-id></element-citation></ref><ref id="bib26"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Menendez</surname><given-names>D</given-names></name><name><surname>Nguyen</surname><given-names>TA</given-names></name><name><surname>Freudenberg</surname><given-names>JM</given-names></name><name><surname>Mathew</surname><given-names>VJ</given-names></name><name><surname>Anderson</surname><given-names>CW</given-names></name><name><surname>Jothi</surname><given-names>R</given-names></name><name><surname>Resnick</surname><given-names>MA</given-names></name></person-group><year>2013</year><article-title>Diverse stresses dramatically alter genome-wide p53 binding and transactivation landscape in human cancer cells</article-title><source>Nucleic Acids Research</source><volume>41</volume><fpage>7286</fpage><lpage>7301</lpage><pub-id pub-id-type="doi">10.1093/nar/gkt504</pub-id></element-citation></ref><ref id="bib27"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Momand</surname><given-names>J</given-names></name><name><surname>Zambetti</surname><given-names>GP</given-names></name><name><surname>Olson</surname><given-names>DC</given-names></name><name><surname>George</surname><given-names>D</given-names></name><name><surname>Levine</surname><given-names>AJ</given-names></name></person-group><year>1992</year><article-title>The mdm-2 oncogene product forms a complex with the p53 protein and inhibits p53-mediated transactivation</article-title><source>Cell</source><volume>69</volume><fpage>1237</fpage><lpage>1245</lpage><pub-id pub-id-type="doi">10.1016/0092-8674(92)90644-R</pub-id></element-citation></ref><ref id="bib28"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Morachis</surname><given-names>JM</given-names></name><name><surname>Murawsky</surname><given-names>CM</given-names></name><name><surname>Emerson</surname><given-names>BM</given-names></name></person-group><year>2010</year><article-title>Regulation of the p53 transcriptional response by structurally diverse core promoters</article-title><source>Genes & Development</source><volume>24</volume><fpage>135</fpage><lpage>147</lpage><pub-id pub-id-type="doi">10.1101/gad.1856710</pub-id></element-citation></ref><ref id="bib29"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Nikulenkov</surname><given-names>F</given-names></name><name><surname>Spinnler</surname><given-names>C</given-names></name><name><surname>Li</surname><given-names>H</given-names></name><name><surname>Tonelli</surname><given-names>C</given-names></name><name><surname>Shi</surname><given-names>Y</given-names></name><name><surname>Turunen</surname><given-names>M</given-names></name><name><surname>Kivioja</surname><given-names>T</given-names></name><name><surname>Ignatiev</surname><given-names>I</given-names></name><name><surname>Kel</surname><given-names>A</given-names></name><name><surname>Taipale</surname><given-names>J</given-names></name><name><surname>Selivanova</surname><given-names>G</given-names></name></person-group><year>2012</year><article-title>Insights into p53 transcriptional function via genome-wide chromatin occupancy and gene expression analysis</article-title><source>Cell Death and Differentiation</source><volume>19</volume><fpage>1992</fpage><lpage>2002</lpage><pub-id pub-id-type="doi">10.1038/cdd.2012.89</pub-id></element-citation></ref><ref id="bib30"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ohkubo</surname><given-names>S</given-names></name><name><surname>Tanaka</surname><given-names>T</given-names></name><name><surname>Taya</surname><given-names>Y</given-names></name><name><surname>Kitazato</surname><given-names>K</given-names></name><name><surname>Prives</surname><given-names>C</given-names></name></person-group><year>2006</year><article-title>Excess HDM2 impacts cell cycle and apoptosis and has a selective effect on p53-dependent transcription</article-title><source>The Journal of Biological Chemistry</source><volume>281</volume><fpage>16943</fpage><lpage>16950</lpage><pub-id pub-id-type="doi">10.1074/jbc.M601388200</pub-id></element-citation></ref><ref id="bib31"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Oliner</surname><given-names>JD</given-names></name><name><surname>Pietenpol</surname><given-names>JA</given-names></name><name><surname>Thiagalingam</surname><given-names>S</given-names></name><name><surname>Gyuris</surname><given-names>J</given-names></name><name><surname>Kinzler</surname><given-names>KW</given-names></name><name><surname>Vogelstein</surname><given-names>B</given-names></name></person-group><year>1993</year><article-title>Oncoprotein MDM2 conceals the activation domain of tumour suppressor p53</article-title><source>Nature</source><volume>362</volume><fpage>857</fpage><lpage>860</lpage><pub-id pub-id-type="doi">10.1038/362857a0</pub-id></element-citation></ref><ref id="bib32"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Picksley</surname><given-names>SM</given-names></name><name><surname>Vojtesek</surname><given-names>B</given-names></name><name><surname>Sparks</surname><given-names>A</given-names></name><name><surname>Lane</surname><given-names>DP</given-names></name></person-group><year>1994</year><article-title>Immunochemical analysis of the interaction of p53 with MDM2;–fine mapping of the MDM2 binding site on p53 using synthetic peptides</article-title><source>Oncogene</source><volume>9</volume><fpage>2523</fpage><lpage>2529</lpage></element-citation></ref><ref id="bib33"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Polyak</surname><given-names>K</given-names></name><name><surname>Xia</surname><given-names>Y</given-names></name><name><surname>Zweier</surname><given-names>JL</given-names></name><name><surname>Kinzler</surname><given-names>KW</given-names></name><name><surname>Vogelstein</surname><given-names>B</given-names></name></person-group><year>1997</year><article-title>A model for p53-induced apoptosis</article-title><source>Nature</source><volume>389</volume><fpage>300</fpage><lpage>305</lpage><pub-id pub-id-type="doi">10.1038/38525</pub-id></element-citation></ref><ref id="bib34"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Poyurovsky</surname><given-names>MV</given-names></name><name><surname>Katz</surname><given-names>C</given-names></name><name><surname>Laptenko</surname><given-names>O</given-names></name><name><surname>Beckerman</surname><given-names>R</given-names></name><name><surname>Lokshin</surname><given-names>M</given-names></name><name><surname>Ahn</surname><given-names>J</given-names></name><name><surname>Byeon</surname><given-names>IJ</given-names></name><name><surname>Gabizon</surname><given-names>R</given-names></name><name><surname>Mattia</surname><given-names>M</given-names></name><name><surname>Zupnick</surname><given-names>A</given-names></name><name><surname>Brown</surname><given-names>LM</given-names></name><name><surname>Friedler</surname><given-names>A</given-names></name><name><surname>Prives</surname><given-names>C</given-names></name></person-group><year>2010</year><article-title>The C terminus of p53 binds the N-terminal domain of MDM2</article-title><source>Nature Structural & Molecular Biology</source><volume>17</volume><fpage>982</fpage><lpage>989</lpage><pub-id pub-id-type="doi">10.1038/nsmb.1872</pub-id></element-citation></ref><ref id="bib35"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Riemenschneider</surname><given-names>MJ</given-names></name><name><surname>Buschges</surname><given-names>R</given-names></name><name><surname>Wolter</surname><given-names>M</given-names></name><name><surname>Reifenberger</surname><given-names>J</given-names></name><name><surname>Bostrom</surname><given-names>J</given-names></name><name><surname>Kraus</surname><given-names>JA</given-names></name><name><surname>Schlegel</surname><given-names>U</given-names></name><name><surname>Reifenberger</surname><given-names>G</given-names></name></person-group><year>1999</year><article-title>Amplification and overexpression of the MDM4 (MDMX) gene from 1q32 in a subset of malignant gliomas without TP53 mutation or MDM2 amplification</article-title><source>Cancer Research</source><volume>59</volume><fpage>6091</fpage><lpage>6096</lpage></element-citation></ref><ref id="bib36"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Riley</surname><given-names>T</given-names></name><name><surname>Sontag</surname><given-names>E</given-names></name><name><surname>Chen</surname><given-names>P</given-names></name><name><surname>Levine</surname><given-names>A</given-names></name></person-group><year>2008</year><article-title>Transcriptional control of human p53-regulated genes</article-title><source>Nature Reviews Molecular Cell Biology</source><volume>9</volume><fpage>402</fpage><lpage>412</lpage><pub-id pub-id-type="doi">10.1038/nrm2395</pub-id></element-citation></ref><ref id="bib37"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Schlereth</surname><given-names>K</given-names></name><name><surname>Heyl</surname><given-names>C</given-names></name><name><surname>Krampitz</surname><given-names>AM</given-names></name><name><surname>Mernberger</surname><given-names>M</given-names></name><name><surname>Finkernagel</surname><given-names>F</given-names></name><name><surname>Scharfe</surname><given-names>M</given-names></name><name><surname>Jarek</surname><given-names>M</given-names></name><name><surname>Leich</surname><given-names>E</given-names></name><name><surname>Rosenwald</surname><given-names>A</given-names></name><name><surname>Stiewe</surname><given-names>T</given-names></name></person-group><year>2013</year><article-title>Characterization of the p53 cistrome–DNA binding cooperativity dissects p53's tumor suppressor functions</article-title><source>PLOS Genetics</source><volume>9</volume><fpage>e1003726</fpage><pub-id pub-id-type="doi">10.1371/journal.pgen.1003726</pub-id></element-citation></ref><ref id="bib37a"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Smeenk</surname><given-names>L</given-names></name><name><surname>van Heeringen</surname><given-names>SJ</given-names></name><name><surname>Koeppel</surname><given-names>M</given-names></name><name><surname>van Driel</surname><given-names>MA</given-names></name><name><surname>Bartels</surname><given-names>SJ</given-names></name><name><surname>Akkers</surname><given-names>RC</given-names></name><name><surname>Denissov</surname><given-names>S</given-names></name><name><surname>Stunnenberg</surname><given-names>HG</given-names></name><name><surname>Lohrum</surname><given-names>M</given-names></name></person-group><year>2008</year><article-title>Characterization of genome-wide p53-binding sites upon stress response</article-title><source>Nucleic Acids Research</source><volume>36</volume><fpage>3639</fpage><lpage>3654</lpage><pub-id pub-id-type="doi">10.1093/nar/gkn232</pub-id></element-citation></ref><ref id="bib38"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Smeenk</surname><given-names>L</given-names></name><name><surname>van Heeringen</surname><given-names>SJ</given-names></name><name><surname>Koeppel</surname><given-names>M</given-names></name><name><surname>Gilbert</surname><given-names>B</given-names></name><name><surname>Janssen-Megens</surname><given-names>E</given-names></name><name><surname>Stunnenberg</surname><given-names>HG</given-names></name><name><surname>Lohrum</surname><given-names>M</given-names></name></person-group><year>2011</year><article-title>Role of p53 serine 46 in p53 target gene regulation</article-title><source>PLOS ONE</source><volume>6</volume><fpage>e17574</fpage><pub-id pub-id-type="doi">10.1371/journal.pone.0017574</pub-id></element-citation></ref><ref id="bib39"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sullivan</surname><given-names>KD</given-names></name><name><surname>Padilla-Just</surname><given-names>N</given-names></name><name><surname>Henry</surname><given-names>RE</given-names></name><name><surname>Porter</surname><given-names>CC</given-names></name><name><surname>Kim</surname><given-names>J</given-names></name><name><surname>Tentler</surname><given-names>JJ</given-names></name><name><surname>Eckhardt</surname><given-names>SG</given-names></name><name><surname>Tan</surname><given-names>AC</given-names></name><name><surname>DeGregori</surname><given-names>J</given-names></name><name><surname>Espinosa</surname><given-names>JM</given-names></name></person-group><year>2012</year><article-title>ATM and MET kinases are synthetic lethal with nongenotoxic activation of p53</article-title><source>Nature Chemical Biology</source><volume>8</volume><fpage>646</fpage><lpage>654</lpage><pub-id pub-id-type="doi">10.1038/nchembio.965</pub-id></element-citation></ref><ref id="bib40"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Szak</surname><given-names>ST</given-names></name><name><surname>Mays</surname><given-names>D</given-names></name><name><surname>Pietenpol</surname><given-names>JA</given-names></name></person-group><year>2001</year><article-title>Kinetics of p53 binding to promoter sites in vivo</article-title><source>Molecular and Cellular Biology</source><volume>21</volume><fpage>3375</fpage><lpage>3386</lpage><pub-id pub-id-type="doi">10.1128/MCB.21.10.3375-3386.2001</pub-id></element-citation></ref><ref id="bib41"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tang</surname><given-names>HY</given-names></name><name><surname>Zhao</surname><given-names>K</given-names></name><name><surname>Pizzolato</surname><given-names>JF</given-names></name><name><surname>Fonarev</surname><given-names>M</given-names></name><name><surname>Langer</surname><given-names>JC</given-names></name><name><surname>Manfredi</surname><given-names>JJ</given-names></name></person-group><year>1998</year><article-title>Constitutive expression of the cyclin-dependent kinase inhibitor p21 is transcriptionally regulated by the tumor suppressor protein p53</article-title><source>The Journal of Biological Chemistry</source><volume>273</volume><fpage>29156</fpage><lpage>29163</lpage><pub-id pub-id-type="doi">10.1074/jbc.273.44.29156</pub-id></element-citation></ref><ref id="bib42"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Thut</surname><given-names>CJ</given-names></name><name><surname>Goodrich</surname><given-names>JA</given-names></name><name><surname>Tjian</surname><given-names>R</given-names></name></person-group><year>1997</year><article-title>Repression of p53-mediated transcription by MDM2: a dual mechanism</article-title><source>Genes & Development</source><volume>11</volume><fpage>1974</fpage><lpage>1986</lpage><pub-id pub-id-type="doi">10.1101/gad.11.15.1974</pub-id></element-citation></ref><ref id="bib43"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Tovar</surname><given-names>C</given-names></name><name><surname>Rosinski</surname><given-names>J</given-names></name><name><surname>Filipovic</surname><given-names>Z</given-names></name><name><surname>Higgins</surname><given-names>B</given-names></name><name><surname>Kolinsky</surname><given-names>K</given-names></name><name><surname>Hilton</surname><given-names>H</given-names></name><name><surname>Zhao</surname><given-names>X</given-names></name><name><surname>Vu</surname><given-names>BT</given-names></name><name><surname>Qing</surname><given-names>W</given-names></name><name><surname>Packman</surname><given-names>K</given-names></name><name><surname>Myklebost</surname><given-names>O</given-names></name><name><surname>Heimbrook</surname><given-names>DC</given-names></name><name><surname>Vassilev</surname><given-names>LT</given-names></name></person-group><year>2006</year><article-title>Small-molecule MDM2 antagonists reveal aberrant p53 signaling in cancer: implications for therapy</article-title><source>Proceedings of the National Academy of Sciences of the United States of America</source><volume>103</volume><fpage>1888</fpage><lpage>1893</lpage><pub-id pub-id-type="doi">10.1073/pnas.0507493103</pub-id></element-citation></ref><ref id="bib44"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Vassilev</surname><given-names>LT</given-names></name><name><surname>Vu</surname><given-names>BT</given-names></name><name><surname>Graves</surname><given-names>B</given-names></name><name><surname>Carvajal</surname><given-names>D</given-names></name><name><surname>Podlaski</surname><given-names>F</given-names></name><name><surname>Filipovic</surname><given-names>Z</given-names></name><name><surname>Kong</surname><given-names>N</given-names></name><name><surname>Kammlott</surname><given-names>U</given-names></name><name><surname>Lukacs</surname><given-names>C</given-names></name><name><surname>Klein</surname><given-names>C</given-names></name><name><surname>Fotouhi</surname><given-names>N</given-names></name><name><surname>Liu</surname><given-names>EA</given-names></name></person-group><year>2004</year><article-title>In vivo activation of the p53 pathway by small-molecule antagonists of MDM2</article-title><source>Science</source><volume>303</volume><fpage>844</fpage><lpage>848</lpage><pub-id pub-id-type="doi">10.1126/science.1092472</pub-id></element-citation></ref><ref id="bib45"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Vousden</surname><given-names>KH</given-names></name><name><surname>Prives</surname><given-names>C</given-names></name></person-group><year>2009</year><article-title>Blinded by the light: the growing complexity of p53</article-title><source>Cell</source><volume>137</volume><fpage>413</fpage><lpage>431</lpage><pub-id pub-id-type="doi">10.1016/j.cell.2009.04.037</pub-id></element-citation></ref><ref id="bib46"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname><given-names>B</given-names></name><name><surname>Niu</surname><given-names>D</given-names></name><name><surname>Lam</surname><given-names>TH</given-names></name><name><surname>Xiao</surname><given-names>Z</given-names></name><name><surname>Ren</surname><given-names>EC</given-names></name></person-group><year>2013</year><article-title>Mapping the p53 transcriptome universe using p53 natural polymorphs</article-title><source>Cell Death and Differentiation</source><volume>21</volume><fpage>521</fpage><lpage>532</lpage><pub-id pub-id-type="doi">10.1038/cdd.2013.132</pub-id></element-citation></ref><ref id="bib47"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wei</surname><given-names>CL</given-names></name><name><surname>Wu</surname><given-names>Q</given-names></name><name><surname>Vega</surname><given-names>VB</given-names></name><name><surname>Chiu</surname><given-names>KP</given-names></name><name><surname>Ng</surname><given-names>P</given-names></name><name><surname>Zhang</surname><given-names>T</given-names></name><name><surname>Shahab</surname><given-names>A</given-names></name><name><surname>Yong</surname><given-names>HC</given-names></name><name><surname>Fu</surname><given-names>Y</given-names></name><name><surname>Weng</surname><given-names>Z</given-names></name><name><surname>Liu</surname><given-names>J</given-names></name><name><surname>Zhao</surname><given-names>XD</given-names></name><name><surname>Chew</surname><given-names>JL</given-names></name><name><surname>Lee</surname><given-names>YL</given-names></name><name><surname>Kuznetsov</surname><given-names>VA</given-names></name><name><surname>Sung</surname><given-names>WK</given-names></name><name><surname>Miller</surname><given-names>LD</given-names></name><name><surname>Lim</surname><given-names>B</given-names></name><name><surname>Liu</surname><given-names>ET</given-names></name><name><surname>Yu</surname><given-names>Q</given-names></name><name><surname>Ng</surname><given-names>HH</given-names></name><name><surname>Ruan</surname><given-names>Y</given-names></name></person-group><year>2006</year><article-title>A global map of p53 transcription-factor binding sites in the human genome</article-title><source>Cell</source><volume>124</volume><fpage>207</fpage><lpage>219</lpage><pub-id pub-id-type="doi">10.1016/j.cell.2005.10.043</pub-id></element-citation></ref><ref id="bib48"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>White</surname><given-names>DE</given-names></name><name><surname>Talbott</surname><given-names>KE</given-names></name><name><surname>Arva</surname><given-names>NC</given-names></name><name><surname>Bargonetti</surname><given-names>J</given-names></name></person-group><year>2006</year><article-title>Mouse double minute 2 associates with chromatin in the presence of p53 and is released to facilitate activation of transcription</article-title><source>Cancer Research</source><volume>66</volume><fpage>3463</fpage><lpage>3470</lpage><pub-id pub-id-type="doi">10.1158/0008-5472.CAN-05-1381</pub-id></element-citation></ref><ref id="bib49"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yan</surname><given-names>J</given-names></name><name><surname>Jiang</surname><given-names>J</given-names></name><name><surname>Lim</surname><given-names>CA</given-names></name><name><surname>Wu</surname><given-names>Q</given-names></name><name><surname>Ng</surname><given-names>HH</given-names></name><name><surname>Chin</surname><given-names>KC</given-names></name></person-group><year>2007</year><article-title>BLIMP1 regulates cell growth through repression of p53 transcription</article-title><source>Proceedings of the National Academy of Sciences of the United States of America</source><volume>104</volume><fpage>1841</fpage><lpage>1846</lpage><pub-id pub-id-type="doi">10.1073/pnas.0605562104</pub-id></element-citation></ref><ref id="bib50"><element-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhao</surname><given-names>R</given-names></name><name><surname>Gish</surname><given-names>K</given-names></name><name><surname>Murphy</surname><given-names>M</given-names></name><name><surname>Yin</surname><given-names>Y</given-names></name><name><surname>Notterman</surname><given-names>D</given-names></name><name><surname>Hoffman</surname><given-names>WH</given-names></name><name><surname>Tom</surname><given-names>E</given-names></name><name><surname>Mack</surname><given-names>DH</given-names></name><name><surname>Levine</surname><given-names>AJ</given-names></name></person-group><year>2000</year><article-title>Analysis of p53-regulated gene expression patterns using oligonucleotide arrays</article-title><source>Genes & Development</source><volume>14</volume><fpage>981</fpage><lpage>993</lpage><pub-id pub-id-type="doi">10.1101/gad.827700</pub-id></element-citation></ref></ref-list></back><sub-article article-type="article-commentary" id="SA1"><front-stub><article-id pub-id-type="doi">10.7554/eLife.02200.018</article-id><title-group><article-title>Decision letter</article-title></title-group><contrib-group content-type="section"><contrib contrib-type="editor"><name><surname>Prives</surname><given-names>Carol</given-names></name><role>Reviewing editor</role><aff><institution>Columbia University</institution>, <country>United States</country></aff></contrib></contrib-group></front-stub><body><boxed-text><p>eLife posts the editorial decision letter and author response on a selection of the published articles (subject to the approval of the authors). An edited version of the letter sent to the authors after peer review is shown, indicating the substantive concerns or comments; minor concerns are not usually shown. Reviewers have the opportunity to discuss the decision before the letter is sent (see <ext-link ext-link-type="uri" xlink:href="http://elifesciences.org/review-process">review process</ext-link>). Similarly, the author response typically shows only responses to the major concerns raised by the reviewers.</p></boxed-text><p>Thank you for sending your work entitled “Global analysis of p53-regulated transcription identifies its direct targets and unexpected regulatory mechanisms” for consideration at <italic>eLife</italic>. Your article has been favorably evaluated by a Senior editor and 3 reviewers, one of whom is a member of our Board of Reviewing Editors.</p><p>The Reviewing editor and the other reviewers discussed their comments before we reached this decision, and the Reviewing editor has assembled the following comments to help you prepare a revised submission.</p><p>The three referees agree that this is an interesting study that could be appropriate in principle for publication in <italic>eLife</italic>. It was felt, however, that although the study has some very valuable components, some more indepth analysis is needed. The two most interesting findings to address with more information and validation relate to the observation that Mdm2 is constitutively repressing p53 activation of some of its targets, and that GRO-seq analysis reveals several new targets of p53. The authors are requested to respond to the following comments.</p><p>Major comments:</p><p>1) It is possible that Nutlin does increase p53 levels in a small percent of the cells by 1 hour and it is those few cells with high p53 that have the increased transcription of p53 targets. Can the authors rule this out? Have they performed immunofluorescence analysis to see whether there are some ’jackpot‘ cells where p53 protein is rapidly induced by Nutlin?</p><p>2) While the Nutlin treatment allows direct ability to assay p53 function, a second method of p53 induction would greatly enhance these findings. While not direct necessarily, does an oncogenic or DNA-damage insult induce some of these new target genes (even if a few were assayed by qPCR on a gene by gene basis)? This would also address the survival vs. apoptotic regulation, since Nutlin treatment may have bias towards survival.</p><p>3) The authors should perform an Mdm2 ChIP at a few select targets such as p21 to see whether the Mdm2 signal is down by 1 hour. This is an important experiment because it would extend and confirm the most interesting aspect of this study. Endogenously expressed Mdm2 ChIPs are not very commonly reported in the literature, but have been accomplished. See for example, papers by J. Bargonetti and colleagues (White et al., Cancer Research 2006; Arva et al., JBC 2005).</p><p>4) Relatedly, in the Results section, the authors conclude that p53 acts as a repressor at a subset of genes under basal conditions (also shown in supplemental <xref ref-type="fig" rid="fig3">figure 3B</xref>). This is based on the comparison of expression in the p53 wt vs. null cells. Even though these genes are considered 'direct' targets from GRO-seq analysis after nutlin treatment, it is unclear whether they have evidence of p53 binding this small subset of genes, and since this is not an immediate knockdown experiment, should the authors consider secondary effects as they did in the first part of the paper? To support a direct role for p53 in repressing these promoters, the authors could ask whether p53 and/or MDM2 are at these basally regulated promoters by performing ChIP-qPCR or by using published ChIP-seq data.</p><p>5) It would be very interesting to choose several genes of interest (out of the genes represented in the figures) and check the Pol II CTD status across these genes. This will reveal that status of transcription and could possibly provide an explanation of why the mature mRNA does not correlate with the GRO-seq observations. This should also be compared to the Mdm2 ChIP(s) requested above to see if genes to which Mdm2 binds cause the Pol II CTD to have different phosphorylation patterns.</p><p>6) The expression of the newly identified genes, while some were of low abundance, were not confirmed at the protein level. Is there a functional readout of these RNAs?</p><p>7) The authors identify the primary response genes upon p53 activation. Is there a better way to leverage the valuable list of directly regulated genes to infer anything about the p53 activation/repression network? For instance, since they find that p53 does not directly down regulate many genes, can they see evidence of upregulation of known or putative repressors that might be responsible for down-regulation of genes identified as such in microarray studies?</p></body></sub-article><sub-article article-type="reply" id="SA2"><front-stub><article-id pub-id-type="doi">10.7554/eLife.02200.019</article-id><title-group><article-title>Author response</article-title></title-group></front-stub><body><p><italic>1) It is possible that Nutlin does increase p53 levels in a small percent of the cells by 1 hour and it is those few cells with high p53 that have the increased transcription of p53 targets. Can the authors rule this out? Have they performed immunofluorescence analysis to see whether there are some ’jackpot’ cells where p53 protein</italic> <italic>is rapidly induced by Nutlin?</italic></p><p>We are grateful for this comment, which inspired us to include new experiments, including the requested immunofluorescence, investigating the mechanism of action of Nutlin at the very short time point used for GRO-seq. First of all, the immunofluorescence discards the notion of ‘jackpot’ cells: although ∼3% cells show high p53 staining at the 1 hour time point of Nutlin treatment, this number is not significantly different from that observed in control p53 +/+ cells (<xref ref-type="fig" rid="fig1s1">Figure 1–figure supplement 1G–H</xref>).</p><p>The greater question behind this comment is: how could we observe such drastic transactivation of p53 target genes at 1 hour of Nutlin treatment without any noticeable increase in total p53 levels? Although it is well established that MDM2 represses p53 both by masking its transactivation domain and by targeting it for degradation, the contribution of each mechanism to the regulation of p53 target genes has been very difficult to tease out experimentally. This is not a trivial issue, because previous studies employing steady state mRNA measurements concluded that prolonged p53 activation and higher p53 levels were differentially required for activation of apoptotic genes (<xref ref-type="bibr" rid="bib5">Chen et al., 1996</xref>; <xref ref-type="bibr" rid="bib8">Das et al., 2007</xref>; <xref ref-type="bibr" rid="bib11">Espinosa et al., 2003</xref>; <xref ref-type="bibr" rid="bib40">Szak et al., 2001</xref>; <xref ref-type="bibr" rid="bib50">Zhao et al., 2000</xref>). However, GRO-seq demonstrates that a 1 hour time point of Nutlin treatment induces transcription of genes in every major pathway downstream of p53 (<xref ref-type="fig" rid="fig1">Figure 1 G</xref>). The observation that key survival (e.g. p21) and apoptotic genes (e.g. PIG3) show >6-fold increase in transcription at a time point preceding a proportional increase in total p53 levels as measured by Western blot, indicates that the mere unmasking of the p53 transactivation domain suffices to activate a multifaceted transcriptional program.</p><p>To further test this notion, we performed flow cytometry analyses using a monoclonal antibody (DO-1) that recognizes an epitope in the p53 N-terminal transactivation domain 1 (TAD1) that overlaps with the MDM2-binding surface competed by Nutlin (<xref ref-type="bibr" rid="bib32">Picksley et al., 1994</xref>). In fact, DO-1 is able to compete the p53-MDM2 interaction in vitro (<xref ref-type="bibr" rid="bib6">Cohen et al., 1998</xref>), in analogous fashion to Nutlin. Under the denaturing conditions of a Western Blot assay, where the p53-MDM2 complexes are fully disrupted, this antibody reveals no significant increase in total p53 levels at 1 hour Nutlin treatment (<xref ref-type="fig" rid="fig1">Figure 1C</xref>). However, we posited that DO-1 reactivity should increase significantly upon Nutlin treatment under the fixed conditions employed in flow cytometry. Flow cytometry quantitation shows that, even before Nutlin treatment, p53 +/+ cells have significantly more DO-1 reactivity than p53 -/- cells (<xref ref-type="fig" rid="fig1">Figure 1F</xref>). The functional importance of this ‘basal p53 activity’ is investigated at length later in our report (<xref ref-type="fig" rid="fig3">Figure 3</xref>). Interestingly, the p53 +/+ cell population shifts to significantly higher DO-1 reactivity at the 1 hour time point, as expected from epitope unmasking. A further increase is observed at 12 hours of Nutlin treatment, when total p53 levels have risen considerably as measured by Western blots (<xref ref-type="fig" rid="fig1">Figure 1C,G</xref>).</p><p>Altogether, these results indicate that the low levels of p53 present in proliferating cancer cells suffice to directly activate a multifunctional transcriptional program, including many canonical apoptotic genes, upon unmasking of the p53 transactivation domain by Nutlin. However, as discussed later in the paper (<xref ref-type="fig" rid="fig4">Figure 4</xref>), this conclusion does not necessarily conflict with previous reports showing differential timing of mRNA induction between arrest and apoptotic genes as seen by steady state RNA measurements.</p><p><italic>2) While the Nutlin treatment allows direct ability to assay p53 function, a second method of p53 induction would greatly enhance these findings. While not direct necessarily, does an oncogenic or DNA-damage insult induce some of these new target genes? (even if a few were assayed by qPCR on a gene by gene basis). This would also address the survival vs. apoptotic regulation, since Nutlin treatment may have bias towards survival</italic>.</p><p>We are grateful for this comment and the revised manuscript includes Q-RT-PCR for 14 of the novel genes discovered by GRO-seq showing that 12 of them are also induced by doxorubicin, a well established genotoxic p53 activating agent. This is now shown in <xref ref-type="fig" rid="fig1s1">Figure 1–figure supplement 1E</xref>.</p><p><italic>3) The authors should perform an Mdm2 ChIP at a few select targets such as p21 to see whether the Mdm2 signal is down by one hour. This is an important experiment because it would extend and confirm the most interesting aspect of this study. Endogenously expressed Mdm2 ChIPs are not very commonly reported in the literature, but have been accomplished. See for example papers by J. Bargonetti and colleagues (White et al., Cancer Research 2006; Arva et al., JBC 2005)</italic>.</p><p><italic>4) Relatedly, in the Results section, the authors conclude that p53 acts as a repressor at a subset of genes under basal conditions (also shown in supplemental</italic> <xref ref-type="fig" rid="fig3"><italic>figure 3B</italic></xref><italic>). This is based on the comparison of expression in the p53 wt vs. null cells. Even though these genes are considered 'direct' targets from GRO-seq analysis after nutlin treatment, it is unclear whether they have evidence of p53 binding this small subset of genes, and since this is not an immediate knockdown experiment, shouldn't the authors consider secondary effects as they did in the first part of the paper. To support a direct role for p53 in repressing these promoters, the authors could ask whether p53 and/or MDM2 are at these basally regulated promoters by performing ChIP-qPCR or by using published ChIP-seq data</italic>.</p><p>We thank the reviewers for these comment and we now include MDM2 ChIP data and p53 ChIP data for CDKN1A/p21 (basally activated) and two p53 target genes that are ‘basally repressed’ by low levels of p53 (<italic>PTP4A1</italic> and <italic>HES2</italic>) (<xref ref-type="fig" rid="fig3">Figure 3E</xref> and <xref ref-type="fig" rid="fig3s1">Figure 3–figure supplement 1B</xref>). The MDM2 ChIPs were technically challenging in HCT116 cells as MDM2 does not contact DNA directly and we could not confidently define MDM2 chromatin binding above background levels. However, we successfully managed to obtain MDM2 ChIP signals that are significantly enriched over IgG controls by employing the SJSA cell line, which expresses high MDM2 levels due to MDM2 gene amplification. We also made slight modifications to our ChIP protocol washing conditions. These technical modifications are explained in the revised Methods section.</p><p>These results are now explained in the revised text as follows:</p><p>“Overall, these results indicate that p53 acts as a repressor at a subset of its targets in a manner that is relieved by Nutlin, suggesting that MDM2 recruitment by basal levels of p53 may repress transcription at specific loci. To test this hypothesis, we performed ChIP experiments for p53 and MDM2 under conditions matching the GRO-seq experiment. […] Of note, although Nutlin disrupts the interaction between the p53 N-terminus and the hydrophobic pocket in the N-terminal domain of MDM2, a second molecular interaction occurs between the p53 C-terminus and the MDM2 N-terminus that is not competed by Nutlin in vitro (<xref ref-type="bibr" rid="bib34">Poyurovsky et al., 2010</xref>), which may explain why the MDM2 signal is not completely abrogated upon a short time point of Nutlin treatment. Western blots demonstrating specific MDM2 immunoprecipitation under the ChIP conditions utilized are shown in <xref ref-type="fig" rid="fig3s1">Figure 3–figure supplement 1C</xref>.”</p><p><italic>5) It would be very interesting to choose several genes of interest (out of the genes represented in the figures) and check the Pol II CTD status across these genes. This will reveal that status of transcription and could possibly provide an explanation of why the mature mRNA does not correlate with the GRO-seq observations. This should also be compared to the Mdm2 ChIP(s) requested above ip to see if genes to which Mdm2 binds cause the Pol II CTD to have different phosphorylation patterns</italic>.</p><p>As requested, we performed ChIP assays for Serine 5- and Serine 2-phosphorylated RNAPII matching our GRO-seq conditions. We now show ChIPs for <italic>GDF15</italic> (a ‘basally activated’ gene) and <italic>PTP4A1</italic> (a ‘basally repressed’ gene). Briefly, the phospho-RNAPII ChIPs reproduce the pattern of RNAPII behavior revealed by GRO-seq. The ‘basally repressed’ genes have lower RNAPII ChIP signals (and GRO-seq signals), when MDM2 is bound to the corresponding p53REs as shown by ChIP.</p><p>On a related note, we have now expanded our explanation of why GRO-seq and microarray measurements produce different pictures of p53 action, including a more detailed explanation of how microarrays fail to detect low abundance mRNAs.</p><p><italic>6. The expression of the newly identified genes, while some were of low abundance, were not confirmed at the protein level. Is there</italic> <italic>a functional readout of these RNAs?</italic></p><p>We now provide Western blot data showing induction at the protein level for three novel p53 target genes: GRIN2C, PTCHD4 and RINL. These data are in <xref ref-type="fig" rid="fig1s1">Figure 1–figure supplement 1F</xref>. Furthermore, we have expanded our explanation of why GRO-seq enabled us to discover so many novel genes not predicted by ChIP-seq/microarray. More specifically, we have included new p53 ChIP data for six of these novel genes (<xref ref-type="fig" rid="fig2s2">Figure 2–figure supplement 2</xref>), which shows that some of them have very distal p53 binding sites (i.e., >30kb from TSS), which would explain why they have been missed by ChIP-seq analysis requiring closer proximity to the TSS.</p><p><italic>7) The authors identify the primary response genes upon p53 activation. Is there a better way to leverage the valuable list of directly regulated genes to infer anything about the p53 activation/repression network? For instance, since they find that p53 does not directly down regulate many genes, can they see evidence of upregulation</italic> <italic>of known or putative repressors that might be responsible for downregulation of genes identified as such in microarray studies?</italic></p><p>As explained in our manuscript, Ingenuity Pathway Analysis of the ‘repressed’ genes identified E2F4, RB and CDKN1A as the top regulators of this gene set. Since CDKN1A (p21) is the most actively transcribed p53 target gene according to GRO-seq, its strong upregulation upon Nutlin treatment would lead to dephosphorylation of Rb and transcriptional repression by E2F4 complexes, thus explaining a large fraction of the ‘repression’ observed. In the manuscript we also show that another large fraction of the repression can be explained by upregulation of miR-34a, a p53-inducible miRNA detected by GRO-seq as a primary target. Analysis of our ‘repressed’ gene set led to the observation that up to 72% of them are downregulated in this same cell type upon miR-34a overexpression. Therefore, two direct targets or p53 transactivation (CDKN1A and miR-34a) can drive most of the indirect gene repression seen at a late time point by microarrays.</p><p>We searched for other transcriptional repressors within our list of immediate direct p53 targets and identified PRDM1 (also known as BLIMP1). PRDM1 is as a DNA binding protein involved in gene repression that was previously characterized as p53 target (<xref ref-type="bibr" rid="bib49">Yan et al., 2007</xref>). However, IPA analysis failed to identify PRDM1 as a top regulator of our ‘repressed’ gene set. In any case, it would be interesting to test in the future whether PRDM1 contributes to indirect gene repression downstream of p53. The revised manuscript mentions PRDM1 as another possible mediator of indirect gene repression downstream of p53.</p><p>Finally, our list of immediate p53 targets contains four other miRNAs beyond miR-34a (miR- 1204, miR-3189, miR4679-1 and miR-4692), which could also contribute to indirect post- transcriptional repression. Unfortunately, there is no experimental genomics data on the targets of these four miRNAs (only bioinformatic target predictions). The revised manuscript mentions these miRNAs and the possibility that they may also contribute to indirect gene repression downstream of p53.</p></body></sub-article></article> |