A polymorphism in IRF4 affects human pigmentation through a tyrosinase-dependent MITF/TFAP2A pathway.

Sequence polymorphisms linked to human diseases and phenotypes in genome-wide association studies often affect noncoding regions. A SNP within an intron of the gene encoding Interferon Regulatory Factor 4 (IRF4), a transcription factor with no known role in melanocyte biology, is strongly associated with sensitivity of skin to sun exposure, freckles, blue eyes, and brown hair color. Here, we demonstrate that this SNP lies within an enhancer of IRF4 transcription in melanocytes. The allele associated with this pigmentation phenotype impairs binding of the TFAP2A transcription factor that, together with the melanocyte master regulator MITF, regulates activity of the enhancer. Assays in zebrafish and mice reveal that IRF4 cooperates with MITF to activate expression of Tyrosinase (TYR), an essential enzyme in melanin synthesis. Our findings provide a clear example of a noncoding polymorphism that affects a phenotype by modulating a developmental gene regulatory network.

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.

Shared genetic risk between major orofacial cleft phenotypes in an African population.

Nonsyndromic orofacial clefts (NSOFCs) represent a large proportion (70%-80%) of all OFCs. They can be broadly categorized into nonsyndromic cleft lip with or without cleft palate (NSCL/P) and nonsyndromic cleft palate only (NSCPO). Although NSCL/P and NSCPO are considered etiologically distinct, recent evidence suggests the presence of shared genetic risks. Thus, we investigated the genetic overlap between NSCL/P and NSCPO using African genome-wide association study (GWAS) data on NSOFCs. These data consist of 814 NSCL/P, 205 NSCPO cases, and 2159 unrelated controls. We generated common single-nucleotide variants (SNVs) association summary statistics separately for each phenotype (NSCL/P and NSCPO) under an additive genetic model. Subsequently, we employed the pleiotropic analysis under the composite null (PLACO) method to test for genetic overlap. Our analysis identified two loci with genome-wide significance (rs181737795 [p = 2.58E-08] and rs2221169 [p = 4.5E-08]) and one locus with marginal significance (rs187523265 [p = 5.22E-08]). Using mouse transcriptomics data and information from genetic phenotype databases, we identified MDN1, MAP3k7, KMT2A, ARCN1, and VADC2 as top candidate genes for the associated SNVs. These findings enhance our understanding of genetic variants associated with NSOFCs and identify potential candidate genes for further exploration.

TFAP2 paralogs facilitate chromatin access for MITF at pigmentation and cell proliferation genes.

In developing melanocytes and in melanoma cells, multiple paralogs of the Activating-enhancer-binding Protein 2 family of transcription factors (TFAP2) contribute to expression of genes encoding pigmentation regulators, but their interaction with Microphthalmia transcription factor (MITF), a master regulator of these cells, is unclear. Supporting the model that TFAP2 facilitates MITF's ability to activate expression of pigmentation genes, single-cell seq analysis of zebrafish embryos revealed that pigmentation genes are only expressed in the subset of mitfa-expressing cells that also express tfap2 paralogs. To test this model in SK-MEL-28 melanoma cells we deleted the two TFAP2 paralogs with highest expression, TFAP2A and TFAP2C, creating TFAP2 knockout (TFAP2-KO) cells. We then assessed gene expression, chromatin accessibility, binding of TFAP2A and of MITF, and the chromatin marks H3K27Ac and H3K27Me3 which are characteristic of active enhancers and silenced chromatin, respectively. Integrated analyses of these datasets indicate TFAP2 paralogs directly activate enhancers near genes enriched for roles in pigmentation and proliferation, and directly repress enhancers near genes enriched for roles in cell adhesion. Consistently, compared to WT cells, TFAP2-KO cells proliferate less and adhere to one another more. TFAP2 paralogs and MITF co-operatively activate a subset of enhancers, with the former necessary for MITF binding and chromatin accessibility. By contrast, TFAP2 paralogs and MITF do not appear to co-operatively inhibit enhancers. These studies reveal a mechanism by which TFAP2 profoundly influences the set of genes activated by MITF, and thereby the phenotype of pigment cells and melanoma cells.

Transcriptome analyses of murine right and left maxilla-mandibular complex.

