Multimodal spatiotemporal transcriptomic resolution of embryonic palate osteogenesis.

The terminal differentiation of osteoblasts and subsequent formation of bone marks an important phase in palate development that leads to the separation of the oral and nasal cavities. While the morphogenetic events preceding palatal osteogenesis are well explored, major gaps remain in our understanding of the molecular mechanisms driving the formation of this bony union of the fusing palate. Through bulk, single-nucleus, and spatially resolved RNA-sequencing analyses of the developing secondary palate, we identify a shift in transcriptional programming between embryonic days 14.5 and 15.5 pinpointing the onset of osteogenesis. We define spatially restricted expression patterns of key osteogenic marker genes that are differentially expressed between these developmental timepoints. Finally, we identify genes in the palate highly expressed by palate nasal epithelial cells, also enriched within palatal osteogenic mesenchymal cells. This investigation provides a relevant framework to advance palate-specific diagnostic and therapeutic biomarker discovery.

Integrated spatiotemporal transcriptomic resolution of embryonic palate osteogenesis.

The differentiation of osteoblasts and the subsequent formation of bone marks an important terminal phase in palate formation that leads to the separation of the oral and nasal cavities. While the developmental events that precede palatal osteogenesis are well explored, major gaps remain in our understanding of the molecular mechanisms that lead to the bony union of fusing palatal shelves. Herein, the timeline of osteogenic transcriptional programming is unveiled in the embryonic palate by way of integrated bulk, single-cell, and spatially resolved RNA-seq analyses. We define spatially restricted expression patterns of key marker genes, both regulatory and structural, that are differentially expressed during palatal fusion, including the identification of several novel genes ( Deup1, Dynlrb2, Lrrc23 ) spatially restricted in expression to the palate, providing a relevant framework for future studies that identify new candidate genes for cleft palate anomalies in humans as well as the timing of mammalian embryonic palatal osteogenesis.

Absence of TRIC-B from type XIV Osteogenesis Imperfecta osteoblasts alters cell adhesion and mitochondrial function - A multi-omics study.

Osteogenesis Imperfecta (OI) is a heritable collagen-related bone dysplasia characterized by bone fractures, growth deficiency and skeletal deformity. Type XIV OI is a recessive OI form caused by null mutations in TMEM38B, which encodes the ER membrane intracellular cation channel TRIC-B. Previously, we showed that absence of TMEM38B alters calcium flux in the ER of OI patient osteoblasts and fibroblasts, which further disrupts collagen synthesis and secretion. How the absence of TMEM38B affects osteoblast function is still poorly understood. Here we further investigated the role of TMEM38B in human osteoblast differentiation and mineralization. TMEM38B-null osteoblasts showed altered expression of osteoblast marker genes and decreased mineralization. RNA-Seq analysis revealed that cell-cell adhesion was one of the most downregulated pathways in TMEM38B-null osteoblasts, with further validation by real-time PCR and Western blot. Gap and tight junction proteins were also decreased by TRIC-B absence, both in patient osteoblasts and in calvarial osteoblasts of Tmem38b-null mice. Disrupted cell adhesion decreased mutant cell proliferation and cell cycle progression. An important novel finding was that TMEM38B-null osteoblasts had elongated mitochondria with altered fusion and fission markers, MFN2 and DRP1. In addition, TMEM38B-null osteoblasts exhibited a significant increase in superoxide production in mitochondria, further supporting mitochondrial dysfunction. Together these results emphasize the novel role of TMEM38B/TRIC-B in osteoblast differentiation, affecting cell-cell adhesion processes, gap and tight junction, proliferation, cell cycle, and mitochondrial function.

Enhancer-promoter interactions can bypass CTCF-mediated boundaries and contribute to phenotypic robustness.

