Severe nasal clefting and abnormal embryonic apoptosis in Alx3/Alx4 double mutant mice

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.

Clinical and Genetic Correlation in Neurocristopathies: Bridging a Precision Medicine Gap

Neurocristopathies (NCPs) encompass a spectrum of disorders arising from issues during the formation and migration of neural crest cells (NCCs). NCCs undergo epithelial-mesenchymal transition (EMT) and upon key developmental gene deregulation, fetuses and neonates are prone to exhibit diverse manifestations depending on the affected area. These conditions are generally rare and often have a genetic basis, with many following Mendelian inheritance patterns, thus making them perfect candidates for precision medicine. Examples include cranial NCPs, like Goldenhar syndrome and Axenfeld-Rieger syndrome; cardiac-vagal NCPs, such as DiGeorge syndrome; truncal NCPs, like congenital central hypoventilation syndrome and Waardenburg syndrome; and enteric NCPs, such as Hirschsprung disease. Additionally, NCCs' migratory and differentiating nature makes their derivatives prone to tumors, with various cancer types categorized based on their NCC origin. Representative examples include schwannomas and pheochromocytomas. This review summarizes current knowledge of diseases arising from defects in NCCs' specification and highlights the potential of precision medicine to remedy a clinical phenotype by targeting the genotype, particularly important given that those affected are primarily infants and young children.

Unraveling the transcriptomic landscape of eye migration and visual adaptations during flatfish metamorphosis

Flatfish undergo a remarkable metamorphosis from symmetrical pelagic larvae to fully asymmetrical benthic juveniles. The most distinctive features of this transformation is the migration of one eye. The molecular role of thyroid hormone in the metamorphosis process in flatfishes is well established. However, the regulatory network that facilitates eye movement remains enigmatic. This paper presents a morphological investigation of the metamorphic process in turbot eyes, using advanced imaging techniques and a global view of gene expression. The study covers migrant and non-migrant eyes and aims to identify the genes that are active during ocular migration. Our transcriptomic analysis shows a significant up-regulation of immune-related genes. The analysis of eye-specific genes reveals distinct patterns during the metamorphic process. Myosin is highlighted in the non-migrant eye, while ependymin is highlighted in the migrant eye, possibly involved in optic nerve regeneration. Furthermore, a potential association between the alx3 gene and cranial restructuring has been identified. Additionally, it confirmed simultaneous adaptation to low light in both eyes, as described by changes in opsins expression during the metamorphic process. The study also revealed that ocular migration activates systems asynchronously in both eyes, providing insight into multifaceted reorganization processes during metamorphosis of flatfish.

Loss of Neogenin alters branchial arch development and leads to craniofacial skeletal defects

The formation of complex structures, such as the craniofacial skeleton, requires precise and intricate two-way signalling between populations of cells of different embryonic origins. For example, the lower jaw, or mandible, arises from cranial neural crest cells (CNCCs) in the mandibular portion of the first branchial arch (mdBA1) of the embryo, and its development is regulated by signals from the ectoderm and cranial mesoderm (CM) within this structure. The molecular mechanisms underlying CM cell influence on CNCC development in the mdBA1 remain poorly defined. Herein we identified the receptor Neogenin as a key regulator of craniofacial development. We found that ablation of Neogenin expression via gene-targeting resulted in several craniofacial skeletal defects, including reduced size of the CNCC-derived mandible. Loss of Neogenin did not affect the formation of the mdBA1 CM core but resulted in altered Bmp4 and Fgf8 expression, increased apoptosis, and reduced osteoblast differentiation in the mdBA1 mesenchyme. Reduced BMP signalling in the mdBA1 of Neogenin mutant embryos was associated with alterations in the gene regulatory network, including decreased expression of transcription factors of the Hand, Msx, and Alx families, which play key roles in the patterning and outgrowth of the mdBA1. Tissue-specific Neogenin loss-of-function studies revealed that Neogenin expression in mesodermal cells contributes to mandible formation. Thus, our results identify Neogenin as a novel regulator of craniofacial skeletal formation and demonstrates it impinges on CNCC development via a non-cell autonomous mechanism.

DNA-guided transcription factor cooperativity shapes face and limb mesenchyme

Transcription factors (TFs) can define distinct cellular identities despite nearly identical DNA-binding specificities. One mechanism for achieving regulatory specificity is DNA-guided TF cooperativity. Although in vitro studies suggest that it may be common, examples of such cooperativity remain scarce in cellular contexts. Here, we demonstrate how "Coordinator," a long DNA motif composed of common motifs bound by many basic helix-loop-helix (bHLH) and homeodomain (HD) TFs, uniquely defines the regulatory regions of embryonic face and limb mesenchyme. Coordinator guides cooperative and selective binding between the bHLH family mesenchymal regulator TWIST1 and a collective of HD factors associated with regional identities in the face and limb. TWIST1 is required for HD binding and open chromatin at Coordinator sites, whereas HD factors stabilize TWIST1 occupancy at Coordinator and titrate it away from HD-independent sites. This cooperativity results in the shared regulation of genes involved in cell-type and positional identities and ultimately shapes facial morphology and evolution.

TFAP2 paralogs regulate midfacial development in part through a conserved ALX genetic pathway

Cranial neural crest development is governed by positional gene regulatory networks (GRNs). Fine-tuning of the GRN components underlies facial shape variation, yet how those networks in the midface are connected and activated remain poorly understood. Here, we show that concerted inactivation of Tfap2a and Tfap2b in the murine neural crest, even during the late migratory phase, results in a midfacial cleft and skeletal abnormalities. Bulk and single-cell RNA-seq profiling reveal that loss of both TFAP2 family members dysregulates numerous midface GRN components involved in midface morphogenesis, patterning and differentiation. Notably, Alx1, Alx3 and Alx4 (ALX) transcript levels are reduced, whereas ChIP-seq analyses suggest TFAP2 family members directly and positively regulate ALX gene expression. Tfap2a, Tfap2b and ALX co-expression in midfacial neural crest cells of both mouse and zebrafish implies conservation of this regulatory axis across vertebrates. Consistent with this notion, tfap2a zebrafish mutants present with abnormal alx3 expression patterns, Tfap2a binds ALX loci and tfap2a-alx3 genetic interactions are observed. Together, these data demonstrate TFAP2 paralogs regulate vertebrate midfacial development in part by activating expression of ALX transcription factor genes.

Developmental origin of the mammalian premaxilla

The evolution of jaws has played a major role in the success of vertebrate expansion into a wide variety of ecological niches. A fundamental, yet unresolved, question in craniofacial biology is about the origin of the premaxilla, the most distal bone present in the upper jaw of all amniotes. Recent reports have suggested that the mammalian premaxilla is derived from embryonic maxillary prominences rather than the frontonasal ectomesenchyme as previously shown in studies of chicken embryos. However, whether mammalian embryonic frontonasal ectomesenchyme contributes to the premaxillary bone has not been investigated and a tool to trace the contributions of the frontonasal ectomesenchyme to facial structures in mammals is lacking. The expression of the Alx3 gene is activated highly specifically in the frontonasal ectomesenchyme, but not in the maxillary mesenchyme, from the beginning of facial morphogenesis in mice. Here, we report the generation and characterization of a novel Alx3CreERT2 knock-in mouse line that express tamoxifen-inducible Cre DNA recombinase from the Alx3 locus. Tamoxifen treatment of Alx3CreERT2/+;Rosa26mTmG/+ embryos at E7.5, E8.5, E9.5, and E10.5, each induced specific labeling of the embryonic medial nasal and lateral nasal mesenchyme but not the maxillary mesenchyme. Lineage tracing of Alx3CreERT2-labeled frontonasal mesenchyme from E9.5 to E16.5 clearly showed that the frontonasal mesenchyme cells give rise to the osteoblasts generating the premaxillary bone. Furthermore, we characterize a Dlx1-Cre BAC transgenic mouse line that expresses Cre activity in the embryonic maxillary but not the frontonasal mesenchyme and show that the Dlx1-Cre labeled embryonic maxillary mesenchyme cells contribute to the maxillary bone as well as the soft tissues lateral to both the premaxillary and maxillary bones but not to the premaxillary bone. These results clearly demonstrate the developmental origin of the premaxillary bone from embryonic frontonasal ectomesenchyme cells in mice and confirm the evolutionary homology of the premaxilla across amniotes.