OBJECTIVE: The objective of the study was to investigate differential gene expression between murine right and left maxilla-mandibular (MxMn) complexes. SETTING AND SAMPLE POPULATION: Wild-type (WT) C57BL/6 embryonic (E) day 14.5 (n = 3) and 18.5 (n = 3) murine embryos. METHODS: The E14.5 and 18.5 embryos were harvested and hemi-sectioned the MxMn complexes into right and left halves in the mid-sagittal plane. We isolated total RNA using Trizol reagent and further purified using the RNA-easy kit (QIAGEN). We confirmed equal expression of house-keeping genes in right and left halves using RT-PCR and then performed paired-end whole mRNA sequencing in LC Sciences (Houston, TX) followed by differential transcript analyses (>1 or <-1 log fold change; p < .05; q < .05; and FPKM >0.5 in 2/3 samples). The Mouse Genome Informatics and Online Mendelian Inheritance in Man databases as well as gnomAD constraint scores were used to prioritize differentially expressed transcripts. RESULTS: There were 19 upregulated and 19 downregulated transcripts at E14.5 and 8 upregulated and 17 downregulated transcripts at E18.5 time-points. These differentially expressed transcripts were statistically significant and shown to be associated with craniofacial phenotypes in mouse models. These transcripts also have significant gnomAD constraint scores and are enriched in biological processes critical for embryogenesis. CONCLUSIONS: We identified significant differential expression of transcripts between E14.5 and 18.5 murine right and left MxMn complexes. These findings when extrapolated to humans, they may provide a biological basis for facial asymmetry. Further experiments are required to validate these findings in murine models with craniofacial asymmetry.

Anomalous incisor morphology indicates tissue-specific roles for Tfap2a and Tfap2b in tooth development.

Mice possess two types of teeth that differ in their cusp patterns; incisors have one cusp and molars have multiple cusps. The patterning of these two types of teeth relies on fine-tuning of the reciprocal molecular signaling between dental epithelial and mesenchymal tissues during embryonic development. The AP-2 transcription factors, particularly Tfap2a and Tfap2b, are essential components of such epithelial-mesenchymal signaling interactions that coordinate craniofacial development in mice and other vertebrates, but little is known about their roles in the regulation of tooth development and shape. Here we demonstrate that incisors and molars differ in their temporal and spatial expression of Tfap2a and Tfap2b. At the bud stage, Tfap2a is expressed in both the epithelium and mesenchyme of the incisors and molars, but Tfap2b expression is restricted to the molar mesenchyme, only later appearing in the incisor epithelium. Tissue-specific deletions show that loss of the epithelial domain of Tfap2a and Tfap2b affects the number and spatial arrangement of the incisors, notably resulting in duplicated lower incisors. In contrast, deletion of these two genes in the mesenchymal domain has little effect on tooth development. Collectively these results implicate epithelial expression of Tfap2a and Tfap2b in regulating the extent of the dental lamina associated with patterning the incisors and suggest that these genes contribute to morphological differences between anterior (incisor) and posterior (molar) teeth within the mammalian dentition.

AP-2α and AP-2ß cooperatively orchestrate homeobox gene expression during branchial arch patterning.

The evolution of a hinged moveable jaw with variable morphology is considered a major factor behind the successful expansion of the vertebrates. DLX homeobox transcription factors are crucial for establishing the positional code that patterns the mandible, maxilla and intervening hinge domain, but how the genes encoding these proteins are regulated remains unclear. Herein, we demonstrate that the concerted action of the AP-2α and AP-2ß transcription factors within the mouse neural crest is essential for jaw patterning. In the absence of these two proteins, the hinge domain is lost and there are alterations in the size and patterning of the jaws correlating with dysregulation of homeobox gene expression, with reduced levels of Emx, Msx and Dlx paralogs accompanied by an expansion of Six1 expression. Moreover, detailed analysis of morphological features and gene expression changes indicate significant overlap with various compound Dlx gene mutants. Together, these findings reveal that the AP-2 genes have a major function in mammalian neural crest development, influencing patterning of the craniofacial skeleton via the DLX code, an effect that has implications for vertebrate facial evolution, as well as for human craniofacial disorders.

TFAP2 paralogs regulate melanocyte differentiation in parallel with MITF.

