TF and craniosynostosis: These TFs were determined based on over-representation of TF binding sites in regulatory regions of craniosynostosis DEGs, measured through Enrichr (Figure 4A, FDR < 0.05).