Quantitative analyses link modulation of sonic hedgehog signaling to continuous variation in facial growth and shape

PBX1 and PBX3 transcription factors regulate SHH expression in the Frontonasal Ectodermal Zone through complementary mechanisms

Sonic hedgehog (SHH) signaling from the frontonasal ectodermal zone (FEZ) is a key regulator of craniofacial morphogenesis. Along with SHH, pre-B-cell leukemia homeobox (PBX) transcription factors regulate midfacial development. PBXs act in the epithelium during fusion of facial primordia, but their specific interactions with SHH have not been fully investigated. We hypothesized that PBX1/3 regulate SHH expression in the FEZ by activating or repressing transcription. The hypothesis was tested by manipulating PBX1/3 expression in chick embryos and profiling epigenomic landscapes at early developmental stages. PBX1/3 expression was perturbed in the chick face beginning at stage 10 (HH10) using RCAS viruses, and the resulting SHH expression was assessed at HH22. Overexpressing PBX1 expanded SHH expression, while overexpressing PBX3 decreased SHH expression. Conversely, reducing PBX1 expression decreased SHH expression, but reducing PBX3 induced ectopic SHH expression. We performed ATAC-seq and mapped binding of PBX1 and PBX3 with ChIP-seq on the FEZ at HH22 to assess direct interactions of PBX1/3 with the SHH locus. These multi-omics approaches uncovered a 400 bp PBX1-enriched element within intron 1 of SHH (chr2:8,173,222-8,173,621). Enhancer activity of this element was demonstrated by electroporation of reporter constructs in ovo and luciferase reporter assays in vitro . When bound by PBX1, this element upregulates transcription, while it downregulates transcription when bound by PBX3. The present study identifies a cis- regulatory element, named SFE1, that interacts with PBX1/3 to modulate SHH expression in the FEZ and establishes that PBX1 and PBX3 play complementary roles in SHH regulation during embryonic development.

Shaping faces: genetic and epigenetic control of craniofacial morphogenesis

Major differences in facial morphology distinguish vertebrate species. Variation of facial traits underlies the uniqueness of human individuals, and abnormal craniofacial morphogenesis during development leads to birth defects that significantly affect quality of life. Studies during the past 40 years have advanced our understanding of the molecular mechanisms that establish facial form during development, highlighting the crucial roles in this process of a multipotent cell type known as the cranial neural crest cell. In this Review, we discuss recent advances in multi-omics and single-cell technologies that enable genes, transcriptional regulatory networks and epigenetic landscapes to be closely linked to the establishment of facial patterning and its variation, with an emphasis on normal and abnormal craniofacial morphogenesis. Advancing our knowledge of these processes will support important developments in tissue engineering, as well as the repair and reconstruction of the abnormal craniofacial complex.

Zebrafish crocc2 mutants exhibit divergent craniofacial shape, misregulated variability, and aberrant cartilage morphogenesis

BACKGROUND

Phenotypic variation is of paramount importance in development, evolution, and human health; however, the molecular mechanisms that influence organ shape and shape variability are not well understood. During craniofacial development, the behavior of skeletal precursors is regulated by both biochemical and environmental inputs, and the primary cilia play critical roles in transducing both types of signals. Here, we examine a gene that encodes a key constituent of the ciliary rootlets, crocc2, and its role in cartilage morphogenesis in larval zebrafish.

RESULTS

Geometric morphometric analysis of crocc2 mutants revealed altered craniofacial shapes and expanded variation. At the cellular level, we observed altered chondrocyte shapes and planar cell polarity across multiple stages in crocc2 mutants. Notably, cellular defects were specific to areas that experience direct mechanical input. Cartilage cell number, apoptosis, and bone patterning were not affected in crocc2 mutants.

CONCLUSIONS

Whereas "regulatory" genes are widely implicated in patterning the craniofacial skeleton, genes that encode "structural" aspects of the cell are increasingly implicated in shaping the face. Our results add crocc2 to this list, and demonstrate that it affects craniofacial geometry and canalizes phenotypic variation. We propose that it does so via mechanosensing, possibly through the ciliary rootlet. If true, this would implicate a new organelle in skeletal development and evolution.

Fgf8 dosage regulates jaw shape and symmetry through pharyngeal-cardiac tissue relationships

BACKGROUND

Asymmetries in craniofacial anomalies are commonly observed. In the facial skeleton, the left side is more commonly and/or severely affected than the right. Such asymmetries complicate treatment options. Mechanisms underlying variation in disease severity between individuals as well as within individuals (asymmetries) are still relatively unknown.

RESULTS

Developmental reductions in fibroblast growth factor 8 (Fgf8) have a dosage dependent effect on jaw size, shape, and symmetry. Further, Fgf8 mutants have directionally asymmetric jaws with the left side being more affected than the right. Defects in lower jaw development begin with disruption to Meckel's cartilage, which is discontinuous. All skeletal elements associated with the proximal condensation are dysmorphic, exemplified by a malformed and misoriented malleus. At later stages, Fgf8 mutants exhibit syngnathia, which falls into two broad categories: bony fusion of the maxillary and mandibular alveolar ridges and zygomatico-mandibular fusion. All of these morphological defects exhibit both inter- and intra-specimen variation.

CONCLUSIONS

We hypothesize that these asymmetries are linked to heart development resulting in higher levels of Fgf8 on the right side of the face, which may buffer the right side to developmental perturbations. This mouse model may facilitate future investigations of mechanisms underlying human syngnathia and facial asymmetry.

MusMorph, a database of standardized mouse morphology data for morphometric meta-analyses

Complex morphological traits are the product of many genes with transient or lasting developmental effects that interact in anatomical context. Mouse models are a key resource for disentangling such effects, because they offer myriad tools for manipulating the genome in a controlled environment. Unfortunately, phenotypic data are often obtained using laboratory-specific protocols, resulting in self-contained datasets that are difficult to relate to one another for larger scale analyses. To enable meta-analyses of morphological variation, particularly in the craniofacial complex and brain, we created MusMorph, a database of standardized mouse morphology data spanning numerous genotypes and developmental stages, including E10.5, E11.5, E14.5, E15.5, E18.5, and adulthood. To standardize data collection, we implemented an atlas-based phenotyping pipeline that combines techniques from image registration, deep learning, and morphometrics. Alongside stage-specific atlases, we provide aligned micro-computed tomography images, dense anatomical landmarks, and segmentations (if available) for each specimen (N = 10,056). Our workflow is open-source to encourage transparency and reproducible data collection. The MusMorph data and scripts are available on FaceBase ( www.facebase.org , https://doi.org/10.25550/3-HXMC ) and GitHub ( https://github.com/jaydevine/MusMorph ).

Concepts in Multifactorial Etiology of Developmental Disorders: Gene-Gene and Gene-Environment Interactions in Holoprosencephaly

Many common developmental disorders are thought to arise from a complex set of genetic and environmental risk factors. These factors interact with each other to affect the strength and duration of key developmental signaling pathways, thereby increasing the possibility that they fail to achieve the thresholds required for normal embryonic patterning. One such disorder, holoprosencephaly (HPE), serves as a useful model system in understanding various forms of multifactorial etiology. Genomic analysis of HPE cases, epidemiology, and mechanistic studies of animal models have illuminated multiple potential ways that risk factors interact to produce adverse developmental outcomes. Among these are: 1) interactions between driver and modifier genes; 2) oligogenic inheritance, wherein each parent provides predisposing variants in one or multiple distinct loci; 3) interactions between genetic susceptibilities and environmental risk factors that may be insufficient on their own; and 4) interactions of multiple genetic variants with multiple non-genetic risk factors. These studies combine to provide concepts that illuminate HPE and are also applicable to additional disorders with complex etiology, including neural tube defects, congenital heart defects, and oro-facial clefting.

