Craniofacial genetics: Where have we been and where are we going?

The level of protein in the maternal murine diet modulates the facial appearance of the offspring via mTORC1 signaling

The development of craniofacial skeletal structures is fascinatingly complex and elucidation of the underlying mechanisms will not only provide novel scientific insights, but also help develop more effective clinical approaches to the treatment and/or prevention of the numerous congenital craniofacial malformations. To this end, we performed a genome-wide analysis of RNA transcription from non-coding regulatory elements by CAGE-sequencing of the facial mesenchyme of human embryos and cross-checked the active enhancers thus identified against genes, identified by GWAS for the normal range human facial appearance. Among the identified active cis-enhancers, several belonged to the components of the PI3/AKT/mTORC1/autophagy pathway. To assess the functional role of this pathway, we manipulated it both genetically and pharmacologically in mice and zebrafish. These experiments revealed that mTORC1 signaling modulates craniofacial shaping at the stage of skeletal mesenchymal condensations, with subsequent fine-tuning during clonal intercalation. This ability of mTORC1 pathway to modulate facial shaping, along with its evolutionary conservation and ability to sense external stimuli, in particular dietary amino acids, indicate that the mTORC1 pathway may play a role in facial phenotypic plasticity. Indeed, the level of protein in the diet of pregnant female mice influenced the activity of mTORC1 in fetal craniofacial structures and altered the size of skeletogenic clones, thus exerting an impact on the local geometry and craniofacial shaping. Overall, our findings indicate that the mTORC1 signaling pathway is involved in the effect of environmental conditions on the shaping of craniofacial structures.

Jaw size variation is associated with a novel craniofacial function for galanin receptor 2 in an adaptive radiation of pupfishes

Understanding the genetic basis of novel adaptations in new species is a fundamental question in biology. Here we demonstrate a new role for galr2 in vertebrate craniofacial development using an adaptive radiation of trophic specialist pupfishes endemic to San Salvador Island, Bahamas. We confirmed the loss of a putative Sry transcription factor binding site upstream of galr2 in scale-eating pupfish and found significant spatial differences in galr2 expression among pupfish species in Meckel's cartilage using in situ hybridization chain reaction (HCR). We then experimentally demonstrated a novel role for Galr2 in craniofacial development by exposing embryos to Garl2-inhibiting drugs. Galr2-inhibition reduced Meckel's cartilage length and increased chondrocyte density in both trophic specialists but not in the generalist genetic background. We propose a mechanism for jaw elongation in scale-eaters based on the reduced expression of galr2 due to the loss of a putative Sry binding site. Fewer Galr2 receptors in the scale-eater Meckel's cartilage may result in their enlarged jaw lengths as adults by limiting opportunities for a circulating Galr2 agonist to bind to these receptors during development. Our findings illustrate the growing utility of linking candidate adaptive SNPs in non-model systems with highly divergent phenotypes to novel vertebrate gene functions.

Integrative analysis of transcriptome dynamics during human craniofacial development identifies candidate disease genes

Craniofacial disorders arise in early pregnancy and are one of the most common congenital defects. To fully understand how craniofacial disorders arise, it is essential to characterize gene expression during the patterning of the craniofacial region. To address this, we performed bulk and single-cell RNA-seq on human craniofacial tissue from 4-8 weeks post conception. Comparisons to dozens of other human tissues revealed 239 genes most strongly expressed during craniofacial development. Craniofacial-biased developmental enhancers were enriched +/- 400 kb surrounding these craniofacial-biased genes. Gene co-expression analysis revealed that regulatory hubs are enriched for known disease causing genes and are resistant to mutation in the normal healthy population. Combining transcriptomic and epigenomic data we identified 539 genes likely to contribute to craniofacial disorders. While most have not been previously implicated in craniofacial disorders, we demonstrate this set of genes has increased levels of de novo mutations in orofacial clefting patients warranting further study.

Jaw size variation is associated with a novel craniofacial function for galanin receptor 2 in an adaptive radiation of pupfishes