Mutations in the gene encoding transcription factor TFAP2A result in pigmentation anomalies in model organisms and premature hair graying in humans. However, the pleiotropic functions of TFAP2A and its redundantly-acting paralogs have made the precise contribution of TFAP2-type activity to melanocyte differentiation unclear. Defining this contribution may help to explain why TFAP2A expression is reduced in advanced-stage melanoma compared to benign nevi. To identify genes with TFAP2A-dependent expression in melanocytes, we profile zebrafish tissue and mouse melanocytes deficient in Tfap2a, and find that expression of a small subset of genes underlying pigmentation phenotypes is TFAP2A-dependent, including Dct, Mc1r, Mlph, and Pmel. We then conduct TFAP2A ChIP-seq in mouse and human melanocytes and find that a much larger subset of pigmentation genes is associated with active regulatory elements bound by TFAP2A. These elements are also frequently bound by MITF, which is considered the "master regulator" of melanocyte development. For example, the promoter of TRPM1 is bound by both TFAP2A and MITF, and we show that the activity of a minimal TRPM1 promoter is lost upon deletion of the TFAP2A binding sites. However, the expression of Trpm1 is not TFAP2A-dependent, implying that additional TFAP2 paralogs function redundantly to drive melanocyte differentiation, which is consistent with previous results from zebrafish. Paralogs Tfap2a and Tfap2b are both expressed in mouse melanocytes, and we show that mouse embryos with Wnt1-Cre-mediated deletion of Tfap2a and Tfap2b in the neural crest almost completely lack melanocytes but retain neural crest-derived sensory ganglia. These results suggest that TFAP2 paralogs, like MITF, are also necessary for induction of the melanocyte lineage. Finally, we observe a genetic interaction between tfap2a and mitfa in zebrafish, but find that artificially elevating expression of tfap2a does not increase levels of melanin in mitfa hypomorphic or loss-of-function mutants. Collectively, these results show that TFAP2 paralogs, operating alongside lineage-specific transcription factors such as MITF, directly regulate effectors of terminal differentiation in melanocytes. In addition, they suggest that TFAP2A activity, like MITF activity, has the potential to modulate the phenotype of melanoma cells.

The old and new face of craniofacial research: How animal models inform human craniofacial genetic and clinical data.

The craniofacial skeletal structures that comprise the human head develop from multiple tissues that converge to form the bones and cartilage of the face. Because of their complex development and morphogenesis, many human birth defects arise due to disruptions in these cellular populations. Thus, determining how these structures normally develop is vital if we are to gain a deeper understanding of craniofacial birth defects and devise treatment and prevention options. In this review, we will focus on how animal model systems have been used historically and in an ongoing context to enhance our understanding of human craniofacial development. We do this by first highlighting "animal to man" approaches; that is, how animal models are being utilized to understand fundamental mechanisms of craniofacial development. We discuss emerging technologies, including high throughput sequencing and genome editing, and new animal repository resources, and how their application can revolutionize the future of animal models in craniofacial research. Secondly, we highlight "man to animal" approaches, including the current use of animal models to test the function of candidate human disease variants. Specifically, we outline a common workflow deployed after discovery of a potentially disease causing variant based on a select set of recent examples in which human mutations are investigated in vivo using animal models. Collectively, these topics will provide a pipeline for the use of animal models in understanding human craniofacial development and disease for clinical geneticist and basic researchers alike.

MEMO1 drives cranial endochondral ossification and palatogenesis.

The cranial base is a component of the neurocranium and has a central role in the structural integration of the face, brain and vertebral column. Consequently, alteration in the shape of the human cranial base has been intimately linked with primate evolution and defective development is associated with numerous human facial abnormalities. Here we describe a novel recessive mutant mouse strain that presented with a domed head and fully penetrant cleft secondary palate coupled with defects in the formation of the underlying cranial base. Mapping and non-complementation studies revealed a specific mutation in Memo1 - a gene originally associated with cell migration. Expression analysis of Memo1 identified robust expression in the perichondrium and periosteum of the developing cranial base, but only modest expression in the palatal shelves. Fittingly, although the palatal shelves failed to elevate in Memo1 mutants, expression changes were modest within the shelves themselves. In contrast, the cranial base, which forms via endochondral ossification had major reductions in the expression of genes responsible for bone formation, notably matrix metalloproteinases and markers of the osteoblast lineage, mirrored by an increase in markers of cartilage and extracellular matrix development. Concomitant with these changes, mutant cranial bases showed an increased zone of hypertrophic chondrocytes accompanied by a reduction in both vascular invasion and mineralization. Finally, neural crest cell-specific deletion of Memo1 caused a failure of anterior cranial base ossification indicating a cell autonomous role for MEMO1 in the development of these neural crest cell derived structures. However, palate formation was largely normal in these conditional mutants, suggesting a non-autonomous role for MEMO1 in palatal closure. Overall, these findings assign a new function to MEMO1 in driving endochondral ossification in the cranium, and also link abnormal development of the cranial base with more widespread effects on craniofacial shape relevant to human craniofacial dysmorphology.

