The Interaction of Genetic Background and Mutational Effects in Regulation of Mouse Craniofacial Shape

New CRISPR/Cas9-based Fgfr2C361Y/+ mouse model of Crouzon syndrome exhibits skull and behavioral abnormalities

Crouzon syndrome (CS), a syndromic craniosynostosis, is a craniofacial developmental deformity caused by mutations in fibroblast growth factor receptor 2 (FGFR2). Previous CS mouse models constructed using traditional gene editing techniques faced issues such as low targeting efficiency, extended lineage cycles, and inconsistent and unstable phenotypes. In this study, a CRISPR/Cas9-mediated strategy was employed to induce a functional augmentation of the Fgfr2 point mutation in mice. Various techniques, including bone staining, micro-CT, histological methods, and behavioral experiments, were employed to systematically examine and corroborate phenotypic disparities between mutant mice (Fgfr2C361Y/+) and their wild-type littermates. Confirmed via PCR-Sanger sequencing, we successfully induced the p.Cys361Tyr missense mutation in the Fgfr2 IIIc isoform of the extracellular domain (corresponding to the p.Cys342Tyr mutation in humans) based on Fgfr2-215 transcript (ENSMUST00000122054.8). Fgfr2C361Y/+ mice exhibited characteristics consistent with the phenotypic features associated with CS, including skull-vault craniosynostosis, skull deformity, shallow orbits accompanied by exophthalmos, midface hypoplasia with malocclusion, and shortened skull base, notably without any apparent limb defects. Furthermore, mutant mice displayed behavioral abnormalities encompassing deficits in learning and memory, social interaction, and motor dysfunction, without anxiety-related disorders. Histopathological examination of the hippocampal region revealed structural abnormalities, suggesting possible brain development impairment secondary to craniosynostosis. In conclusion, we constructed a novel gene-edited Fgfr2C361Y/+ mice strain based on CRISPR/Cas9, which displayed skull and behavioral abnormalities, serving as a new model for studying genetic molecular mechanisms and exploring treatments for CS. KEY MESSAGES: CRISPR/Cas9 crafted a Crouzon model by enhancing Fgfr2-C361Y in mice. Fgfr2C361Y/+ mice replicate CS phenotypes-craniosynostosis and midface anomalies. Mutant mice show diverse behavioral abnormalities, impacting learning and memory. Fgfr2C361Y/+ mice offer a novel model for cranial suture studies and therapeutic exploration.

Foraging-induced craniofacial plasticity is associated with an early, robust and dynamic transcriptional response

Phenotypic plasticity is the ability of a single genotype to vary its phenotype in response to the environment. Plasticity of the skeletal system in response to mechanical input is widely studied, but the timing of its transcriptional regulation is not well understood. Here, we used the cichlid feeding apparatus to examine the transcriptional dynamics of skeletal plasticity over time. Using three closely related species that vary in their ability to remodel bone and a panel of 11 genes, including well-studied skeletal differentiation markers and newly characterized environmentally sensitive genes, we examined plasticity at one, two, four and eight weeks following the onset of alternate foraging challenges. We found that the plastic species exhibited environment-specific bursts in gene expression beginning at one week, followed by a sharp decline in levels, while the species with more limited plasticity exhibited consistently low levels of gene expression. This trend held across nearly all genes, suggesting that it is a hallmark of the larger plasticity regulatory network. We conclude that plasticity of the cichlid feeding apparatus is not the result of slowly accumulating gene expression difference over time, but rather is stimulated by early bursts of environment-specific gene expression followed by a return to homeostatic levels.

Craniosynostosis, inner ear, and renal anomalies in a child with complete loss of SPRY1 (sprouty homolog 1) function

INTRODUCTION

SPRY1 encodes protein sprouty homolog 1 (Spry-1), a negative regulator of receptor tyrosine kinase signalling. Null mutant mice display kidney/urinary tract abnormalities and altered size of the skull; complete loss-of-function of Spry-1 in humans has not been reported.

METHODS

