AP-2α and AP-2β cooperatively function in the craniofacial surface ectoderm to regulate chromatin and gene expression dynamics during facial development

AP-2α/AP-2β Transcription Factors Are Key Regulators of Epidermal Homeostasis

AP-2 transcription factors regulate ectodermal development, but their roles in epidermal homeostasis in adult skin are unknown. We find that AP-2α is the predominant AP-2 family member in adult epidermis, followed by AP-2β. Through inactivation of AP-2α, AP-2β, or both in keratinocytes, we assessed the effects of a gradient of epidermal AP-2 activity on skin function. We find that (i) loss of AP-2β in keratinocytes is compensated for by AP-2α, (ii) loss of AP-2α impairs terminal keratinocyte differentiation and hair morphogenesis, and (iii) the combined loss of AP-2α/AP-2β results in more severe skin and hair abnormalities. Keratinocyte differentiation defects precede progressive neutrophilic skin inflammation. Inducible inactivation of AP-2α/AP-2β in the adult phenocopies these manifestations. Transcriptomic analyses of epidermis lacking AP-2α or AP-2α/AP-2β in keratinocytes demonstrate a terminal keratinocyte differentiation defect with upregulation of alarmin keratins and of several immune pathway regulators. Moreover, our analyses suggest a key role of reduced AP-2α-dependent gene expression of CXCL14 and the keratin 15 gene K15 as an early pathogenic event toward the manifestation of skin inflammation. Thus, AP-2α and AP-2β are critical regulators of epidermal homeostasis in adult skin.

HOXA-AS2 Epigenetically Inhibits HBV Transcription by Recruiting the MTA1-HDAC1/2 Deacetylase Complex to cccDNA Minichromosome

Persistent transcription of HBV covalently closed circular DNA (cccDNA) is critical for chronic HBV infection. Silencing cccDNA transcription through epigenetic mechanisms offers an effective strategy to control HBV. Long non-coding RNAs (lncRNAs), as important epigenetic regulators, have an unclear role in cccDNA transcription regulation. In this study, lncRNA sequencing (lncRNA seq) is conducted on five pairs of HBV-positive and HBV-negative liver tissue. Through analysis, HOXA-AS2 (HOXA cluster antisense RNA 2) is identified as a significantly upregulated lncRNA in HBV-infected livers. Further experiments demonstrate that HBV DNA polymerase (DNA pol) induces HOXA-AS2 after establishing persistent high-level HBV replication. Functional studies reveal that HOXA-AS2 physically binds to cccDNA and significantly inhibits its transcription. Mechanistically, HOXA-AS2 recruits the MTA1-HDAC1/2 deacetylase complex to cccDNA minichromosome by physically interacting with metastasis associated 1 (MTA1) subunit, resulting in reduced acetylation of histone H3 at lysine 9 (H3K9ac) and lysine 27 (H3K27ac) associated with cccDNA and subsequently suppressing cccDNA transcription. Altogether, the study reveals a mechanism to self-limit HBV replication, wherein the upregulation of lncRNA HOXA-AS2, induced by HBV DNA pol, can epigenetically suppress cccDNA transcription.

Whole genome sequencing identifies associations for nonsyndromic sagittal craniosynostosis with the intergenic region of BMP2 and noncoding RNA gene LINC01428