Systematic identification of conditionally folded intrinsically disordered regions by AlphaFold2

The AlphaFold Protein Structure Database contains predicted structures for millions of proteins. For the majority of human proteins that contain intrinsically disordered regions (IDRs), which do not adopt a stable structure, it is generally assumed that these regions have low AlphaFold2 confidence scores that reflect low-confidence structural predictions. Here, we show that AlphaFold2 assigns confident structures to nearly 15% of human IDRs. By comparison to experimental NMR data for a subset of IDRs that are known to conditionally fold (i.e., upon binding or under other specific conditions), we find that AlphaFold2 often predicts the structure of the conditionally folded state. Based on databases of IDRs that are known to conditionally fold, we estimate that AlphaFold2 can identify conditionally folding IDRs at a precision as high as 88% at a 10% false positive rate, which is remarkable considering that conditionally folded IDR structures were minimally represented in its training data. We find that human disease mutations are nearly fivefold enriched in conditionally folded IDRs over IDRs in general and that up to 80% of IDRs in prokaryotes are predicted to conditionally fold, compared to less than 20% of eukaryotic IDRs. These results indicate that a large majority of IDRs in the proteomes of human and other eukaryotes function in the absence of conditional folding, but the regions that do acquire folds are more sensitive to mutations. We emphasize that the AlphaFold2 predictions do not reveal functionally relevant structural plasticity within IDRs and cannot offer realistic ensemble representations of conditionally folded IDRs.

aristaless-like homeobox-3 is wound induced and promotes a low-Wnt environment required for planarian head regeneration

The planarian Schmidtea mediterranea is a well-established model of adult regeneration, which is dependent on a large population of adult stem cells called neoblasts. Upon amputation, planarians undergo transcriptional wounding programs and coordinated stem cell proliferation to give rise to missing tissues. Interestingly, the Wnt signaling pathway is key to guiding what tissues are regenerated, yet less known are the transcriptional regulators that ensure proper activation and timing of signaling pathway components. Here, we have identified an aristaless-like homeobox transcription factor, alx-3, that is enriched in a population of putative neural-fated progenitor cells at homeostasis, and is also upregulated in stem cells and muscle cells at anterior-facing wounds upon amputation. Knockdown of alx-3 results in failure of head regeneration and patterning defects in amputated tail fragments. alx-3 is required for the expression of several early wound-induced genes, including the Wnt inhibitor notum, which is required to establish anterior polarity during regeneration. Together, these findings reveal a role for alx-3 as an early wound-response transcriptional regulator in both muscle cells and stem cells that is required for anterior regeneration by promoting a low-Wnt environment.

Linking Vertebrate Gene Duplications to the New Head Hypothesis

Vertebrates have diverse morphologies and various anatomical novelties that set them apart from their closest invertebrate relatives. A conspicuous head housing a large brain, paired sense organs, and protected by a skeleton of cartilage and bone is unique to vertebrates and is a defining feature of this taxon. Gans and Northcutt (1980s) proposed that the evolution of this "new head" was dependent on two key developmental innovations: neural crest cells (NCCs) and ectodermal placodes. NCCs are migratory embryonic cells that form bone, cartilage, and neurons in the new head. Based on genome size, Ohno (1970s) proposed a separate hypothesis, stating that vertebrate genome content was quadrupled via two rounds (2R) of whole genome duplications (WGDs), and the surplus of genetic material potentiated vertebrate morphological diversification. While both hypotheses offer explanations for vertebrate success, it is unclear if, and how, the "new head" and "2R" hypotheses are linked. Here, we consider both hypotheses and evaluate the experimental evidence connecting the two. Overall, evidence suggests that while the origin of the NC GRN predates the vertebrate WGDs, these genomic events may have potentiated the evolution of distinct genetic subnetworks in different neural crest subpopulations. We describe the general composition of the NC GRN and posit that its increased developmental modularity facilitated the independent evolution of NC derivatives and the diversification of the vertebrate head skeleton. Lastly, we discuss experimental strategies needed to test whether gene duplications drove the diversification of neural crest derivatives and the "new head".

Creating morphological diversity in reptilian temporal skull region: A review of potential developmental mechanisms

Reptilian skull morphology is highly diverse and broadly categorized into three categories based on the number and position of the temporal fenestrations: anapsid, synapsid, and diapsid. According to recent phylogenetic analysis, temporal fenestrations evolved twice independently in amniotes, once in Synapsida and once in Diapsida. Although functional aspects underlying the evolution of tetrapod temporal fenestrations have been well investigated, few studies have investigated the developmental mechanisms responsible for differences in the pattern of temporal skull region. To determine what these mechanisms might be, we first examined how the five temporal bones develop by comparing embryonic cranial osteogenesis between representative extant reptilian species. The pattern of temporal skull region may depend on differences in temporal bone growth rate and growth direction during ontogeny. Next, we compared the histogenesis patterns and the expression of two key osteogenic genes, Runx2 and Msx2, in the temporal region of the representative reptilian embryos. Our comparative analyses suggest that the embryonic histological condition of the domain where temporal fenestrations would form predicts temporal skull morphology in adults and regulatory modifications of Runx2 and Msx2 expression in osteogenic mesenchymal precursor cells are likely involved in generating morphological diversity in the temporal skull region of reptiles.

Genome Wide Association Study with Imputed Whole Genome Sequence Data Identifies a 431 kb Risk Haplotype on CFA18 for Congenital Laryngeal Paralysis in Alaskan Sled Dogs

Congenital laryngeal paralysis (CLP) is an inherited disorder that affects the ability of the dog to exercise and precludes it from functioning as a working sled dog. Though CLP is known to occur in Alaskan sled dogs (ASDs) since 1986, the genetic mutation underlying the disease has not been reported. Using a genome-wide association study (GWAS), we identified a 708 kb region on CFA 18 harboring 226 SNPs to be significantly associated with CLP. The significant SNPs explained 47.06% of the heritability of CLP. We narrowed the region to 431 kb through autozygosity mapping and found 18 of the 20 cases to be homozygous for the risk haplotype. Whole genome sequencing of two cases and a control ASD, and comparison with the genome of 657 dogs from various breeds, confirmed the homozygous status of the risk haplotype to be unique to the CLP cases. Most of the dogs that were homozygous for the risk allele had blue eyes. Gene annotation and a gene-based association study showed that the risk haplotype encompasses genes implicated in developmental and neurodegenerative disorders. Pathway analysis showed enrichment of glycoproteins and glycosaminoglycans biosynthesis, which play a key role in repairing damaged nerves. In conclusion, our results suggest an important role for the identified candidate region in CLP.

Single-cell transcriptomic signatures and gene regulatory networks modulated by Wls in mammalian midline facial formation and clefts

