Dynamic enhancer landscapes in human craniofacial development.

The genetic basis of human facial variation and craniofacial birth defects remains poorly understood. Distant-acting transcriptional enhancers control the fine-tuned spatiotemporal expression of genes during critical stages of craniofacial development. However, a lack of accurate maps of the genomic locations and cell type-resolved activities of craniofacial enhancers prevents their systematic exploration in human genetics studies. Here, we combine histone modification, chromatin accessibility, and gene expression profiling of human craniofacial development with single-cell analyses of the developing mouse face to define the regulatory landscape of facial development at tissue- and single cell-resolution. We provide temporal activity profiles for 14,000 human developmental craniofacial enhancers. We find that 56% of human craniofacial enhancers share chromatin accessibility in the mouse and we provide cell population- and embryonic stage-resolved predictions of their in vivo activity. Taken together, our data provide an expansive resource for genetic and developmental studies of human craniofacial development.

Single-cell, whole-embryo phenotyping of mammalian developmental disorders.

Mouse models are a critical tool for studying human diseases, particularly developmental disorders1. However, conventional approaches for phenotyping may fail to detect subtle defects throughout the developing mouse2. Here we set out to establish single-cell RNA sequencing of the whole embryo as a scalable platform for the systematic phenotyping of mouse genetic models. We applied combinatorial indexing-based single-cell RNA sequencing3 to profile 101 embryos of 22 mutant and 4 wild-type genotypes at embryonic day 13.5, altogether profiling more than 1.6 million nuclei. The 22 mutants represent a range of anticipated phenotypic severities, from established multisystem disorders to deletions of individual regulatory regions4,5. We developed and applied several analytical frameworks for detecting differences in composition and/or gene expression across 52 cell types or trajectories. Some mutants exhibit changes in dozens of trajectories whereas others exhibit changes in only a few cell types. We also identify differences between widely used wild-type strains, compare phenotyping of gain- versus loss-of-function mutants and characterize deletions of topological associating domain boundaries. Notably, some changes are shared among mutants, suggesting that developmental pleiotropy might be 'decomposable' through further scaling of this approach. Overall, our findings show how single-cell profiling of whole embryos can enable the systematic molecular and cellular phenotypic characterization of mouse mutants with unprecedented breadth and resolution.

Cell Type- and Tissue-specific Enhancers in Craniofacial Development.

The genetic basis of craniofacial birth defects and general variation in human facial shape remains poorly understood. Distant-acting transcriptional enhancers are a major category of non-coding genome function and have been shown to control the fine-tuned spatiotemporal expression of genes during critical stages of craniofacial development1-3. However, a lack of accurate maps of the genomic location and cell type-specific in vivo activities of all craniofacial enhancers prevents their systematic exploration in human genetics studies. Here, we combined histone modification and chromatin accessibility profiling from different stages of human craniofacial development with single-cell analyses of the developing mouse face to create a comprehensive catalogue of the regulatory landscape of facial development at tissue- and single cell-resolution. In total, we identified approximately 14,000 enhancers across seven developmental stages from weeks 4 through 8 of human embryonic face development. We used transgenic mouse reporter assays to determine the in vivo activity patterns of human face enhancers predicted from these data. Across 16 in vivo validated human enhancers, we observed a rich diversity of craniofacial subregions in which these enhancers are active in vivo. To annotate the cell type specificities of human-mouse conserved enhancers, we performed single-cell RNA-seq and single-nucleus ATAC-seq of mouse craniofacial tissues from embryonic days e11.5 to e15.5. By integrating these data across species, we find that the majority (56%) of human craniofacial enhancers are functionally conserved in mice, providing cell type- and embryonic stage-resolved predictions of their in vivo activity profiles. Using retrospective analysis of known craniofacial enhancers in combination with single cell-resolved transgenic reporter assays, we demonstrate the utility of these data for predicting the in vivo cell type specificity of enhancers. Taken together, our data provide an expansive resource for genetic and developmental studies of human craniofacial development.

Topologically associating domain boundaries are required for normal genome function.

Topologically associating domain (TAD) boundaries partition the genome into distinct regulatory territories. Anecdotal evidence suggests that their disruption may interfere with normal gene expression and cause disease phenotypes1-3, but the overall extent to which this occurs remains unknown. Here we demonstrate that targeted deletions of TAD boundaries cause a range of disruptions to normal in vivo genome function and organismal development. We used CRISPR genome editing in mice to individually delete eight TAD boundaries (11-80 kb in size) from the genome. All deletions examined resulted in detectable molecular or organismal phenotypes, which included altered chromatin interactions or gene expression, reduced viability, and anatomical phenotypes. We observed changes in local 3D chromatin architecture in 7 of 8 (88%) cases, including the merging of TADs and altered contact frequencies within TADs adjacent to the deleted boundary. For 5 of 8 (63%) loci examined, boundary deletions were associated with increased embryonic lethality or other developmental phenotypes. For example, a TAD boundary deletion near Smad3/Smad6 caused complete embryonic lethality, while a deletion near Tbx5/Lhx5 resulted in a severe lung malformation. Our findings demonstrate the importance of TAD boundary sequences for in vivo genome function and reinforce the critical need to carefully consider the potential pathogenicity of noncoding deletions affecting TAD boundaries in clinical genetics screening.

