p63 consensus DNA-binding site: identification, analysis and application into a p63MH algorithm

SERPINE1: Role in Cholangiocarcinoma Progression and a Therapeutic Target in the Desmoplastic Microenvironment

A heterogenous population of inflammatory elements, other immune and nonimmune cells and cancer-associated fibroblasts (CAFs) are evident in solid malignancies where they coexist with the growing tumor mass. In highly desmoplastic malignancies, CAFs are the prominent mesenchymal cell type in the tumor microenvironment (TME), where their presence and abundance signal a poor prognosis. CAFs play a major role in the progression of various cancers by remodeling the supporting stroma into a dense, fibrotic matrix while secreting factors that promote the maintenance of cancer stem-like characteristics, tumor cell survival, aggressive growth and metastasis and reduced sensitivity to chemotherapeutics. Tumors with high stromal fibrotic signatures are more likely to be associated with drug resistance and eventual relapse. Identifying the molecular underpinnings for such multidirectional crosstalk among the various normal and neoplastic cell types in the TME may provide new targets and novel opportunities for therapeutic intervention. This review highlights recent concepts regarding the complexity of CAF biology in cholangiocarcinoma, a highly desmoplastic cancer. The discussion focuses on CAF heterogeneity, functionality in drug resistance, contributions to a progressively fibrotic tumor stroma, the involved signaling pathways and the participating genes.

Transcription Factor ATF3 Participates in DeltaNp63-Mediated Proliferation of Corneal Epithelial Cells

Understanding the regulatory mechanisms underlying corneal epithelial cell (CEC) proliferation in vitro may provide the means to boost CEC production in cell therapy for ocular disorders. The transcription factor ΔNp63 plays a crucial role in the proliferation of CECs, but the underlying mechanisms is yet to be elucidated. TP63 and ΔNp63 are encoded by the TP63 gene via alternative promoters. We previously reported that both ΔNp63 and activating transcription factor (ATF3) are substantially expressed in cultured CECs, but the regulatory relationship between ΔNp63 and ATF3 is unknown. In the present study, we found that ΔNp63 increased ATF3 expression and ATF3 promoter activity in cultured CECs. The deletion of the p63 binding core site reduced ATF3 promoter activity. CECs overexpressing ATF3 exhibited significantly greater proliferation than control CECs. ATF3 knockdown suppressed the ΔNp63-induced increase in cell proliferation. Overexpression of ATF3 in CECs significantly elevated protein and mRNA levels of cyclin D. The protein levels of keratin 3/14, integrin β1, and involucrin did not differ between ATF3-overexpressing CECs, ATF3-downregulated CECs, and control cells. In conclusion, our results suggest that ΔNp63 increases CEC proliferation via the ΔNp63/ATF3/CDK pathway.

A luminal intermediate cell state maintains long-term prostate homeostasis and contributes to tumorigenesis

Cellular heterogeneity poses tremendous challenges for developing cell-targeted therapies and biomarkers of clinically significant prostate cancer. The origins of this heterogeneity within normal adult and aging tissue remain unknown, leaving cellular states and transcriptional programs that allow expansions of malignant clones unidentified. To define cell states that contribute to early cancer development, we performed clonal analyses and single cell transcriptomics of normal prostate from genetically-engineered mouse models. We uncovered a luminal transcriptional state with a unique "basal-like" Wnt/p63 signaling ( luminal intermediate , LumI) which contributes to the maintenance of long-term prostate homeostasis. Moreover, LumI cells greatly expand during early stages of tumorigenesis in several mouse models of prostate cancer. Genetic ablation of p63 in vivo in luminal cells reduced the formation of aggressive clones in mouse prostate tumor models. Finally, the LumI cells and Wnt signaling appear to significantly increase in human aging prostate and prostate cancer samples, highlighting the importance of this hybrid cell state for human pathologies with potential translational impact.

Overview of frequent pattern mining

Various methods of frequent pattern mining have been applied to genetic problems, specifically, to the combined association of two genotypes (a genotype pattern, or diplotype) at different DNA variants with disease. These methods have the ability to come up with a selection of genotype patterns that are more common in affected than unaffected individuals, and the assessment of statistical significance for these selected patterns poses some unique problems, which are briefly outlined here.

Mutual exclusivity of ESR1 and TP53 mutations in endocrine resistant metastatic breast cancer

Both TP53 and ESR1 mutations occur frequently in estrogen receptor positive (ER+) metastatic breast cancers (MBC) and their distinct roles in breast cancer tumorigenesis and progression are well appreciated. Recent clinical studies discovered mutual exclusivity between TP53 and ESR1 mutations in metastatic breast cancers; however, mechanisms underlying this intriguing clinical observation remain largely understudied and unknown. Here, we explored the interplay between TP53 and ESR1 mutations using publicly available clinical and experimental data sets. We first confirmed the robust mutational exclusivity using six independent cohorts with 1,056 ER+ MBC samples and found that the exclusivity broadly applies to all ER+ breast tumors regardless of their clinical and distinct mutational features. ESR1 mutant tumors do not exhibit differential p53 pathway activity, whereas we identified attenuated ER activity and expression in TP53 mutant tumors, driven by a p53-associated E2 response gene signature. Further, 81% of these p53-associated E2 response genes are either direct targets of wild-type (WT) p53-regulated transactivation or are mutant p53-associated microRNAs, representing bimodal mechanisms of ER suppression. Lastly, we analyzed the very rare cases with co-occurrences of TP53 and ESR1 mutations and found that their simultaneous presence was also associated with reduced ER activity. In addition, tumors with dual mutations showed higher levels of total and PD-L1 positive macrophages. In summary, our study utilized multiple publicly available sources to explore the mechanism underlying the mutual exclusivity between ESR1 and TP53 mutations, providing further insights and testable hypotheses of the molecular interplay between these two pivotal genes in ER+ MBC.

p63, a key regulator of Ago2, links to the microRNA-144 cluster

Abstract

As a key component of the RNA-induced silencing complex (RISC), Argonaute2 (Ago2) exhibits a dual function regulatory role in tumor progression. However, the mechanistic basis of differential regulation remains elusive. p63 is a homolog of the tumor suppressor p53. p63 isoforms play a critical role in tumorigenesis and metastasis. Herein, we show that p63 isoforms physically interact with and stabilize Ago2. Expression of p63 isoforms increases the levels of Ago2 protein, while depletion of p63 isoforms by shRNA decreases Ago2 protein levels. p63 strongly guides Ago2 dual functions in vitro and in vivo. Ectopic expression of the miR-144/451 cluster increases p63 protein levels; TAp63 transactivates the miR-144/451 cluster, forming a positive feedback loop. Notably, miR-144 activates p63 by directly targeting Itch, an E3 ligase of p63. Ectopic expression of miR-144 induces apoptosis in H1299 cells. miR-144 enhances TAp63 tumor suppressor function and inhibits cell invasion. Our findings uncover a novel function of p63 linking the miRNA-144 cluster and the Ago2 pathway.

Facts and questions

Identification of Ago2 as a p63 target.

Ago2 exhibits a dual function regulatory role in tumor progression; however, the molecular mechanism of Ago2 regulation remains unknown.

p63 strongly guides Ago2 dual functions in vitro and in vivo.

Unraveling a novel function of p63 links the miRNA-144 cluster and the Ago2 pathway.

Evolutionary history of the p53 family DNA-binding domain: insights from an Alvinella pompejana homolog