Formation of highly unique and complex facial structures is controlled by genetic programs that are responsible for the precise coordination of three-dimensional tissue morphogenesis. However, the underlying mechanisms governing these processes remain poorly understood. We combined mouse genetic and genomic approaches to define the mechanisms underlying normal and defective midfacial morphogenesis. Conditional inactivation of the Wnt secretion protein Wls in Pax3-expressing lineage cells disrupted frontonasal primordial patterning, cell survival and directional outgrowth, resulting in altered facial structures, including midfacial hypoplasia and midline facial clefts. Single-cell RNA sequencing revealed unique transcriptomic atlases of mesenchymal subpopulations in the midfacial primordia, which are disrupted in the conditional Wls mutants. Differentially expressed genes and cis-regulatory sequence analyses uncovered that Wls modulates and integrates a core gene regulatory network, consisting of key midfacial regulatory transcription factors (including Msx1, Pax3 and Pax7) and their downstream targets (including Wnt, Shh, Tgfβ and retinoic acid signaling components), in a mesenchymal subpopulation of the medial nasal prominences that is responsible for midline facial formation and fusion. These results reveal fundamental mechanisms underlying mammalian midfacial morphogenesis and related defects at single-cell resolution.

ALX1-transcribed LncRNA AC132217.4 promotes osteogenesis and bone healing via IGF-AKT signaling in mesenchymal stem cells

The osteogenic potential of bone marrow mesenchymal stem cells (BMSCs) is critical for bone formation and regeneration. A high non-/delayed-union rate of fracture healing still occurs in specific populations, implying an urgent need to discover novel targets for promoting osteogenesis and bone regeneration. Long non-coding (lnc)RNAs are emerging regulators of multiple physiological processes, including osteogenesis. Based on differential expression analysis of RNA sequencing data, we found that lncRNA AC132217.4, a 3'UTR-overlapping lncRNA of insulin growth factor 2 (IGF2), was highly induced during osteogenic differentiation of BMSCs. Afterward, both gain-of-function and loss-of-function experiments proved that AC132217.4 promotes osteoblast development from BMSCs. As for its molecular mechanism, we found that AC132217.4 binds with IGF2 mRNA to regulate its expression and downstream AKT activation to control osteoblast maturation and function. Furthermore, we identified two splicing factors, splicing component 35 KDa (SC35) and heterogeneous nuclear ribonucleoprotein A1 (HNRNPA1), which regulate the biogenesis of AC132217.4 at the post-transcriptional level. We also identified a transcription factor, ALX1, which regulates AC132217.7 expression at the transcriptional level to promote osteogenesis. Importantly, in-vivo over-expression of AC132217.4 essentially promotes the bone healing process in a murine tibial drill-hole model. Our study demonstrates that lncRNA AC132217.4 is a novel anabolic regulator of BMSC osteogenesis and could be a plausible therapeutic target for improving bone regeneration.

Zebrafish models of alx-linked frontonasal dysplasia reveal a role for Alx1 and Alx3 in the anterior segment and vasculature of the developing eye

The cellular and genetic mechanisms that coordinate formation of facial sensory structures with surrounding skeletal and soft tissue elements remain poorly understood. Alx1, a homeobox transcription factor, is a key regulator of midfacial morphogenesis. ALX1 mutations in humans are linked to severe congenital anomalies of the facial skeleton (frontonasal dysplasia, FND) with malformation or absence of eyes and orbital contents (micro- and anophthalmia). Zebrafish with loss-of-function alx1 mutations develop with craniofacial and ocular defects of variable penetrance, likely due to compensatory upregulation in expression of a paralogous gene, alx3. Here we show that zebrafish alx1;alx3 mutants develop with highly penetrant cranial and ocular defects that resemble human ALX1-linked FND. alx1 and alx3 are expressed in anterior cranial neural crest (aCNC), which gives rise to the anterior neurocranium (ANC), anterior segment structures of the eye and vascular pericytes. Consistent with a functional requirement for alx genes in aCNC, alx1; alx3 mutants develop with nearly absent ANC and grossly aberrant hyaloid vasculature and ocular anterior segment, but normal retina. In vivo lineage labeling identified a requirement for alx1 and alx3 during aCNC migration, and transcriptomic analysis suggested oxidative stress response as a key target mechanism of this function. Oxidative stress is a hallmark of fetal alcohol toxicity, and we found increased penetrance of facial and ocular malformations in alx1 mutants exposed to ethanol, consistent with a protective role for alx1 against ethanol toxicity. Collectively, these data demonstrate a conserved role for zebrafish alx genes in controlling ocular and facial development, and a novel role in protecting these key midfacial structures from ethanol toxicity during embryogenesis. These data also reveal novel roles for alx genes in ocular anterior segment formation and vascular development and suggest that retinal deficits in alx mutants may be secondary to aberrant ocular vascularization and anterior segment defects. This study establishes robust zebrafish models for interrogating conserved genetic mechanisms that coordinate facial and ocular development, and for exploring gene--environment interactions relevant to fetal alcohol syndrome.

Transcriptional Regulation of Jaw Osteoblasts: Development to Pathology

Craniofacial and jaw bones have unique physiological specificities when compared to axial and appendicular bones. However, the molecular profile of the jaw osteoblast (OB) remains incomplete. The present study aimed to decipher the bone site-specific profiles of transcription factors (TFs) expressed in OBs in vivo. Using RNA sequencing analysis, we mapped the transcriptome of confirmed OBs from 2 different skeletal sites: mandible (Md) and tibia (Tb). The OB transcriptome contains 709 TF genes: 608 are similarly expressed in Md-OB and Tb-OB, referred to as "OB-core"; 54 TF genes are upregulated in Md-OB, referred to as "Md-set"; and 18 TF genes are upregulated in Tb-OB, referred to as "Tb-set." Notably, the expression of 29 additional TF genes depends on their RNA transcript variants. TF genes with no previously known role in OBs and bone were identified. Bioinformatics analysis combined with review of genetic disease databases and a comprehensive literature search showed a significant contribution of anatomical origin to the OB signatures. Md-set and Tb-set are enriched with site-specific TF genes associated with development and morphogenesis (neural crest vs. mesoderm), and this developmental imprint persists during growth and homeostasis. Jaw and tibia site-specific OB signatures are associated with craniofacial and appendicular skeletal disorders as well as neurocristopathies, dental disorders, and digit malformations. The present study demonstrates the feasibility of a new method to isolate pure OB populations and map their gene expression signature in the context of OB physiological environment, avoiding in vitro culture and its associated biases. Our results provide insights into the site-specific developmental pathways governing OBs and identify new major OB regulators of bone physiology. We also established the importance of the OB transcriptome as a prognostic tool for human rare bone diseases to explore the hidden pathophysiology of craniofacial malformations, among the most prevalent congenital defects in humans.

Alx1 Deficient Mice Recapitulate Craniofacial Phenotype and Reveal Developmental Basis of ALX1-Related Frontonasal Dysplasia

Loss of ALX1 function causes the frontonasal dysplasia syndrome FND3, characterized by severe facial clefting and microphthalmia. Whereas the laboratory mouse has been the preeminent animal model for studying developmental mechanisms of human craniofacial birth defects, the roles of ALX1 in mouse frontonasal development have not been well characterized because the only previously reported Alx1 mutant mouse line exhibited acrania due to a genetic background-dependent failure of cranial neural tube closure. Using CRISPR/Cas9-mediated genome editing, we have generated an Alx1-deletion mouse model that recapitulates the FND craniofacial malformations, including median orofacial clefting and disruption of development of the eyes and alae nasi. In situ hybridization analysis showed that Alx1 is strongly expressed in frontonasal neural crest cells that give rise to periocular and frontonasal mesenchyme. Alx1del/del embryos exhibited increased apoptosis of periocular mesenchyme and decreased expression of ocular developmental regulators Pitx2 and Lmxb1 in the periocular mesenchyme, followed by defective optic stalk morphogenesis. Moreover, Alx1del/del embryos exhibited disruption of frontonasal mesenchyme identity, with loss of expression of Pax7 and concomitant ectopic expression of the jaw mesenchyme regulators Lhx6 and Lhx8 in the developing lateral nasal processes. The function of ALX1 in patterning the frontonasal mesenchyme is partly complemented by ALX4, a paralogous ALX family transcription factor whose loss-of-function causes a milder and distinctive FND. Together, these data uncover previously unknown roles of ALX1 in periocular mesenchyme development and frontonasal mesenchyme patterning, providing novel insights into the pathogenic mechanisms of ALX1-related FND.

