Diverse stresses dramatically alter genome-wide p53 binding and transactivation landscape in human cancer cells

Mechanisms of transcriptional regulation by p53

p53 is a transcription factor that suppresses tumor growth through regulation of dozens of target genes with diverse biological functions. The activity of this master transcription factor is inactivated in nearly all tumors, either by mutations in the TP53 locus or by oncogenic events that decrease the activity of the wild-type protein, such as overexpression of the p53 repressor MDM2. However, despite decades of intensive research, our collective understanding of the p53 signaling cascade remains incomplete. In this review, we focus on recent advances in our understanding of mechanisms of p53-dependent transcriptional control as they relate to five key areas: (1) the functionally distinct N-terminal transactivation domains, (2) the diverse regulatory roles of its C-terminal domain, (3) evidence that p53 is solely a direct transcriptional activator, not a direct repressor, (4) the ability of p53 to recognize many of its enhancers across diverse chromatin environments, and (5) mechanisms that modify the p53-dependent transcriptional program in a context-dependent manner.

p53 controls expression of the DNA deaminase APOBEC3B to limit its potential mutagenic activity in cancer cells

Cancer genome sequencing has implicated the cytosine deaminase activity of apolipoprotein B mRNA editing enzyme catalytic polypeptide-like (APOBEC) genes as an important source of mutations in diverse cancers, with APOBEC3B (A3B) expression especially correlated with such cancer mutations. To better understand the processes directing A3B over-expression in cancer, and possible therapeutic avenues for targeting A3B, we have investigated the regulation of A3B gene expression. Here, we show that A3B expression is inversely related to p53 status in different cancer types and demonstrate that this is due to a direct and pivotal role for p53 in repressing A3B expression. This occurs through the induction of p21 (CDKN1A) and the recruitment of the repressive DREAM complex to the A3B gene promoter, such that loss of p53 through mutation, or human papilloma virus-mediated inhibition, prevents recruitment of the complex, thereby causing elevated A3B expression and cytosine deaminase activity in cancer cells. As p53 is frequently mutated in cancer, our findings provide a mechanism by which p53 loss can promote cancer mutagenesis.

p53 regulates enhancer accessibility and activity in response to DNA damage

The tumor suppressor p53 is a well-characterized transcription factor that can bind gene promoters and regulate target gene transcription in response to DNA damage. Recent studies, however, have revealed that p53 binding events occur predominantly within regulatory enhancer elements. The effect of p53 binding on enhancer function has not been systematically evaluated. Here, we perform a genome-scale analysis of enhancer activity from p53-bound sequences using a series of massively parallel reporter assays (MPRAs) coupled with the assay for transposase-accessible chromatin (ATAC-Seq). We find that the majority of sequences examined display p53-dependent enhancer activity during the DNA damage response. Furthermore, we observe that p53 is bound to enhancer elements in healthy fibroblasts and poised for rapid activation in response to DNA damage. Surprisingly, our analyses revealed that most p53-bound enhancers are located within regions of inaccessible chromatin. A large subset of these enhancers become accessible following DNA damage indicating that p53 regulates their activity, in part, by modulating chromatin accessibility. The recognition and activation of enhancer elements located within inaccessible chromatin may contribute to the ability of the p53 network to function across the diverse chromatin landscapes of different tissues and cell types.

p53 binding sites in normal and cancer cells are characterized by distinct chromatin context

The tumor suppressor protein p53 interacts with DNA in a sequence-dependent manner. Thousands of p53 binding sites have been mapped genome-wide in normal and cancer cells. However, the way p53 selectively binds its cognate sites in different types of cells is not fully understood. Here, we performed a comprehensive analysis of 25 published p53 cistromes and identified 3,551 and 6,039 ‘high-confidence’ binding sites in normal and cancer cells, respectively. Our analysis revealed 2 distinct epigenetic features underlying p53-DNA interactions in vivo. First, p53 binding sites are associated with transcriptionally active histone marks (H3K4me3 and H3K36me3) in normal-cell chromatin, but with repressive histone marks (H3K27me3) in cancer-cell chromatin. Second, p53 binding sites in cancer cells are characterized by a lower level of DNA methylation than their counterparts in normal cells, probably related to global hypomethylation in cancers. Intriguingly, regardless of the cell type, p53 sites are highly enriched in the endogenous retroviral elements of the ERV1 family, highlighting the importance of this repeat family in shaping the transcriptional network of p53. Moreover, the p53 sites exhibit an unusual combination of chromatin patterns: high nucleosome occupancy and, at the same time, high sensitivity to DNase I. Our results suggest that p53 can access its target sites in a chromatin environment that is non-permissive to most DNA-binding transcription factors, which may allow p53 to act as a pioneer transcription factor in the context of chromatin.

Diverse p53/DNA binding modes expand the repertoire of p53 response elements

The tumor suppressor protein p53 acts as a transcription factor, binding sequence-specifically to defined DNA sites, thereby activating the expression of genes leading to diverse cellular outcomes. Canonical p53 response elements (REs) are made of two decameric half-sites separated by a variable number of base pairs (spacers). Fifty percent of all validated p53 REs contain spacers between 1 and 18 bp; however, their functional significance is unclear at present. Here, we show that p53 forms two different tetrameric complexes with consensus or natural REs, both with long spacers: a fully specific complex where two p53 dimers bind to two specific half-sites, and a hemispecific complex where one dimer binds to a specific half-site and the second binds to an adjacent spacer sequence. The two types of complexes have comparable binding affinity and specificity, as judged from binding competition against bulk genomic DNA. Structural analysis of the p53 REs in solution shows that these sites are not bent in both their free and p53-bound states when the two half-sites are either abutting or separated by spacers. Cell-based assay supports the physiological relevance of our findings. We propose that p53 REs with long spacers comprise separate specific half-sites that can lead to several different tetrameric complexes. This finding expands the universe of p53 binding sites and demonstrates that even isolated p53 half-sites can form tetrameric complexes. Moreover, it explains the manner in which p53 binds to clusters of more than one canonical binding site, common in many natural REs.

Identification of a core TP53 transcriptional program with highly distributed tumor suppressive activity

The tumor suppressor TP53 is the most frequently mutated gene product in human cancer. Close to half of all solid tumors carry inactivating mutations in the TP53 gene, while in the remaining cases, TP53 activity is abrogated by other oncogenic events, such as hyperactivation of its endogenous repressors MDM2 or MDM4. Despite identification of hundreds of genes regulated by this transcription factor, it remains unclear which direct target genes and downstream pathways are essential for the tumor suppressive function of TP53. We set out to address this problem by generating multiple genomic data sets for three different cancer cell lines, allowing the identification of distinct sets of TP53-regulated genes, from early transcriptional targets through to late targets controlled at the translational level. We found that although TP53 elicits vastly divergent signaling cascades across cell lines, it directly activates a core transcriptional program of ∼100 genes with diverse biological functions, regardless of cell type or cellular response to TP53 activation. This core program is associated with high-occupancy TP53 enhancers, high levels of paused RNA polymerases, and accessible chromatin. Interestingly, two different shRNA screens failed to identify a single TP53 target gene required for the anti-proliferative effects of TP53 during pharmacological activation in vitro. Furthermore, bioinformatics analysis of thousands of cancer genomes revealed that none of these core target genes are frequently inactivated in tumors expressing wild-type TP53. These results support the hypothesis that TP53 activates a genetically robust transcriptional program with highly distributed tumor suppressive functions acting in diverse cellular contexts.

Expression of SESN1, UHRF1BP1, and miR-377-3p as prognostic markers in mutated TP53 squamous cell carcinoma of the head and neck

The tumor suppressor gene TP53 is the most frequently mutated gene in human papillomavirus (HPV)-negative head and neck squamous cell carcinoma (HNSCC). It represents a known transcription factor that controls different microRNAs (miRNA) and target genes involved in the regulation of cellular stress, apoptosis and response to DNA damage. We used The Cancer Genome Atlas database to investigate the difference in transcriptome and proteome levels between mutated and wild-type TP53 HPV-negative HNSCC. Using different databases and an extensive literature review, we built the transcriptional and post-transcriptional network regulated by TP53. TP53 mutation was associated with poor overall survival in 203 HPV-negative patients compared to 40 patients with TP53 wild-type tumors. Using the enrichment analysis, we found that UHRF1BP1 and SESN1 mRNA were linked to prognosis in the TP53 mutated group. This is also the case for miR-377-3p, an important miRNA regulator of SESN1. Our study shows that SESN1 mRNA, UHRF1BP11 mRNA and miRNA-377-3p levels are prognostically relevant in HPV-negative HNSCC patients. This finding may help with patient stratification and the development of potential new therapeutic targets to treat patients with HNSCC.