The extremophile Alvinella pompejana, an annelid worm living on the edge of hydrothermal vents in the Pacific Ocean, is an excellent model system for studying factors that govern protein stability. Low intrinsic stability is a crucial factor for the susceptibility of the transcription factor p53 to inactivating mutations in human cancer. Understanding its molecular basis may facilitate the design of novel therapeutic strategies targeting mutant p53. By analyzing expressed sequence tag (EST) data, we discovered a p53 family gene in A. pompejana. Protein crystallography and biophysical studies showed that it has a p53/p63-like DNA-binding domain (DBD) that is more thermostable than all vertebrate p53 DBDs tested so far, but not as stable as that of human p63. We also identified features associated with its increased thermostability. In addition, the A. pompejana homolog shares DNA-binding properties with human p53 family DBDs, despite its evolutionary distance, consistent with a potential role in maintaining genome integrity. Through extensive structural and phylogenetic analyses, we could further trace key evolutionary events that shaped the structure, stability, and function of the p53 family DBD over time, leading to a potent but vulnerable tumor suppressor in humans.

YAP1 activation by human papillomavirus E7 promotes basal cell identity in squamous epithelia

Persistent human papillomavirus (HPV) infection of stratified squamous epithelial cells causes nearly 5% of cancer cases worldwide. HPV-positive oropharyngeal cancers harbor few mutations in the Hippo signaling pathway compared to HPV-negative cancers at the same anatomical site, prompting the hypothesis that an HPV-encoded protein inactivates the Hippo pathway and activates the Hippo effector yes-associated protein (YAP1). The HPV E7 oncoprotein is required for HPV infection and for HPV-mediated oncogenic transformation. We investigated the effects of HPV oncoproteins on YAP1 and found that E7 activates YAP1, promoting YAP1 nuclear localization in basal epithelial cells. YAP1 activation by HPV E7 required that E7 binds and degrades the tumor suppressor protein tyrosine phosphatase non-receptor type 14 (PTPN14). E7 required YAP1 transcriptional activity to extend the lifespan of primary keratinocytes, indicating that YAP1 activation contributes to E7 carcinogenic activity. Maintaining infection in basal cells is critical for HPV persistence, and here we demonstrate that YAP1 activation causes HPV E7 expressing cells to be retained in the basal compartment of stratified epithelia. We propose that YAP1 activation resulting from PTPN14 inactivation is an essential, targetable activity of the HPV E7 oncoprotein relevant to HPV infection and carcinogenesis.

The ‘epithelial’ cells that cover our bodies are in a constant state of turnover. Every few weeks, the outermost layers die and are replaced by new cells from the layers below. For scientists, this raises a difficult question. Cells infected by human papillomaviruses, often known as HPV, can become cancerous over years or even decades. How do infected cells survive while the healthy cells around them mature and get replaced?

One clue could lie in PTPN14, a human protein which many papillomaviruses eliminate using their viral E7 protein; this mechanism could be essential for the virus to replicate and cause cancer. To find out the impact of losing PTPN14, Hatterschide et al. used human epithelial cells to make three-dimensional models of infected tissues. These experiments showed that, when papillomaviruses destroy PTPN14, a human protein called YAP1 turns on in the lowest, most long-lived layer of the tissue. Cells in which YAP1 is activated survive while those that carry the inactive version mature and die. This suggests that papillomaviruses turn on YAP1 to remain in tissues for long periods.

Papillomaviruses cause about five percent of all human cancers. Finding ways to stop them from activating YAP1 has the potential to prevent disease. Overall, the research by Hatterschide et al. also sheds light on other epithelial cancers which are not caused by viruses.

Genomic gain of RRS1 promotes hepatocellular carcinoma through reducing the RPL11-MDM2-p53 signaling

Hepatocellular carcinomas (HCCs) are characterized by frequent somatic genomic copy number alterations (CNAs), with most of them biologically unexplored. Here, we performed integrative analyses combining CNAs with the transcriptomic data to reveal the cis- and trans-effects of CNAs in HCC. We identified recurrent genomic gains of chromosome 8q, which exhibit strong trans-effects and are broadly associated with ribosome biogenesis activity. Furthermore, 8q gain-driven overexpression of ribosome biogenesis regulator (RRS1) promotes growth of HCC cells in vitro and in vivo. Mechanistically, RRS1 attenuates ribosomal stress through retaining RPL11 in the nucleolus, which, in turn, potentiates MDM2-mediated ubiquitination and degradation of p53. Clinically, higher RRS1 expression levels predict poor clinical outcomes for patients with HCC, especially in those with intact p53 Our findings established that the chromosome 8q oncogene RRS1 promotes HCC development through attenuating the RPL11-MDM2-p53 pathway and provided new potential targets for treatment of this malignancy.

p63 and p53: Collaborative Partners or Dueling Rivals?

The tumor suppressor p53 and its oncogenic sibling p63 (ΔNp63) direct opposing fates in tumor development. These paralog proteins are transcription factors that elicit their tumor suppressive and oncogenic capacity through the regulation of both shared and unique target genes. Both proteins predominantly function as activators of transcription, leading to a paradigm shift away from ΔNp63 as a dominant negative to p53 activity. The discovery of p53 and p63 as pioneer transcription factors regulating chromatin structure revealed new insights into how these paralogs can both positively and negatively influence each other to direct cell fate. The previous view of a strict rivalry between the siblings needs to be revisited, as p53 and p63 can also work together toward a common goal.

Hsp70 acts as a fine-switch that controls E3 ligase CHIP-mediated TAp63 and ΔNp63 ubiquitination and degradation

The major clinical problem in human cancer is metastasis. Metastases are the cause of 90% of human cancer deaths. TAp63 is a critical suppressor of tumorigenesis and metastasis. ΔNp63 acts as a dominant-negative inhibitor to block the function of p53 and TAp63. Although several ubiquitin E3 ligases have been reported to regulate p63 stability, the mechanism of p63 regulation remains partially understood. Herein, we show that CHIP, an E3 ligase with a U-box domain, physically interacts with p63 and promotes p63 degradation. Notably, Hsp70 depletion by siRNA stabilizes TAp63 in H1299 cells and destabilizes ΔNp63 in SCC9 cells. Loss of Hsp70 results in a reduction in the TAp63-CHIP interaction in H1299 cells and an increase in the interaction between ΔNp63 and CHIP in SCC9 cells. Our results reveal that Hsp70 acts as a molecular switch to control CHIP-mediated ubiquitination and degradation of p63 isoforms. Furthermore, regulation of p63 by the Hsp70-CHIP axis contributes to the migration and invasion of tumor cells. Hence, our findings demonstrate that Hsp70 is a crucial regulator of CHIP-mediated ubiquitination and degradation of p63 isoforms and identify a new pathway for maintaining TAp63 or ΔNp63 stability in cancers.

Loss-of-Function Variants in the Tumor-Suppressor Gene PTPN14 Confer Increased Cancer Risk