Understanding the genetic basis of novel adaptations in new species is a fundamental question in biology that also provides an opportunity to uncover new genes and regulatory networks with potential clinical relevance. Here we demonstrate a new role for galr2 in vertebrate craniofacial development using an adaptive radiation of trophic specialist pupfishes endemic to San Salvador Island in the Bahamas. We confirmed the loss of a putative Sry transcription factor binding site in the upstream region of galr2 in scale-eating pupfish and found significant spatial differences in galr2 expression among pupfish species in Meckel's cartilage and premaxilla using in situ hybridization chain reaction (HCR). We then experimentally demonstrated a novel function for Galr2 in craniofacial development and jaw elongation by exposing embryos to drugs that inhibit Galr2 activity. Galr2-inhibition reduced Meckel's cartilage length and increased chondrocyte density in both trophic specialists but not in the generalist genetic background. We propose a mechanism for jaw elongation in scale-eaters based on the reduced expression of galr2 due to the loss of a putative Sry binding site. Fewer Galr2 receptors in the scale-eater Meckel's cartilage may result in their enlarged jaw lengths as adults by limiting opportunities for a postulated Galr2 agonist to bind to these receptors during development. Our findings illustrate the growing utility of linking candidate adaptive SNPs in non-model systems with highly divergent phenotypes to novel vertebrate gene functions.

Surgical Strategy for the Treatment of Facial Clefts

Craniofacial clefts have an incidence of 1/700 [...]

Retrospective cumulative dietary risk assessment of craniofacial alterations by residues of pesticides

EFSA established cumulative assessment groups and conducted retrospective cumulative risk assessments for two types of craniofacial alterations (alterations due to abnormal skeletal development, head soft tissue alterations and brain neural tube defects) for 14 European populations of women in childbearing age. Cumulative acute exposure calculations were performed by probabilistic modelling using monitoring data collected by Member States in 2017, 2018 and 2019. A rigorous uncertainty analysis was performed using expert knowledge elicitation. Considering all sources of uncertainty, their dependencies and differences between populations, it was concluded with varying degrees of certainty that the MOET resulting from cumulative exposure is above 100 for the two types of craniofacial alterations. The threshold for regulatory consideration established by risk managers is therefore not exceeded. Considering the severity of the effects under consideration, it was also assessed whether the MOET is above 500. This was the case with varying levels of certainty for the head soft tissue alterations and brain neural tube defects. However, for the alterations due to abnormal skeletal development, it was found about as likely as not that the MOET is above 500 in most populations. For two populations, it was even found more likely that the MOET is below 500. These results were discussed in the light of the conservatism of the methodological approach.

The Medaka Inbred Kiyosu-Karlsruhe (MIKK) panel

Background

Unraveling the relationship between genetic variation and phenotypic traits remains a fundamental challenge in biology. Mapping variants underlying complex traits while controlling for confounding environmental factors is often problematic. To address this, we establish a vertebrate genetic resource specifically to allow for robust genotype-to-phenotype investigations. The teleost medaka (Oryzias latipes) is an established genetic model system with a long history of genetic research and a high tolerance to inbreeding from the wild.

Results

Here we present the Medaka Inbred Kiyosu-Karlsruhe (MIKK) panel: the first near-isogenic panel of 80 inbred lines in a vertebrate model derived from a wild founder population. Inbred lines provide fixed genomes that are a prerequisite for the replication of studies, studies which vary both the genetics and environment in a controlled manner, and functional testing. The MIKK panel will therefore enable phenotype-to-genotype association studies of complex genetic traits while allowing for careful control of interacting factors, with numerous applications in genetic research, human health, drug development, and fundamental biology.

Conclusions

Here we present a detailed characterization of the genetic variation across the MIKK panel, which provides a rich and unique genetic resource to the community by enabling large-scale experiments for mapping complex traits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02623-z.

Comparison of Presurgical Anthropometric Measures of Right and Left Complete Unilateral Cleft Lip and/or Palate

BACKGROUND

Left clefts occur twice as frequently as right ones. The sidedness has been suggested to influence certain outcomes. Some surgeons consider a right cleft more challenging to repair. This is often attributed to their reduced prevalence. The authors question whether this may be caused by morphologic differences. The authors' hypothesis is that there are anthropometric differences between left and right complete cleft lips.

METHODS

Patients with complete unilateral cleft lip, with or without cleft palate, operated on at the age of 3 to 6 months, between 2000 and 2018, by a single surgeon, were included. Eight standardized anthropometric measurements of the cleft lip, collected just before cleft lip repair, compare lip and vermillion dimensions and ratios between left and right clefts.

RESULTS

One hundred thirty-nine left and 80 right unilateral cleft lips were compared. A significant difference was found between left and right clefts for cleft-side to non-cleft-side ratios comparing the lateral lip element vertical heights and vermillion heights.

CONCLUSIONS