Retinoic Acid Deficiency Underlies the Etiology of Midfacial Defects

Embryonic craniofacial development depends on the coordinated outgrowth and fusion of multiple facial primordia, which are populated with cranial neural crest cells and covered by the facial ectoderm. Any disturbance in these developmental events, their progenitor tissues, or signaling pathways can result in craniofacial deformities such as orofacial clefts, which are among the most common birth defects in humans. In the present study, we show that Rdh10 loss of function leads to a substantial reduction in retinoic acid (RA) signaling in the developing frontonasal process during early embryogenesis, which results in a variety of craniofacial anomalies, including midfacial cleft and ectopic chondrogenic nodules. Elevated apoptosis and perturbed cell proliferation in postmigratory cranial neural crest cells and a substantial reduction in Alx1 and Alx3 transcription in the developing frontonasal process were associated with midfacial cleft in Rdh10-deficient mice. More important, expanded Shh signaling in the ventral forebrain, as well as partial abrogation of midfacial defects in Rdh10 mutants via inhibition of Hh signaling, indicates that misregulation of Shh signaling underlies the pathogenesis of reduced RA signaling-associated midfacial defects. Taken together, these data illustrate the precise spatiotemporal function of Rdh10 and RA signaling during early embryogenesis and their importance in orchestrating molecular and cellular events essential for normal midfacial development.

Analysis of the DNA-binding properties of Alx1, an evolutionarily conserved regulator of skeletogenesis in echinoderms

Alx1, a homeodomain-containing transcription factor, is a highly conserved regulator of skeletogenesis in echinoderms. In sea urchins, Alx1 plays a central role in the differentiation of embryonic primary mesenchyme cells (PMCs) and positively regulates the transcription of most biomineralization genes expressed by these cells. The alx1 gene arose via duplication and acquired a skeletogenic function distinct from its paralog (alx4) through the exonization of a 41–amino acid motif (the D2 domain). Alx1 and Alx4 contain glutamine-50 paired-type homeodomains, which interact preferentially with palindromic binding sites in vitro. Chromatin immunoprecipitation sequencing (ChIP-seq) studies have shown, however, that Alx1 binds both to palindromic and half sites in vivo. To address this apparent discrepancy and explore the function of the D2 domain, we used an endogenous cis-regulatory module associated with Sp-mtmmpb, a gene that encodes a PMC-specific metalloprotease, to analyze the DNA-binding properties of Alx1. We find that Alx1 forms dimeric complexes on TAAT-containing half sites by a mechanism distinct from the well-known mechanism of dimerization on palindromic sites. We used transgenic reporter assays to analyze the functional roles of half sites in vivo and demonstrate that two sites with partially redundant functions are essential for the PMC-specific activity of the Sp-mtmmpb cis-regulatory module. Finally, we show that the D2 domain influences the DNA-binding properties of Alx1 in vitro, suggesting that the exonization of this motif may have facilitated the acquisition of new transcriptional targets and consequently a novel developmental function.

The alx3 gene shapes the zebrafish neurocranium by regulating frontonasal neural crest cell differentiation timing

During craniofacial development, different populations of cartilage- and bone-forming cells develop in precise locations in the head. Most of these cells are derived from pluripotent cranial neural crest cells and differentiate with distinct developmental timing and cellular morphologies. The mechanisms that divide neural crest cells into discrete populations are not fully understood. Here, we use single-cell RNA sequencing to transcriptomically define different populations of cranial neural crest cells. We discovered that the gene family encoding the Alx transcription factors is enriched in the frontonasal population of neural crest cells. Genetic mutant analyses indicate that alx3 functions to regulate the distinct differentiation timing and cellular morphologies among frontonasal neural crest cell subpopulations. This study furthers our understanding of how genes controlling developmental timing shape craniofacial skeletal elements.

SATB2: A versatile transcriptional regulator of craniofacial and skeleton development, neurogenesis and tumorigenesis, and its applications in regenerative medicine

SATB2 (special AT-rich sequence-binding protein 2) is a member of the special AT-rich binding protein family. As a transcription regulator, SATB2 mainly integrates higher-order chromatin organization. SATB2 expression appears to be tissue- and stage-specific, and is governed by several cellular signaling molecules and mediators. Expressed in branchial arches and osteoblast-lineage cells, SATB2 plays a significant role in craniofacial pattern and skeleton development. In addition to regulating osteogenic differentiation, SATB2 also displays versatile functions in neural development and cancer progression. As an osteoinductive factor, SATB2 holds great promise in improving bone regeneration toward bone defect repair. In this review, we have summarized our current understanding of the physiological and pathological functions of SATB2 in craniofacial and skeleton development, neurogenesis, tumorigenesis and regenerative medicine.

ALX1-related frontonasal dysplasia results from defective neural crest cell development and migration

A pedigree of subjects presented with frontonasal dysplasia (FND). Genome sequencing and analysis identified a p.L165F missense variant in the homeodomain of the transcription factor ALX1 which was imputed to be pathogenic. Induced pluripotent stem cells (iPSC) were derived from the subjects and differentiated to neural crest cells (NCC). NCC derived from ALX1L165F/L165F iPSC were more sensitive to apoptosis, showed an elevated expression of several neural crest progenitor state markers, and exhibited impaired migration compared to wild-type controls. NCC migration was evaluated in vivo using lineage tracing in a zebrafish model, which revealed defective migration of the anterior NCC stream that contributes to the median portion of the anterior neurocranium, phenocopying the clinical presentation. Analysis of human NCC culture media revealed a change in the level of bone morphogenic proteins (BMP), with a low level of BMP2 and a high level of BMP9. Soluble BMP2 and BMP9 antagonist treatments were able to rescue the defective migration phenotype. Taken together, these results demonstrate a mechanistic requirement of ALX1 in NCC development and migration.

Transcriptional analysis of cleft palate in TGFβ3 mutant mice

Cleft palate (CP) is one of the most common craniofacial birth defects, impacting about 1 in 800 births in the USA. Tgf-β3 plays a critical role in regulating murine palate development, and Tgf-β3 null mutants develop cleft palate with 100% penetrance. In this study, we compared global palatal transcriptomes of wild type (WT) and Tgf-β3 −/− homozygous (HM) mouse embryos at the crucial palatogenesis stages of E14.5, and E16.5, using RNA-seq data. We found 1,809 and 2,127 differentially expressed genes at E16.5 vs. E14.5 in the WT and HM groups, respectively (adjusted p < 0.05; |fold change|> 2.0). We focused on the genes that were uniquely up/downregulated in WT or HM at E16.5 vs. E14.5 to identify genes associated with CP. Systems biology analysis relating to cell behaviors and function of WT and HM specific genes identified functional non-Smad pathways and preference of apoptosis to epithelial-mesenchymal transition. We identified 24 HM specific and 11 WT specific genes that are CP-related and/or involved in Tgf-β3 signaling. We validated the expression of 29 of the 35 genes using qRT-PCR and the trend of mRNA expression is similar to that of RNA-seq data . Our results enrich our understanding of genes associated with CP that are directly or indirectly regulated via TGF-β.

Crucial and Overlapping Roles of Six1 and Six2 in Craniofacial Development