The success of genome-wide association studies (GWAS) in identifying common, low-penetrance variant-cancer associations for the past decade is undisputed. However, discovering additional high-penetrance cancer mutations in unknown cancer predisposing genes requires detection of variant-cancer association of ultra-rare coding variants. Consequently, large-scale next-generation sequence data with associated phenotype information are needed. Here, we used genotype data on 166,281 Icelanders, of which, 49,708 were whole-genome sequenced and 408,595 individuals from the UK Biobank, of which, 41,147 were whole-exome sequenced, to test for association between loss-of-function burden in autosomal genes and basal cell carcinoma (BCC), the most common cancer in Caucasians. A total of 25,205 BCC cases and 683,058 controls were tested. Rare germline loss-of-function variants in PTPN14 conferred substantial risks of BCC (OR, 8.0; P = 1.9 × 10^-12), with a quarter of carriers getting BCC before age 70 and over half in their lifetime. Furthermore, common variants at the PTPN14 locus were associated with BCC, suggesting PTPN14 as a new, high-impact BCC predisposition gene. A follow-up investigation of 24 cancers and three benign tumor types showed that PTPN14 loss-of-function variants are associated with high risk of cervical cancer (OR, 12.7, P = 1.6 × 10^-4) and low age at diagnosis. Our findings, using power-increasing methods with high-quality rare variant genotypes, highlight future prospects for new discoveries on carcinogenesis. SIGNIFICANCE: This study identifies the tumor-suppressor gene PTPN14 as a high-impact BCC predisposition gene and indicates that inactivation of PTPN14 by germline sequence variants may also lead to increased risk of cervical cancer.

Evaluating the Influence of a G-Quadruplex Prone Sequence on the Transactivation Potential by Wild-Type and/or Mutant P53 Family Proteins through a Yeast-Based Functional Assay

P53, P63, and P73 proteins belong to the P53 family of transcription factors, sharing a common gene organization that, from the P1 and P2 promoters, produces two groups of mRNAs encoding proteins with different N-terminal regions; moreover, alternative splicing events at C-terminus further contribute to the generation of multiple isoforms. P53 family proteins can influence a plethora of cellular pathways mainly through the direct binding to specific DNA sequences known as response elements (REs), and the transactivation of the corresponding target genes. However, the transcriptional activation by P53 family members can be regulated at multiple levels, including the DNA topology at responsive promoters. Here, by using a yeast-based functional assay, we evaluated the influence that a G-quadruplex (G4) prone sequence adjacent to the p53 RE derived from the apoptotic PUMA target gene can exert on the transactivation potential of full-length and N-terminal truncated P53 family α isoforms (wild-type and mutant). Our results show that the presence of a G4 prone sequence upstream or downstream of the P53 RE leads to significant changes in the relative activity of P53 family proteins, emphasizing the potential role of structural DNA features as modifiers of P53 family functions at target promoter sites.

G9a controls pluripotent-like identity and tumor-initiating function in human colorectal cancer

Colorectal tumors are hierarchically organized and governed by populations of self-renewing cancer stem cells, representing one of the deadliest types of cancers worldwide. Emergence of cancer stemness phenotype depends on epigenetic reprogramming, associated with profound transcriptional changes. As described for pluripotent reprogramming, epigenetic modifiers play a key role in cancer stem cells by establishing embryonic stem-like transcriptional programs, thus impacting the balance between self-renewal and differentiation. We identified overexpression of histone methyltransferase G9a as a risk factor for colorectal cancer, associated with shorter relapse-free survival. Moreover, using human transformed pluripotent cells as a surrogate model for cancer stem cells, we observed that G9a activity is essential for the maintenance of embryonic-like transcriptional signature promoting self-renewal, tumorigenicity, and undifferentiated state. Such a role was also applicable to colorectal cancer, where inhibitors of G9a histone methyltransferase function induced intestinal differentiation while restricting tumor-initiating activity in patient-derived colorectal tumor samples. Finally, by integrating transcriptome profiling with G9a/H3K9me2 loci co-occupancy, we identified the canonical Wnt pathway, epithelial-to-mesenchyme transition, and extracellular matrix organization as potential targets of such a chromatin regulation mechanism in colorectal cancer stem cells. Overall, our findings provide novel insights on the role of G9a as a driver of cancer stem cell phenotype, promoting self-renewal, tumorigenicity, and undifferentiated state.

Dissecting the DNA binding landscape and gene regulatory network of p63 and p53

The transcription factor p53 is the best-known tumor suppressor, but its sibling p63 is a master regulator of epidermis development and a key oncogenic driver in squamous cell carcinomas (SCC). Despite multiple gene expression studies becoming available, the limited overlap of reported p63-dependent genes has made it difficult to decipher the p63 gene regulatory network. Particularly, analyses of p63 response elements differed substantially among the studies. To address this intricate data situation, we provide an integrated resource that enables assessing the p63-dependent regulation of any human gene of interest. We use a novel iterative de novo motif search approach in conjunction with extensive ChIP-seq data to achieve a precise global distinction between p53-and p63-binding sites, recognition motifs, and potential co-factors. We integrate these data with enhancer:gene associations to predict p63 target genes and identify those that are commonly de-regulated in SCC representing candidates for prognosis and therapeutic interventions.

P63 modulates the expression of the WDFY2 gene which is implicated in cancer regulation and limb development

TP63 is a member of the TP53 gene family, sharing a common gene structure that produces two groups of mRNAs' encoding proteins with different N-terminal regions (ΔN and TA isoforms); both transcripts are also subjected to alternative splicing mechanisms at C-terminus, generating a variety of isoforms. p63 is a master regulator of epidermal development and homoeostasis as well as an important player in tumorigenesis and cancer progression with both oncogenic and tumour suppressive roles. A number of studies have aimed at the identification of p63 target genes, allowing the dissection of the molecular pathways orchestrated by the different isoforms. In the present study we investigated in more detail the p63 responsiveness of the WDFY2 (WD repeat and FYVE domain containing 2) gene, encoding for an endosomal protein identified as a binding partner of the PI-3K/AKT signalling pathway. We showed that overexpression of different p63 isoforms was able to induce WDFY2 expression in TP53-null cells. The p63-dependent transcriptional activation was associated with specific response elements (REs) that have been identified by a bioinformatics tool and validated by yeast- and mammal-based assays. Interestingly, to confirm that WDFY2 belongs to the p63 network of cancer regulation, we analysed the impact of WDFY2 alterations, by showing its frequent deletion in different types of tumours and suggesting its expression level as a prognostic biomarker. Lastly, we identified a chromosomal translocation involving the WDFY2 locus in a patient affected by a rare congenital limb anomaly, indicating WDFY2 as a possible susceptibility gene placed downstream p63 in the network of limb development.

ΔNp63α Suppresses TGFB2 Expression and RHOA Activity to Drive Cell Proliferation in Squamous Cell Carcinomas

The transcriptional repressor ΔNp63α is a potent oncogene widely overexpressed in squamous cell carcinomas (SCCs) of diverse tissue origins, where it promotes malignant cell proliferation and survival. We report here the results of a genome-wide CRISPR screen to identify pathways controlling ΔNp63α-dependent cell proliferation, which revealed that the small GTPase RHOA blocks cell division upon ΔNp63α knockdown. After ΔNp63α depletion, RHOA activity is increased, and cells undergo RHOA-dependent proliferation arrest along with transcriptome changes indicative of increased TGF-β signaling. Mechanistically, ΔNp63α represses transcription of TGFB2, which induces a cell cycle arrest that is partially dependent on RHOA. Ectopic TGFB2 activates RHOA and impairs SCC proliferation, and TGFB2 neutralization restores cell proliferation during ΔNp63α depletion. Genomic data from tumors demonstrate inactivation of RHOA and the TGFBR2 receptor and ΔNp63α overexpression in more than 80% of lung SCCs. These results reveal a signaling pathway controlling SCC proliferation that is potentially amenable to pharmacological intervention.

α-Ketoglutarate links p53 to cell fate during tumour suppression