Patients with right cleft lips have a greater degree of lateral lip element hypoplasia, demonstrating greater deficiencies of lateral lip element vertical height and vermillion height when compared to patients with left clefts. This has clinical implications for preoperative assessment, choice of surgical technique, and postoperative and long-term outcomes.

Involvement of homeobox transcription factor Mohawk in palatogenesis

Palatogenesis is affected by many factors, including gene polymorphisms and exposure to toxic chemicals during sensitive developmental periods. Cleft palate is one of the most common congenital anomalies, and ongoing efforts to elucidate the molecular mechanisms underlying palatogenesis are providing useful insights to reduce the risk of this disorder. To identify novel potential regulators of palatogenesis, we analyzed public transcriptome datasets from a mouse model of cleft palate caused by selective deletion of transforming growth factor-β (TGFβ) receptor type 2 in cranial neural crest cells. We identified the homeobox transcription factor Mohawk (Mkx) as a gene downregulated in the maxilla of TGFβ knockout mice compared with wild-type mice. To examine the role of mkx in palatogenesis, we used CRISPR/Cas9 editing to generate zebrafish with impaired expression of mkxa and mkxb, the zebrafish homologs of Mkx. We found that mkx crispants expressed reduced levels of gli1, a critical transcription factor in the Sonic hedgehog (SHH) signaling pathway that plays an important role in the regulation of palatogenesis. Furthermore, we found that mkxa^-/- zebrafish were more susceptible than mkxa^+/+ zebrafish to the deleterious effects of cyclopamine, an inhibitor of SHH signaling, on upper jaw development. These results suggest that Mkx may be involved in palatogenesis regulated by TGFβ and SHH signaling, and that impairment in Mkx function may be related to the etiology of cleft palate.

Complex genetic architecture of three-dimensional craniofacial shape variation in domestic pigeons

Deciphering the genetic basis of vertebrate craniofacial variation is a longstanding biological problem with broad implications in evolution, development, and human pathology. One of the most stunning examples of craniofacial diversification is the adaptive radiation of birds, in which the beak serves essential roles in virtually every aspect of their life histories. The domestic pigeon (Columba livia) provides an exceptional opportunity to study the genetic underpinnings of craniofacial variation because of its unique balance of experimental accessibility and extraordinary phenotypic diversity within a single species. We used traditional and geometric morphometrics to quantify craniofacial variation in an F₂ laboratory cross derived from the straight-beaked Pomeranian Pouter and curved-beak Scandaroon pigeon breeds. Using a combination of genome-wide quantitative trait locus scans and multi-locus modeling, we identified a set of genetic loci associated with complex shape variation in the craniofacial skeleton, including beak shape, braincase shape, and mandible shape. Some of these loci control coordinated changes between different structures, while others explain variation in the size and shape of specific skull and jaw regions. We find that in domestic pigeons, a complex blend of both independent and coupled genetic effects underlie three-dimensional craniofacial morphology.

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.

Duration of the pubertal growth spurt in patients with increased craniofacial growth component in sagittal and vertical planes-retrospective and cross-sectional study

OBJECTIVES

The aim of the study is to assess the skeletal age at the onset and end of the pubertal growth spurt and determine its duration in four growth type groups: (1) normodivergent skeletal Class I (I N), (2) normodivergent skeletal Class III (III N), (3) high-angle skeletal Class III (III H) and (4) high-angle skeletal Class I (I H).

MATERIALS AND METHODS

Two hundred thirteen subjects were selected from 2163 examined files. The cervical vertebral maturation stage was recorded by means of Baccetti's method. The sagittal and vertical skeletal relations were evaluated according to Steiner analysis with Kaminek's modification. The duration of the pubertal growth spurt was calculated from the difference between the means of the chronological age related to CS3 and CS4 maturation stages.

RESULTS

The shortest lasting pubertal growth spurt was observed in group I N (1.1), followed by group III N (1.6). Major differences between arithmetic means CS4-CS3 were seen in groups I H and III H (2.3 and 2.7, respectively).

CONCLUSIONS

The following tendency was observed in the duration of the pubertal growth spurt: I N < III N < I H < III H. This tendency has statistical significance only in high-angle patients in comparison with normodivergent skeletal Class I.

CLINICAL RELEVANCE

Knowledge on the longer pubertal growth spurt in high-angle patients compared to patients with normal anteroposterior and vertical relationships can be useful in the selection of an appropriate therapeutic method and a treatment time.