Identification of a p53 target, CD137L, that mediates growth suppression and immune response of osteosarcoma cells

p53 encodes a transcription factor that transactivates downstream target genes involved in tumour suppression. Although osteosarcoma frequently has p53 mutations, the role of p53 in osteosarcomagenesis is not fully understood. To explore p53-target genes comprehensively in calvarial bone and find out novel druggable p53 target genes for osteosarcoma, we performed RNA sequencing using the calvarial bone and 23 other tissues from p53^+/+ and p53^−/− mice after radiation exposure. Of 23,813 genes, 69 genes were induced more than two-fold in irradiated p53^+/+ calvarial bone, and 127 genes were repressed. Pathway analysis of the p53-induced genes showed that genes associated with cytokine-cytokine receptor interactions were enriched. Three genes, CD137L, CDC42 binding protein kinase gamma and Follistatin, were identified as novel direct p53 target genes that exhibited growth-suppressive effects on osteosarcoma cell lines. Of the three genes, costimulatory molecule Cd137l was induced only in calvarial bone among the 24 tissues tested. CD137L-expressing cells exhibited growth-suppressive effects in vivo. In addition, recombinant Fc-fusion Cd137l protein activated the immune response in vitro and suppressed osteosarcoma cell growth in vivo. We clarified the role of CD137L in osteosarcomagenesis and its potential therapeutic application. Our transcriptome analysis also indicated the regulation of the immune response through p53.

Identification and Analysis of P53-Mediated Competing Endogenous RNA Network in Human Hepatocellular Carcinoma

Recent studies have indicated that long non-coding RNAs (lncRNAs) and mRNA function as competing endogenous RNAs (ceRNAs) that compete to bind to shared microRNA (miRNA) recognition elements (MREs) to perform specific biological functions during tumorigenesis. The tumor suppressor p53 is a master regulator of cancer-related biological processes by acting as a transcription factor to regulate target genes including miRNA and lncRNA. However, the mechanism in human hepatocellular carcinoma and whether p53-mediated RNA targets could form ceRNA network remain unclear. Here, we identified a series of differential expressed miRNAs, lncRNA and mRNA which were potentially regulated by p53 using RNA sequencing in HepG2. Genomic characteristics comparative analysis showed significant differences between mRNAs and lncRNAs. By integrating experimentally confirmed Ago2 and p53 binding sites, we constructed a highly reliable p53-mediated ceRNA network using hypergeometric test. The KEGG pathway enrichment analysis showed that the ceRNA network highly enriched in the cancer or p53-associated signaling pathways. Finally, using betweenness centrality analysis, we identified five master miRNAs (hsa-miR-3620-5p, hsa-miR-3613-3p, hsa-miR-6881-3p, hsa-miR-6087 and hsa-miR-18a-3p) that regulated most of the target RNAs, suggesting these miRNAs play central roles in the whole p53-mediated ceRNAs network. Taken together, our results provide a new regulatory mechanism of p53 networks for future studies in cancer therapeutics.

Effects of distinct drugs on gene transcription in an osteosarcoma cell line

Osteosarcoma (OS) is a common cancerous bone tumor which has a detrimental impact on the lives of patients and their families. The present study aimed at investigating the underlying molecular mechanism of various drug treatments pertaining to OS, including dimethyl sulfoxide (DMSO), doxorubicin (DXP), Nutlin-3, actinomycin D (ActD) and etoposide (Eto). Microarray and p53 chromatin immunoprecipitation combined with sequencing (ChIP-seq) datasets of the OS cell line U2OS treated with distinct drugs were acquired from the Gene Expression Omnibus and differentially-expressed genes (DEGs) were screened for alignment analysis. The p53-binding target genes were identified and ChIP-seq and microarray gene expression data were combined to identify directly and indirectly targeted genes. A regulatory network of p53 was constructed with the acquired data. Finally, the Database for Annotation, Visualization and Integrated Discovery was interrogated for annotation of target genes. A total of 212 p53-binding peaks were obtained in the untreated group, whereas thousands of peaks were obtained in the treated groups. In total, ~1,000 target genes were identified in each of DXP, DMSO, Eto and ActD treatment groups, whereas the Nutlin-3 treatment group identified an increased number, with 5,458 target genes obtained. Several common DEGs including MDM2, TP53I3, RRM2B, FAS and SESN1 were targeted by all the drugs with the exception of DMSO. p53 regulated various genes including EHF, HOXA10 and BHLHE40 in the Nutlin-3 treatment group, whereas p53 regulated EHF, RFX3, TRAF40 and TCF7L2 in the DXR treatment group. The results of the present study indicate that p53 was able to directly regulate target genes including MDM2, TP53I3 and RRM2B or indirectly regulate numerous further genes through several hub genes including EHF and RFX through various drug treatments in U2OS cells. Furthermore, p53 regulated distinct molecular processes in various drug treatments.

Assessment of the DNA damaging potential of environmental chemicals using a quantitative high-throughput screening approach to measure p53 activation

Genotoxicity potential is a critical component of any comprehensive toxicological profile. Compounds that induce DNA or chromosomal damage often activate p53, a transcription factor essential to cell cycle regulation. Thus, within the US Tox21 Program, we screened a library of ∼10,000 (∼8,300 unique) environmental compounds and drugs for activation of the p53-signaling pathway using a quantitative high-throughput screening assay employing HCT-116 cells (p53^+/+ ) containing a stably integrated β-lactamase reporter gene under control of the p53 response element (p53RE). Cells were exposed (-S9) for 16 hr at 15 concentrations (generally 1.2 nM to 92 μM) three times, independently. Excluding compounds that failed analytical chemistry analysis or were suspected of inducing assay interference, 365 (4.7%) of 7,849 unique compounds were concluded to activate p53. As part of an in-depth characterization of our results, we first compared them with results from traditional in vitro genotoxicity assays (bacterial mutation, chromosomal aberration); ∼15% of known, direct-acting genotoxicants in our library activated the p53RE. Mining the Comparative Toxicogenomics Database revealed that these p53 actives were significantly associated with increased expression of p53 downstream genes involved in DNA damage responses. Furthermore, 53 chemical substructures associated with genotoxicity were enriched in certain classes of p53 actives, for example, anthracyclines (antineoplastics) and vinca alkaloids (tubulin disruptors). Interestingly, the tubulin disruptors manifested unusual nonmonotonic concentration response curves suggesting activity through a unique p53 regulatory mechanism. Through the analysis of our results, we aim to define a role for this assay as one component of a comprehensive toxicological characterization of large compound libraries. Environ. Mol. Mutagen. 58:494-507, 2017. © 2017 Wiley Periodicals, Inc.

DDX54 regulates transcriptome dynamics during DNA damage response

The cellular response to genotoxic stress is mediated by a well-characterized network of DNA surveillance pathways. The contribution of post-transcriptional gene regulatory networks to the DNA damage response (DDR) has not been extensively studied. Here, we systematically identified RNA-binding proteins differentially interacting with polyadenylated transcripts upon exposure of human breast carcinoma cells to ionizing radiation (IR). Interestingly, more than 260 proteins, including many nucleolar proteins, showed increased binding to poly(A)+ RNA in IR-exposed cells. The functional analysis of DDX54, a candidate genotoxic stress responsive RNA helicase, revealed that this protein is an immediate-to-early DDR regulator required for the splicing efficacy of its target IR-induced pre-mRNAs. Upon IR exposure, DDX54 acts by increased interaction with a well-defined class of pre-mRNAs that harbor introns with weak acceptor splice sites, as well as by protein-protein contacts within components of U2 snRNP and spliceosomal B complex, resulting in lower intron retention and higher processing rates of its target transcripts. Because DDX54 promotes survival after exposure to IR, its expression and/or mutation rate may impact DDR-related pathologies. Our work indicates the relevance of many uncharacterized RBPs potentially involved in the DDR.

Census and evaluation of p53 target genes

