Functional evolution of the p53 regulatory network through its target response elements

Human HDAC6 senses valine abundancy to regulate DNA damage

As an essential branched amino acid, valine is pivotal for protein synthesis, neurological behaviour, haematopoiesis and leukaemia progression1-3. However, the mechanism by which cellular valine abundancy is sensed for subsequent cellular functions remains undefined. Here we identify that human histone deacetylase 6 (HDAC6) serves as a valine sensor by directly binding valine through a primate-specific SE14 repeat domain. The nucleus and cytoplasm shuttling of human, but not mouse, HDAC6 is tightly controlled by the intracellular levels of valine. Valine deprivation leads to HDAC6 retention in the nucleus and induces DNA damage. Mechanistically, nuclear-localized HDAC6 binds and deacetylates ten-eleven translocation 2 (TET2) to initiate active DNA demethylation, which promotes DNA damage through thymine DNA glycosylase-driven excision. Dietary valine restriction inhibits tumour growth in xenograft and patient-derived xenograft models, and enhances the therapeutic efficacy of PARP inhibitors. Collectively, our study identifies human HDAC6 as a valine sensor that mediates active DNA demethylation and DNA damage in response to valine deprivation, and highlights the potential of dietary valine restriction for cancer treatment.

p53motifDB: integration of genomic information and tumor suppressor p53 binding motifs

The tumor suppressor gene TP53 encodes the DNA binding transcription factor p53 and is one of the most commonly mutated genes in human cancer. Tumor suppressor activity requires binding of p53 to its DNA response elements and subsequent transcriptional activation of a diverse set of target genes. Despite decades of close study, the logic underlying p53 interactions with its numerous potential genomic binding sites and target genes is not yet fully understood. Here, we present a database of DNA and chromatin-based information focused on putative p53 binding sites in the human genome to allow users to generate and test new hypotheses related to p53 activity in the genome. Users can query genomic locations based on experimentally observed p53 binding, regulatory element activity, genetic variation, evolutionary conservation, chromatin modification state, and chromatin structure. We present multiple use cases demonstrating the utility of this database for generating novel biological hypotheses, such as chromatin-based determinants of p53 binding and potential cell type-specific p53 activity. All database information is also available as a precompiled sqlite database for use in local analysis or as a Shiny web application.

TP53 germline pathogenic variants in modern humans were likely originated during recent human history

TP53 is crucial for maintaining genome stability and preventing oncogenesis. Germline pathogenic variation in TP53 damages its function, causing genome instability and increased cancer risk. Despite extensive study in TP53, the evolutionary origin of the human TP53 germline pathogenic variants remains largely unclear. In this study, we applied phylogenetic and archaeological approaches to identify the evolutionary origin of TP53 germline pathogenic variants in modern humans. In the phylogenic analysis, we searched 406 human TP53 germline pathogenic variants in 99 vertebrates distributed in eight clades of Primate, Euarchontoglires, Laurasiatheria, Afrotheria, Mammal, Aves, Sarcopterygii and Fish, but we observed no direct evidence for the cross-species conservation as the origin; in the archaeological analysis, we searched the variants in 5031 ancient human genomes dated between 45045 and 100 years before present, and identified 45 pathogenic variants in 62 ancient humans dated mostly within the last 8000 years; we also identified 6 pathogenic variants in 3 Neanderthals dated 44000 to 38515 years before present and 1 Denisovan dated 158 550 years before present. Our study reveals that TP53 germline pathogenic variants in modern humans were likely originated in recent human history and partially inherited from the extinct Neanderthals and Denisovans.

The complex architecture of p53 binding sites

Sequence-specific protein-DNA interactions are at the heart of the response of the tumor-suppressor p53 to numerous physiological and stress-related signals. Large variability has been previously reported in p53 binding to and transactivating from p53 response elements (REs) due, at least in part, to changes in direct (base) and indirect (shape) readouts of p53 REs. Here, we dissect p53 REs to decipher the mechanism by which p53 optimizes this highly regulated variable level of interaction with its DNA binding sites. We show that hemi-specific binding is more prevalent in p53 REs than previously envisioned. We reveal that sequences flanking the REs modulate p53 binding and activity and show that these effects extend to 4–5 bp from the REs. Moreover, we show here that the arrangement of p53 half-sites within its REs, relative to transcription direction, has been fine-tuned by selection pressure to optimize and regulate the response levels from p53 REs. This directionality in the REs arrangement is at least partly encoded in the structural properties of the REs. Furthermore, we show here that in the p21-5′ RE the orientation of the half-sites is such that the effect of the flanking sequences is minimized and we discuss its advantages.

Tumor suppressor p53: from engaging DNA to target gene regulation

The p53 transcription factor confers its potent tumor suppressor functions primarily through the regulation of a large network of target genes. The recent explosion of next generation sequencing protocols has enabled the study of the p53 gene regulatory network (GRN) and underlying mechanisms at an unprecedented depth and scale, helping us to understand precisely how p53 controls gene regulation. Here, we discuss our current understanding of where and how p53 binds to DNA and chromatin, its pioneer-like role, and how this affects gene regulation. We provide an overview of the p53 GRN and the direct and indirect mechanisms through which p53 affects gene regulation. In particular, we focus on delineating the ubiquitous and cell type-specific network of regulatory elements that p53 engages; reviewing our understanding of how, where, and when p53 binds to DNA and the mechanisms through which these events regulate transcription. Finally, we discuss the evolution of the p53 GRN and how recent work has revealed remarkable differences between vertebrates, which are of particular importance to cancer researchers using mouse models.

Regulating tumor suppressor genes: post-translational modifications

Tumor suppressor genes cooperate with each other in tumors. Three important tumor suppressor proteins, retinoblastoma (Rb), p53, phosphatase, and tensin homolog deleted on chromosome ten (PTEN) are functionally associated and they regulated by post-translational modification (PTMs) as well. PTMs include phosphorylation, SUMOylation, acetylation, and other novel modifications becoming growing appreciated. Because most of PTMs are reversible, normal cells use them as a switch to control the state of cells being the resting or proliferating, and PTMs also involve in cell survival and cell cycle, which may lead to abnormal proliferation and tumorigenesis. Although a lot of studies focus on the importance of each kind of PTM, further discoveries shows that tumor suppressor genes (TSGs) form a complex "network" by the interaction of modification. Recently, there are several promising strategies for TSGs for they change more frequently than carcinogenic genes in cancers. We here review the necessity, characteristics, and mechanisms of each kind of post-translational modification on Rb, p53, PTEN, and its influence on the precise and selective function. We also discuss the current antitumoral therapies of Rb, p53 and PTEN as predictive, prognostic, and therapeutic target in cancer.

A Protein in the Yeast Saccharomyces cerevisiae Presents DNA Binding Homology to the p53 Checkpoint Protein and Tumor Suppressor