How enhancers activate their distal target promoters remains incompletely understood. Here we dissect how CTCF-mediated loops facilitate and restrict such regulatory interactions. Using an allelic series of mouse mutants, we show that CTCF is neither required for the interaction of the Sox2 gene with distal enhancers, nor for its expression. Insertion of various combinations of CTCF motifs, between Sox2 and its distal enhancers, generated boundaries with varying degrees of insulation that directly correlated with reduced transcriptional output. However, in both epiblast and neural tissues, enhancer contacts and transcriptional induction could not be fully abolished, and insertions failed to disrupt implantation and neurogenesis. In contrast, Sox2 expression was undetectable in the anterior foregut of mutants carrying the strongest boundaries, and these animals fully phenocopied loss of SOX2 in this tissue. We propose that enhancer clusters with a high density of regulatory activity can better overcome physical barriers to maintain faithful gene expression and phenotypic robustness.

Decreasing Wapl dosage partially corrects embryonic growth and brain transcriptome phenotypes in Nipbl+/- embryos.

Cohesin rings interact with DNA and modulate the expression of thousands of genes. NIPBL loads cohesin onto chromosomes, and WAPL takes it off. Haploinsufficiency for NIPBL causes a developmental disorder, Cornelia de Lange syndrome (CdLS), that is modeled by Nipbl+/- mice. Mutations in WAPL have not been shown to cause disease or gene expression changes in mammals. Here, we show dysregulation of >1000 genes in WaplΔ/+ embryonic mouse brain. The patterns of dysregulation are highly similar in Wapl and Nipbl heterozygotes, suggesting that Wapl mutations may also cause human disease. Since WAPL and NIPBL have opposite effects on cohesin's association with DNA, we asked whether decreasing Wapl dosage could correct phenotypes seen in Nipbl+/- mice. Gene expression and embryonic growth are partially corrected, but perinatal lethality is not. Our data are consistent with the view that cohesin dynamics play a key role in regulating gene expression.

An epigenome atlas of neural progenitors within the embryonic mouse forebrain.

A comprehensive characterization of epigenomic organization in the embryonic mouse forebrain will enhance our understanding of neurodevelopment and provide insight into mechanisms of neurological disease. Here we collected single-cell chromatin accessibility profiles from four distinct neurogenic regions of the embryonic mouse forebrain using single nuclei ATAC-Seq (snATAC-Seq). We identified thousands of differentially accessible peaks, many restricted to distinct progenitor cell types or brain regions. We integrated snATAC-Seq and single cell transcriptome data to characterize changes of chromatin accessibility at enhancers and promoters with associated transcript abundance. Multi-modal integration of histone modifications (CUT&Tag and CUT&RUN), promoter-enhancer interactions (Capture-C) and high-order chromatin structure (Hi-C) extended these initial observations. This dataset reveals a diverse chromatin landscape with region-specific regulatory mechanisms and genomic interactions in distinct neurogenic regions of the embryonic mouse brain and represents an extensive public resource of a 'ground truth' epigenomic landscape at this critical stage of neurogenesis.

Extensive co-binding and rapid redistribution of NANOG and GATA6 during emergence of divergent lineages.

Fate-determining transcription factors (TFs) can promote lineage-restricted transcriptional programs from common progenitor states. The inner cell mass (ICM) of mouse blastocysts co-expresses the TFs NANOG and GATA6, which drive the bifurcation of the ICM into either the epiblast (Epi) or the primitive endoderm (PrE), respectively. Here, we induce GATA6 in embryonic stem cells-that also express NANOG-to characterize how a state of co-expression of opposing TFs resolves into divergent lineages. Surprisingly, we find that GATA6 and NANOG co-bind at the vast majority of Epi and PrE enhancers, a phenomenon we also observe in blastocysts. The co-bound state is followed by eviction and repression of Epi TFs, and quick remodeling of chromatin and enhancer-promoter contacts thus establishing the PrE lineage while repressing the Epi fate. We propose that co-binding of GATA6 and NANOG at shared enhancers maintains ICM plasticity and promotes the rapid establishment of Epi- and PrE-specific transcriptional programs.