Environmental Thermal Stress Induces Neuronal Cell Death and Developmental Malformations in Reptiles

Every stage of organismal life history is being challenged by global warming. Many species are already experiencing temperatures approaching their physiological limits; this is particularly true for ectothermic species, such as lizards. Embryos are markedly sensitive to thermal insult. Here, we demonstrate that temperatures currently experienced in natural nesting areas can modify gene expression levels and induce neural and craniofacial malformations in embryos of the lizard Anolis sagrei. Developmental abnormalities ranged from minor changes in facial structure to significant disruption of anterior face and forebrain. The first several days of postoviposition development are particularly sensitive to this thermal insult. These results raise new concern over the viability of ectothermic species under contemporary climate change. Herein, we propose and test a novel developmental hypothesis that describes the cellular and developmental origins of those malformations: cell death in the developing forebrain and abnormal facial induction due to disrupted Hedgehog signaling. Based on similarities in the embryonic response to thermal stress among distantly related species, we propose that this developmental hypothesis represents a common embryonic response to thermal insult among amniote embryos. Our results emphasize the importance of adopting a broad, multidisciplinary approach that includes both lab and field perspectives when trying to understand the future impacts of anthropogenic change on animal development.

Relating multivariate shapes to genescapes using phenotype-biological process associations for craniofacial shape

Realistic mappings of genes to morphology are inherently multivariate on both sides of the equation. The importance of coordinated gene effects on morphological phenotypes is clear from the intertwining of gene actions in signaling pathways, gene regulatory networks, and developmental processes underlying the development of shape and size. Yet, current approaches tend to focus on identifying and localizing the effects of individual genes and rarely leverage the information content of high-dimensional phenotypes. Here, we explicitly model the joint effects of biologically coherent collections of genes on a multivariate trait – craniofacial shape – in a sample of n = 1145 mice from the Diversity Outbred (DO) experimental line. We use biological process Gene Ontology (GO) annotations to select skeletal and facial development gene sets and solve for the axis of shape variation that maximally covaries with gene set marker variation. We use our process-centered, multivariate genotype-phenotype (process MGP) approach to determine the overall contributions to craniofacial variation of genes involved in relevant processes and how variation in different processes corresponds to multivariate axes of shape variation. Further, we compare the directions of effect in phenotype space of mutations to the primary axis of shape variation associated with broader pathways within which they are thought to function. Finally, we leverage the relationship between mutational and pathway-level effects to predict phenotypic effects beyond craniofacial shape in specific mutants. We also introduce an online application that provides users the means to customize their own process-centered craniofacial shape analyses in the DO. The process-centered approach is generally applicable to any continuously varying phenotype and thus has wide-reaching implications for complex trait genetics.

Ciliary Hedgehog signaling regulates cell survival to build the facial midline

Craniofacial defects are among the most common phenotypes caused by ciliopathies, yet the developmental and molecular etiology of these defects is poorly understood. We investigated multiple mouse models of human ciliopathies (including Tctn2, Cc2d2a, and Tmem231 mutants) and discovered that each displays hypotelorism, a narrowing of the midface. As early in development as the end of gastrulation, Tctn2 mutants displayed reduced activation of the Hedgehog (HH) pathway in the prechordal plate, the head organizer. This prechordal plate defect preceded a reduction of HH pathway activation and Shh expression in the adjacent neurectoderm. Concomitant with the reduction of HH pathway activity, Tctn2 mutants exhibited increased cell death in the neurectoderm and facial ectoderm, culminating in a collapse of the facial midline. Enhancing HH signaling by decreasing the gene dosage of a negative regulator of the pathway, Ptch1, decreased cell death and rescued the midface defect in both Tctn2 and Cc2d2a mutants. These results reveal that ciliary HH signaling mediates communication between the prechordal plate and the neurectoderm to provide cellular survival cues essential for development of the facial midline.

Novel neuroanatomical integration and scaling define avian brain shape evolution and development

How do large and unique brains evolve? Historically, comparative neuroanatomical studies have attributed the evolutionary genesis of highly encephalized brains to deviations along, as well as from, conserved scaling relationships among brain regions. However, the relative contributions of these concerted (integrated) and mosaic (modular) processes as drivers of brain evolution remain unclear, especially in non-mammalian groups. While proportional brain sizes have been the predominant metric used to characterize brain morphology to date, we perform a high-density geometric morphometric analysis on the encephalized brains of crown birds (Neornithes or Aves) compared to their stem taxa—the non-avialan coelurosaurian dinosaurs and Archaeopteryx. When analyzed together with developmental neuroanatomical data of model archosaurs (Gallus, Alligator), crown birds exhibit a distinct allometric relationship that dictates their brain evolution and development. Furthermore, analyses by neuroanatomical regions reveal that the acquisition of this derived shape-to-size scaling relationship occurred in a mosaic pattern, where the avian-grade optic lobe and cerebellum evolved first among non-avialan dinosaurs, followed by major changes to the evolutionary and developmental dynamics of cerebrum shape after the origin of Avialae. Notably, the brain of crown birds is a more integrated structure than non-avialan archosaurs, implying that diversification of brain morphologies within Neornithes proceeded in a more coordinated manner, perhaps due to spatial constraints and abbreviated growth period. Collectively, these patterns demonstrate a plurality in evolutionary processes that generate encephalized brains in archosaurs and across vertebrates.

Connexin 43 contributes to phenotypic robustness of the mouse skull

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.

Differences in Cell Proliferation and Craniofacial Phenotype of Closely Related Species in the Pupfish Genus Cyprinodon

Understanding the genetic basis for phenotypic differences is fundamental to the study of macroevolutionary patterns of biological diversity. While technological advances in DNA sequencing have made researching genetic variation in wild taxa routine, fully understanding how these variants affect phenotype requires taking the next step to investigate how genetic changes alter cell and tissue interactions that ultimately produce phenotypes. In this article, we investigate a role for cell proliferation as a developmental source of craniofacial diversity in a radiation of 3 species of Cyprinodon from San Salvador Island, Bahamas. Patterns of cell proliferation in the heads of hatching-age fish differ among species of Cyprinodon, and correlate with differences in allometric growth rate among the jaws of 3 distinct species. Regional patterns of cell proliferation in the head are complex, resulting in an unintuitive result in which lower levels of cell proliferation in the posterior head region are associated with the development of relatively larger jaws in one species. We combine these data with previously published morphological and genomic data to show how studying the mechanisms generating phenotype at the cellular and tissue levels of biological organization can help mechanistically link genomic studies with classic morphological studies.

Symmetry and fluctuation of cell movements in neural crest-derived facial mesenchyme

In the face, symmetry is established when bilateral streams of neural crest cells leave the neural tube at the same time, follow identical migration routes and then give rise to the facial prominences. However, developmental instability exists, particularly surrounding the steps of lip fusion. The causes of instability are unknown but inability to cope with developmental fluctuations are a likely cause of congenital malformations, such as non-syndromic orofacial clefts. Here, we tracked cell movements over time in the frontonasal mass, which forms the facial midline and participates in lip fusion, using live-cell imaging of chick embryos. Our mathematical examination of cell velocity vectors uncovered temporal fluctuations in several parameters, including order/disorder, symmetry/asymmetry and divergence/convergence. We found that treatment with a Rho GTPase inhibitor completely disrupted the temporal fluctuations in all measures and blocked morphogenesis. Thus, we discovered that genetic control of symmetry extends to mesenchymal cell movements and that these movements are of the type that could be perturbed in asymmetrical malformations, such as non-syndromic cleft lip.