Saccharomyces cerevisiae does not contain a p53 homolog. Utilizing this yeast as an in vivo test tube model, our aim was to investigate if a yeast protein would show p53 DNA binding homology. Electrophoretic mobility shift analyses revealed the formation of specific DNA-protein complexes consisting of S. cerevisiae nuclear protein(s) and oligonucleotides containing p53 DNA binding sites. A S. cerevisiae p53 binding site factor (Scp53BSF) bound to a p53 synthetic DNA-consensus sequence (SCS) and a p53 binding-site sequence from the MDM2 oncogene. The complexes were of comparable size. Like mammalian p53, the affinity of Scp53BSF for the SCS oligonucleotide was higher than for the MDM2 oligonucleotide. Binding of Scp53BSF to the SCS and MDM2 oligonucleotides was strongly competed by unlabeled oligonucleotides containing mammalian p53 sites, but very little by a mutated site oligonucleotide. Importantly, Scp53BSF-DNA binding activity was significantly induced in extracts from cells with DNA damage. This resulted in dose-dependent coordinated activation of transcription when using p53-binding site reporter constructs. An ancient p53-like DNA binding protein may have been found, and activation of DNA-associated factors to p53 response elements may have functions not yet determined.

Discordant gene responses to radiation in humans and mice and the role of hematopoietically humanized mice in the search for radiation biomarkers

The mouse (Mus musculus) is an extensively used model of human disease and responses to stresses such as ionizing radiation. As part of our work developing gene expression biomarkers of radiation exposure, dose, and injury, we have found many genes are either up-regulated (e.g. CDKN1A, MDM2, BBC3, and CCNG1) or down-regulated (e.g. TCF4 and MYC) in both species after irradiation at ~4 and 8 Gy. However, we have also found genes that are consistently up-regulated in humans and down-regulated in mice (e.g. DDB2, PCNA, GADD45A, SESN1, RRM2B, KCNN4, IFI30, and PTPRO). Here we test a hematopoietically humanized mouse as a potential in vivo model for biodosimetry studies, measuring the response of these 14 genes one day after irradiation at 2 and 4 Gy, and comparing it with that of human blood irradiated ex vivo, and blood from whole body irradiated mice. We found that human blood cells in the hematopoietically humanized mouse in vivo environment recapitulated the gene expression pattern expected from human cells, not the pattern seen from in vivo irradiated normal mice. The results of this study support the use of hematopoietically humanized mice as an in vivo model for screening of radiation response genes relevant to humans.

The p53 mRNA: an integral part of the cellular stress response

A large number of signalling pathways converge on p53 to induce different cellular stress responses that aim to promote cell cycle arrest and repair or, if the damage is too severe, to induce irreversible senescence or apoptosis. The differentiation of p53 activity towards specific cellular outcomes is tightly regulated via a hierarchical order of post-translational modifications and regulated protein-protein interactions. The mechanisms governing these processes provide a model for how cells optimize the genetic information for maximal diversity. The p53 mRNA also plays a role in this process and this review aims to illustrate how protein and RNA interactions throughout the p53 mRNA in response to different signalling pathways control RNA stability, translation efficiency or alternative initiation of translation. We also describe how a p53 mRNA platform shows riboswitch-like features and controls the rate of p53 synthesis, protein stability and modifications of the nascent p53 protein. A single cancer-derived synonymous mutation disrupts the folding of this platform and prevents p53 activation following DNA damage. The role of the p53 mRNA as a target for signalling pathways illustrates how mRNA sequences have co-evolved with the function of the encoded protein and sheds new light on the information hidden within mRNAs.

Shaping the regulation of the p53 mRNA tumour suppressor: the co-evolution of genetic signatures

Structured RNA regulatory motifs exist from the prebiotic stages of the RNA world to the more complex eukaryotic systems. In cases where a functional RNA structure is within the coding sequence a selective pressure drives a parallel co-evolution of the RNA structure and the encoded peptide domain. The p53-MDM2 axis, describing the interactions between the p53 tumor suppressor and the MDM2 E3 ubiquitin ligase, serves as particularly useful model revealing how secondary RNA structures have co-evolved along with corresponding interacting protein motifs, thus having an impact on protein - RNA and protein - protein interactions; and how such structures developed signal-dependent regulation in mammalian systems. The p53(BOX-I) RNA sequence binds the C-terminus of MDM2 and controls p53 synthesis while the encoded peptide domain binds MDM2 and controls p53 degradation. The BOX-I peptide domain is also located within p53 transcription activation domain. The folding of the p53 mRNA structure has evolved from temperature-regulated in pre-vertebrates to an ATM kinase signal-dependent pathway in mammalian cells. The protein - protein interaction evolved in vertebrates and became regulated by the same signaling pathway. At the same time the protein - RNA and protein - protein interactions evolved, the p53 trans-activation domain progressed to become integrated into a range of cellular pathways. We discuss how a single synonymous mutation in the BOX-1, the p53(L22 L), observed in a chronic lymphocyte leukaemia patient, prevents the activation of p53 following DNA damage. The concepts analysed and discussed in this review may serve as a conceptual mechanistic paradigm of the co-evolution and function of molecules having roles in cellular regulation, or the aetiology of genetic diseases and how synonymous mutations can affect the encoded protein.

Conservation and divergence of the p53 gene regulatory network between mice and humans

Understanding the p53 tumor suppressor pathway remains crucial for the design of anticancer strategies. Studies in human tumors and mouse models help to unravel the molecular mechanisms that underlie the p53 signaling pathway. Yet, the p53 gene regulatory network (GRN) is not the same in mice and humans. The comparison of the regulatory networks of p53 in mice and humans reveals that gene up- and down-regulation by p53 are distinctly affected during evolution. Importantly, gene up-regulation by p53 underwent more rapid evolution and gene down-regulation has been evolutionarily constrained. This difference stems from the two major mechanisms employed by p53 to regulate gene expression: up-regulation through direct p53 target gene binding and indirect down-regulation through the p53-p21-DREAM pathway. More than 1000 genes have been identified to differ in their p53-dependent expression between mice and humans. Analysis of p53 gene expression profiles and p53 binding data reveal that turnover of p53 binding sites is the major mechanism underlying extensive variation in p53-dependent gene up-regulation. Only a core set of high-confidence genes appears to be directly regulated by p53 in both species. In contrast to up-regulation, p53-induced down-regulation is well conserved between mice and humans and controls cell cycle genes. Here a curated data set is provided that extends the previously established web-atlas at www.targetgenereg.org to assess the p53 response of any human gene of interest and its mouse ortholog. Taken together, the analysis reveals a limited translation potential from mouse models to humans for the p53 GRN.

DNA Damage Response and Repair, DNA Methylation, and Cell Death in Human Neurons and Experimental Animal Neurons Are Different

Neurological disorders affecting individuals in infancy to old age elude interventions for meaningful protection against neurodegeneration, and preclinical work has not translated to humans. We studied human neuron responses to injury and death stimuli compared to those of animal neurons in culture under similar settings of insult (excitotoxicity, oxidative stress, and DNA damage). Human neurons were differentiated from a cortical neuron cell line and the embryonic stem cell-derived H9 line. Mouse neurons were differentiated from forebrain neural stem cells and embryonic cerebral cortex; pig neurons were derived from forebrain neural stem cells. Mitochondrial morphology was different in human and mouse neurons. Human and mouse neurons challenged with DNA-damaging agent camptothecin showed different chromatin condensation, cell death, and DNA damage sensor activation. DNA damage accumulation and repair kinetics differed among human, mouse, and pig neurons. Promoter CpG island methylation microarrays showed significant differential DNA methylation in human and mouse neurons after injury. Therefore, DNA damage response, DNA repair, DNA methylation, and autonomous cell death mechanisms in human neurons and experimental animal neurons are different.