The tumor suppressor p53 functions primarily as a transcription factor. Mutation of the TP53 gene alters its response pathway, and is central to the development of many cancers. The discovery of a large number of p53 target genes, which confer p53's tumor suppressor function, has led to increasingly complex models of p53 function. Recent meta-analysis approaches, however, are simplifying our understanding of how p53 functions as a transcription factor. In the survey presented here, a total set of 3661 direct p53 target genes is identified that comprise 3509 potential targets from 13 high-throughput studies, and 346 target genes from individual gene analyses. Comparison of the p53 target genes reported in individual studies with those identified in 13 high-throughput studies reveals limited consistency. Here, p53 target genes have been evaluated based on the meta-analysis data, and the results show that high-confidence p53 target genes are involved in multiple cellular responses, including cell cycle arrest, DNA repair, apoptosis, metabolism, autophagy, mRNA translation and feedback mechanisms. However, many p53 target genes are identified only in a small number of studies and have a higher likelihood of being false positives. While numerous mechanisms have been proposed for mediating gene regulation in response to p53, recent advances in our understanding of p53 function show that p53 itself is solely an activator of transcription, and gene downregulation by p53 is indirect and requires p21. Taking into account the function of p53 as an activator of transcription, recent results point to an unsophisticated means of regulation.

The p53 gene with emphasis on its paralogues in mosquitoes

The p53 gene is highly important in human cancers, as it serves as a tumor-suppressor gene. Subsequently, two p53 homologues, i.e., p73 and p63, with high identity of amino acids were identified, leading to construction of the p53 family. The p53 gene is highly important in human cancer because it usually transcribes genes that function by causing apoptosis in mammalian cells. In contrast, p63 and p73 tend to be more important in modulating development than inducing cell death, even though they share similar protein structures. Relatively recently, p53 was also identified in mosquitoes and many other insect species. Uniquely, its structure lacks the sterile alpha motif domain which is a putative protein-protein interaction domain and exclusively exists at the C-terminal region in p73 and p63 in mammals. A phylogenetic analysis revealed that the p53 gene derived from mosquitoes is composed of two paralogues, p53-1 and p53-2. Of these, only p53-2 is responsively upregulated by dengue 2 virus (DENV2) in C6/36 cells which usually survive the infection. This indicates that the p53 gene is closely related to DENV infection in mosquito cells. The specific significance of p53-2's involvement in cell survival from virus-induced stress is described and briefly discussed in this report.

Prosurvival long noncoding RNA PINCR regulates a subset of p53 targets in human colorectal cancer cells by binding to Matrin 3

Thousands of long noncoding RNAs (lncRNAs) have been discovered, yet the function of the vast majority remains unclear. Here, we show that a p53-regulated lncRNA which we named PINCR (p53-induced noncoding RNA), is induced ~100-fold after DNA damage and exerts a prosurvival function in human colorectal cancer cells (CRC) in vitro and tumor growth in vivo. Targeted deletion of PINCR in CRC cells significantly impaired G1 arrest and induced hypersensitivity to chemotherapeutic drugs. PINCR regulates the induction of a subset of p53 targets involved in G1 arrest and apoptosis, including BTG2, RRM2B and GPX1. Using a novel RNA pulldown approach that utilized endogenous S1-tagged PINCR, we show that PINCR associates with the enhancer region of these genes by binding to RNA-binding protein Matrin 3 that, in turn, associates with p53. Our findings uncover a critical prosurvival function of a p53/PINCR/Matrin 3 axis in response to DNA damage in CRC cells.

DOI:http://dx.doi.org/10.7554/eLife.23244.001

Though DNA contains the information needed to build the proteins that keep cells alive, only 2% of the DNA in a human cell codes for proteins. The remaining 98% is referred to as non-coding DNA. The information in some of these non-coding regions can still be copied into molecules of RNA, including long molecules called lncRNAs. Little is known about what lncRNAs actually do, but growing evidence suggests that these molecules are important for a number of vital processes including cell growth and survival.

When the DNA in an animal cell gets damaged, the cell needs to decide whether to pause growth and repair the damage, or to kill itself if the harm is too great. One of the best-studied proteins guiding this decision is the p53 protein, which increases the number of protein-coding genes needed to carry out either option in this decision. That is to say that, p53 regulates the genes needed to kill the cell and the genes needed to temporarily pause its growth and repair the damage, which instead keeps the cell alive. So, how does the p53 protein guide the decision, and are lncRNA molecules involved?

Using human colon cancer cells, Chaudhary et al. now report that when DNA is damaged, the levels of a specific lncRNA increase 100-fold. Further experiments showed that this lncRNA – named PINCR, which refers to p53-induced noncoding RNA – promotes the survival of cells. Chaudhary et al. showed that PINCR molecules do this by recruiting a protein called Matrin 3 to a certain region in the DNA called an enhancer and then links it to promoter region in the DNA of specific genes that temporarily pause cell growth but keep the cell alive. This in turn activates these ‘pro-survival genes’. In further experiments, when the PINCR molecules were essentially deleted, p53 was not able to fully activate these genes and as a result more of the cells died.

Together these findings increase our knowledge of how lncRNAs can work, especially in the context of DNA damage in cancer cells. A next important step will be to uncover other roles for the PINCR molecule in both cancer and healthy cells.

DOI:http://dx.doi.org/10.7554/eLife.23244.002

The Cytidine Deaminase APOBEC3 Family Is Subject to Transcriptional Regulation by p53

The APOBEC3 (A3) family of proteins are DNA cytidine deaminases that act as sentinels in the innate immune response against retroviral infections and are responsive to IFN. Recently, a few A3 genes were identified as potent enzymatic sources of mutations in several human cancers. Using human cancer cells and lymphocytes, we show that under stress conditions and immune challenges, all A3 genes are direct transcriptional targets of the tumor suppressor p53. Although the expression of most A3 genes (including A3C and A3H) was stimulated by the activation of p53, treatment with the DNA-damaging agent doxorubicin or the p53 stabilizer Nutlin led to repression of the A3B gene. Furthermore, p53 could enhance IFN type-I induction of A3 genes. Interestingly, overexpression of a group of tumor-associated p53 mutants in TP53-null cancer cells promoted A3B expression. These findings establish that the "guardian of the genome" role ascribed to p53 also extends to a unique component of the immune system, the A3 genes, thereby integrating human immune and chromosomal stress responses into an A3/p53 immune axis.Implications: Activated p53 can integrate chromosomal stresses and immune responses through its influence on expression of APOBEC3 genes, which are key components of the innate immune system that also influence genomic stability. Mol Cancer Res; 15(6); 735-44. ©2017 AACR.

RECQ1 expression is upregulated in response to DNA damage and in a p53-dependent manner

Sensitivity of cancer cells to DNA damaging chemotherapeutics is determined by DNA repair processes. Consequently, cancer cells may upregulate the expression of certain DNA repair genes as a mechanism to promote chemoresistance. Here, we report that RECQ1, a breast cancer susceptibility gene that encodes the most abundant RecQ helicase in humans, is a p53-regulated gene, potentially acting as a defense against DNA damaging agents. We show that RECQ1 mRNA and protein levels are upregulated upon treatment of cancer cells with a variety of DNA damaging agents including the DNA-alkylating agent methylmethanesulfonate (MMS). The MMS-induced upregulation of RECQ1 expression is p53-dependent as it was observed in p53-proficient but not in isogenic p53-deficient cells. The RECQ1 promoter is bound by endogenous p53 and is responsive to p53 in luciferase reporter assays suggesting that RECQ1 is a direct target of p53. Treatment with the chemotherapeutic drugs temozolomide and fotemustine also increased RECQ1 mRNA levels whereas depletion of RECQ1 enhanced cellular sensitivity to these agents. These results identify a previously unrecognized p53-mediated upregulation of RECQ1 expression in response to DNA damage and implicate RECQ1 in the repair of DNA lesions including those induced by alkylating and other chemotherapeutic agents.

MicroRNAs as Key Effectors in the p53 Network

