Mutant p63 Affects Epidermal Cell Identity through Rewiring the Enhancer Landscape

TRAPT: a multi-stage fused deep learning framework for predicting transcriptional regulators based on large-scale epigenomic data

It is challenging to identify regulatory transcriptional regulators (TRs), which control gene expression via regulatory elements and epigenomic signals, in context-specific studies on the onset and progression of diseases. The use of large-scale multi-omics epigenomic data enables the representation of the complex epigenomic patterns of control of the regulatory elements and the regulators. Herein, we propose Transcription Regulator Activity Prediction Tool (TRAPT), a multi-modality deep learning framework, which infers regulator activity by learning and integrating the regulatory potentials of target gene cis-regulatory elements and genome-wide binding sites. The results of experiments on 570 TR-related datasets show that TRAPT outperformed state-of-the-art methods in predicting the TRs, especially in terms of forecasting transcription co-factors and chromatin regulators. Moreover, we successfully identify key TRs associated with diseases, genetic variations, cell-fate decisions, and tissues. Our method provides an innovative perspective on identifying TRs by using epigenomic data.

Molecular and Cellular Function of p63 in Skin Development and Genetic Diseases

The transcription factor p63 is a master regulator of multiple ectodermal derivatives. During epidermal commitment, p63 interacts with several chromatin remodeling complexes to transactivate epidermal-specific genes and repress transcription of simple epithelial and nonepithelial genes. In the postnatal epidermis, p63 is required to control the proliferative potential of progenitor cells, maintain epidermal integrity, and contribute to epidermal differentiation. Autosomal dominant sequence variant in p63 cause a spectrum of syndromic disorders that affect several tissues, including or derived from stratified epithelia. In this review, we describe the recent studies that have provided novel insights into disease pathogenesis and potential therapeutic targets.

Nucleosome binding by TP53, TP63, and TP73 is determined by the composition, accessibility, and helical orientation of their binding sites

The TP53 family of transcription factors plays key roles in driving development and combating cancer by regulating gene expression. TP53, TP63, and TP73-the three members of the TP53 family-regulate gene expression by binding to their DNA binding sites, many of which are situated within nucleosomes. To thoroughly examine the nucleosome-binding abilities of the TP53 family, we used Pioneer-seq, a technique that assesses a transcription factor's binding affinity to its DNA-binding sites at all possible positions within the nucleosome core particle. Using Pioneer-seq, we analyzed the binding affinities of TP53, TP63, and TP73 to 10 TP53 family binding sites across the nucleosome core particle. We find that the affinities of TP53, TP63, and TP73 for nucleosomes are primarily determined by the positioning of TP53 family binding sites within nucleosomes; TP53 family members bind strongly to the more accessible edges of nucleosomes but weakly to the less accessible centers of nucleosomes. Our results further show that the DNA-helical orientation of TP53 family binding sites within nucleosomal DNA impacts the nucleosome-binding affinities of TP53 family members, with binding-site composition impacting the affinity of each TP53 family member only when the binding-site location is accessible. Taken together, our results show that the accessibility, composition, and helical orientation of TP53 family binding sites collectively determine the nucleosome-binding affinities of TP53, TP63, and TP73. These findings help explain the rules underlying TP53 family-nucleosome binding and thus provide requisite insight into how we may better control gene expression changes involved in development and tumor suppression.

A multi-step completion process model of cell plasticity

Plasticity is the potential for cells or cell populations to change their phenotypes and behaviors in response to internal or external cues. Plasticity is fundamental to many complex biological processes, yet to date there remains a lack of mathematical models that can elucidate and predict molecular behaviors in a plasticity program. Here, we report a new mathematical framework that models cell plasticity as a multi-step completion process, where the system moves from the initial state along a path guided by multiple intermediate attractors until the final state (i.e. a new homeostasis) is reached. Using omics time-series data as model input, we show that our method fits data well; identifies attractor states by their timing and molecular markers which are well-aligned with domain knowledge; and can make quantitative and time-resolved predictions such as the molecular outcomes of blocking a plasticity program from reaching completion, to an R2 of 0.53-0.63. We demonstrate that application of our model to primary patient-derived data can provide quantitative insights and predictions that may be useful in guiding further research and potential biomedical interventions.

Single-nucleus CUT&RUN elucidates the function of intrinsic and genomics-driven epigenetic heterogeneity in head and neck cancer progression

Interrogating regulatory epigenetic alterations during tumor progression at the resolution of single cells has remained an understudied area of research. Here we developed a highly sensitive single-nucleus CUT&RUN (snCUT&RUN) assay to profile histone modifications in isogenic primary, metastatic, and cisplatin-resistant head and neck squamous cell carcinoma (HNSCC) patient-derived tumor cell lines. We find that the epigenome can be involved in diverse modes to contribute toward HNSCC progression. First, we demonstrate that gene expression changes during HNSCC progression can be comodulated by alterations in both copy number and chromatin activity, driving epigenetic rewiring of cell states. Furthermore, intratumor epigenetic heterogeneity (ITeH) may predispose subclonal populations within the primary tumor to adapt to selective pressures and foster the acquisition of malignant characteristics. In conclusion, snCUT&RUN serves as a valuable addition to the existing toolkit of single-cell epigenomic assays and can be used to dissect the functionality of the epigenome during cancer progression.

Crosstalk between paralogs and isoforms influences p63-dependent regulatory element activity