The tumour suppressor TP53 is mutated in the majority of human cancers, and in over 70% of pancreatic ductal adenocarcinoma (PDAC)^1,2. Wild-type p53 accumulates in response to cellular stress, and regulates gene expression to alter cell fate and prevent tumour development². Wild-type p53 is also known to modulate cellular metabolic pathways³, although p53-dependent metabolic alterations that constrain cancer progression remain poorly understood. Here we find that p53 remodels cancer-cell metabolism to enforce changes in chromatin and gene expression that favour a premalignant cell fate. Restoring p53 function in cancer cells derived from KRAS-mutant mouse models of PDAC leads to the accumulation of α-ketoglutarate (αKG, also known as 2-oxoglutarate), a metabolite that also serves as an obligate substrate for a subset of chromatin-modifying enzymes. p53 induces transcriptional programs that are characteristic of premalignant differentiation, and this effect can be partially recapitulated by the addition of cell-permeable αKG. Increased levels of the αKG-dependent chromatin modification 5-hydroxymethylcytosine (5hmC) accompany the tumour-cell differentiation that is triggered by p53, whereas decreased 5hmC characterizes the transition from premalignant to de-differentiated malignant lesions that is associated with mutations in Trp53. Enforcing the accumulation of αKG in p53-deficient PDAC cells through the inhibition of oxoglutarate dehydrogenase-an enzyme of the tricarboxylic acid cycle-specifically results in increased 5hmC, tumour-cell differentiation and decreased tumour-cell fitness. Conversely, increasing the intracellular levels of succinate (a competitive inhibitor of αKG-dependent dioxygenases) blunts p53-driven tumour suppression. These data suggest that αKG is an effector of p53-mediated tumour suppression, and that the accumulation of αKG in p53-deficient tumours can drive tumour-cell differentiation and antagonize malignant progression.

p63 and SOX2 Dictate Glucose Reliance and Metabolic Vulnerabilities in Squamous Cell Carcinomas

Squamous cell carcinoma (SCC), a malignancy arising across multiple anatomical sites, is responsible for significant cancer mortality due to insufficient therapeutic options. Here, we identify exceptional glucose reliance among SCCs dictated by hyperactive GLUT1-mediated glucose influx. Mechanistically, squamous lineage transcription factors p63 and SOX2 transactivate the intronic enhancer cluster of SLC2A1. Elevated glucose influx fuels generation of NADPH and GSH, thereby heightening the anti-oxidative capacity in SCC tumors. Systemic glucose restriction by ketogenic diet and inhibiting renal glucose reabsorption with SGLT2 inhibitor precipitate intratumoral oxidative stress and tumor growth inhibition. Furthermore, reduction of blood glucose lowers blood insulin levels, which suppresses PI3K/AKT signaling in SCC cells. Clinically, we demonstrate a robust correlation between blood glucose concentration and worse survival among SCC patients. Collectively, this study identifies the exceptional glucose reliance of SCC and suggests its candidacy as a highly vulnerable cancer type to be targeted by systemic glucose restriction.

Roles of p53 Family Structure and Function in Non-Canonical Response Element Binding and Activation

The p53 canonical consensus sequence is a 10-bp repeat of PuPuPuC(A/T)(A/T)GPyPyPy, separated by a spacer with up to 13 bases. C(A/T)(A/T)G is the core sequence and purine (Pu) and pyrimidine (Py) bases comprise the flanking sequence. However, in the p53 noncanonical sequences, there are many variations, such as length of consensus sequence, variance of core sequence or flanking sequence, and variance in number of bases making up the spacer or AT gap composition. In comparison to p53, the p53 family members p63 and p73 have been found to have more tolerance to bind and activate several of these noncanonical sequences. The p53 protein forms monomers, dimers, and tetramers, and its nonspecific binding domain is well-defined; however, those for p63 or p73 are still not fully understood. Study of p63 and p73 structure to determine the monomers, dimers or tetramers to bind and regulate noncanonical sequence is a new challenge which is crucial to obtaining a complete picture of structure and function in order to understand how p63 and p73 regulate genes differently from p53. In this review, we will summarize the rules of p53 family non-canonical sequences, especially focusing on the structure of p53 family members in the regulation of specific target genes. In addition, we will compare different software programs for prediction of p53 family responsive elements containing parameters with canonical or non-canonical sequences.

Molecular Mechanisms of p63-Mediated Squamous Cancer Pathogenesis

The p63 gene is a member of the p53/p63/p73 family of transcription factors and plays a critical role in development and homeostasis of squamous epithelium. p63 is transcribed as multiple isoforms; ΔNp63α, the predominant p63 isoform in stratified squamous epithelium, is localized to the basal cells and is overexpressed in squamous cell cancers of multiple organ sites, including skin, head and neck, and lung. Further, p63 is considered a stem cell marker, and within the epidermis, ΔNp63α directs lineage commitment. ΔNp63α has been implicated in numerous processes of skin biology that impact normal epidermal homeostasis and can contribute to squamous cancer pathogenesis by supporting proliferation and survival with roles in blocking terminal differentiation, apoptosis, and senescence, and influencing adhesion and migration. ΔNp63α overexpression may also influence the tissue microenvironment through remodeling of the extracellular matrix and vasculature, as well as by enhancing cytokine and chemokine secretion to recruit pro-inflammatory infiltrate. This review focuses on the role of ΔNp63α in normal epidermal biology and how dysregulation can contribute to cutaneous squamous cancer development, drawing from knowledge also gained by squamous cancers from other organ sites that share p63 overexpression as a defining feature.

Epigenetic evolution and lineage histories of chronic lymphocytic leukaemia

Genetic and epigenetic intra-tumoral heterogeneity cooperate to shape the evolutionary course of cancer¹. Chronic lymphocytic leukaemia (CLL) is a highly informative model for cancer evolution as it undergoes substantial genetic diversification and evolution after therapy^2,3. The CLL epigenome is also an important disease-defining feature^4,5, and growing populations of cells in CLL diversify by stochastic changes in DNA methylation known as epimutations⁶. However, previous studies using bulk sequencing methods to analyse the patterns of DNA methylation were unable to determine whether epimutations affect CLL populations homogeneously. Here, to measure the epimutation rate at single-cell resolution, we applied multiplexed single-cell reduced-representation bisulfite sequencing to B cells from healthy donors and patients with CLL. We observed that the common clonal origin of CLL results in a consistently increased epimutation rate, with low variability in the cell-to-cell epimutation rate. By contrast, variable epimutation rates across healthy B cells reflect diverse evolutionary ages across the trajectory of B cell differentiation, consistent with epimutations serving as a molecular clock. Heritable epimutation information allowed us to reconstruct lineages at high-resolution with single-cell data, and to apply this directly to patient samples. The CLL lineage tree shape revealed earlier branching and longer branch lengths than in normal B cells, reflecting rapid drift after the initial malignant transformation and a greater proliferative history. Integration of single-cell bisulfite sequencing analysis with single-cell transcriptomes and genotyping confirmed that genetic subclones mapped to distinct clades, as inferred solely on the basis of epimutation information. Finally, to examine potential lineage biases during therapy, we profiled serial samples during ibrutinib-associated lymphocytosis, and identified clades of cells that were preferentially expelled from the lymph node after treatment, marked by distinct transcriptional profiles. The single-cell integration of genetic, epigenetic and transcriptional information thus charts the lineage history of CLL and its evolution with therapy.

PTPN14 degradation by high-risk human papillomavirus E7 limits keratinocyte differentiation and contributes to HPV-mediated oncogenesis