The guardian of the genome p53 is embedded in a fine-spun network of MicroRNAs. p53 is able to activate or repress directly the transcription of MicroRNAs that are participating in the tumor-suppressive mission of p53. On the other hand, the expression of p53 is under tight control of MicroRNAs that are either targeting directly p53 or factors that are modifying its protein level or activity. Although the most important function of p53 is suggested to be transcriptional regulation, there are several nontranscriptional functions described. One of those regards the modulation of MicroRNA biogenesis. Wild-type p53 is increasing the maturation of selected MicroRNAs from the primary transcript to the precursor MiRNA by interacting with the Microprocessor complex. Furthermore, p53 is modulating the mRNA accessibility for certain MicroRNAs by association with the RISC complex and transcriptional regulation of RNA-binding proteins. In this way p53 is able to remodel the MiRNA-mRNA interaction network. As wild-type p53 is employing MicroRNAs to suppress cancer development, gain-of-function mutant p53 proteins use MicroRNAs to confer oncogenic properties like chemoresistance and the ability to drive metastasis. Like its wild-type counterpart mutant p53 is able to regulate MicroRNAs transcriptionally and posttranscriptionally. Mutant p53 affects the MiRNA processing at two cleavage steps through interfering with the Microprocessor complex and by downregulating Dicer and KSRP, a modulator of MiRNA biogenesis. Thus, MicroRNAs are essential components in the p53 pathway, contributing substantially to combat or enhance tumor development depending on the wild-type or mutant p53 context.

p53 Maintains Baseline Expression of Multiple Tumor Suppressor Genes

TP53 is the most commonly mutated tumor suppressor gene and its mutation drives tumorigenesis. Using ChIP-seq for p53 in the absence of acute cell stress, we found that wild-type but not mutant p53 binds and activates numerous tumor suppressor genes, including PTEN, STK11(LKB1), miR-34a, KDM6A(UTX), FOXO1, PHLDA3, and TNFRSF10B through consensus binding sites in enhancers and promoters. Depletion of p53 reduced expression of these target genes, and analysis across 18 tumor types showed that mutation of TP53 associated with reduced expression of many of these genes. Regarding PTEN, p53 activated expression of a luciferase reporter gene containing the p53-consensus site in the PTEN enhancer, and homozygous deletion of this region in cells decreased PTEN expression and increased growth and transformation. These findings show that p53 maintains expression of a team of tumor suppressor genes that may together with the stress-induced targets mediate the ability of p53 to suppress cancer development. p53 mutations selected during tumor initiation and progression, thus, inactivate multiple tumor suppressor genes in parallel, which could account for the high frequency of p53 mutations in cancer.Implications: In this study, we investigate the activities of p53 under normal low-stress conditions and discover that p53 is capable of maintaining the expression of a group of important tumor suppressor genes at baseline, many of which are haploinsufficient, which could contribute to p53-mediated tumor suppression. Mol Cancer Res; 15(8); 1051-62. ©2017 AACR.

Human p53 interacts with the elongating RNAPII complex and is required for the release of actinomycin D induced transcription blockage

The p53 tumour suppressor regulates the transcription initiation of selected genes by binding to specific DNA sequences at their promoters. Here we report a novel role of p53 in transcription elongation in human cells. Our data demonstrate that upon transcription elongation blockage, p53 is associated with genes that have not been reported as its direct targets. p53 could be co-immunoprecipitated with active forms of DNA-directed RNA polymerase II subunit 1 (RPB1), highlighting its association with the elongating RNA polymerase II. During a normal transcription cycle, p53 and RPB1 are localised at distinct regions of selected non-canonical p53 target genes and this pattern of localisation was changed upon blockage of transcription elongation. Additionally, transcription elongation blockage induced the proteasomal degradation of RPB1. Our results reveal a novel role of p53 in human cells during transcription elongation blockage that may facilitate the removal of RNA polymerase II from DNA.

Genome-wide analysis of p53-regulated transcription in Myc-driven lymphomas

The tumour suppressor p53 is a transcription factor that controls cellular stress responses. Here, we dissected the transcriptional programmes triggered upon restoration of p53 in Myc-driven lymphomas, based on the integrated analysis of p53 genomic occupancy and gene regulation. p53 binding sites were identified at promoters and enhancers, both characterized by the pre-existence of active chromatin marks. Only a small fraction of these sites showed the 20 base-pair p53 consensus motif, suggesting that p53 recruitment to genomic DNA was primarily mediated through protein-protein interactions in a chromatin context. p53 also targeted distal sites devoid of activation marks, at which binding was prevalently driven by sequence recognition. In all instances, the relevant motif was the canonical unsplit consensus element, with no clear evidence for p53 recruitment by split motifs. At promoters, p53 binding to the consensus motif was associated with gene induction, but not repression, indicating that the latter was most likely indirect. Altogether, our data highlight key features of genome recognition by p53 and provide unprecedented insight into the pathways associated with p53 reactivation and tumour regression, paving the way for their therapeutic application.

p53 Specifically Binds Triplex DNA In Vitro and in Cells

Triplex DNA is implicated in a wide range of biological activities, including regulation of gene expression and genomic instability leading to cancer. The tumor suppressor p53 is a central regulator of cell fate in response to different type of insults. Sequence and structure specific modes of DNA recognition are core attributes of the p53 protein. The focus of this work is the structure-specific binding of p53 to DNA containing triplex-forming sequences in vitro and in cells and the effect on p53-driven transcription. This is the first DNA binding study of full-length p53 and its deletion variants to both intermolecular and intramolecular T.A.T triplexes. We demonstrate that the interaction of p53 with intermolecular T.A.T triplex is comparable to the recognition of CTG-hairpin non-B DNA structure. Using deletion mutants we determined the C-terminal DNA binding domain of p53 to be crucial for triplex recognition. Furthermore, strong p53 recognition of intramolecular T.A.T triplexes (H-DNA), stabilized by negative superhelicity in plasmid DNA, was detected by competition and immunoprecipitation experiments, and visualized by AFM. Moreover, chromatin immunoprecipitation revealed p53 binding T.A.T forming sequence in vivo. Enhanced reporter transactivation by p53 on insertion of triplex forming sequence into plasmid with p53 consensus sequence was observed by luciferase reporter assays. In-silico scan of human regulatory regions for the simultaneous presence of both consensus sequence and T.A.T motifs identified a set of candidate p53 target genes and p53-dependent activation of several of them (ABCG5, ENOX1, INSR, MCC, NFAT5) was confirmed by RT-qPCR. Our results show that T.A.T triplex comprises a new class of p53 binding sites targeted by p53 in a DNA structure-dependent mode in vitro and in cells. The contribution of p53 DNA structure-dependent binding to the regulation of transcription is discussed.

The novel p53 target TNFAIP8 variant 2 is increased in cancer and offsets p53-dependent tumor suppression

Tumor necrosis factor-α-induced protein 8 (TNFAIP8) is a stress-response gene that has been associated with cancer, but no studies have differentiated among or defined the regulation or function of any of its several recently described expression variants. We found that TNFAIP8 variant 2 (v2) is overexpressed in multiple human cancers, whereas other variants are commonly downregulated in cancer (v1) or minimally expressed in cancer or normal tissue (v3-v6). Silencing v2 in cancer cells induces p53-independent inhibition of DNA synthesis, widespread binding of p53, and induction of target genes and p53-dependent cell cycle arrest and DNA damage sensitization. Cell cycle arrest induced by v2 silencing requires p53-dependent induction of p21. In response to the chemotherapeutic agent doxorubicin, p53 regulates v2 through binding to an intragenic enhancer, together indicating that p53 and v2 engage in complex reciprocal regulation. We propose that TNFAIP8 v2 promotes human cancer by broadly repressing p53 function, in essence offsetting p53-dependent tumor suppression.

The Tail That Wags the Dog: How the Disordered C-Terminal Domain Controls the Transcriptional Activities of the p53 Tumor-Suppressor Protein

The p53 tumor suppressor is a transcription factor (TF) that exerts antitumor functions through its ability to regulate the expression of multiple genes. Within the p53 protein resides a relatively short unstructured C-terminal domain (CTD) that remarkably participates in virtually every aspect of p53 performance as a TF. Because these aspects are often interdependent and it is not always possible to dissect them experimentally, there has been a great deal of controversy about the CTD. In this review we evaluate the significance and key features of this interesting region of p53 and its impact on the many aspects of p53 function in light of previous and more recent findings.

Identification of cellular senescence-specific genes by comparative transcriptomics

