A gene expression atlas of early craniofacial development

Genetic architecture of trait variance in craniofacial morphology

The genetic architecture of trait variance has long been of interest in genetics and evolution. One of the earliest attempts to understand this architecture was presented in Lerner's Genetic Homeostasis (1954). Lerner proposed that heterozygotes should be better able to tolerate environmental perturbations because of functional differences between the alleles at a given locus, with each allele optimal for slightly different environments. This greater robustness to environmental variance, he argued, would result in smaller trait variance for heterozygotes. The evidence for Lerner's hypothesis has been inconclusive. To address this question using modern genomic methods, we mapped loci associated with differences in trait variance (vQTL) on 1,101 individuals from the F34 of an advanced intercross between LG/J and SM/J mice. We also mapped epistatic interactions for these vQTL in order to understand the influence of epistasis for the architecture of trait variance. We did not find evidence supporting Lerner's hypothesis, that heterozygotes tend to have smaller trait variances than homozygotes. We further show that the effects of most mapped loci on trait variance are produced by epistasis affecting trait means and that those epistatic effects account for about a half of the differences in genotypic-specific trait variances. Finally, we propose a model where the different interactions between the additive and dominance effects of the vQTL and their epistatic partners can explain Lerner's original observations but can also be extended to include other conditions where heterozygotes are not the least variable genotype.

A. Epigenetic regulation of craniofacial development and disease

BACKGROUND

The formation of the craniofacial complex relies on proper neural crest development. The gene regulatory networks (GRNs) and signaling pathways orchestrating this process have been extensively studied. These GRNs and signaling cascades are tightly regulated as alterations to any stage of neural crest development can lead to common congenital birth defects, including multiple syndromes affecting facial morphology as well as nonsyndromic facial defects, such as cleft lip with or without cleft palate. Epigenetic factors add a hierarchy to the regulation of transcriptional networks and influence the spatiotemporal activation or repression of specific gene regulatory cascades; however less is known about their exact mechanisms in controlling precise gene regulation.

AIMS

In this review, we discuss the role of epigenetic factors during neural crest development, specifically during craniofacial development and how compromised activities of these regulators contribute to congenital defects that affect the craniofacial complex.

Determining zebrafish dorsal organizer size by a negative feedback loop between canonical/non-canonical Wnts and Tlr4/NFκB

In vertebrate embryos, the canonical Wnt ligand primes the formation of dorsal organizers that govern dorsal-ventral patterns by secreting BMP antagonists. In contrast, in Drosophila embryos, Toll-like receptor (Tlr)-mediated NFκB activation initiates dorsal-ventral patterning, wherein Wnt-mediated negative feedback regulation of Tlr/NFκB generates a BMP antagonist-secreting signalling centre to control the dorsal-ventral pattern. Although both Wnt and BMP antagonist are conserved among species, the involvement of Tlr/NFκB and feedback regulation in vertebrate organizer formation remains unclear. By imaging and genetic modification, we reveal that a negative feedback loop between canonical and non-canonical Wnts and Tlr4/NFκB determines the size of zebrafish organizer, and that Tlr/NFκB and Wnts switch initial cue and feedback mediator roles between Drosophila and zebrafish. Here, we show that canonical Wnt signalling stimulates the expression of the non-canonical Wnt5b ligand, activating the Tlr4 receptor to stimulate NFκB-mediated transcription of the Wnt antagonist frzb, restricting Wnt-dependent dorsal organizer formation.

Combined genome-wide association study of 136 quantitative ear morphology traits in multiple populations reveal 8 novel loci

Human ear morphology, a complex anatomical structure represented by a multidimensional set of correlated and heritable phenotypes, has a poorly understood genetic architecture. In this study, we quantitatively assessed 136 ear morphology traits using deep learning analysis of digital face images in 14,921 individuals from five different cohorts in Europe, Asia, and Latin America. Through GWAS meta-analysis and C-GWASs, a recently introduced method to effectively combine GWASs of many traits, we identified 16 genetic loci involved in various ear phenotypes, eight of which have not been previously associated with human ear features. Our findings suggest that ear morphology shares genetic determinants with other surface ectoderm-derived traits such as facial variation, mono eyebrow, and male pattern baldness. Our results enhance the genetic understanding of human ear morphology and shed light on the shared genetic contributors of different surface ectoderm-derived phenotypes. Additionally, gene editing experiments in mice have demonstrated that knocking out the newly ear-associated gene (Intu) and a previously ear-associated gene (Tbx15) causes deviating mouse ear morphology.

Induced pluripotent stem cell technology in bone biology

Technologies on the development and differentiation of human induced pluripotent stem cells (hiPSCs) are rapidly improving, and have been applied to create cell types relevant to the bone field. Differentiation protocols to form bona fide bone-forming cells from iPSCs are available, and can be used to probe details of differentiation and function in depth. When applied to iPSCs bearing disease-causing mutations, the pathogenetic mechanisms of diseases of the skeleton can be elucidated, along with the development of novel therapeutics. These cells can also be used for development of cell therapies for cell and tissue replacement.

Syntaxin 4 is essential for hearing in human and zebrafish

Congenital hearing impairment (HI) is a genetically highly heterogeneous disorder in which prompt recognition and intervention are crucial to optimize outcomes. In this study, we used exome sequencing to investigate a large consanguineous Pakistani family with eight affected individuals showing bilateral severe-to-profound HI. This identified a homozygous splice region variant in STX4 (c.232 + 6T>C), which causes exon skipping and a frameshift, that segregated with HI (two-point logarithm of odds (LOD) score = 5.9). STX4, a member of the syntaxin family, is a component of the SNARE machinery involved in several vesicle transport and recycling pathways. In silico analysis showed that murine orthologue Stx4a is highly and widespread expressed in the developing and adult inner ear. Immunofluorescent imaging revealed localization of STX4A in the cell body, cell membrane and stereocilia of inner and outer hair cells. Furthermore, a morpholino-based knockdown of stx4 in zebrafish showed an abnormal startle response, morphological and developmental defects, and a disrupted mechanotransduction function in neuromast hair cells measured via FM1-43 uptake. Our findings indicate that STX4 dysfunction leads to HI in humans and zebrafish and supports the evolutionary conserved role of STX4 in inner ear development and hair cell functioning.

Identity, lineage and fates of a temporally distinct progenitor population in the embryonic olfactory epithelium

We defined a temporally and transcriptionally divergent precursor cohort in the medial olfactory epithelium (OE) shortly after it differentiates as a distinct tissue at mid-gestation in the mouse. This temporally distinct population of Ascl1+ cells in the dorsomedial OE is segregated from Meis1+/Pax7+ progenitors in the lateral OE, and does not appear to be generated by Pax7+ lateral OE precursors. The medial Ascl1+ precursors do not yield a substantial number of early-generated ORNs. Instead, they first generate additional proliferative precursors as well as a distinct population of frontonasal mesenchymal cells associated with the migratory mass that surrounds the nascent olfactory nerve. Parallel to these in vivo distinctions, isolated medial versus lateral OE precursors in vitro retain distinct proliferative capacities and modes of division that reflect their in vivo identities. At later fetal stages, these early dorsomedial Ascl1+ precursors cells generate spatially restricted subsets of ORNs as well as other non-neuronal cell classes. Accordingly, the initial compliment of ORNs and other OE cell types is derived from at least two distinct early precursor populations: lateral Meis1/Pax7+ precursors that generate primarily early ORNs, and a temporally, spatially, and transcriptionally distinct subset of medial Ascl1+ precursors that initially generate additional OE progenitors and apparent migratory mass cells before yielding a subset of ORNs and likely supporting cell classes.

