p63 establishes epithelial enhancers at critical craniofacial development genes

Functional analysis of ESRP1/2 gene variants and CTNND1 isoforms in orofacial cleft pathogenesis

Orofacial cleft (OFC) is a common human congenital anomaly. Epithelial-specific RNA splicing regulators ESRP1 and ESRP2 regulate craniofacial morphogenesis and their disruption result in OFC in zebrafish, mouse and humans. Using esrp1/2 mutant zebrafish and murine Py2T cell line models, we functionally tested the pathogenicity of human ESRP1/2 gene variants. We found that many variants predicted by in silico methods to be pathogenic were functionally benign. Esrp1 also regulates the alternative splicing of Ctnnd1 and these genes are co-expressed in the embryonic and oral epithelium. In fact, over-expression of ctnnd1 is sufficient to rescue morphogenesis of epithelial-derived structures in esrp1/2 zebrafish mutants. Additionally, we identified 13 CTNND1 variants from genome sequencing of OFC cohorts, confirming CTNND1 as a key gene in human OFC. This work highlights the importance of functional assessment of human gene variants and demonstrates the critical requirement of Esrp-Ctnnd1 acting in the embryonic epithelium to regulate palatogenesis.

Functional analysis of ESRP1/2 gene variants and CTNND1 isoforms in orofacial cleft pathogenesis

Orofacial cleft (OFC) is a common human congenital anomaly. Epithelial-specific RNA splicing regulators ESRP1 and ESRP2 regulate craniofacial morphogenesis and their disruption result in OFC in zebrafish, mouse and humans. Using esrp1/2 mutant zebrafish and murine Py2T cell line models, we functionally tested the pathogenicity of human ESRP1/2 gene variants. We found that many variants predicted by in silico methods to be pathogenic were functionally benign. Esrp1 also regulates the alternative splicing of Ctnnd1 and these genes are co-expressed in the embryonic and oral epithelium. In fact, over-expression of ctnnd1 is sufficient to rescue morphogenesis of epithelial-derived structures in esrp1/2 zebrafish mutants. Additionally, we identified 13 CTNND1 variants from genome sequencing of OFC cohorts, confirming CTNND1 as a key gene in human OFC. This work highlights the importance of functional assessment of human gene variants and demonstrates the critical requirement of Esrp-Ctnnd1 acting in the embryonic epithelium to regulate palatogenesis.

Context dependent activity of p63-bound gene regulatory elements

The p53 family of transcription factors regulate numerous organismal processes including the development of skin and limbs, ciliogenesis, and preservation of genetic integrity and tumor suppression. p53 family members control these processes and gene expression networks through engagement with DNA sequences within gene regulatory elements. Whereas p53 binding to its cognate recognition sequence is strongly associated with transcriptional activation, p63 can mediate both activation and repression. How the DNA sequence of p63-bound gene regulatory elements is linked to these varied activities is not yet understood. Here, we use massively parallel reporter assays (MPRA) in a range of cellular and genetic contexts to investigate the influence of DNA sequence on p63-mediated transcription. Most regulatory elements with a p63 response element motif (p63RE) activate transcription, with those sites bound by p63 more frequently or adhering closer to canonical p53 family response element sequences driving higher transcriptional output. The most active regulatory elements are those also capable of binding p53. Elements uniquely bound by p63 have varied activity, with p63RE-mediated repression associated with lower overall GC content in flanking sequences. Comparison of activity across cell lines suggests differential activity of elements may be regulated by a combination of p63 abundance or context-specific cofactors. Finally, changes in p63 isoform expression dramatically alters regulatory element activity, primarily shifting inactive elements towards a strong p63-dependent activity. Our analysis of p63-bound gene regulatory elements provides new insight into how sequence, cellular context, and other transcription factors influence p63-dependent transcription. These studies provide a framework for understanding how p63 genomic binding locally regulates transcription. Additionally, these results can be extended to investigate the influence of sequence content, genomic context, chromatin structure on the interplay between p63 isoforms and p53 family paralogs.

Dynamic enhancer landscapes in human craniofacial development

The genetic basis of human facial variation and craniofacial birth defects remains poorly understood. Distant-acting transcriptional enhancers control the fine-tuned spatiotemporal expression of genes during critical stages of craniofacial development. However, a lack of accurate maps of the genomic locations and cell type-resolved activities of craniofacial enhancers prevents their systematic exploration in human genetics studies. Here, we combine histone modification, chromatin accessibility, and gene expression profiling of human craniofacial development with single-cell analyses of the developing mouse face to define the regulatory landscape of facial development at tissue- and single cell-resolution. We provide temporal activity profiles for 14,000 human developmental craniofacial enhancers. We find that 56% of human craniofacial enhancers share chromatin accessibility in the mouse and we provide cell population- and embryonic stage-resolved predictions of their in vivo activity. Taken together, our data provide an expansive resource for genetic and developmental studies of human craniofacial development.