High-risk human papillomavirus (HPV) E7 proteins enable oncogenic transformation of HPV-infected cells by inactivating host cellular proteins. High-risk but not low-risk HPV E7 target PTPN14 for proteolytic degradation, suggesting that PTPN14 degradation may be related to their oncogenic activity. HPV infects human keratinocytes but the role of PTPN14 in keratinocytes and the consequences of PTPN14 degradation are unknown. Using an HPV16 E7 variant that can inactivate retinoblastoma tumor suppressor (RB1) but cannot degrade PTPN14, we found that high-risk HPV E7-mediated PTPN14 degradation impairs keratinocyte differentiation. Deletion of PTPN14 from primary human keratinocytes decreased keratinocyte differentiation gene expression. Related to oncogenic transformation, both HPV16 E7-mediated PTPN14 degradation and PTPN14 deletion promoted keratinocyte survival following detachment from a substrate. PTPN14 degradation contributed to high-risk HPV E6/E7-mediated immortalization of primary keratinocytes and HPV⁺ but not HPV^- cancers exhibit a gene-expression signature consistent with PTPN14 inactivation. We find that PTPN14 degradation impairs keratinocyte differentiation and propose that this contributes to high-risk HPV E7-mediated oncogenic activity independent of RB1 inactivation.

IL-1β Inflammatory Cytokine-Induced TP63 Isoform ∆NP63α Signaling Cascade Contributes to Cisplatin Resistance in Human Breast Cancer Cells

The mechanisms behind the induction of malignancy and chemoresistance in breast cancer cells are still not completely understood. Inflammation is associated with the induction of malignancy in different types of cancer and is highlighted as an important factor for chemoresistance. In previous work, we demonstrated that the inflammatory cytokine interleukin 1β (IL-1β)-induced upregulation of genes was associated with chemoresistance in breast cancer cells. Here, we evaluated the participation and the expression profile of TP63 in the induction of resistance to cisplatin. By loss-of-function assays, we identified that IL-1β particularly upregulates the expression of the tumor protein 63 (TP63) isoform ΔNP63α, through the activation of the IL-1β/IL-1RI/β-catenin signaling pathway. Upregulation of ΔNP63α leads to an increase in the expression of the cell survival factors epidermal growth factor receptor (EGFR) and phosphatase 1D (Wip1), and a decrease in the DNA damage sensor, ataxia-telangiectasia mutated (ATM). The participation of these processes in the increase of resistance to cisplatin was confirmed by silencing TP63 expression or by inhibition of the phosphoinositide 3-kinase (PI3K)/protein kinase B (AKT) activity in the IL-1β/IL-1RI/β-catenin signaling pathway. These data reinforced the importance of an inflammatory environment in the induction of drug resistance in cancer cells and uncovered a molecular mechanism where the IL-1β signaling pathway potentiates the acquisition of cisplatin resistance in breast cancer cells.

ΔNp63 Inhibits Oxidative Stress-Induced Cell Death, Including Ferroptosis, and Cooperates with the BCL-2 Family to Promote Clonogenic Survival

The BCL-2 family proteins are central regulators of apoptosis. However, cells deficient for BAX and BAK or overexpressing BCL-2 still succumb to oxidative stress upon DNA damage or matrix detachment. Here, we show that ΔNp63α overexpression protects cells from oxidative stress induced by oxidants, DNA damage, anoikis, or ferroptosis-inducing agents. Conversely, ΔNp63α deficiency increases oxidative stress. Mechanistically, ΔNp63α orchestrates redox homeostasis through transcriptional control of glutathione biogenesis, utilization, and regeneration. Analysis of a lung squamous cell carcinoma dataset from The Cancer Genome Atlas (TCGA) reveals that TP63 amplification/overexpression upregulates the glutathione metabolism pathway in primary human tumors. Strikingly, overexpression of ΔNp63α promotes clonogenic survival of p53^-/-Bax^-/-Bak^-/- cells against DNA damage. Furthermore, co-expression of BCL-2 and ΔNp63α confers clonogenic survival against matrix detachment, disrupts the luminal clearance of mammary acini, and promotes cancer metastasis. Our findings highlight the need for a simultaneous blockade of apoptosis and oxidative stress to promote long-term cellular well-being.

Motifome comparison between modern human, Neanderthal and Denisovan

BACKGROUND

The availability of the genomes of two archaic humans, Neanderthal and Denisovan, and that of modern humans provides researchers an opportunity to investigate genetic differences between these three subspecies on a genome-wide scale. Here we describe an algorithm that predicts statistically significant motifs based on the difference between a given motif's actual and expected distributions. The algorithm was previously applied to plants but was modified for this work.

RESULTS

The result of applying the algorithm to the human, Neanderthal, and Denisovan genomes is a catalog of potential regulatory motifs in these three human subspecies. We examined the distributions of these motifs in genetic elements including human retroviruses, human accelerated regions, and human accelerated conserved noncoding sequences regions. Differences in these distributions could be the origin of differences in phenotype between the three subspecies. Twenty significant motifs common to all three genomes were found; thirty-three were found in endogenous retroviruses in Neanderthal and Denisovan. Ten of these motifs mapped to the 22 bp core of MiR-1304. The core of this genetic element regulates the ENAM and AMTN genes, which take part in odontogenesis and whose 3' UTRs contained significant motifs. The introns of 20 genes were found to contain a large number of significant motifs, which were also overrepresented in 49 human accelerated regions. These genes include NAV2, SorCS2, TRAPPC9, GRID1, PRDM16, CAMTA1, and ASIC which are all involved in neuroregulation. Further analysis of these genes using the GO database indicates that many are associated with neurodevelopment. Also, varying numbers of significant motifs were found to occur in regions of the Neanderthal and Denisovan genomes that are missing from the human genome, suggesting further functional differences between modern and archaic humans.

CONCLUSION

Although Neanderthal and Denisovan are now extinct, detailed examination of elements from their genomes can shed light on possible phenotypic and cognitive differences between these two archaic human subspecies and modern humans. Genetic similarities and differences between these three subspecies and other fossil hominids would also be of interest.

ΔNp63α down-regulates c-Myc modulator MM1 via E3 ligase HERC3 in the regulation of cell senescence

p63 and c-Myc are key transcription factors controlling genes involved in the cell cycle and cellular senescence. We previously reported that p63α can destabilize MM1 protein to derepress c-Myc, resulting in cell cycle progress and tumorigenesis. However, how the proteasomal degradation of MM1 is facilitated remains unclear. In the present study, we identified a novel E3 ligase, HERC3, which can mediate ubiquitination of MM1 and promote its proteasome-dependent degradation. We found that ΔNp63α transcriptionally up-regulates HERC3 and knockdown of HERC3 abrogates ΔNp63α-induced down-regulation of MM1. Either overexpression of MM1 or ablation of HERC3 induces cell senescence, while knockdown of MM1 rescues cell senescence induced by deficiency of either ΔNp63α or HERC3, implicating the involvement of the ΔNp63α/HERC3/MM1/c-Myc axis in the modulation of cell senescence. Additionally, our Oncomine analysis indicates activation of the ΔNp63α/HERC3/MM1/c-Myc axis in invasive breast carcinoma. Together, our data illuminate a novel axis regulating cell senescence: ΔNp63α stimulates transcription of E3 ligase HERC3, which mediates ubiquitination of c-Myc modulator MM1 and targets it to proteasomal degradation; subsequently, c-Myc is derepressed by ΔNp63α, thereby cell senescence is modulated by this axis. Our work provides a new interpretation of crosstalk between p63 and c-Myc, and also sheds new light on ΔNp63α-controlled cell senescence and tumorigenesis.

A double dealing tale of p63: an oncogene or a tumor suppressor