This article has an associated ‘The people behind the papers’ interview.

Highlighted Article: Live imaging of the chick embryo face followed by mathematical analysis of mesenchymal cell tracks captures novel fluctuations between states of order/disorder as well as symmetry/asymmetry, revealing developmental instabilities that are part of normal morphogenesis.

Genetic variations at 10q26 regions near FGFR2 gene and its association with non-syndromic cleft lip with or without cleft palate

OBJECTIVES

In our study, we focussed on three SNPs in the non-coding regions near FGFR2 gene, as studies on non-coding variants in the genome are the novel trends to identify the susceptible loci for nonsyndromic cleft lip with or without cleft palate (NSCL/P). FGFR2 gene is selected as a candidate gene based on knock out animal models and the role played in syndromic forms of clefting. FGFR2 gene also plays an important role in FGF signaling pathway during craniofacial development.

METHODS

In the present study 148 case-parent triads were assessed for three SNPs rs10749408, rs11199874 and rs10788165 near FGFR2 gene by using TaqMan allelic discrimination method. Transmission disequilibrium test (TDT) was used to find the allelic association. Linkage disequilibrium (LD) between the markers was analysed using Haploview program 4.2. Haplotype transmission effects were estimated using FAMHAP package. The possible parent-of-origin effects were assessed by likelihood based approach.

RESULTS

TDT analysis of three SNPs failed to show significant transmission disortion from heterozygous parents to the affected child and are not associated with NSCL/P. Linkage disequilibrium analysis showed strong LD between rs11199874 and rs10788165 SNPs. In the haplotype TDT analysis, GG haplotype of rs11199874-rs10788165 showed significant undertransmission to affected child. No significant parent-of-origin effects were observed.

CONCLUSION

The present study on noncoding variants near FGFR2 gene is not associated with NSCL/P. As the numbers of triads included in the study are less, further studies are needed including large sample size to find association with NSCL/P.

Human-chimpanzee fused cells reveal cis-regulatory divergence underlying skeletal evolution

Gene regulatory divergence is thought to play a central role in determining human-specific traits. However, our ability to link divergent regulation to divergent phenotypes is limited. Here, we utilized human-chimpanzee hybrid induced pluripotent stem cells to study gene expression separating these species. The tetraploid hybrid cells allowed us to separate cis- from trans-regulatory effects, and to control for non-genetic confounding factors. We differentiated these cells into cranial neural crest cells (CNCCs), the primary cell type giving rise to the face. We discovered evidence of lineage-specific selection on the hedgehog signaling pathway, including a human-specific 6-fold down-regulation of EVC2 (LIMBIN), a key hedgehog gene. Inducing a similar down-regulation of EVC2 substantially reduced hedgehog signaling output. Mice and humans lacking functional EVC2 show striking phenotypic parallels to human-chimpanzee craniofacial differences, suggesting that the regulatory divergence of hedgehog signaling may have contributed to the unique craniofacial morphology of humans.

MicroRNA function in craniofacial bone formation, regeneration and repair

Bone formation in the craniofacial complex is regulated by cranial neural crest (CNC) and mesoderm-derived cells. Different elements of the developing skull, face, mandible, maxilla (jaws) and nasal bones are regulated by an array of transcription factors, signaling molecules and microRNAs (miRs). miRs are molecular modulators of these factors and act to restrict their expression in a temporal-spatial mechanism. miRs control the different genetic pathways that form the craniofacial complex. By understanding how miRs function in vivo during development they can be adapted to regenerate and repair craniofacial genetic anomalies as well as bone diseases and defects due to traumatic injuries. This review will highlight some of the new miR technologies and functions that form new bone or inhibit bone regeneration.

Gli3 utilizes Hand2 to synergistically regulate tissue-specific transcriptional networks

Despite a common understanding that Gli TFs are utilized to convey a Hh morphogen gradient, genetic analyses suggest craniofacial development does not completely fit this paradigm. Using the mouse model (Mus musculus), we demonstrated that rather than being driven by a Hh threshold, robust Gli3 transcriptional activity during skeletal and glossal development required interaction with the basic helix-loop-helix TF Hand2. Not only did genetic and expression data support a co-factorial relationship, but genomic analysis revealed that Gli3 and Hand2 were enriched at regulatory elements for genes essential for mandibular patterning and development. Interestingly, motif analysis at sites co-occupied by Gli3 and Hand2 uncovered mandibular-specific, low-affinity, 'divergent' Gli-binding motifs (dGBMs). Functional validation revealed these dGBMs conveyed synergistic activation of Gli targets essential for mandibular patterning and development. In summary, this work elucidates a novel, sequence-dependent mechanism for Gli transcriptional activity within the craniofacial complex that is independent of a graded Hh signal.

Functional analysis of Sonic Hedgehog variants associated with holoprosencephaly in humans using a CRISPR/Cas9 zebrafish model

Genetic variation in the highly conserved Sonic Hedgehog (SHH) gene is one of the most common genetic causes for the malformations of the brain and face in humans described as the holoprosencephaly clinical spectrum. However, only a minor fraction of known SHH variants have been experimentally proven to lead to abnormal function. Employing a phenotypic rescue assay with synthetic human messenger RNA variant constructs in shha^-/- knockout zebrafish, we evaluated 104 clinically reported in-frame and missense SHH variants. Our data helped us to classify them into loss of function variants (31), hypomorphic variants (33), and nonpathogenic variants (40). We discuss the strengths and weaknesses of currently accepted predictors of variant deleteriousness and the American College of Medical Genetics and Genomics guidelines for variant interpretation in the context of this functional model; furthermore, we demonstrate the robustness of model systems such as zebrafish as a rapid method to resolve variants of uncertain significance.

miR-199 family contributes to regulation of sonic hedgehog expression during craniofacial development

BACKGROUND

The frontonasal ectodermal zone (FEZ) is a signaling center that regulates patterned development of the upper jaw, and Sonic hedgehog (SHH) mediates FEZ activity. Induction of SHH expression in the FEZ results from SHH-dependent signals from the brain and neural crest cells. Given the role of miRNAs in modulating gene expression, we investigated the extent to which miRNAs regulate SHH expression and FEZ signaling.

RESULTS

In the FEZ, the miR-199 family appears to be regulated by SHH-dependent signals from the brain; expression of this family increased from HH18 to HH22, and upon activation of SHH signaling in the brain. However, the miR-199 family is more broadly expressed in the mesenchyme of the frontonasal process and adjacent neuroepithelium. Downregulating the miR-199 genes expanded SHH expression in the FEZ, resulting in wider faces, while upregulating miR-199 genes resulted in decreased SHH expression and narrow faces. Hypoxia inducible factor 1 alpha (HIF1A) and mitogen-activated protein kinase kinase kinase 4 (MAP3K4) appear to be potential targets of miR-199b. Reduction of MAP3K4 altered beak development but increased apoptosis, while reducing HIF1A reduced expression of SHH in the FEZ and produced malformations independent of apoptosis.

CONCLUSIONS

Our results demonstrate that this miRNA family appears to participate in regulating SHH expression in the FEZ; however, specific molecular mechanisms remain unknown.

Facial shape and allometry quantitative trait locus intervals in the Diversity Outbred mouse are enriched for known skeletal and facial development genes

