Single Amino Acid Change Underlies Distinct Roles of H2A.Z Subtypes in Human Syndrome.

The developmental disorder Floating-Harbor syndrome (FHS) is caused by heterozygous truncating mutations in SRCAP, a gene encoding a chromatin remodeler mediating incorporation of histone variant H2A.Z. Here, we demonstrate that FHS-associated mutations result in loss of SRCAP nuclear localization, alter neural crest gene programs in human in vitro models and Xenopus embryos, and cause craniofacial defects. These defects are mediated by one of two H2A.Z subtypes, H2A.Z.2, whose knockdown mimics and whose overexpression rescues the FHS phenotype. Selective rescue by H2A.Z.2 is conferred by one of the three amino acid differences between the H2A.Z subtypes, S38/T38. We further show that H2A.Z.1 and H2A.Z.2 genomic occupancy patterns are qualitatively similar, but quantitatively distinct, and H2A.Z.2 incorporation at AT-rich enhancers and expression of their associated genes are both sensitized to SRCAP truncations. Altogether, our results illuminate the mechanism underlying a human syndrome and uncover selective functions of H2A.Z subtypes during development.

Integrated Control of Predatory Hunting by the Central Nucleus of the Amygdala.

Superior predatory skills led to the evolutionary triumph of jawed vertebrates. However, the mechanisms by which the vertebrate brain controls predation remain largely unknown. Here, we reveal a critical role for the central nucleus of the amygdala in predatory hunting. Both optogenetic and chemogenetic stimulation of central amygdala of mice elicited predatory-like attacks upon both insect and artificial prey. Coordinated control of cervical and mandibular musculatures, which is necessary for accurately positioning lethal bites on prey, was mediated by a central amygdala projection to the reticular formation in the brainstem. In contrast, prey pursuit was mediated by projections to the midbrain periaqueductal gray matter. Targeted lesions to these two pathways separately disrupted biting attacks upon prey versus the initiation of prey pursuit. Our findings delineate a neural network that integrates distinct behavioral modules and suggest that central amygdala neurons instruct predatory hunting across jawed vertebrates.

FGF signaling regulates development by processes beyond canonical pathways.

FGFs are key developmental regulators that engage a signal transduction cascade through receptor tyrosine kinases, prominently engaging ERK1/2 but also other pathways. However, it remains unknown whether all FGF activities depend on this canonical signal transduction cascade. To address this question, we generated allelic series of knock-in Fgfr1 and Fgfr2 mouse strains, carrying point mutations that disrupt binding of signaling effectors, and a kinase dead allele of Fgfr2 that broadly phenocopies the null mutant. When interrogated in cranial neural crest cells, we identified discrete functions for signaling pathways in specific craniofacial contexts, but point mutations, even when combined, failed to recapitulate the single or double null mutant phenotypes. Furthermore, the signaling mutations abrogated established FGF-induced signal transduction pathways, yet FGF functions such as cell-matrix and cell-cell adhesion remained unaffected, though these activities did require FGFR kinase activity. Our studies establish combinatorial roles of Fgfr1 and Fgfr2 in development and uncouple novel FGFR kinase-dependent cell adhesion properties from canonical intracellular signaling.

Mesenchymal Wnts are required for morphogenetic movements of calvarial osteoblasts during apical expansion.

Apical expansion of calvarial osteoblast progenitors from the cranial mesenchyme (CM) above the eye is integral to calvarial growth and enclosure of the brain. The cellular behaviors and signals underlying the morphogenetic process of calvarial expansion are unknown. Time-lapse light-sheet imaging of mouse embryos revealed calvarial progenitors intercalate in 3D in the CM above the eye, and exhibit protrusive and crawling activity more apically. CM cells express non-canonical Wnt/planar cell polarity (PCP) core components and calvarial osteoblasts are bidirectionally polarized. We found non-canonical ligand Wnt5a-/- mutants have less dynamic cell rearrangements and protrusive activity. Loss of CM-restricted Wntless (CM-Wls), a gene required for secretion of all Wnt ligands, led to diminished apical expansion of Osx+ calvarial osteoblasts in the frontal bone primordia in a non-cell autonomous manner without perturbing proliferation or survival. Calvarial osteoblast polarization, progressive cell elongation and enrichment for actin along the baso-apical axis were dependent on CM-Wnts. Thus, CM-Wnts regulate cellular behaviors during calvarial morphogenesis for efficient apical expansion of calvarial osteoblasts. These findings also offer potential insights into the etiologies of calvarial dysplasias.