Transcriptional heterogeneity of ventricular zone cells in the ganglionic eminences of the mouse forebrain.

The ventricular zone (VZ) of the nervous system contains radial glia cells that were originally considered relatively homogenous in their gene expression, but a detailed characterization of transcriptional diversity in these VZ cells has not been reported. Here, we performed single-cell RNA sequencing to characterize transcriptional heterogeneity of neural progenitors within the VZ and subventricular zone (SVZ) of the ganglionic eminences (GEs), the source of all forebrain GABAergic neurons. By using a transgenic mouse line to enrich for VZ cells, we characterize significant transcriptional heterogeneity, both between GEs and within spatial subdomains of specific GEs. Additionally, we observe differential gene expression between E12.5 and E14.5 VZ cells, which could provide insights into temporal changes in cell fate. Together, our results reveal a previously unknown spatial and temporal genetic diversity of VZ cells in the ventral forebrain that will aid our understanding of initial fate decisions in the forebrain.

The histone demethylase Lsd1 regulates multiple repressive gene programs during T cell development.

Analysis of the transcriptional profiles of developing thymocytes has shown that T lineage commitment is associated with loss of stem cell and early progenitor gene signatures and the acquisition of T cell gene signatures. Less well understood are the epigenetic alterations that accompany or enable these transcriptional changes. Here, we show that the histone demethylase Lsd1 (Kdm1a) performs a key role in extinguishing stem/progenitor transcriptional programs in addition to key repressive gene programs during thymocyte maturation. Deletion of Lsd1 caused a block in late T cell development and resulted in overexpression of interferon response genes as well as genes regulated by the Gfi1, Bcl6, and, most prominently, Bcl11b transcriptional repressors in CD4+CD8+ thymocytes. Transcriptional overexpression in Lsd1-deficient thymocytes was not always associated with increased H3K4 trimethylation at gene promoters, indicating that Lsd1 indirectly affects the expression of many genes. Together, these results identify a critical function for Lsd1 in the epigenetic regulation of multiple repressive gene signatures during T cell development.

NMDARs Drive the Expression of Neuropsychiatric Disorder Risk Genes Within GABAergic Interneuron Subtypes in the Juvenile Brain.

Medial ganglionic eminence (MGE)-derived parvalbumin (PV)+, somatostatin (SST)+and Neurogliaform (NGFC)-type cortical and hippocampal interneurons, have distinct molecular, anatomical, and physiological properties. However, the molecular mechanisms regulating their maturation remain poorly understood. Here, via single-cell transcriptomics, we show that the obligate NMDA-type glutamate receptor (NMDAR) subunit gene Grin1 mediates transcriptional regulation of gene expression in specific subtypes of MGE-derived interneurons, leading to altered subtype abundances. Notably, MGE-specific early developmental Grin1 loss results in a broad downregulation of diverse transcriptional, synaptogenic and membrane excitability regulatory programs in the juvenile brain. These widespread gene expression abnormalities mirror aberrations that are typically associated with neurodevelopmental disorders. Our study hence provides a road map for the systematic examination of NMDAR signaling in interneuron subtypes, revealing potential MGE-specific genetic targets that could instruct future therapies of psychiatric disorders.

Cardiac pathologies in mouse loss of imprinting models are due to misexpression of H19 long noncoding RNA.

Maternal loss of imprinting (LOI) at the H19/IGF2 locus results in biallelic IGF2 and reduced H19 expression and is associated with Beckwith--Wiedemann syndrome (BWS). We use mouse models for LOI to understand the relative importance of Igf2 and H19 mis-expression in BWS phenotypes. Here we focus on cardiovascular phenotypes and show that neonatal cardiomegaly is exclusively dependent on increased Igf2. Circulating IGF2 binds cardiomyocyte receptors to hyperactivate mTOR signaling, resulting in cellular hyperplasia and hypertrophy. These Igf2-dependent phenotypes are transient: cardiac size returns to normal once Igf2 expression is suppressed postnatally. However, reduced H19 expression is sufficient to cause progressive heart pathologies including fibrosis and reduced ventricular function. In the heart, H19 expression is primarily in endothelial cells (ECs) and regulates EC differentiation both in vivo and in vitro. Finally, we establish novel mouse models to show that cardiac phenotypes depend on H19 lncRNA interactions with Mirlet7 microRNAs.