Ocular manifestations in Koolen-de Vries syndrome: an international study

OBJECTIVES

Koolen-de Vries Syndrome (KdVS) is a rare multisystem neurodevelopmental disorder. Ocular manifestations, including strabismus, ptosis, and hyperopia, have been reported in KdVS patients, but detailed clinical data are limited. This study aims to investigate the already known ocular malformations and their frequency while uncovering novel ocular associations.

METHODS

This was an international cross-sectional study. An anonymous questionnaire was sent to 237 KdVS patients registered in the GenIDA database. The questionnaire inquired about demographic data, ocular symptoms, findings reported by ophthalmologists, and ophthalmologic surgical interventions. The main outcome measures included ocular findings and surgical interventions.

RESULTS

Sixty-seven respondents worldwide completed the questionnaire, most (n = 53; 79%) under 18 years of age. Ophthalmologic abnormalities, noted in 79% of patients, included refractive errors (n = 35; 52.2%), strabismus (n = 23; 34.3%), amblyopia (n = 13; 19.5%), and eyelid ptosis (n = 9; 13.4%). Lacrimal disorders were present (n = 6; 9.0%), as were retinal findings (n = 7; 10.4%), including retinal hyperpigmentation or hypopigmentation (n = 4; 7.5%), Sjögren's pigment epithelial reticular dystrophy (n = 1; 1.5%), and macular chorioretinal coloboma (n = 1; 1.5%). Other manifestations included ocular surface disorders (n = 5; 7.5%), cataracts (n = 3; 4.5%), Brown syndrome (n = 1; 1.5%), glaucoma (n = 1; 1.5%), cerebral visual impairment (n = 1; 1.5%), and optic atrophy (n = 1; 1.5%). Fourteen patients (20.8%) had undergone surgical interventions.

CONCLUSIONS

KdVS is associated with various ophthalmic findings, such as amblyopia, refractive errors, strabismus, and eyelid ptosis. We describe, for the first time, a high rate of nasolacrimal disorders and retinal abnormalities consisting mainly of pigmentary findings, including a rare case of Sjögren's pigment epithelial reticular dystrophy. A comprehensive ophthalmic evaluation is therefore recommended for all KdVS patients at initial diagnosis or at 4-6 months of age for diagnosed newborns.

Epigenetic priming of an epithelial enhancer by p63 and CTCF controls expression of a skin-restricted gene XP33

The transcription factor p63 is a renowned master regulator of gene expression of stratified epithelia. While multiple proteins have been identified as p63 bona fide targets, little is known about non-coding RNAs (ncRNAs) whose transcription is controlled by p63. Here, we describe a skin-specific non-coding RNA XP33 as a novel target of p63. XP33 levels are increased during keratinocyte differentiation in vitro, while its depletion results in decreased expression of late cornified gene LCE2D. By using publicly available multi-omics data, we show that CTCF and p63 establish an epithelial enhancer to prime XP33 transcription in a tissue-restricted manner. XP33 promoter and enhancer form a chromatin loop exclusively in keratinocytes but not in other cell types. Moreover, the XP33 enhancer is occupied by differentiation-specific factors that control XP33 transcription. Altogether, we identify a tissue-specific non-coding RNA whose expression is epigenetically regulated by p63 and CTCF.

ΔNp63α facilitates proliferation and migration, and modulates the chromatin landscape in intrahepatic cholangiocarcinoma cells

p63 plays a crucial role in epithelia-originating tumours; however, its role in intrahepatic cholangiocarcinoma (iCCA) has not been completely explored. Our study revealed the oncogenic properties of p63 in iCCA and identified the major expressed isoform as ΔNp63α. We collected iCCA clinical data from The Cancer Genome Atlas database and analyzed p63 expression in iCCA tissue samples. We further established genetically modified iCCA cell lines in which p63 was overexpressed or knocked down to study the protein function/function of p63 in iCCA. We found that cells overexpressing p63, but not p63 knockdown counterparts, displayed increased proliferation, migration, and invasion. Transcriptome analysis showed that p63 altered the iCCA transcriptome, particularly by affecting cell adhesion-related genes. Moreover, chromatin accessibility decreased at p63 target sites when p63 binding was lost and increased when p63 binding was gained. The majority of the p63 bound sites were located in the distal intergenic regions and showed strong enhancer marks; however, active histone modifications around the Transcription Start Site changed as p63 expression changed. We also detected an interaction between p63 and the chromatin structural protein YY1. Taken together, our results suggest an oncogenic role for p63 in iCCA.