Cellular senescence is defined as permanent cell cycle arrest induced by various stresses. Although the p53 transcriptional activity is essential for senescence induction, the downstream genes that are crucial for senescence remain unsolved. Here, by using a developed experimental system in which cellular senescence or apoptosis is induced preferentially by altering concentration of etoposide, a DNA-damaging drug, we compared gene expression profiles of senescent and apoptotic cells by microarray analysis. Subtraction of the expression profile of apoptotic cells identified 20 genes upregulated specifically in senescent cells. Furthermore, 6 out of 20 genes showed p53-dependent upregulation by comparing gene expression between p53-proficient and -deficient cells. These 6 genes were also upregulated during replicative senescence of normal human diploid fibroblasts, suggesting that upregulation of these genes is a general phenomenon in senescence. Among these genes, 2 genes (PRODH and DAO) were found to be directly regulated by p53, and ectopic expression of 4 genes (PRODH, DAO, EPN3, and GPR172B) affected senescence phenotypes induced by etoposide treatment. Collectively, our results identified several proteins as novel downstream effectors of p53-mediated senescence and provided new clues for further research on the complex signalling networks underlying the induction and maintenance of senescence.

Integration of TP53, DREAM, MMB-FOXM1 and RB-E2F target gene analyses identifies cell cycle gene regulatory networks

Cell cycle (CC) and TP53 regulatory networks are frequently deregulated in cancer. While numerous genome-wide studies of TP53 and CC-regulated genes have been performed, significant variation between studies has made it difficult to assess regulation of any given gene of interest. To overcome the limitation of individual studies, we developed a meta-analysis approach to identify high confidence target genes that reflect their frequency of identification in independent datasets. Gene regulatory networks were generated by comparing differential expression of TP53 and CC-regulated genes with chromatin immunoprecipitation studies for TP53, RB1, E2F, DREAM, B-MYB, FOXM1 and MuvB. RNA-seq data from p21-null cells revealed that gene downregulation by TP53 generally requires p21 (CDKN1A). Genes downregulated by TP53 were also identified as CC genes bound by the DREAM complex. The transcription factors RB, E2F1 and E2F7 bind to a subset of DREAM target genes that function in G1/S of the CC while B-MYB, FOXM1 and MuvB control G2/M gene expression. Our approach yields high confidence ranked target gene maps for TP53, DREAM, MMB-FOXM1 and RB-E2F and enables prediction and distinction of CC regulation. A web-based atlas at www.targetgenereg.org enables assessing the regulation of any human gene of interest.

Genome-wide analysis of p53 transcriptional programs in B cells upon exposure to genotoxic stress in vivo

The tumor suppressor p53 is a transcription factor that coordinates the cellular response to DNA damage. Here we provide an integrated analysis of p53 genomic occupancy and p53-dependent gene regulation in the splenic B and non-B cell compartments of mice exposed to whole-body ionizing radiation, providing insight into general principles of p53 activity in vivo. In unstressed conditions, p53 bound few genomic targets; induction of p53 by ionizing radiation increased the number of p53 bound sites, leading to highly overlapping profiles in the different cell types. Comparison of these profiles with chromatin features in unstressed B cells revealed that, upon activation, p53 localized at active promoters, distal enhancers, and a smaller set of unmarked distal regions. At promoters, recognition of the canonical p53 motif as well as binding strength were associated with p53-dependent transcriptional activation, but not repression, indicating that the latter was most likely indirect. p53-activated targets constituted the core of a cell type-independent response, superimposed onto a cell type-specific program. Core response genes included most of the known p53-regulated genes, as well as many new ones. Our data represent a unique characterization of the p53-regulated response to ionizing radiation in vivo.

Developing a Gene Biomarker at the Tipping Point of Adaptive and Adverse Responses in Human Bronchial Epithelial Cells

Determining mechanism-based biomarkers that distinguish adaptive and adverse cellular processes is critical to understanding the health effects of environmental exposures. Shifting from in vivo, low-throughput toxicity studies to high-throughput screening (HTS) paradigms and risk assessment based on in vitro and in silico testing requires utilizing toxicity pathway information to distinguish adverse outcomes from recoverable adaptive events. Little work has focused on oxidative stresses in human airway for the purposes of predicting adverse responses. We hypothesize that early gene expression-mediated molecular changes could be used to delineate adaptive and adverse responses to environmentally-based perturbations. Here, we examined cellular responses of the tracheobronchial airway to zinc (Zn) exposure, a model oxidant. Airway derived BEAS-2B cells exposed to 2–10 μM Zn²⁺ elicited concentration- and time-dependent cytotoxicity. Normal, adaptive, and cytotoxic Zn²⁺ exposure conditions were determined with traditional apical endpoints, and differences in global gene expression around the tipping point of the responses were used to delineate underlying molecular mechanisms. Bioinformatic analyses of differentially expressed genes indicate early enrichment of stress signaling pathways, including those mediated by the transcription factors p53 and NRF2. After 4 h, 154 genes were differentially expressed (p < 0.01) between the adaptive and cytotoxic Zn²⁺ concentrations. Nearly 40% of the biomarker genes were related to the p53 signaling pathway with 30 genes identified as likely direct targets using a database of p53 ChIP-seq studies. Despite similar p53 activation profiles, these data revealed widespread dampening of p53 and NRF2-related genes as early as 4 h after exposure at higher, unrecoverable Zn²⁺ exposures. Thus, in our model early increased activation of stress response pathways indicated a recoverable adaptive event. Overall, this study highlights the importance of characterizing molecular mechanisms around the tipping point of adverse responses to better inform HTS paradigms.

Multiplex enhancer-reporter assays uncover unsophisticated TP53 enhancer logic

Transcription factors regulate their target genes by binding to regulatory regions in the genome. Although the binding preferences of TP53 are known, it remains unclear what distinguishes functional enhancers from nonfunctional binding. In addition, the genome is scattered with recognition sequences that remain unoccupied. Using two complementary techniques of multiplex enhancer-reporter assays, we discovered that functional enhancers could be discriminated from nonfunctional binding events by the occurrence of a single TP53 canonical motif. By combining machine learning with a meta-analysis of TP53 ChIP-seq data sets, we identified a core set of more than 1000 responsive enhancers in the human genome. This TP53 cistrome is invariably used between cell types and experimental conditions, whereas differences among experiments can be attributed to indirect nonfunctional binding events. Our data suggest that TP53 enhancers represent a class of unsophisticated cell-autonomous enhancers containing a single TP53 binding site, distinct from complex developmental enhancers that integrate signals from multiple transcription factors.

p53 amplifies Toll-like receptor 5 response in human primary and cancer cells through interaction with multiple signal transduction pathways

The p53 tumor suppressor regulates transcription of genes associated with diverse cellular functions including apoptosis, growth arrest, DNA repair and differentiation. Recently, we established that p53 can modulate expression of Toll-like receptor (TLR) innate immunity genes but the degree of cross-talk between p53 and TLR pathways remained unclear. Here, using gene expression profiling we characterize the global effect of p53 on the TLR5-mediated transcription in MCF7 cells. We found that combined activation of p53 and TLR5 pathways synergistically increases expression of over 200 genes, mostly associated with immunity and inflammation. The synergy was observed in several human cancer cells and primary lymphocytes. The p53-dependent amplification of transcriptional response to TLR5 activation required expression of NFκB subunit p65 and was mediated by several molecular mechanisms including increased phosphorylation of p38 MAP kinase, PI3K and STAT3 signaling. Additionally, p53 induction increased cytokine expression in response to TNFα, another activator of NFκB and MAP kinase pathways, suggesting a broad interaction between p53 and these signaling pathways. The expression of many synergistically induced genes is elevated in breast cancer patients responsive to chemotherapy. We suggest that p53's capacity to enhance immune response could be exploited to increase antitumor immunity and to improve cancer treatment.

Biphasic and Stage-Associated Expression of CPEB4 in Hepatocellular Carcinoma