The p53 family of transcription factors (p53, p63 and p73) regulate diverse organismal processes including tumor suppression, maintenance of genome integrity and the development of skin and limbs. Crosstalk between transcription factors with highly similar DNA binding profiles, like those in the p53 family, can dramatically alter gene regulation. While p53 is primarily associated with transcriptional activation, p63 mediates both activation and repression. The specific mechanisms controlling p63-dependent gene regulatory activity are not well understood. Here, we use massively parallel reporter assays (MPRA) to investigate how local DNA sequence context influences p63-dependent transcriptional activity. Most regulatory elements with a p63 response element motif (p63RE) activate transcription, although binding of the p63 paralog, p53, drives a substantial proportion of that activity. p63RE sequence content and co-enrichment with other known activating and repressing transcription factors, including lineage-specific factors, correlates with differential p63RE-mediated activities. p63 isoforms dramatically alter transcriptional behavior, primarily shifting inactive regulatory elements towards high p63-dependent activity. Our analysis provides novel insight into how local sequence and cellular context influences p63-dependent behaviors and highlights the key, yet still understudied, role of transcription factor paralogs and isoforms in controlling gene regulatory element activity.

Fibroblast-Mediated Macrophage Recruitment Supports Acute Wound Healing

Epithelial and immune cells have long been appreciated for their contribution to the early immune response after injury; however, much less is known about the role of mesenchymal cells. Using single-nuclei RNA sequencing, we defined changes in gene expression associated with inflammation 1 day after wounding in mouse skin. Compared with those in keratinocytes and myeloid cells, we detected enriched expression of proinflammatory genes in fibroblasts associated with deeper layers of the skin. In particular, SCA1+ fibroblasts were enriched for numerous chemokines, including CCL2, CCL7, and IL-33, compared with SCA1- fibroblasts. Genetic deletion of Ccl2 in fibroblasts resulted in fewer wound-bed macrophages and monocytes during injury-induced inflammation, with reduced revascularization and re-epithelialization during the proliferation phase of healing. These findings highlight the important contribution of fibroblast-derived factors to injury-induced inflammation and the impact of immune cell dysregulation on subsequent tissue repair.

The Aryl Hydrocarbon Receptor Regulates Epidermal Differentiation through Transient Activation of TFAP2A

The aryl hydrocarbon receptor (AHR) is an evolutionary conserved environmental sensor identified as an indispensable regulator of epithelial homeostasis and barrier organ function. Molecular signaling cascade and target genes upon AHR activation and their contribution to cell and tissue function are however not fully understood. Multiomics analyses using human skin keratinocytes revealed that upon ligand activation, AHR binds open chromatin to induce expression of transcription factors, for example, TFAP2A, as a swift response to environmental stimuli. The terminal differentiation program, including upregulation of barrier genes, FLG and keratins, was mediated by TFAP2A as a secondary response to AHR activation. The role of AHR-TFAP2A axis in controlling keratinocyte terminal differentiation for proper barrier formation was further confirmed using CRISPR/Cas9 in human epidermal equivalents. Overall, the study provides additional insights into the molecular mechanism behind AHR-mediated barrier function and identifies potential targets for the treatment of skin barrier diseases.

ΔNp63 bookmarks and creates an accessible epigenetic environment for TGFβ-induced cancer cell stemness and invasiveness

BACKGROUND

p63 is a transcription factor with intrinsic pioneer factor activity and pleiotropic functions. Transforming growth factor β (TGFβ) signaling via activation and cooperative action of canonical, SMAD, and non-canonical, MAP-kinase (MAPK) pathways, elicits both anti- and pro-tumorigenic properties, including cell stemness and invasiveness. TGFβ activates the ΔNp63 transcriptional program in cancer cells; however, the link between TGFβ and p63 in unmasking the epigenetic landscape during tumor progression allowing chromatin accessibility and gene transcription, is not yet reported.

METHODS

Small molecule inhibitors, including protein kinase inhibitors and RNA-silencing, provided loss of function analyses. Sphere formation assays in cancer cells, chromatin immunoprecipitation and mRNA expression assays were utilized in order to gain mechanistic evidence. Mass spectrometry analysis coupled to co-immunoprecipitation assays revealed novel p63 interactors and their involvement in p63-dependent transcription.

RESULTS

The sphere-forming capacity of breast cancer cells was enhanced upon TGFβ stimulation and significantly decreased upon ΔNp63 depletion. Activation of TGFβ signaling via p38 MAPK signaling induced ΔNp63 phosphorylation at Ser 66/68 resulting in stabilized ΔNp63 protein with enhanced DNA binding properties. TGFβ stimulation altered the ratio of H3K27ac and H3K27me3 histone modification marks, pointing towards higher H3K27ac and increased p300 acetyltransferase recruitment to chromatin. By silencing the expression of ΔNp63, the TGFβ effect on chromatin remodeling was abrogated. Inhibition of H3K27me3, revealed the important role of TGFβ as the upstream signal for guiding ΔNp63 to the TGFβ/SMAD gene loci, as well as the indispensable role of ΔNp63 in recruiting histone modifying enzymes, such as p300, to these genomic regions, regulating chromatin accessibility and gene transcription. Mechanistically, TGFβ through SMAD activation induced dissociation of ΔNp63 from NURD or NCOR/SMRT histone deacetylation complexes, while promoted the assembly of ΔNp63-p300 complexes, affecting the levels of histone acetylation and the outcome of ΔNp63-dependent transcription.