ΔNp63 silencing, DNA methylation shifts, and epithelial-mesenchymal transition resulted from TAp63 genome editing in squamous cell carcinoma

TP63 (p63) is strongly expressed in lower-grade carcinomas of the head and neck, skin, breast, and urothelium to maintain a well-differentiated phenotype. TP63 has two transcription start sites at exons 1 and 3' that produce TAp63 and ΔNp63 isoforms, respectively. The major protein, ΔNp63α, epigenetically activates genes essential for epidermal/craniofacial differentiation, including ΔNp63 itself. To examine the specific role of weakly expressed TAp63, we disrupted exon 1 using CRISPR-Cas9 homology-directed repair in a head and neck squamous cell carcinoma (SCC) line. Surprisingly, TAp63 knockout cells having either monoallelic GFP cassette insertion paired with a frameshift deletion allele or biallelic GFP cassette insertion exhibited ΔNp63 silencing. Loss of keratinocyte-specific gene expression, switching of intermediate filament genes from KRT(s) to VIM, and suppression of cell-cell and cell-matrix adhesion components indicated the core events of epithelial-mesenchymal transition. Many of the positively and negatively affected genes, including ΔNp63, displayed local DNA methylation changes. Furthermore, ΔNp63 expression was partially rescued by transfection of the TAp63 knockout cells with TAp63α and application of DNA methyltransferase inhibitor zebularine. These results suggest that TAp63, a minor part of the TP63 gene, may be involved in the auto-activation mechanism of ΔNp63 by which the keratinocyte-specific epigenome is maintained in SCC.

An integrated genomic approach identifies follistatin as a target of the p63-epidermal growth factor receptor oncogenic network in head and neck squamous cell carcinoma

Although numerous putative oncogenes have been associated with the etiology of head and neck squamous cell carcinoma (HNSCC), the mechanisms by which these oncogenes and their downstream targets mediate tumor progression have not been fully elucidated. We performed an integrative analysis to identify a crucial set of targets of the oncogenic transcription factor p63 that are common across multiple transcriptomic datasets obtained from HNSCC patients, and representative cell line models. Notably, our analysis revealed FST which encodes follistatin, a secreted glycoprotein that inhibits the transforming growth factor TGFβ/activin signaling pathways, to be a direct transcriptional target of p63. In addition, we found that FST expression is also driven by epidermal growth factor receptor EGFR signaling, thus mediating a functional link between the TGF-β and EGFR pathways. We show through loss- and gain-of-function studies that FST predominantly imparts a tumor-growth and migratory phenotype in HNSCC cells. Furthermore, analysis of single-cell RNA sequencing data from HNSCC patients unveiled cancer cells as the dominant source of FST within the tumor microenvironment and exposed a correlation between the expression of FST and its regulators with immune infiltrates. We propose FST as a prognostic biomarker for patient survival and a compelling candidate mediating the broad effects of p63 on the tumor and its associated microenvironment.

Cell Type- and Tissue-specific Enhancers in Craniofacial Development

The genetic basis of craniofacial birth defects and general variation in human facial shape remains poorly understood. Distant-acting transcriptional enhancers are a major category of non-coding genome function and have been shown to control the fine-tuned spatiotemporal expression of genes during critical stages of craniofacial development1-3. However, a lack of accurate maps of the genomic location and cell type-specific in vivo activities of all craniofacial enhancers prevents their systematic exploration in human genetics studies. Here, we combined histone modification and chromatin accessibility profiling from different stages of human craniofacial development with single-cell analyses of the developing mouse face to create a comprehensive catalogue of the regulatory landscape of facial development at tissue- and single cell-resolution. In total, we identified approximately 14,000 enhancers across seven developmental stages from weeks 4 through 8 of human embryonic face development. We used transgenic mouse reporter assays to determine the in vivo activity patterns of human face enhancers predicted from these data. Across 16 in vivo validated human enhancers, we observed a rich diversity of craniofacial subregions in which these enhancers are active in vivo. To annotate the cell type specificities of human-mouse conserved enhancers, we performed single-cell RNA-seq and single-nucleus ATAC-seq of mouse craniofacial tissues from embryonic days e11.5 to e15.5. By integrating these data across species, we find that the majority (56%) of human craniofacial enhancers are functionally conserved in mice, providing cell type- and embryonic stage-resolved predictions of their in vivo activity profiles. Using retrospective analysis of known craniofacial enhancers in combination with single cell-resolved transgenic reporter assays, we demonstrate the utility of these data for predicting the in vivo cell type specificity of enhancers. Taken together, our data provide an expansive resource for genetic and developmental studies of human craniofacial development.

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.

