LBX1 and craniosynostosis: The key findings of this study were that, (1) each craniosynostosis phenotype had at least 20 unique differentially expressed genes (DEGs), (2) more DEGs were identified in the male-stratified analysis compared to the female-stratified analysis, and (3) that homeobox genes were prevalent as DEGs across the primary and male-stratified models.