CONCLUSIONS

ΔNp63, phosphorylated and recruited by TGFβ to the TGFβ/SMAD/ΔNp63 gene loci, promotes chromatin accessibility and transcription of target genes related to stemness and cell invasion.

Deep skin fibroblast-mediated macrophage recruitment supports acute wound healing

Epithelial and immune cells have long been appreciated for their contribution to the early immune response after injury; however, much less is known about the role of mesenchymal cells. Using single nuclei RNA-sequencing, we defined changes in gene expression associated with inflammation at 1-day post-wounding (dpw) in mouse skin. Compared to keratinocytes and myeloid cells, we detected enriched expression of pro-inflammatory genes in fibroblasts associated with deeper layers of the skin. In particular, SCA1+ fibroblasts were enriched for numerous chemokines, including CCL2, CCL7, and IL33 compared to SCA1- fibroblasts. Genetic deletion of Ccl2 in fibroblasts resulted in fewer wound bed macrophages and monocytes during injury-induced inflammation with reduced revascularization and re-epithelialization during the proliferation phase of healing. These findings highlight the important contribution of deep skin fibroblast-derived factors to injury-induced inflammation and the impact of immune cell dysregulation on subsequent tissue repair.

Identification of the regulatory circuit governing corneal epithelial fate determination and disease

The transparent corneal epithelium in the eye is maintained through the homeostasis regulated by limbal stem cells (LSCs), while the nontransparent epidermis relies on epidermal keratinocytes for renewal. Despite their cellular similarities, the precise cell fates of these two types of epithelial stem cells, which give rise to functionally distinct epithelia, remain unknown. We performed a multi-omics analysis of human LSCs from the cornea and keratinocytes from the epidermis and characterized their molecular signatures, highlighting their similarities and differences. Through gene regulatory network analyses, we identified shared and cell type-specific transcription factors (TFs) that define specific cell fates and established their regulatory hierarchy. Single-cell RNA-seq (scRNA-seq) analyses of the cornea and the epidermis confirmed these shared and cell type-specific TFs. Notably, the shared and LSC-specific TFs can cooperatively target genes associated with corneal opacity. Importantly, we discovered that FOSL2, a direct PAX6 target gene, is a novel candidate associated with corneal opacity, and it regulates genes implicated in corneal diseases. By characterizing molecular signatures, our study unveils the regulatory circuitry governing the LSC fate and its association with corneal opacity.

Integrated single-cell chromatin and transcriptomic analyses of human scalp identify gene-regulatory programs and critical cell types for hair and skin diseases

Genome-wide association studies have identified many loci associated with hair and skin disease, but identification of causal variants requires deciphering of gene-regulatory networks in relevant cell types. We generated matched single-cell chromatin profiles and transcriptomes from scalp tissue from healthy controls and patients with alopecia areata, identifying diverse cell types of the hair follicle niche. By interrogating these datasets at multiple levels of cellular resolution, we infer 50-100% more enhancer-gene links than previous approaches and show that aggregate enhancer accessibility for highly regulated genes predicts expression. We use these gene-regulatory maps to prioritize cell types, genes and causal variants implicated in the pathobiology of androgenetic alopecia (AGA), eczema and other complex traits. AGA genome-wide association studies signals are enriched in dermal papilla regulatory regions, supporting the role of these cells as drivers of AGA pathogenesis. Finally, we train machine learning models to nominate single-nucleotide polymorphisms that affect gene expression through disruption of transcription factor binding, predicting candidate functional single-nucleotide polymorphism for AGA and eczema.

The aryl hydrocarbon receptor regulates epidermal differentiation through transient activation of TFAP2A

The aryl hydrocarbon receptor (AHR) is an evolutionary conserved environmental sensor identified as indispensable regulator of epithelial homeostasis and barrier organ function. Molecular signaling cascade and target genes upon AHR activation and their contribution to cell and tissue function are however not fully understood. Multi-omics analyses using human skin keratinocytes revealed that, upon ligand activation, AHR binds open chromatin to induce expression of transcription factors (TFs), e.g., Transcription Factor AP-2α (TFAP2A), as a swift response to environmental stimuli. The terminal differentiation program including upregulation of barrier genes, filaggrin and keratins, was mediated by TFAP2A as a secondary response to AHR activation. The role of AHR-TFAP2A axis in controlling keratinocyte terminal differentiation for proper barrier formation was further confirmed using CRISPR/Cas9 in human epidermal equivalents. Overall, the study provides novel insights into the molecular mechanism behind AHR-mediated barrier function and potential novel targets for the treatment of skin barrier diseases.

Convergence of YAP/TAZ, TEAD and TP63 activity is associated with bronchial premalignant severity and progression

BACKGROUND

Bronchial premalignant lesions (PMLs) are composed primarily of cells resembling basal epithelial cells of the airways, which through poorly understood mechanisms have the potential to progress to lung squamous cell carcinoma (LUSC). Despite ongoing efforts that have mapped gene expression and cell diversity across bronchial PML pathologies, signaling and transcriptional events driving malignancy are poorly understood. Evidence has suggested key roles for the Hippo pathway effectors YAP and TAZ and associated TEAD and TP63 transcription factor families in bronchial basal cell biology and LUSC. In this study we examine the functional association of YAP/TAZ, TEADs and TP63 in bronchial epithelial cells and PMLs.

METHODS