Craniosynostosis (CS) is a major birth defect resulting from premature fusion of cranial sutures. Nonsyndromic CS occurs more frequently than syndromic CS, with sagittal nonsyndromic craniosynostosis (sNCS) presenting as the most common CS phenotype. Previous genome-wide association and targeted sequencing analyses of sNCS have identified multiple associated loci, with the strongest association on chromosome 20. Herein, we report the first whole-genome sequencing study of sNCS using 63 proband-parent trios. Sequencing data for these trios were analyzed using the transmission disequilibrium test (TDT) and rare variant TDT (rvTDT) to identify high-risk rare gene variants. Sequencing data were also examined for copy number variants (CNVs) and de novo variants. TDT analysis identified a highly significant locus at 20p12.3, localized to the intergenic region between BMP2 and the noncoding RNA gene LINC01428. Three variants (rs6054763, rs6054764, rs932517) were identified as potential causal variants due to their probability of being transcription factor binding sites, deleterious combined annotation dependent depletion scores, and high minor allele enrichment in probands. Morphometric analysis of cranial vault shape in an unaffected cohort validated the effect of these three single nucleotide variants (SNVs) on dolichocephaly. No genome-wide significant rare variants, de novo loci, or CNVs were identified. Future efforts to identify risk variants for sNCS should include sequencing of larger and more diverse population samples and increased omics analyses, such as RNA-seq and ATAC-seq.

[Role of Fibroblast Growth Factor 7 in Craniomaxillofacial Development]

Craniomaxillofacial development involves a series of highly ordered temporal-spatial cellular differentiation processes in which a variety of cell signaling factors, such as fibroblast growth factors, play important regulatory roles. As a classic fibroblast growth factor, fibroblast growth factor 7 (FGF7) serves a wide range of regulatory functions. Previous studies have demonstrated that FGF7 regulates the proliferation and migration of epithelial cells, protects them, and promotes their repair. Furthermore, recent findings indicate that epithelial cells are not the only ones subjected to the broad and powerful regulatory capacity of FGF7. It has potential effects on skeletal system development as well. In addition, FGF7 plays an important role in the development of craniomaxillofacial organs, such as the palate, the eyes, and the teeth. Nonetheless, the role of FGF7 in oral craniomaxillofacial development needs to be further elucidated. In this paper, we summarized the published research on the role of FGF7 in oral craniomaxillofacial development to demonstrate the overall understanding of FGF7 and its potential functions in oral craniomaxillofacial development.

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.

Temporally resolved early BMP-driven transcriptional cascade during human amnion specification

Amniogenesis, a process critical for continuation of healthy pregnancy, is triggered in a collection of pluripotent epiblast cells as the human embryo implants. Previous studies have established that BMP signaling is a major driver of this lineage specifying process, but the downstream BMP-dependent transcriptional networks that lead to successful amniogenesis remain to be identified. This is, in part, due to the current lack of a robust and reproducible model system that enables mechanistic investigations exclusively into amniogenesis. Here, we developed an improved model of early amnion specification, using a human pluripotent stem cell-based platform in which the activation of BMP signaling is controlled and synchronous. Uniform amniogenesis is seen within 48 hours after BMP activation, and the resulting cells share transcriptomic characteristics with amnion cells of a gastrulating human embryo. Using detailed time-course transcriptomic analyses, we established a previously uncharacterized BMP-dependent amniotic transcriptional cascade, and identified markers that represent five distinct stages of amnion fate specification; the expression of selected markers was validated in early post-implantation macaque embryos. Moreover, a cohort of factors that could potentially control specific stages of amniogenesis was identified, including the transcription factor TFAP2A. Functionally, we determined that, once amniogenesis is triggered by the BMP pathway, TFAP2A controls the progression of amniogenesis. This work presents a temporally resolved transcriptomic resource for several previously uncharacterized amniogenesis states and demonstrates a critical intermediate role for TFAP2A during amnion fate specification.

Potential prognosis and immunotherapy predictor TFAP2A in pan-cancer

BACKGROUND

TFAP2A is critical in regulating the expression of various genes, affecting various biological processes and driving tumorigenesis and tumor development. However, the significance of TFAP2A in carcinogenesis processes remains obscure.

METHODS

In our study, we explored multiple databases including TCGA, GTEx, HPA, cBioPortal, TCIA, and other well-established databases for further analysis to expound TFAP2A expression, genetic alternations, and their relationship with the prognosis and cellular signaling network alternations. GO term and KEGG pathway enrichment analysis as well as GSEA were conducted to examine the common functions of TFAP2A. RT-qPCR, Western Blot and Dual Luciferase Reporter assay were employed to perform experimental validation.

