A common cis-regulatory variant impacts normal-range and disease-associated human facial shape through regulation of PKDCC during chondrogenesis.

Genome-wide association studies (GWAS) identified thousands of genetic variants linked to phenotypic traits and disease risk. However, mechanistic understanding of how GWAS variants influence complex morphological traits and can, in certain cases, simultaneously confer normal-range phenotypic variation and disease predisposition, is still largely lacking. Here, we focus on rs6740960, a single nucleotide polymorphism (SNP) at the 2p21 locus, which in GWAS studies has been associated both with normal-range variation in jaw shape and with an increased risk of non-syndromic orofacial clefting. Using in vitro derived embryonic cell types relevant for human facial morphogenesis, we show that this SNP resides in an enhancer that regulates chondrocytic expression of PKDCC - a gene encoding a tyrosine kinase involved in chondrogenesis and skeletal development. In agreement, we demonstrate that the rs6740960 SNP is sufficient to confer chondrocyte-specific differences in PKDCC expression. By deploying dense landmark morphometric analysis of skull elements in mice, we show that changes in Pkdcc dosage are associated with quantitative changes in the maxilla, mandible, and palatine bone shape that are concordant with the facial phenotypes and disease predisposition seen in humans. We further demonstrate that the frequency of the rs6740960 variant strongly deviated among different human populations, and that the activity of its cognate enhancer diverged in hominids. Our study provides a mechanistic explanation of how a common SNP can mediate normal-range and disease-associated morphological variation, with implications for the evolution of human facial features.

Palatal segment contributions to midfacial growth in three inbred mouse strains.

Following facial prominence fusion, anterior-posterior (A-P) elongation of the palate is a critical aspect of palatogenesis and integrated midfacial elongation. Reciprocal epithelial-mesenchymal interactions drive secondary palate elongation and periodic signaling center formation within the rugae growth zone (RGZ). However, the relationship between RGZ dynamics and the morphogenetic behavior of underlying palatal bone mesenchymal precursors has remained enigmatic. Our results indicate that cellular activity at the RGZ simultaneously drives rugae formation and elongation of the maxillary bone primordium within the anterior secondary palate, which more than doubles in length prior to palatal shelf fusion. The first formed palatal ruga, found just posterior to the RGZ, represents a consistent morphological boundary between anterior and posterior secondary palate bone precursors, being found at the future maxillary-palatine suture. These results suggest a model where changes in RGZ-driven A-P growth of the anterior secondary palate may produce interspecies and intraspecies differences in facial prognathism and differences in the proportional contribution of palatal segment-associated bones to total palate length. An ontogenetic comparison of three inbred mouse strains indicated that while RGZ-driven growth of the anterior secondary palate is critical for early midfacial outgrowth, subtle strain-specific bony contributions to adult palate length are not present until after this initial palatal growth period. This multifaceted illustration of normal midfacial growth dynamics confirms a one-to-one relationship between palatal segments and upper jaw bones during the earliest stages of palatal growth, which may serve as the basis for evolutionary change in upper jaw morphology. Additionally, identified mouse strain-specific differences in palatal segment elongation provide a useful foundation for understanding the impact of background genetic effects on facial morphogenesis.

Loss of Extreme Long-Range Enhancers in Human Neural Crest Drives a Craniofacial Disorder.

Non-coding mutations at the far end of a large gene desert surrounding the SOX9 gene result in a human craniofacial disorder called Pierre Robin sequence (PRS). Leveraging a human stem cell differentiation model, we identify two clusters of enhancers within the PRS-associated region that regulate SOX9 expression during a restricted window of facial progenitor development at distances up to 1.45 Mb. Enhancers within the 1.45 Mb cluster exhibit highly synergistic activity that is dependent on the Coordinator motif. Using mouse models, we demonstrate that PRS phenotypic specificity arises from the convergence of two mechanisms: confinement of Sox9 dosage perturbation to developing facial structures through context-specific enhancer activity and heightened sensitivity of the lower jaw to Sox9 expression reduction. Overall, we characterize the longest-range human enhancers involved in congenital malformations, directly demonstrate that PRS is an enhanceropathy, and illustrate how small changes in gene expression can lead to morphological variation.

FaceBase 3: analytical tools and FAIR resources for craniofacial and dental research.

The FaceBase Consortium was established by the National Institute of Dental and Craniofacial Research in 2009 as a 'big data' resource for the craniofacial research community. Over the past decade, researchers have deposited hundreds of annotated and curated datasets on both normal and disordered craniofacial development in FaceBase, all freely available to the research community on the FaceBase Hub website. The Hub has developed numerous visualization and analysis tools designed to promote integration of multidisciplinary data while remaining dedicated to the FAIR principles of data management (findability, accessibility, interoperability and reusability) and providing a faceted search infrastructure for locating desired data efficiently. Summaries of the datasets generated by the FaceBase projects from 2014 to 2019 are provided here. FaceBase 3 now welcomes contributions of data on craniofacial and dental development in humans, model organisms and cell lines. Collectively, the FaceBase Consortium, along with other NIH-supported data resources, provide a continuously growing, dynamic and current resource for the scientific community while improving data reproducibility and fulfilling data sharing requirements.

Pbx loss in cranial neural crest, unlike in epithelium, results in cleft palate only and a broader midface.

Orofacial clefting represents the most common craniofacial birth defect. Cleft lip with or without cleft palate (CL/P) is genetically distinct from cleft palate only (CPO). Numerous transcription factors (TFs) regulate normal development of the midface, comprising the premaxilla, maxilla and palatine bones, through control of basic cellular behaviors. Within the Pbx family of genes encoding Three Amino-acid Loop Extension (TALE) homeodomain-containing TFs, we previously established that in the mouse, Pbx1 plays a preeminent role in midfacial morphogenesis, and Pbx2 and Pbx3 execute collaborative functions in domains of coexpression. We also reported that Pbx1 loss from cephalic epithelial domains, on a Pbx2- or Pbx3-deficient background, results in CL/P via disruption of a regulatory network that controls apoptosis at the seam of frontonasal and maxillary process fusion. Conversely, Pbx1 loss in cranial neural crest cell (CNCC)-derived mesenchyme on a Pbx2-deficient background results in CPO, a phenotype not yet characterized. In this study, we provide in-depth analysis of PBX1 and PBX2 protein localization from early stages of midfacial morphogenesis throughout development of the secondary palate. We further establish CNCC-specific roles of PBX TFs and describe the developmental abnormalities resulting from their loss in the murine embryonic secondary palate. Additionally, we compare and contrast the phenotypes arising from PBX1 loss in CNCC with those caused by its loss in the epithelium and show that CNCC-specific Pbx1 deletion affects only later secondary palate morphogenesis. Moreover, CNCC mutants exhibit perturbed rostro-caudal organization and broadening of the midfacial complex. Proliferation defects are pronounced in CNCC mutants at gestational day (E)12.5, suggesting altered proliferation of mutant palatal progenitor cells, consistent with roles of PBX factors in maintaining progenitor cell state. Although the craniofacial skeletal abnormalities in CNCC mutants do not result from overt patterning defects, osteogenesis is delayed, underscoring a critical role of PBX factors in CNCC morphogenesis and differentiation. Overall, the characterization of tissue-specific Pbx loss-of-function mouse models with orofacial clefting establishes these strains as unique tools to further dissect the complexities of this congenital craniofacial malformation. This study closely links PBX TALE homeodomain proteins to the variation in maxillary shape and size that occurs in pathological settings and during evolution of midfacial morphology.