A highly specific SpCas9 variant is identified by in vivo screening in yeast

Despite the utility of CRISPR-Cas9 nucleases for genome editing, the potential for off-target activity limits their application, especially for therapeutic purposes. We developed a yeast-based assay to identify optimized Streptococcus pyogenes Cas9 (SpCas9) variants that enables simultaneous evaluation of on- and off-target activity. We screened a library of SpCas9 variants carrying random mutations in the REC3 domain and identified mutations that increased editing accuracy while maintaining editing efficiency. We combined four beneficial mutations to generate evoCas9, a variant that has fidelity exceeding both wild-type (79-fold improvement) and rationally designed Cas9 variants (fourfold average improvement), while maintaining near wild-type on-target editing efficiency (90% median residual activity). Evaluating evoCas9 on endogenous genomic loci, we demonstrated a substantially improved specificity and observed no off-target sites for four of the eight single guide RNAs (sgRNAs) tested. Finally, we showed that following long-term expression (40 d), evoCas9 strongly limited the nonspecific cleavage of a difficult-to-discriminate off-target site and fully abrogated the cleavage of two additional off-target sites.

Conservation and Divergence of p53 Oscillation Dynamics across Species

The tumor-suppressing transcription factor p53 is highly conserved at the protein level and plays a key role in the DNA damage response. One important aspect of p53 regulation is its dynamics in response to DNA damage, which include oscillations. Here, we observe that, while the qualitative oscillatory nature of p53 dynamics is conserved across cell lines derived from human, monkey, dog, mouse, and rat, the oscillation period is variable. Specifically, rodent cells exhibit rapid p53 oscillations, whereas dog, monkey, and human cells show slower oscillations. Computational modeling and experiments identify stronger negative feedback between p53 and MDM2 as the driver of faster oscillations in rodents, suggesting that the period of oscillation is a network-level property. In total, our study shows that despite highly conserved signaling, the quantitative features of p53 oscillations can diverge across evolution. We caution that strong amino acid conservation of proteins and transcriptional network similarity do not necessarily imply conservation of time dynamics.

Human transcription factors in yeast: the fruitful examples of P53 and NF-кB

The observation that human transcription factors (TFs) can function when expressed in yeast cells has stimulated the development of various functional assays to investigate (i) the role of binding site sequences (herein referred to as response elements, REs) in transactivation specificity, (ii) the impact of polymorphic nucleotide variants on transactivation potential, (iii) the functional consequences of mutations in TFs and (iv) the impact of cofactors or small molecules. These approaches have found applications in basic as well as applied research, including the identification and the characterisation of mutant TF alleles from clinical samples. The ease of genome editing of yeast cells and the availability of regulated systems for ectopic protein expression enabled the development of quantitative reporter systems, integrated at a chosen chromosomal locus in isogenic yeast strains that differ only at the level of a specific RE targeted by a TF or for the expression of distinct TF alleles. In many cases, these assays were proven predictive of results in higher eukaryotes. The potential to work in small volume formats and the availability of yeast strains with modified chemical uptake have enhanced the scalability of these approaches. Next to well-established one-, two-, three-hybrid assays, the functional assays with non-chimeric human TFs enrich the palette of opportunities for functional characterisation. We review ∼25 years of research on human sequence-specific TFs expressed in yeast, with an emphasis on the P53 and NF-кB family of proteins, highlighting outcomes, advantages, challenges and limitations of these heterologous assays.

Identification of p53-target genes in Danio rerio

To orchestrate the genomic response to cellular stress signals, p53 recognizes and binds to DNA containing specific and well-characterized p53-responsive elements (REs). Differences in RE sequences can strongly affect the p53 transactivation capacity and occur even between closely related species. Therefore, the identification and characterization of a species-specific p53 Binding sistes (BS) consensus sequence and of the associated target genes may help to provide new insights into the evolution of the p53 regulatory networks across different species. Although p53 functions were studied in a wide range of species, little is known about the p53-mediated transcriptional signature in Danio rerio. Here, we designed and biochemically validated a computational approach to identify novel p53 target genes in Danio rerio genome. Screening all the Danio rerio genome by pattern-matching-based analysis, we found p53 RE-like patterns proximal to 979 annotated Danio rerio genes. Prioritization analysis identified a subset of 134 candidate pattern-related genes, 31 of which have been investigated in further biochemical assays. Our study identified runx1, axin1, traf4a, hspa8, col4a5, necab2, and dnajc9 genes as novel direct p53 targets and 12 additional p53-controlled genes in Danio rerio genome. The proposed combinatorial approach resulted to be highly sensitive and robust for identifying new p53 target genes also in additional animal species.

The DNA damage-induced cell death response: a roadmap to kill cancer cells

Upon massive DNA damage cells fail to undergo productive DNA repair and trigger the cell death response. Resistance to cell death is linked to cellular transformation and carcinogenesis as well as radio- and chemoresistance, making the underlying signaling pathways a promising target for therapeutic intervention. Diverse DNA damage-induced cell death pathways are operative in mammalian cells and finally culminate in the induction of programmed cell death via activation of apoptosis or necroptosis. These signaling routes affect nuclear, mitochondria- and plasma membrane-associated key molecules to activate the apoptotic or necroptotic response. In this review, we highlight the main signaling pathways, molecular players and mechanisms guiding the DNA damage-induced cell death response.

Two p53 tetramers bind one consensus DNA response element

p53 tumor suppressor is a transcription factor that controls cell cycle and genetic integrity. In response to genotoxic stress p53 activates DNA repair, cell cycle arrest, apoptosis or senescence, which are initiated via p53 binding to its specific DNA response elements (RE). The consensus p53 DNA RE consists of two decameric palindromic half-site sequences. Crystallographic studies have demonstrated that two isolated p53 DNA-binding core domains interact with one half-site of the p53 DNA REs suggesting that one p53 tetramer is bound to one RE. However, our recent 3D cryo-EM studies showed that the full-length p53 tetramer is bound to only one half-site of RE.

Here, we have used biochemical and electron microscopy (EM) methods to analyze DNA-binding of human and murine p53 tetramers to various p53 DNA REs. Our new results demonstrate that two p53 tetramers can interact sequence-specifically with one DNA RE at the same time. In particular, the EM structural analysis revealed that two p53 tetramers bind one DNA RE simultaneously with DNA positioned between them. These results demonstrate a mode different from that assumed previously for the p53-DNA interaction and suggest important biological implications on p53 activity as a transcriptional regulator of cellular response to stress.

TLR2∆22 (-196-174) significantly increases the risk of breast cancer in females carrying proline allele at codon 72 of TP53 gene: a case-control study from four ethnic groups of North Eastern region of India