Linkage-specific deubiquitylation by OTUD5 defines an embryonic pathway intolerant to genomic variation.

Reversible modification of proteins with linkage-specific ubiquitin chains is critical for intracellular signaling. Information on physiological roles and underlying mechanisms of particular ubiquitin linkages during human development are limited. Here, relying on genomic constraint scores, we identify 10 patients with multiple congenital anomalies caused by hemizygous variants in OTUD5, encoding a K48/K63 linkage-specific deubiquitylase. By studying these mutations, we find that OTUD5 controls neuroectodermal differentiation through cleaving K48-linked ubiquitin chains to counteract degradation of select chromatin regulators (e.g., ARID1A/B, histone deacetylase 2, and HCF1), mutations of which underlie diseases that exhibit phenotypic overlap with OTUD5 patients. Loss of OTUD5 during differentiation leads to less accessible chromatin at neuroectodermal enhancers and aberrant gene expression. Our study describes a previously unidentified disorder we name LINKED (LINKage-specific deubiquitylation deficiency-induced Embryonic Defects) syndrome and reveals linkage-specific ubiquitin cleavage from chromatin remodelers as an essential signaling mode that coordinates chromatin remodeling during embryogenesis.

KRAB-zinc finger protein gene expansion in response to active retrotransposons in the murine lineage.

The Krüppel-associated box zinc finger protein (KRAB-ZFP) family diversified in mammals. The majority of human KRAB-ZFPs bind transposable elements (TEs), however, since most TEs are inactive in humans it is unclear whether KRAB-ZFPs emerged to suppress TEs. We demonstrate that many recently emerged murine KRAB-ZFPs also bind to TEs, including the active ETn, IAP, and L1 families. Using a CRISPR/Cas9-based engineering approach, we genetically deleted five large clusters of KRAB-ZFPs and demonstrate that target TEs are de-repressed, unleashing TE-encoded enhancers. Homozygous knockout mice lacking one of two KRAB-ZFP gene clusters on chromosome 2 and chromosome 4 were nonetheless viable. In pedigrees of chromosome 4 cluster KRAB-ZFP mutants, we identified numerous novel ETn insertions with a modest increase in mutants. Our data strongly support the current model that recent waves of retrotransposon activity drove the expansion of KRAB-ZFP genes in mice and that many KRAB-ZFPs play a redundant role restricting TE activity.

CD5 dynamically calibrates basal NF-κB signaling in T cells during thymic development and peripheral activation.

Immature T cells undergo a process of positive selection in the thymus when their new T cell receptor (TCR) engages and signals in response to self-peptides. As the T cell matures, a slew of negative regulatory molecules, including the inhibitory surface glycoprotein CD5, are up-regulated in proportion to the strength of the self-peptide signal. Together these regulators dampen TCR-proximal signaling and help avoid any subsequent peripheral activation of T cells by self-peptides. Paradoxically, antigen-specific T cells initially expressing more CD5 (CD5hi) have been found to better persist as effector/memory cells after a peripheral challenge. The molecular mechanisms underlying such a duality in CD5 function is not clear. We found that CD5 alters the basal activity of the NF-κB signaling in resting peripheral T cells. When CD5 was conditionally ablated, T cells were unable to maintain higher expression of the cytoplasmic NF-κB inhibitor IκBα. Consistent with this, resting CD5hi T cells expressed more of the NF-κB p65 protein than CD5lo cells, without significant increases in transcript levels, in the absence of TCR signals. This posttranslationally stabilized cellular NF-κB depot potentially confers a survival advantage to CD5hi T cells over CD5lo ones. Taken together, these data suggest a two-step model whereby the strength of self-peptide-induced TCR signal lead to the up-regulation of CD5, which subsequently maintains a proportional reserve of NF-κB in peripheral T cells poised for responding to agonistic antigen-driven T cell activation.