RESULTS

TFAP2A mRNA expression level was upregulated and its genetic alternations were frequently present in most cancer types. The enrichment analysis results prompted us to investigate the changes in the tumor immune microenvironment further. We discovered that the expression of TFAP2A was significantly associated with the expression of immune checkpoint genes, immune subtypes, ESTIMATE scores, tumor-infiltrating immune cells, and the possible role of TFAP2A in predicting immunotherapy efficacy. In addition, high TFAP2A expression significantly correlated with several ICP genes, and promoted the expression of PD-L1 on mRNA and protein levels through regulating its expression at the transcriptional level. TFAP2A protein level was upregulated in fresh colon tumor tissue samples compared to that in the adjacent normal tissues, which essentially positively correlated with the expression of PD-L1.

CONCLUSIONS

Our study suggests that targeting TFAP2A may provide a novel and effective strategy for cancer treatment.

TFAP2 paralogs regulate midfacial development in part through a conserved ALX genetic pathway

Cranial neural crest development is governed by positional gene regulatory networks (GRNs). Fine-tuning of the GRN components underlies facial shape variation, yet how those networks in the midface are connected and activated remain poorly understood. Here, we show that concerted inactivation of Tfap2a and Tfap2b in the murine neural crest, even during the late migratory phase, results in a midfacial cleft and skeletal abnormalities. Bulk and single-cell RNA-seq profiling reveal that loss of both TFAP2 family members dysregulates numerous midface GRN components involved in midface morphogenesis, patterning and differentiation. Notably, Alx1, Alx3 and Alx4 (ALX) transcript levels are reduced, whereas ChIP-seq analyses suggest TFAP2 family members directly and positively regulate ALX gene expression. Tfap2a, Tfap2b and ALX co-expression in midfacial neural crest cells of both mouse and zebrafish implies conservation of this regulatory axis across vertebrates. Consistent with this notion, tfap2a zebrafish mutants present with abnormal alx3 expression patterns, Tfap2a binds ALX loci and tfap2a-alx3 genetic interactions are observed. Together, these data demonstrate TFAP2 paralogs regulate vertebrate midfacial development in part by activating expression of ALX transcription factor genes.

Computational identification of key transcription factors for embryonic and postnatal Sox2+ dental epithelial stem cell

While many reptiles can replace their tooth throughout life, human loss the tooth replacement capability after formation of the permanent teeth. It was thought that the difference in tooth regeneration capability depends on the persistence of a specialized dental epithelial structure, the dental lamina that contains dental epithelial stem cells (DESC). Currently, we know very little about DESC such as what genes are expressed or its chromatin accessibility profile. Multiple markers of DESC have been proposed such as Sox2 and Lgr5 . Few single cell RNA-seq experiments have been performed previously, but no obvious DESC cluster was identified in these scRNA-seq datasets, possible due to that the expression level of DESC markers such as Sox2 and Lgr5 is too low or the percentage of DESC is too low in whole tooth. We utilize a mouse line Sox2-GFP to enrich Sox2+ DESC and use Smart-Seq2 protocol and ATAC-seq protocol to generate transcriptome profile and chromatin accessibility profile of P2 Sox2+ DESC. Additionally, we generate transcriptome profile and chromatin accessibility profile of E11.5 Sox2+ dental lamina cells. With transcriptome profile and chromatin accessibility profile, we systematically identify potential key transcription factors for E11.5 Sox2+ cells and P2 Sox2+ cells. We identified transcription factors including Pitx2, Id3, Pitx1, Tbx1, Trp63, Nkx2-3, Grhl3, Dlx2, Runx1, Nfix, Zfp536 , etc potentially formed the core transcriptional regulatory networks of Sox2+ DESC in both embryonic and postnatal stages.

AP-2α/AP-2β transcription factors are key regulators of epidermal homeostasis