Breast cancer (BC) is the second most common cancer in women. In the North Eastern Region (NER) of India, BC is emerging as an important concern as evidenced by the data available from population and hospital-based cancer registries. Studies on genetic susceptibility to BC are important to understand the increase in the incidence of BC in NER. The present case control study was conducted to investigate the association between tumour suppressor gene TP53 codon 72 polymorphism and innate immune pathway gene TLR2∆22 (-196-174) polymorphism with BC in females of NER of India for the identification of novel biomarker of BC. Four hundred sixty-two histopathologically confirmed BC cases from four states of NER of India, and 770 healthy controls were included by organizing community surveys from the neighbourhood of cases. In our study, no significant association between TP53 codon 72 polymorphisms and the risk of BC was found. However, our study has shown that TP53 codon 72 polymorphism is an important effect modifier. In the present study it was found that females carrying 22 base-pair deletion in the promoter region of their TLR2 gene had two times (AOR= 2.18, 95 % CI 1.13-4.21, p=0.019 in dominant model; AOR= 2.17, 95 % CI 1.09-4.34, p=0.027 in co-dominant model) increased risk of BC whwn they also carry proline allele at codon 72 of their TP53 gene.

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.

Evolution of p53 transactivation specificity through the lens of a yeast-based functional assay

Co-evolution of transcription factors (TFs) with their respective cis-regulatory network enhances functional diversity in the course of evolution. We present a new approach to investigate transactivation capacity of sequence-specific TFs in evolutionary studies. Saccharomyces cerevisiae was used as an in vivo test tube and p53 proteins derived from human and five commonly used animal models were chosen as proof of concept. p53 is a highly conserved master regulator of environmental stress responses. Previous reports indicated conserved p53 DNA binding specificity in vitro, even for evolutionary distant species. We used isogenic yeast strains where p53-dependent transactivation was measured towards chromosomally integrated p53 response elements (REs). Ten REs were chosen to sample a wide range of DNA binding affinity and transactivation capacity for human p53 and proteins were expressed at two levels using an inducible expression system. We showed that the assay is amenable to study thermo-sensitivity of frog p53, and that chimeric constructs containing an ectopic transactivation domain could be rapidly developed to enhance the activity of proteins, such as fruit fly p53, that are poorly effective in engaging the yeast transcriptional machinery. Changes in the profile of relative transactivation towards the ten REs were measured for each p53 protein and compared to the profile obtained with human p53. These results, which are largely independent from relative p53 protein levels, revealed widespread evolutionary divergence of p53 transactivation specificity, even between human and mouse p53. Fruit fly and human p53 exhibited the largest discrimination among REs while zebrafish p53 was the least selective.

∆N-P63α and TA-P63α exhibit intrinsic differences in transactivation specificities that depend on distinct features of DNA target sites

TP63 is a member of the TP53 gene family that encodes for up to ten different TA and ΔN isoforms through alternative promoter usage and alternative splicing. Besides being a master regulator of gene expression for squamous epithelial proliferation, differentiation and maintenance, P63, through differential expression of its isoforms, plays important roles in tumorigenesis. All P63 isoforms share an immunoglobulin-like folded DNA binding domain responsible for binding to sequence-specific response elements (REs), whose overall consensus sequence is similar to that of the canonical p53 RE. Using a defined assay in yeast, where P63 isoforms and RE sequences are the only variables, and gene expression assays in human cell lines, we demonstrated that human TA- and ΔN-P63α proteins exhibited differences in transactivation specificity not observed with the corresponding P73 or P53 protein isoforms. These differences 1) were dependent on specific features of the RE sequence, 2) could be related to intrinsic differences in their oligomeric state and cooperative DNA binding, and 3) appeared to be conserved in evolution. Since genotoxic stress can change relative ratio of TA- and ΔN-P63α protein levels, the different transactivation specificity of each P63 isoform could potentially influence cellular responses to specific stresses.

Codon 104 variation of p53 gene provides adaptive apoptotic responses to extreme environments in mammals of the Tibet plateau

Mutational changes in p53 correlate well with tumorigenesis. Remarkably, however, relatively little is known about the role that p53 variations may play in environmental adaptation. Here we report that codon asparagine-104 (104N) and glutamic acid-104 (104E), respectively, of the p53 gene in the wild zokor (Myospalax baileyi) and root vole (Microtus oeconomus) are adaptively variable, meeting the environmental stresses of the Tibetan plateau. They differ from serine-104 (104S) seen in other rodents, including the lowland subterranean zokor Myospalax cansus, and from serine 106 (106S) in humans. Based on site-directed mutational analysis in human cell lines, the codon 104N variation in M. baileyi is responsible for the adaptive balance of the transactivation of apoptotic genes under hypoxia, cold, and acidic stresses. The 104E p53 variant in Microtus oeconomus suppresses apoptotic gene transactivation and cell apoptosis. Neither 104N nor 104E affects the cell-cycle genes. We propose that these variations in p53 codon 104 are an outcome of environmental adaptation and evolutionary selection that enhance cellular strategies for surviving the environmental stresses of hypoxia and cold (in M. baileyi and M. oeconomus) and hypercapnia (in M. baileyi) in the stressful environments of the Qinghai-Tibet plateau.

Identification of new p53 target microRNAs by bioinformatics and functional analysis

BACKGROUND

The tumor suppressor p53 is a sequence-specific transcription factor that regulates an extensive network of coding genes, long non-coding RNAs and microRNAs, that establish intricate gene regulatory circuits influencing many cellular responses beyond the prototypical control of cell cycle, apoptosis and DNA repair.

METHODS

Using bioinformatic approaches, we identified an additional group of candidate microRNAs (miRs) under direct p53 transcriptional control. To validate p53 family-mediated responsiveness of the newly predicted target miRs we first evaluated the potential for wild type p53, p63β and p73β to transactivate from p53 response elements (REs) mapped in the miR promoters, using an established yeast-based assay.

RESULTS

The REs found in miR-10b, -23b, -106a, -151a, -191, -198, -202, -221, -320, -1204, -1206 promoters were responsive to p53 and 8 of them were also responsive to p63β or p73β. The potential for germline p53 mutations to drive transactivation at selected miR-associated REs was also examined. Chromatin Immuno-Precipitation (ChIP) assays conducted in doxorubicin-treated MCF7 cells and HCT116 p53+/+ revealed moderate induction of p53 occupancy at the miR-202, -1204, -1206, -10b RE-containing sites, while weak occupancy was observed for the miR-23b-associated RE only in MCF7 cells. RT-qPCR analyses cells showed modest doxorubicin- and/or Nutlin-dependent induction of the levels of mature miR-10b, -23b, -151a in HCT116 p53+/+ and MCF7 cells. The long noncoding RNA PVT1 comprising miR-1204 and -1206 was weakly induced only in HCT116 p53+/+ cells, but the mature miRs were not detected. miR-202 expression was not influenced by p53-activating stimuli in our cell systems.

CONCLUSIONS

Our study reveals additional miRs, particularly miR-10b and miR-151a, that could be directly regulated by the p53-family of transcription factors and contribute to the tuning of p53-induced responses.

Rotational positioning of nucleosomes facilitates selective binding of p53 to response elements associated with cell cycle arrest