Zebrafish raptor mutation inhibits the activity of mTORC1, inducing craniofacial defects due to autophagy-induced neural crest cell death.

The mechanistic target of rapamycin (mTOR) coordinates metabolism and cell growth with environmental inputs. mTOR forms two functional complexes: mTORC1 and mTORC2. Proper development requires both complexes but mTORC1 has unique roles in numerous cellular processes, including cell growth, survival and autophagy. Here, we investigate the function of mTORC1 in craniofacial development. We created a zebrafish raptor mutant via CRISPR/Cas9, to specifically disrupt mTORC1. The entire craniofacial skeleton and eyes were reduced in size in mutants; however, overall body length and developmental timing were not affected. The craniofacial phenotype associates with decreased chondrocyte size and increased neural crest cell death. We found that autophagy is elevated in raptor mutants. Chemical inhibition of autophagy reduced cell death and improved craniofacial phenotypes in raptor mutants. Genetic inhibition of autophagy, via mutation of the autophagy gene atg7, improved facial phenotypes in atg7;raptor double mutants, relative to raptor single mutants. We conclude that finely regulated levels of autophagy, via mTORC1, are crucial for craniofacial development.

BRD4 binds to active cranial neural crest enhancers to regulate RUNX2 activity during osteoblast differentiation.

Cornelia de Lange syndrome (CdLS) is a congenital disorder featuring facial dysmorphism, postnatal growth deficits, cognitive disability and upper limb abnormalities. CdLS is genetically heterogeneous, with cases arising from mutation of BRD4, a bromodomain protein that binds and reads acetylated histones. In this study, we have modeled CdLS facial pathology through mouse neural crest cell (NCC)-specific mutation of BRD4 to characterize cellular and molecular function in craniofacial development. Mice with BRD4 NCC loss of function died at birth with severe facial hypoplasia, cleft palate, mid-facial clefting and exencephaly. Following migration, BRD4 mutant NCCs initiated RUNX2 expression for differentiation to osteoblast lineages but failed to induce downstream RUNX2 targets required for lineage commitment. BRD4 bound to active enhancers to regulate expression of osteogenic transcription factors and extracellular matrix components integral for bone formation. RUNX2 physically interacts with a C-terminal domain in the long isoform of BRD4 and can co-occupy osteogenic enhancers. This BRD4 association is required for RUNX2 recruitment and appropriate osteoblast differentiation. We conclude that BRD4 controls facial bone development through osteoblast enhancer regulation of the RUNX2 transcriptional program.

Grainyhead-like 2 interacts with noggin to regulate tissue fusion in mouse.

Defective tissue fusion during mammalian embryogenesis results in congenital anomalies, such as exencephaly, spina bifida and cleft lip and/or palate. The highly conserved transcription factor grainyhead-like 2 (Grhl2) is a crucial regulator of tissue fusion, with mouse models lacking GRHL2 function presenting with a fully penetrant open cranial neural tube, facial and abdominal clefting (abdominoschisis), and an open posterior neuropore. Here, we show that GRHL2 interacts with the soluble morphogen protein and bone morphogenetic protein (BMP) inhibitor noggin (NOG) to impact tissue fusion during development. The maxillary prominence epithelium in embryos lacking Grhl2 shows substantial morphological abnormalities and significant upregulation of NOG expression, together with aberrantly distributed pSMAD5-positive cells within the neural crest cell-derived maxillary prominence mesenchyme, indicative of disrupted BMP signalling. Reducing this elevated NOG expression (by generating Grhl2-/-;Nog+/- embryos) results in delayed embryonic lethality, partial tissue fusion rescue, and restoration of tissue form within the craniofacial epithelia. These data suggest that aberrant epithelial maintenance, partially regulated by noggin-mediated regulation of BMP-SMAD pathways, may underpin tissue fusion defects in Grhl2-/- mice.