AP-2 transcription factors regulate ectodermal development but their roles for epidermal homeostasis in the adult skin are unknown. We find that AP-2α is the predominant AP-2 family member in adult epidermis, followed by AP-2β. Through inactivation of AP-2α, AP-2β, or both in keratinocytes we assessed the effects of a gradient of epidermal AP-2 activity on skin function. We find that (1) loss of AP-2β in keratinocytes is compensated for by AP-2α, (2) loss of AP-2α impairs terminal keratinocyte differentiation and hair morphogenesis, and (3) the combined loss of AP-2α/AP-2β results in more severe skin and hair abnormalities. Keratinocyte differentiation defects precede a progressive neutrophilic skin inflammation. Inducible inactivation of AP-2α/AP-2β in the adult phenocopies these manifestations. Transcriptomic analyses of epidermis lacking AP-2α or AP-2α/AP-2β in keratinocytes demonstrate a terminal keratinocyte differentiation defect with upregulation of alarmin keratins and of several immune pathway regulators. Moreover, our analyses suggest a key role of loss of AP-2α-dependent gene expression of CXCL14 and KRT15 as an early pathogenic event towards the manifestation of skin inflammation. Thus, AP-2α/AP-2β are critical regulators of epidermal homeostasis in the adult skin.

Shaping faces: genetic and epigenetic control of craniofacial morphogenesis

Major differences in facial morphology distinguish vertebrate species. Variation of facial traits underlies the uniqueness of human individuals, and abnormal craniofacial morphogenesis during development leads to birth defects that significantly affect quality of life. Studies during the past 40 years have advanced our understanding of the molecular mechanisms that establish facial form during development, highlighting the crucial roles in this process of a multipotent cell type known as the cranial neural crest cell. In this Review, we discuss recent advances in multi-omics and single-cell technologies that enable genes, transcriptional regulatory networks and epigenetic landscapes to be closely linked to the establishment of facial patterning and its variation, with an emphasis on normal and abnormal craniofacial morphogenesis. Advancing our knowledge of these processes will support important developments in tissue engineering, as well as the repair and reconstruction of the abnormal craniofacial complex.

Structural basis for specific DNA sequence motif recognition by the TFAP2 transcription factors

The TFAP2 family regulates gene expression during differentiation, development, and organogenesis, and includes five homologs in humans. They all possess a highly conserved DNA binding domain (DBD) followed by a helix-span-helix (HSH) domain. The DBD-HSH tandem domain specifically binds to a GCC(N3)GGC consensus sequence, but the precise recognition mechanisms remain unclear. Here, we found that TFAP2 preferred binding to the GCC(N3)GGC sequence, and the pseudo-palindromic GCC and GGC motifs and the length of the central spacer between the two motifs determined their binding specificity. Structural studies revealed that the two flat amphipathic α-helical HSH domains of TFAP2A stacked with each other to form a dimer via hydrophobic interactions, while the stabilized loops from both DBD domains inserted into two neighboring major grooves of the DNA duplex to form base-specific interactions. This specific DNA binding mechanism controlled the length of the central spacer and determined the DNA sequence specificity of TFAP2. Mutations of the TFAP2 proteins are implicated in various diseases. We illustrated that reduction or disruption of the DNA binding ability of the TFAP2 proteins is the primary cause of TFAP2 mutation-associated diseases. Thus, our findings also offer valuable insights into the pathogenesis of disease-associated mutations in TFAP2 proteins.

Ap-2β regulates cranial osteogenic potential via the activation of Wnt/β-catenin signaling pathway