Analysis of whole-genome sequencing data from individuals with craniosynostosis enrolled in the 100,000 Genomes Project identified a likely pathogenic variant within SPRY1. Reverse-transcriptase PCR and western blot analysis were used to investigate the effect of the variant on SPRY1 mRNA and protein, in lymphoblastoid cell lines from the patient and both parents.

RESULTS

A nonsense variant in SPRY1, encoding p.(Leu27*), was confirmed to be heterozygous in the unaffected parents and homozygous in the child. The child's phenotype, which included sagittal craniosynostosis, subcutaneous cystic lesions overlying the lambdoid sutures, hearing loss associated with bilateral cochlear and vestibular dysplasia and a unilateral renal cyst, overlapped the features reported in Spry1-/- null mice. Functional studies supported escape from nonsense-mediated decay, but western blot analysis demonstrated complete absence of full-length protein in the affected child and a marked reduction in both parents.

CONCLUSION

This is the first report of complete loss of Spry-1 function in humans, associated with abnormalities of the cranial sutures, inner ear, and kidneys.

HDAC9 structural variants disrupting TWIST1 transcriptional regulation lead to craniofacial and limb malformations

Structural variants (SVs) can affect protein-coding sequences as well as gene regulatory elements. However, SVs disrupting protein-coding sequences that also function as cis-regulatory elements remain largely uncharacterized. Here, we show that craniosynostosis patients with SVs containing the histone deacetylase 9 (HDAC9) protein-coding sequence are associated with disruption of TWIST1 regulatory elements that reside within the HDAC9 sequence. Based on SVs within the HDAC9-TWIST1 locus, we defined the 3'-HDAC9 sequence as a critical TWIST1 regulatory region, encompassing craniofacial TWIST1 enhancers and CTCF sites. Deletions of either Twist1 enhancers (eTw5-7Δ/Δ) or CTCF site (CTCF-5Δ/Δ) within the Hdac9 protein-coding sequence led to decreased Twist1 expression and altered anterior/posterior limb expression patterns of SHH pathway genes. This decreased Twist1 expression results in a smaller sized and asymmetric skull and polydactyly that resembles Twist1+/- mouse phenotype. Chromatin conformation analysis revealed that the Twist1 promoter interacts with Hdac9 sequences that encompass Twist1 enhancers and a CTCF site, and that interactions depended on the presence of both regulatory regions. Finally, a large inversion of the entire Hdac9 sequence (Hdac9 INV/+) in mice that does not disrupt Hdac9 expression but repositions Twist1 regulatory elements showed decreased Twist1 expression and led to a craniosynostosis-like phenotype and polydactyly. Thus, our study elucidates essential components of TWIST1 transcriptional machinery that reside within the HDAC9 sequence. It suggests that SVs encompassing protein-coding sequences could lead to a phenotype that is not attributed to its protein function but rather to a disruption of the transcriptional regulation of a nearby gene.

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 ).

Spatiotemporal Expression Patterns of Critical Genes Involved in FGF Signaling During Morphogenesis and Odontogenesis of Deciduous Molars in Miniature Pigs

The fibroblast growth factor (FGF) pathway plays an important role in epithelial-mesenchymal interactions during tooth development. Nevertheless, how the ligands, receptors, and antagonists of the FGF pathway are involved in epithelial-mesenchymal interactions remains largely unknown. Miniature pigs exhibit tooth anatomy and replacement patterns like those in humans and hence can serve as large animal models. The present study investigated the spatiotemporal expression patterns of critical genes encoding FGF ligands (FGF3, FGF4, FGF7, and FGF9), antagonists (SPRY2 and SPRY4) and receptors (FGFR1, FGFR2, and FGFR3) in the third deciduous molars of miniature pigs at the cap (embryonic day 40, E40), early bell (E50), and late bell (E60) stages. The results of in situ hybridization (ISH) with tyramide signal amplification and of qRT-PCR analysis revealed increased expression of FGF7, FGFR1, FGFR2, and SPRY4 in dental epithelium and of FGF7 and FGFR1 in mesenchyme from E40 to E50. In contrast, the results revealed decreased expression of FGF3, FGF4, FGF9, and FGFR3 in dental epithelium and of FGF4, FGF9, FGFR2, and FGFR3 in the mesenchyme from E40 to E60. Mesenchyme signals of FGF3, FGF4, FGF7, SPRY2, FGFR2, and FGFR3 were concentrated in the odontoblast layer from E50 to E60. The distinct expression patterns of these molecules indicated elaborate regulation during dental morphogenesis. Our results provide a foundation for further investigation into fine-tuning dental morphogenesis and odontogenesis by controlling interactions between dental epithelium and mesenchyme, thus promoting tooth regeneration in large mammals.

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.