A novel autosomal dominant GREB1L variant associated with non-syndromic hearing impairment in Ghana

BACKGROUND

Childhood hearing impairment (HI) is genetically heterogeneous with many implicated genes, however, only a few of these genes are reported in African populations.

METHODS

This study used exome and Sanger sequencing to resolve the possible genetic cause of non-syndromic HI in a Ghanaian family.

RESULTS

We identified a novel variant c.3041G > A: p.(Gly1014Glu) in GREB1L (DFNA80) in the index case. The GREB1L: p.(Gly1014Glu) variant had a CADD score of 26.5 and was absent from human genomic databases such as TopMed and gnomAD. In silico homology protein modeling approaches displayed major structural differences between the wildtype and mutant proteins. Additionally, the variant was predicted to probably affect the secondary protein structure that may impact its function. Publicly available expression data shows a higher expression of Greb1L in the inner ear of mice during development and a reduced expression in adulthood, underscoring its importance in the development of the inner ear structures.

CONCLUSION

This report on an African individual supports the association of GREB1L variant with non-syndromic HI and extended the evidence of the implication of GREB1L variants in HI in diverse populations.

The RhoGAP-myosin Myo9b regulates ocular lens pit morphogenesis

BACKGROUND

During eye development the lens placode invaginates to form the lens pit. Further bending of lens epithelium and separation from ectoderm leads eventually to a spherical lens vesicle with enclosed extracellular fluid. Changes in epithelial morphology involve the actin cytoskeleton and its regulators. The myosin Myo9b is simultaneously an actin-based motor and Rho GTPase-activating protein that regulates actin cytoskeleton organization. Myo9b-deficient adult mice and embryos were analyzed for eye malformations and alterations in lens development.

RESULTS

Myo9b-deficient mice showed a high incidence of microphthalmia and cataracts with occasional blepharitis. Formation of the lens vesicle during embryonic lens development was disordered in virtually all embryos. Lens placode invagination was less deep and gave rise to a conical structure instead of a spherical pit. At later stages either no lens vesicle was formed or a significantly smaller one that was not enclosed by the optic cup. Expression of the cell fate marker Pax6 was not altered. Staining of adherens junctions and F-actin was most intense at the tip of conical invaginations, suggesting that mechanical forces are not properly coordinated between epithelial cells that form the pit.

CONCLUSIONS

Myo9b is a critical regulator of ocular lens vesicle morphogenesis during eye development.

Exome sequencing of families from Ghana reveals known and candidate hearing impairment genes

We investigated hearing impairment (HI) in 51 families from Ghana with at least two affected members that were negative for GJB2 pathogenic variants. DNA samples from 184 family members underwent whole-exome sequencing (WES). Variants were found in 14 known non-syndromic HI (NSHI) genes [26/51 (51.0%) families], five genes that can underlie either syndromic HI or NSHI [13/51 (25.5%)], and one syndromic HI gene [1/51 (2.0%)]. Variants in CDH23 and MYO15A contributed the most to HI [31.4% (16/51 families)]. For DSPP, an autosomal recessive mode of inheritance was detected. Post-lingual expression was observed for a family segregating a MARVELD2 variant. To our knowledge, seven novel candidate HI genes were identified (13.7%), with six associated with NSHI (INPP4B, CCDC141, MYO19, DNAH11, POTEI, and SOX9); and one (PAX8) with Waardenburg syndrome. MYO19 and DNAH11 were replicated in unrelated Ghanaian probands. Six of the novel genes were expressed in mouse inner ear. It is known that Pax8-/- mice do not respond to sound, and depletion of Sox9 resulted in defective vestibular structures and abnormal utricle development. Most variants (48/60; 80.0%) have not previously been associated with HI. Identifying seven candidate genes in this study emphasizes the potential of novel HI genes discovery in Africa.

High prevalence of deleterious mutations in concomitant nonsyndromic cleft and outflow tract heart defects

Our previous work demonstrating enrichment of outflow tract (OFT) congenital heart disease (CHD) in children with cleft lip and/or palate (CL/P) suggests derangements in common underlying developmental pathways. The current pilot study examines the underlying genetics of concomitant nonsyndromic CL/P and OFT CHD phenotype. Of 575 patients who underwent CL/P surgery at Children's Hospital Los Angeles, seven with OFT CHD, negative chromosomal microarray analysis, and no recognizable syndromic association were recruited with their parents (as available). Whole genome sequencing of blood samples paired with whole‐blood‐based RNA sequencing for probands was performed. A pathogenic or potentially pathogenic variant was identified in 6/7 (85.7%) probands. A total of seven candidate genes were mutated (CHD7, SMARCA4, MED12, APOB, RNF213, SETX, and JAG1). Gene ontology analysis of variants predicted involvement in binding (100%), regulation of transcription (42.9%), and helicase activity (42.9%). Four patients (57.1%) expressed gene variants (CHD7, SMARCA4, MED12, and RNF213) previously involved in the Wnt signaling pathway. Our pilot analysis of a small cohort of patients with combined CL/P and OFT CHD phenotype suggests a potentially significant prevalence of deleterious mutations. In our cohort, an overrepresentation of mutations in molecules associated with Wnt‐signaling was found. These variants may represent an expanded phenotypic heterogeneity within known monogenic disease genes or provide novel evidence of shared developmental pathways. The mechanistic implications of these mutations and subsequent developmental derangements resulting in the CL/P and OFT CHD phenotype require further analysis in a larger cohort of patients.

Targeted Re-Sequencing of the 2p21 Locus Identifies Non-Syndromic Cleft Lip Only Novel Susceptibility Gene ZFP36L2

rs7590268 present on the 2p21 locus was identified to be associated with non-syndromic cleft lip with or without cleft palate (NSCL/P) in several populations, including the Chinese Han population, indicating that 2p21 was a susceptibility locus for NSCL/P. However, previous studies have only identified common single-nucleotide polymorphism (SNP) within the THADA gene, neglecting the rare variants and other genes in 2p21; thus, this study was designed to investigate additional variants and novel susceptibility genes in 2p21. A total of 159 NSCL/P patients and 542 controls were recruited in the discovery phase, whereas 1830 NSCL/P patients and 2,436 controls were recruited in the replication phase. After targeted region sequencing, we performed association and burden analyses for the common and rare variants, respectively. Furthermore, RNA-seq, proliferation assay and cell cycle analysis were performed to clarify the possible function of the candidate gene ZFP36L2. Association analysis showed that four SNPs were specifically associated with non-syndromic cleft lip only (NSCLO) and two SNPs were associated with both NSCLO and NSCL/P. Burden analysis indicated that ZFP36L2 was associated with NSCLO (p = .0489, OR = 2.41, 95% CI: 0.98–5.90). Moreover, SNPs in the ZFP36L2 targeted gene JUP were also associated with NSCLO. ZFP36L2 also inhibited cell proliferation and induced G2 phase arrest in the GMSM-K cell line. Therefore, we proposed that ZFP36L2 is a novel susceptibility gene of NSCLO in the 2p21 locus, which could lead to NSCLO by modulating cell proliferation and cycle.

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

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