We performed RNA-seq in primary human bronchial epithelial cells following small interfering RNA (siRNA)-mediated depletion of YAP/TAZ, TEADs or TP63, and combined these data with ChIP-seq analysis of these factors. Directly activated or repressed genes were identified and overlapping genes were profiled across gene expression data obtained from progressive or regressive human PMLs and across lung single cell RNA-seq data sets.

RESULTS

Analysis of genes regulated by YAP/TAZ, TEADs, and TP63 in human bronchial epithelial cells revealed a converged transcriptional network that is strongly associated with the pathological progression of bronchial PMLs. Our observations suggest that YAP/TAZ-TEAD-TP63 associate to cooperatively promote basal epithelial cell proliferation and repress signals associated with interferon responses and immune cell communication. Directly repressed targets we identified include the MHC Class II transactivator CIITA, which is repressed in progressive PMLs and associates with adaptive immune responses in the lung. Our findings provide molecular insight into the control of gene expression events driving PML progression, including those contributing to immune evasion, offering potential new avenues for lung cancer interception.

CONCLUSIONS

Our study identifies important gene regulatory functions for YAP/TAZ-TEAD-TP63 in the early stages of lung cancer development, which notably includes immune-suppressive roles, and suggest that an assessment of the activity of this transcriptional complex may offer a means to identify immune evasive bronchial PMLs and serve as a potential therapeutic target.

Disease-related p63 DBD mutations impair DNA binding by distinct mechanisms and varying degree

The transcription factor p63 shares a high sequence identity with the tumour suppressor p53 which manifests itself in high structural similarity and preference for DNA sequences. Mutations in the DNA binding domain (DBD) of p53 have been studied in great detail, enabling a general mechanism-based classification. In this study we provide a detailed investigation of all currently known mutations in the p63 DBD, which are associated with developmental syndromes, by measuring their impact on transcriptional activity, DNA binding affinity, zinc binding capacity and thermodynamic stability. Some of the mutations we have further characterized with respect to their ability to convert human dermal fibroblasts into induced keratinocytes. Here we propose a classification of the p63 DBD mutations based on the four different mechanisms of DNA binding impairment which we identified: direct DNA contact, zinc finger region, H2 region, and dimer interface mutations. The data also demonstrate that, in contrast to p53 cancer mutations, no p63 mutation induces global unfolding and subsequent aggregation of the domain. The dimer interface mutations that affect the DNA binding affinity by disturbing the interaction between the individual DBDs retain partial DNA binding capacity which correlates with a milder patient phenotype.

ΔNp63 drives dysplastic alveolar remodeling and restricts epithelial plasticity upon severe lung injury

The lung exhibits a robust, multifaceted regenerative response to severe injuries such as influenza infection, during which quiescent lung-resident epithelial progenitors participate in two distinct reparative pathways: functionally beneficial regeneration via alveolar type 2 (AT2) cell proliferation and differentiation, and dysplastic tissue remodeling via intrapulmonary airway-resident basal p63+ progenitors. Here we show that the basal cell transcription factor ΔNp63 is required for intrapulmonary basal progenitors to participate in dysplastic alveolar remodeling following injury. We find that ΔNp63 restricts the plasticity of intrapulmonary basal progenitors by maintaining either active or repressive histone modifications at key differentiation gene loci. Following loss of ΔNp63, intrapulmonary basal progenitors are capable of either airway or alveolar differentiation depending on their surrounding environment both in vitro and in vivo. Uncovering these regulatory mechanisms of dysplastic repair and lung basal cell fate choice highlight potential therapeutic targets to promote functional alveolar regeneration following severe lung injuries.

Molecular Pathway-Based Classification of Ectodermal Dysplasias: First Five-Yearly Update

To keep pace with the rapid advancements in molecular genetics and rare diseases research, we have updated the list of ectodermal dysplasias based on the latest classification approach that was adopted in 2017 by an international panel of experts. For this purpose, we searched the databases PubMed and OMIM for the term "ectodermal dysplasia", referring mainly to changes in the last 5 years. We also tried to obtain information about those diseases on which the last scientific report appeared more than 15 years ago by contacting the authors of the most recent publication. A group of experts, composed of researchers who attended the 8th International Conference on Ectodermal Dysplasias and additional members of the previous classification panel, reviewed the proposed amendments and agreed on a final table listing all 49 currently known ectodermal dysplasias for which the molecular genetic basis has been clarified, including 15 new entities. A newly reported ectodermal dysplasia, linked to the gene LRP6, is described here in more detail. These ectodermal dysplasias, in the strict sense, should be distinguished from syndromes with features of ectodermal dysplasia that are related to genes extraneous to the currently known pathways involved in ectodermal development. The latter group consists of 34 syndromes which had been placed on the previous list of ectodermal dysplasias, but most if not all of them could actually be classified elsewhere. This update should streamline the classification of ectodermal dysplasias, provide guidance to the correct diagnosis of rare disease entities, and facilitate the identification of individuals who could benefit from novel treatment options.

Identification of transcription factors dictating blood cell development using a bidirectional transcription network-based computational framework