As a member of tumor suppressor p53 family, p63, a gene encoding versatile protein variant, has been documented to correlate with cancer formation and progression, though it is rarely mutated in cancer patients. However, it has long been controversial on whether p63 is an oncogene or a tumor suppressor. Here, we comprehensively reviewed reports on roles of p63 in development, tumorigenesis and tumor progression. According to data from molecular cell biology, genetic models and clinic research, we conclude that p63 may act as either an oncogene or a tumor suppressor gene in different scenarios: TA isoforms of p63 gene are generally tumor-suppressive through repressing cell proliferation, survival and metastasis; ΔN isoforms, however, may initiate tumorigenesis via promoting cell proliferation and survival, but inhibit tumor metastasis and progression; effects of p63 on tumor formation and progression depend on the context of the whole p53 family, and either amplification or loss of p63 gene locus can break the balance to cause tumorigenesis.

Tumor-Suppressor Inactivation of GDF11 Occurs by Precursor Sequestration in Triple-Negative Breast Cancer

Triple-negative breast cancer (TNBC) is an aggressive and heterogeneous carcinoma in which various tumor-suppressor genes are lost by mutation, deletion, or silencing. Here we report a tumor-suppressive mode of action for growth-differentiation factor 11 (GDF11) and an unusual mechanism of its inactivation in TNBC. GDF11 promotes an epithelial, anti-invasive phenotype in 3D triple-negative cultures and intraductal xenografts by sustaining expression of E-cadherin and inhibitor of differentiation 2 (ID2). Surprisingly, clinical TNBCs retain the GDF11 locus and expression of the protein itself. GDF11 bioactivity is instead lost because of deficiencies in its convertase, proprotein convertase subtilisin/kexin type 5 (PCSK5), causing inactive GDF11 precursor to accumulate intracellularly. PCSK5 reconstitution mobilizes the latent TNBC reservoir of GDF11 in vitro and suppresses triple-negative mammary cancer metastasis to the lung of syngeneic hosts. Intracellular GDF11 retention adds to the concept of tumor-suppressor inactivation and reveals a cell-biological vulnerability for TNBCs lacking therapeutically actionable mutations.

Copy number gains at chr3p25 and chr11p11 are associated with lymph node involvement and survival in muscle-invasive bladder tumors

Patients with muscle-invasive bladder cancer (MIBC) have poorer prognoses if cancer has metastasized to the lymph nodes. Genomic markers of lymph node involvement (LNI) would be useful for treatment planning, especially if measured at the biopsy stage, but large-scale studies of tumor tissue at any stage are needed to discover robust markers of LNI. We performed a genome-wide query of copy number alterations (CNA) in 237 MIBC surgical tumor specimens from patients in The Cancer Genome Atlas who had radical cystectomy and lymphadenectomy without neoadjuvant treatment. Pathology reports were independently reviewed to confirm LNI, and copy number data was analyzed to confirm gene-level gains and losses while adjusting for tumor purity and ploidy. Using logistic regression and elastic net models, we identified the CNA most significantly associated with LNI. Multivariable logistic regression was used to describe these CNA associations while adjusting for clinical variables. Kaplan-Meier and Cox regression were used to describe their association with overall survival. Gains in 26 genes were identified as having strong associations with LNI. After adjusting for age, gender, race, pathological tumor stage, histology, and number of nodes examined, gains in 22 genes on chr3p25 or chr11p11 remained significantly associated with LNI (p<0.01) and improved model discrimination over clinical variables alone (p = 0.04). They were also associated with shorter overall survival (adjusted p = 0.02). These results suggest that a simple genomic test for gains in chr3p25 and chr11p11 could inform adjuvant treatment or clinical trial decisions if validated in external cohorts. Additional studies will also be needed to determine if these CNA are detectible in biopsy tissue and can inform clinical decisions at the preoperative stage.

Evolutionary re-wiring of p63 and the epigenomic regulatory landscape in keratinocytes and its potential implications on species-specific gene expression and phenotypes

Although epidermal keratinocyte development and differentiation proceeds in similar fashion between humans and mice, evolutionary pressures have also wrought significant species-specific physiological differences. These differences between species could arise in part, by the rewiring of regulatory network due to changes in the global targets of lineage-specific transcriptional master regulators such as p63. Here we have performed a systematic and comparative analysis of the p63 target gene network within the integrated framework of the transcriptomic and epigenomic landscape of mouse and human keratinocytes. We determined that there exists a core set of ∼1600 genomic regions distributed among enhancers and super-enhancers, which are conserved and occupied by p63 in keratinocytes from both species. Notably, these DNA segments are typified by consensus p63 binding motifs under purifying selection and are associated with genes involved in key keratinocyte and skin-centric biological processes. However, the majority of the p63-bound mouse target regions consist of either murine-specific DNA elements that are not alignable to the human genome or exhibit no p63 binding in the orthologous syntenic regions, typifying an occupancy lost subset. Our results suggest that these evolutionarily divergent regions have undergone significant turnover of p63 binding sites and are associated with an underlying inactive and inaccessible chromatin state, indicative of their selective functional activity in the transcriptional regulatory network in mouse but not human. Furthermore, we demonstrate that this selective targeting of genes by p63 correlates with subtle, but measurable transcriptional differences in mouse and human keratinocytes that converges on major metabolic processes, which often exhibit species-specific trends. Collectively our study offers possible molecular explanation for the observable phenotypic differences between the mouse and human skin and broadly informs on the prevailing principles that govern the tug-of-war between evolutionary forces of rigidity and plasticity over transcriptional regulatory programs.

The Regulation of Tumor Suppressor p63 by the Ubiquitin-Proteasome System

The protein p63 has been identified as a homolog of the tumor suppressor protein p53 and is capable of inducing apoptosis, cell cycle arrest, or senescence. p63 has at least six isoforms, which can be divided into two major groups: the TAp63 variants that contain the N-terminal transactivation domain and the ΔNp63 variants that lack the N-terminal transactivation domain. The TAp63 variants are generally considered to be tumor suppressors involved in activating apoptosis and suppressing metastasis. ΔNp63 variants cannot induce apoptosis but can act as dominant negative inhibitors to block the function of TAp53, TAp73, and TAp63. p63 is rarely mutated in human tumors and is predominately regulated at the post-translational level by phosphorylation and ubiquitination. This review focuses primarily on regulation of p63 by the ubiquitin E-3 ligase family of enzymes via ubiquitination and proteasome-mediated degradation, and introduces a new key regulator of the p63 protein.

Neomorphic mutations create therapeutic challenges in cancer

Oncogenesis is a pathologic process driven by genomic aberrations, including changes in nucleotide sequences. The majority of these mutational events fall into two broad categories: inactivation of tumor suppressor genes (hypomorph, antimorph or amorph) or activation of oncogenes (hypermorph). The recent surge in genome sequence data and functional genomics research has ushered in the discovery of aberrations in a third category: gain-of-novel-function mutation (neomorph). These neomorphic mutations, which can be found in both tumor suppressor genes and oncogenes, produce proteins with entirely different functions from their respective wild-type (WT) proteins and the other morphs. The unanticipated phenotypic outcomes elicited by neomorphic mutations imply that tumors with the neomorphic mutations may not respond to therapies designed to target the WT protein. Therefore, understanding the functional activities of each genomic aberration to be targeted is crucial in devising effective treatment strategies that will benefit specific cancer patients.

Solution structure and binding specificity of the p63 DNA binding domain