Etiology of craniofacial and cardiac malformations in a mouse model of SF3B4-related syndromes.

Pathogenic variants in SF3B4, a component of the U2 snRNP complex important for branchpoint sequence recognition and splicing, are responsible for the acrofacial disorders Nager and Rodriguez Syndrome, also known as SF3B4-related syndromes. Patients exhibit malformations in the head, face, limbs, vertebrae as well as the heart. To uncover the etiology of craniofacial malformations found in SF3B4-related syndromes, mutant mouse lines with homozygous deletion of Sf3b4 in neural crest cells (NCC) were generated. Like in human patients, these embryos had craniofacial and cardiac malformations with variable expressivity and penetrance. The severity and survival of Sf3b4 NCC mutants was modified by the level of Sf3b4 in neighboring non-NCC. RNA sequencing analysis of heads of embryos prior to morphological abnormalities revealed significant changes in expression of genes forming the NCC regulatory network, as well as an increase in exon skipping. Additionally, several key histone modifiers involved in craniofacial and cardiac development showed increased exon skipping. Increased exon skipping was also associated with use of a more proximal branch point, as well as an enrichment in thymidine bases in the 50 bp around the branch points. We propose that decrease in Sf3b4 causes changes in the expression and splicing of transcripts required for proper craniofacial and cardiac development, leading to abnormalities.

De novo variants implicate chromatin modification, transcriptional regulation, and retinoic acid signaling in syndromic craniosynostosis.

Craniosynostosis (CS) is the most common congenital cranial anomaly. Several Mendelian forms of syndromic CS are well described, but a genetic etiology remains elusive in a substantial fraction of probands. Analysis of exome sequence data from 526 proband-parent trios with syndromic CS identified a marked excess (observed 98, expected 33, p = 4.83 × 10-20) of damaging de novo variants (DNVs) in genes highly intolerant to loss-of-function variation (probability of LoF intolerance > 0.9). 30 probands harbored damaging DNVs in 21 genes that were not previously implicated in CS but are involved in chromatin modification and remodeling (4.7-fold enrichment, p = 1.1 × 10-11). 17 genes had multiple damaging DNVs, and 13 genes (CDK13, NFIX, ADNP, KMT5B, SON, ARID1B, CASK, CHD7, MED13L, PSMD12, POLR2A, CHD3, and SETBP1) surpassed thresholds for genome-wide significance. A recurrent gain-of-function DNV in the retinoic acid receptor alpha (RARA; c.865G>A [p.Gly289Arg]) was identified in two probands with similar CS phenotypes. CS risk genes overlap with those identified for autism and other neurodevelopmental disorders, are highly expressed in cranial neural crest cells, and converge in networks that regulate chromatin modification, gene transcription, and osteoblast differentiation. Our results identify several CS loci and have major implications for genetic testing and counseling.

Genome-wide analysis of copy-number variation in humans with cleft lip and/or cleft palate identifies COBLL1, RIC1, and ARHGEF38 as clefting genes.

Cleft lip with or without cleft palate (CL/P) is a common birth defect with a complex, heterogeneous etiology. It is well established that common and rare sequence variants contribute to the formation of CL/P, but the contribution of copy-number variants (CNVs) to cleft formation remains relatively understudied. To fill this knowledge gap, we conducted a large-scale comparative analysis of genome-wide CNV profiles of 869 individuals from the Philippines and 233 individuals of European ancestry with CL/P with three primary goals: first, to evaluate whether differences in CNV number, amount of genomic content, or amount of coding genomic content existed within clefting subtypes; second, to assess whether CNVs in our cohort overlapped with known Mendelian clefting loci; and third, to identify unestablished Mendelian clefting genes. Significant differences in CNVs across cleft types or in individuals with non-syndromic versus syndromic clefts were not observed; however, several CNVs in our cohort overlapped with known syndromic and non-syndromic Mendelian clefting loci. Moreover, employing a filtering strategy relying on population genetics data that rare variants are on the whole more deleterious than common variants, we identify several CNV-associated gene losses likely driving non-syndromic clefting phenotypes. By prioritizing genes deleted at a rare frequency across multiple individuals with clefts yet enriched in our cohort of individuals with clefts compared to control subjects, we identify COBLL1, RIC1, and ARHGEF38 as clefting genes. CRISPR-Cas9 mutagenesis of these genes in Xenopus laevis and Danio rerio yielded craniofacial dysmorphologies, including clefts analogous to those seen in human clefting disorders.