The skull is a fundamental bone that protects the development of brain and consists of several bony elements, such as the frontal and parietal bones. Frontal bone exhibited superior in osteogenic potential and regeneration of cranial defects compared to parietal bone. However, how this regional difference is regulated remains largely unknown. In this study, we identified an Ap-2β transcriptional factor with a higher expression in frontal bone, but its molecular function in osteoblasts needs to be elucidated. We found that Ap-2β knockdown in preosteoblasts leads to reduced proliferation, increased cell death and impaired differentiation. Through RNA-seq analysis, we found that Ap-2β influences multiple signaling pathways including the Wnt pathway, and overexpression of Ap-2β showed increased nuclear β-catenin and its target genes expressions in osteoblasts. Pharmacological activation of Wnt/β-catenin signaling using LiCl treatment cannot rescue the reduced luciferase activities of the β-catenin/TCF/LEF reporter in Ap-2β knockdown preosteoblasts. Besides, transient expression of Ap-2β via the lentivirus system could sufficiently rescue the inferior osteogenic potential in parietal osteoblasts, while Ap-2β knockdown in frontal osteoblasts resulted in reduced osteoblast activity, reduced active β-catenin and target genes expressions. Taken together, our data demonstrated that Ap-2β modulates osteoblast proliferation and differentiation through the regulation of Wnt/β-catenin signaling pathway and plays an important role in regulating regional osteogenic potential in frontal and parietal bone.

TFAP2 paralogs regulate midfacial development in part through a conserved ALX genetic pathway

Cranial neural crest development is governed by positional gene regulatory networks (GRNs). Fine-tuning of the GRN components underly facial shape variation, yet how those in the midface are connected and activated remain poorly understood. Here, we show that concerted inactivation of Tfap2a and Tfap2b in the murine neural crest even during the late migratory phase results in a midfacial cleft and skeletal abnormalities. Bulk and single-cell RNA-seq profiling reveal that loss of both Tfap2 members dysregulated numerous midface GRN components involved in midface fusion, patterning, and differentiation. Notably, Alx1/3/4 (Alx) transcript levels are reduced, while ChIP-seq analyses suggest TFAP2 directly and positively regulates Alx gene expression. TFAP2 and ALX co-expression in midfacial neural crest cells of both mouse and zebrafish further implies conservation of this regulatory axis across vertebrates. Consistent with this notion, tfap2a mutant zebrafish present abnormal alx3 expression patterns, and the two genes display a genetic interaction in this species. Together, these data demonstrate a critical role for TFAP2 in regulating vertebrate midfacial development in part through ALX transcription factor gene expression.

Fgf8 regulates first pharyngeal arch segmentation through pouch-cleft interactions

Introduction: The pharyngeal arches are transient developmental structures that, in vertebrates, give rise to tissues of the head and neck. A critical process underlying the specification of distinct arch derivatives is segmentation of the arches along the anterior-posterior axis. Formation of ectodermal-endodermal interfaces is a key mediator of this process, and although it is essential, mechanisms regulating the establishment of these interfaces vary between pouches and between taxa. Methods: Here, we focus on the patterning and morphogenesis of epithelia associated with the first pharyngeal arch, the first pharyngeal pouch (pp1) and the first pharyngeal cleft (pc1), and the role of Fgf8 dosage in these processes in the mouse model system. Results: We find that severe reductions of Fgf8 levels disrupt both pp1 and pc1 development. Notably, out-pocketing of pp1 is largely robust to Fgf8 reductions, however, pp1 extension along the proximal-distal axis fails when Fgf8 is low. Our data indicate that Fgf8 is required for specification of regional identity in both pp1 and pc1, for localized changes in cell polarity, and for elongation and extension of both pp1 and pc1. Discussion: Based on Fgf8-mediated changes in tissue relationships between pp1 and pc1, we hypothesize that extension of pp1 requires physical interaction with pc1. Overall, our data indicate a critical role for the lateral surface ectoderm in segmentation of the first pharyngeal arch that has previously been under-appreciated.

GRHL2 and AP2a coordinate early surface ectoderm lineage commitment during development