The tumor suppressor protein p53 exhibits high affinity to the response elements regulating cell cycle arrest genes (CCA-sites), but relatively low affinity to the sites associated with apoptosis (Apo-sites). This in vivo tendency cannot be explained solely by the p53-DNA binding constants measured in vitro. Since p53 can bind nucleosomal DNA, we sought to understand if the two groups of p53 sites differ in their accessibility when embedded in nucleosomes. To this aim, we analyzed the sequence-dependent bending anisotropy of human genomic DNA containing p53 sites. For the 20 CCA-sites, we calculated rotational positioning patterns predicting that most of the sites are exposed on the nucleosomal surface. This is consistent with experimentally observed positioning of human nucleosomes. Remarkably, the sequence-dependent DNA anisotropy of both the p53 sites and flanking DNA work in concert producing strong positioning signals. By contrast, both the predicted and observed rotational settings of the 38 Apo-sites in nucleosomes suggest that many of these sites are buried inside, thus preventing immediate p53 recognition and delaying gene induction. The distinct chromatin organization of the CCA response elements appears to be one of the key factors facilitating p53-DNA binding and subsequent activation of genes associated with cell cycle arrest.

Characterization of the p53 cistrome--DNA binding cooperativity dissects p53's tumor suppressor functions

p53 protects us from cancer by transcriptionally regulating tumor suppressive programs designed to either prevent the development or clonal expansion of malignant cells. How p53 selects target genes in the genome in a context- and tissue-specific manner remains largely obscure. There is growing evidence that the ability of p53 to bind DNA in a cooperative manner prominently influences target gene selection with activation of the apoptosis program being completely dependent on DNA binding cooperativity. Here, we used ChIP-seq to comprehensively profile the cistrome of p53 mutants with reduced or increased cooperativity. The analysis highlighted a particular relevance of cooperativity for extending the p53 cistrome to non-canonical binding sequences characterized by deletions, spacer insertions and base mismatches. Furthermore, it revealed a striking functional separation of the cistrome on the basis of cooperativity; with low cooperativity genes being significantly enriched for cell cycle and high cooperativity genes for apoptotic functions. Importantly, expression of high but not low cooperativity genes was correlated with superior survival in breast cancer patients. Interestingly, in contrast to most p53-activated genes, p53-repressed genes did not commonly contain p53 binding elements. Nevertheless, both the degree of gene activation and repression were cooperativity-dependent, suggesting that p53-mediated gene repression is largely indirect and mediated by cooperativity-dependently transactivated gene products such as CDKN1A, E2F7 and non-coding RNAs. Since both activation of apoptosis genes with non-canonical response elements and repression of pro-survival genes are crucial for p53's apoptotic activity, the cistrome analysis comprehensively explains why p53-induced apoptosis, but not cell cycle arrest, strongly depends on the intermolecular cooperation of p53 molecules as a possible safeguard mechanism protecting from accidental cell killing.

The tumor suppressor gene p53 counteracts tumor growth by activating genes that prevent cell proliferation or induce cell death. How p53 selects genes in the genome to direct cell fate specifically into one or the other direction remains unclear. We show that the ability of p53 molecules to interact and thereby cooperate, influences which genes in the genome p53 is regulating. In the absence of cooperation, p53 only binds and regulates a limited ‘default’ set of genes that is proficient to stop cell proliferation but insufficient to induce cell death. Cooperation increases p53's DNA binding and enables context-dependent activation of apoptosis genes and repression of pro-survival genes which together triggers cell death. As the concerted effort of p53 molecules is needed, the threshold for cell killing is raised possibly to protect us from accidental cell loss. Thus, by shaping the genomic binding pattern, p53 cooperation fine-tunes the gene activity pattern to steer cell fate into the most appropriate, context-dependent direction. The genome-wide binding patterns of cooperating and non-cooperating p53 proteins generated in this study provide a comprehensive list of p53 binding sites as a resource for the scientific community to further explore mechanisms of tumor suppression by p53.

Transactivation specificity is conserved among p53 family proteins and depends on a response element sequence code

Structural and biochemical studies have demonstrated that p73, p63 and p53 recognize DNA with identical amino acids and similar binding affinity. Here, measuring transactivation activity for a large number of response elements (REs) in yeast and human cell lines, we show that p53 family proteins also have overlapping transactivation profiles. We identified mutations at conserved amino acids of loops L1 and L3 in the DNA-binding domain that tune the transactivation potential nearly equally in p73, p63 and p53. For example, the mutant S139F in p73 has higher transactivation potential towards selected REs, enhanced DNA-binding cooperativity in vitro and a flexible loop L1 as seen in the crystal structure of the protein–DNA complex. By studying, how variations in the RE sequence affect transactivation specificity, we discovered a RE-transactivation code that predicts enhanced transactivation; this correlation is stronger for promoters of genes associated with apoptosis.

P53 family members modulate the expression of PRODH, but not PRODH2, via intronic p53 response elements

The tumor suppressor p53 was previously shown to markedly up-regulate the expression of the PRODH gene, encoding the proline dehydrogenase (PRODH) enzyme, which catalyzes the first step in proline degradation. Also PRODH2, which degrades 4-hydroxy-L-proline, a product of protein (e.g. collagen) catabolism, was recently described as a p53 target. Here, we confirmed p53-dependent induction of endogenous PRODH in response to genotoxic damage in cell lines of different histological origin. We established that over-expression of TAp73β or TAp63β is sufficient to induce PRODH expression in p53-null cells and that PRODH expression parallels the modulation of endogenous p73 by genotoxic drugs in several cell lines. The p53, p63, and p73-dependent transcriptional activation was linked to specific intronic response elements (REs), among those predicted by bioinformatics tools and experimentally validated by a yeast-based transactivation assay. p53 occupancy measurements were validated in HCT116 and MCF7 human cell lines. Conversely, PRODH2 was not responsive to p63 nor p73 and, at best, could be considered a weak p53 target. In fact, minimal levels of PRODH2 transcript induction by genotoxic stress was observed exclusively in one of four p53 wild-type cell lines tested. Consistently, all predicted p53 REs in PRODH2 were poor matches to the p53 RE consensus and showed very weak responsiveness, only to p53, in the functional assay. Taken together, our results highlight that PRODH, but not PRODH2, expression is under the control of p53 family members, specifically p53 and p73. This supports a deeper link between proteins of the p53-family and metabolic pathways, as PRODH modulates the balance of proline and glutamate levels and those of their derivative alpha-keto-glutarate (α-KG) under normal and pathological (tumor) conditions.

Gene entropy-fractal dimension informatics with application to mouse-human translational medicine

DNA informatics represented by Shannon entropy and fractal dimension have been used to form 2D maps of related genes in various mammals. The distance between points on these maps for corresponding mRNA sequences in different species is used to study evolution. By quantifying the similarity of genes between species, this distance might be indicated when studies on one species (mouse) would tend to be valid in the other (human). The hypothesis that a small distance from mouse to human could facilitate mouse to human translational medicine success is supported by the studied ESR-1, LMNA, Myc, and RNF4 sequences. ID1 and PLCZ1 have larger separation. The collinearity of displacement vectors is further analyzed with a regression model, and the ID1 result suggests a mouse-chimp-human translational medicine approach. Further inference was found in the tumor suppression gene, p53, with a new hypothesis of including the bovine PKM2 pathways for targeting the glycolysis preference in many types of cancerous cells, consistent with quantum metabolism models. The distance between mRNA and protein coding CDS is proposed as a measure of the pressure associated with noncoding processes. The Y-chromosome DYS14 in fetal micro chimerism that could offer protection from Alzheimer's disease is given as an example.