Cytoplasmic polyadenylation element binding protein 4 (CPEB4) is a sequence-specific RNA-binding protein and translational regulator, with expression associated with tumor progression. Nevertheless, CPEB4 seems to play paradoxical roles in cancers–an oncogenic promoter in pancreatic ductal adenocarcinoma (PDA) and glioblastomas but a tumor suppressor in hepatocellular carcinoma (HCC). To assess whether CPEB4-regulated carcinogenesis is tissue-specific, we reevaluated the role of CPEB4 in HCC. Although proliferation of hepatocytes appeared normal in CPEB4-knockout (KO) mice after partial hepatectomy, knockdown (KD) of CPEB4 in HepG2 liver cancer cells promoted colony formation in vitro. Moreover, the growth of CPEB4-KD cells was accelerated in an in vivo xenograft mouse model. In tumorous and adjacent non-tumorous paired liver specimens from 49 HCC patients, the protein level of CPEB4 was significantly upregulated in early-stage HCC but decreased toward late-stage HCC. This finding agrees with changes in CPEB4 mRNA level from analysis of two sets of HCC microarray data from the Gene Expression Omnibus (GEO) database. Taken together, downregulation of CPEB4 likely occurs at the late cancer stage to facilitate HCC progression. Biphasic alteration of CPEB4 expression during HCC progression suggests its complicated role in tumorigenesis.

The common stress responsive transcription factor ATF3 binds genomic sites enriched with p300 and H3K27ac for transcriptional regulation

Background

Dysregulation of the common stress responsive transcription factor ATF3 has been causally linked to many important human diseases such as cancer, atherosclerosis, infections, and hypospadias. Although it is believed that the ATF3 transcription activity is central to its cellular functions, how ATF3 regulates gene expression remains largely unknown. Here, we employed ATF3 wild-type and knockout isogenic cell lines to carry out the first comprehensive analysis of global ATF3-binding profiles in the human genome under basal and stressed (DNA damage) conditions.

Results

Although expressed at a low basal level, ATF3 was found to bind a large number of genomic sites that are often associated with genes involved in cellular stress responses. Interestingly, ATF3 appears to bind a large portion of genomic sites distal to transcription start sites and enriched with p300 and H3K27ac. Global gene expression profiling analysis indicates that genes proximal to these genomic sites were often regulated by ATF3. While DNA damage elicited by camptothecin dramatically altered the ATF3 binding profile, most of the genes regulated by ATF3 upon DNA damage were pre-bound by ATF3 before the stress. Moreover, we demonstrated that ATF3 was co-localized with the major stress responder p53 at genomic sites, thereby collaborating with p53 to regulate p53 target gene expression upon DNA damage.

Conclusions

These results suggest that ATF3 likely bookmarks genomic sites and interacts with other transcription regulators to control gene expression.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2664-8) contains supplementary material, which is available to authorized users.

The p53 C terminus controls site-specific DNA binding and promotes structural changes within the central DNA binding domain

DNA binding by numerous transcription factors including the p53 tumor suppressor protein constitutes a vital early step in transcriptional activation. While the role of the central core DNA binding domain (DBD) of p53 in site-specific DNA binding has been established, the contribution of the sequence-independent C-terminal domain (CTD) is still not well understood. We investigated the DNA-binding properties of a series of p53 CTD variants using a combination of in vitro biochemical analyses and in vivo binding experiments. Our results provide several unanticipated and interconnected findings. First, the CTD enables DNA binding in a sequence-dependent manner that is drastically altered by either its modification or deletion. Second, dependence on the CTD correlates with the extent to which the p53 binding site deviates from the canonical consensus sequence. Third, the CTD enables stable formation of p53-DNA complexes to divergent binding sites via DNA-induced conformational changes within the DBD itself.

The RNA Response to DNA Damage

Multicellular organisms must ensure genome integrity to prevent accumulation of mutations, cell death, and cancer. The DNA damage response (DDR) is a complex network that senses, signals, and executes multiple programs including DNA repair, cell cycle arrest, senescence, and apoptosis. This entails regulation of a variety of cellular processes: DNA replication and transcription, RNA processing, mRNA translation and turnover, and post-translational modification, degradation, and relocalization of proteins. Accumulated evidence over the past decades has shown that RNAs and RNA metabolism are both regulators and regulated actors of the DDR. This review aims to present a comprehensive overview of the current knowledge on the many interactions between the DNA damage and RNA fields.

Transglutaminase 2 contributes to a TP53-induced autophagy program to prevent oncogenic transformation

Genetic alterations which impair the function of the TP53 signaling pathway in TP53 wild-type human tumors remain elusive. To identify new components of this pathway, we performed a screen for genes whose loss-of-function debilitated TP53 signaling and enabled oncogenic transformation of human mammary epithelial cells. We identified transglutaminase 2 (TGM2) as a putative tumor suppressor in the TP53 pathway. TGM2 suppressed colony formation in soft agar and tumor formation in a xenograft mouse model. The depletion of growth supplements induced both TGM2 expression and autophagy in a TP53-dependent manner, and TGM2 promoted autophagic flux by enhancing autophagic protein degradation and autolysosome clearance. Reduced expression of both CDKN1A, which regulates the cell cycle downstream of TP53, and TGM2 synergized to promote oncogenic transformation. Our findings suggest that TGM2-mediated autophagy and CDKN1A-mediated cell cycle arrest are two important barriers in the TP53 pathway that prevent oncogenic transformation.

DOI:http://dx.doi.org/10.7554/eLife.07101.001

Cancers grow from rogue cells that manage to defy the strict rules that normally stop a cell from dividing when it should not. Each cell contains many proteins that are responsible for implementing these rules, and thus help to prevent tumors from forming. One of these proteins – p53 (which is also called TP53) – plays a central role in this process. Information about many processes within and around a cell filters through the p53 protein, before being passed on to a range of different proteins.

The proteins that are alerted by p53 are commonly referred to as its 'downstream effectors', and it is these proteins that stop cells from dividing too much. For example, the protein p21 (also called CDKN1A) – which is the best understood of p53’s downstream effectors – hinders the machinery that causes cells to divide. Other p53 effectors can cause cells to kill themselves to prevent cancer growth. However, recent experiments with mice predicted that there may be other p53’s effectors that are important too.

Yeo, Itahana et al. have now depleted the proteins that potentially work in p53’s network, one by one, in human cells called mammary epithelial cells, to test if these cells can become cancerous in the laboratory. The experiments showed that another downstream effector protein of p53 – an enzyme called transglutaminase 2 – contributes to prevent these mammary epithelial cells from becoming cancerous. Transglutaminase 2 promotes a process known as autophagy, which recycles damaged and old components of the cell, and therefore normally helps to keep cells healthy.

Yeo, Itahana et al. also demonstrated that the effects of both p21 and transglutaminase 2 are critical to stop human mammary epithelial cells grown in the laboratory from dividing too much and from forming tumors when injected into mice.

These experiments provide a deeper understanding of how most cells manage to remain healthy rather than becoming cancerous and reveal a potential new target for the early detection of cancer. Further investigations could now explore whether therapies could re-activate this enzyme to prevent or treat cancer.

DOI:http://dx.doi.org/10.7554/eLife.07101.002

A/T gap tolerance in the core sequence and flanking sequence requirements of non-canonical p53 response elements

The canonical core sequence of the p53 response element, CATG, has a two-base A/T gap. Previously, we found that p53 can also activate a non-canonical four-base A/T gap CATATG core sequence. In this study, we investigated the possible number of A/T bases used by p53 and showed that a six-base A/T gap CATATATG core sequence was the maximum A/T gap in the p53 response element that could be upregulated by p53 and p63. Canonical and non-canonical p53 response elements also have three-base flanking sequences. A/T bases could be substituted by G/C bases, including CACACG and CGTGTG, but not CGCGCG. We found that the SV40 promoter with functional six- and two-base A/T gap core sequences could be activated by TAp63γ and that TAp63γ could upregulate SV40 small and large T antigens expression in COS7 cells. We also found that the distal region of PUMA promoter with functional two six-base A/T gap core sequences could be activated by TAp63γ in 293T cells. These new findings could provide novel rules for the non-canonical p53 family response element and could extend the entire p53 family regulation network.

The Establishment of a Hyperactive Structure Allows the Tumour Suppressor Protein p53 to Function through P-TEFb during Limited CDK9 Kinase Inhibition

CDK9 is the catalytic subunit of positive elongation factor b (P-TEFb) that controls the transition of RNA polymerase II (RNAPII) into elongation. CDK9 inhibitors block mRNA synthesis and trigger activation of the stress-sensitive p53 protein. This in turn induces transcription of CDKN1A (p21) and other cell cycle control genes. It is presently unclear if and how p53 circumvents a general P-TEFb-requirement when it activates its target genes. Our investigations using a panel of specific inhibitors reason for a critical role of CDK9 also in the case of direct inhibition of the kinase. At the prototypic p21 gene, the activator p53 initially accumulates at the pre-bound upstream enhancer followed—with significant delay—by de novo binding to a secondary enhancer site within the first intron of p21. This is accompanied by recruitment of the RNAPII initiation machinery to both elements. ChIP and functional analyses reason for a prominent role of CDK9 itself and elongation factor complexes PAF1c and SEC involved in pause and elongation control. It appears that the strong activation potential of p53 facilitates gene activation in the situation of global repression of RNAPII transcription. The data further underline the fundamental importance of CDK9 for class II gene transcription.

