Expansion of the sagittal suture induces proliferation of skeletal stem cells and sustains endogenous calvarial bone regeneration.

In newborn humans, and up to approximately 2 y of age, calvarial bone defects can naturally regenerate. This remarkable regeneration potential is also found in newborn mice and is absent in adult mice. Since previous studies showed that the mouse calvarial sutures are reservoirs of calvarial skeletal stem cells (cSSCs), which are the cells responsible for calvarial bone regeneration, here we hypothesized that the regenerative potential of the newborn mouse calvaria is due to a significant amount of cSSCs present in the newborn expanding sutures. Thus, we tested whether such regenerative potential can be reverse engineered in adult mice by artificially inducing an increase of the cSSCs resident within the adult calvarial sutures. First, we analyzed the cellular composition of the calvarial sutures in newborn and in older mice, up to 14-mo-old mice, showing that the sutures of the younger mice are enriched in cSSCs. Then, we demonstrated that a controlled mechanical expansion of the functionally closed sagittal sutures of adult mice induces a significant increase of the cSSCs. Finally, we showed that if a calvarial critical size bone defect is created simultaneously to the mechanical expansion of the sagittal suture, it fully regenerates without the need for additional therapeutic aids. Using a genetic blockade system, we further demonstrate that this endogenous regeneration is mediated by the canonical Wnt signaling. This study shows that controlled mechanical forces can harness the cSSCs and induce calvarial bone regeneration. Similar harnessing strategies may be used to develop novel and more effective bone regeneration autotherapies.

Damaging variants in FOXI3 cause microtia and craniofacial microsomia.

PURPOSE: Craniofacial microsomia (CFM) represents a spectrum of craniofacial malformations, ranging from isolated microtia with or without aural atresia to underdevelopment of the mandible, maxilla, orbit, facial soft tissue, and/or facial nerve. The genetic causes of CFM remain largely unknown. METHODS: We performed genome sequencing and linkage analysis in patients and families with microtia and CFM of unknown genetic etiology. The functional consequences of damaging missense variants were evaluated through expression of wild-type and mutant proteins in vitro. RESULTS: We studied a 5-generation kindred with microtia, identifying a missense variant in FOXI3 (p.Arg236Trp) as the cause of disease (logarithm of the odds = 3.33). We subsequently identified 6 individuals from 3 additional kindreds with microtia-CFM spectrum phenotypes harboring damaging variants in FOXI3, a regulator of ectodermal and neural crest development. Missense variants in the nuclear localization sequence were identified in cases with isolated microtia with aural atresia and found to affect subcellular localization of FOXI3. Loss of function variants were found in patients with microtia and mandibular hypoplasia (CFM), suggesting dosage sensitivity of FOXI3. CONCLUSION: Damaging variants in FOXI3 are the second most frequent genetic cause of CFM, causing 1% of all cases, including 13% of familial cases in our cohort.

Isolated frontosphenoidal craniosynostosis: An argument for genetic testing.

Isolated frontosphenoidal craniosynostosis (IFSC) is a rare congenital defect defined as premature fusion of the frontosphenoidal suture in the absence of other suture fusion. Until now, IFSC was regarded as a phenomenon with an unclear genetic etiology. We have identified three cases with IFSC with underlying syndromic diagnoses that were attributable to pathogenic mutations involving FGFR3 and MN1, as well as 22q11.2 deletion syndrome. These findings suggest a genetic predisposition to IFSC may exist, thereby justifying the recommendation for genetic evaluation and testing in this population. Furthermore, due to improved imaging resolution, cases of IFSC are now readily identified. With the identification of IFSC with underlying genetic diagnoses, in combination with significant improvements in imaging resolution, we recommend genetic evaluation in children with IFSC.

Unexpected role of SIX1 variants in craniosynostosis: expanding the phenotype of SIX1-related disorders.