It Takes 15 to Tango: Making Sense of the Many Ubiquitin Ligases of p53

The transcription factor p53 regulates numerous cellular processes to guard against tumorigenesis. Cell-cycle inhibition, apoptosis, and autophagy are all regulated by p53 in a cell- and context-specific manner, underscoring the need for p53 activity to be kept low in most circumstances. p53 is kept in check primarily through its regulated ubiquitination and degradation by a number of different factors, whose contributions may reflect complex context-specific needs to restrain p53 activity. Chief among these E3 ubiquitin ligases in p53 homeostasis is the ubiquitously expressed proto-oncogene MDM2, whose loss renders vertebrates unable to limit p53 activity, resulting in early embryonic lethality. MDM2 has been validated as a critical, universal E3 ubiquitin ligase for p53 in numerous tissues and organisms to date, but additional E3 ligases have also been identified for p53 whose contribution to p53 activity is unclear. In this review, we summarize the recent advances in our knowledge regarding how p53 activity is apparently controlled by a multitude of ubiquitin ligases beyond MDM2.

Comparative epigenomic annotation of regulatory DNA

Despite the explosive growth of genomic data, functional annotation of regulatory sequences remains difficult. Here, we introduce "comparative epigenomics"-interspecies comparison of DNA and histone modifications-as an approach for annotation of the regulatory genome. We measured in human, mouse, and pig pluripotent stem cells the genomic distributions of cytosine methylation, H2A.Z, H3K4me1/2/3, H3K9me3, H3K27me3, H3K27ac, H3K36me3, transcribed RNAs, and P300, TAF1, OCT4, and NANOG binding. We observed that epigenomic conservation was strong in both rapidly evolving and slowly evolving DNA sequences, but not in neutrally evolving sequences. In contrast, evolutionary changes of the epigenome and the transcriptome exhibited a linear correlation. We suggest that the conserved colocalization of different epigenomic marks can be used to discover regulatory sequences. Indeed, seven pairs of epigenomic marks identified exhibited regulatory functions during differentiation of embryonic stem cells into mesendoderm cells. Thus, comparative epigenomics reveals regulatory features of the genome that cannot be discerned from sequence comparisons alone.

The human TLR innate immune gene family is differentially influenced by DNA stress and p53 status in cancer cells

The transcription factor p53 regulates genes associated with a wide range of functions, including the Toll-like receptor (TLR) set of innate immunity genes, suggesting that p53 also modulates the human immune response. The TLR family comprises membrane glycoproteins that recognize pathogen-associated molecular patterns (PAMP) and mediate innate immune responses, and TLR agonists are being used as adjuvants in cancer treatments. Here, we show that doxorubicin, 5-fluorouracil, and UV and ionizing radiation elicit changes in TLR expression that are cell line- and damage-specific. Specifically, treatment-induced expression changes led to increased downstream cytokine expression in response to ligand stimulation. The effect of DNA stressors on TLR expression was mainly mediated by p53, and several p53 cancer-associated mutants dramatically altered the pattern of TLR gene expression. In all cell lines tested, TLR3 induction was p53-dependent, whereas induction of TLR9, the most stress-responsive family member, was less dependent on status of p53. In addition, each of the 10 members of the innate immune TLR gene family tested was differentially inducible. Our findings therefore show that the matrix of p53 status, chromosome stress, and responsiveness of individual TLRs should be considered in TLR-based cancer therapies.

The transcriptional response of mammalian cancer cells to irradiation is dominated by a cell cycle signature which is strongly attenuated in non-cancer cells and tissues

PURPOSE

Our goal was to identify genes showing a general transcriptional response to irradiation in mammalian cells and to analyze their response in function of dose, time and quality of irradiation and of cell type.

MATERIALS AND METHODS

We used a modified MIAME (Minimal Information About Microarray Experiments) protocol to import microarray data from 177 different irradiation conditions in the Radiation Genes database and performed cut-off-based selections and hierarchical gene clustering.

RESULTS

We identified a set of 29 genes which respond to a wide range of irradiation conditions in different cell types and tissues. Functional analysis of the negatively modulated genes revealed a dominant signature of mitotic cell cycle regulation which appears both dose and time-dependent. This signature is prominent in cancer cells and highly proliferating tissues but it is strongly attenuated in non cancer cells.

CONCLUSIONS

The transcriptional response of mammalian cancer cells to irradiation is dominated by a mitotic cell cycle signature both dose and time-dependent. This core response, which is present in cancer cells and highly proliferating tissues such as skin, blood and lymph node, is weaker or absent in non-cancer cells and in liver and spleen. CDKN1A (cyclin-dependent kinase inhibitor 1A) appears as the most generally induced mammalian gene and its response (mostly dose- and time-independent) seems to go beyond the typical DNA damage response.

Structure of p73 DNA-binding domain tetramer modulates p73 transactivation

The transcription factor p73 triggers developmental pathways and overlaps stress-induced p53 transcriptional pathways. How p53-family response elements determine and regulate transcriptional specificity remains an unsolved problem. In this work, we have determined the first crystal structures of p73 DNA-binding domain tetramer bound to response elements with spacers of different length. The structure and function of the adaptable tetramer are determined by the distance between two half-sites. The structures with zero and one base-pair spacers show compact p73 DNA-binding domain tetramers with large tetramerization interfaces; a two base-pair spacer results in DNA unwinding and a smaller tetramerization interface, whereas a four base-pair spacer hinders tetramerization. Functionally, p73 is more sensitive to spacer length than p53, with one base-pair spacer reducing 90% of transactivation activity and longer spacers reducing transactivation to basal levels. Our results establish the quaternary structure of the p73 DNA-binding domain required as a scaffold to promote transactivation.

Functional analyses of human DNA repair proteins important for aging and genomic stability using yeast genetics

Model systems have been extremely useful for studying various theories of aging. Studies of yeast have been particularly helpful to explore the molecular mechanisms and pathways that affect aging at the cellular level in the simple eukaryote. Although genetic analysis has been useful to interrogate the aging process, there has been both interest and debate over how functionally conserved the mechanisms of aging are between yeast and higher eukaryotes, especially mammalian cells. One area of interest has been the importance of genomic stability for age-related processes, and the potential conservation of proteins and pathways between yeast and human. Translational genetics have been employed to examine the functional roles of mammalian proteins using yeast as a pliable model system. In the current review recent advancements made in this area are discussed, highlighting work which shows that the cellular functions of human proteins in DNA repair and maintenance of genomic stability can be elucidated by genetic rescue experiments performed in yeast.

Targeting p53 as a therapeutic strategy in sensitizing TRAIL-induced apoptosis in cancer cells