SIX1 and SIX2 encode closely related transcription factors of which disruptions have been associated with distinct craniofacial syndromes, with mutations in SIX1 associated with branchiootic syndrome 3 (BOS3) and heterozygous deletions of SIX2 associated with frontonasal dysplasia defects. Whereas mice deficient in Six1 recapitulated most of the developmental defects associated with BOS3, mice lacking Six2 function had no obvious frontonasal defects. We show that Six1 and Six2 exhibit partly overlapping patterns of expression in the developing mouse embryonic frontonasal, maxillary, and mandibular processes. We found that Six1 ^-/- Six2 ^-/- double-mutant mice were born with severe craniofacial deformity not seen in the Six1 ^-/- or Six2 ^-/- single mutants, including skull bone agenesis, midline facial cleft, and syngnathia. Moreover, whereas Six1 ^-/- mice exhibited partial transformation of maxillary zygomatic bone into a mandibular condyle-like structure, Six1 ^-/-Six2 ^+/- mice exhibit significantly increased penetrance of the maxillary malformation. In addition to ectopic Dlx5 expression at the maxillary-mandibular junction as recently reported in E10.5 Six1 ^-/- embryos, the E10.5 Six1 ^-/- Six2 ^+/- embryos showed ectopic expression of Bmp4, Msx1, and Msx2 messenger RNAs in the maxillary-mandibular junction. Genetically inactivating 1 allele of either Ednra or Bmp4 significantly reduced the penetrance of maxillary malformation in both Six1 ^-/- and Six1 ^-/- Six2 ^+/- embryos, indicating that Six1 and Six2 regulate both endothelin and bone morphogenetic protein-4 signaling pathways to pattern the facial structures. Furthermore, we show that neural crest-specific inactivation of Six1 in Six2 ^-/- embryos resulted in midline facial cleft and frontal bone agenesis. We show that Six1 ^-/- Six2 ^-/- embryos exhibit significantly reduced expression of key frontonasal development genes Alx1 and Alx3 as well as increased apoptosis in the developing frontonasal mesenchyme. Together, these results indicate that Six1 and Six2 function partly redundantly to control multiple craniofacial developmental processes and play a crucial neural crest cell-autonomous role in frontonasal morphogenesis.

Transcriptome analysis of Xenopus orofacial tissues deficient in retinoic acid receptor function

BACKGROUND

Development of the face and mouth is orchestrated by a large number of transcription factors, signaling pathways and epigenetic regulators. While we know many of these regulators, our understanding of how they interact with each other and implement changes in gene expression during orofacial development is still in its infancy. Therefore, this study focuses on uncovering potential cooperation between transcriptional regulators and one important signaling pathway, retinoic acid, during development of the midface.

RESULTS

Transcriptome analyses was performed on facial tissues deficient for retinoic acid receptor function at two time points in development; early (35 hpf) just after the neural crest migrates and facial tissues are specified and later (60 hpf) when the mouth has formed and facial structures begin to differentiate. Functional and network analyses revealed that retinoic acid signaling could cooperate with novel epigenetic factors and calcium-NFAT signaling during early orofacial development. At the later stage, retinoic acid may work with WNT and BMP and regulate homeobox containing transcription factors. Finally, there is an overlap in genes dysregulated in Xenopus embryos with median clefts with human genes associated with similar orofacial defects.

CONCLUSIONS

This study uncovers novel signaling pathways required for orofacial development as well as pathways that could interact with retinoic acid signaling during the formation of the face. We show that frog faces are an important tool for studying orofacial development and birth defects.

Derivation and characterization of putative craniofacial mesenchymal progenitor cells from human induced pluripotent stem cells

The introduction and widespread adoption of induced pluripotent stem cell (iPSC) technology has opened new avenues for craniofacial regenerative medicine. Neural crest cells (NCCs) are the precursor population to many craniofacial structures, including dental and periodontal structures, and iPSC-derived NCCs may, in the near future, offer an unlimited supply of patient-specific cells for craniofacial repair interventions. Here, we used an established protocol involving simultaneous Wnt signaling activation and TGF-β signaling inhibition to differentiate three human iPSC lines to cranial NCCs. We then derived a mesenchymal progenitor cell (NCC-MPCs) population with chondrogenic and osteogenic potential from cranial NCCs and investigated their similarity to widely studied human postnatal dental or periodontal stem/progenitor cells. NCC-MPCs were quite distinct from both their precursor cells (NCCs) and bone-marrow mesenchymal stromal cells, a stromal population of mesodermal origin. Despite their similarity with dental stem/progenitor cells, NCC-MPCs were clearly differentiated by a core set of 43 genes, including ACKR3 (CXCR7), whose expression (both at transcript and protein level) appear to be specific to NCC-MPCs. Altogether, our data demonstrate the feasibility of craniofacial mesenchymal progenitor derivation from human iPSCs through a neural crest-intermediate and set the foundation for future studies regarding their full differentiation repertoire and their in vivo existence.

Ectomesoderm and epithelial-mesenchymal transition-related genes in spiralian development

BACKGROUND

Spiralians (e.g., annelids, molluscs, and flatworms) possess two sources of mesoderm. One is from endodermal precursors (endomesoderm), which is considered to be the ancestral source in metazoans. The second is from ectoderm (ectomesoderm) and may represent a novel cell type in the Spiralia. In the mollusc Crepidula fornicata, ectomesoderm is derived from micromere daughters within the A and B cell quadrants. Their progeny lie along the anterolateral edges of the blastopore. There they undergo epithelial-mesenchymal transition (EMT), become rounded and undergo delamination/ingression. Subsequently, they assume the mesenchymal phenotype, and migrate beneath the surface ectoderm to differentiate various cell types, including muscles and pigment cells.

RESULTS

We examined expression of several genes whose homologs are known to regulate Type 1 EMT in other metazoans. Most of these genes were expressed within spiralian ectomesoderm during EMT.

CONCLUSIONS

We propose that spiralian ectomesoderm, which exhibits analogous cellular behaviors to other populations of mesenchymal cells, may be controlled by the same genes that drive EMT in other metazoans. Perhaps these genes comprise a conserved metazoan EMT gene regulatory network (GRN). This study represents the first step in elucidating the GRN controlling the development of a novel spiralian cell type (ectomesoderm). Developmental Dynamics 247:1097-1120, 2018. © 2018 Wiley Periodicals, Inc.

Direct-to-consumer DNA testing of 6,000 dogs reveals 98.6-kb duplication associated with blue eyes and heterochromia in Siberian Huskies

Consumer genomics enables genetic discovery on an unprecedented scale by linking very large databases of personal genomic data with phenotype information voluntarily submitted via web-based surveys. These databases are having a transformative effect on human genomics research, yielding insights on increasingly complex traits, behaviors, and disease by including many thousands of individuals in genome-wide association studies (GWAS). The promise of consumer genomic data is not limited to human research, however. Genomic tools for dogs are readily available, with hundreds of causal Mendelian variants already characterized, because selection and breeding have led to dramatic phenotypic diversity underlain by a simple genetic structure. Here, we report the results of the first consumer genomics study ever conducted in a non-human model: a GWAS of blue eyes based on more than 3,000 customer dogs with validation panels including nearly 3,000 more, the largest canine GWAS to date. We discovered a novel association with blue eyes on chromosome 18 (P = 1.3x10-68) and used both sequence coverage and microarray probe intensity data to identify the putative causal variant: a 98.6-kb duplication directly upstream of the Homeobox gene ALX4, which plays an important role in mammalian eye development. This duplication is largely restricted to Siberian Huskies, is strongly associated with the blue-eyed phenotype (chi-square P = 5.2x10-290), and is highly, but not completely, penetrant. These results underscore the power of consumer-data-driven discovery in non-human species, especially dogs, where there is intense owner interest in the personal genomic information of their pets, a high level of engagement with web-based surveys, and an underlying genetic architecture ideal for mapping studies.

Pre-transplantational Control of the Post-transplantational Fate of Human Pluripotent Stem Cell-Derived Cartilage