BACKGROUND: Pathogenic heterozygous SIX1 variants (predominantly missense) occur in branchio-otic syndrome (BOS), but an association with craniosynostosis has not been reported. METHODS: We investigated probands with craniosynostosis of unknown cause using whole exome/genome (n=628) or RNA (n=386) sequencing, and performed targeted resequencing of SIX1 in 615 additional patients. Expression of SIX1 protein in embryonic cranial sutures was examined in the Six1nLacZ/+ reporter mouse. RESULTS: From 1629 unrelated cases with craniosynostosis we identified seven different SIX1 variants (three missense, including two de novo mutations, and four nonsense, one of which was also present in an affected twin). Compared with population data, enrichment of SIX1 loss-of-function variants was highly significant (p=0.00003). All individuals with craniosynostosis had sagittal suture fusion; additionally four had bilambdoid synostosis. Associated BOS features were often attenuated; some carrier relatives appeared non-penetrant. SIX1 is expressed in a layer basal to the calvaria, likely corresponding to the dura mater, and in the mid-sagittal mesenchyme. CONCLUSION: Craniosynostosis is associated with heterozygous SIX1 variants, with possible enrichment of loss-of-function variants compared with classical BOS. We recommend screening of SIX1 in craniosynostosis, particularly when sagittal±lambdoid synostosis and/or any BOS phenotypes are present. These findings highlight the role of SIX1 in cranial suture homeostasis.

Transcriptomic Signatures of Single-Suture Craniosynostosis Phenotypes.

Craniosynostosis is a birth defect where calvarial sutures close prematurely, as part of a genetic syndrome or independently, with unknown cause. This study aimed to identify differences in gene expression in primary calvarial cell lines derived from patients with four phenotypes of single-suture craniosynostosis, compared to controls. Calvarial bone samples (N = 388 cases/85 controls) were collected from clinical sites during reconstructive skull surgery. Primary cell lines were then derived from the tissue and used for RNA sequencing. Linear models were fit to estimate covariate adjusted associations between gene expression and four phenotypes of single-suture craniosynostosis (lambdoid, metopic, sagittal, and coronal), compared to controls. Sex-stratified analysis was also performed for each phenotype. Differentially expressed genes (DEGs) included 72 genes associated with coronal, 90 genes associated with sagittal, 103 genes associated with metopic, and 33 genes associated with lambdoid craniosynostosis. The sex-stratified analysis revealed more DEGs in males (98) than females (4). There were 16 DEGs that were homeobox (HOX) genes. Three TFs (SUZ12, EZH2, AR) significantly regulated expression of DEGs in one or more phenotypes. Pathway analysis identified four KEGG pathways associated with at least one phenotype of craniosynostosis. Together, this work suggests unique molecular mechanisms related to craniosynostosis phenotype and fetal sex.

A genome-wide association study implicates the BMP7 locus as a risk factor for nonsyndromic metopic craniosynostosis.

Our previous genome-wide association study (GWAS) for sagittal nonsyndromic craniosynostosis (sNCS) provided important insights into the genetics of midline CS. In this study, we performed a GWAS for a second midline NCS, metopic NCS (mNCS), using 215 non-Hispanic white case-parent triads. We identified six variants with genome-wide significance (P ≤ 5 × 10-8): rs781716 (P = 4.71 × 10-9; odds ratio [OR] = 2.44) intronic to SPRY3; rs6127972 (P = 4.41 × 10-8; OR = 2.17) intronic to BMP7; rs62590971 (P = 6.22 × 10-9; OR = 0.34), located ~ 155 kb upstream from TGIF2LX; and rs2522623, rs2573826, and rs2754857, all intronic to PCDH11X (P = 1.76 × 10-8, OR = 0.45; P = 3.31 × 10-8, OR = 0.45; P = 1.09 × 10-8, OR = 0.44, respectively). We performed a replication study of these variants using an independent non-Hispanic white sample of 194 unrelated mNCS cases and 333 unaffected controls; only the association for rs6127972 (P = 0.004, OR = 1.45; meta-analysis P = 1.27 × 10-8, OR = 1.74) was replicated. Our meta-analysis examining single nucleotide polymorphisms common to both our mNCS and sNCS studies showed the strongest association for rs6127972 (P = 1.16 × 10-6). Our imputation analysis identified a linkage disequilibrium block encompassing rs6127972, which contained an enhancer overlapping a CTCF transcription factor binding site (chr20:55,798,821-55,798,917) that was significantly hypomethylated in mesenchymal stem cells derived from fused metopic compared to open sutures from the same probands. This study provides additional insights into genetic factors in midline CS.