The biology of how faces are built and come to differ from one another is complex. Discovering normal variants that contribute to differences in facial morphology is one key to untangling this complexity, with important implications for medicine and evolutionary biology. This study maps quantitative trait loci (QTL) for skeletal facial shape using Diversity Outbred (DO) mice. The DO is a randomly outcrossed population with high heterozygosity that captures the allelic diversity of eight inbred mouse lines from three subspecies. The study uses a sample of 1147 DO animals (the largest sample yet employed for a shape QTL study in mouse), each characterized by 22 three-dimensional landmarks, 56,885 autosomal and X-chromosome markers, and sex and age classifiers. We identified 37 facial shape QTL across 20 shape principal components (PCs) using a mixed effects regression that accounts for kinship among observations. The QTL include some previously identified intervals as well as new regions that expand the list of potential targets for future experimental study. Three QTL characterized shape associations with size (allometry). Median support interval size was 3.5 Mb. Narrowing additional analysis to QTL for the five largest magnitude shape PCs, we found significant overrepresentation of genes with known roles in growth, skeletal and facial development, and sensory organ development. For most intervals, one or more of these genes lies within 0.25 Mb of the QTL's peak. QTL effect sizes were small, with none explaining more than 0.5% of facial shape variation. Thus, our results are consistent with a model of facial diversity that is influenced by key genes in skeletal and facial development and, simultaneously, is highly polygenic.

Aberrant cell segregation in the craniofacial primordium and the emergence of facial dysmorphology in craniofrontonasal syndrome

Craniofrontonasal syndrome (CFNS) is a rare X-linked disorder characterized by craniofacial, skeletal, and neurological anomalies and is caused by mutations in EFNB1. Heterozygous females are more severely affected by CFNS than hemizygous males, a phenomenon called cellular interference that results from EPHRIN-B1 mosaicism. In Efnb1 heterozygous mice, mosaicism for EPHRIN-B1 results in cell sorting and more severe phenotypes than Efnb1 hemizygous males, but how craniofacial dysmorphology arises from cell segregation is unknown and CFNS etiology therefore remains poorly understood. Here, we couple geometric morphometric techniques with temporal and spatial interrogation of embryonic cell segregation in mouse mutant models to elucidate mechanisms underlying CFNS pathogenesis. By generating EPHRIN-B1 mosaicism at different developmental timepoints and in specific cell populations, we find that EPHRIN-B1 regulates cell segregation independently in early neural development and later in craniofacial development, correlating with the emergence of quantitative differences in face shape. Whereas specific craniofacial shape changes are qualitatively similar in Efnb1 heterozygous and hemizygous mutant embryos, heterozygous embryos are quantitatively more severely affected, indicating that Efnb1 mosaicism exacerbates loss of function phenotypes rather than having a neomorphic effect. Notably, neural tissue-specific disruption of Efnb1 does not appear to contribute to CFNS craniofacial dysmorphology, but its disruption within neural crest cell-derived mesenchyme results in phenotypes very similar to widespread loss. EPHRIN-B1 can bind and signal with EPHB1, EPHB2, and EPHB3 receptor tyrosine kinases, but the signaling partner(s) relevant to CFNS are unknown. Geometric morphometric analysis of an allelic series of Ephb1; Ephb2; Ephb3 mutant embryos indicates that EPHB2 and EPHB3 are key receptors mediating Efnb1 hemizygous-like phenotypes, but the complete loss of EPHB1-3 does not fully recapitulate the severity of CFNS-like Efnb1 heterozygosity. Finally, by generating Efnb1+/Δ; Ephb1; Ephb2; Ephb3 quadruple knockout mice, we determine how modulating cumulative receptor activity influences cell segregation in craniofacial development and find that while EPHB2 and EPHB3 play an important role in craniofacial cell segregation, EPHB1 is more important for cell segregation in the brain; surprisingly, complete loss of EPHB1-EPHB3 does not completely abrogate cell segregation. Together, these data advance our understanding of the etiology and signaling interactions underlying CFNS dysmorphology.

How the diversity of the faces arises

BACKGROUND

The evolution of the face is crucial for each species to adapt to different diets, environments, and in some species, to promote social interaction. The diversity in the shapes of the face results from divergence in the process of facial development that begins during early embryonic development.

HIGHLIGHTS

Here we review the recent advancements in the understanding of the genetic, epigenetic, molecular, and cellular basis of facial diversity. We also review the robustness of facial development and how it relates to the evolution of the face. Finally, we discuss the current challenges in achieving a deeper understanding of facial diversity.

CONCLUSION

We have gained much knowledge with respect to cis-regulatory elements, gene expression, cellular behavior, and the physical forces in facial development in the past two decades. Significant interdisciplinary work is needed to integrate these varied pieces of information into a complete picture of how the diversity of faces arises.

Nonlinear gene expression-phenotype relationships contribute to variation and clefting in the A/WySn mouse

BACKGROUND

Cleft lip and palate is one of the most common human birth defects, but the underlying etiology is poorly understood. The A/WySn mouse is a spontaneously occurring model of multigenic clefting in which 20% to 30% of individuals develop an orofacial cleft. Recent work has shown altered methylation at a specific retrotransposon insertion downstream of the Wnt9b locus in clefting animals, which results in decreased Wnt9b expression.

RESULTS

Using a newly developed protocol that allows us to measure morphology, gene expression, and DNA methylation in the same embryo, we relate gene expression in an individual embryo directly to its three-dimensional morphology for the first time. We find that methylation at the retrotransposon relates to Wnt9b expression and morphology. IAP methylation relates to shape of the nasal process in a manner consistent with clefting. Embryos with low IAP methylation exhibit increased among-individual variance in facial shape.

CONCLUSIONS

Methylation and gene expression relate nonlinearly to nasal process morphology. Individuals at one end of a continuum of phenotypic states display a clinical phenotype and increased phenotypic variation. Variable penetrance and expressivity in this model is likely determined both by among-individual variation in methylation and changes in phenotypic robustness along the underlying liability distribution for orofacial clefting.

Molecular and cellular mechanisms underlying the evolution of form and function in the amniote jaw

The amniote jaw complex is a remarkable amalgamation of derivatives from distinct embryonic cell lineages. During development, the cells in these lineages experience concerted movements, migrations, and signaling interactions that take them from their initial origins to their final destinations and imbue their derivatives with aspects of form including their axial orientation, anatomical identity, size, and shape. Perturbations along the way can produce defects and disease, but also generate the variation necessary for jaw evolution and adaptation. We focus on molecular and cellular mechanisms that regulate form in the amniote jaw complex, and that enable structural and functional integration. Special emphasis is placed on the role of cranial neural crest mesenchyme (NCM) during the species-specific patterning of bone, cartilage, tendon, muscle, and other jaw tissues. We also address the effects of biomechanical forces during jaw development and discuss ways in which certain molecular and cellular responses add adaptive and evolutionary plasticity to jaw morphology. Overall, we highlight how variation in molecular and cellular programs can promote the phenomenal diversity and functional morphology achieved during amniote jaw evolution or lead to the range of jaw defects and disease that affect the human condition.

Are endocasts good proxies for brain size and shape in archosaurs throughout ontogeny?

Cranial endocasts, or the internal molds of the braincase, are a crucial correlate for investigating the neuroanatomy of extinct vertebrates and tracking brain evolution through deep time. Nevertheless, the validity of such studies pivots on the reliability of endocasts as a proxy for brain morphology. Here, we employ micro-computed tomography imaging, including diffusible iodine-based contrast-enhanced CT, and a three-dimensional geometric morphometric framework to examine both size and shape differences between brains and endocasts of two exemplar archosaur taxa - the American alligator (Alligator mississippiensis) and the domestic chicken (Gallus gallus). With ontogenetic sampling, we quantitatively evaluate how endocasts differ from brains and whether this deviation changes during development. We find strong size and shape correlations between brains and endocasts, divergent ontogenetic trends in the brain-to-endocast correspondence between alligators and chickens, and a comparable magnitude between brain-endocast shape differences and intraspecific neuroanatomical variation. The results have important implications for paleoneurological studies in archosaurs. Notably, we demonstrate that the pattern of endocranial shape variation closely reflects brain shape variation. Therefore, analyses of endocranial morphology are unlikely to generate spurious conclusions about large-scale trends in brain size and shape. To mitigate any artifacts, however, paleoneurological studies should consider the lower brain-endocast correspondence in the hindbrain relative to the forebrain; higher size and shape correspondences in chickens than alligators throughout postnatal ontogeny; artificially 'pedomorphic' shape of endocasts relative to their corresponding brains; and potential biases in both size and shape data due to the lack of control for ontogenetic stages in endocranial sampling.