DNA methylation differences in monozygotic twins with Van der Woude syndrome

Introduction

Van der Woude Syndrome (VWS) is an autosomal dominant disorder responsible for 2% of all syndromic orofacial clefts (OFCs) with IRF6 being the primary causal gene (70%). Cases may present with lip pits and either cleft lip, cleft lip with cleft palate, or cleft palate, with marked phenotypic discordance even among individuals carrying the same mutation. This suggests that genetic or epigenetic modifiers may play additional roles in the syndrome's etiology and variability in expression. We report the first DNA methylation profiling of 2 pairs of monozygotic twins with VWS. Our goal is to explore epigenetic contributions to VWS etiology and variable phenotypic expressivity by comparing DNAm profiles in both twin pairs. While the mutations that cause VWS in these twins are known, the additional mechanism behind their phenotypic risk and variability in expression remains unclear.

Methods

We generated whole genome DNAm data for both twin pairs. Differentially methylated positions (DMPs) were selected based on: (1) a coefficient of variation in DNAm levels in unaffected individuals < 20%, and (2) intra-twin pair absolute difference in DNAm levels >5% (delta beta > | 0.05|). We then divided the DMPs in two subgroups for each twin pair for further analysis: (1) higher methylation levels in twin A (Twin A > Twin B); and (2) higher methylation levels in twin B (Twin B >Twin A).

Results and Discussion

Gene ontology analysis revealed a list of enriched genes that showed significant differential DNAm, including clef-associated genes. Among the cleft-associated genes, TP63 was the most significant hit (p=7.82E-12). Both twin pairs presented differential DNAm levels in CpG sites in/near TP63 (Twin 1A > Twin 1B and Twin 2A < Twin 2B). The genes TP63 and IRF6 function in a biological regulatory loop to coordinate epithelial proliferation and differentiation in a process that is critical for palatal fusion. The effects of the causal mutations in IRF6 can be further impacted by epigenetic dysregulation of IRF6 itself, or genes in its pathway. Our data shows evidence that changes in DNAm is a plausible mechanism that can lead to markedly distinct phenotypes, even among individuals carrying the same mutation.

Analysis of candidate genes for cleft lip ± cleft palate using murine single-cell expression data

Introduction: Cleft lip ± cleft palate (CL/P) is one of the most common birth defects. Although research has identified multiple genetic risk loci for different types of CL/P (i.e., syndromic or non-syndromic forms), determining the respective causal genes and understanding the relevant functional networks remain challenging. The recent introduction of single-cell RNA sequencing (scRNA-seq) has provided novel opportunities to study gene expression patterns at cellular resolution. The aims of our study were to: (i) aggregate available scRNA-seq data from embryonic mice and provide this as a resource for the craniofacial community; and (ii) demonstrate the value of these data in terms of the investigation of the gene expression patterns of CL/P candidate genes. Methods and Results: First, two published scRNA-seq data sets from embryonic mice were re-processed, i.e., data representing the murine time period of craniofacial development: (i) facial data from embryonic day (E) E11.5; and (ii) whole embryo data from E9.5-E13.5 from the Mouse Organogenesis Cell Atlas (MOCA). Marker gene expression analyses demonstrated that at E11.5, the facial data were a high-resolution representation of the MOCA data. Using CL/P candidate gene lists, distinct groups of genes with specific expression patterns were identified. Among others we identified that a co-expression network including Irf6, Grhl3 and Tfap2a in the periderm, while it was limited to Irf6 and Tfap2a in palatal epithelia, cells of the ectodermal surface, and basal cells at the fusion zone. The analyses also demonstrated that additional CL/P candidate genes (e.g., Tpm1, Arid3b, Ctnnd1, and Wnt3) were exclusively expressed in Irf6+ facial epithelial cells (i.e., as opposed to Irf6- epithelial cells). The MOCA data set was finally used to investigate differences in expression profiles for candidate genes underlying different types of CL/P. These analyses showed that syndromic CL/P genes (syCL/P) were expressed in significantly more cell types than non-syndromic CL/P candidate genes (nsCL/P). Discussion: The present study illustrates how scRNA-seq data can empower research on craniofacial development and disease.

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

