GATA3 is essential for separating patterning domains during facial morphogenesis

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.

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.

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.

FOXI3 pathogenic variants cause one form of craniofacial microsomia

Craniofacial microsomia (CFM; also known as Goldenhar syndrome), is a craniofacial developmental disorder of variable expressivity and severity with a recognizable set of abnormalities. These birth defects are associated with structures derived from the first and second pharyngeal arches, can occur unilaterally and include ear dysplasia, microtia, preauricular tags and pits, facial asymmetry and other malformations. The inheritance pattern is controversial, and the molecular etiology of this syndrome is largely unknown. A total of 670 patients belonging to unrelated pedigrees with European and Chinese ancestry with CFM, are investigated. We identify 18 likely pathogenic variants in 21 probands (3.1%) in FOXI3. Biochemical experiments on transcriptional activity and subcellular localization of the likely pathogenic FOXI3 variants, and knock-in mouse studies strongly support the involvement of FOXI3 in CFM. Our findings indicate autosomal dominant inheritance with reduced penetrance, and/or autosomal recessive inheritance. The phenotypic expression of the FOXI3 variants is variable. The penetrance of the likely pathogenic variants in the seemingly dominant form is reduced, since a considerable number of such variants in affected individuals were inherited from non-affected parents. Here we provide suggestive evidence that common variation in the FOXI3 allele in trans with the pathogenic variant could modify the phenotypic severity and accounts for the incomplete penetrance.

Exploration of Novel Genetic Evidence and Clinical Significance Into Hemifacial Microsomia Pathogenesis

The authors browsed through past genetic findings in hemifacial microsomia along with our previously identified mutations in ITGB4 and PDE4DIP from whole genome sequencing of hemifacial microsomia patients. Wondering whether these genes influence mandibular bone modeling by regulation on osteogenesis, the authors approached mechanisms of hemifacial microsomia through this investigation into gene knockdown effects in vitro. MC3T3E1 cells were divided into 5 groups: the negative control group without osteogenesis induction or siRNA, the positive control group with only osteogenesis induction, and 3 gene silenced groups with both osteogenesis induction and siRNA. Validation of transfection was through fluorescence microscopy and quantitative real-time Polymerase chain reaction on knockdown efficiency. Changes in expression levels of the 3 genes during osteogenesis and impact of Itgb4 and Pde4dip knockdown on osteogenesis were examined by quantitative real-time Polymerase chain reaction, alkaline phosphatase, and alizarin red staining. Elevation of osteogenic genes Alpl, Col1a1, Bglap, Spp1, and Runx2 verified successful osteogenesis. Both genes were upregulated under osteogenic induction, while they had different trends over time. Intracellular fluorophores under microscope validated successful transfection and si-m-Itgb4_003, si-m-Pde4dip_002 had satisfactory knockdown effects. During osteogenesis, Pde4dip knockdown enhanced Spp1 expression (1.95±0.13 folds, P =0.045). The authors speculated that these genes may have different involvements in osteogenesis. Stimulated expression of Spp1 by Pde4dip knockdown may suggest that Pde4dip inhibits osteogenesis.

Predicting Modifiers of Genotype-Phenotype Correlations in Craniofacial Development

Most human birth defects are phenotypically variable even when they share a common genetic basis. Our understanding of the mechanisms of this variation is limited, but they are thought to be due to complex gene-environment interactions. Loss of the transcription factor Gata3 associates with the highly variable human birth defects HDR syndrome and microsomia, and can lead to disruption of the neural crest-derived facial skeleton. We have demonstrated that zebrafish gata3 mutants model the variability seen in humans, with genetic background and candidate pathways modifying the resulting phenotype. In this study, we sought to use an unbiased bioinformatic approach to identify environmental modifiers of gata3 mutant craniofacial phenotypes. The LINCs L1000 dataset identifies chemicals that generate differential gene expression that either positively or negatively correlates with an input gene list. These chemicals are predicted to worsen or lessen the mutant phenotype, respectively. We performed RNA-seq on neural crest cells isolated from zebrafish across control, Gata3 loss-of-function, and Gata3 rescue groups. Differential expression analyses revealed 551 potential targets of gata3. We queried the LINCs database with the 100 most upregulated and 100 most downregulated genes. We tested the top eight available chemicals predicted to worsen the mutant phenotype and the top eight predicted to lessen the phenotype. Of these, we found that vinblastine, a microtubule inhibitor, and clofibric acid, a PPAR-alpha agonist, did indeed worsen the gata3 phenotype. The Topoisomerase II and RNA-pol II inhibitors daunorubicin and triptolide, respectively, lessened the phenotype. GO analysis identified Wnt signaling and RNA polymerase function as being enriched in our RNA-seq data, consistent with the mechanism of action of some of the chemicals. Our study illustrates multiple potential pathways for Gata3 function, and demonstrates a systematic, unbiased process to identify modifiers of genotype-phenotype correlations.

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.

Gibbin mesodermal regulation patterns epithelial development

Proper ectodermal patterning during human development requires previously identified transcription factors such as GATA3 and p63, as well as positional signalling from regional mesoderm1-6. However, the mechanism by which ectoderm and mesoderm factors act to stably pattern gene expression and lineage commitment remains unclear. Here we identify the protein Gibbin, encoded by the Xia-Gibbs AT-hook DNA-binding-motif-containing 1 (AHDC1) disease gene7-9, as a key regulator of early epithelial morphogenesis. We find that enhancer- or promoter-bound Gibbin interacts with dozens of sequence-specific zinc-finger transcription factors and methyl-CpG-binding proteins to regulate the expression of mesoderm genes. The loss of Gibbin causes an increase in DNA methylation at GATA3-dependent mesodermal genes, resulting in a loss of signalling between developing dermal and epidermal cell types. Notably, Gibbin-mutant human embryonic stem-cell-derived skin organoids lack dermal maturation, resulting in p63-expressing basal cells that possess defective keratinocyte stratification. In vivo chimeric CRISPR mouse mutants reveal a spectrum of Gibbin-dependent developmental patterning defects affecting craniofacial structure, abdominal wall closure and epidermal stratification that mirror patient phenotypes. Our results indicate that the patterning phenotypes seen in Xia-Gibbs and related syndromes derive from abnormal mesoderm maturation as a result of gene-specific DNA methylation decisions.

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.