Sperm Transcriptional State Associated with Paternal Transmission of Stress Phenotypes

Paternal stress can induce long-lasting changes in germ cells potentially propagating heritable changes across generations. To date, no studies have investigated differences in transmission patterns between stress-resilient and stress-susceptible mice. We tested the hypothesis that transcriptional alterations in sperm during chronic social defeat stress (CSDS) transmit increased susceptibility to stress phenotypes to the next generation. We demonstrate differences in offspring from stressed fathers that depend on paternal category (resilient vs susceptible) and offspring sex. Importantly, artificial insemination (AI) reveals that sperm mediates some of the behavioral phenotypes seen in offspring. Using RNA-sequencing (RNA-seq), we report substantial and distinct changes in the transcriptomic profiles of sperm following CSDS in susceptible versus resilient fathers, with alterations in long noncoding RNAs (lncRNAs) predominating especially in susceptibility. Correlation analysis revealed that these alterations were accompanied by a loss of regulation of protein-coding genes by lncRNAs in sperm of susceptible males. We also identify several co-expression gene modules that are enriched in differentially expressed genes (DEGs) in sperm from either resilient or susceptible fathers. Taken together, these studies advance our understanding of intergenerational epigenetic transmission of behavioral experience.SIGNIFICANCE STATEMENT This manuscript contributes to the complex factors that influence the paternal transmission of stress phenotypes. By leveraging the segregation of males exposed to chronic social defeat stress (CSDS) into either resilient or susceptible categories we were able to identify the phenotypic differences in the paternal transmission of stress phenotypes across generations between the two lineages. Importantly, this work also alludes to the significance of both long noncoding RNAs (lncRNAs) and protein coding genes (PCGs) mediating the paternal transmission of stress. The knowledge gained from these data are of particular interest in understanding the risk for the development of psychiatric disorders such as anxiety and depression.

Cadherins in early neural development

During early neural development, changes in signalling inform the expression of transcription factors that in turn instruct changes in cell identity. At the same time, switches in adhesion molecule expression result in cellular rearrangements that define the morphology of the emerging neural tube. It is becoming increasingly clear that these two processes influence each other; adhesion molecules do not simply operate downstream of or in parallel with changes in cell identity but rather actively feed into cell fate decisions. Why are differentiation and adhesion so tightly linked? It is now over 60 years since Conrad Waddington noted the remarkable "Constancy of the Wild Type" (Waddington in Nature 183: 1654-1655, 1959) yet we still do not fully understand the mechanisms that make development so reproducible. Conversely, we do not understand why directed differentiation of cells in a dish is sometimes unpredictable and difficult to control. It has long been suggested that cells make decisions as 'local cooperatives' rather than as individuals (Gurdon in Nature 336: 772-774, 1988; Lander in Cell 144: 955-969, 2011). Given that the cadherin family of adhesion molecules can simultaneously influence morphogenesis and signalling, it is tempting to speculate that they may help coordinate cell fate decisions between neighbouring cells in the embryo to ensure fidelity of patterning, and that the uncoupling of these processes in a culture dish might underlie some of the problems with controlling cell fate decisions ex-vivo. Here we review the expression and function of cadherins during early neural development and discuss how and why they might modulate signalling and differentiation as neural tissues are formed.

Crosstalk among Calcium ATPases: PMCA, SERCA and SPCA in Mental Diseases

Calcium in mammalian neurons is essential for developmental processes, neurotransmitter release, apoptosis, and signal transduction. Incorrectly processed Ca²⁺ signal is well-known to trigger a cascade of events leading to altered response to variety of stimuli and persistent accumulation of pathological changes at the molecular level. To counterbalance potentially detrimental consequences of Ca²⁺, neurons are equipped with sophisticated mechanisms that function to keep its concentration in a tightly regulated range. Calcium pumps belonging to the P-type family of ATPases: plasma membrane Ca²⁺-ATPase (PMCA), sarco/endoplasmic Ca²⁺-ATPase (SERCA) and secretory pathway Ca²⁺-ATPase (SPCA) are considered efficient line of defense against abnormal Ca²⁺ rises. However, their role is not limited only to Ca²⁺ transport, as they present tissue-specific functionality and unique sensitive to the regulation by the main calcium signal decoding protein—calmodulin (CaM). Based on the available literature, in this review we analyze the contribution of these three types of Ca²⁺-ATPases to neuropathology, with a special emphasis on mental diseases.

GREB1L, CRELD2 and ITGA10 expression in the human developmental and postnatal kidneys: an immunohistochemical study

BACKGROUND

Aim of our study is to provide an insight into the genetic expression landscape of GREB1L, ITGA10 and CRELD2 which are important in human genitourinary tract development which might help elucidate the critical stages for the onset of kidney anomalies.

METHODS

Morphological parameters were analyzed using immunohistochemistry on human foetal (13-38 w) and postnatal (1.5 and 7.5y) human kidney samples.

RESULTS

GREB1L marker had a strong intensity and the highest rate in proximal tubules (PTC) of 1.5 years' kidney (90.25%). In the distal tubules (DCT) there were statistically significant differences in 13 w, 15 w, 16 w, 21 w, 38 w and 7.5y regarding 1.5y (Kruskal-Wallis test, p < 0.001). There was significantly more GREB1L in the glomeruli at 21 w and 38 w in regard to all other stages (Kruskal-Wallis test, p < 0.01). ITGA10 staining intensity was strongest in PCT with the highest rate in 13 w (92.75%), while the lowest rate was found in glomeruli and DCT (Kruskal-Wallis test, p < 0.001). CRELD2 had the strongest staining intensity in PCT with the highest rate in 13 w and 1.5y (92.25%) and lowest in the glomeruli of 7.5 years (24.3 %). In DCT there were statistically significant differences in CRELD2 positive cells in 13 w, 15 w, 16 w, 21 w, 38 w and 7.5y regarding 1.5y (Kruskal-Wallis test, p < 0.01). ITGA10 and CRELD2 co-localised in the postnatal period in DCT.

CONCLUSION

High kidney expressions of GREB1L, ITGA10 and CRELD2 even in the postnatal period implicate their importance not only for the onset of CAKUT in the case of their mutation but also for maintenance of kidney homeostasis.

An Alternative Splicing Program for Mouse Craniofacial Development

