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.

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.

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.

Identification of differentially expressed proteins between fused and open sutures in sagittal nonsyndromic craniosynostosis during suture development by quantitative proteomic analysis.

PURPOSE: Nonsyndromic craniosynostosis (NCS), the premature fusion of cranial sutures, results in an abnormal skull shape and is associated with a significant morbidity. Proteomics is a promising tool for disease characterization and biomarker discovery; we aimed to identify biologically relevant differentially expressed proteins for NCS. EXPERIMENTAL DESIGN: Label-based quantitative proteomic profiling using TMT was performed on protein extracted from mesenchymal stem cells, osteoblasts and bone tissue of five open and five fused sutures of sagittal NCS (sNCS) and analyzed using quantitative LC-MS/MS based bottom-up proteomics. Differential protein abundance between open and fused sutures was determined to identify biologically relevant proteins of interest. Proteins were validated in an independent sample set by western blot and immunohistochemistry. RESULTS: We observed 838 differentially expressed proteins between open and fused sutures of sNCS. Decorin, lumican, and asporin were significantly downregulated while COL4A1 and TGFß1|1 were upregulated in fused compared to open sutures. CONCLUSIONS AND CLINICAL RELEVANCE: The majority of significantly differentially expressed proteins between open and fused sutures were observed in the proteomes of osteoblasts suggesting that protein changes contributing to premature sagittal suture fusion occur predominantly at the osteoblast level. Our findings suggest a possible ineffective ECM deposition at the osteoblast cell stage.

Gain-of-function variants and overexpression of RUNX2 in patients with nonsyndromic midline craniosynostosis.

Craniosynostosis (CS), the premature fusion of one or more cranial sutures, is a relatively common congenital anomaly, occurring in 3-5 per 10,000 live births. Nonsyndromic CS (NCS) accounts for up to 80% of all CS cases, yet the genetic factors contributing to the disorder remain largely unknown. The RUNX2 gene, encoding a transcription factor critical for bone and skull development, is a well known CS candidate gene, as copy number variations of this gene locus have been found in patients with syndromic craniosynostosis. In the present study, we aimed to characterize RUNX2 to better understand its role in the genetic etiology and in the molecular mechanisms underlying midline suture ossification in NCS. We report four nonsynonymous variants, one intronic variant and one 18 bp in-frame deletion in RUNX2 not found in our study control population. Significant difference in allele frequency (AF) for the deletion variant RUNX2 p.Ala84-Ala89del (ClinVar 257,095; dbSNP rs11498192) was observed in our sagittal NCS cohort when compared to the general population (P = 1.28 × 10-6), suggesting a possible role in the etiology of NCS. Dual-luciferase assays showed that three of four tested RUNX2 variants conferred a gain-of-function effect on RUNX2, further suggesting their putative pathogenicity in the tested NCS cases. Downregulation of RUNX2 expression was observed in prematurely ossified midline sutures. Metopic sites showed significant downregulation of promoter 1-specific isoforms compared to sagittal sites. Suture-derived mesenchymal stromal cells showed an increased expression of RUNX2 over matched unfused suture derived cells. This demonstrates that RUNX2, and particularly the distal promoter 1-isoform group, are overexpressed in the osteogenic precursors within the pathological suture sites.