Cartilage pellets generated from ectomesenchymal progeny of human pluripotent stem cells (hPSCs) in vitro eventually show signs of commitment of chondrocytes to hypertrophic differentiation. When transplanted subcutaneously, most of the surviving pellets were fully mineralized by 8 weeks. In contrast, treatment with the adenylyl cyclase activator, forskolin, in vitro resulted in slightly enlarged cartilage pellets containing an increased proportion of proliferating immature chondrocytes that expressed very low levels of hypertrophic/terminally matured chondrocyte-specific genes. Forskolin treatment also enhanced hyaline cartilage formation by reducing type I collagen gene expression and increasing sulfated glycosaminoglycan accumulation in the developed cartilage. Chondrogenic mesoderm from hPSCs and dedifferentiated nasal chondrocytes responded similarly to forskolin. Furthermore, forskolin treatment in vitro increased the frequency at which the cartilage pellets maintained unmineralized chondrocytes after subcutaneous transplantation. Thus, the post-transplantational fate of chondrocytes originating from hPSC-derived chondroprogenitors can be controlled during their genesis in vitro.

•

Forskolin/cAMP suppresses/delays BMP-induced chondrocyte maturation in vitro

•

Forskolin supports chondrocyte proliferation and hyaline chondrogenesis in vitro

•

Forskolin suppresses osteogenesis and BMP signaling gene expression in cartilage

•

In vitro forskolin treatment improves in vivo maintenance of uncalcified cartilage

FACS-Seq analysis of Pax3-derived cells identifies non-myogenic lineages in the embryonic forelimb

Skeletal muscle in the forelimb develops during embryonic and fetal development and perinatally. While much is known regarding the molecules involved in forelimb myogenesis, little is known about the specific mechanisms and interactions. Migrating skeletal muscle precursor cells express Pax3 as they migrate into the forelimb from the dermomyotome. To compare gene expression profiles of the same cell population over time, we isolated lineage-traced Pax3⁺ cells (Pax3^EGFP) from forelimbs at different embryonic days. We performed whole transcriptome profiling via RNA-Seq of Pax3⁺ cells to construct gene networks involved in different stages of embryonic and fetal development. With this, we identified genes involved in the skeletal, muscular, vascular, nervous and immune systems. Expression of genes related to the immune, skeletal and vascular systems showed prominent increases over time, suggesting a non-skeletal myogenic context of Pax3-derived cells. Using co-expression analysis, we observed an immune-related gene subnetwork active during fetal myogenesis, further implying that Pax3-derived cells are not a strictly myogenic lineage, and are involved in patterning and three-dimensional formation of the forelimb through multiple systems.

Gene datasets associated with mouse cleft palate

This article presents data on genes associated with cleft palate (CP), retrieved through both a full-text systematic review and a mouse genome informatics (MGI) database search. In order to group CP-associated genes according to function, pathway, biological process, and cellular component, the genes were analyzed using category enrichment bioinformatics tools, the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO). This approach provides invaluable opportunities for the identification of candidate pathways and genes in CP research.

Exome sequencing revealed a novel nonsense variant in ALX3 gene underlying frontorhiny

Frontorhiny is one of the two forms of mid-facial malformations characterized by ocular hypertelorism, wide and short nasal ridge, bifid nasal tip, broad columella, widely separated nares, long and wide philtrum and V-shaped hairline. Sometimes these phenotypes are associated with ptosis and midline dermoid cysts. Frontorhiny inherits in an autosomal recessive pattern. Sequence variants in the Aristaless-like homeobox 3 (ALX3) gene underlying frontorhiny have been reported previously. Here, in the present study, we have investigated four patients in a consanguineous family of Pakistani origin segregating frontorhiny in autosomal recessive manner. Genome scan using 250k Nsp1 array followed by exome and Sanger sequence analysis revealed a novel homozygous nonsense variant (c.604C>T, p.Gln202*) in the ALX3 gene resulting in frontorhiny in the family. This is the first mutation in the ALX3 gene, underlying frontorhiny, in Pakistani population.

Alx4 relays sequential FGF signaling to induce lacrimal gland morphogenesis

The sequential use of signaling pathways is essential for the guidance of pluripotent progenitors into diverse cell fates. Here, we show that Shp2 exclusively mediates FGF but not PDGF signaling in the neural crest to control lacrimal gland development. In addition to preventing p53-independent apoptosis and promoting the migration of Sox10-expressing neural crests, Shp2 is also required for expression of the homeodomain transcription factor Alx4, which directly controls Fgf10 expression in the periocular mesenchyme that is necessary for lacrimal gland induction. We show that Alx4 binds an Fgf10 intronic element conserved in terrestrial but not aquatic animals, underlying the evolutionary emergence of the lacrimal gland system in response to an airy environment. Inactivation of ALX4/Alx4 causes lacrimal gland aplasia in both human and mouse. These results reveal a key role of Alx4 in mediating FGF-Shp2-FGF signaling in the neural crest for lacrimal gland development.

The dry eye disease caused by lacrimal gland dysgenesis is one of the most common ocular ailments. In this study, we show that Shp2 mediates the sequential use of FGF signaling in lacrimal gland development. Our study identifies Alx4 as a novel target of Shp2 signaling and a causal gene for lacrimal gland aplasia in humans. Given this result, there may also be a potential role for Alx4 in guiding pluripotent stem cells to produce lacrimal gland tissue. Finally, our data reveals an Alx4-Fgf10 regulatory unit broadly conserved in the diverse array of terrestrial animals from humans to reptiles, but not in aquatic animals such as amphibians and fish, which sheds light on how the lacrimal gland arose as an evolutionary innovation of terrestrial animals to adapt to their newfound exposure to an airy environment.

Mllt10 knockout mouse model reveals critical role of Af10-dependent H3K79 methylation in midfacial development

Epigenetic regulation is required to ensure the precise spatial and temporal pattern of gene expression that is necessary for embryonic development. Although the roles of some epigenetic modifications in embryonic development have been investigated in depth, the role of methylation at lysine 79 (H3K79me) is poorly understood. Dot1L, a unique methyltransferase for H3K79, forms complexes with distinct sets of co-factors. To further understand the role of H3K79me in embryogenesis, we generated a mouse knockout of Mllt10, the gene encoding Af10, one Dot1L complex co-factor. We find homozygous Mllt10 knockout mutants (Mllt10-KO) exhibit midline facial cleft. The midfacial defects of Mllt10-KO embryos correspond to hyperterolism and are associated with reduced proliferation of mesenchyme in developing nasal processes and adjacent tissue. We demonstrate that H3K79me level is significantly decreased in nasal processes of Mllt10-KO embryos. Importantly, we find that expression of AP2α, a gene critical for midfacial development, is directly regulated by Af10-dependent H3K79me, and expression AP2α is reduced specifically in nasal processes of Mllt10-KO embryos. Suppression of H3K79me completely mimicked the Mllt10-KO phenotype. Together these data are the first to demonstrate that Af10-dependent H3K79me is essential for development of nasal processes and adjacent tissues, and consequent midfacial formation.

Generation of an immortalized mouse embryonic palatal mesenchyme cell line

Palatogenesis is a complex morphogenetic process, disruptions in which result in highly prevalent birth defects in humans. In recent decades, the use of model systems such as genetically-modified mice, mouse palatal organ cultures and primary mouse embryonic palatal mesenchyme (MEPM) cultures has provided significant insight into the molecular and cellular defects underlying cleft palate. However, drawbacks in each of these systems have prevented high-throughput, large-scale studies of palatogenesis in vitro. Here, we report the generation of an immortalized MEPM cell line that maintains the morphology, migration ability, transcript expression and responsiveness to exogenous growth factors of primary MEPM cells, with increased proliferative potential over primary cultures. The immortalization method described in this study will facilitate the generation of palatal mesenchyme cells with an unlimited capacity for expansion from a single genetically-modified mouse embryo and enable mechanistic studies of palatogenesis that have not been possible using primary culture.