p63 is a close homologue of p53 and, together with p73, is grouped into the p53 family of transcription factors. p63 is known to be involved in the induction of controlled apoptosis important for differentiation processes, germ line integrity and development. Despite its high homology to p53, especially within the DNA binding domain (DBD), p63-DBD does not show cooperative DNA binding properties and is significantly more stable against thermal and chemical denaturation. Here, we determined the solution structure of p63-DBD and show that it is markedly less dynamic than p53-DBD. In addition, we also investigate the effect of a double salt bridge present in p53-DBD, but not in p63-DBD on the cooperative binding behavior and specificity to various DNA sites. Restoration of the salt bridges in p63-DBD by mutagenesis leads to enhanced binding affinity to p53-specific, but not p63-specific response elements. Furthermore, we show that p63-DBD is capable of binding to anti-apoptotic BclxL via its DNA binding interface, a feature that has only been shown for p53 so far. These data suggest that all p53 family members - despite alterations in the specificity and binding affinity - are capable of activating pro-apoptotic pathways in a tissue specific manner.

Identification of a Functional Risk Variant for Pemphigus Vulgaris in the ST18 Gene

Pemphigus vulgaris (PV) is a life-threatening autoimmune mucocutaneous blistering disease caused by disruption of intercellular adhesion due to auto-antibodies directed against epithelial components. Treatment is limited to immunosuppressive agents, which are associated with serious adverse effects. The propensity to develop the disease is in part genetically determined. We therefore reasoned that the delineation of PV genetic basis may point to novel therapeutic strategies. Using a genome-wide association approach, we recently found that genetic variants in the vicinity of the ST18 gene confer a significant risk for the disease. Here, using targeted deep sequencing, we identified a PV-associated variant residing within the ST18 promoter region (p<0.0002; odds ratio = 2.03). This variant was found to drive increased gene transcription in a p53/p63-dependent manner, which may explain the fact that ST18 is up-regulated in the skin of PV patients. We then discovered that when overexpressed, ST18 stimulates PV serum-induced secretion of key inflammatory molecules and contributes to PV serum-induced disruption of keratinocyte cell-cell adhesion, two processes previously implicated in the pathogenesis of PV. Thus, the present findings indicate that ST18 may play a direct role in PV and consequently represents a potential target for the treatment of this disease.

p73 Is Required for Multiciliogenesis and Regulates the Foxj1-Associated Gene Network

We report that p73 is expressed in multiciliated cells (MCCs), is required for MCC differentiation, and directly regulates transcriptional modulators of multiciliogenesis. Loss of ciliary biogenesis provides a unifying mechanism for many phenotypes observed in p73 knockout mice including hydrocephalus; hippocampal dysgenesis; sterility; and chronic inflammation/infection of lung, middle ear, and sinus. Through p73 and p63 ChIP-seq using murine tracheal cells, we identified over 100 putative p73 target genes that regulate MCC differentiation and homeostasis. We validated Foxj1, a transcriptional regulator of multiciliogenesis, and many other cilia-associated genes as direct target genes of p73 and p63. We show p73 and p63 are co-expressed in a subset of basal cells and suggest that p73 marks these cells for MCC differentiation. In summary, p73 is essential for MCC differentiation, functions as a critical regulator of a transcriptome required for MCC differentiation, and, like p63, has an essential role in development of tissues.

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.

SETD1A modulates cell cycle progression through a miRNA network that regulates p53 target genes

Expression of the p53-inducible antiproliferative gene BTG2 is suppressed in many cancers in the absence of inactivating gene mutations, suggesting alternative mechanisms of silencing. Using a shRNA screen targeting 43 histone lysine methyltransferases (KMTs), we show that SETD1A suppresses BTG2 expression through its induction of several BTG2-targeting miRNAs. This indirect but highly specific mechanism, by which a chromatin regulator that mediates transcriptional activating marks can lead to the downregulation of a critical effector gene, is shared with multiple genes in the p53 pathway. Through such miRNA-dependent effects, SETD1A regulates cell cycle progression in vitro and modulates tumorigenesis in mouse xenograft models. Together, these observations help explain the remarkably specific genetic consequences associated with alterations in generic chromatin modulators in cancer.

The p53-inducible antiproliferative gene BTG2 is suppressed in many cancers, in the absence of inactivating gene mutations. Here the authors show that the histone lysine methyltransferase SETD1A suppresses the expression of several p53 target genes including BTG2 by inducing a network of microRNAs.

Gene regulatory mechanisms orchestrated by p63 in epithelial development and related disorders

The transcription factor p63 belongs to the p53 family and is a key regulator in epithelial commitment and development. Mutations in p63 give rise to several epithelial related disorders with defects in skin, limb and orofacial structures. Since the discovery of p63, efforts have been made to identify its target genes using individual gene approaches and to understand p63 function in normal epithelial development and related diseases. Recent genome-wide approaches have identified tens of thousands of potential p63-regulated target genes and regulatory elements, and reshaped the concept of gene regulation orchestrated by p63. These data also provide insights into p63-related disease mechanisms. In this review, we discuss the regulatory role of p63 in normal and diseased epithelial development in light of these novel findings. We also propose future perspectives for dissecting the molecular mechanism of p63-mediated epithelial development and related disorders as well as for potential therapeutic strategies.

p53 and ΔNp63α Coregulate the Transcriptional and Cellular Response to TGFβ and BMP Signals

The TGFβ superfamily regulates a broad range of cellular processes, including proliferation, cell-fate specification, differentiation, and migration. Molecular mechanisms underlying this high degree of pleiotropy and cell-type specificity are not well understood. The TGFβ family is composed of two branches: (i) TGFβs, activins, and nodals, which signal through SMAD2/3, and (ii) bone morphogenetic proteins (BMP), which signal through SMAD1/5/8. SMADs have weak DNA-binding affinity and rely on coactivators and corepressors to specify their transcriptional outputs. This report reveals that p53 and ΔNp63α act as transcriptional partners for SMAD proteins and thereby influence cellular responses to TGFβ and BMPs. Suppression of p53 or overexpression of ΔNp63α synergistically enhance BMP-induced transcription. Mechanistically, p53 and ΔNp63α physically interact with SMAD1/5/8 proteins and co-occupy the promoter region of inhibitor of differentiation (ID2), a prosurvival BMP target gene. Demonstrating further convergence of these pathways, TGFβ-induced canonical BMP regulated transcription in a ΔNp63α- and p53-dependent manner. Furthermore, bioinformatic analyses revealed that SMAD2/3 and ΔNp63α coregulate a significant number of transcripts involved in the regulation of epithelial-to-mesenchymal transition. Thus, p53 and ΔNp63α are transcriptional partners for a subset of TGFβ- and BMP-regulated SMAD target genes in the mammary epithelium. Collectively, these results establish an integrated gene network of SMADs, p53, and ΔNp63α that contribute to EMT and metastasis.

IMPLICATIONS

This study identifies aberrant BMP activation as a result of p53 mutation or ΔNp63α expression.

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

Locally disordered methylation forms the basis of intratumor methylome variation in chronic lymphocytic leukemia

Intratumoral heterogeneity plays a critical role in tumor evolution. To define the contribution of DNA methylation to heterogeneity within tumors, we performed genome-scale bisulfite sequencing of 104 primary chronic lymphocytic leukemias (CLLs). Compared with 26 normal B cell samples, CLLs consistently displayed higher intrasample variability of DNA methylation patterns across the genome, which appears to arise from stochastically disordered methylation in malignant cells. Transcriptome analysis of bulk and single CLL cells revealed that methylation disorder was linked to low-level expression. Disordered methylation was further associated with adverse clinical outcome. We therefore propose that disordered methylation plays a similar role to that of genetic instability, enhancing the ability of cancer cells to search for superior evolutionary trajectories.

Role of chromatin and transcriptional co-regulators in mediating p63-genome interactions in keratinocytes