POLR1A variants underlie phenotypic heterogeneity in craniofacial, neural, and cardiac anomalies.

Heterozygous pathogenic variants in POLR1A, which encodes the largest subunit of RNA Polymerase I, were previously identified as the cause of acrofacial dysostosis, Cincinnati-type. The predominant phenotypes observed in the cohort of 3 individuals were craniofacial anomalies reminiscent of Treacher Collins syndrome. We subsequently identified 17 additional individuals with 12 unique heterozygous variants in POLR1A and observed numerous additional phenotypes including neurodevelopmental abnormalities and structural cardiac defects, in combination with highly prevalent craniofacial anomalies and variable limb defects. To understand the pathogenesis of this pleiotropy, we modeled an allelic series of POLR1A variants in vitro and in vivo. In vitro assessments demonstrate variable effects of individual pathogenic variants on ribosomal RNA synthesis and nucleolar morphology, which supports the possibility of variant-specific phenotypic effects in affected individuals. To further explore variant-specific effects in vivo, we used CRISPR-Cas9 gene editing to recapitulate two human variants in mice. Additionally, spatiotemporal requirements for Polr1a in developmental lineages contributing to congenital anomalies in affected individuals were examined via conditional mutagenesis in neural crest cells (face and heart), the second heart field (cardiac outflow tract and right ventricle), and forebrain precursors in mice. Consistent with its ubiquitous role in the essential function of ribosome biogenesis, we observed that loss of Polr1a in any of these lineages causes cell-autonomous apoptosis resulting in embryonic malformations. Altogether, our work greatly expands the phenotype of human POLR1A-related disorders and demonstrates variant-specific effects that provide insights into the underlying pathogenesis of ribosomopathies.

Successful therapeutic intervention in new mouse models of frizzled 2-associated congenital malformations.

Frizzled 2 (FZD2) is a transmembrane Wnt receptor. We previously identified a pathogenic human FZD2 variant in individuals with FZD2-associated autosomal dominant Robinow syndrome. The variant encoded a protein with a premature stop and loss of 17 amino acids, including a region of the consensus dishevelled-binding sequence. To model this variant, we used zygote microinjection and i-GONAD-based CRISPR/Cas9-mediated genome editing to generate a mouse allelic series. Embryos mosaic for humanized Fzd2W553* knock-in exhibited cleft palate and shortened limbs, consistent with patient phenotypes. We also generated two germline mouse alleles with small deletions: Fzd2D3 and Fzd2D4. Homozygotes for each allele exhibit a highly penetrant cleft palate phenotype, shortened limbs compared with wild type and perinatal lethality. Fzd2D4 craniofacial tissues indicated decreased canonical Wnt signaling. In utero treatment with IIIC3a (a DKK inhibitor) normalized the limb lengths in Fzd2D4 homozygotes. The in vivo replication represents an approach for further investigating the mechanism of FZD2 phenotypes and demonstrates the utility of CRISPR knock-in mice as a tool for investigating the pathogenicity of human genetic variants. We also present evidence for a potential therapeutic intervention.

Identification of a heterogeneous and dynamic ciliome during embryonic development and cell differentiation.