Correction: SMAD6 variants in craniosynostosis: genotype and phenotype evaluation.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

SMAD6 variants in craniosynostosis: genotype and phenotype evaluation.

PURPOSE: Enrichment of heterozygous missense and truncating SMAD6 variants was previously reported in nonsyndromic sagittal and metopic synostosis, and interaction of SMAD6 variants with a common polymorphism nearBMP2 (rs1884302) was proposed to contribute to inconsistent penetrance. We determined the occurrence of SMAD6 variants in all types of craniosynostosis, evaluated the impact of different missense variants on SMAD6 function, and tested independently whether rs1884302 genotype significantly modifies the phenotype. METHODS: We performed resequencing of SMAD6 in 795 unsolved patients with any type of craniosynostosis and genotyped rs1884302 in SMAD6-positive individuals and relatives. We examined the inhibitory activity and stability of SMAD6 missense variants. RESULTS: We found 18 (2.3%) different rare damaging SMAD6 variants, with the highest prevalence in metopic synostosis (5.8%) and an 18.3-fold enrichment of loss-of-function variants comparedwith gnomAD data (P < 10-7). Combined with eight additional variants, ≥20/26 were transmitted from an unaffected parent but rs1884302 genotype did not predict phenotype. CONCLUSION: Pathogenic SMAD6 variants substantially increase the risk of both nonsyndromic and syndromic presentations of craniosynostosis, especially metopic synostosis. Functional analysis is important to evaluate missense variants. Genotyping of rs1884302 is not clinically useful. Mechanisms to explain the remarkable diversity of phenotypes associated with SMAD6 variants remain obscure.

Motor Function in School-Aged Children With Positional Plagiocephaly or Brachycephaly.

OBJECTIVE: To determine whether children with a history of positional plagiocephaly/brachycephaly (PPB) show persistent deficits in motor development. METHODS: In a longitudinal cohort study, we completed follow-up assessments with 187 school-aged children with PPB and 149 participants without PPB who were originally enrolled in infancy. Primary outcomes were the Bruininks-Oseretsky Test of Motor Proficiency-Second Edition (BOT-2) composite scores. RESULTS: Children with PPB scored lower than controls on the BOT-2. Stratified analyses indicated that differences were restricted to children who had moderate-severe PPB. No consistent differences were observed in children who had mild PPB. CONCLUSION: Children who had moderate-severe PPB in infancy show persistent differences in motor function. We suggest close developmental monitoring and early intervention to address motor deficits.

A genotype-specific surgical approach for patients with Pfeiffer syndrome due to W290C pathogenic variant in FGFR2 is associated with improved developmental outcomes and reduced mortality.

PURPOSE: Among children with FGFR2-associated Pfeiffer syndrome, those with the W290C pathogenic variant (PV) are reported to have the worst clinical outcomes. Mortality is high, and severe neurocognitive impairment has been reported in all surviving patients. However, it is unclear whether these poor outcomes are an unavoidable consequence of the PV itself, or could be improved with a genotype-specific treatment approach. The purpose of this report is to describe the more intensive surgical approach used for each of the three patients with W290C PV in FGFR2 at our center, all of whom survived and have normal neurocognitive functioning. METHODS: Retrospective chart review. RESULTS: In contrast to other patients with Pfeiffer syndrome at our center, all three patients who were subsequently found to have a W290C PV required a similar and more aggressive approach based on early cephalocranial disproportion. In contrast to previously reported W290C cases, each of our three patients survived and demonstrate normal neurocognitive functioning. CONCLUSION: While previously reported outcomes in W290C-associated Pfeiffer syndrome have been extremely poor, we present three patients who underwent an intensive surgical approach and have normal development. This suggests that a personalized and aggressive surgical approach for children with W290C PV may dramatically improve clinical outcome.

