Craniosynostosis (CRS) is characterized by the development of abnormal cranial suture ossification and premature fusion. Despite the identification of several associated genetic disorders, the genetic determinants of CRS remain poorly understood. In this study, we conducted integrative analyses on 225 exomes, comprising 121 CRS probands and 104 parental exomes (52 trios). These analyses encompassed de novo and pathogenic variants, and digenic combinations within haploinsufficient genes harboring rare variants. Our analysis unveils a shared molecular network between genes associated with CRS and those linked to skeletal and neurodevelopmental disorders, with a notable enrichment of deleterious variants within haploinsufficient genes. Additionally, we identified a unique digenic pair (IL6ST and TRPS1) within haploinsufficient genes that was present in 2 patients with nonsyndromic CRS but absent in parents or 1,048 population controls. In vitro experiments provided evidence that the identified missense variants were hypomorphs, and accelerated bone mineralization could result from the additive effects of diminished IL6ST and TRPS1 activities in osteoblasts. Overall, our study underscores the important role of rare variations in haploinsufficient genes and suggests that in a subset of undiagnosed patients, the CRS phenotype may arise from multiple genetic variations.