Extensive morphological variation found in mammals reflects the wide spectrum of their ecological adaptations. The highest morphological diversity is present in the craniofacial region, where geometry is mainly dictated by the bony skull. Mammalian craniofacial development represents complex multistep processes governed by numerous conserved genes that require precise spatio-temporal control. A central question in contemporary evolutionary biology is how a defined set of conserved genes can orchestrate formation of fundamentally different structures, and therefore how morphological variability arises. In principle, differential gene expression patterns during development are the source of morphological variation. With the emergence of multicellular organisms, precise regulation of gene expression in time and space is attributed to cis-regulatory elements. These elements contribute to higher-order chromatin structure and together with trans-acting factors control transcriptional landscapes that underlie intricate morphogenetic processes. Consequently, divergence in cis-regulation is believed to rewire existing gene regulatory networks and form the core of morphological evolution. This review outlines the fundamental principles of the genetic code and genomic regulation interplay during development. Recent work that deepened our comprehension of cis-regulatory element origin, divergence and function is presented here to illustrate the state-of-the-art research that uncovered the principles of morphological novelty. This article is part of the theme issue \'The mammalian skull: development, structure and function\'.

BACKGROUND: We compared skull shape and variation among genetically modified mice that exhibit different levels of connexin43 (Cx43) channel function, to determine whether Cx43 contributes to craniofacial phenotypic robustness. Specifically, we used two heterozygous mutant mouse models (G60S/+ and I130T/+) that, when compared to their wildtype counterparts, have an ~80% and ~50% reduction in Cx43 function, respectively.
RESULTS: Both mutant strains showed significant differences in skull shape compared to wildtype littermates and while these differences were more severe in the G60S/+ mouse, shape differences were localized to similar regions of the skull in both mutants. However, increased skull shape variation was observed in G60S/+ mutants only. Additionally, covariation of skull structures was disrupted in the G60S/+ mutants only, indicating that while a 50% reduction in Cx43 function is sufficient to cause a shift in mean skull shape, the threshold for Cx43 function for disrupting craniofacial phenotypic robustness is lower.
CONCLUSIONS: Collectively, our results indicate Cx43 can contribute to phenotypic robustness of the skull through a nonlinear relationship between Cx43 gap junctional function and phenotypic outcomes.

Mutations in the gene encoding the gap-junctional protein connexin43 (Cx43) are the cause of the human disease oculodentodigital dysplasia (ODDD). The mandible is often affected in this disease, with clinical reports describing both mandibular overgrowth and conversely, retrognathia. These seemingly opposing observations underscore our relative lack of understanding of how ODDD affects mandibular morphology. Using two mutant mouse models that mimic the ODDD phenotype (I130T/+ and G60S/+), we sought to uncover how altered Cx43 function may affect mandibular development. Specifically, mandibles of newborn mice were imaged using micro-CT, to enable statistical comparisons of shape. Tissue-level comparisons of key regions of the mandible were conducted using histomorphology, and we quantified the mRNA expression of several cartilage and bone cell differentiation markers. Both G60S/+ and I130T/+ mutant mice had altered mandibular morphology compared to their wildtype counterparts, and the morphological effects were similarly localized for both mutants. Specifically, the biggest phenotypic differences in mutant mice were focused in regions exposed to mechanical forces, such as alveolar bone, muscular attachment sites, and articular surfaces. Histological analyses revealed differences in ossification of the intramembranous bone of the mandibles of both mutant mice compared to their wildtype littermates. However, chondrocyte organization within the secondary cartilages of the mandible was unaffected in the mutant mice. Overall, our results suggest that the morphological differences seen in G60S/+ and I130T/+ mouse mandibles are due to delayed ossification and suggest that mechanical forces may exacerbate the effects of ODDD on the skeleton.