Systems biology of facial development: contributions of ectoderm and mesenchyme

The rapid increase in gene-centric biological knowledge coupled with analytic approaches for genomewide data integration provides an opportunity to develop systems-level understanding of facial development. Experimental analyses have demonstrated the importance of signaling between the surface ectoderm and the underlying mesenchyme are coordinating facial patterning. However, current transcriptome data from the developing vertebrate face is dominated by the mesenchymal component, and the contributions of the ectoderm are not easily identified. We have generated transcriptome datasets from critical periods of mouse face formation that enable gene expression to be analyzed with respect to time, prominence, and tissue layer. Notably, by separating the ectoderm and mesenchyme we considerably improved the sensitivity compared to data obtained from whole prominences, with more genes detected over a wider dynamic range. From these data we generated a detailed description of ectoderm-specific developmental programs, including pan-ectodermal programs, prominence- specific programs and their temporal dynamics. The genes and pathways represented in these programs provide mechanistic insights into several aspects of ectodermal development. We also used these data to identify co-expression modules specific to facial development. We then used 14 co-expression modules enriched for genes involved in orofacial clefts to make specific mechanistic predictions about genes involved in tongue specification, in nasal process patterning and in jaw development. Our multidimensional gene expression dataset is a unique resource for systems analysis of the developing face; our co-expression modules are a resource for predicting functions of poorly annotated genes, or for predicting roles for genes that have yet to be studied in the context of facial development; and our analytic approaches provide a paradigm for analysis of other complex developmental programs.

Transcriptome profiling reveals expression signatures of cranial neural crest cells arising from different axial levels

Background

Cranial neural crest cells (NCCs) are a unique embryonic cell type which give rise to a diverse array of derivatives extending from neurons and glia through to bone and cartilage. Depending on their point of origin along the antero-posterior axis cranial NCCs are rapidly sorted into distinct migratory streams that give rise to axial specific structures. These migratory streams mirror the underlying segmentation of the brain with NCCs exiting the diencephalon and midbrain following distinct paths compared to those exiting the hindbrain rhombomeres (r). The genetic landscape of cranial NCCs arising at different axial levels remains unknown.

Results

Here we have used RNA sequencing to uncover the transcriptional profiles of mouse cranial NCCs arising at different axial levels. Whole transcriptome analysis identified over 120 transcripts differentially expressed between NCCs arising anterior to r3 (referred to as r1-r2 migratory stream for simplicity) and the r4 migratory stream. Eight of the genes differentially expressed between these populations were validated by RT-PCR with 2 being further validated by in situ hybridisation. We also explored the expression of the Neuropilins (Nrp1 and Nrp2) and their co-receptors and show that the A-type Plexins are differentially expressed in different cranial NCC streams.

Conclusions

Our analyses identify a large number of genes differentially regulated between cranial NCCs arising at different axial levels. This data provides a comprehensive description of the genetic landscape driving diversity of distinct cranial NCC streams and provides novel insight into the regulatory networks controlling the formation of specific skeletal elements and the mechanisms promoting migration along different paths.

Electronic supplementary material

The online version of this article (doi:10.1186/s12861-017-0147-z) contains supplementary material, which is available to authorized users.

Embryonic defence mechanisms against glucose-dependent oxidative stress require enhanced expression of Alx3 to prevent malformations during diabetic pregnancy

Oxidative stress constitutes a major cause for increased risk of congenital malformations associated to severe hyperglycaemia during pregnancy. Mutations in the gene encoding the transcription factor ALX3 cause congenital craniofacial and neural tube defects. Since oxidative stress and lack of ALX3 favour excessive embryonic apoptosis, we investigated whether ALX3-deficiency further increases the risk of embryonic damage during gestational hyperglycaemia in mice. We found that congenital malformations associated to ALX3-deficiency are enhanced in diabetic pregnancies. Increased expression of genes encoding oxidative stress-scavenging enzymes in embryos from diabetic mothers was blunted in the absence of ALX3, leading to increased oxidative stress. Levels of ALX3 increased in response to glucose, but ALX3 did not activate oxidative stress defence genes directly. Instead, ALX3 stimulated the transcription of Foxo1, a master regulator of oxidative stress-scavenging genes, by binding to a newly identified binding site located in the Foxo1 promoter. Our data identify ALX3 as an important component of the defence mechanisms against the occurrence of developmental malformations during diabetic gestations, stimulating the expression of oxidative stress-scavenging genes in a glucose-dependent manner via Foxo1 activation. Thus, ALX3 deficiency provides a novel molecular mechanism for developmental defects arising from maternal hyperglycaemia.

Developmental mechanisms of stripe patterns in rodents

Mammalian colour patterns are among the most recognizable characteristics found in nature and can have a profound impact on fitness. However, little is known about the mechanisms underlying the formation and subsequent evolution of these patterns. Here we show that, in the African striped mouse (Rhabdomys pumilio), periodic dorsal stripes result from underlying differences in melanocyte maturation, which give rise to spatial variation in hair colour. We identify the transcription factor ALX3 as a regulator of this process. In embryonic dorsal skin, patterned expression of Alx3 precedes pigment stripes and acts to directly repress Mitf, a master regulator of melanocyte differentiation, thereby giving rise to light-coloured hair. Moreover, Alx3 is upregulated in the light stripes of chipmunks, which have independently evolved a similar dorsal pattern. Our results show a previously undescribed mechanism for modulating spatial variation in hair colour and provide insights into how phenotypic novelty evolves.

Understanding Mechanisms of GLI-Mediated Transcription during Craniofacial Development and Disease Using the Ciliopathic Mutant, talpid2

The primary cilium is a ubiquitous, microtubule-based organelle that cells utilize to transduce molecular signals. Ciliopathies are a group of diseases that are caused by a disruption in the structure or function of the primary cilium. Over 30% of all ciliopathies are primarily defined by their craniofacial phenotypes, which typically include midfacial defects, cleft lip/palate, micrognathia, aglossia, and craniosynostosis. The frequency and severity of craniofacial phenotypes in ciliopathies emphasizes the importance of the cilium during development of the craniofacial complex. Molecularly, many ciliopathic mutants, including the avian talpid² (ta²), report pathologically high levels of full-length GLI3 (GLI3_FL), which can go on to function as an activator (GLI_A), and reduced production of truncated GLI3 (GLI3_T), which can go on to function as a repressor (GLI_R). These observations suggest that the craniofacial phenotypes of ciliary mutants like ta² are caused either by excessive activity of the GLI_A or reduced activity of GLI_R. To decipher between these two scenarios, we examined GLI3 occupation at the regulatory regions of target genes and subsequent target gene expression. Using in silico strategies we identified consensus GLI binding regions (GBRs) in the avian genome and confirmed GLI3 binding to the regulatory regions of its targets by chromatin immunoprecipitation (ChIP). In ta² mutants, there was a strikingly low number of GLI3 target genes that had significantly increased expression in facial prominences compared to the control embryo and GLI3 occupancy at GBRs associated with target genes was largely reduced. In vitro DNA binding assays, further supported ChIP results, indicated that the excessive GLI3_FL generated in ta² mutants did not bind to GBRs. In light of these results, we explored the possibility of GLI co-regulator proteins playing a role in regulatory mechanism of GLI-mediated transcription. Taken together our studies suggest that craniofacial ciliopathic phenotypes are produced via reduced GLI_T production, allowing for target gene transcription to be mediated by the combinatorial code of GLI co-regulators.

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

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

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

Exome sequencing revealed a novel splice site variant in the ALX1 gene underlying frontonasal dysplasia