Advanced computational methods exploit gene expression and epigenetic datasets to predict gene regulatory networks controlled by transcription factors (TFs). These methods have identified cell fate determining TFs but require large amounts of reference data and experimental expertise. Here, we present an easy to use network-based computational framework that exploits enhancers defined by bidirectional transcription, using as sole input CAGE sequencing data to correctly predict TFs key to various human cell types. Next, we applied this Analysis Algorithm for Networks Specified by Enhancers based on CAGE (ANANSE-CAGE) to predict TFs driving red and white blood cell development, and THP-1 leukemia cell immortalization. Further, we predicted TFs that are differentially important to either cell line- or primary- associated MLL-AF9-driven gene programs, and in primary MLL-AF9 acute leukemia. Our approach identified experimentally validated as well as thus far unexplored TFs in these processes. ANANSE-CAGE will be useful to identify transcription factors that are key to any cell fate change using only CAGE-seq data as input.

Structural diversity of p63 and p73 isoforms

Abstract

The p53 protein family is the most studied protein family of all. Sequence analysis and structure determination have revealed a high similarity of crucial domains between p53, p63 and p73. Functional studies, however, have shown a wide variety of different tasks in tumor suppression, quality control and development. Here we review the structure and organization of the individual domains of p63 and p73, the interaction of these domains in the context of full-length proteins and discuss the evolutionary origin of this protein family.

Facts

Open questions

Enhanced pro-apoptosis gene signature following the activation of TAp63α in oocytes upon γ irradiation

Specialized surveillance mechanisms are essential to maintain the genetic integrity of germ cells, which are not only the source of all somatic cells but also of the germ cells of the next generation. DNA damage and chromosomal aberrations are, therefore, not only detrimental for the individual but affect the entire species. In oocytes, the surveillance of the structural integrity of the DNA is maintained by the p53 family member TAp63α. The TAp63α protein is highly expressed in a closed and inactive state and gets activated to the open conformation upon the detection of DNA damage, in particular DNA double-strand breaks. To understand the cellular response to DNA damage that leads to the TAp63α triggered oocyte death we have investigated the RNA transcriptome of oocytes following irradiation at different time points. The analysis shows enhanced expression of pro-apoptotic and typical p53 target genes such as CDKn1a or Mdm2, concomitant with the activation of TAp63α. While DNA repair genes are not upregulated, inflammation-related genes become transcribed when apoptosis is initiated by activation of STAT transcription factors. Furthermore, comparison with the transcriptional profile of the ΔNp63α isoform from other studies shows only a minimal overlap, suggesting distinct regulatory programs of different p63 isoforms.

Transcriptional and epigenetic regulation in thymic epithelial cells

The thymus is required for the development of both adaptive and innate-like T cell subsets. There is keen interest in manipulating thymic function for therapeutic purposes in circumstances of autoimmunity, immunodeficiency, and for purposes of immunotherapy. Within the thymus, thymic epithelial cells play essential roles in directing T cell development. Several transcription factors are known to be essential for thymic epithelial cell development and function, and a few transcription factors have been studied in considerable detail. However, the role of many other transcription factors is less well understood. Further, it is likely that roles exist for other transcription factors not yet known to be important in thymic epithelial cells. Recent progress in understanding of thymic epithelial cell heterogeneity has provided some new insight into transcriptional requirements in subtypes of thymic epithelial cells. However, it is unknown whether progenitors of thymic epithelial cells exist in the adult thymus, and consequently, developmental relationships linking putative precursors with differentiated cell types are poorly understood. While we do not presently possess a clear understanding of stage-specific requirements for transcription factors in thymic epithelial cells, new single-cell transcriptomic and epigenomic technologies should enable rapid progress in this field. Here, we review our current knowledge of transcription factors involved in the development, maintenance, and function of thymic epithelial cells, and the mechanisms by which they act.

ANANSE: an enhancer network-based computational approach for predicting key transcription factors in cell fate determination

Proper cell fate determination is largely orchestrated by complex gene regulatory networks centered around transcription factors. However, experimental elucidation of key transcription factors that drive cellular identity is currently often intractable. Here, we present ANANSE (ANalysis Algorithm for Networks Specified by Enhancers), a network-based method that exploits enhancer-encoded regulatory information to identify the key transcription factors in cell fate determination. As cell type-specific transcription factors predominantly bind to enhancers, we use regulatory networks based on enhancer properties to prioritize transcription factors. First, we predict genome-wide binding profiles of transcription factors in various cell types using enhancer activity and transcription factor binding motifs. Subsequently, applying these inferred binding profiles, we construct cell type-specific gene regulatory networks, and then predict key transcription factors controlling cell fate transitions using differential networks between cell types. This method outperforms existing approaches in correctly predicting major transcription factors previously identified to be sufficient for trans-differentiation. Finally, we apply ANANSE to define an atlas of key transcription factors in 18 normal human tissues. In conclusion, we present a ready-to-implement computational tool for efficient prediction of transcription factors in cell fate determination and to study transcription factor-mediated regulatory mechanisms. ANANSE is freely available at https://github.com/vanheeringen-lab/ANANSE.

Terminal keratinocyte differentiation in vitro is associated with a stable DNA methylome

The epidermal compartment of the skin is regenerated constantly by proliferation of epidermal keratinocytes. Differentiation of a subset of these keratinocytes allows the epidermis to retain its barrier properties. Regulation of keratinocyte fate-whether to remain proliferative or terminally differentiate-is complex and not fully understood. The objective of our study was to assess if DNA methylation changes contribute to the regulation of keratinocyte fate. We employed genome-wide MethylationEPIC beadchip array measuring approximately 850 000 probes combined with RNA sequencing of in vitro cultured non-differentiated and terminally differentiated adult human primary keratinocytes. We did not observe a correlation between methylation status and transcriptome changes. Moreover, only two differentially methylated probes were detected, of which one was located in the TRIM29 gene. Although TRIM29 knock-down resulted in lower expression levels of terminal differentiation genes, these changes were minor. From these results, we conclude that-in our in vitro experimental setup-it is unlikely that changes in DNA methylation have an important regulatory role in terminal keratinocyte differentiation.

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.