BACKGROUND

The Transcription Factor (TF) p63 is a master regulator of epidermal development and differentiation as evident from the remarkable skin phenotype of p63 mouse knockouts. Furthermore, ectopic expression of p63 alone is sufficient to convert simple epithelium into stratified epithelial tissues in vivo and p63 is required for efficient transdifferentiation of fibroblasts into keratinocytes. However, little is known about the molecular mechanisms of p63 function, in particular how it selects its target sites in the genome. p63, which acts both as an activator and repressor of transcription, recognizes a canonical binding motif that occurs over 1 million times in the human genome. But, in human keratinocytes less than 12,000 of these sites are bound in vivo suggesting that underlying chromatin architecture and cooperating TFs mediate p63-genome interactions.

RESULTS

We find that the chromatin architecture at p63-bound targets possess distinctive features and can be used to categorize p63 targets into proximal promoters (1%), enhancers (59%) and repressed or inactive (40%) regulatory elements. Our analysis shows that the chromatin modifications H3K4me1, H3K27me3, along with overall chromatin accessibility status can accurately predict bonafide p63-bound sites without a priori DNA sequence information. Interestingly, however there exists a qualitative correlation between the p63 binding motif and accessibility and H3K4me1 levels. Furthermore, we use a comprehensive in silico approach that leverages ENCODE data to identify several known TFs such as AP1, AP2 and novel TFs (RFX5 for e.g.) that can potentially cooperate with p63 to modulate its myriad biological functions in keratinocytes.

CONCLUSIONS

Our analysis shows that p63 bound genomic locations in keratinocytes are accessible, marked by active histone modifications, and co-targeted by other developmentally important transcriptional regulators. Collectively, our results suggest that p63 might actively remodel and/or influence chromatin dynamics at its target sites and in the process dictate its own DNA binding and possibly that of adjacent TFs.

The associations between immunity-related genes and breast cancer prognosis in Korean women

We investigated the role of common genetic variation in immune-related genes on breast cancer disease-free survival (DFS) in Korean women. 107 breast cancer patients of the Seoul Breast Cancer Study (SEBCS) were selected for this study. A total of 2,432 tag single nucleotide polymorphisms (SNPs) in 283 immune-related genes were genotyped with the GoldenGate Oligonucleotide pool assay (OPA). A multivariate Cox-proportional hazard model and polygenic risk score model were used to estimate the effects of SNPs on breast cancer prognosis. Harrell’s C index was calculated to estimate the predictive accuracy of polygenic risk score model. Subsequently, an extended gene set enrichment analysis (GSEA-SNP) was conducted to approximate the biological pathway. In addition, to confirm our results with current evidence, previous studies were systematically reviewed. Sixty-two SNPs were statistically significant at p-value less than 0.05. The most significant SNPs were rs1952438 in SOCS4 gene (hazard ratio (HR) = 11.99, 95% CI = 3.62–39.72, P = 4.84E-05), rs2289278 in TSLP gene (HR = 4.25, 95% CI = 2.10–8.62, P = 5.99E-05) and rs2074724 in HGF gene (HR = 4.63, 95% CI = 2.18–9.87, P = 7.04E-05). In the polygenic risk score model, the HR of women in the 3^rd tertile was 6.78 (95% CI = 1.48–31.06) compared to patients in the 1^st tertile of polygenic risk score. Harrell’s C index was 0.813 with total patients and 0.924 in 4-fold cross validation. In the pathway analysis, 18 pathways were significantly associated with breast cancer prognosis (P<0.1). The IL-6R, IL-8, IL-10RB, IL-12A, and IL-12B was associated with the prognosis of cancer in data of both our study and a previous study. Therefore, our results suggest that genetic polymorphisms in immune-related genes have relevance to breast cancer prognosis among Korean women.

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.

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.

Regulators of global genome repair do not respond to DNA damaging therapy but correlate with survival in melanoma

Nucleotide excision repair (NER) orchestrates the repair of helix distorting DNA damage, induced by both ultraviolet radiation (UVR) and cisplatin. There is evidence that the global genome repair (GGR) arm of NER is dysfunctional in melanoma and it is known to have limited induction in melanoma cell lines after cisplatin treatment. The aims of this study were to examine mRNA transcript levels of regulators of GGR and to investigate the downstream effect on global transcript expression in melanoma cell lines after cisplatin treatment and in melanoma tumours. The GGR regulators, BRCA1 and PCNA, were induced in melanocytes after cisplatin, but not in melanoma cell lines. Transcripts associated with BRCA1, BRCA2, ATM and CHEK2 showed altered expression in melanoma cell lines after cisplatin treatment. In melanoma tumour tissue BRCA1 transcript expression correlated with poor survival and XPB expression correlated with solar elastosis levels. Taken together, these findings provide evidence of the mechanisms underlying NER deficiency in melanoma.

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.

Delineating Molecular Mechanisms of Squamous Tissue Homeostasis and Neoplasia: Focus on p63

Mouse models have informed us that p63 is critical for normal epidermal development and homeostasis. The p53/p63/p73 family is expressed as multiple protein isoforms due to a combination of alternative promoter usage and C-terminal alternative splicing. These isoforms can mimic or interfere with one another, and their balance ultimately determines biological outcome in a context-dependent manner. While not frequently mutated, p63, and in particular the ΔNp63 subclass, is commonly overexpressed in human squamous cell cancers. In vitro keratinocytes and murine transgenic and transplantation models have been invaluable in elucidating the contribution of altered p63 levels to cancer development, and studies have identified the roles for ΔNp63 isoforms in keratinocyte survival and malignant progression, likely due in part to their transcriptional regulatory function. These findings can be extended to human cancers; for example, the novel recognition of NFκB/c-Rel as a downstream effector of p63 has identified a role for NFκB/c-Rel in human squamous cell cancers. These models will be critical in enhancing the understanding of the specific molecular mechanisms of cancer development and progression.

Smad4-Irf6 genetic interaction and TGFβ-mediated IRF6 signaling cascade are crucial for palatal fusion in mice

Cleft palate is one of the most common human birth defects and is associated with multiple genetic and environmental risk factors. Although mutations in the genes encoding transforming growth factor beta (TGFβ) signaling molecules and interferon regulatory factor 6 (Irf6) have been identified as genetic risk factors for cleft palate, little is known about the relationship between TGFβ signaling and IRF6 activity during palate formation. Here, we show that TGFβ signaling regulates expression of Irf6 and the fate of the medial edge epithelium (MEE) during palatal fusion in mice. Haploinsufficiency of Irf6 in mice with basal epithelial-specific deletion of the TGFβ signaling mediator Smad4 (Smad4(fl/fl);K14-Cre;Irf6(+/R84C)) results in compromised p21 expression and MEE persistence, similar to observations in Tgfbr2(fl/fl);K14-Cre mice, although the secondary palate of Irf6(+/R84C) and Smad4(fl/fl);K14-Cre mice form normally. Furthermore, Smad4(fl/fl);K14-Cre;Irf6(+/R84C) mice show extra digits that are consistent with abnormal toe and nail phenotypes in individuals with Van der Woude and popliteal pterygium syndromes, suggesting that the TGFβ/SMAD4/IRF6 signaling cascade might be a well-conserved mechanism in regulating multiple organogenesis. Strikingly, overexpression of Irf6 rescued p21 expression and MEE degeneration in Tgfbr2(fl/fl);K14-Cre mice. Thus, IRF6 and SMAD4 synergistically regulate the fate of the MEE, and TGFβ-mediated Irf6 activity is responsible for MEE degeneration during palatal fusion in mice.