Loss of Sprouty Produces a Ciliopathic Skeletal Phenotype in Mice Through Upregulation of Hedgehog Signaling

The Sprouty family is a highly conserved group of intracellular modulators of receptor tyrosine kinase (RTK)-signaling pathways, which have been recently linked to primary cilia. Disruptions in the structure and function of primary cilia cause inherited disorders called ciliopathies. We aimed to evaluate Sprouty2 and Sprouty4 gene-dependent alterations of ciliary structure and to focus on the determination of its association with Hedgehog signaling defects in chondrocytes. Analysis of the transgenic mice phenotype with Sprouty2 and Sprouty4 deficiency revealed several defects, including improper endochondral bone formation and digit patterning, or craniofacial and dental abnormalities. Moreover, reduced bone thickness and trabecular bone mass, skull deformities, or chondroma-like lesions were revealed. All these pathologies might be attributed to ciliopathies. Elongation of the ciliary axonemes in embryonic and postnatal growth plate chondrocytes was observed in Sprouty2^-/- and Sprouty2^+/- /Sprouty4^-/- mutants compared with corresponding littermate controls. Also, cilia-dependent Hedgehog signaling was upregulated in Sprouty2/4 mutant animals. Ptch1 and Ihh expression were upregulated in the autopodium and the proximal tibia of Sprouty2^-/- /Sprouty4^-/- mutants. Increased levels of the GLI3 repressor (GLI3R) form were detected in Sprouty2/4 mutant primary fibroblast embryonic cell cultures and tissues. These findings demonstrate that mouse lines deficient in Sprouty proteins manifest phenotypic features resembling ciliopathic phenotypes in multiple aspects and may serve as valuable models to study the association between overactivation of RTK and dysfunction of primary cilia during skeletogenesis. © 2021 American Society for Bone and Mineral Research (ASBMR).

The Mafb cleft-associated variant H131Q is not required for palatogenesis in the mouse

BACKGROUND

Orofacial clefts (OFCs) are common birth defects with complex etiology. Genome wide association studies for OFC have identified SNPs in and near MAFB. MAFB is a transcription factor critical for structural development of digits, kidneys, skin, and brain. MAFB is also expressed in the craniofacial region. Previous sequencing of MAFB in a Filipino population revealed a novel missense variant significantly associated with an increased risk for OFC. This MAFB variant, leading to the amino acid change H131Q, was knocked into the mouse Mafb, resulting in the MafbH131Q allele. The MafbH131Q construct was engineered to allow for deletion of Mafb ("Mafbdel ").

RESULTS

Mafbdel/del animals died shortly after birth. Conversely, MafbH131Q/H131Q mice survived into adulthood at Mendelian ratios. Mafbdel/del and MafbH131Q/H131Q heads exhibited normal macroscopic and histological appearance at all embryonic time points evaluated. The periderm was intact based on expression of keratin 6, p63, and E-cadherin. Despite no effect on craniofacial morphogenesis, H131Q inhibited the Mafb-dependent promoter activation of Arhgap29 in palatal mesenchymal, but not ectodermal-derived epithelial cells in a luciferase assay.

CONCLUSIONS

Mafb is dispensable for murine palatogenesis in vivo, and the cleft-associated variant H131Q, despite its lack of morphogenic effect, altered the expression of Arhgap29 in a cell-dependent context.

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.

Variation in phenotypes from a Bmp-Gata3 genetic pathway is modulated by Shh signaling