Developmental processes regulate craniofacial variation in disease and evolution

Variation in development mediates phenotypic differences observed in evolution and disease. Although the mechanisms underlying phenotypic variation are still largely unknown, recent research suggests that variation in developmental processes may play a key role. Developmental processes mediate genotype-phenotype relationships and consequently play an important role regulating phenotypes. In this review, we provide an example of how shared and interacting developmental processes may explain convergence of phenotypes in spliceosomopathies and ribosomopathies. These data also suggest a shared pathway to disease treatment. We then discuss three major mechanisms that contribute to variation in developmental processes: genetic background (gene-gene interactions), gene-environment interactions, and developmental stochasticity. Finally, we comment on evolutionary alterations to developmental processes, and the evolution of disease buffering mechanisms.

FoxO6 regulates Hippo signaling and growth of the craniofacial complex

The mechanisms that regulate post-natal growth of the craniofacial complex and that ultimately determine the size and shape of our faces are not well understood. Hippo signaling is a general mechanism to control tissue growth and organ size, and although it is known that Hippo signaling functions in neural crest specification and patterning during embryogenesis and before birth, its specific role in postnatal craniofacial growth remains elusive. We have identified the transcription factor FoxO6 as an activator of Hippo signaling regulating neonatal growth of the face. During late stages of mouse development, FoxO6 is expressed specifically in craniofacial tissues and FoxO6-/- mice undergo expansion of the face, frontal cortex, olfactory component and skull. Enlargement of the mandible and maxilla and lengthening of the incisors in FoxO6-/- mice are associated with increases in cell proliferation. In vitro and in vivo studies demonstrated that FoxO6 activates Lats1 expression, thereby increasing Yap phosphorylation and activation of Hippo signaling. FoxO6-/- mice have significantly reduced Hippo Signaling caused by a decrease in Lats1 expression and decreases in Shh and Runx2 expression, suggesting that Shh and Runx2 are also linked to Hippo signaling. In vitro, FoxO6 activates Hippo reporter constructs and regulates cell proliferation. Furthermore PITX2, a regulator of Hippo signaling is associated with Axenfeld-Rieger Syndrome causing a flattened midface and we show that PITX2 activates FoxO6 expression. Craniofacial specific expression of FoxO6 postnatally regulates Hippo signaling and cell proliferation. Together, these results identify a FoxO6-Hippo regulatory pathway that controls skull growth, odontogenesis and face morphology.

Quantitative Morphological Variation in the Developing Drosophila Wing

Quantitative genetic variation in morphology is pervasive in all species and is the basis for the evolution of differences among species. The measurement of morphological form in adults is now beginning to be combined with comparable measurements of form during development. Here we compare the shape of the developing wing to its adult form in a holometabolous insect, Drosophila melanogaster. We used protein expression patterns to measure shape in the developing precursors of the final adult wing. Three developmental stages were studied: late larval third instar, post-pupariation and in the adult fly. We studied wild-type animals in addition to mutants of two genes (shf and ds) that have known effects on adult wing shape and size. Despite experimental noise related to the difficulty of comparing developing structures, we found consistent differences in wing shape and size at each developmental stage between genotypes. Quantitative comparisons of variation arising at different developmental stages with the variation in the final structure enable us to determine when variation arises, and to generate hypotheses about the causes of that variation. In addition we provide linear rules allowing us to link wing morphology in the larva, with wing morphology in the pupa. Our approach provides a framework to analyze quantitative morphological variation in the developing fly wing. This framework should help to characterize the natural variation of the larval and pupal wing shape, and to measure the contribution of the processes occurring during these developmental stages to the natural variation in adult wing morphology.

Neural crest and the origin of species-specific pattern

For well over half of the 150 years since the discovery of the neural crest, the special ability of these cells to function as a source of species-specific pattern has been clearly recognized. Initially, this observation arose in association with chimeric transplant experiments among differentially pigmented amphibians, where the neural crest origin for melanocytes had been duly noted. Shortly thereafter, the role of cranial neural crest cells in transmitting species-specific information on size and shape to the pharyngeal arch skeleton as well as in regulating the timing of its differentiation became readily apparent. Since then, what has emerged is a deeper understanding of how the neural crest accomplishes such a presumably difficult mission, and this includes a more complete picture of the molecular and cellular programs whereby neural crest shapes the face of each species. This review covers studies on a broad range of vertebrates and describes neural-crest-mediated mechanisms that endow the craniofacial complex with species-specific pattern. A major focus is on experiments in quail and duck embryos that reveal a hierarchy of cell-autonomous and non-autonomous signaling interactions through which neural crest generates species-specific pattern in the craniofacial integument, skeleton, and musculature. By controlling size and shape throughout the development of these systems, the neural crest underlies the structural and functional integration of the craniofacial complex during evolution.

The developmental-genetics of canalization

Canalization, or robustness to genetic or environmental perturbations, is fundamental to complex organisms. While there is strong evidence for canalization as an evolved property that varies among genotypes, the developmental and genetic mechanisms that produce this phenomenon are very poorly understood. For evolutionary biology, understanding how canalization arises is important because, by modulating the phenotypic variation that arises in response to genetic differences, canalization is a determinant of evolvability. For genetics of disease in humans and for economically important traits in agriculture, this subject is important because canalization is a potentially significant cause of missing heritability that confounds genomic prediction of phenotypes. We review the major lines of thought on the developmental-genetic basis for canalization. These fall into two groups. One proposes specific evolved molecular mechanisms while the other deals with robustness or canalization as a more general feature of development. These explanations for canalization are not mutually exclusive and they overlap in several ways. General explanations for canalization are more likely to involve emergent features of development than specific molecular mechanisms. Disentangling these explanations is also complicated by differences in perspectives between genetics and developmental biology. Understanding canalization at a mechanistic level will require conceptual and methodological approaches that integrate quantitative genetics and developmental biology.

Developmental nonlinearity drives phenotypic robustness

Robustness to perturbation is a fundamental feature of complex organisms. Mutations are the raw material for evolution, yet robustness to their effects is required for species survival. The mechanisms that produce robustness are poorly understood. Nonlinearities are a ubiquitous feature of development that may link variation in development to phenotypic robustness. Here, we manipulate the gene dosage of a signaling molecule, Fgf8, a critical regulator of vertebrate development. We demonstrate that variation in Fgf8 expression has a nonlinear relationship to phenotypic variation, predicting levels of robustness among genotypes. Differences in robustness are not due to gene expression variance or dysregulation, but emerge from the nonlinearity of the genotype–phenotype curve. In this instance, embedded features of development explain robustness differences. How such features vary in natural populations and relate to genetic variation are key questions for unraveling the origin and evolvability of this feature of organismal development.

Developmental processes often involve nonlinearities, but the consequences for translating genotype to phenotype are not well characterized. Here, Green et al. vary Fgf8 signaling across allelic series of mice and show that phenotypic robustness in craniofacial shape is explained by a nonlinear effect of Fgf8 expression.