Genome-Wide Association Study Reveals Multiple Loci Influencing Normal Human Facial Morphology.

Numerous lines of evidence point to a genetic basis for facial morphology in humans, yet little is known about how specific genetic variants relate to the phenotypic expression of many common facial features. We conducted genome-wide association meta-analyses of 20 quantitative facial measurements derived from the 3D surface images of 3118 healthy individuals of European ancestry belonging to two US cohorts. Analyses were performed on just under one million genotyped SNPs (Illumina OmniExpress+Exome v1.2 array) imputed to the 1000 Genomes reference panel (Phase 3). We observed genome-wide significant associations (p < 5 x 10-8) for cranial base width at 14q21.1 and 20q12, intercanthal width at 1p13.3 and Xq13.2, nasal width at 20p11.22, nasal ala length at 14q11.2, and upper facial depth at 11q22.1. Several genes in the associated regions are known to play roles in craniofacial development or in syndromes affecting the face: MAFB, PAX9, MIPOL1, ALX3, HDAC8, and PAX1. We also tested genotype-phenotype associations reported in two previous genome-wide studies and found evidence of replication for nasal ala length and SNPs in CACNA2D3 and PRDM16. These results provide further evidence that common variants in regions harboring genes of known craniofacial function contribute to normal variation in human facial features. Improved understanding of the genes associated with facial morphology in healthy individuals can provide insights into the pathways and mechanisms controlling normal and abnormal facial morphogenesis.

Activation of the IGF1 pathway mediates changes in cellular contractility and motility in single-suture craniosynostosis.

Insulin growth factor 1 (IGF1) is a major anabolic signal that is essential during skeletal development, cellular adhesion and migration. Recent transcriptomic studies have shown that there is an upregulation in IGF1 expression in calvarial osteoblasts derived from patients with single-suture craniosynostosis (SSC). Upregulation of the IGF1 signaling pathway is known to induce increased expression of a set of osteogenic markers that previously have been shown to be correlated with contractility and migration. Although the IGF1 signaling pathway has been implicated in SSC, a correlation between IGF1, contractility and migration has not yet been investigated. Here, we examined the effect of IGF1 activation in inducing cellular contractility and migration in SSC osteoblasts using micropost arrays and time-lapse microscopy. We observed that the contractile forces and migration speeds of SSC osteoblasts correlated with IGF1 expression. Moreover, both contractility and migration of SSC osteoblasts were directly affected by the interaction of IGF1 with IGF1 receptor (IGF1R). Our results suggest that IGF1 activity can provide valuable insight for phenotype-genotype correlation in SSC osteoblasts and might provide a target for therapeutic intervention.

Associations between laterality of orofacial clefts and medical and academic outcomes.

Patients with oral clefts have an increased risk of other malformations, syndromes, and lower academic performance in school. Few studies have investigated if laterality of clefts is associated with medical and academic outcomes. Oral clefts have nonrandom laterality, with left-sided clefts occurring approximately twice as often as right-sided clefts. Using a retrospective study design, we examined potential associations of cleft attributes and outcomes in patients with cleft lip with or without cleft palate (CL/P) born in 2003-2010 who were treated at the Seattle Children's Craniofacial Center. The following variables were extracted from medical records: cleft type, medical history, maternal hyperglycemia, other malformations, and the need for academic support at school. We used logistic regression to examine risk of associations with outcomes of interest. Relative to patients with left-sided clefts, patients with bilateral CL/P were more likely to have a syndrome. Patients with nonsyndromic right-sided CL/P had a higher risk (OR and 95%CI: 3.5, 1.3-9.5, and 5.5, 1.9-16.0, respectively) of having other malformations and requiring academic support at school, when compared to patients with left-sided CL/P. Understanding the etiology of oral clefts is complicated, in part because both genetic and environmental factors contribute to the risk of developing a cleft. However, the different outcomes associated with cleft laterality suggest that right-sided clefts may have a distinct etiology. Using laterality to study cleft subgroups may advance our understanding of the etiology of this common birth defect.