Somatic SMAD3-activating mutations cause melorheostosis by up-regulating the TGF-ß/SMAD pathway.

Melorheostosis is a rare sclerosing dysostosis characterized by asymmetric exuberant bone formation. Recently, we reported that somatic mosaicism for MAP2K1-activating mutations causes radiographical "dripping candle wax" melorheostosis. We now report somatic SMAD3 mutations in bone lesions of four unrelated patients with endosteal pattern melorheostosis. In vitro, the SMAD3 mutations stimulated the TGF-ß pathway in osteoblasts, enhanced nuclear translocation and target gene expression, and inhibited proliferation. Osteoblast differentiation and mineralization were stimulated by the SMAD3 mutation, consistent with higher mineralization in affected than in unaffected bone, but differing from MAP2K1 mutation-positive melorheostosis. Conversely, osteoblast differentiation and mineralization were inhibited when osteogenesis of affected osteoblasts was driven in the presence of BMP2. Transcriptome profiling displayed that TGF-ß pathway activation and ossification-related processes were significantly influenced by the SMAD3 mutation. Co-expression clustering illuminated melorheostosis pathophysiology, including alterations in ECM organization, cell growth, and interferon signaling. These data reveal antagonism of TGF-ß/SMAD3 activation by BMP signaling in SMAD3 mutation-positive endosteal melorheostosis, which may guide future therapies.

Ldb1 is required for Lmo2 oncogene-induced thymocyte self-renewal and T-cell acute lymphoblastic leukemia.

Prolonged or enhanced expression of the proto-oncogene Lmo2 is associated with a severe form of T-cell acute lymphoblastic leukemia (T-ALL), designated early T-cell precursor ALL, which is characterized by the aberrant self-renewal and subsequent oncogenic transformation of immature thymocytes. It has been suggested that Lmo2 exerts these effects by functioning as component of a multi-subunit transcription complex that includes the ubiquitous adapter Ldb1 along with b-HLH and/or GATA family transcription factors; however, direct experimental evidence for this mechanism is lacking. In this study, we investigated the importance of Ldb1 for Lmo2-induced T-ALL by conditional deletion of Ldb1 in thymocytes in an Lmo2 transgenic mouse model of T-ALL. Our results identify a critical requirement for Ldb1 in Lmo2-induced thymocyte self-renewal and thymocyte radiation resistance and for the transition of preleukemic thymocytes to overt T-ALL. Moreover, Ldb1 was also required for acquisition of the aberrant preleukemic ETP gene expression signature in immature Lmo2 transgenic thymocytes. Co-binding of Ldb1 and Lmo2 was detected at the promoters of key upregulated T-ALL driver genes (Hhex, Lyl1, and Nfe2) in preleukemic Lmo2 transgenic thymocytes, and binding of both Ldb1 and Lmo2 at these sites was reduced following Cre-mediated deletion of Ldb1. Together, these results identify a key role for Ldb1, a nonproto-oncogene, in T-ALL and support a model in which Lmo2-induced T-ALL results from failure to downregulate Ldb1/Lmo2-nucleated transcription complexes which normally function to enforce self-renewal in bone marrow hematopoietic progenitors.

Body Condition Scoring for Adult Zebrafish (Danio rerio).