Unmodern Synthesis: Developmental Hierarchies and the Origin of Phenotypes

The question of whether the modern evolutionary synthesis requires an extension has recently become a topic of discussion, and a source of controversy. We suggest that this debate is, for the most part, not about the modern synthesis at all. Rather, it is about the extent to which genetic mechanisms can be regarded as the primary determinants of phenotypic characters. The modern synthesis has been associated with the idea that phenotypes are the result of gene products, while supporters of the extended synthesis have suggested that environmental factors, along with processes such as epigenetic inheritance, and niche construction play an important role in character formation. We argue that the methodology of the modern evolutionary synthesis has been enormously successful, but does not provide an accurate characterization of the origin of phenotypes. For its part, the extended synthesis has yet to be transformed into a testable theory, and accordingly, has yielded few results. We conclude by suggesting that the origin of phenotypes can only be understood by integrating findings from all levels of the organismal hierarchy. In most cases, parts and processes from a single level fail to accurately explain the presence of a given phenotypic trait.

Zebrafish zic2 controls formation of periocular neural crest and choroid fissure morphogenesis

The vertebrate retina develops in close proximity to the forebrain and neural crest-derived cartilages of the face and jaw. Coloboma, a congenital eye malformation, is associated with aberrant forebrain development (holoprosencephaly) and with craniofacial defects (frontonasal dysplasia) in humans, suggesting a critical role for cross-lineage interactions during retinal morphogenesis. ZIC2, a zinc-finger transcription factor, is linked to human holoprosencephaly. We have previously used morpholino assays to show zebrafish zic2 functions in the developing forebrain, retina and craniofacial cartilage. We now report that zebrafish with genetic lesions in zebrafish zic2 orthologs, zic2a and zic2b, develop with retinal coloboma and craniofacial anomalies. We demonstrate a requirement for zic2 in restricting pax2a expression and show evidence that zic2 function limits Hh signaling. RNA-seq transcriptome analysis identified an early requirement for zic2 in periocular neural crest as an activator of alx1, a transcription factor with essential roles in craniofacial and ocular morphogenesis in human and zebrafish. Collectively, these data establish zic2 mutant zebrafish as a powerful new genetic model for in-depth dissection of cell interactions and genetic controls during craniofacial complex development.

The Cyprinodon variegatus genome reveals gene expression changes underlying differences in skull morphology among closely related species

BACKGROUND

Understanding the genetic and developmental origins of phenotypic novelty is central to the study of biological diversity. In this study we identify modifications to the expression of genes at four developmental stages that may underlie jaw morphological differences among three closely related species of pupfish (genus Cyprinodon) from San Salvador Island, Bahamas. Pupfishes on San Salvador Island are trophically differentiated and include two endemic species that have evolved jaw morphologies unlike that of any other species in the genus Cyprinodon.

RESULTS

We find that gene expression differs significantly across recently diverged species of pupfish. Genes such as Bmp4 and calmodulin, previously implicated in jaw diversification in African cichlid fishes and Galapagos finches, were not found to be differentially expressed among species of pupfish. Instead we find multiple growth factors and cytokine/chemokine genes to be differentially expressed among these pupfish taxa. These include both genes and pathways known to affect craniofacial development, such as Wnt signaling, as well as novel genes and pathways not previously implicated in craniofacial development. These data highlight both shared and potentially unique sources of jaw diversity in pupfish and those identified in other evolutionary model systems such as Galapagos finches and African cichlids.

CONCLUSIONS

We identify modifications to the expression of genes involved in Wnt signaling, Igf signaling, and the inflammation response as promising avenues for future research. Our project provides insight into the magnitude of gene expression changes contributing to the evolution of morphological novelties, such as jaw structure, in recently diverged pupfish species.

Facial shape differences between rats selected for tame and aggressive behaviors

Domestication has been consistently accompanied by a suite of traits called the domestication syndrome. These include increased docility, changes in coat coloration, prolonged juvenile behaviors, modified function of adrenal glands and reduced craniofacial dimensions. Wilkins et al recently proposed that the mechanistic factor underlying traits that encompass the domestication syndrome was altered neural crest cell (NCC) development. NCC form the precursors to a large number of tissue types including pigment cells, adrenal glands, teeth and the bones of the face. The hypothesis that deficits in NCC development can account for the domestication syndrome was partly based on the outcomes of Dmitri Belyaev's domestication experiments initially conducted on silver foxes. After generations of selecting for tameness, the foxes displayed phenotypes observed in domesticated species. Belyaev also had a colony of rats selected over 64 generations for either tameness or defensive aggression towards humans. Here we focus on the facial morphology of Belyaev's tame, 'domesticated' rats to test whether: 1) tameness in rats causes craniofacial changes similar to those observed in the foxes; 2) facial shape, i.e. NCC-derived region, is distinct in the tame and aggressive rats. We used computed-tomography scans of rat skulls and landmark-based geometric morphometrics to quantify and analyze the facial skeleton. We found facial shape differences between the tame and aggressive rats that were independent of size and which mirrored changes seen in domesticated animals compared to their wild counterparts. However, there was no evidence of reduced sexual dimorphism in the face of the tame rats. This indicates that not all morphological changes in NCC-derived regions in the rats follow the pattern of shape change reported in domesticated animals or the silver foxes. Thus, certain phenotypic trends that are part of the domestication syndrome might not be consistently present in all experimental animal models.

Craniofacial Ciliopathies Reveal Specific Requirements for GLI Proteins during Development of the Facial Midline

Ciliopathies represent a broad class of disorders that affect multiple organ systems. The craniofacial complex is among those most severely affected when primary cilia are not functional. We previously reported that loss of primary cilia on cranial neural crest cells, via a conditional knockout of the intraflagellar transport protein KIF3a, resulted in midfacial widening due to a gain of Hedgehog (HH) activity. Here, we examine the molecular mechanism of how a loss of primary cilia can produce facial phenotypes associated with a gain of HH function. We show that loss of intraflagellar transport proteins (KIF3a or IFT88) caused aberrant GLI processing such that the amount of GLI3FL and GLI2FL was increased, thus skewing the ratio of GLIFL to GLIR in favor of the FL isoform. Genetic addition of GLI3R partially rescued the ciliopathic midfacial widening. Interestingly, despite several previous studies suggesting midfacial development relies heavily on GLI3R activity, the conditional loss of GLI3 alone did not reproduce the ciliopathic phenotype. Only the combined loss of both GLI2 and GLI3 was able to phenocopy the ciliopathic midfacial appearance. Our findings suggest that ciliopathic facial phenotypes are generated via loss of both GLI3R and GLI2R and that this pathology occurs via a de-repression mechanism. Furthermore, these studies suggest a novel role for GLI2R in craniofacial development.

Primary cilia are ubiquitous organelles that serve to transduce molecular signals within a cell. Loss of functional primary cilia results in a disease class called ciliopathies. Ciliopathies have a broad range of phenotypes; however, severe facial anomalies are commonly associated with this disease class. The facial midline is particularly sensitive to loss of primary cilia, frequently undergoing a significant widening. This phenotype is similar to that which occurs when there are gain-of-function defects in the Sonic Hedgehog pathway. This manuscript addresses the molecular basis for midfacial widening in ciliopathies. Importantly, we determine mechanisms to both rescue and phenocopy the ciliopathic midfacial phenotype. In sum, this work provides novel insight into the molecular mechanisms of midfacial patterning and the extent to which loss of cilia impact that process.

Genome-Wide Association Study Reveals Multiple Loci Influencing Normal Human Facial Morphology