ΔNp63 is a pioneer factor that binds inaccessible chromatin and elicits chromatin remodeling

BACKGROUND

ΔNp63 is a master transcriptional regulator playing critical roles in epidermal development and other cellular processes. Recent studies suggest that ΔNp63 functions as a pioneer factor that can target its binding sites within inaccessible chromatin and induce chromatin remodeling.

METHODS

In order to examine if ΔNp63 can bind to inaccessible chromatin and to determine if specific histone modifications are required for binding, we induced ΔNp63 expression in two p63-naïve cell lines. ΔNp63 binding was then examined by ChIP-seq and the chromatin at ΔNp63 targets sites was examined before and after binding. Further analysis with competitive nucleosome binding assays was used to determine how ΔNp63 directly interacts with nucleosomes.

RESULTS

Our results show that before ΔNp63 binding, targeted sites lack histone modifications, indicating ΔNp63's capability to bind at unmodified chromatin. Moreover, the majority of the sites that are bound by ectopic ΔNp63 expression exist in an inaccessible state. Once bound, ΔNp63 induces acetylation of the histone and the repositioning of nucleosomes at its binding sites. Further analysis with competitive nucleosome binding assays reveal that ΔNp63 can bind directly to nucleosome edges with significant binding inhibition occurring within 50 bp of the nucleosome dyad.

CONCLUSION

Overall, our results demonstrate that ΔNp63 is a pioneer factor that binds nucleosome edges at inaccessible and unmodified chromatin sites and induces histone acetylation and nucleosome repositioning.

The foggy world(s) of p63 isoform regulation in normal cells and cancer

The p53 family member p63 exists as two major protein variants (TAp63 and ΔNp63) with distinct expression patterns and functional properties. Whilst downstream target genes of p63 have been studied intensively, how p63 variants are themselves controlled has been relatively neglected. Here, we review advances in understanding ΔNp63 and TAp63 regulation, highlighting their distinct pathways. TAp63 has roles in senescence and metabolism, and in germ cell genome maintenance, where it is activated post-transcriptionally by phosphorylation cascades after DNA damage. The function and regulation of TAp63 in mesenchymal and haematopoietic cells is less clear but may involve epigenetic control through DNA methylation. ΔNp63 functions to maintain stem/progenitor cells in various epithelia and is overexpressed in squamous and certain other cancers. ΔNp63 is transcriptionally regulated through multiple enhancers in concert with chromatin modifying proteins. Many signalling pathways including growth factors, morphogens, inflammation, and the extracellular matrix influence ΔNp63 levels, with inconsistent results reported. There is also evidence for reciprocal regulation, including ΔNp63 activating its own transcription. ΔNp63 is downregulated during cell differentiation through transcriptional regulation, while post-transcriptional events cause proteasomal degradation. Throughout the review, we identify knowledge gaps and highlight discordances, providing potential explanations including cell-context and cell-matrix interactions. Identifying individual p63 variants has roles in differential diagnosis and prognosis, and understanding their regulation suggests clinically approved agents for targeting p63 that may be useful combination therapies for selected cancer patients. © 2021 The Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

An NKX2-1GFP and TP63tdTomato dual fluorescent reporter for the investigation of human lung basal cell biology

Basal cells are multipotent stem cells responsible for the repair and regeneration of all the epithelial cell types present in the proximal lung. In mice, the elusive origins of basal cells and their contribution to lung development were recently revealed by high-resolution, lineage tracing studies. It however remains unclear if human basal cells originate and participate in lung development in a similar fashion, particularly with mounting evidence for significant species-specific differences in this process. To address this outstanding question, in the last several years differentiation protocols incorporating human pluripotent stem cells (hPSC) have been developed to produce human basal cells in vitro with varying efficiencies. To facilitate this endeavour, we introduced tdTomato into the human TP63 gene, whose expression specifically labels basal cells, in the background of a previously described hPSC line harbouring an NKX2-1GFP reporter allele. The functionality and specificity of the NKX2-1GFP;TP63tdTomato hPSC line was validated by directed differentiation into lung progenitors as well as more specialised lung epithelial subtypes using an organoid platform. This dual fluorescent reporter hPSC line will be useful for tracking, isolating and expanding basal cells from heterogenous differentiation cultures for further study.

SPT6 promotes epidermal differentiation and blockade of an intestinal-like phenotype through control of transcriptional elongation

In adult tissue, stem and progenitor cells must tightly regulate the balance between proliferation and differentiation to sustain homeostasis. How this exquisite balance is achieved is an area of active investigation. Here, we show that epidermal genes, including ~30% of induced differentiation genes already contain stalled Pol II at the promoters in epidermal stem and progenitor cells which is then released into productive transcription elongation upon differentiation. Central to this process are SPT6 and PAF1 which are necessary for the elongation of these differentiation genes. Upon SPT6 or PAF1 depletion there is a loss of human skin differentiation and stratification. Unexpectedly, loss of SPT6 also causes the spontaneous transdifferentiation of epidermal cells into an intestinal-like phenotype due to the stalled transcription of the master regulator of epidermal fate P63. Our findings suggest that control of transcription elongation through SPT6 plays a prominent role in adult somatic tissue differentiation and the inhibition of alternative cell fate choices.