Indirect p53-dependent transcriptional repression of Survivin, CDC25C, and PLK1 genes requires the cyclin-dependent kinase inhibitor p21/CDKN1A and CDE/CHR promoter sites binding the DREAM complex

The transcription factor p53 is central to cell cycle control by downregulation of cell cycle-promoting genes upon cell stress such as DNA damage. Survivin (BIRC5), CDC25C, and PLK1 encode important cell cycle regulators that are repressed following p53 activation. Here, we provide evidence that p53-dependent repression of these genes requires activation of p21 (CDKN1A, WAF1, CIP1). Chromatin immunoprecipitation (ChIP) data indicate that promoter binding of B-MYB switches to binding of E2F4 and p130 resulting in a replacement of the MMB (Myb-MuvB) by the DREAM complex. We demonstrate that this replacement depends on p21. Furthermore, transcriptional repression by p53 requires intact DREAM binding sites in the target promoters. The CDE and CHR cell cycle promoter elements are the sites for DREAM binding. These elements as well as the p53 response of Survivin, CDC25C, and PLK1 are evolutionarily conserved. No binding of p53 to these genes is detected by ChIP and mutation of proposed p53 binding sites does not alter the p53 response. Thus, a mechanism for direct p53-dependent transcriptional repression is not supported by the data. In contrast, repression by DREAM is consistent with most previous findings and unifies models based on p21-, E2F4-, p130-, and CDE/CHR-dependent repression by p53. In conclusion, the presented data suggest that the p53-p21-DREAM-CDE/CHR pathway regulates p53-dependent repression of Survivin, CDC25C, and PLK1.

Sequence Variation in the Response Element Determines Binding by the Transcription Factor p73

How the sequence of a response element affects the binding of a transcription factor and, ultimately, the differential rate of transcription of genes under its control is not well-understood. In the case of the p73 transcription factor, it binds to >200 response elements to trigger developmental, cell arrest, and apoptotic pathways. The p73 response elements match the 20 bp consensus sequence of the p53 response elements that are formed by two 10 bp half-sites, where each half-site is an inverted repeat of two 5 bp quarter-sites. Using sedimentation velocity and fluorescence anisotropy experiments, we studied how systematic variations in the sequence of a half-site response element modify the DNA binding affinity of the p73 DNA-binding domain. We observed that each nucleotide position in the response element has a different influence in determining the binding of the p73 DNA-binding domain. The cytosine in the fourth position of each quarter-site is the largest determinant of DNA binding, followed by the nucleotide in the fifth position, and last, the first three positions show a slight regulatory preference for purines. Together with previous structural and functional results, our data suggest a hierarchical model of binding in which some nucleotide positions in the response element are more important than others in determining the binding of the transcription factor.

The Role of Structural Elements of the 5'-Terminal Region of p53 mRNA in Translation under Stress Conditions Assayed by the Antisense Oligonucleotide Approach

The p53 protein is one of the major factors responsible for cell cycle regulation and stress response. In the 5'-terminal region of p53 mRNA, an IRES element has been found which takes part in the translational regulation of p53 expression. Two characteristic hairpin motifs are present in this mRNA region: G56-C169, with the first AUG codon, and U180-A218, which interacts with the Hdm2 protein (human homolog of mouse double minute 2 protein). 2'-OMe modified antisense oligomers hybridizing to the 5'-terminal region of p53 mRNA were applied to assess the role of these structural elements in translation initiation under conditions of cellular stress. Structural changes in the RNA target occurring upon oligomers' binding were monitored by the Pb2+-induced cleavage method. The impact of antisense oligomers on the synthesis of two proteins, the full-length p53 and its isoform Δ40p53, was analysed in HT-29, MCF-7 and HepG2 cells, under normal conditions and under stress, as well as in vitro conditions. The results revealed that the hairpin U180-A218 and adjacent single-stranded region A219-A228 were predominantly responsible for high efficacy of IRES-mediated translation in the presence of stress factors. These motifs play a role of cis-acting elements which are able to modulate IRES activity, likely via interactions with protein factors.

Cooperative interactions between p53 and NFκB enhance cell plasticity

The p53 and NFκB sequence-specific transcription factors play crucial roles in cell proliferation and survival with critical, even if typically opposite, effects on cancer progression. To investigate a possible crosstalk between p53 and NFκB driven by chemotherapy-induced responses in the context of an inflammatory microenvironment, we performed a proof of concept study using MCF7 cells. Transcriptome analyses upon single or combined treatments with doxorubicin (Doxo, 1.5μM) and the NFκB inducer TNF-alpha (TNF⍺, 5ng/ml) revealed 432 up-regulated (log₂ FC> 2), and 390 repressed genes (log₂ FC< -2) for the Doxo+TNF⍺ treatment. 239 up-regulated and 161 repressed genes were synergistically regulated by the double treatment. Annotation and pathway analyses of Doxo+TNF⍺ selectively up-regulated genes indicated strong enrichment for cell migration terms. A panel of genes was examined by qPCR coupled to p53 activation by Doxo, 5-Fluoruracil and Nutlin-3a, or to p53 or NFκB inhibition. Transcriptome data were confirmed for 12 of 15 selected genes and seven (PLK3, LAMP3, ETV7, UNC5B, NTN1, DUSP5, SNAI1) were synergistically up-regulated after Doxo+TNF⍺ and dependent both on p53 and NFκB. Migration assays consistently showed an increase in motility for MCF7 cells upon Doxo+TNF⍺. A signature of 29 Doxo+TNF⍺ highly synergistic genes exhibited prognostic value for luminal breast cancer patients, with adverse outcome correlating with higher relative expression. We propose that the crosstalk between p53 and NFκB can lead to the activation of specific gene expression programs that may impact on cancer phenotypes and potentially modify the efficacy of cancer therapy.

The p53-p21-DREAM-CDE/CHR pathway regulates G2/M cell cycle genes

The tumor suppressor p53 functions predominantly as a transcription factor by activating and downregulating gene expression, leading to cell cycle arrest or apoptosis. p53 was shown to indirectly repress transcription of the CCNB2, KIF23 and PLK4 cell cycle genes through the recently discovered p53-p21-DREAM-CDE/CHR pathway. However, it remained unclear whether this pathway is commonly used. Here, we identify genes regulated by p53 through this pathway in a genome-wide computational approach. The bioinformatic analysis is based on genome-wide DREAM complex binding data, p53-depedent mRNA expression data and a genome-wide definition of phylogenetically conserved CHR promoter elements. We find 210 target genes that are expected to be regulated by the p53-p21-DREAM-CDE/CHR pathway. The target gene list was verified by detailed analysis of p53-dependent repression of the cell cycle genes B-MYB (MYBL2), BUB1, CCNA2, CCNB1, CHEK2, MELK, POLD1, RAD18 and RAD54L. Most of the 210 target genes are essential regulators of G₂ phase and mitosis. Thus, downregulation of these genes through the p53-p21-DREAM-CDE/CHR pathway appears to be a principal mechanism for G₂/M cell cycle arrest by p53.

The transcription factor p53: not a repressor, solely an activator

The predominant function of the tumor suppressor p53 is transcriptional regulation. It is generally accepted that p53-dependent transcriptional activation occurs by binding to a specific recognition site in promoters of target genes. Additionally, several models for p53-dependent transcriptional repression have been postulated. Here, we evaluate these models based on a computational meta-analysis of genome-wide data. Surprisingly, several major models of p53-dependent gene regulation are implausible. Meta-analysis of large-scale data is unable to confirm reports on directly repressed p53 target genes and falsifies models of direct repression. This notion is supported by experimental re-analysis of representative genes reported as directly repressed by p53. Therefore, p53 is not a direct repressor of transcription, but solely activates its target genes. Moreover, models based on interference of p53 with activating transcription factors as well as models based on the function of ncRNAs are also not supported by the meta-analysis. As an alternative to models of direct repression, the meta-analysis leads to the conclusion that p53 represses transcription indirectly by activation of the p53-p21-DREAM/RB pathway.