Activation of AcvR1-Mediated Signaling Results in Semilunar Valve Defects.

Calcific aortic valve disease (CAVD) is a common cardiac defect, particularly in the aging population. While several risk factors, such as bi-leaflet valve structure and old age, have been identified in CAVD pathogenesis, molecular mechanisms resulting in this condition are still under active investigation. Bone morphogenetic protein signaling via the activin type I receptor (AcvRI) plays an important role during physiological and pathological processes involving calcification, e.g., bone formation and heterotopic ossification. In addition, AcvRI is required for normal cardiac valve development, yet its role in aortic valve disease, if any, is currently unknown. Here, we induced the expression of constitutively active AcvRI in developing mouse embryos in the endocardium and in cells at the valve leaflet-wall junction that are not of endocardium origin using the Nfac1Cre transgene. The mutant mice were born alive, but showed thickened aortic and pulmonary valve leaflets during the early postnatal period. Adult mutant mice developed aortic stenosis with high frequency, sclerotic aortic valves, and displayed Alcian Blue-positive hypertrophic chondrocyte-like cells at the leaflet-wall junction. Calcification was only seen with low penetrance. In addition, we observed that the expression levels of gene sets associated with inflammation-related cytokine signaling, smooth muscle cell contraction, and cGMP signaling were altered in the mutants when compared with those of the controls. This work shows that, in a mouse model, such continuous AcvRI activity in the Nfatc1Cre recombination domain results in pathological changes in the aortic valve structure and function.

Trim33 is required for appropriate development of pre-cardiogenic mesoderm.

Congenital cardiac malformations are among the most common birth defects in humans. Here we show that Trim33, a member of the Tif1 subfamily of tripartite domain containing transcriptional cofactors, is required for appropriate differentiation of the pre-cardiogenic mesoderm during a narrow time window in late gastrulation. While mesoderm-specific Trim33 mutants did not display noticeable phenotypes, epiblast-specific Trim33 mutant embryos developed ventricular septal defects, showed sparse trabeculation and abnormally thin compact myocardium, and died as a result of cardiac failure during late gestation. Differentiating embryoid bodies deficient in Trim33 showed an enrichment of gene sets associated with cardiac differentiation and contractility, while the total number of cardiac precursor cells was reduced. Concordantly, cardiac progenitor cell proliferation was reduced in Trim33-deficient embryos. ChIP-Seq performed using antibodies against Trim33 in differentiating embryoid bodies revealed more than 4000 peaks, which were significantly enriched close to genes implicated in stem cell maintenance and mesoderm development. Nearly half of the Trim33 peaks overlapped with binding sites of the Ctcf insulator protein. Our results suggest that Trim33 is required for appropriate differentiation of precardiogenic mesoderm during late gastrulation and that it will likely mediate some of its functions via multi-protein complexes, many of which include the chromatin architectural and insulator protein Ctcf.

Trim33 regulates early maturation of mouse embryoid bodies in vitro.

Embryonic stem cells (ESCs) are an established model for investigating developmental processes, disease conditions, tissue regeneration and therapeutic targets. Previous studies have shown that tripartite motif-containing 33 protein (Trim33) functions as a chromatin reader during Nodal-induced mesoderm induction. Here we report that despite reduced proliferation, mouse ESCs deficient in Trim33 remained pluripotent when cultured under non-differentiating conditions. However, when induced to differentiate to embryoid bodies (EBs), the mutant cultures showed increased cell shedding and apoptosis at day 3 of differentiation. Gene set enrichment analysis (GSEA) indicated that several molecular functions associated with cell survival, transcriptional/translational activity and growth factor signaling were affected already by the second day of differentiation in Trim33-deficient EBs. Consistent with increased apoptosis, expression of Rac1, a critical factor for EB cell survival, was reduced in Trim33 mutant EBs. In addition, a set of genes involved in regulation of pluripotency was upregulated in mutant EBs. Our results suggest that Trim33 regulates early maturation of mouse embryoid bodies in vitro.

Elevated Fibroblast Growth Factor Signaling Is Critical for the Pathogenesis of the Dwarfism in Evc2/Limbin Mutant Mice.

Ellis-van Creveld (EvC) syndrome is a skeletal dysplasia, characterized by short limbs, postaxial polydactyly, and dental abnormalities. EvC syndrome is also categorized as a ciliopathy because of ciliary localization of proteins encoded by the two causative genes, EVC and EVC2 (aka LIMBIN). While recent studies demonstrated important roles for EVC/EVC2 in Hedgehog signaling, there is still little known about the pathophysiological mechanisms underlying the skeletal dysplasia features of EvC patients, and in particular why limb development is affected, but not other aspects of organogenesis that also require Hedgehog signaling. In this report, we comprehensively analyze limb skeletogenesis in Evc2 mutant mice and in cell and tissue cultures derived from these mice. Both in vivo and in vitro data demonstrate elevated Fibroblast Growth Factor (FGF) signaling in Evc2 mutant growth plates, in addition to compromised but not abrogated Hedgehog-PTHrP feedback loop. Elevation of FGF signaling, mainly due to increased Fgf18 expression upon inactivation of Evc2 in the perichondrium, critically contributes to the pathogenesis of limb dwarfism. The limb dwarfism phenotype is partially rescued by inactivation of one allele of Fgf18 in the Evc2 mutant mice. Taken together, our data uncover a novel pathogenic mechanism to understand limb dwarfism in patients with Ellis-van Creveld syndrome.