Single suture craniosynostosis: Identification of rare variants in genes associated with syndromic forms.

We report RNA-Sequencing results on a cohort of patients with single suture craniosynostosis and demonstrate significant enrichment of heterozygous, rare, and damaging variants among key craniosynostosis-related genes. Genetic burden analysis identified a significant increase in damaging variants in ATR, EFNA4, ERF, MEGF8, SCARF2, and TGFBR2. Of 391 participants, 15% were found to have damaging and potentially causal variants in 29 genes. We observed transmission in 96% of the affected individuals, and thus penetrance, epigenetics, and oligogenic factors need to be considered when recommending genetic testing in patients with nonsyndromic craniosynostosis.

Exome sequencing identifies a recurrent de novo ZSWIM6 mutation associated with acromelic frontonasal dysostosis.

Acromelic frontonasal dysostosis (AFND) is a rare disorder characterized by distinct craniofacial, brain, and limb malformations, including frontonasal dysplasia, interhemispheric lipoma, agenesis of the corpus callosum, tibial hemimelia, preaxial polydactyly of the feet, and intellectual disability. Exome sequencing of one trio and two unrelated probands revealed the same heterozygous variant (c.3487C>T [p. Arg1163Trp]) in a highly conserved protein domain of ZSWIM6; this variant has not been seen in the 1000 Genomes data, dbSNP, or the Exome Sequencing Project. Sanger validation of the three trios confirmed that the variant was de novo and was also present in a fourth isolated proband. In situ hybridization of early zebrafish embryos at 24 hr postfertilization (hpf) demonstrated telencephalic expression of zswim6 and onset of midbrain, hindbrain, and retinal expression at 48 hpf. Immunohistochemistry of later-stage mouse embryos demonstrated tissue-specific expression in the derivatives of all three germ layers. qRT-PCR expression analysis of osteoblast and fibroblast cell lines available from two probands was suggestive of Hedgehog pathway activation, indicating that the ZSWIM6 mutation associated with AFND may lead to the craniofacial, brain and limb malformations through the disruption of Hedgehog signaling.

The 3D Facial Norms Database: Part 1. A Web-Based Craniofacial Anthropometric and Image Repository for the Clinical and Research Community.

With the current widespread use of three-dimensional (3D) facial surface imaging in clinical and research environments, there is a growing demand for high-quality craniofacial norms based on 3D imaging technology. The principal goal of the 3D Facial Norms (3DFN) project was to create an interactive, Web-based repository of 3D facial images and measurements. Unlike other repositories, users can gain access to both summary-level statistics and individual-level data, including 3D facial landmark coordinates, 3D-derived anthropometric measurements, 3D facial surface images, and genotypes from every individual in the dataset. The 3DFN database currently consists of 2454 male and female participants ranging in age from 3 to 40 years. The subjects were recruited at four US sites and screened for a history of craniofacial conditions. The goal of this article is to introduce readers to the 3DFN repository by providing a general overview of the project, explaining the rationale behind the creation of the database, and describing the methods used to collect the data. Sex- and age-specific summary statistics (means and standard deviations) and growth curves for every anthropometric measurement in the 3DFN dataset are provided as a supplement available online. These summary statistics and growth curves can aid clinicians in the assessment of craniofacial dysmorphology.

A human homeotic transformation resulting from mutations in PLCB4 and GNAI3 causes auriculocondylar syndrome.