Alternative splicing acts as a fundamental mechanism to increase the number of functional transcripts that can be derived from the genome - and its appropriate regulation is required to direct normal development, differentiation, and physiology, in many species. Recent studies have highlighted that mutation of splicing factors, resulting in the disruption of alternative splicing, can have profound consequences for mammalian craniofacial development. However, there has been no systematic analysis of the dynamics of differential splicing during the critical period of face formation with respect to age, tissue layer, or prominence. Here we used deep RNA sequencing to profile transcripts expressed in the developing mouse face for both ectodermal and mesenchymal tissues from the three facial prominences at critical ages for facial development, embryonic days 10.5, 11.5, and 12.5. We also derived separate expression data from the nasal pit relating to the differentiation of the olfactory epithelium for a total of 60 independent datasets. Analysis of these datasets reveals the differential expression of multiple genes, but we find a similar number of genes are regulated only via differential splicing, indicating that alternative splicing is a major source of transcript diversity during facial development. Notably, splicing changes between tissue layers and over time are more prevalent than between prominences, with exon skipping the most common event. We next examined how the variation in splicing correlated with the expression of RNA binding proteins across the various datasets. Further, we assessed how binding sites for splicing regulatory molecules mapped with respect to intron exon boundaries. Overall these studies help define an alternative splicing regulatory program that has important consequences for facial development.

Lineage-specific differentiation of osteogenic progenitors from pluripotent stem cells reveals the FGF1-RUNX2 association in neural crest-derived osteoprogenitors

Human pluripotent stem cells (hPSCs) can provide a platform to model bone organogenesis and disease. To reflect the developmental process of the human skeleton, hPSC differentiation methods should include osteogenic progenitors (OPs) arising from three distinct embryonic lineages: the paraxial mesoderm, lateral plate mesoderm, and neural crest. Although OP differentiation protocols have been developed, the lineage from which they are derived, as well as characterization of their genetic and molecular differences, has not been well reported. Therefore, to generate lineage-specific OPs from human embryonic stem cells and human induced pluripotent stem cells, we employed stepwise differentiation of paraxial mesoderm-like cells, lateral plate mesoderm-like cells, and neural crest-like cells toward their respective OP subpopulation. Successful differentiation, confirmed through gene expression and in vivo assays, permitted the identification of transcriptomic signatures of all three cell populations. We also report, for the first time, high FGF1 levels in neural crest-derived OPs-a notable finding given the critical role of fibroblast growth factors (FGFs) in osteogenesis and mineral homeostasis. Our results indicate that FGF1 influences RUNX2 levels, with concomitant changes in ERK1/2 signaling. Overall, our study further validates hPSCs' power to model bone development and disease and reveals new, potentially important pathways influencing these processes.

Autosomal Dominantly Inherited GREB1L Variants in Individuals with Profound Sensorineural Hearing Impairment

Congenital hearing impairment is a sensory disorder that is genetically highly heterogeneous. By performing exome sequencing in two families with congenital nonsyndromic profound sensorineural hearing loss (SNHL), we identified autosomal dominantly inherited missense variants [p.(Asn283Ser); p.(Thr116Ile)] in GREB1L, a neural crest regulatory molecule. The p.(Thr116Ile) variant was also associated with bilateral cochlear aplasia and cochlear nerve aplasia upon temporal bone imaging, an ultra-rare phenotype previously seen in patients with de novo GREB1L variants. An important role of GREB1L in normal ear development has also been demonstrated by greb1l^-/- zebrafish, which show an abnormal sensory epithelia innervation. Last, we performed a review of all disease-associated variation described in GREB1L, as it has also been implicated in renal, bladder and genital malformations. We show that the spectrum of features associated with GREB1L is broad, variable and with a high level of reduced penetrance, which is typically characteristic of neurocristopathies. So far, seven GREB1L variants (14%) have been associated with ear-related abnormalities. In conclusion, these results show that autosomal dominantly inherited variants in GREB1L cause profound SNHL. Furthermore, we provide an overview of the phenotypic spectrum associated with GREB1L variants and strengthen the evidence of the involvement of GREB1L in human hearing.

SpatialDB: a database for spatially resolved transcriptomes

Spatially resolved transcriptomic techniques allow the characterization of spatial organization of cells in tissues, which revolutionize the studies of tissue function and disease pathology. New strategies for detecting spatial gene expression patterns are emerging, and spatially resolved transcriptomic data are accumulating rapidly. However, it is not convenient for biologists to exploit these data due to the diversity of strategies and complexity in data analysis. Here, we present SpatialDB, the first manually curated database for spatially resolved transcriptomic techniques and datasets. The current version of SpatialDB contains 24 datasets (305 sub-datasets) from 5 species generated by 8 spatially resolved transcriptomic techniques. SpatialDB provides a user-friendly web interface for visualization and comparison of spatially resolved transcriptomic data. To further explore these data, SpatialDB also provides spatially variable genes and their functional enrichment annotation. SpatialDB offers a repository for research community to investigate the spatial cellular structure of tissues, and may bring new insights into understanding the cellular microenvironment in disease. SpatialDB is freely available at https://www.spatialomics.org/SpatialDB.

An internal deletion of ADAR rescued by MAVS deficiency leads to a minute phenotype

The RNA-editing protein ADAR is essential for early development in the mouse. Genetic evidence suggests that A to I editing marks endogenous RNAs as ‘self’. Today, different Adar knockout alleles have been generated that show a common phenotype of apoptosis, liver disintegration, elevated immune response and lethality at E12.5. All the Adar knockout alleles can be rescued by a concomitant deletion of the innate immunity genes Mavs or Ifih1 (MDA5), albeit to different extents. This suggests multiple functions of ADAR. We analyze Adar^Δ7-9 mice that show a unique growth defect phenotype when rescued by Mavs. We show that Adar^Δ7-9 can form a truncated, unstable, editing deficient protein that is mislocalized. Histological and hematologic analysis of these mice indicate multiple tissue- and hematopoietic defects. Gene expression profiling shows dysregulation of Rps3a1 and Rps3a3 in rescued Adar^Δ7-9. Consistently, a distortion in 40S and 60S ribosome ratios is observed in liver cells. This dysregulation is also seen in Adar^Δ2-13; Mavs^−/− but not in Adar^E861A/E861A; Ifih1^−/− mice, suggesting editing-independent functions of ADAR in regulating expression levels of Rps3a1 and Rps3a3. In conclusion, our study demonstrates the importance of ADAR in post-natal development which cannot be compensated by ADARB1.

Noncoding regions underpin avian bill shape diversification at macroevolutionary scales