Frontonasal dysplasia (FND) is a heterogeneous group of disorders characterized by hypertelorism, telecanthus, broad nasal root, wide prominent nasal bridge, short and wide nasal ridge, broad columella and smooth philtrum. To date one X-linked and three autosomal recessive forms of FND have been reported in different ethnic groups. We sought to identify the gene responsible for FND in a consanguineous Pakistani family segregating the disorder in autosomal recessive pattern. Genome-wide homozygosity mapping using 250KNsp array revealed five homozygous regions in the selected affected individuals. Exome sequencing found a novel splice acceptor site variant (c.661-1G>C: NM_006982.2) in ALX1. Sanger sequencing confirmed the correct segregation of the pathogenic variant in the whole family. Our study concludes that the splice site variant identified in the ALX1 gene causes mild form of FND.

Glucose-dependent downregulation of glucagon gene expression mediated by selective interactions between ALX3 and PAX6 in mouse alpha cells

AIMS/HYPOTHESIS

The stimulation of glucagon secretion in response to decreased glucose levels has been studied extensively. In contrast, little is known about the regulation of glucagon gene expression in response to fluctuations in glucose concentration. Paired box 6 (PAX6) is a key transcription factor that regulates the glucagon promoter by binding to the G1 and G3 elements. Here, we investigated the role of the transcription factor aristaless-like homeobox 3 (ALX3) as a glucose-dependent modulator of PAX6 activity in alpha cells.

METHODS

Experiments were performed in wild-type or Alx3-deficient islets and alphaTC1 cells. We used chromatin immunoprecipitations and electrophoretic mobility shift assays for DNA binding, immunoprecipitations and pull-down assays for protein interactions, transfected cells for promoter activity, and small interfering RNA and quantitative RT-PCR for gene expression.

RESULTS

Elevated glucose concentration resulted in stimulated expression of Alx3 and decreased glucagon gene expression in wild-type islets. In ALX3-deficient islets, basal glucagon levels were non-responsive to changes in glucose concentration. In basal conditions ALX3 bound to the glucagon promoter at G3, but not at G1. ALX3 could form heterodimers with PAX6 that were permissive for binding to G3 but not to G1. Thus, increasing the levels of ALX3 in response to glucose resulted in the sequestration of PAX6 by ALX3 for binding to G1, thus reducing glucagon promoter activation and glucagon gene expression.

CONCLUSIONS/INTERPRETATION

Glucose-stimulated expression of ALX3 in alpha cells provides a regulatory mechanism for the downregulation of glucagon gene expression by interfering with PAX6-mediated transactivation on the glucagon G1 promoter element.

Using frogs faces to dissect the mechanisms underlying human orofacial defects

In this review I discuss how Xenopus laevis is an effective model to dissect the mechanisms underlying orofacial defects. This species has been particularly useful in studying the understudied structures of the developing face including the embryonic mouth and primary palate. The embryonic mouth is the first opening between the foregut and the environment and is critical for adult mouth development. The final step in embryonic mouth formation is the perforation of a thin layer of tissue covering the digestive tube called the buccopharyngeal membrane. When this tissue does not perforate in humans it can pose serious health risks for the fetus and child. The primary palate forms just dorsal to the embryonic mouth and in non-amniotes it functions as the roof of the adult mouth. Defects in the primary palate result in a median oral cleft that appears similar across the vertebrates. In humans, these median clefts are often severe and surgically difficult to repair. Xenopus has several qualities that make it advantageous for craniofacial research. The free living embryo has an easily accessible face and we have also developed several new tools to analyze the development of the region. Further, Xenopus is readily amenable to chemical screens allowing us to uncover novel gene-environment interactions during orofacial development, as well as to define underlying mechanisms governing such interactions. In conclusion, we are utilizing Xenopus in new and innovative ways to contribute to craniofacial research.

Molecular mechanisms of midfacial developmental defects

The morphogenesis of midfacial processes requires the coordination of a variety of cellular functions of both mesenchymal and epithelial cells to develop complex structures. Any failure or delay in midfacial development as well as any abnormal fusion of the medial and lateral nasal and maxillary prominences will result in developmental defects in the midface with a varying degree of severity, including cleft, hypoplasia, and midline expansion. Despite the advances in human genome sequencing technology, the causes of nearly 70% of all birth defects, which include midfacial development defects, remain unknown. Recent studies in animal models have highlighted the importance of specific signaling cascades and genetic-environmental interactions in the development of the midfacial region. This review will summarize the current understanding of the morphogenetic processes and molecular mechanisms underlying midfacial birth defects based on mouse models with midfacial developmental abnormalities.

Common mechanisms in development and disease: BMP signaling in craniofacial development

BMP signaling is one of the key pathways regulating craniofacial development. It is involved in the early patterning of the head, the development of cranial neural crest cells, and facial patterning. It regulates development of its mineralized structures, such as cranial bones, maxilla, mandible, palate, and teeth. Targeted mutations in the mouse have been instrumental to delineate the functional involvement of this signaling network in different aspects of craniofacial development. Gene polymorphisms and mutations in BMP pathway genes have been associated with various non-syndromic and syndromic human craniofacial malformations. The identification of intricate cellular interactions and underlying molecular pathways illustrate the importance of local fine-regulation of Bmp signaling to control proliferation, apoptosis, epithelial-mesenchymal interactions, and stem/progenitor differentiation during craniofacial development. Thus, BMP signaling contributes both to shape and functionality of our facial features. BMP signaling also regulates postnatal craniofacial growth and is associated with dental structures life-long. A more detailed understanding of BMP function in growth, homeostasis, and repair of postnatal craniofacial tissues will contribute to our ability to rationally manipulate this signaling network in the context of tissue engineering.

A new frontonasal dysplasia syndrome associated with deletion of the SIX2 gene

The frontonasal dysplasias are a group of craniofacial phenotypes characterized by hypertelorism, nasal clefting, frontal bossing, and abnormal hairline. These conditions are caused by recessive mutations in members of the aristaless gene family, resulting in abnormal cranial neural crest migration and differentiation. We report a family with a dominantly inherited craniofacial phenotype comprised of frontal bossing with high hairline, ptosis, hypertelorism, broad nasal tip, large anterior fontanelle, cranial base anomalies, and sagittal synostosis. Chromosomal microarray identified a heterozygous 108.3 kilobase deletion of chromosome 2p21 segregating with phenotype and limited to the sine oculis homeobox gene SIX2 and surrounding noncoding DNA. Similar to the human SIX2 deletion phenotype, one mouse model of frontonasal dysplasia, brachyrrhine, exhibits dominant inheritance and impaired cranial base chondrogenesis associated with reduced Six2 expression. We report the first human autosomal dominant frontonasal dysplasia syndrome associated with SIX2 deletion and with phenotypic similarities to murine models of Six2 Loss-of-function.

A novel allele of Alx4 results in reduced Fgf10 expression and failure of eyelid fusion in mice

Normal fusion of developing eyelids requires coordination of inductive signals from the eyelid mesenchyme with migration of the periderm cell layer and constriction of the eyelids across the eye. Failure of this process results in an eyelids open at birth (EOB) phenotype in mice. We have identified a novel spontaneous allele of Alx4 that displays EOB, in addition to polydactyly and cranial malformations. Alx4 is expressed in the eyelid mesenchyme prior to and during eyelid fusion in a domain overlapping the expression of genes that also play a role in normal eyelid development. We show that Alx4 mutant mice have reduced expression of Fgf10, a key factor expressed in the mesenchyme that is required for initiation of eyelid fusion by the periderm. This is accompanied by a reduced number of periderm cells expressing phosphorylated c-Jun, consistent with the incomplete ablation of Fgf10 expression. Together, these data demonstrate that eyelid fusion in mice requires the expression of Alx4, accompanied by the loss of normal expression of essential components of the eyelid fusion pathway.