A cis-regulatory lexicon of DNA motif combinations mediating cell-type-specific gene regulation

Gene expression is controlled by transcription factors (TFs) that bind cognate DNA motif sequences in cis-regulatory elements (CREs). The combinations of DNA motifs acting within homeostasis and disease, however, are unclear. Gene expression, chromatin accessibility, TF footprinting, and H3K27ac-dependent DNA looping data were generated and a random-forest-based model was applied to identify 7,531 cell-type-specific cis-regulatory modules (CRMs) across 15 diploid human cell types. A co-enrichment framework within CRMs nominated 838 cell-type-specific, recurrent heterotypic DNA motif combinations (DMCs), which were functionally validated using massively parallel reporter assays. Cancer cells engaged DMCs linked to neoplasia-enabling processes operative in normal cells while also activating new DMCs only seen in the neoplastic state. This integrative approach identifies cell-type-specific cis-regulatory combinatorial DNA motifs in diverse normal and diseased human cells and represents a general framework for deciphering cis-regulatory sequence logic in gene regulation.

Novel Candidate Genes for Non-Syndromic Tooth Agenesis Identified Using Targeted Next-Generation Sequencing

Non-syndromic tooth agenesis (ns-TA) is one of the most common dental anomalies characterized by the congenital absence of at least one permanent tooth (excluding third molars). Regarding the essential role of genetic factors in ns-TA aetiology, the present study aimed to identify novel pathogenic variants underlying hypodontia and oligodontia. In a group of 65 ns-TA patients and 127 healthy individuals from the genetically homogenous Polish population, the coding sequences of 423 candidate genes were screened using targeted next-generation sequencing. Pathogenic and likely pathogenic variants were identified in 37 (56.92%) patients, including eight nucleotide alternations of genes not previously implicated in ns-TA (CHD7, CREBBP, EVC, LEF1, ROR2, TBX22 and TP63). However, since only single variants were detected, future research is required to confirm and fully understand their role in the aetiology of ns-TA. Additionally, our results support the importance of already known ns-TA candidate genes (AXIN2, EDA, EDAR, IRF6, LAMA3, LRP6, MSX1, PAX9 and WNT10A) and provide additional evidence that ns-TA might be an oligogenic condition involving the cumulative effect of rare variants in two or more distinct genes.

Targeted re-sequencing on 1p22 among non-syndromic orofacial clefts from Han Chinese population

Rs560426 at 1p22 was proved to be associated with NSCL/P (non-syndromic cleft lip with or without the palate) in several populations, including Han Chinese population. Here, we conducted a deep sequencing around rs560426 to locate more susceptibility variants in this region. In total, 2,293 NSCL/P cases and 3,235 normal controls were recruited. After sequencing, association analysis was performed. Western blot, RT-qPCR, HE, immunofluorescence staining, and RNA sequencing were conducted for functional analyses of the selected variants. Association analysis indicated that rs77179923 was the only SNP associated with NSCLP specifically (p = 4.70E-04, OR = 1.84), and rs12071152 was uniquely associated with LCLO (p = 4.00E-04, OR = 1.30, 95%CI: 1.12–1.51). Moreover, de novo harmful rare variant NM_004815.3, NP_004806.3; c.1652G>C, p.R551T in ARHGAP29 resulted in a decreased expression level of ARHGAP29, which in turn affected NSCL/P-related biological processes; however, no overt cleft palate (CP) phenotype was observed. In conclusion, rs12071152 was a new susceptible variant, which is specifically associated with LCLO among the Han Chinese population. Allele A of it could increase the risk of having a cleft baby. Rs77179923 and rare variant NM_004815.3, NP_004806.3; c.1652G>C, p.R551T at 1p22 were both associated with NSCLP among the Han Chinese population. However, this missense variation contributes to no overt CP phenotype due to dosage insufficiency or compensation from other genes.

Distinct interactors define the p63 transcriptional signature in epithelial development or cancer

The TP63 is an indispensable transcription factor for development and homeostasis of epithelia and its derived glandular tissue. It is also involved in female germline cell quality control, muscle and thymus development. It is expressed as multiple isoforms transcribed by two independent promoters, in addition to alternative splicing occurring at the mRNA 3′-UTR. Expression of the TP63 gene, specifically the amino-deleted p63 isoform, ΔNp63, is required to regulate numerous biological activities, including lineage specification, self-renewal capacity of epithelial stem cells, proliferation/expansion of basal keratinocytes, differentiation of stratified epithelia. In cancer, ΔNp63 is implicated in squamous cancers pathogenesis of different origin including skin, head and neck and lung and in sustaining self-renewal of cancer stem cells. How this transcription factor can control such a diverse set of biological pathways is central to the understanding of the molecular mechanisms through which p63 acquires oncogenic activity, profoundly changing its down-stream transcriptional signature. Here, we highlight how different proteins interacting with p63 allow it to regulate the transcription of several central genes. The interacting proteins include transcription factors/regulators, epigenetic modifiers, and post-transcriptional modifiers. Moreover, as p63 depends on its interactome, we discuss the hypothesis to target the protein interactors to directly affect p63 oncogenic activities and p63-related diseases.