Recent progress has been made in identifying genomic regions implicated in trait evolution on a microevolutionary scale in many species, but whether these are relevant over macroevolutionary time remains unclear. Here, we directly address this fundamental question using bird beak shape, a key evolutionary innovation linked to patterns of resource use, divergence, and speciation, as a model trait. We integrate class-wide geometric-morphometric analyses with evolutionary sequence analyses of 10,322 protein-coding genes as well as 229,001 genomic regions spanning 72 species. We identify 1434 protein-coding genes and 39,806 noncoding regions for which molecular rates were significantly related to rates of bill shape evolution. We show that homologs of the identified protein-coding genes as well as genes in close proximity to the identified noncoding regions are involved in craniofacial embryo development in mammals. They are associated with embryonic stem cell pathways, including BMP and Wnt signaling, both of which have repeatedly been implicated in the morphological development of avian beaks. This suggests that identifying genotype-phenotype association on a genome-wide scale over macroevolutionary time is feasible. Although the coding and noncoding gene sets are associated with similar pathways, the actual genes are highly distinct, with significantly reduced overlap between them and bill-related phenotype associations specific to noncoding loci. Evidence for signatures of recent diversifying selection on our identified noncoding loci in Darwin finch populations further suggests that regulatory rather than coding changes are major drivers of morphological diversification over macroevolutionary times.

A Novel Missense Variant in the ALX4 Gene Underlies Mild to Severe Frontonasal Dysplasia in a Consanguineous Family

Background: Frontonasal dysplasia (FND) is a rare developmental disorder characterized by mild to severe changes in skull and brain structures. It is a phenotypically variable and heterogeneous disorder. This study was designed to provide a clinical and genetic analysis of FND in a consanguineous family of Pakistani origin. Methodology and Results: Affected individuals in the family showed characteristic features of frontonasal dysplasia type-2 (FND2), such as nasal bone hypoplasia, hypertelorism, and alopecia. Skull and brain imaging of affected members revealed ossification defects and various types of brain structural anomalies that created a split-brain. Sanger sequencing of the ALX4 gene revealed a homozygous missense variant [NM_021926.4: c.652C>T; p.(Arg218Trp)] in three affected members who demonstrated severe craniofacial anomalies. Heterozygous carriers in the family showed mild FND2 phenotypes. Conclusion: Clinical and genetic analysis of a family, exhibiting FND2 phenotypes, revealed several previously unreported clinical features and a novel missense variant in the ALX4 gene. These results will facilitate diagnosis and genetic counseling of the FND patients in the Pakistani population.

Investigating the gene expression profiles of cells in seven embryonic stages with machine learning algorithms

The development of embryonic cells involves several continuous stages, and some genes are related to embryogenesis. To date, few studies have systematically investigated changes in gene expression profiles during mammalian embryogenesis. In this study, a computational analysis using machine learning algorithms was performed on the gene expression profiles of mouse embryonic cells at seven stages. First, the profiles were analyzed through a powerful Monte Carlo feature selection method for the generation of a feature list. Second, increment feature selection was applied on the list by incorporating two classification algorithms: support vector machine (SVM) and repeated incremental pruning to produce error reduction (RIPPER). Through SVM, we extracted several latent gene biomarkers, indicating the stages of embryonic cells, and constructed an optimal SVM classifier that produced a nearly perfect classification of embryonic cells. Furthermore, some interesting rules were accessed by the RIPPER algorithm, suggesting different expression patterns for different stages.

Generation of Retinaldehyde for Retinoic Acid Biosynthesis

The concentration of all-trans-retinoic acid, the bioactive derivative of vitamin A, is critically important for the optimal performance of numerous physiological processes. Either too little or too much of retinoic acid in developing or adult tissues is equally harmful. All-trans-retinoic acid is produced by the irreversible oxidation of all-trans-retinaldehyde. Thus, the concentration of retinaldehyde as the immediate precursor of retinoic acid has to be tightly controlled. However, the enzymes that produce all-trans-retinaldehyde for retinoic acid biosynthesis and the mechanisms responsible for the control of retinaldehyde levels have not yet been fully defined. The goal of this review is to summarize the current state of knowledge regarding the identities of physiologically relevant retinol dehydrogenases, their enzymatic properties, and tissue distribution, and to discuss potential mechanisms for the regulation of the flux from retinol to retinaldehyde.

High-Resolution Epigenomic Atlas of Human Embryonic Craniofacial Development

Defects in patterning during human embryonic development frequently result in craniofacial abnormalities. The gene regulatory programs that build the craniofacial complex are likely controlled by information located between genes and within intronic sequences. However, systematic identification of regulatory sequences important for forming the human face has not been performed. Here, we describe comprehensive epigenomic annotations from human embryonic craniofacial tissues and systematic comparisons with multiple tissues and cell types. We identified thousands of tissue-specific craniofacial regulatory sequences and likely causal regions for rare craniofacial abnormalities. We demonstrate significant enrichment of common variants associated with orofacial clefting in enhancers active early in embryonic development, while those associated with normal facial variation are enriched near the end of the embryonic period. These data are provided in easily accessible formats for both craniofacial researchers and clinicians to aid future experimental design and interpretation of noncoding variation in those affected by craniofacial abnormalities.

The molecular anatomy of mammalian upper lip and primary palate fusion at single cell resolution

The mammalian lip and primary palate form when coordinated growth and morphogenesis bring the nasal and maxillary processes into contact, and the epithelia co-mingle, remodel and clear from the fusion site to allow mesenchyme continuity. Although several genes required for fusion have been identified, an integrated molecular and cellular description of the overall process is lacking. Here, we employ single cell RNA sequencing of the developing mouse face to identify ectodermal, mesenchymal and endothelial populations associated with patterning and fusion of the facial prominences. This analysis indicates that key cell populations at the fusion site exist within the periderm, basal epithelial cells and adjacent mesenchyme. We describe the expression profiles that make each population unique, and the signals that potentially integrate their behaviour. Overall, these data provide a comprehensive high-resolution description of the various cell populations participating in fusion of the lip and primary palate, as well as formation of the nasolacrimal groove, and they furnish a powerful resource for those investigating the molecular genetics of facial development and facial clefting that can be mined for crucial mechanistic information concerning this prevalent human birth defect.

Integrating gene regulatory pathways into differential network analysis of gene expression data

The advent of next-generation sequencing has introduced new opportunities in analyzing gene expression data. Research in systems biology has taken advantage of these opportunities by gleaning insights into gene regulatory networks through the analysis of gene association networks. Contrasting networks from different populations can reveal the many different roles genes fill, which can lead to new discoveries in gene function. Pathologies can also arise from aberrations in these gene-gene interactions. Exposing these network irregularities provides a new avenue for understanding and treating diseases. A general framework for integrating known gene regulatory pathways into a differential network analysis between two populations is proposed. The framework importantly allows for any gene-gene association measure to be used, and inference is carried out through permutation testing. A simulation study investigates the performance in identifying differentially connected genes when incorporating known pathways, even if the pathway knowledge is partially inaccurate. Another simulation study compares the general framework with four state-of-the-art methods. Two RNA-seq datasets are analyzed to illustrate the use of this framework in practice. In both examples, the analysis reveals genes and pathways that are known to be biologically significant along with potentially novel findings that may be used to motivate future research.

Laser Capture Microdissection of Murine Embryonic Neural Crest Cells