Body condition scoring (BCS) is a simple, rapid, noninvasive tool used to assess body condition in animals. In this study, wedeveloped and validated a diagram-based BCS for adult zebrafish (Danio rerio), a popular research model. After receiving 20min of hands-on training regarding the scoring system, 5 people each rated 95 adult zebrafish. The fish then were euthanizedand measured to establish body condition indices (BMI and the Fulton K factor). Both condition indices were highly correlatedwith fish width. Using correlation data and observed trends in fish width, we established expected BCS definitions. Wevalidated the BCS definitions in 2 ways. First, we calculated the Pearson correlation coefficient between the average observedBCS and expected BCS; this statistic revealed very strong correlation between observed and expected BCS. In addition, weassessed the predictive power of BCS by using multinomial logistic regression and then applied the fitted model to evaluatethe accuracy of the predictions (BCS compared with BMI, 85%; BCS compared with K factor, 61%). Finally, to determine therobustness of BCS to variation among raters, we calculated the intraclass correlation coefficient and demonstrated high interrater reliability. In conclusion, adult zebrafish BCS can be used to quickly identify animals with different body conditionindices (thin to obese). In addition, the diagram-based chart is easy to use and implement accurately, with minimal training.

Circular RNAs and competing endogenous RNA (ceRNA) networks.

In recent years, advances in bioinformatics approaches have allowed a systematic characterization of circular RNAs (circRNAs) across a variety of cell types. Demonstration of cell type specificity, regulated expression, and conservation between species all suggest that circRNAs have functional importance. Especially, investigators have begun focusing on the possibility that circRNAs operate as part of competing endogenous RNA (ceRNA) regulatory networks that are proposed to play critical roles in normal development and in pathologic conditions like cancer.

Conditional ablation and conditional rescue models for Casq2 elucidate the role of development and of cell-type specific expression of Casq2 in the CPVT2 phenotype.

Cardiac calsequestrin (Casq2) associates with the ryanodine receptor 2 channel in the junctional sarcoplasmic reticulum to regulate Ca2+ release into the cytoplasm. Patients carrying mutations in CASQ2 display low resting heart rates under basal conditions and stress-induced polymorphic ventricular tachycardia (CPVT). In this study, we generate and characterize novel conditional deletion and conditional rescue mouse models to test the influence of developmental programs on the heart rate and CPVT phenotypes. We also compare the requirements for Casq2 function in the cardiac conduction system (CCS) and in working cardiomyocytes. Our study shows that the CPVT phenotype is dependent upon concurrent loss of Casq2 function in both the CCS and in working cardiomyocytes. Accordingly, restoration of Casq2 in only the CCS prevents CPVT. In addition, occurrence of CPVT is independent of the developmental history of Casq2-deficiency. In contrast, resting heart rate depends upon Casq2 gene activity only in the CCS and upon developmental history. Finally, our data support a model where low basal heart rate is a significant risk factor for CPVT.

Loss of imprinting mutations define both distinct and overlapping roles for misexpression of IGF2 and of H19 lncRNA.

Imprinted genes occur in discrete clusters that are coordinately regulated by shared DNA elements called Imprinting Control Regions. H19 and Igf2 are linked imprinted genes that play critical roles in development. Loss of imprinting (LOI) at the IGF2/H19 locus on the maternal chromosome is associated with the developmental disorder Beckwith Wiedemann Syndrome (BWS) and with several cancers. Here we use comprehensive genetic and genomic analyses to follow muscle development in a mouse model of BWS to dissect the separate and shared roles for misexpression of Igf2 and H19 in the disease phenotype. We show that LOI results in defects in muscle differentiation and hypertrophy and identify primary downstream targets: Igf2 overexpression results in over-activation of MAPK signaling while loss of H19 lncRNA prevents normal down regulation of p53 activity and therefore results in reduced AKT/mTOR signaling. Moreover, we demonstrate instances where H19 and Igf2 misexpression work separately, cooperatively, and antagonistically to establish the developmental phenotype. This study thus identifies new biochemical roles for the H19 lncRNA and underscores that LOI phenotypes are multigenic so that complex interactions will contribute to disease outcomes.