Target sequencing reveals the association between variants in VAX1 and NSCL/P in Chinese population

OBJECTIVE

A significant genetic association between rs7078160 in VAX1 and NSCL/P has been established through genome-wide association studies (GWAS), and we previously replicated the association in the Chinese population. The critical issue in the post-GWAS era is to identify functional variations that have a real impact on disease in the susceptible regions highlighted by GWAS. This study aimed to elucidate functional variants in VAX1 fully.

MATERIALS AND METHODS

Firstly, target sequencing was performed on 159 NSCL/P patients, followed by association analysis to discover disease-associated single nucleotide polymorphisms (SNPs); we then replicated the findings using a larger sample (1626 cases, 2255 controls) and investigated how candidate SNPs affect disease occurrence using extensive annotation databases. Additionally, we compared the genetic profiles of NSCL/P subtypes.

RESULTS

In this study, 6 SNPs in VAX1 were identified to be associated with NSCL/P in the Western Han Chinese population. Five of them were predicted to influence transcriptional factor biding ability and were expression quantitative trait loci (eQTLs) of nearby genes in multiple tissues.

CONCLUSION

The previously reported association between rs7078160 and NSCL/P was successfully replicated. Moreover, our findings firstly revealed that 5 SNPs in VAX1 are associated with NSCL/P in the Western Han Chinese population.

BrewerIX enables allelic expression analysis of imprinted and X-linked genes from bulk and single-cell transcriptomes

Genomic imprinting and X chromosome inactivation (XCI) are two prototypical epigenetic mechanisms whereby a set of genes is expressed mono-allelically in order to fine-tune their expression levels. Defects in genomic imprinting have been observed in several neurodevelopmental disorders, in a wide range of tumours and in induced pluripotent stem cells (iPSCs). Single Nucleotide Variants (SNVs) are readily detectable by RNA-sequencing allowing the determination of whether imprinted or X-linked genes are aberrantly expressed from both alleles, although standardised analysis methods are still missing. We have developed a tool, named BrewerIX, that provides comprehensive information about the allelic expression of a large, manually-curated set of imprinted and X-linked genes. BrewerIX does not require programming skills, runs on a standard personal computer, and can analyze both bulk and single-cell transcriptomes of human and mouse cells directly from raw sequencing data. BrewerIX confirmed previous observations regarding the bi-allelic expression of some imprinted genes in naive pluripotent cells and extended them to preimplantation embryos. BrewerIX also identified misregulated imprinted genes in breast cancer cells and in human organoids and identified genes escaping XCI in human somatic cells. We believe BrewerIX will be useful for the study of genomic imprinting and XCI during development and reprogramming, and for detecting aberrations in cancer, iPSCs and organoids. Due to its ease of use to non-computational biologists, its implementation could become standard practice during sample assessment, thus raising the robustness and reproducibility of future studies.

Targeted Re-Sequencing of the 2p21 Locus Identifies Non-Syndromic Cleft Lip Only Novel Susceptibility Gene ZFP36L2

rs7590268 present on the 2p21 locus was identified to be associated with non-syndromic cleft lip with or without cleft palate (NSCL/P) in several populations, including the Chinese Han population, indicating that 2p21 was a susceptibility locus for NSCL/P. However, previous studies have only identified common single-nucleotide polymorphism (SNP) within the THADA gene, neglecting the rare variants and other genes in 2p21; thus, this study was designed to investigate additional variants and novel susceptibility genes in 2p21. A total of 159 NSCL/P patients and 542 controls were recruited in the discovery phase, whereas 1830 NSCL/P patients and 2,436 controls were recruited in the replication phase. After targeted region sequencing, we performed association and burden analyses for the common and rare variants, respectively. Furthermore, RNA-seq, proliferation assay and cell cycle analysis were performed to clarify the possible function of the candidate gene ZFP36L2. Association analysis showed that four SNPs were specifically associated with non-syndromic cleft lip only (NSCLO) and two SNPs were associated with both NSCLO and NSCL/P. Burden analysis indicated that ZFP36L2 was associated with NSCLO (p = .0489, OR = 2.41, 95% CI: 0.98–5.90). Moreover, SNPs in the ZFP36L2 targeted gene JUP were also associated with NSCLO. ZFP36L2 also inhibited cell proliferation and induced G2 phase arrest in the GMSM-K cell line. Therefore, we proposed that ZFP36L2 is a novel susceptibility gene of NSCLO in the 2p21 locus, which could lead to NSCLO by modulating cell proliferation and cycle.

