Craniosynostosis is the early fusion of one or more sutures of the infant skull and is a common defect occurring in approximately 1 of every 2,500 live births. Non-syndromic craniosynostosis accounts for approximately 80% of all cases and is thought to have strong genetic determinants that are yet to be identified. ALX4 is a homeodomain transcription factor with known involvement in osteoblast regulation. By direct sequencing of the ALX4 coding region in sagittal or sagittal-suture-involved non-syndromic craniosynostosis probands, we identified novel, nonsynonymous, familial variants in three of 203 individuals with NSC. Using dual-luciferase assay we show that two of these variants (V7F and K211E) confer a significant gain-of-function effect on ALX4. Our results suggest that ALX4 variants may have an impact on the genetic etiology of NSC.

Craniosynostosis (CS) is a major birth defect resulting from premature fusion of cranial sutures. Nonsyndromic CS occurs more frequently than syndromic CS, with sagittal nonsyndromic craniosynostosis (sNCS) presenting as the most common CS phenotype. Previous genome-wide association and targeted sequencing analyses of sNCS have identified multiple associated loci, with the strongest association on chromosome 20. Herein, we report the first whole-genome sequencing study of sNCS using 63 proband-parent trios. Sequencing data for these trios were analyzed using the transmission disequilibrium test (TDT) and rare variant TDT (rvTDT) to identify high-risk rare gene variants. Sequencing data were also examined for copy number variants (CNVs) and de novo variants. TDT analysis identified a highly significant locus at 20p12.3, localized to the intergenic region between BMP2 and the noncoding RNA gene LINC01428. Three variants (rs6054763, rs6054764, rs932517) were identified as potential causal variants due to their probability of being transcription factor binding sites, deleterious combined annotation dependent depletion scores, and high minor allele enrichment in probands. Morphometric analysis of cranial vault shape in an unaffected cohort validated the effect of these three single nucleotide variants (SNVs) on dolichocephaly. No genome-wide significant rare variants, de novo loci, or CNVs were identified. Future efforts to identify risk variants for sNCS should include sequencing of larger and more diverse population samples and increased omics analyses, such as RNA-seq and ATAC-seq.

Sagittal craniosynostosis is the most common form of craniosynostosis, affecting approximately one in 5,000 newborns. We conducted, to our knowledge, the first genome-wide association study for nonsyndromic sagittal craniosynostosis (sNSC) using 130 non-Hispanic case-parent trios of European ancestry (NHW). We found robust associations in a 120-kb region downstream of BMP2 flanked by rs1884302 (P = 1.13 × 10(-14), odds ratio (OR) = 4.58) and rs6140226 (P = 3.40 × 10(-11), OR = 0.24) and within a 167-kb region of BBS9 between rs10262453 (P = 1.61 × 10(-10), OR = 0.19) and rs17724206 (P = 1.50 × 10(-8), OR = 0.22). We replicated the associations to both loci (rs1884302, P = 4.39 × 10(-31) and rs10262453, P = 3.50 × 10(-14)) in an independent NHW population of 172 unrelated probands with sNSC and 548 controls. Both BMP2 and BBS9 are genes with roles in skeletal development that warrant functional studies to further understand the etiology of sNSC.

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.