Primary cilia are nearly ubiquitous organelles that transduce molecular and mechanical signals. Although the basic structure of the cilium and the cadre of genes that contribute to ciliary formation and function (the ciliome) are believed to be evolutionarily conserved, the presentation of ciliopathies with narrow, tissue-specific phenotypes and distinct molecular readouts suggests that an unappreciated heterogeneity exists within this organelle. Here, we provide a searchable transcriptomic resource for a curated primary ciliome, detailing various subgroups of differentially expressed genes within the ciliome that display tissue and temporal specificity. Genes within the differentially expressed ciliome exhibited a lower level of functional constraint across species, suggesting organism and cell-specific function adaptation. The biological relevance of ciliary heterogeneity was functionally validated by using Cas9 gene-editing to disrupt ciliary genes that displayed dynamic gene expression profiles during osteogenic differentiation of multipotent neural crest cells. Collectively, this novel primary cilia-focused resource will allow researchers to explore longstanding questions related to how tissue and cell-type specific functions and ciliary heterogeneity may contribute to the range of phenotypes associated with ciliopathies.

Craniofacial developmental biology in the single-cell era.

The evolution of a unique craniofacial complex in vertebrates made possible new ways of breathing, eating, communicating and sensing the environment. The head and face develop through interactions of all three germ layers, the endoderm, ectoderm and mesoderm, as well as the so-called fourth germ layer, the cranial neural crest. Over a century of experimental embryology and genetics have revealed an incredible diversity of cell types derived from each germ layer, signaling pathways and genes that coordinate craniofacial development, and how changes to these underlie human disease and vertebrate evolution. Yet for many diseases and congenital anomalies, we have an incomplete picture of the causative genomic changes, in particular how alterations to the non-coding genome might affect craniofacial gene expression. Emerging genomics and single-cell technologies provide an opportunity to obtain a more holistic view of the genes and gene regulatory elements orchestrating craniofacial development across vertebrates. These single-cell studies generate novel hypotheses that can be experimentally validated in vivo. In this Review, we highlight recent advances in single-cell studies of diverse craniofacial structures, as well as potential pitfalls and the need for extensive in vivo validation. We discuss how these studies inform the developmental sources and regulation of head structures, bringing new insights into the etiology of structural birth anomalies that affect the vertebrate head.

Endothelin signaling in development.

Since the discovery of endothelin 1 (EDN1) in 1988, the role of endothelin ligands and their receptors in the regulation of blood pressure in normal and disease states has been extensively studied. However, endothelin signaling also plays crucial roles in the development of neural crest cell-derived tissues. Mechanisms of endothelin action during neural crest cell maturation have been deciphered using a variety of in vivo and in vitro approaches, with these studies elucidating the basis of human syndromes involving developmental differences resulting from altered endothelin signaling. In this Review, we describe the endothelin pathway and its functions during the development of neural crest-derived tissues. We also summarize how dysregulated endothelin signaling causes developmental differences and how this knowledge may lead to potential treatments for individuals with gene variants in the endothelin pathway.

Craniofacial dysmorphology in Down syndrome is caused by increased dosage of Dyrk1a and at least three other genes.

Down syndrome (DS), trisomy of human chromosome 21 (Hsa21), occurs in 1 in 800 live births and is the most common human aneuploidy. DS results in multiple phenotypes, including craniofacial dysmorphology, which is characterised by midfacial hypoplasia, brachycephaly and micrognathia. The genetic and developmental causes of this are poorly understood. Using morphometric analysis of the Dp1Tyb mouse model of DS and an associated mouse genetic mapping panel, we demonstrate that four Hsa21-orthologous regions of mouse chromosome 16 contain dosage-sensitive genes that cause the DS craniofacial phenotype, and identify one of these causative genes as Dyrk1a. We show that the earliest and most severe defects in Dp1Tyb skulls are in bones of neural crest (NC) origin, and that mineralisation of the Dp1Tyb skull base synchondroses is aberrant. Furthermore, we show that increased dosage of Dyrk1a results in decreased NC cell proliferation and a decrease in size and cellularity of the NC-derived frontal bone primordia. Thus, DS craniofacial dysmorphology is caused by an increased dosage of Dyrk1a and at least three other genes.

Reassessing the embryonic origin and potential of craniofacial ectomesenchyme.