Isoform-Specific Roles of Mutant p63 in Human Diseases

Simple Summary

The protein p63 belongs to the family of the p53 tumor suppressor. Mouse models have, however, shown that it is not a classical tumor suppressor but instead involved in developmental processes. Mutations in the p63 gene cause several developmental defects in human patients characterized by limb deformation, cleft lip/palate, and ectodermal dysplasia due to p63’s role as a master regulator of epidermal development. In addition, p63 plays a key role as a quality control factor in oocytes and p63 mutations can result either in compromised genetic quality control or premature cell death of all oocytes.

Abstract

The p63 gene encodes a master regulator of epidermal commitment, development, and differentiation. Heterozygous mutations in the DNA binding domain cause Ectrodactyly, Ectodermal Dysplasia, characterized by limb deformation, cleft lip/palate, and ectodermal dysplasia while mutations in in the C-terminal domain of the α-isoform cause Ankyloblepharon-Ectodermal defects-Cleft lip/palate (AEC) syndrome, a life-threatening disorder characterized by skin fragility, severe, long-lasting skin erosions, and cleft lip/palate. The molecular disease mechanisms of these syndromes have recently become elucidated and have enhanced our understanding of the role of p63 in epidermal development. Here we review the molecular cause and functional consequences of these p63-mutations for skin development and discuss the consequences of p63 mutations for female fertility.

Distinct p63 and p73 Protein Interactions Predict Specific Functions in mRNA Splicing and Polyploidy Control in Epithelia

Epithelial organs are the first barrier against microorganisms and genotoxic stress, in which the p53 family members p63 and p73 have both overlapping and distinct functions. Intriguingly, p73 displays a very specific localization to basal epithelial cells in human tissues, while p63 is expressed in both basal and differentiated cells. Here, we analyse systematically the literature describing p63 and p73 protein-protein interactions to reveal distinct functions underlying the aforementioned distribution. We have found that p73 and p63 cooperate in the genome stability surveillance in proliferating cells; p73 specific interactors contribute to the transcriptional repression, anaphase promoting complex and spindle assembly checkpoint, whereas p63 specific interactors play roles in the regulation of mRNA processing and splicing in both proliferating and differentiated cells. Our analysis reveals the diversification of the RNA and DNA specific functions within the p53 family.

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.

TP63 links chromatin remodeling and enhancer reprogramming to epidermal differentiation and squamous cell carcinoma development

Squamous cell carcinoma (SCC) is an aggressive malignancy that can originate from various organs. TP63 is a master regulator that plays an essential role in epidermal differentiation. It is also a lineage-dependent oncogene in SCC. ΔNp63α is the prominent isoform of TP63 expressed in epidermal cells and SCC, and overexpression promotes SCC development through a variety of mechanisms. Recently, ΔNp63α was highlighted to act as an epidermal-specific pioneer factor that binds closed chromatin and enhances chromatin accessibility at epidermal enhancers. ΔNp63α coordinates chromatin-remodeling enzymes to orchestrate the tissue-specific enhancer landscape and three-dimensional high-order architecture of chromatin. Moreover, ΔNp63α establishes squamous-like enhancer landscapes to drive oncogenic target expression during SCC development. Importantly, ΔNp63α acts as an upstream regulator of super enhancers to activate a number of oncogenic transcripts linked to poor prognosis in SCC. Mechanistically, ΔNp63α activates genes transcription through physically interacting with a number of epigenetic modulators to establish enhancers and enhance chromatin accessibility. In contrast, ΔNp63α also represses gene transcription via interacting with repressive epigenetic regulators. ΔNp63α expression is regulated at multiple levels, including transcriptional, post-transcriptional, and post-translational levels. In this review, we summarize recent advances of p63 in epigenomic and transcriptional control, as well as the mechanistic regulation of p63.

Attenuation of p53 mutant as an approach for treatment Her2-positive cancer

Breast cancer is one of the world's leading causes of oncological disease-related death. It is characterized by a high degree of heterogeneity on the clinical, morphological, and molecular levels. Based on molecular profiling breast carcinomas are divided into several subtypes depending on the expression of a number of cell surface receptors, e.g., ER, PR, and HER2. The Her2-positive subtype occurs in ~10-15% of all cases of breast cancer, and is characterized by a worse prognosis of patient survival. This is due to a high and early relapse rate, as well as an increased level of metastases. Several FDA-approved drugs for the treatment of Her2-positive tumors have been developed, although eventually cancer cells develop drug resistance. These drugs target either the homo- or heterodimerization of Her2 receptors or the receptors' RTK activity, both of them being critical for the proliferation of cancer cells. Notably, Her2-positive cancers also frequently harbor mutations in the TP53 tumor suppressor gene, which exacerbates the unfavorable prognosis. In this review, we describe the molecular mechanisms of RTK-specific drugs and discuss new perspectives of combinatorial treatment of Her2-positive cancers through inhibition of the mutant form of p53.

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.

p63-related signaling at a glance