The purpose of this chapter is to provide a step-by-step protocol to enable performance of laser capture microdissection (LCM) on tissue sections from mammalian embryos or postnatal organism stages in order to collect pure populations of neural crest cells from which sufficient amounts of nucleic acids and/or protein can be obtained for quantitative analysis. The methods (1) define a strategy to genetically and indelibly label mammalian neural crest-derived cells with a fluorescent marker, thus enabling their isolation throughout the pre- and postnatal life span of the organism, and (2) describe subsequent isolation by LCM of the labeled neural crest cells, or their derivatives, from embryonic/postnatal tissue cryosections. Details are provided for using the Arcturus PixCell^®IIe Laser Capture Microdissection System (Arcturus) and CapSure LCM Caps (Thermo Fisher Scientific), to which the selected cells adhere upon laser-mediated capture. The protocol outlined herein can be applied in any situation wherein limited cellular samples are available for isolation by LCM. Nucleic acids or proteins can be extracted from LCM-isolated cells and processed for high-density gene expression profiling analyses (microarrays or RNA sequencing), Real-Time PCR (q-PCR) for specific candidate gene expression, investigation of DNA methylation, as well as for varied protein analyses.

Transcriptomics reveals complex kinetics of dorsal-ventral patterning gene expression in the mandibular arch

The mandibular or first pharyngeal arch forms the upper and lower jaws in all gnathostomes. A gene regulatory network that defines ventral, intermediate, and dorsal domains along the dorsal-ventral (D-V) axis of the arch has emerged from studies in zebrafish and mice, but the temporal dynamics of this process remain unclear. To define cell fate trajectories in the arches we have performed quantitative gene expression analyses of D-V patterning genes in pharyngeal arch primordia in zebrafish and mice. Using NanoString technology to measure transcript numbers per cell directly we show that, in many cases, genes expressed in similar D-V domains and induced by similar signals vary dramatically in their temporal profiles. This suggests that cellular responses to D-V patterning signals are likely shaped by the baseline kinetics of target gene expression. Furthermore, similarities in the temporal dynamics of genes that occupy distinct pathways suggest novel shared modes of regulation. Incorporating these gene expression kinetics into our computational models for the mandibular arch improves the accuracy of patterning, and facilitates temporal comparisons between species. These data suggest that the magnitude and timing of target gene expression help diversify responses to patterning signals during craniofacial development.

Secretory pathway calcium ATPase 1 (SPCA1) controls mouse neural tube closure by regulating cytoskeletal dynamics

Neural tube closure relies on the apical constriction of neuroepithelial cells. Research in frog and fly embryos has found links between the levels of intracellular calcium, actomyosin dynamics and apical constriction. However, genetic evidence for a role of calcium in apical constriction during mammalian neurulation is still lacking. Secretory pathway calcium ATPase (SPCA1) regulates calcium homeostasis by pumping cytosolic calcium into the Golgi apparatus. Loss of function in Spca1 causes cranial exencephaly and spinal cord defects in mice, phenotypes previously ascribed to apoptosis. However, our characterization of a novel allele of Spca1 revealed that neurulation defects in Spca1 mutants are not due to cell death, but rather to a failure of neuroepithelial cells to apically constrict. We show that SPCA1 influences cell contractility by regulating myosin II localization. Furthermore, we found that loss of Spca1 disrupts actin dynamics and the localization of the actin remodeling protein cofilin 1. Taken together, our results provide evidence that SPCA1 promotes neurulation by regulating the cytoskeletal dynamics that promote apical constriction and identify cofilin 1 as a downstream effector of SPCA1 function.

De novo variants in GREB1L are associated with non-syndromic inner ear malformations and deafness

Congenital inner ear malformations affecting both the osseous and membranous labyrinth can have a devastating impact on hearing and language development. With the exception of an enlarged vestibular aqueduct, non-syndromic inner ear malformations are rare, and their underlying molecular biology has thus far remained understudied. To identify molecular factors that might be important in the developing inner ear, we adopted a family-based trio exome sequencing approach in young unrelated subjects with severe inner ear malformations. We identified two previously unreported de novo loss-of-function variants in GREB1L [c.4368G>T;p.(Glu1410fs) and c.982C>T;p.(Arg328*)] in two affected subjects with absent cochleae and eighth cranial nerve malformations. The cochlear aplasia in these affected subjects suggests that a developmental arrest or problem at a very early stage of inner ear development exists, e.g., during the otic pit formation. Craniofacial Greb1l RNA expression peaks in mice during this time frame (E8.5). It also peaks in the developing inner ear during E13-E16, after which it decreases in adulthood. The crucial function of Greb1l in craniofacial development is also evidenced in knockout mice, which develop severe craniofacial abnormalities. In addition, we show that Greb1l-/- zebrafish exhibit a loss of abnormal sensory epithelia innervation. An important role for Greb1l in sensory epithelia innervation development is supported by the eighth cranial nerve deficiencies seen in both affected subjects. In conclusion, we demonstrate that GREB1L is a key player in early inner ear and eighth cranial nerve development. Abnormalities in cochleovestibular anatomy can provide challenges for cochlear implantation. Combining a molecular diagnosis with imaging techniques might aid the development of individually tailored therapeutic interventions in the future.

Neural crest and the origin of species-specific pattern

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

Developmental constraint through negative pleiotropy in the zygomatic arch

BACKGROUND

Previous analysis suggested that the relative contribution of individual bones to regional skull lengths differ between inbred mouse strains. If the negative correlation of adjacent bone lengths is associated with genetic variation in a heterogeneous population, it would be an example of negative pleiotropy, which occurs when a genetic factor leads to opposite effects in two phenotypes. Confirming negative pleiotropy and determining its basis may reveal important information about the maintenance of overall skull integration and developmental constraint on skull morphology.

RESULTS

We identified negative correlations between the lengths of the frontal and parietal bones in the midline cranial vault as well as the zygomatic bone and zygomatic process of the maxilla, which contribute to the zygomatic arch. Through gene association mapping of a large heterogeneous population of Diversity Outbred (DO) mice, we identified a quantitative trait locus on chromosome 17 driving the antagonistic contribution of these two zygomatic arch bones to total zygomatic arch length. Candidate genes in this region were identified and real-time PCR of the maxillary processes of DO founder strain embryos indicated differences in the RNA expression levels for two of the candidate genes, Camkmt and Six2.

CONCLUSIONS

A genomic region underlying negative pleiotropy of two zygomatic arch bones was identified, which provides a mechanism for antagonism in component bone lengths while constraining overall zygomatic arch length. This type of mechanism may have led to variation in the contribution of individual bones to the zygomatic arch noted across mammals. Given that similar genetic and developmental mechanisms may underlie negative correlations in other parts of the skull, these results provide an important step toward understanding the developmental basis of evolutionary variation and constraint in skull morphology.

Inhibition of MAPK-Erk pathway in vivo attenuates aortic valve disease processes in Emilin1-deficient mouse model