Of all the cell types arising from the neural crest, ectomesenchyme is likely the most unusual. In contrast to the neuroglial cells generated by neural crest throughout the embryo, consistent with its ectodermal origin, cranial neural crest-derived cells (CNCCs) generate many connective tissue and skeletal cell types in common with mesoderm. Whether this ectoderm-derived mesenchyme (ectomesenchyme) potential reflects a distinct developmental origin from other CNCC lineages, and/or epigenetic reprogramming of the ectoderm, remains debated. Whereas decades of lineage tracing studies have defined the potential of CNCC ectomesenchyme, these are being revisited by modern genetic techniques. Recent work is also shedding light on the extent to which intrinsic and extrinsic cues determine ectomesenchyme potential, and whether maintenance or reacquisition of CNCC multipotency influences craniofacial repair.

Dyrk1a is required for craniofacial development in Xenopus laevis.

Loss of function variations in the dual specificity tyrosine-phosphorylation-regulated kinase 1 A (DYRK1A) gene are associated with craniofacial malformations in humans. Here we characterized the effects of deficient DYRK1A in craniofacial development using a developmental model, Xenopus laevis. Dyrk1a mRNA and protein were expressed throughout the developing head and both were enriched in the branchial arches which contribute to the face and jaw. Consistently, reduced Dyrk1a function, using dyrk1a morpholinos and pharmacological inhibitors, resulted in orofacial malformations including hypotelorism, altered mouth shape, slanted eyes, and narrower face accompanied by smaller jaw cartilage and muscle. Inhibition of Dyrk1a function resulted in misexpression of key craniofacial regulators including transcription factors and members of the retinoic acid signaling pathway. Two such regulators, sox9 and pax3 are required for neural crest development and their decreased expression corresponds with smaller neural crest domains within the branchial arches. Finally, we determined that the smaller size of the faces, jaw elements and neural crest domains in embryos deficient in Dyrk1a could be explained by increased cell death and decreased proliferation. This study is the first to provide insight into why craniofacial birth defects might arise in humans with variants of DYRK1A.

Exploring the roles of noncoding RNAs in craniofacial abnormalities: A systematic review.

Congenital craniofacial abnormalities are congenital anomalies of variable expressivity and severity with a recognizable set of abnormalities, which are derived from five identifiable primordial structures. They can occur unilaterally or bilaterally and include various malformations such as cleft lip with/without palate, craniosynostosis, and craniofacial microsomia. To date, the molecular etiology of craniofacial abnormalities is largely unknown. Noncoding RNAs (ncRNAs), including microRNAs, long ncRNAs, circular RNAs and PIWI-interacting RNAs, function as major regulators of cellular epigenetic hallmarks via regulation of various molecular and cellular processes. Recently, aberrant expression of ncRNAs has been implicated in many diseases, including craniofacial abnormalities. Consequently, this review focuses on the role and mechanism of ncRNAs in regulating craniofacial development in the hope of providing clues to identify potential therapeutic targets.

Shared genetic risk between major orofacial cleft phenotypes in an African population.

Nonsyndromic orofacial clefts (NSOFCs) represent a large proportion (70%-80%) of all OFCs. They can be broadly categorized into nonsyndromic cleft lip with or without cleft palate (NSCL/P) and nonsyndromic cleft palate only (NSCPO). Although NSCL/P and NSCPO are considered etiologically distinct, recent evidence suggests the presence of shared genetic risks. Thus, we investigated the genetic overlap between NSCL/P and NSCPO using African genome-wide association study (GWAS) data on NSOFCs. These data consist of 814 NSCL/P, 205 NSCPO cases, and 2159 unrelated controls. We generated common single-nucleotide variants (SNVs) association summary statistics separately for each phenotype (NSCL/P and NSCPO) under an additive genetic model. Subsequently, we employed the pleiotropic analysis under the composite null (PLACO) method to test for genetic overlap. Our analysis identified two loci with genome-wide significance (rs181737795 [p = 2.58E-08] and rs2221169 [p = 4.5E-08]) and one locus with marginal significance (rs187523265 [p = 5.22E-08]). Using mouse transcriptomics data and information from genetic phenotype databases, we identified MDN1, MAP3k7, KMT2A, ARCN1, and VADC2 as top candidate genes for the associated SNVs. These findings enhance our understanding of genetic variants associated with NSOFCs and identify potential candidate genes for further exploration.