Generation of Evc2/Limbin global and conditional KO mice and its roles during mineralized tissue formation.

Ellis-van Creveld (EvC) syndrome (OMIM 225500) is an autosomal recessive disease characterized with chondrodysplastic dwarfism in association with abnormalities in oral cavity. Ciliary proteins EVC and EVC2 have been identified as causative genes and they play an important role on Hedgehog signal transduction. We have also identified a causative gene LIMBIN for bovine chondrodysplastic dwarfism (bcd) that is later identified as the bovine ortholog of EVC2. Here, we report generation of conventional and conditional mutant Evc2/Limbin alleles that mimics mutations found in EvC patients and bcd cattle. Resulted homozygous mice showed no ciliary localization of EVC2 and EVC and displayed reduced Hedgehog signaling activity in association with skeletal and oral defects similar to the EvC patients. Cartilage-specific disruption of Evc2/Limbin resulted in similar but milder skeletal defects, whereas osteoblast-specific disruption did not cause overt changes in skeletal system. Neural crest-specific disruption of Evc2/Limbin resulted in defective incisor growth similar to that seen in conventional knockouts; however, differentiation of amelobolasts was relatively normal in the conditional knockouts. These results showcased functions of EVC2/LIMBIN during formation of mineralized tissues. Availability of the conditional allele for this gene should facilitate further detailed analyses of the role of EVC2/LIMBIN in pathogenesis of EvC syndrome. genesis 53:612-626, 2015. © 2015 Wiley Periodicals, Inc.

AcvR1-mediated BMP signaling in second heart field is required for arterial pole development: implications for myocardial differentiation and regional identity.

BMP signaling plays an essential role in second heart field-derived heart and arterial trunk development, including myocardial differentiation, right ventricular growth, and interventricular, outflow tract and aortico-pulmonary septation. It is mediated by a number of different BMP ligands, and receptors, many of which are present simultaneously. The mechanisms by which they regulate morphogenetic events and degree of redundancy amongst them have still to be elucidated. We therefore assessed the role of BMP Type I receptor AcvR1 in anterior second heart field-derived cell development, and compared it with that of BmpR1a. By removing Acvr1 using the driver Mef2c[AHF]-Cre, we show that AcvR1 plays an essential role in arterial pole morphogenesis, identifying defects in outflow tract wall and cushion morphology that preceded a spectrum of septation defects from double outlet right ventricle to common arterial trunk in mutants. Its absence caused dysregulation in gene expression important for myocardial differentiation (Isl1, Fgf8) and regional identity (Tbx2, Tbx3, Tbx20, Tgfb2). Although these defects resemble to some degree those in the equivalent Bmpr1a mutant, a novel gene knock-in model in which Bmpr1a was expressed in the Acvr1 locus only partially restored septation in Acvr1 mutants. These data show that both BmpR1a and AcvR1 are needed for normal heart development, in which they play some non-redundant roles, and refine our understanding of the genetic and morphogenetic processes underlying Bmp-mediated heart development important in human congenital heart disease.

Enamel malformations associated with a defined dentin sialophosphoprotein mutation in two families.

Dentin sialophosphoprotein (DSPP) mutations cause dentin dysplasia type II (DD-II) and dentinogenesis imperfecta types II and III (DGI-II and DGI-III, respectively). We identified two kindreds with DGI-II who exhibited vertical bands of hypoplastic enamel. Both families had a previously reported DSPP mutation that segregated with the disease phenotype. Oral photographs and dental radiographs of four affected and one unaffected participant in one family and of the proband in the second family were used to document the dental phenotypes. We aligned the 33 unique allelic DSPP sequences showing variable patterns of insertions and deletions (indels), generated a merged dentin phosphoprotein (DPP) sequence that includes sequences from all DSPP length haplotypes, and mapped the known DSPP mutations in this context. Analyses of the DSPP sequence changes and their probable effects on protein expression, as well as published findings of the dental phenotype in Dspp null mice, support the hypothesis that all DSPP mutations cause pathology through dominant-negative effects. Noting that Dspp is transiently expressed by mouse pre-ameloblasts during formation of the dentino-enamel junction, we hypothesize that DSPP dominant-negative effects potentially cause cellular pathology in pre-ameloblasts that, in turn, causes enamel defects. We conclude that enamel defects can be part of the dental phenotype caused by DSPP mutations, although DSPP is not critical for dental enamel formation.