Auriculocondylar syndrome (ACS) is a rare, autosomal-dominant craniofacial malformation syndrome characterized by variable micrognathia, temporomandibular joint ankylosis, cleft palate, and a characteristic "question-mark" ear malformation. Careful phenotypic characterization of severely affected probands in our cohort suggested the presence of a mandibular patterning defect resulting in a maxillary phenotype (i.e., homeotic transformation). We used exome sequencing of five probands and identified two novel (exclusive to the patient and/or family studied) missense mutations in PLCB4 and a shared mutation in GNAI3 in two unrelated probands. In confirmatory studies, three additional novel PLCB4 mutations were found in multigenerational ACS pedigrees. All mutations were confirmed by Sanger sequencing, were not present in more than 10,000 control chromosomes, and resulted in amino-acid substitutions located in highly conserved protein domains. Additionally, protein-structure modeling demonstrated that all ACS substitutions disrupt the catalytic sites of PLCB4 and GNAI3. We suggest that PLCB4 and GNAI3 are core signaling molecules of the endothelin-1-distal-less homeobox 5 and 6 (EDN1-DLX5/DLX6) pathway. Functional studies demonstrated a significant reduction in downstream DLX5 and DLX6 expression in ACS cases in assays using cultured osteoblasts from probands and controls. These results support the role of the previously implicated EDN1-DLX5/6 pathway in regulating mandibular specification in other species, which, when disrupted, results in a maxillary phenotype. This work defines the molecular basis of ACS as a homeotic transformation (mandible to maxilla) in humans.

Expanding the genetic and phenotypic spectrum of popliteal pterygium disorders.

The popliteal pterygia syndromes are a distinct subset of the hundreds of Mendelian orofacial clefting syndromes. Popliteal pterygia syndromes have considerable variability in severity and in the associated phenotypic features but are all characterized by cutaneous webbing across one or more major joints, cleft lip and/or palate, syndactyly, and genital malformations. Heterozygous mutations in IRF6 cause popliteal pterygium syndrome (PPS) while homozygous mutations in RIPK4 or CHUK (IKKA) cause the more severe Bartsocas-Papas syndrome (BPS) and Cocoon syndrome, respectively. In this study, we report mutations in six pedigrees with children affected with PPS or BPS. Using a combination of Sanger and exome sequencing, we report the first case of an autosomal recessive popliteal pterygium syndrome caused by homozygous mutation of IRF6 and the first case of uniparental disomy of chromosome 21 leading to a recessive disorder. We also demonstrate that mutations in RIPK4 can cause features with a range of severity along the PPS-BPS spectrum and that mutations in IKKA can cause a range of features along the BPS-Cocoon spectrum. Our findings have clinical implications for genetic counseling of families with pterygia syndromes and further implicate IRF6, RIPK4, and CHUK (IKKA) in potentially interconnected pathways governing epidermal and craniofacial development.

Clinical delineation and natural history of the PIK3CA-related overgrowth spectrum.

Somatic mutations in the phosphatidylinositol/AKT/mTOR pathway cause segmental overgrowth disorders. Diagnostic descriptors associated with PIK3CA mutations include fibroadipose overgrowth (FAO), Hemihyperplasia multiple Lipomatosis (HHML), Congenital Lipomatous Overgrowth, Vascular malformations, Epidermal nevi, Scoliosis/skeletal and spinal (CLOVES) syndrome, macrodactyly, and the megalencephaly syndrome, Megalencephaly-Capillary malformation (MCAP) syndrome. We set out to refine the understanding of the clinical spectrum and natural history of these phenotypes, and now describe 35 patients with segmental overgrowth and somatic PIK3CA mutations. The phenotypic data show that these previously described disease entities have considerable overlap, and represent a spectrum. While this spectrum overlaps with Proteus syndrome (sporadic, mosaic, and progressive) it can be distinguished by the absence of cerebriform connective tissue nevi and a distinct natural history. Vascular malformations were found in 15/35 (43%) and epidermal nevi in 4/35 (11%) patients, lower than in Proteus syndrome. Unlike Proteus syndrome, 31/35 (89%) patients with PIK3CA mutations had congenital overgrowth, and in 35/35 patients this was asymmetric and disproportionate. Overgrowth was mild with little postnatal progression in most, while in others it was severe and progressive requiring multiple surgeries. Novel findings include: adipose dysregulation present in all patients, unilateral overgrowth that is predominantly left-sided, overgrowth that affects the lower extremities more than the upper extremities and progresses in a distal to proximal pattern, and in the most severely affected patients is associated with marked paucity of adipose tissue in unaffected areas. While the current data are consistent with some genotype-phenotype correlation, this cannot yet be confirmed.