Impact of low-frequency coding variants on human facial shape

The contribution of low-frequency variants to the genetic architecture of normal-range facial traits is unknown. We studied the influence of low-frequency coding variants (MAF < 1%) in 8091 genes on multi-dimensional facial shape phenotypes in a European cohort of 2329 healthy individuals. Using three-dimensional images, we partitioned the full face into 31 hierarchically arranged segments to model facial morphology at multiple levels, and generated multi-dimensional phenotypes representing the shape variation within each segment. We used MultiSKAT, a multivariate kernel regression approach to scan the exome for face-associated low-frequency variants in a gene-based manner. After accounting for multiple tests, seven genes (AR, CARS2, FTSJ1, HFE, LTB4R, TELO2, NECTIN1) were significantly associated with shape variation of the cheek, chin, nose and mouth areas. These genes displayed a wide range of phenotypic effects, with some impacting the full face and others affecting localized regions. The missense variant rs142863092 in NECTIN1 had a significant effect on chin morphology and was predicted bioinformatically to have a deleterious effect on protein function. Notably, NECTIN1 is an established craniofacial gene that underlies a human syndrome that includes a mandibular phenotype. We further showed that nectin1a mutations can affect zebrafish craniofacial development, with the size and shape of the mandibular cartilage altered in mutant animals. Findings from this study expanded our understanding of the genetic basis of normal-range facial shape by highlighting the role of low-frequency coding variants in several novel genes.

Insights into the genetic architecture of the human face

The human face is complex and multipartite, and characterization of its genetic architecture remains challenging. Using a multivariate genome-wide association study meta-analysis of 8,246 European individuals, we identified 203 genome-wide-significant signals (120 also study-wide significant) associated with normal-range facial variation. Follow-up analyses indicate that the regions surrounding these signals are enriched for enhancer activity in cranial neural crest cells and craniofacial tissues, several regions harbor multiple signals with associations to different facial phenotypes, and there is evidence for potential coordinated actions of variants. In summary, our analyses provide insights into the understanding of how complex morphological traits are shaped by both individual and coordinated genetic actions.

Chemical-Induced Cleft Palate Is Caused and Rescued by Pharmacological Modulation of the Canonical Wnt Signaling Pathway in a Zebrafish Model

Cleft palate is one of the most frequent birth defects worldwide. It causes severe problems regarding eating and speaking and requires long-term treatment. Effective prenatal treatment would contribute to reducing the risk of cleft palate. The canonical Wnt signaling pathway is critically involved in palatogenesis, and genetic or chemical disturbance of this signaling pathway leads to cleft palate. Presently, preventative treatment for cleft palate during prenatal development has limited efficacy, but we expect that zebrafish will provide a useful high-throughput chemical screening model for effective prevention. To achieve this, the zebrafish model should recapitulate cleft palate development and its rescue by chemical modulation of the Wnt pathway. Here, we provide proof of concept for a zebrafish chemical screening model. Zebrafish embryos were treated with 12 chemical reagents known to induce cleft palate in mammals, and all 12 chemicals induced cleft palate characterized by decreased proliferation and increased apoptosis of palatal cells. The cleft phenotype was enhanced by combinatorial treatment with Wnt inhibitor and teratogens. Furthermore, the expression of tcf7 and lef1 as a readout of the pathway was decreased. Conversely, cleft palate was prevented by Wnt agonist and the cellular defects were also prevented. In conclusion, we provide evidence that chemical-induced cleft palate is caused by inhibition of the canonical Wnt pathway. Our results indicate that this zebrafish model is promising for chemical screening for prevention of cleft palate as well as modulation of the Wnt pathway as a therapeutic target.

Innovative Molecular and Cellular Therapeutics in Cleft Palate Tissue Engineering