Ectodermal dysplasias including skin abnormalities and cleft lip/palate result from improper surface ectoderm (SE) patterning. However, the connection between SE gene regulatory networks and disease remains poorly understood. Here, we dissect human SE differentiation with multiomics and establish GRHL2 as a key mediator of early SE commitment, which acts by skewing cell fate away from the neural lineage. GRHL2 and master SE regulator AP2a balance early cell fate output, with GRHL2 facilitating AP2a binding to SE loci. In turn, AP2a restricts GRHL2 DNA binding away from de novo chromatin contacts. Integration of these regulatory sites with ectodermal dysplasia-associated genomic variants annotated within the Biomedical Data Commons identifies 55 loci previously implicated in craniofacial disorders. These include ABCA4/ARHGAP29 and NOG regulatory regions where disease-linked variants directly affect GRHL2/AP2a binding and gene transcription. These studies elucidate the logic underlying SE commitment and deepen our understanding of human oligogenic disease pathogenesis.

Fgf8 regulates first pharyngeal arch segmentation through pouch-cleft interactions

The pharyngeal arches are transient developmental structures that, in vertebrates, give rise to tissues of the head and neck. A critical process underlying the specification of distinct arch derivatives is segmentation of the arches along the anterior-posterior axis. Out-pocketing of the pharyngeal endoderm between the arches is a key mediator of this process, and although it is essential, mechanisms regulating out-pocketing vary between pouches and between taxa. Here, we focus on the patterning and morphogenesis of epithelia associated with the first pharyngeal arch, the first pharyngeal pouch (pp1) and the first pharyngeal cleft (pc1), and the role of Fgf8 dosage in these processes. We find that severe reductions of Fgf8 levels disrupt both pp1 and pc1 development. Notably, out-pocketing of pp1 is largely robust to Fgf8 reductions, however, pp1 extension along the proximal-distal axis fails when Fgf8 is low. Our data indicate that extension of pp1 requires physical interaction with pc1, and that multiple aspects of pc1 morphogenesis require Fgf8 . In particular, Fgf8 is required for specification of regional identity in both pp1 and pc1, for localized changes in cell polarity, and for elongation and extension of both pp1 and pc1. Overall, our data indicate a critical role for the lateral surface ectoderm in segmentation of the first pharyngeal arch that has previously been under-appreciated.

Understanding the development of oral epithelial organs through single cell transcriptomic analysis

During craniofacial development, the oral epithelium begins as a morphologically homogeneous tissue that gives rise to locally complex structures, including the teeth, salivary glands and taste buds. How the epithelium is initially patterned and specified to generate diverse cell types remains largely unknown. To elucidate the genetic programs that direct the formation of distinct oral epithelial populations, we mapped the transcriptional landscape of embryonic day 12 mouse mandibular epithelia at single cell resolution. Our analysis identified key transcription factors and gene regulatory networks that define different epithelial cell types. By examining the spatiotemporal patterning process along the oral-aboral axis, our results propose a model in which the dental field is progressively confined to its position by the formation of the aboral epithelium anteriorly and the non-dental oral epithelium posteriorly. Using our data, we also identified Ntrk2 as a proliferation driver in the forming incisor, contributing to its invagination. Together, our results provide a detailed transcriptional atlas of the embryonic mandibular epithelium, and unveil new genetic markers and regulators that are present during the specification of various oral epithelial structures.

Cis-regulatory chromatin loops analysis identifies GRHL3 as a master regulator of surface epithelium commitment

Understanding the regulatory network of cell fate acquisition remains a major challenge. Using the induction of surface epithelium (SE) from human embryonic stem cells as a paradigm, we show that the dynamic changes in morphology-related genes (MRGs) closely correspond to SE fate transitions. The marked remodeling of cytoskeleton indicates the initiation of SE differentiation. By integrating promoter interactions, epigenomic features, and transcriptome, we delineate an SE-specific cis-regulatory network and identify grainyhead-like 3 (GRHL3) as an initiation factor sufficient to drive SE commitment. Mechanically, GRHL3 primes the SE chromatin accessibility landscape and activates SE-initiating gene expression. In addition, the evaluation of GRHL3-mediated promoter interactions unveils a positive feedback loop of GRHL3 and bone morphogenetic protein 4 on SE fate decisions. Our work proposes a concept that MRGs could be used to identify cell fate transitions and provides insights into regulatory principles of SE lineage development and stem cell-based regenerative medicine.