Aortic valve disease (AVD) is a common condition with a progressive natural history, and presently, there are no pharmacologic treatment strategies. Elastic fiber fragmentation (EFF) is a hallmark of AVD, and increasing evidence implicates developmental elastic fiber assembly defects. Emilin1 is a glycoprotein necessary for elastic fiber assembly that is present in both developing and mature human and mouse aortic valves. The Emilin1‐deficient mouse (Emilin1^−/−) is a model of latent AVD, characterized by activated TGFβ/MEK/p‐Erk signaling and upregulated elastase activity. Emilin1^−/− aortic valves demonstrate early EFF and aberrant angiogenesis followed by late neovascularization and fibrosis. The objective of this study was to test the effectiveness of three different targeted therapies. Aged (12–14 months) Emilin1^−/− mice were treated with refametinib (RDEA‐119, MEK1/2 inhibitor), doxycycline (elastase inhibitor), or G6‐31 (anti‐VEGF‐A mouse antibody) for 4 weeks. Refametinib‐ and doxycycline‐treated Emilin1^−/− mice markedly reduced MEK/p‐Erk activation in valve tissue. Furthermore, both refametinib and doxycycline attenuated elastolytic cathepsin K, L, MMP‐2, and MMP‐9 activation, and abrogated macrophage and neutrophil infiltration in Emilin1^−/− aortic valves. RNAseq analysis was performed in aortic valve tissue from adult (4 months) and aged (14 months) Emilin1^−/− and age‐matched wild‐type control mice, and demonstrated upregulation of genes associated with MAPK/MEK/p‐Erk signaling and elastases at the adult stage and inflammatory pathways at the aged stage controlling for age. These results suggest that Erk1/2 signaling is an important modulator of early elastase activation, and pharmacological inhibition using refametinib may be a promising treatment to halt AVD progression

Genome-wide analysis of facial skeletal regionalization in zebrafish

Patterning of the facial skeleton involves the precise deployment of thousands of genes in distinct regions of the pharyngeal arches. Despite the significance for craniofacial development, how genetic programs drive this regionalization remains incompletely understood. Here we use combinatorial labeling of zebrafish cranial neural crest-derived cells (CNCCs) to define global gene expression along the dorsoventral axis of the developing arches. Intersection of region-specific transcriptomes with expression changes in response to signaling perturbations demonstrates complex roles for Endothelin 1 (Edn1) signaling in the intermediate joint-forming region, yet a surprisingly minor role in ventralmost regions. Analysis of co-variance across multiple sequencing experiments further reveals clusters of co-regulated genes, with in situ hybridization confirming the domain-specific expression of novel genes. We then created loss-of-function alleles for 12 genes and uncovered antagonistic functions of two new Edn1 targets, follistatin a (fsta) and emx2, in regulating cartilaginous joints in the hyoid arch. Our unbiased discovery and functional analysis of genes with regional expression in zebrafish arch CNCCs reveals complex regulation by Edn1 and points to novel candidates for craniofacial disorders.

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.

TFAP2A is a component of the ZEB1/2 network that regulates TGFB1-induced epithelial to mesenchymal transition

Background

The transition between epithelial and mesenchymal phenotypes (EMT) occurs in a variety of contexts. It is critical for mammalian development and it is also involved in tumor initiation and progression. Master transcription factor (TF) regulators of this process are conserved between mouse and human.

Methods

From a computational analysis of a variety of high-throughput sequencing data sets we initially inferred that TFAP2A is connected to the core EMT network in both species. We then analysed publicly available human breast cancer data for TFAP2A expression and also studied the expression (by mRNA sequencing), activity (by monitoring the expression of its predicted targets), and binding (by electrophoretic mobility shift assay and chromatin immunoprecipitation) of this factor in a mouse mammary gland EMT model system (NMuMG) cell line.

Results

We found that upon induction of EMT, the activity of TFAP2A, reflected in the expression level of its predicted targets, is up-regulated in a variety of systems, both murine and human, while TFAP2A’s expression is increased in more “stem-like” cancers. We provide strong evidence for the direct interaction between the TFAP2A TF and the ZEB2 promoter and we demonstrate that this interaction affects ZEB2 expression. Overexpression of TFAP2A from an exogenous construct perturbs EMT, however, in a manner similar to the downregulation of endogenous TFAP2A that takes place during EMT.

Conclusions

Our study reveals that TFAP2A is a conserved component of the core network that regulates EMT, acting as a repressor of many genes, including ZEB2.

Reviewers

This article has been reviewed by Dr. Martijn Huynen and Dr. Nicola Aceto.

Electronic supplementary material

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

Regulation of pituitary stem cells by epithelial to mesenchymal transition events and signaling pathways

The anterior pituitary gland is comprised of specialized cell-types that produce and secrete polypeptide hormones in response to hypothalamic input and feedback from target organs. These specialized cells arise from stem cells that express SOX2 and the pituitary transcription factor PROP1, which is necessary to establish the stem cell pool and promote an epithelial to mesenchymal-like transition, releasing progenitors from the niche. The adult anterior pituitary responds to physiological challenge by mobilizing the SOX2-expressing progenitor pool and producing additional hormone-producing cells. Knowledge of the role of signaling pathways and extracellular matrix components in these processes may lead to improvements in the efficiency of differentiation of embryonic stem cells or induced pluripotent stem cells into hormone producing cells in vitro. Advances in our basic understanding of pituitary stem cell regulation and differentiation may lead to improved diagnosis and treatment for patients with hypopituitarism.

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.

Implications of Vertebrate Craniodental Evo-Devo for Human Oral Health

Highly processed diets eaten by postindustrial modern human populations coincide with higher frequencies of third molar impaction, malocclusion, and temporomandibular joint disorders that affect millions of people worldwide each year. Current treatments address symptoms, not causes, because the multifactorial etiologies of these three concerns mask which factors incline certain people to malocclusion, impaction, and/or joint issues. Deep scientific curiosity about the origins of jaws and dentitions continues to yield rich insights about the developmental genetic mechanisms that underpin healthy craniodental morphogenesis and integration. Mounting evidence from evolution and development (Evo-Devo) studies suggests that function is another mechanism important to healthy craniodental integration and fit. Starting as early as weaning, softer diets and thus lower bite forces appear to relax or disrupt integration of oral tissues, alter development and growth, and catalyze impaction, malocclusion, and jaw joint disorders. How developing oral tissues respond to bite forces remains poorly understood, but biomechanical feedback seems to alter balances of local bone resorption and deposition at the tooth-bone interface as well as affect tempos and amounts of facial outgrowth. Also, behavioral changes in jaw function and parafunction contribute to degeneration and pain in joint articular cartilages and masticatory muscles. The developmental genetic contribution to craniodental misfits and disorders is undeniable but still unclear; however, at present, human diet and jaw function remain important and much more actionable clinical targets. New Evo-Devo studies are needed to explain how function interfaces with craniodental phenotypic plasticity, variation, and evolvability to yield a spectrum of healthy and mismatched dentitions and jaws.

An integrated genomic approach for the study of mandibular prognathism in the European seabass (Dicentrarchus labrax)