p53-Regulated Networks of Protein, mRNA, miRNA, and lncRNA Expression Revealed by Integrated Pulsed Stable Isotope Labeling With Amino Acids in Cell Culture (pSILAC) and Next Generation Sequencing (NGS) Analyses

We determined the effect of p53 activation on de novo protein synthesis using quantitative proteomics (pulsed stable isotope labeling with amino acids in cell culture/pSILAC) in the colorectal cancer cell line SW480. This was combined with mRNA and noncoding RNA expression analyses by next generation sequencing (RNA-, miR-Seq). Furthermore, genome-wide DNA binding of p53 was analyzed by chromatin-immunoprecipitation (ChIP-Seq). Thereby, we identified differentially regulated proteins (542 up, 569 down), mRNAs (1258 up, 415 down), miRNAs (111 up, 95 down) and lncRNAs (270 up, 123 down). Changes in protein and mRNA expression levels showed a positive correlation (r = 0.50, p < 0.0001). In total, we detected 133 direct p53 target genes that were differentially expressed and displayed p53 occupancy in the vicinity of their promoter. More transcriptionally induced genes displayed occupied p53 binding sites (4.3% mRNAs, 7.2% miRNAs, 6.3% lncRNAs, 5.9% proteins) than repressed genes (2.4% mRNAs, 3.2% miRNAs, 0.8% lncRNAs, 1.9% proteins), suggesting indirect mechanisms of repression. Around 50% of the down-regulated proteins displayed seed-matching sequences of p53-induced miRNAs in the corresponding 3'-UTRs. Moreover, proteins repressed by p53 significantly overlapped with those previously shown to be repressed by miR-34a. We confirmed up-regulation of the novel direct p53 target genes LINC01021, MDFI, ST14 and miR-486 and showed that ectopic LINC01021 expression inhibits proliferation in SW480 cells. Furthermore, KLF12, HMGB1 and CIT mRNAs were confirmed as direct targets of the p53-induced miR-34a, miR-205 and miR-486-5p, respectively. In line with the loss of p53 function during tumor progression, elevated expression of KLF12, HMGB1 and CIT was detected in advanced stages of cancer. In conclusion, the integration of multiple omics methods allowed the comprehensive identification of direct and indirect effectors of p53 that provide new insights and leads into the mechanisms of p53-mediated tumor suppression.

Quantitative Analysis of NF-κB Transactivation Specificity Using a Yeast-Based Functional Assay

The NF-κB transcription factor family plays a central role in innate immunity and inflammation processes and is frequently dysregulated in cancer. We developed an NF-κB functional assay in yeast to investigate the following issues: transactivation specificity of NF-κB proteins acting as homodimers or heterodimers; correlation between transactivation capacity and in vitro DNA binding measurements; impact of co-expressed interacting proteins or of small molecule inhibitors on NF-κB-dependent transactivation. Full-length p65 and p50 cDNAs were cloned into centromeric expression vectors under inducible GAL1 promoter in order to vary their expression levels. Since p50 lacks a transactivation domain (TAD), a chimeric construct containing the TAD derived from p65 was also generated (p50TAD) to address its binding and transactivation potential. The p50TAD and p65 had distinct transactivation specificities towards seventeen different κB response elements (κB-REs) where single nucleotide changes could greatly impact transactivation. For four κB-REs, results in yeast were predictive of transactivation potential measured in the human MCF7 cell lines treated with the NF-κB activator TNFα. Transactivation results in yeast correlated only partially with in vitro measured DNA binding affinities, suggesting that features other than strength of interaction with naked DNA affect transactivation, although factors such as chromatin context are kept constant in our isogenic yeast assay. The small molecules BAY11-7082 and ethyl-pyruvate as well as expressed IkBα protein acted as NF-κB inhibitors in yeast, more strongly towards p65. Thus, the yeast-based system can recapitulate NF-κB features found in human cells, thereby providing opportunities to address various NF-κB functions, interactions and chemical modulators.

Whole-genome cartography of p53 response elements ranked on transactivation potential

Background

Many recent studies using ChIP-seq approaches cross-referenced to trascriptome data and also to potentially unbiased in vitro DNA binding selection experiments are detailing with increasing precision the p53-directed gene regulatory network that, nevertheless, is still expanding. However, most experiments have been conducted in established cell lines subjected to specific p53-inducing stimuli, both factors potentially biasing the results.

Results

We developed p53retriever, a pattern search algorithm that maps p53 response elements (REs) and ranks them according to predicted transactivation potentials in five classes. Besides canonical, full site REs, we developed specific pattern searches for non-canonical half sites and 3/4 sites and show that they can mediate p53-dependent responsiveness of associated coding sequences. Using ENCODE data, we also mapped p53 REs in about 44,000 distant enhancers and identified a 16-fold enrichment for high activity REs within those sites in the comparison with genomic regions near transcriptional start sites (TSS). Predictions from our pattern search were cross-referenced to ChIP-seq, ChIP-exo, expression, and various literature data sources. Based on the mapping of predicted functional REs near TSS, we examined expression changes of thirteen genes as a function of different p53-inducing conditions, providing further evidence for PDE2A, GAS6, E2F7, APOBEC3H, KCTD1, TRIM32, DICER, HRAS, KITLG and TGFA p53-dependent regulation, while MAP2K3, DNAJA1 and potentially YAP1 were identified as new direct p53 target genes.

Conclusions

We provide a comprehensive annotation of canonical and non-canonical p53 REs in the human genome, ranked on predicted transactivation potential. We also establish or corroborate direct p53 transcriptional control of thirteen genes. The entire list of identified and functionally classified p53 REs near all UCSC-annotated genes and within ENCODE mapped enhancer elements is provided. Our approach is distinct from, and complementary to, existing methods designed to identify p53 response elements. p53retriever is available as an R package at: http://tomateba.github.io/p53retriever.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1643-9) contains supplementary material, which is available to authorized users.

Microarray and ChIP-seq data analysis revealed changes in p53-mediated transcriptional regulation in Nutlin-3-treated U2OS cells

Integrative analysis of chromatin immunoprecipitation-sequencing (ChIP-seq) data and microarray data was performed to illustrate the effect of Nutlin‑3 on promoter selectivity and transcriptional regulation by the tumor suppressor p53 in U2OS human osteosarcoma cells. Raw data (accession number, GSE46642) were downloaded from Gene Expression Omnibus. Differential analyses were performed using package limma of R software. Gene ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed for the differentially expressed genes (DEGs) using the Database for Annotation, Visualization and Integration Discovery. Integrative analysis of ChIP‑seq data and microarray data were confirmed with ChIP‑Array. A total of 565 DEGs were identified, including 373 upregulated genes and 192 downregulated genes. Genes involved in the p53 signaling pathway, cell cycle, DNA replication, cytokine‑cytokine receptor interaction and melanoma were markedly over‑represented in the DEGs. A total of 39 DEGs were directly regulated by p53 and two were the transcription factors (TFs), E2F2 and HOXA1. E2F2 regulated 25 DEGs, while HOXA1 regulated one DEG. The cell cycle, p53 signaling pathway, melanoma and pathways involved in cancer were enriched in the direct and indirect target genes. Changes in the p53‑binding pattern induced by Nutlin‑3 were described in the present study, which may advance the understanding of the regulatory network of p53 in osteosarcoma and aid in the development of novel therapies.

Integrative genomic analysis reveals widespread enhancer regulation by p53 in response to DNA damage

The tumor suppressor p53 has been studied extensively as a direct transcriptional activator of protein-coding genes. Recent studies, however, have shed light on novel regulatory functions of p53 within noncoding regions of the genome. Here, we use a systematic approach that integrates transcriptome-wide expression analysis, genome-wide p53 binding profiles and chromatin state maps to characterize the global regulatory roles of p53 in response to DNA damage. Notably, our approach identified conserved features of the p53 network in both human and mouse primary fibroblast models. In addition to known p53 targets, we identify many previously unappreciated mRNAs and long noncoding RNAs that are regulated by p53. Moreover, we find that p53 binding occurs predominantly within enhancers in both human and mouse model systems. The ability to modulate enhancer activity offers an additional layer of complexity to the p53 network and greatly expands the diversity of genomic elements directly regulated by p53.