Numerous lines of evidence point to a genetic basis for facial morphology in humans, yet little is known about how specific genetic variants relate to the phenotypic expression of many common facial features. We conducted genome-wide association meta-analyses of 20 quantitative facial measurements derived from the 3D surface images of 3118 healthy individuals of European ancestry belonging to two US cohorts. Analyses were performed on just under one million genotyped SNPs (Illumina OmniExpress+Exome v1.2 array) imputed to the 1000 Genomes reference panel (Phase 3). We observed genome-wide significant associations (p < 5 x 10⁻⁸) for cranial base width at 14q21.1 and 20q12, intercanthal width at 1p13.3 and Xq13.2, nasal width at 20p11.22, nasal ala length at 14q11.2, and upper facial depth at 11q22.1. Several genes in the associated regions are known to play roles in craniofacial development or in syndromes affecting the face: MAFB, PAX9, MIPOL1, ALX3, HDAC8, and PAX1. We also tested genotype-phenotype associations reported in two previous genome-wide studies and found evidence of replication for nasal ala length and SNPs in CACNA2D3 and PRDM16. These results provide further evidence that common variants in regions harboring genes of known craniofacial function contribute to normal variation in human facial features. Improved understanding of the genes associated with facial morphology in healthy individuals can provide insights into the pathways and mechanisms controlling normal and abnormal facial morphogenesis.

There is a great deal of evidence that genes influence facial appearance. This is perhaps most apparent when we look at our own families, since we are more likely to share facial features in common with our close relatives than with unrelated individuals. Nevertheless, little is known about how variation in specific regions of the genome relates to the kinds of distinguishing facial characteristics that give us our unique identities, e.g., the size and shape of our nose or how far apart our eyes are spaced. In this paper, we investigate this question by examining the association between genetic variants across the whole genome and a set of measurements designed to capture key aspects of facial form. We found evidence of genetic associations involving measures of eye, nose, and facial breadth. In several cases, implicated regions contained genes known to play roles in embryonic face formation or in syndromes in which the face is affected. Our ability to connect specific genetic variants to ubiquitous facial traits can inform our understanding of normal and abnormal craniofacial development, provide potential predictive models of evolutionary changes in human facial features, and improve our ability to create forensic facial reconstructions from DNA.

Hypertelorism and Orofacial Clefting Revisited: An Anthropometric Investigation

OBJECTIVE

Since the 1960s, multiple studies have reported a tendency toward hypertelorism in individuals with nonsyndromic orofacial clefts (OFCs). However, the association between specific cleft types and increased interorbital distance has been inconsistent. Using three-dimensional (3D) surface imaging, we tested whether different forms of clefting showed evidence of increased interorbital distance.

METHODS

Intercanthal and outercanthal distances and intercanthal indices were calculated from 3D facial surface images of 287 individuals with repaired OFCs. Raw measurements were converted to sex and age-normalized Z-scores. Mean Z-scores for individuals with cleft lip (CL), cleft lip and palate (CLP), and cleft palate (CP) were compared with reference normative values (controls) and one another directly using t tests and analysis of variance.

RESULTS

The CLP group showed a significant increase in intercanthal width (P = .001) and intercanthal index (P < .001) compared with reference norms. The CP group showed a significant decrease (P < .001) in outercanthal width. The CL group showed no difference from reference norms. The proportion of clinically hyperteloric individuals was generally low but highest in the CLP group (7.4%). Cleft severity had little effect on interorbital spacing.

CONCLUSIONS

Individuals with CLP exhibited on average a tendency toward mild hypertelorism, driven primarily by an increase in intercanthal distance. This tendency was not seen in CL or CP.

Cichlid fishes as a model to understand normal and clinical craniofacial variation

We have made great strides towards understanding the etiology of craniofacial disorders, especially for 'simple' Mendelian traits. However, the facial skeleton is a complex trait, and the full spectrum of genetic, developmental, and environmental factors that contribute to its final geometry remain unresolved. Forward genetic screens are constrained with respect to complex traits due to the types of genes and alleles commonly identified, developmental pleiotropy, and limited information about the impact of environmental interactions. Here, we discuss how studies in an evolutionary model - African cichlid fishes - can complement traditional approaches to understand the genetic and developmental origins of complex shape. Cichlids exhibit an unparalleled range of natural craniofacial morphologies that model normal human variation, and in certain instances mimic human facial dysmorphologies. Moreover, the evolutionary history and genomic architecture of cichlids make them an ideal system to identify the genetic basis of these phenotypes via quantitative trait loci (QTL) mapping and population genomics. Given the molecular conservation of developmental genes and pathways, insights from cichlids are applicable to human facial variation and disease. We review recent work in this system, which has identified lbh as a novel regulator of neural crest cell migration, determined the Wnt and Hedgehog pathways mediate species-specific bone morphologies, and examined how plastic responses to diet modulate adult facial shapes. These studies have not only revealed new roles for existing pathways in craniofacial development, but have identified new genes and mechanisms involved in shaping the craniofacial skeleton. In all, we suggest that combining work in traditional laboratory and evolutionary models offers significant potential to provide a more complete and comprehensive picture of the myriad factors that are involved in the development of complex traits.

The shapes of bird beaks are highly controlled by nondietary factors

Bird beaks are textbook examples of ecological adaptation to diet, but their shapes are also controlled by genetic and developmental histories. To test the effects of these factors on the avian craniofacial skeleton, we conducted morphometric analyses on raptors, a polyphyletic group at the base of the landbird radiation. Despite common perception, we find that the beak is not an independently targeted module for selection. Instead, the beak and skull are highly integrated structures strongly regulated by size, with axes of shape change linked to the actions of recently identified regulatory genes. Together, size and integration account for almost 80% of the shape variation seen between different species to the exclusion of morphological dietary adaptation. Instead, birds of prey use size as a mechanism to modify their feeding ecology. The extent to which shape variation is confined to a few major axes may provide an advantage in that it facilitates rapid morphological evolution via changes in body size, but may also make raptors especially vulnerable when selection pressures act against these axes. The phylogenetic position of raptors suggests that this constraint is prevalent in all landbirds and that breaking the developmental correspondence between beak and braincase may be the key novelty in classic passerine adaptive radiations.

The human brain and face: mechanisms of cranial, neurological and facial development revealed through malformations of holoprosencephaly, cyclopia and aberrations in chromosome 18

The study of inborn genetic errors can lend insight into mechanisms of normal human development and congenital malformations. Here, we present the first detailed comparison of cranial and neuro pathology in two exceedingly rare human individuals with cyclopia and alobar holoprosencephaly (HPE) in the presence and absence of aberrant chromosome 18 (aCh18). The aCh18 fetus contained one normal Ch18 and one with a pseudo-isodicentric duplication of chromosome 18q and partial deletion of 18p from 18p11.31 where the HPE gene, TGIF, resides, to the p terminus. In addition to synophthalmia, the aCh18 cyclopic malformations included a failure of induction of most of the telencephalon - closely approximating anencephaly, unchecked development of brain stem structures, near absence of the sphenoid bone and a malformed neurocranium and viscerocranium that constitute the median face. Although there was complete erasure of the olfactory and superior nasal structures, rudiments of nasal structures derived from the maxillary bone were evident, but with absent pharyngeal structures. The second non-aCh18 cyclopic fetus was initially classified as a true Cyclops, as it appeared to have a proboscis and one median eye with a single iris, but further analysis revealed two eye globes as expected for synophthalmic cyclopia. Furthermore, the proboscis was associated with the medial ethmoid ridge, consistent with an incomplete induction of these nasal structures, even as the nasal septum and paranasal sinuses were apparently developed. An important conclusion of this study is that it is the brain that predicts the overall configuration of the face, due to its influence on the development of surrounding skeletal structures. The present data using a combination of macroscopic, computed tomography (CT) and magnetic resonance imaging (MRI) techniques provide an unparalleled analysis on the extent of the effects of median defects, and insight into normal development and patterning of the brain, face and their skeletal support.

