A coupled reaction–diffusion–strain model predicts cranial vault formation in development and disease

How cells utilize instructions provided by genes and integrate mechanical forces generated by tissue growth to produce morphology is a fundamental question of biology. Dermal bones of the vertebrate cranial vault are formed through the direct differentiation of mesenchymal cells on the neural surface into osteoblasts through intramembranous ossification. Here we join a self-organizing Turing mechanism, computational biomechanics, and experimental data to produce a 3D representative model of the growing cerebral surface, cranial vault bones, and sutures. We show how changes in single parameters regulating signaling during osteoblast differentiation and bone formation may explain cranial vault shape variation in craniofacial disorders. A key result is that toggling a parameter in our model results in closure of a cranial vault suture, an event that occurred during evolution of the cranial vault and that occurs in craniofacial disorders. Our approach provides an initial and important step toward integrating biomechanics into the genotype phenotype map to explain the production of variation in head morphology by developmental mechanisms.