Clefts of the lip and/or palate are the most prevalent orofacial birth defects occurring in about 1:700 live human births worldwide. Early postnatal surgical interventions are extensive and staged to bring about optimal growth and fusion of palatal shelves. Severe cleft defects pose a challenge to correct with surgery alone, resulting in complications and sequelae requiring life-long, multidisciplinary care. Advances made in materials science innovation, including scaffold-based delivery systems for precision tissue engineering, now offer new avenues for stimulating bone formation at the site of surgical correction for palatal clefts. In this study, we review the present scientific literature on key developmental events that can go awry in palate development and the common surgical practices and challenges faced in correcting cleft defects. How key osteoinductive pathways implicated in palatogenesis inform the design and optimization of constructs for cleft palate correction is discussed within the context of translation to humans. Finally, we highlight new osteogenic agents and innovative delivery systems with the potential to be adopted in engineering-based therapeutic approaches for the correction of palatal defects. Impact statement Tissue-engineered scaffolds supplemented with osteogenic growth factors have attractive, largely unexplored possibilities to modulate molecular signaling networks relevant to driving palatogenesis in the context of congenital anomalies (e.g., cleft palate). Constructs that address this need may obviate current use of autologous bone grafts, thereby avoiding donor-site morbidity and other regenerative challenges in patients afflicted with palatal clefts. Combinations of biomaterials and drug delivery of diverse regenerative cues and biologics are currently transforming strategies exploited by engineers, scientists, and clinicians for palatal cleft repair.

Molecular Diagnostics and In Utero Therapeutics for Orofacial Clefts

Orofacial clefts and their management impose a substantial burden on patients, on their families, and on the health system. Under the current standard of care, affected patients are subjected to a lifelong journey of corrective surgeries and multidisciplinary management to replace bone and soft tissues, as well as restore esthetics and physiologic functions while restoring self-esteem and psychological health. Hence, a better understanding of the dynamic interplay of molecular signaling pathways at critical phases of palate development is necessary to pioneer novel prenatal interventions. Such pathways include transforming growth factor-β (Tgfβ), sonic hedgehog (Shh), wingless-integrated site (Wnt)/β-catenin, bone morphogenetic protein (Bmp), and fibroblast growth factor (Fgf) and its associated receptors, among others. Here, we summarize commonly used surgical methods used to correct cleft defects postnatally. We also review the advances made in prenatal diagnostics of clefts through imaging and genomics and the various in utero surgical corrections that have been attempted thus far. An overview of how key mediators of signaling that drive palatogenesis are emphasized in the context of the framework and rationale for the development and testing of therapeutics in animal model systems and in humans is provided. The pros and cons of in utero therapies that can potentially restore molecular homeostasis needed for the proper growth and fusion of palatal shelves are presented. The theme advanced throughout this review is the need to develop preclinical molecular therapies that could ultimately be translated into human trials that can correct orofacial clefts at earlier stages of development.

Chemical-induced craniofacial anomalies caused by disruption of neural crest cell development in a zebrafish model

Background

Craniofacial anomalies are among the most frequent birth defects worldwide, and are thought to be caused by gene‐environment interactions. Genetically manipulated zebrafish simulate human diseases and provide great advantages for investigating the etiology and pathology of craniofacial anomalies. Although substantial advances have been made in understanding genetic factors causing craniofacial disorders, limited information about the etiology by which environmental factors, such as teratogens, induce craniofacial anomalies is available in zebrafish.

Results

Zebrafish embryos displayed craniofacial malformations after teratogen treatments. Further observations revealed characteristic disruption of chondrocyte number, shape and stacking. These findings suggested aberrant development of cranial neural crest (CNC) cells, which was confirmed by gene expression analysis of the CNC. Notably, these observations suggested conserved etiological pathways between zebrafish and mammals including human. Furthermore, several of these chemicals caused malformations of the eyes, otic vesicle, and/or heart, representing a phenocopy of neurocristopathy, and these chemicals altered the expression levels of the responsible genes.

Conclusions

Our results demonstrate that chemical‐induced craniofacial malformation is caused by aberrant development of neural crest. This study indicates that zebrafish provide a platform for investigating contributions of environmental factors as causative agents of craniofacial anomalies and neurocristopathy.

What Do Animal Models Teach Us About Congenital Craniofacial Defects?

The formation of the head and face is a complex process which involves many different signaling cues regulating the migration, differentiation, and proliferation of the neural crest. This highly complex process is very error-prone, resulting in craniofacial defects in nearly 10,000 births in the United States annually. Due to the highly conserved mechanisms of craniofacial development, animal models are widely used to understand the pathogenesis of various human diseases and assist in the diagnosis and generation of preventative therapies and treatments. Here, we provide a brief background of craniofacial development and discuss several rare diseases affecting craniofacial bone development. We focus on rare congenital diseases of the cranial bone, facial jaw bones, and two classes of diseases, ciliopathies and RASopathies. Studying the animal models of these rare diseases sheds light not only on the etiology and pathology of each disease, but also provides meaningful insights towards the mechanisms which regulate normal development of the head and face.