Skeletal anomalies in farmed fish are a relevant issue affecting animal welfare and health and causing significant economic losses. Here, a high-density genetic map of European seabass for QTL mapping of jaw deformity was constructed and a genome-wide association study (GWAS) was carried out on a total of 298 juveniles, 148 of which belonged to four full-sib families. Out of 298 fish, 107 were affected by mandibular prognathism (MP). Three significant QTLs and two candidate SNPs associated with MP were identified. The two GWAS candidate markers were located on ChrX and Chr17, both in close proximity with the peaks of the two most significant QTLs. Notably, the SNP marker on Chr17 was positioned within the Sobp gene coding region, which plays a pivotal role in craniofacial development. The analysis of differentially expressed genes in jaw-deformed animals highlighted the “nervous system development” as a crucial pathway in MP. In particular, Zic2, a key gene for craniofacial morphogenesis in model species, was significantly down-regulated in MP-affected animals. Gene expression data revealed also a significant down-regulation of Sobp in deformed larvae. Our analyses, integrating transcriptomic and GWA methods, provide evidence for putative mechanisms underlying seabass jaw deformity.

Genetics of Combined Pituitary Hormone Deficiency: Roadmap into the Genome Era

The genetic basis for combined pituitary hormone deficiency (CPHD) is complex, involving 30 genes in a variety of syndromic and nonsyndromic presentations. Molecular diagnosis of this disorder is valuable for predicting disease progression, avoiding unnecessary surgery, and family planning. We expect that the application of high throughput sequencing will uncover additional contributing genes and eventually become a valuable tool for molecular diagnosis. For example, in the last 3 years, six new genes have been implicated in CPHD using whole-exome sequencing. In this review, we present a historical perspective on gene discovery for CPHD and predict approaches that may facilitate future gene identification projects conducted by clinicians and basic scientists. Guidelines for systematic reporting of genetic variants and assigning causality are emerging. We apply these guidelines retrospectively to reports of the genetic basis of CPHD and summarize modes of inheritance and penetrance for each of the known genes. In recent years, there have been great improvements in databases of genetic information for diverse populations. Some issues remain that make molecular diagnosis challenging in some cases. These include the inherent genetic complexity of this disorder, technical challenges like uneven coverage, differing results from variant calling and interpretation pipelines, the number of tolerated genetic alterations, and imperfect methods for predicting pathogenicity. We discuss approaches for future research in the genetics of CPHD.

MORN5 Expression during Craniofacial Development and Its Interaction with the BMP and TGFβ Pathways

MORN5 (MORN repeat containing 5) is encoded by a locus positioned on chromosome 17 in the chicken genome. The MORN motif is found in multiple copies in several proteins including junctophilins or phosphatidylinositol phosphate kinase family and the MORN proteins themselves are found across the animal and plant kingdoms. MORN5 protein has a characteristic punctate pattern in the cytoplasm in immunofluorescence imaging. Previously, MORN5 was found among differentially expressed genes in a microarray profiling experiment of the chicken embryo head. Here, we provided in situ hybridization to analyse, in detail, the MORN5 expression in chick craniofacial structures. The expression of MORN5 was first observed at stage HH17-18 (E2.5). MORN5 expression gradually appeared on either side of the primitive oral cavity, within the maxillary region. At stage HH20 (E3), prominent expression was localized in the mandibular prominences lateral to the midline. From stage HH20 up to HH29 (E6), there was strong expression in restricted regions of the maxillary and mandibular prominences. The frontonasal mass (in the midline of the face) expressed MORN5, starting at HH27 (E5). The expression was concentrated in the corners or globular processes, which will ultimately fuse with the cranial edges of the maxillary prominences. MORN5 expression was maintained in the fusion zone up to stage HH29. In sections MORN5 expression was localized preferentially in the mesenchyme. Previously, we examined signals that regulate MORN5 expression in the face based on a previous microarray study. Here, we validated the array results with in situ hybridization and QPCR. MORN5 was downregulated 24 h after Noggin and/or RA treatment. We also determined that BMP pathway genes are downstream of MORN5 following siRNA knockdown. Based on these results, we conclude that MORN5 is both regulated by and required for BMP signaling. The restricted expression of MORN5 in the lip fusion zone shown here supports the human genetic data in which MORN5 variants were associated with increased risk of non-syndromic cleft lip with or without cleft palate.

Remission for Loss of Odontogenic Potential in a New Micromilieu In Vitro

During embryonic organogenesis, the odontogenic potential resides in dental mesenchyme from the bud stage until birth. Mouse dental mesenchymal cells (mDMCs) isolated from the inductive dental mesenchyme of developing molars are frequently used in the context of tooth development and regeneration. We wondered if and how the odontogenic potential could be retained when mDMCs were cultured in vitro. In the present study, we undertook to test the odontogenic potential of cultured mDMCs and attempted to maintain the potential during culturing. We found that cultured mDMCs could retain the odontogenic potential for 24 h with a ratio of 60% for tooth formation, but mDMCs were incapable of supporting tooth formation after more than 24 h in culture. This loss of odontogenic potential was accompanied by widespread transcriptomic alteration and, specifically, the downregulation of some dental mesenchyme-specific genes, such as Pax9, Msx1, and Pdgfrα. To prolong the odontogenic potential of mDMCs in vitro, we then cultured mDMCs in a serum-free medium with Knockout Serum Replacement (KSR) and growth factors (fibroblastic growth factor 2 and epidermal growth factor). In this new micromilieu, mDMCs could maintain the odontogenic potential for 48 h with tooth formation ratio of 50%. Moreover, mDMCs cultured in KSR-supplemented medium gave rise to tooth-like structures when recombined with non-dental second-arch epithelium. Among the supplements, KSR is essential for the survival and adhesion of mDMCs, and both Egf and Fgf2 induced the expression of certain dental mesenchyme-related genes. Taken together, our results demonstrated that the transcriptomic changes responded to the alteration of odontogenic potential in cultured mDMCs and a new micromilieu partly retained this potential in vitro, providing insight into the long-term maintenance of odontogenic potential in mDMCs.

The old and new face of craniofacial research: How animal models inform human craniofacial genetic and clinical data

The craniofacial skeletal structures that comprise the human head develop from multiple tissues that converge to form the bones and cartilage of the face. Because of their complex development and morphogenesis, many human birth defects arise due to disruptions in these cellular populations. Thus, determining how these structures normally develop is vital if we are to gain a deeper understanding of craniofacial birth defects and devise treatment and prevention options. In this review, we will focus on how animal model systems have been used historically and in an ongoing context to enhance our understanding of human craniofacial development. We do this by first highlighting "animal to man" approaches; that is, how animal models are being utilized to understand fundamental mechanisms of craniofacial development. We discuss emerging technologies, including high throughput sequencing and genome editing, and new animal repository resources, and how their application can revolutionize the future of animal models in craniofacial research. Secondly, we highlight "man to animal" approaches, including the current use of animal models to test the function of candidate human disease variants. Specifically, we outline a common workflow deployed after discovery of a potentially disease causing variant based on a select set of recent examples in which human mutations are investigated in vivo using animal models. Collectively, these topics will provide a pipeline for the use of animal models in understanding human craniofacial development and disease for clinical geneticist and basic researchers alike.