Tumor necrosis factor (TNF)-related apoptosis-inducing ligand (TRAIL) has been intensively studied as a cancer therapeutic agent due to its unique ability to induce apoptosis in malignant cells but not in normal cells. However, as more human cancer cells are reported to be resistant to TRAIL treatment, it is important to develop new therapeutic strategies to overcome this resistance. p53 is an important tumor suppressor that is widely involved in cellular responses to various stresses. In this mini-review, we aim to provide an overview of the intricate relationship between p53 and the TRAIL-mediated apoptosis pathway, and to summarize the current approaches of targeting p53 as a therapeutic strategy to sensitize TRAIL-induced apoptosis in human cancer cells. Although in some cases TRAIL kills cancer cells in a p53-independent manner, it is believed that in cancers with wild-type and functional p53, targeting p53 may be an important strategy for overcoming TRAIL-resistance in cancer therapy.

p53 transactivation and the impact of mutations, cofactors and small molecules using a simplified yeast-based screening system

BACKGROUND

The p53 tumor suppressor, which is altered in most cancers, is a sequence-specific transcription factor that is able to modulate the expression of many target genes and influence a variety of cellular pathways. Inactivation of the p53 pathway in cancer frequently occurs through the expression of mutant p53 protein. In tumors that retain wild type p53, the pathway can be altered by upstream modulators, particularly the p53 negative regulators MDM2 and MDM4.

METHODOLOGY/PRINCIPAL FINDINGS

Given the many factors that might influence p53 function, including expression levels, mutations, cofactor proteins and small molecules, we expanded our previously described yeast-based system to provide the opportunity for efficient investigation of their individual and combined impacts in a miniaturized format. The system integrates i) variable expression of p53 proteins under the finely tunable GAL1,10 promoter, ii) single copy, chromosomally located p53-responsive and control luminescence reporters, iii) enhanced chemical uptake using modified ABC-transporters, iv) small-volume formats for treatment and dual-luciferase assays, and v) opportunities to co-express p53 with other cofactor proteins. This robust system can distinguish different levels of expression of WT and mutant p53 as well as interactions with MDM2 or 53BP1.

CONCLUSIONS/SIGNIFICANCE

We found that the small molecules Nutlin and RITA could both relieve the MDM2-dependent inhibition of WT p53 transactivation function, while only RITA could impact p53/53BP1 functional interactions. PRIMA-1 was ineffective in modifying the transactivation capacity of WT p53 and missense p53 mutations. This dual-luciferase assay can, therefore, provide a high-throughput assessment tool for investigating a matrix of factors that can influence the p53 network, including the effectiveness of newly developed small molecules, on WT and tumor-associated p53 mutants as well as interacting proteins.

The Toll-like receptor gene family is integrated into human DNA damage and p53 networks

In recent years the functions that the p53 tumor suppressor plays in human biology have been greatly extended beyond "guardian of the genome." Our studies of promoter response element sequences targeted by the p53 master regulatory transcription factor suggest a general role for this DNA damage and stress-responsive regulator in the control of human Toll-like receptor (TLR) gene expression. The TLR gene family mediates innate immunity to a wide variety of pathogenic threats through recognition of conserved pathogen-associated molecular motifs. Using primary human immune cells, we have examined expression of the entire TLR gene family following exposure to anti-cancer agents that induce the p53 network. Expression of all TLR genes, TLR1 to TLR10, in blood lymphocytes and alveolar macrophages from healthy volunteers can be induced by DNA metabolic stressors. However, there is considerable inter-individual variability. Most of the TLR genes respond to p53 via canonical as well as noncanonical promoter binding sites. Importantly, the integration of the TLR gene family into the p53 network is unique to primates, a recurrent theme raised for other gene families in our previous studies. Furthermore, a polymorphism in a TLR8 response element provides the first human example of a p53 target sequence specifically responsible for endogenous gene induction. These findings-demonstrating that the human innate immune system, including downstream induction of cytokines, can be modulated by DNA metabolic stress-have many implications for health and disease, as well as for understanding the evolution of damage and p53 responsive networks.

p53 Superfamily proteins in marine bivalve cancer and stress biology

The human p53 tumour suppressor protein is inactivated in many cancers and is also a major player in apoptotic responses to cellular stress. The p53 protein and the two other members of this protein family (p63, p73) are encoded by distinct genes and their functions have been extensively documented for humans and some other vertebrates. The structure and relative expression levels for members of the p53 superfamily have also been reported for most major invertebrate taxa. The functions of homologous proteins have been investigated for only a few invertebrates (specifically, p53 in flies, nematodes and recently a sea anemone). These studies of classical model organisms all suggest that the gene family originally evolved to mediate apoptosis of damaged germ cells or to protect germ cells from genotoxic stress. Here, we have correlated data from a number of molluscan and other invertebrate sequencing projects to provide a framework for understanding p53 signalling pathways in marine bivalve cancer and stress biology. These data suggest that (a) the two identified p53 and p63/73-like proteins in soft shell clam (Mya arenaria), blue mussel (Mytilus edulis) and Northern European squid (Loligo forbesi) have identical core sequences and may be splice variants of a single gene, while some molluscs and most other invertebrates have two or more distinct genes expressing different p53 family members; (b) transcriptional activation domains (TADs) in bivalve p53 and p63/73-like protein sequences are 67-69% conserved with human p53, while those in ecdysozoan, cnidarian, placozoan and choanozoan eukaryotes are ≤33% conserved; (c) the Mdm2 binding site in the transcriptional activation domain is 100% conserved in all sequenced bivalve p53 proteins (e.g. Mya, Mytilus, Crassostrea and Spisula) but is not present in other non-deuterostome invertebrates; (d) an Mdm2 homologue has been cloned for Mytilus trossulus; (e) homologues for both human p53 upstream regulatory and transcriptional target genes exist in molluscan genomes (missing are ARF, CIP1 and BH3 only proteins) and (f) p53 is demonstrably involved in bivalve haemocyte and germinoma cancers. We usually do not know enough about the molecular biology of marine invertebrates to address molecular mechanisms that characterize particular diseases. Understanding the molecular basis of naturally occurring diseases in marine bivalves is a virtually unexplored aspect of toxicoproteomics and genomics and related drug discovery. Additionally, increases in coastal development and concomitant increases in aquatic pollutants have driven interest in developing models appropriate for evaluating potential hazardous compounds or conditions found in the aquatic environment. Data reviewed in this study are coupled with recent developments in our understanding the molecular biology of the marine bivalve p53 superfamily. Taken together, they suggest that both structurally and functionally, bivalve p53 family proteins are the most highly conserved members of this gene superfamily so far identified outside of higher vertebrates and invertebrate chordates. Marine bivalves provide some of the most relevant and best understood models currently available for experimental studies by biomedical and marine environmental researchers.

The coordinated p53 and estrogen receptor cis-regulation at an FLT1 promoter SNP is specific to genotoxic stress and estrogenic compound

BACKGROUND

Recently, we established that a C>T single nucleotide polymorphism (SNP) in the promoter of the VEGF receptor FLT1 gene generates a (1/2) site p53 response element (RE-T) that results in p53 responsiveness of the promoter. The transcriptional control required an estrogen receptor (ER) (1/2) site response element (ERE1) 225 nt upstream to the RE-T.

METHODOLOGY/PRINCIPAL FINDINGS