Novel Tfap2-mediated control of soxE expression facilitated the evolutionary emergence of the neural crest.

Gene duplication has been proposed to drive the evolution of novel morphologies. After gene duplication, it is unclear whether changes in the resulting paralogs' coding-regions, or in their cis-regulatory elements, contribute most significantly to the assembly of novel gene regulatory networks. The Transcription Factor Activator Protein 2 (Tfap2) was duplicated in the chordate lineage and is essential for development of the neural crest, a tissue that emerged with vertebrates. Using a tfap2-depleted zebrafish background, we test the ability of available gnathostome, agnathan, cephalochordate and insect tfap2 paralogs to drive neural crest development. With the exception of tfap2d (lamprey and zebrafish), all are able to do so. Together with expression analyses, these results indicate that sub-functionalization has occurred among Tfap2 paralogs, but that neo-functionalization of the Tfap2 protein did not drive the emergence of the neural crest. We investigate whether acquisition of novel target genes for Tfap2 might have done so. We show that in neural crest cells Tfap2 directly activates expression of sox10, which encodes a transcription factor essential for neural crest development. The appearance of this regulatory interaction is likely to have coincided with that of the neural crest, because AP2 and SoxE are not co-expressed in amphioxus, and because neural crest enhancers are not detected proximal to amphioxus soxE. We find that sox10 has limited ability to restore the neural crest in Tfap2-deficient embryos. Together, these results show that mutations resulting in novel Tfap2-mediated regulation of sox10 and other targets contributed to the evolution of the neural crest.

Rare Variants Analyses Suggest Novel Cleft Genes in the African Population.

Non-syndromic orofacial clefts (NSOFCs) are common birth defects with a complex etiology. While over 60 common risk loci have been identified, they explain only a small proportion of the heritability for NSOFC. Rare variants have been implicated in the missing heritability. Thus, our study aimed to identify genes enriched with nonsynonymous rare coding variants associated with NSOFCs. Our sample included 814 non-syndromic cleft lip with or without palate (NSCL/P), 205 non-syndromic cleft palate only (NSCPO), and 2150 unrelated control children from Nigeria, Ghana, and Ethiopia. We conducted a gene-based analysis separately for each phenotype using three rare-variants collapsing models: (1) protein-altering (PA), (2) missense variants only (MO); and (3) loss of function variants only (LOFO). Subsequently, we utilized relevant transcriptomics data to evaluate associated gene expression and examined their mutation constraint using the gnomeAD database. In total, 13 genes showed suggestive associations (p = E-04). Among them, eight genes (ABCB1, ALKBH8, CENPF, CSAD, EXPH5, PDZD8, SLC16A9, and TTC28) were consistently expressed in relevant mouse and human craniofacial tissues during the formation of the face, and three genes (ABCB1, TTC28, and PDZD8) showed statistically significant mutation constraint. These findings underscore the role of rare variants in identifying candidate genes for NSOFCs.

Whole genome sequencing of a family with autosomal dominant features within the oculoauriculovertebral spectrum.

