The RUNT-related transcription factor RUNX2 plays a critical role in osteoblast differentiation, and alterations to gene dosage cause distinct craniofacial anomalies. Uniquely amongst the RUNT-related family, vertebrate RUNX2 encodes a polyglutamine/polyalanine repeat (Gln23-Glu-Ala17 in humans), with the length of the polyalanine component completely conserved in great apes. Surprisingly, a frequent 6-amino acid deletion polymorphism, p.(Ala84_Ala89)del, occurs in humans (termed 11A allele), and a previous association study (Cuellar et al. Bone 137:115395;2020) reported that the 11A variant was significantly more frequent in non-syndromic sagittal craniosynostosis (nsSag; allele frequency [AF] = 0.156; 95% confidence interval [CI] 0.126-0.189) compared to non-syndromic metopic craniosynostosis (nsMet; AF = 0.068; 95% CI 0.045-0.098). However, the gnomAD v.2.1.1 control population used by Cuellar et al. did not display Hardy-Weinberg equilibrium, hampering interpretation. To re-examine this association, we genotyped the RUNX2 11A polymorphism in 225 individuals with sporadic nsSag as parent-child trios and 164 singletons with sporadic nsMet, restricting our analysis to individuals of European ancestry. We compared observed allele frequencies to the non-transmitted alleles in the parent-child trios, and to the genome sequencing data from gnomAD v.4, which display Hardy-Weinberg equilibrium. Observed AFs (and 95% CI) were 0.076 (0.053-0.104) in nsSag and 0.082 (0.055-0.118) in nsMet, compared with 0.062 (0.042-0.089) in non-transmitted parental alleles and 0.065 (0.063-0.067) in gnomAD v.4.0.0 non-Finnish European control genomes. In summary, we observed a non-significant excess, compared to gnomAD data, of 11A alleles in both nsSag (relative risk 1.18, 95% CI 0.83-1.67) and nsMet (relative risk 1.29, 95% CI 0.87-1.92), but we did not replicate the much higher excess of RUNX2 11A alleles in nsSag previously reported (p = 0.0001).

Craniosynostosis is a condition with neurologic and aesthetic sequelae requiring invasive surgery. Understanding its pathobiology requires familiarity with the processes underlying physiologic suture closure. Animal studies have shown that cyclical strain from chewing and suckling influences the closure of cranial vault sutures, especially the metopic, an important locus of craniosynostosis. However, there are no human data correlating strain patterns during chewing and suckling with the physiologically early closure pattern of the metopic suture. Furthermore, differences in craniofacial morphology make it challenging to directly extrapolate animal findings to humans. Eight finite-element analysis (FEA) models were built from craniofacial computer tomography (CT) scans at varying stages of metopic suture closure, including two with isolated non-syndromic metopic craniosynostosis. Muscle forces acting on the cranium during chewing and suckling were simulated using subject-specific jaw muscle cross-sectional areas. Chewing and suckling induced tension at the metopic and sagittal sutures, and compressed the coronal, lambdoid, and squamous sutures. Relative to other cranial vault sutures, the metopic suture experienced larger magnitudes of axial strain across the suture and a lower magnitude of shear strain. Strain across the metopic suture decreased during suture closure, but other sutures were unaffected. Strain patterns along the metopic suture mirrored the anterior to posterior sequence of closure: strain magnitudes were highest at the glabella and decreased posteriorly, with minima at the nasion and the anterior fontanelle. In models of physiologic suture closure, increased degree of metopic suture closure correlated with higher maximum principal strains across the frontal bone and mid-face, a strain regime not observed in models of severe metopic craniosynostosis. In summary, our work provides human evidence that bone strain patterns from chewing and suckling correlate with the physiologically early closure pattern of the metopic suture, and that deviations from physiologic strain regimes may contribute to clinically observed craniofacial dysmorphism.

Craniofacial dysmorphology varies significantly along a wide spectrum of severity in metopic cranial synostosis (MCS). This study aimed to quantify craniofacial changes, in MCS, to investigate their relationships with the severity of trigonocephaly. By combining the metopic ridge and interfrontal angles, we identified three groups of trigonocephaly severity (mild group n.14, moderate group n.19, severe group n.18). We perform a quantitative analysis using high-resolution CT images evaluating (1) cranial fossae dimensions; (2) vault indices and ratios: interparietal/ intercoronal (IPD/ICD), interparietal/intertemporal (IPD/ITD), cephalic index, vertico-longitudinal index; (3) orbito-facial distances (midfacial depth, maxillary height, upper facial index, orbital distances, globe protrusions), maxilla and orbital volumes; (4) supratentorial (ICV) and infratentorial (PCFV) cranial volumes and supratentorial (WBV) and infratentorial (PCFBV) brain volumes. In all groups, middle skull base lengths and upper midface index were increased. In moderate and severe groups: anterior hemifossa lengths were reduced, IPD/ICD and vertico-longitudinal index were changed; midfacial depth, anterior, mild, and lateral interorbital distances were reduced; globe protrusions were increased. The comparison between moderate and severe groups showed an increase of both globe protrusions and IPD/ICD. Among all groups, ICV and WBV were reduced in the severe group. This morpho-volumetric study provides new insights in understanding the craniofacial changes occurring in infants at different severity of trigonocephaly. The increase of globe protrusions and the reduction of supratentorial volumes found in the severe group reflect the severity of trigonocephaly; these findings might have a clinical and surgical relevance.