Here we report the identification of a second ER (1/2) site (ERE2) located 145 bp downstream of the RE-T and establish that both EREs can impact p53-mediated transactivation of FLT1-T in a manner that is cell type and ER level dependent. Gene reporter assays and ChIP experiments conducted in the breast cancer-derived MCF7 cells revealed that the ERE2 site was sufficient for p53-mediated ERalpha recruitment and transactivation of the FLT1-T promoter/reporter construct. Surprisingly, unlike the case for other p53 target promoters, p53-mediated transactivation of FLT1-T constructs or expression of the endogenous FLT1 gene, as well as binding of p53 and ER at the promoter constructs, was inducible by doxorubicin but not by 5-fluorouracil. Furthermore, ER activity at FLT1-T was differentially affected by ER ligands, compared to a control TFF1/pS2 ER target promoter. The p53-related transcription factors (TFs) p73 and p63 had no effect on FLT1 transactivation.

CONCLUSIONS/SIGNIFICANCE

We establish a new dimension to the p53 master regulatory network where p53-mediated transcription from a (1/2) site RE can be determined by ER binding at one or more cis-acting EREs in manner that is dependent on level of ER protein, the type of ER ligand and the specific p53-inducing agent.

Foxa1 functions as a pioneer transcription factor at transposable elements to activate Afp during differentiation of embryonic stem cells

Epigenetic control of genes that are silent in embryonic stem cells, but destined for expression during differentiation, includes distinctive hallmarks, such as simultaneous activating/repressing (bivalent) modifications of chromatin and DNA hypomethylation at enhancers of gene expression. Although alpha-fetoprotein (Afp) falls into this class of genes, as it is silent in pluripotent stem cells and activated during differentiation of endoderm, we find that Afp chromatin lacks bivalent histone modifications. However, critical regulatory sites for Afp activation, overlapping Foxa1/p53/Smad-binding elements, are located within a 300-bp region lacking DNA methylation, due to transposed elements underrepresented in CpG sequences: a short interspersed transposable element and a medium reiterated sequence 1 element. Forkhead family member Foxa1 is activated by retinoic acid treatment of embryonic stem cells, binds its DNA consensus site within the short interspersed transposable/medium reiterated sequence 1 elements, and displaces linker histone H1 from silent Afp chromatin. Small interfering RNA depletion of Foxa1 showed that Foxa1 is essential in providing chromatin access to transforming growth factor beta-activated Smad2 and Smad4 and their subsequent DNA binding. Together these transcription factors establish highly acetylated chromatin and promote expression of Afp. Foxa1 acts as a pioneer transcription factor in de novo activation of Afp, by exploiting a lack of methylation at juxtaposed transposed elements, to bind and poise chromatin for intersection with transforming growth factor beta signaling during differentiation of embryonic stem cells.

Estrogen receptor acting in cis enhances WT and mutant p53 transactivation at canonical and noncanonical p53 target sequences

p53 is a master regulatory, sequence-specific transcription factor that directly controls expression of over 100 genes in response to various stress signals. Transactivation is generally considered to occur through p53 binding to a consensus response element (RE) composed of two 5'-RRRCWWGYYY-3' decamers. Recently, studying the human angiogenesis-related gene FLT1 we discovered that p53 can mediate limited transactivation at a noncanonical 1/2 site and could synergize with the estrogen receptor (ER) acting in cis at a nearby ER 1/2 site. To address the generality of concerted transactivation by p53 and ER, the 1/2 site in the FLT1 promoter was replaced with a variety of 1/2 sites, as well as canonical weak and strong p53 REs of human target genes. The p53 transactivation of all tested sequences was greatly enhanced by ligand-activated ER acting in cis. Furthermore, enhanced transactivation extends to several cancer-associated p53 mutants with altered function, suggesting ER-dependent mutant p53 activity for at least some REs. The enhanced transactivation was also found with p63 and p73. We propose a general synergistic relationship between p53 family and ER master regulators in transactivation of p53 target canonical and noncanonical REs, which might be poorly responsive to p53 on their own. This relationship greatly expands the transcriptional master network regulated by p53 in terms of genes affected and levels of expression and has implications for the appearance and possible treatments of cancer.

Conservation of DNA-binding specificity and oligomerisation properties within the p53 family

Background

Transcription factors activate their target genes by binding to specific response elements. Many transcription factor families evolved from a common ancestor by gene duplication and subsequent divergent evolution. Members of the p53 family, which play key roles in cell-cycle control and development, share conserved DNA binding and oligomerisation domains but exhibit distinct functions. In this study, the molecular basis of the functional divergence of related transcription factors was investigated.

Results

We characterised the DNA-binding specificity and oligomerisation properties of human p53, p63 and p73, as well as p53 from other organisms using novel biophysical approaches. All p53 family members bound DNA cooperatively as tetramers with high affinity. Despite structural differences in the oligomerisation domain, the dissociation constants of the tetramers was in the low nanomolar range for all family members, indicating that the strength of tetramerisation was evolutionarily conserved. However, small differences in the oligomerisation properties were observed, which may play a regulatory role. Intriguingly, the DNA-binding specificity of p53 family members was highly conserved even for evolutionarily distant species. Additionally, DNA recognition was only weakly affected by CpG methylation. Prediction of p53/p63/p73 binding sites in the genome showed almost complete overlap between the different homologs.

Conclusion

Diversity of biological function of p53 family members is not reflected in differences in sequence-specific DNA binding. Hence, additional specificity factors must exist, which allowed the acquisition of novel functions during evolution while preserving original roles.

The expanding universe of p53 targets

The p53 tumour suppressor is modified through mutation or changes in expression in most cancers, leading to the altered regulation of hundreds of genes that are directly influenced by this sequence-specific transcription factor. Central to the p53 master regulatory network are the target response element (RE) sequences. The extent of p53 transactivation and transcriptional repression is influenced by many factors, including p53 levels, cofactors and the specific RE sequences, all of which contribute to the role that p53 has in the aetiology of cancer. This Review describes the identification and functionality of REs and highlights the inclusion of non-canonical REs that expand the universe of genes and regulation mechanisms in the p53 tumour suppressor network.

A regulatory loop composed of RAP80-HDM2-p53 provides RAP80-enhanced p53 degradation by HDM2 in response to DNA damage

The ubiquitin interaction motif-containing protein RAP80 plays a key role in DNA damage response signaling. Using genomic and functional analysis, we established that the expression of the RAP80 gene is regulated in a DNA damage-responsive manner by the master regulator p53. This regulation occurs at the transcriptional level through a noncanonical p53 response element in the RAP80 promoter. Although it is inducible by p53, RAP80 is also able to regulate p53 through an association with both p53 and the E3 ubiquitin ligase HDM2, providing HDM2-dependent enhancement of p53 polyubiquitination. Depletion of RAP80 by small interfering RNA stabilizes p53, which, following DNA damage, results in an increased transactivation of several p53 target genes as well as greater apoptosis. Consistent with these observations, exogenous expression of RAP80 selectively inhibits p53-dependent transactivation of target genes in an mdm2-dependent manner in MEF cells. Thus, we identify a new DNA damage-associated role for RAP80. It can function in an autoregulatory loop consisting of RAP80, HDM2, and the p53 master regulatory network, implying an important role for this loop in genome stability and oncogenesis.