BACKGROUND: Oculoauriculovertebral Spectrum (OAVS) encompasses abnormalities on derivatives from the first and second pharyngeal arches including macrostomia, hemifacial microsomia, micrognathia, preauricular tags, ocular and vertebral anomalies. We present genetic findings on a three-generation family affected with macrostomia, preauricular tags and uni- or bilateral ptosis following an autosomal dominant pattern. METHODS: We generated whole genome sequencing data for the proband, affected parent and unaffected paternal grandparent followed by Sanger sequencing on 23 family members for the top 10 candidate genes: KCND2, PDGFRA, CASP9, NCOA3, WNT10A, SIX1, MTF1, KDR/VEGFR2, LRRK1, and TRIM2 We performed parent and sibling-based transmission disequilibrium tests and burden analysis via a penalized linear mixed model, for segregation and mutation burden respectively. Next, via bioinformatic tools we predicted protein function, mutation pathogenicity and pathway enrichment to investigate the biological relevance of mutations identified. RESULTS: Rare missense mutations in SIX1, KDR/VEGFR2, and PDGFRA showed the best segregation with the OAV phenotypes in this family. When considering any of the 3 OAVS phenotypes as an outcome, SIX1 had the strongest associations in parent-TDTs and sib-TDTs (p=0.025, p=0.052) (unadjusted p-values). Burden analysis identified SIX1 (RC=0.87) and PDGFRA (RC=0.98) strongly associated with OAVS severity. Using phenotype-specific outcomes, sib-TDTs identified SIX1 with uni- or bilateral ptosis (p=0.049) and ear tags (p=0.01), and PDGFRA and KDR/VEGFR2 with ear tags (both p<0.01). CONCLUSION: SIX1, PDGFRA, and KDR/VEGFR2 are strongly associated to OAVS phenotypes. SIX1 has been previously associated with OAVS ear malformations and is co-expressed with EYA1 during ear development. Efforts to strengthen the genotype-phenotype co-relation underlying the OAVS are key to discover etiology, family counseling and prevention.

Rare variants analyses suggest novel cleft genes in the African population.

Non-syndromic orofacial clefts (NSOFCs) are common birth defects with a complex etiology. While over 60 common risk loci have been identified, they explain only a small proportion of the heritability for NSOFCs. Rare variants have been implicated in the missing heritability. Thus, our study aimed to identify genes enriched with nonsynonymous rare coding variants associated with NSOFCs. Our sample included 814 non-syndromic cleft lip with or without palate (NSCL/P), 205 non-syndromic cleft palate only (NSCPO), and 2150 unrelated control children from Nigeria, Ghana, and Ethiopia. We conducted a gene-based analysis separately for each phenotype using three rare-variants collapsing models: (1) protein-altering (PA), (2) missense variants only (MO); and (3) loss of function variants only (LOFO). Subsequently, we utilized relevant transcriptomics data to evaluate associated gene expression and examined their mutation constraint using the gnomeAD database. In total, 13 genes showed suggestive associations (p = E-04). Among them, eight genes (ABCB1, ALKBH8, CENPF, CSAD, EXPH5, PDZD8, SLC16A9, and TTC28) were consistently expressed in relevant mouse and human craniofacial tissues during the formation of the face, and three genes (ABCB1, TTC28, and PDZD8) showed statistically significant mutation constraint. These findings underscore the role of rare variants in identifying candidate genes for NSOFCs.

Gene regulatory evolution and the origin of macroevolutionary novelties: insights from the neural crest.

The appearance of novel anatomic structures during evolution is driven by changes to the networks of transcription factors, signaling pathways, and downstream effector genes controlling development. The nature of the changes to these developmental gene regulatory networks (GRNs) is poorly understood. A striking test case is the evolution of the GRN controlling development of the neural crest (NC). NC cells emerge from the neural plate border (NPB) and contribute to multiple adult structures. While all chordates have a NPB, only in vertebrates do NPB cells express all the genes constituting the neural crest GRN (NC-GRN). Interestingly, invertebrate chordates express orthologs of NC-GRN components in other tissues, revealing that during vertebrate evolution new regulatory connections emerged between transcription factors primitively expressed in the NPB and genes primitively expressed in other tissues. Such interactions could have evolved by two mechanisms. First, transcription factors primitively expressed in the NPB may have evolved new DNA and/or cofactor binding properties (protein neofunctionalization). Alternately, cis-regulatory elements driving NPB expression may have evolved near genes primitively expressed in other tissues (cis-regulatory neofunctionalization). Here we discuss how gene duplication can, in principle, promote either form of neofunctionalization. We review recent published examples of interspecies gene-swap, or regulatory-element-swap, experiments that test both models. Such experiments have yielded little evidence to support the importance of protein neofunctionalization in the emergence of the NC-GRN, but do support the importance of novel cis-regulatory elements in this process. The NC-GRN is an excellent model for the study of gene regulatory and macroevolutionary innovation.