A Cross-Sectional Study to Understand 3D Facial Differences in a Population of African Americans and Caucasians

Objective The purpose of this cross-sectional retrospective study was to use three-dimensional surface imaging to determine gender dimorphism and facial morphological changes from adolescence to adulthood in African American and Caucasian populations.

Materials and Methods Three-dimensional images were captured and the total sample size included 371 subjects. Images were combined using Rapidform 2006 Plus Pack 2 software to produce a male and female facial average for each population. Comparisons were conducted within the following categories: (1) gender comparison within each race, (2) adult and adolescent comparison within each race, and (3) adult and adolescent comparison between the races.

Results Adolescent gender comparisons within each race showed high percentages of similarity. However, adult females in both races showed more prominent periorbital, malar, and nasolabial regions and less prominent lower forehead, nose, and lower face compared with adult males of the same race. African American adult females showed increase in length and width of the face, increased nasal tip projection, and decreased periorbital regions compared with African American adolescent females. Welsh adult females had an increase in the nose and chin projection compared with Welsh adolescent females. Adult males of both races had increase in nose and chin projection, increase in length and width of the face, and decreased periorbital, malar, and nasolabial regions compared with adolescent males of the same race. African American adolescents had a wider alar base, more protrusive lips, and periorbital regions, and less prominent nose and chin compared with the Welsh adolescents. African American adults also had a wider alar base; more protrusive lips and periorbital regions; a broader face; and more retrusive chin, nose, nasolabial region; and lower forehead compared with Welsh adults.

Conclusions Few differences were noted between genders within the same racial groups during adolescence. However, changes became more distinct in adulthood. From adolescence to adulthood, facial morphologies were similarly matched within the gender for females; however, there were significant changes for males. Lastly, facial morphology patterns tend to be established early in life.

Vitamin D Receptor Signaling Regulates Craniofacial Cartilage Development in Zebrafish

Vitamin D plays essential roles in supporting the skeletal system. The active form of vitamin D functions through the vitamin D receptor (VDR). A hereditary vitamin-D-resistant rickets with facial dysmorphism has been reported, but the involvement of VDR signaling during early stages of craniofacial development remains to be elucidated. The present study investigated whether VDR signaling is implicated in zebrafish craniofacial cartilage development using a morpholino-based knockdown approach. Two paralogous VDR genes, vdra and vdrb, have been found in zebrafish embryos. Loss-of-vdra has no discernible effect on cartilage elements, whereas loss-of-vdrb causes reduction and malformation of craniofacial cartilages. Disrupting both vdra and vdrb leads to more severe defects or complete loss of cartilage. Notably, knockdown of vdrb results in elevated expression of follistatin a (fsta), a bone morphogenetic protein (BMP) antagonist, in the adjacent pharyngeal endoderm. Taken together, these findings strongly indicate that VDR signaling is required for early craniofacial cartilage development in zebrafish.

Mouse Skull Mean Shape and Shape Robustness Rely on Different Genetic Architectures and Different Loci

The genetic architecture of skull shape has been extensively studied in mice and the results suggest a highly polygenic and additive basis. In contrast few studies have explored the genetic basis of the skull variability. Canalization and developmental stability are the two components of phenotypic robustness. They have been proposed to be emergent properties of the genetic networks underlying the development of the trait itself, but this hypothesis has been rarely tested empirically. Here we use outbred mice to investigate the genetic architecture of canalization of the skull shape by implementing a genome-wide marginal epistatic test on 3D geometric morphometric data. The same data set had been used previously to explore the genetic architecture of the skull mean shape and its developmental stability. Here, we address two questions: (1) Are changes in mean shape and changes in shape variance associated with the same genomic regions? and (2) Do canalization and developmental stability rely on the same loci and genetic architecture and do they involve the same patterns of shape variation? We found that unlike skull mean shape, among-individual shape variance and fluctuating asymmetry (FA) show a total lack of additive effects. They are both associated with complex networks of epistatic interactions involving many genes (protein-coding and regulatory elements). Remarkably, none of the genomic loci affecting mean shape contribute these networks despite their enrichment for genes involved in craniofacial variation and diseases. We also found that the patterns of shape FA and individual variation are largely similar and rely on similar multilocus epistatic genetic networks, suggesting that the processes channeling variation within and among individuals are largely common. However, the loci involved in these two networks are completely different. This in turn underlines the difference in the origin of the variation at these two levels, and points at buffering processes that may be specific to each level.