Mouse models in palate development and orofacial cleft research: Understanding the crucial role and regulation of epithelial integrity in facial and palate morphogenesis

Cleft lip and cleft palate are common birth defects resulting from genetic and/or environmental perturbations of facial development in utero. Facial morphogenesis commences during early embryogenesis, with cranial neural crest cells interacting with the surface ectoderm to form initially partly separate facial primordia consisting of the medial and lateral nasal prominences, and paired maxillary and mandibular processes. As these facial primordia grow around the primitive oral cavity and merge toward the ventral midline, the surface ectoderm undergoes a critical differentiation step to form an outer layer of flattened and tightly connected periderm cells with a non-stick apical surface that prevents epithelial adhesion. Formation of the upper lip and palate requires spatiotemporally regulated inter-epithelial adhesions and subsequent dissolution of the intervening epithelial seam between the maxillary and medial/lateral nasal processes and between the palatal shelves. Proper regulation of epithelial integrity plays a paramount role during human facial development, as mutations in genes encoding epithelial adhesion molecules and their regulators have been associated with syndromic and non-syndromic orofacial clefts. In this chapter, we summarize mouse genetic studies that have been instrumental in unraveling the mechanisms regulating epithelial integrity and periderm differentiation during facial and palate development. Since proper epithelial integrity also plays crucial roles in wound healing and cancer, understanding the mechanisms regulating epithelial integrity during facial development have direct implications for improvement in clinical care of craniofacial patients.

Programmed Cell Death Not as Sledgehammer but as Chisel: Apoptosis in Normal and Abnormal Craniofacial Patterning and Development

Coordination of craniofacial development involves an complex, intricate, genetically controlled and tightly regulated spatiotemporal series of reciprocal inductive and responsive interactions among the embryonic cephalic epithelia (both endodermal and ectodermal) and the cephalic mesenchyme — particularly the cranial neural crest (CNC). The coordinated regulation of these interactions is critical both ontogenetically and evolutionarily, and the clinical importance and mechanistic sensitivity to perturbation of this developmental system is reflected by the fact that one-third of all human congenital malformations affect the head and face. Here, we focus on one element of this elaborate process, apoptotic cell death, and its role in normal and abnormal craniofacial development. We highlight four themes in the temporospatial elaboration of craniofacial apoptosis during development, namely its occurrence at (1) positions of epithelial-epithelial apposition, (2) within intra-epithelial morphogenesis, (3) during epithelial compartmentalization, and (4) with CNC metameric organization. Using the genetic perturbation of Satb2, Pbx1/2, Fgf8, and Foxg1 as exemplars, we examine the role of apoptosis in the elaboration of jaw modules, the evolution and elaboration of the lambdoidal junction, the developmental integration at the mandibular arch hinge, and the control of upper jaw identity, patterning and development. Lastly, we posit that apoptosis uniquely acts during craniofacial development to control patterning cues emanating from core organizing centres.

Integrative approaches generate insights into the architecture of non-syndromic cleft lip with or without cleft palate

Non-syndromic cleft lip with or without cleft palate (nsCL/P) is a common congenital facial malformation with a multifactorial etiology. Genome-wide association studies (GWASs) have identified multiple genetic risk loci. However, functional interpretation of these loci is hampered by the underrepresentation in public resources of systematic functional maps representative of human embryonic facial development. To generate novel insights into the etiology of nsCL/P, we leveraged published GWAS data on nsCL/P as well as available chromatin modification and expression data on mid-facial development. Our analyses identified five novel risk loci, prioritized candidate target genes within associated regions, and highlighted distinct pathways. Furthermore, the results suggest the presence of distinct regulatory effects of nsCL/P risk variants throughout mid-facial development and shed light on its regulatory architecture. Our integrated data provide a platform to advance hypothesis-driven molecular investigations of nsCL/P and other human facial defects.

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.

Endogenous Mechanisms of Craniomaxillofacial Repair: Toward Novel Regenerative Therapies