Hedgehog receptor function during craniofacial development

The Hedgehog signalling pathway plays a fundamental role in orchestrating normal craniofacial development in vertebrates. In particular, Sonic hedgehog (Shh) is produced in three key domains during the early formation of the head; neuroectoderm of the ventral forebrain, facial ectoderm and the pharyngeal endoderm; with signal transduction evident in both ectodermal and mesenchymal tissue compartments. Shh signalling from the prechordal plate and ventral midline of the diencephalon is required for appropriate division of the eyefield and forebrain, with mutation in a number of pathway components associated with Holoprosencephaly, a clinically heterogeneous developmental defect characterized by a failure of the early forebrain vesicle to divide into distinct halves. In addition, signalling from the pharyngeal endoderm and facial ectoderm plays an essential role during development of the face, influencing cranial neural crest cells that migrate into the early facial processes. In recent years, the complexity of Shh signalling has been highlighted by the identification of multiple novel proteins that are involved in regulating both the release and reception of this protein. Here, we review the contributions of Shh signalling during early craniofacial development, focusing on Hedgehog receptor function and describing the consequences of disruption for inherited anomalies of this region in both mouse models and human populations.

Genetics of murine craniofacial morphology: diallel analysis of the eight founders of the Collaborative Cross

Using eight inbred founder strains of the mouse Collaborative Cross (CC) project and their reciprocal F1 hybrids, we quantified variation in craniofacial morphology across mouse strains, explored genetic contributions to craniofacial variation that distinguish the founder strains, and tested whether specific or summary measures of craniofacial shape display stronger additive genetic contributions. This study thus provides critical information about phenotypic diversity among CC founder strains and about the genetic contributions to this phenotypic diversity, which is relevant to understanding the basis of variation in standard laboratory strains and natural populations. Craniofacial shape was quantified as a series of size-adjusted linear dimensions (RDs) and by principal components (PC) analysis of morphological landmarks captured from computed tomography images from 62 of the 64 reciprocal crosses of the CC founder strains. We first identified aspects of skull morphology that vary between these phenotypically 'normal' founder strains and that are defining characteristics of these strains. We estimated the contributions of additive and various non-additive genetic factors to phenotypic variation using diallel analyses of a subset of these strongly differing RDs and the first eight PCs of skull shape variation. We find little difference in the genetic contributions to RD measures and PC scores, suggesting fundamental similarities in the magnitude of genetic contributions to both specific and summary measures of craniofacial phenotypes. Our results indicate that there are stronger additive genetic effects associated with defining phenotypic characteristics of specific founder strains, suggesting these distinguishing measures are good candidates for use in genotype-phenotype association studies of CC mice. Our results add significantly to understanding of genotype-phenotype associations in the skull, which serve as a foundation for modeling the origins of medically and evolutionarily relevant variation.

Developmental mechanisms underlying variation in craniofacial disease and evolution

Craniofacial disease phenotypes exhibit significant variation in penetrance and severity. Although many genetic contributions to phenotypic variation have been identified, genotype-phenotype correlations remain imprecise. Recent work in evolutionary developmental biology has exposed intriguing developmental mechanisms that potentially explain incongruities in genotype-phenotype relationships. This review focuses on two observations from work in comparative and experimental animal model systems that highlight how development structures variation. First, multiple genetic inputs converge on relatively few developmental processes. Investigation of when and how variation in developmental processes occurs may therefore help predict potential genetic interactions and phenotypic outcomes. Second, genetic mutation is typically associated with an increase in phenotypic variance. Several models outlining developmental mechanisms underlying mutational increases in phenotypic variance are discussed using Satb2-mediated variation in jaw size as an example. These data highlight development as a critical mediator of genotype-phenotype correlations. Future research in evolutionary developmental biology focusing on tissue-level processes may help elucidate the "black box" between genotype and phenotype, potentially leading to novel treatment, earlier diagnoses, and better clinical consultations for individuals affected by craniofacial anomalies.

Chronic up-regulation of sonic hedgehog has little effect on postnatal craniofacial morphology of euploid and trisomic mice

BACKGROUND

In Ts65Dn, a mouse model of Down syndrome (DS), brain and craniofacial abnormalities that parallel those in people with DS are linked to an attenuated cellular response to sonic hedgehog (SHH) signaling. If a similarly reduced response to SHH occurs in all trisomic cells, then chronic up-regulation of the pathway might have a positive effect on development in trisomic mice, resulting in amelioration of the craniofacial anomalies.

RESULTS

We crossed Ts65Dn with Ptch1(tm1Mps/+) mice and quantified the craniofacial morphology of Ts65Dn;Ptch(+/-) offspring to assess whether a chronic up-regulation of the SHH pathway rescued DS-related anomalies. Ts65Dn;Ptch1(+/-) mice experience a chronic increase in SHH in SHH-receptive cells due to haploinsufficiency of the pathway suppressor, Ptch1. Chronic up-regulation had minimal effect on craniofacial shape and did not correct facial abnormalities in Ts65Dn;Ptch(+/-) mice. We further compared effects of this chronic up-regulation of SHH with acute pathway stimulation in mice treated on the day of birth with a SHH pathway agonist, SAG. We found that SHH affects facial morphology differently based on chronic vs. acute postnatal pathway up-regulation.

CONCLUSIONS

Our findings have implications for understanding the function of SHH in craniofacial development and for the potential use of SHH-based agonists to treat DS-related abnormalities.

Neural crest-mediated bone resorption is a determinant of species-specific jaw length

Precise control of jaw length during development is crucial for proper form and function. Previously we have shown that in birds, neural crest mesenchyme (NCM) confers species-specific size and shape to the beak by regulating molecular and histological programs for the induction and deposition of cartilage and bone. Here we reveal that a hitherto unrecognized but similarly essential mechanism for establishing jaw length is the ability of NCM to mediate bone resorption. Osteoclasts are considered the predominant cells that resorb bone, although osteocytes have also been shown to participate in this process. In adults, bone resorption is tightly coupled to bone deposition as a means to maintain skeletal homeostasis. Yet, the role and regulation of bone resorption during growth of the embryonic skeleton have remained relatively unexplored. We compare jaw development in short-beaked quail versus long-billed duck and find that quail have substantially higher levels of enzymes expressed by bone-resorbing cells including tartrate-resistant acid phosphatase (TRAP), Matrix metalloproteinase 13 (Mmp13), and Mmp9. Then, we transplant NCM destined to form the jaw skeleton from quail to duck and generate chimeras in which osteocytes arise from quail donor NCM and osteoclasts come exclusively from the duck host. Chimeras develop quail-like jaw skeletons coincident with dramatically elevated expression of TRAP, Mmp13, and Mmp9. To test for a link between bone resorption and jaw length, we block resorption using a bisphosphonate, osteoprotegerin protein, or an MMP13 inhibitor, and this significantly lengthens the jaw. Conversely, activating resorption with RANKL protein shortens the jaw. Finally, we find that higher resorption in quail presages their relatively lower adult jaw bone mineral density (BMD) and that BMD is also NCM-mediated. Thus, our experiments suggest that NCM not only controls bone resorption by its own derivatives but also modulates the activity of mesoderm-derived osteoclasts, and in so doing enlists bone resorption as a key patterning mechanism underlying the functional morphology and evolution of the jaw.