We sought to understand how perturbation of signaling pathways and their targets generates variable phenotypes. In humans, GATA3 associates with highly variable defects, such as HDR syndrome, microsomia and choanal atresia. We previously characterized a zebrafish point mutation in gata3 with highly variable craniofacial defects to the posterior palate. This variability could be due to residual Gata3 function, however, we observe the same phenotypic variability in gata3 null mutants. Using hsp:GATA3-GFP transgenics, we demonstrate that Gata3 function is required between 24 and 30 hpf. At this time maxillary neural crest cells fated to generate the palate express gata3. Transplantation experiments show that neural crest cells require Gata3 function for palatal development. Via a candidate approach, we determined if Bmp signaling was upstream of gata3 and if this pathway explained the mutant's phenotypic variation. Using BRE:d2EGFP transgenics, we demonstrate that maxillary neural crest cells are Bmp responsive by 24 hpf. We find that gata3 expression in maxillary neural crest requires Bmp signaling and that blocking Bmp signaling, in hsp:DN-Bmpr1a-GFP embryos, can phenocopy gata3 mutants. Palatal defects are rescued in hsp:DN-Bmpr1a-GFP;hsp:GATA3-GFP double transgenic embryos, collectively demonstrating that gata3 is downstream of Bmp signaling. However, Bmp attenuation does not alter phenotypic variability in gata3 loss-of-function embryos, implicating a different pathway. Due to phenotypes observed in hypomorphic shha mutants, the Sonic Hedgehog (Shh) pathway was a promising candidate for this pathway. Small molecule activators and inhibitors of the Shh pathway lessen and exacerbate, respectively, the phenotypic severity of gata3 mutants. Importantly, inhibition of Shh can cause gata3 haploinsufficiency, as observed in humans. We find that gata3 mutants in a less expressive genetic background have a compensatory upregulation of Shh signaling. These results demonstrate that the level of Shh signaling can modulate the phenotypes observed in gata3 mutants.

ADAMTS9 and ADAMTS20 are differentially affected by loss of B3GLCT in mouse model of Peters plus syndrome