p63 (also known as TP63) is a transcription factor of the p53 family, along with p73. Multiple isoforms of p63 have been discovered and these have diverse functions encompassing a wide array of cell biology. p63 isoforms are implicated in lineage specification, proliferative potential, differentiation, cell death and survival, DNA damage response and metabolism. Furthermore, p63 is linked to human disease states including cancer. p63 is critical to many aspects of cell signaling, and in this Cell science at a glance article and the accompanying poster, we focus on the signaling cascades regulating TAp63 and ΔNp63 isoforms and those that are regulated by TAp63 and ΔNp63, as well the role of p63 in disease.

RUNX transcription factors: orchestrators of development

RUNX transcription factors orchestrate many different aspects of biology, including basic cellular and developmental processes, stem cell biology and tumorigenesis. In this Primer, we introduce the molecular hallmarks of the three mammalian RUNX genes, RUNX1, RUNX2 and RUNX3, and discuss the regulation of their activities and their mechanisms of action. We then review their crucial roles in the specification and maintenance of a wide array of tissues during embryonic development and adult homeostasis.

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.

p63 cooperates with CTCF to modulate chromatin architecture in skin keratinocytes

The transcription factor p63 regulates epidermal genes and the enhancer landscape in skin keratinocytes. Its molecular function in controlling the chromatin structure is, however, not yet completely understood. Here, we integrated multi-omics profiles, including the transcriptome, transcription factor DNA-binding and chromatin accessibility, in skin keratinocytes isolated from EEC syndrome patients carrying p63 mutations, to examine the role of p63 in shaping the chromatin architecture. We found decreased chromatin accessibility in p63- and CTCF-bound open chromatin regions that potentially contributed to gene deregulation in mutant keratinocytes. Cooperation of p63 and CTCF seemed to assist chromatin interactions between p63-bound enhancers and gene promoters in skin keratinocytes. Our study suggests an intriguing model where cell type-specific transcription factors such as p63 cooperate with the genome organizer CTCF in the three-dimensional chromatin space to regulate the transcription program important for the proper cell identity.

Control of p53-dependent transcription and enhancer activity by the p53 family member p63

Transcriptional activation by p53 provides powerful, organism-wide tumor suppression. We hypothesized that the local chromatin environment, including differential enhancer activities, contributes to various p53-dependent transcriptional activities in different cell types during stress-induced signaling. In this work, using ChIP-sequencing, immunoblotting, quantitative PCR, and computational analyses across various mammalian cell lines, we demonstrate that the p53-induced transcriptome varies by cell type, reflects cell type-specific activities, and is considerably broader than previously anticipated. We found that these molecular events are strongly influenced by p53's engagement with differentially active cell type-specific enhancers and promoters. We also observed that p53 activity depends on the p53 family member tumor protein p63 in epithelial cell types. Notably, we demonstrate that p63 is required for epithelial enhancer identity, including enhancers used by p53 during stress-dependent signaling. Loss of p63, but not p53, caused site-specific depletion of enhancer-associated chromatin modifications, suggesting that p63 functions as an enhancer maintenance factor in epithelial cells. Additionally, a subset of epithelial-specific enhancers depends on the activity of p63 providing a direct link between lineage determination and enhancer structure. These results suggest that a broad, cell-intrinsic mechanism controls p53-dependent cellular stress response through differential regulation of cis-regulatory elements.

p63 establishes epithelial enhancers at critical craniofacial development genes

The transcription factor p63 is a key mediator of epidermal development. Point mutations in p63 in patients lead to developmental defects, including orofacial clefting. To date, knowledge on how pivotal the role of p63 is in human craniofacial development is limited. Using an inducible transdifferentiation model, combined with epigenomic sequencing and multicohort meta-analysis of genome-wide association studies data, we show that p63 establishes enhancers at craniofacial development genes to modulate their transcription. Disease-specific substitution mutation in the DNA binding domain or sterile alpha motif protein interaction domain of p63, respectively, eliminates or reduces establishment of these enhancers. We show that enhancers established by p63 are highly enriched for single-nucleotide polymorphisms associated with nonsyndromic cleft lip ± cleft palate (CL/P). These orthogonal approaches indicate a strong molecular link between p63 enhancer function and CL/P, illuminating molecular mechanisms underlying this developmental defect and revealing vital regulatory elements and new candidate causative genes.

C-terminal α Domain of p63 Binds to p300 to Coactivate β-Catenin

TP63 (p63), a member of the tumor suppressor TP53 (p53) gene family, is essential for ectodermal tissue development and suppresses malignant progression of carcinomas. The most abundant isoform, ΔNp63α (referred to as p63), lacks the N-terminal transactivation (TA) domain, and was originally characterized as a dominant-negative type suppressor against p53 family proteins. It also binds to TCF/LEF to inhibit β-catenin. Nevertheless, transcriptional activation by p63 has also been observed in varied systems. To understand the puzzling results, we analyzed the structure–function relationship of p63 in the control of β-catenin-dependent transcription. p63 acted as a suppressor of moderately induced β-catenin. However, when nuclear targeted S33Y β-catenin was applied to cause the maximum enhancer activation, p63 displayed a β-catenin-coactivating function. The DNA-binding domain of p63 and the target sequence facilitated it. Importantly, we newly found that, despite the absence of TA domain, p63 was associated with p300, a general adaptor protein and chromatin modifier causing transcriptional activation. C-terminal α domain of p63 was essential for p300-binding and for the coactivator function. These results were related to endogenous p63-p300 complex formation and Wnt/β-catenin-responsive gene regulation by p63 in squamous cell carcinoma lines. The novel p63-p300 interaction may be involved in positive regulation of gene expression in tissue development and carcinogenesis.