Differentiation of zebrafish melanophores depends on transcription factors AP2 alpha and AP2 epsilon.

A model of the gene-regulatory-network (GRN), governing growth, survival, and differentiation of melanocytes, has emerged from studies of mouse coat color mutants and melanoma cell lines. In this model, Transcription Factor Activator Protein 2 alpha (TFAP2A) contributes to melanocyte development by activating expression of the gene encoding the receptor tyrosine kinase Kit. Next, ligand-bound Kit stimulates a pathway activating transcription factor Microphthalmia (Mitf), which promotes differentiation and survival of melanocytes by activating expression of Tyrosinase family members, Bcl2, and other genes. The model predicts that in both Tfap2a and Kit null mutants there will be a phenotype of reduced melanocytes and that, because Tfap2a acts upstream of Kit, this phenotype will be more severe, or at least as severe as, in Tfap2a null mutants in comparison to Kit null mutants. Unexpectedly, this is not the case in zebrafish or mouse. Because many Tfap2 family members have identical DNA-binding specificity, we reasoned that another Tfap2 family member may work redundantly with Tfap2a in promoting Kit expression. We report that tfap2e is expressed in melanoblasts and melanophores in zebrafish embryos and that its orthologue, TFAP2E, is expressed in human melanocytes. We provide evidence that Tfap2e functions redundantly with Tfap2a to maintain kita expression in zebrafish embryonic melanophores. Further, we show that, in contrast to in kita mutants where embryonic melanophores appear to differentiate normally, in tfap2a/e doubly-deficient embryonic melanophores are small and under-melanized, although they retain expression of mitfa. Interestingly, forcing expression of mitfa in tfap2a/e doubly-deficient embryos partially restores melanophore differentiation. These findings reveal that Tfap2 activity, mediated redundantly by Tfap2a and Tfap2e, promotes melanophore differentiation in parallel with Mitf by an effector other than Kit. This work illustrates how analysis of single-gene mutants may fail to identify steps in a GRN that are affected by the redundant activity of related proteins.

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.

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.

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

The facial surface ectoderm is essential for normal development of the underlying cranial neural crest cell populations, providing signals that direct appropriate growth, patterning, and morphogenesis. Despite the importance of the ectoderm as a signaling center, the molecular cues and genetic programs implemented within this tissue are understudied. Here, we show that removal of two members of the AP-2 transcription factor family, AP-2α and AP-2ß, within the early embryonic ectoderm of the mouse leads to major alterations in the craniofacial complex. Significantly, there are clefts in both the upper face and mandible, accompanied by fusion of the upper and lower jaws in the hinge region. Comparison of ATAC-seq and RNA-seq analyses between controls and mutants revealed significant changes in chromatin accessibility and gene expression centered on multiple AP-2 binding motifs associated with enhancer elements within these ectodermal lineages. In particular, loss of these AP-2 proteins affects both skin differentiation as well as multiple signaling pathways, most notably the WNT pathway. We also determined that the mutant clefting phenotypes that correlated with reduced WNT signaling could be rescued by Wnt1 ligand overexpression in the ectoderm. Collectively, these findings highlight a conserved ancestral function for AP-2 transcription factors in ectodermal development and signaling, and provide a framework from which to understand the gene regulatory network operating within this tissue that directs vertebrate craniofacial development.

The Skull's Girder: A Brief Review of the Cranial Base.

The cranial base is a multifunctional bony platform within the core of the cranium, spanning rostral to caudal ends. This structure provides support for the brain and skull vault above, serves as a link between the head and the vertebral column below, and seamlessly integrates with the facial skeleton at its rostral end. Unique from the majority of the cranial skeleton, the cranial base develops from a cartilage intermediate-the chondrocranium-through the process of endochondral ossification. Owing to the intimate association of the cranial base with nearly all aspects of the head, congenital birth defects impacting these structures often coincide with anomalies of the cranial base. Despite this critical importance, studies investigating the genetic control of cranial base development and associated disorders lags in comparison to other craniofacial structures. Here, we highlight and review developmental and genetic aspects of the cranial base, including its transition from cartilage to bone, dual embryological origins, and vignettes of transcription factors controlling its formation.