In the fields of oral and craniomaxillofacial surgery, regeneration of multiple tissue types-including bone, skin, teeth, and mucosal soft tissue-is often a desired outcome. However, limited endogenous capacity for regeneration, as well as predisposition of many tissues to fibrotic healing, may prevent recovery of normal form and function for patients. Recent basic science research has advanced our understanding of molecular and cellular pathways of repair in the oral/craniofacial region and how these are influenced by local microenvironment and embryonic origin. Here, we review the current state of knowledge in oral and craniomaxillofacial tissue repair/regeneration in four key areas: bone (in the context of calvarial defects and mandibular regeneration during distraction osteogenesis); skin (in the context of cleft lip/palate surgery); oral mucosa (in the context of minimally scarring repair of mucosal injuries); and teeth (in the context of dental disease/decay). These represent four distinct healing processes and outcomes. We will discuss both divergent and conserved pathways of repair in these contexts, with an eye toward fundamental mechanisms of regeneration vs. fibrosis as well as translational research directions. Ultimately, this knowledge can be leveraged to develop new cell-based and molecular treatment strategies to encourage bone and soft tissue regeneration in oral and craniomaxillofacial surgery.

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

Lyophilized powder of mesenchymal stem cell supernatant attenuates acute lung injury through the IL-6-p-STAT3-p63-JAG2 pathway

BACKGROUND

Acute lung injury (ALI) and acute respiratory distress syndrome (ARDS) are syndromes of acute respiratory failure with extremely high mortality and few effective treatments. Mesenchymal stem cells (MSCs) may reportedly contribute to tissue repair in ALI and ARDS. However, applications of MSCs have been restricted due to safety considerations and limitations in terms of large-scale production and industrial delivery. Alternatively, the MSC secretome has been considered promising for use in therapeutic approaches and has been advanced in pre-clinical and clinical trials. Furthermore, the MSC secretome can be freeze-dried into a stable and ready-to-use supernatant lyophilized powder (SLP) form. Currently, there are no studies on the role of MSC SLP in ALI.

METHODS

Intratracheal bleomycin was used to induce ALI in mice, and intratracheal MSC SLP was administered as a treatment. Histopathological assessment was performed by hematoxylin and eosin, immunohistochemistry, and immunofluorescence staining. Apoptosis, inflammatory infiltration, immunological cell counts, cytokine levels, and mRNA- and protein-expression levels of relevant targets were measured by performing terminal deoxynucleotidyl transferase dUTP nick-end labeling assays, determining total cell and protein levels in bronchoalveolar lavage fluids, flow cytometry, multiple cytokine-detection techniques, and reverse transcriptase-quantitative polymerase chain reaction and western blot analysis, respectively.

RESULTS

We found that intratracheal MSC SLP considerably promoted cell survival, inhibited epithelial cell apoptosis, attenuated inflammatory cell recruitment, and reversed immunological imbalances induced by bleomycin. MSC SLP inhibited the interleukin 6-phosphorylated signal transducer and activator of transcription signaling pathway to activate tumor protein 63-jagged 2 signaling in basal cells, suppress T helper 17 cell differentiation, promote p63⁺ cell proliferation and lung damage repair, and attenuate inflammatory responses.

CONCLUSIONS

MSC SLP ameliorated ALI by activating p63 and promoting p63⁺ cell proliferation and the repair of damaged epithelial cells. The findings of this study also shed insight into ALI pathogenesis and imply that MSC SLP shows considerable therapeutic promise for treating ALI and ARDS.

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.

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.

Cranial Neural Crest Cells and Their Role in the Pathogenesis of Craniofacial Anomalies and Coronal Craniosynostosis

Craniofacial anomalies are among the most common of birth defects. The pathogenesis of craniofacial anomalies frequently involves defects in the migration, proliferation, and fate of neural crest cells destined for the craniofacial skeleton. Genetic mutations causing deficient cranial neural crest migration and proliferation can result in Treacher Collins syndrome, Pierre Robin sequence, and cleft palate. Defects in post-migratory neural crest cells can result in pre- or post-ossification defects in the developing craniofacial skeleton and craniosynostosis (premature fusion of cranial bones/cranial sutures). The coronal suture is the most frequently fused suture in craniosynostosis syndromes. It exists as a biological boundary between the neural crest-derived frontal bone and paraxial mesoderm-derived parietal bone. The objective of this review is to frame our current understanding of neural crest cells in craniofacial development, craniofacial anomalies, and the pathogenesis of coronal craniosynostosis. We will also discuss novel approaches for advancing our knowledge and developing prevention and/or treatment strategies for craniofacial tissue regeneration and craniosynostosis.

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.