Peters plus syndrome (MIM #261540 PTRPLS), characterized by defects in eye development, prominent forehead, hypertelorism, short stature and brachydactyly, is caused by mutations in the β3-glucosyltransferase (B3GLCT) gene. Protein O-fucosyltransferase 2 (POFUT2) and B3GLCT work sequentially to add an O-linked glucose β1-3fucose disaccharide to properly folded thrombospondin type 1 repeats (TSRs). Forty-nine proteins are predicted to be modified by POFUT2, and nearly half are members of the ADAMTS superfamily. Previous studies suggested that O-linked fucose is essential for folding and secretion of POFUT2-modified proteins and that B3GLCT-mediated extension to the disaccharide is essential for only a subset of targets. To test this hypothesis and gain insight into the origin of PTRPLS developmental defects, we developed and characterized two mouse B3glct knockout alleles. Using these models, we tested the role of B3GLCT in enabling function of ADAMTS9 and ADAMTS20, two highly conserved targets whose functions are well characterized in mouse development. The mouse B3glct mutants developed craniofacial and skeletal abnormalities comparable to PTRPLS. In addition, we observed highly penetrant hydrocephalus, white spotting and soft tissue syndactyly. We provide strong genetic and biochemical evidence that hydrocephalus and white spotting in B3glct mutants resulted from loss of ADAMTS20, eye abnormalities from partial reduction of ADAMTS9 and cleft palate from loss of ADAMTS20 and partially reduced ADAMTS9 function. Combined, these results provide compelling evidence that ADAMTS9 and ADAMTS20 were differentially sensitive to B3GLCT inactivation and suggest that the developmental defects in PTRPLS result from disruption of a subset of highly sensitive POFUT2/B3GLCT targets such as ADAMTS20.

Selective breeding modifies mef2ca mutant incomplete penetrance by tuning the opposing Notch pathway

Deleterious genetic mutations allow developmental biologists to understand how genes control development. However, not all loss of function genetic mutants develop phenotypic changes. Many deleterious mutations only produce a phenotype in a subset of mutant individuals, a phenomenon known as incomplete penetrance. Incomplete penetrance can confound analyses of gene function and our understanding of this widespread phenomenon remains inadequate. To better understand what controls penetrance, we capitalized on the zebrafish mef2ca mutant which produces craniofacial phenotypes with variable penetrance. Starting with a characterized mef2ca loss of function mutant allele, we used classical selective breeding methods to generate zebrafish strains in which mutant-associated phenotypes consistently appear with low or high penetrance. Strikingly, our selective breeding for low penetrance converted the mef2ca mutant allele behavior from homozygous lethal to homozygous viable. Meanwhile, selective breeding for high penetrance converted the mef2ca mutant allele from fully recessive to partially dominant. Comparing the selectively-bred low- and high-penetrance strains revealed that the strains initially respond similarly to the mutation, but then gene expression differences between strains emerge during development. Thus, altered temporal genetic circuitry can manifest through selective pressure to modify mutant penetrance. Specifically, we demonstrate differences in Notch signaling between strains, and further show that experimental manipulation of the Notch pathway phenocopies penetrance changes occurring through selective breeding. This study provides evidence that penetrance is inherited as a liability-threshold trait. Our finding that vertebrate animals can overcome a deleterious mutation by tuning genetic circuitry complements other reported mechanisms of overcoming deleterious mutations such as transcriptional adaptation of compensatory genes, alternative mRNA splicing, and maternal deposition of wild-type transcripts, which are not observed in our system. The selective breeding approach and the resultant genetic circuitry change we uncovered advances and expands our current understanding of genetic and developmental resilience.

Some deleterious gene mutations only affect a subset of genetically mutant animals. This widespread phenomenon, known as mutant incomplete penetrance, complicates discovery of causative gene mutations in both model organisms and human disease. This study utilized the zebrafish mef2ca transcription factor mutant that produces craniofacial skeleton defects with incomplete penetrance. Selectively breeding zebrafish families for low- or high-penetrance mutants for many generations created different zebrafish strains with consistently low or high penetrance. Comparing these strains allowed us to gain insight into the mechanisms that control penetrance. Specifically, genes under the control of mef2ca are initially similarly expressed between the two strains, but differences between strains emerge during development. We found that genetic manipulation of these downstream genes mimics the effects of our selective breeding. Thus, selective breeding for penetrance can change the genetic circuitry downstream of the mutated gene. We propose that small differences in gene circuitry between individuals is one mechanism underlying susceptibility or resilience to genetic mutations.

Gpr63 is a modifier of microcephaly in Ttc21b mouse mutants

The primary cilium is a signaling center critical for proper embryonic development. Previous studies have demonstrated that mice lacking Ttc21b have impaired retrograde trafficking within the cilium and multiple organogenesis phenotypes, including microcephaly. Interestingly, the severity of the microcephaly in Ttc21b^aln/aln homozygous null mutants is considerably affected by the genetic background and mutants on an FVB/NJ (FVB) background develop a forebrain significantly smaller than mutants on a C57BL/6J (B6) background. We performed a Quantitative Trait Locus (QTL) analysis to identify potential genetic modifiers and identified two regions linked to differential forebrain size: modifier of alien QTL1 (Moaq1) on chromosome 4 at 27.8 Mb and Moaq2 on chromosome 6 at 93.6 Mb. These QTLs were validated by constructing congenic strains. Further analysis of Moaq1 identified an orphan G-protein coupled receptor (GPCR), Gpr63, as a candidate gene. We identified a SNP that is polymorphic between the FVB and B6 strains in Gpr63 and creates a missense mutation predicted to be deleterious in the FVB protein. We used CRISPR-Cas9 genome editing to create two lines of FVB congenic mice: one with the B6 sequence of Gpr63 and the other with a deletion allele leading to a truncation of the GPR63 C-terminal tail. We then demonstrated that Gpr63 can localize to the cilium in vitro. These alleles affect ciliary localization of GPR63 in vitro and genetically interact with Ttc21b^aln/aln as Gpr63;Ttc21b double mutants show unique phenotypes including spina bifida aperta and earlier embryonic lethality. This validated Gpr63 as a modifier of multiple Ttc21b neural phenotypes and strongly supports Gpr63 as a causal gene (i.e., a quantitative trait gene, QTG) within the Moaq1 QTL.

TTC21B in humans is a known ciliopathy gene and contributes to the pathophysiology of a number of ciliopathies. Mice homozygous for a null allele of Ttc21b also have a spectrum of ciliopathy phenotypes, including microcephaly (small brain). Further work has shown that the severity of the microcephaly significantly depends on the genetic background of the mouse model. The genetic mechanisms of the Ttc21b pathophysiology and the interacting gene network remain far from understood. As an initial attempt to understand the underlying mechanism(s) underlying the variable effects on brain size, we performed a quantitative trait locus (QTL) analysis and found two regions of genomic significance that correlated with smaller brain size. We confirmed both QTLs with congenic lines. One of the two regions was small enough that we considered candidate genes and hypothesized Gpr63 might be a contributing locus for a number of reasons. We evaluated this hypothesis directly with precise variant creation using genome editing and provide evidence that Ttc21b and Gpr63 do indeed genetically interact. Thus, we have been able to combine classical QTL analysis and genome editing to directly test the resulting hypothesis.

Downregulation of FGF Signaling by Spry4 Overexpression Leads to Shape Impairment, Enamel Irregularities, and Delayed Signaling Center Formation in the Mouse Molar

FGF signaling plays a critical role in tooth development, and mutations in modulators of this pathway produce a number of striking phenotypes. However, many aspects of the role of the FGF pathway in regulating the morphological features and the mineral quality of the dentition remain unknown. Here, we used transgenic mice overexpressing the FGF negative feedback regulator Sprouty4 under the epithelial keratin 14 promoter (K14‐Spry4) to achieve downregulation of signaling in the epithelium. This led to highly penetrant defects affecting both cusp morphology and the enamel layer. We characterized the phenotype of erupted molars, identified a developmental delay in K14‐Spry4 transgenic embryos, and linked this with changes in the tooth developmental sequence. These data further delineate the role of FGF signaling in the development of the dentition and implicate the pathway in the regulation of tooth mineralization. © 2019 The Authors. JBMR Plus is published by Wiley Periodicals, Inc. on behalf of American Society for Bone and Mineral Research.

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.

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.

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.

Genetic backgrounds and hidden trait complexity in natural populations

Dissecting the genetic basis of natural phenotypic variation is a major goal in biology. We know that most traits are strongly heritable. However, their genetic architecture is a long-standing question, which is unfortunately confounded by the lack of complete knowledge of the genetic components as well as their phenotypic effect in a specific genetic background. Many genetic variants are known to affect phenotypes but the same functional variant can have a different effect on the phenotype in different individuals of the same species. Understanding the impact of genetic background on the expressivity of a given phenotype is essential because this effect complicates our ability to predict phenotype from genotype. Here, we briefly review recent progress on the exploration of the effect of genetic background and we discuss how a deeper characterization of the inheritance, expressivity and genetic interactions hidden behind the phenotypic landscape of natural variation could provide a better understanding of the relationship between genotype and phenotype.

Feedback regulation of RTK signaling in development

Precise regulation of the amplitude and duration of receptor tyrosine kinase (RTK) signaling is critical for the execution of cellular programs and behaviors. Understanding these control mechanisms has important implications for the field of developmental biology, and in recent years, the question of how augmentation or attenuation of RTK signaling via feedback loops modulates development has become of increasing interest. RTK feedback regulation is also important for human disease research; for example, germline mutations in genes that encode RTK signaling pathway components cause numerous human congenital syndromes, and somatic alterations contribute to the pathogenesis of diseases such as cancers. In this review, we survey regulators of RTK signaling that tune receptor activity and intracellular transduction cascades, with a focus on the roles of these genes in the developing embryo. We detail the diverse inhibitory mechanisms utilized by negative feedback regulators that, when lost or perturbed, lead to aberrant increases in RTK signaling. We also discuss recent biochemical and genetic insights into positive regulators of RTK signaling and how these proteins function in tandem with negative regulators to guide embryonic development.