High incidence and geographic distribution of cleft palate cases in Finland are associated with a regulatory variant in IRF6.

In Finland the frequency of isolated cleft palate (CP) is higher than that of isolated cleft lip with or without cleft palate (CL/P). This trend contrasts to that in other European countries but its genetic underpinnings are unknown. We performed a genome-wide association study for orofacial clefts, which include CL/P and CP, in the Finnish population. We identified rs570516915, a single nucleotide polymorphism that is highly enriched in Finns and Estonians, as being strongly associated with CP ( P = 5.25 × 10 -34 , OR = 8.65, 95% CI 6.11-12.25), but not with CL/P ( P = 7.2 × 10 -5 ), with genome-wide significance. The risk allele frequency of rs570516915 parallels the regional variation of CP prevalence in Finland, and the association was replicated in independent cohorts of CP cases from Finland ( P = 8.82 × 10 -28 ) and Estonia ( P = 1.25 × 10 -5 ). The risk allele of rs570516915 disrupts a conserved binding site for the transcription factor IRF6 within a previously characterized enhancer upstream of the IRF6 gene. Through reporter assay experiments we found that the risk allele of rs570516915 diminishes the enhancer activity. Oral epithelial cells derived from CRISPR-Cas9 edited induced pluripotent stem cells demonstrate that the CP-associated allele of rs570516915 concomitantly decreases the binding of IRF6 and the expression level of IRF6 , suggesting impaired IRF6 autoregulation as a molecular mechanism underlying the risk for CP.

A conserved molecular logic for neurogenesis to gliogenesis switch in the cerebral cortex.

During development, neural stem cells in the cerebral cortex, also known as radial glial cells (RGCs), generate excitatory neurons, followed by production of cortical macroglia and inhibitory neurons that migrate to the olfactory bulb (OB). Understanding the mechanisms for this lineage switch is fundamental for unraveling how proper numbers of diverse neuronal and glial cell types are controlled. We and others recently showed that Sonic Hedgehog (Shh) signaling promotes the cortical RGC lineage switch to generate cortical oligodendrocytes and OB interneurons. During this process, cortical RGCs generate intermediate progenitor cells that express critical gliogenesis genes Ascl1, Egfr, and Olig2. The increased Ascl1 expression and appearance of Egfr+ and Olig2+ cortical progenitors are concurrent with the switch from excitatory neurogenesis to gliogenesis and OB interneuron neurogenesis in the cortex. While Shh signaling promotes Olig2 expression in the developing spinal cord, the exact mechanism for this transcriptional regulation is not known. Furthermore, the transcriptional regulation of Olig2 and Egfr has not been explored. Here, we show that in cortical progenitor cells, multiple regulatory programs, including Pax6 and Gli3, prevent precocious expression of Olig2, a gene essential for production of cortical oligodendrocytes and astrocytes. We identify multiple enhancers that control Olig2 expression in cortical progenitors and show that the mechanisms for regulating Olig2 expression are conserved between the mouse and human. Our study reveals evolutionarily conserved regulatory logic controlling the lineage switch of cortical neural stem cells.

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.

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.

Uncovering the mesendoderm gene regulatory network through multi-omic data integration.

Mesendodermal specification is one of the earliest events in embryogenesis, where cells first acquire distinct identities. Cell differentiation is a highly regulated process that involves the function of numerous transcription factors (TFs) and signaling molecules, which can be described with gene regulatory networks (GRNs). Cell differentiation GRNs are difficult to build because existing mechanistic methods are low throughput, and high-throughput methods tend to be non-mechanistic. Additionally, integrating highly dimensional data composed of more than two data types is challenging. Here, we use linked self-organizing maps to combine chromatin immunoprecipitation sequencing (ChIP-seq)/ATAC-seq with temporal, spatial, and perturbation RNA sequencing (RNA-seq) data from Xenopus tropicalis mesendoderm development to build a high-resolution genome scale mechanistic GRN. We recover both known and previously unsuspected TF-DNA/TF-TF interactions validated through reporter assays. Our analysis provides insights into transcriptional regulation of early cell fate decisions and provides a general approach to building GRNs using highly dimensional multi-omic datasets.

Anterograde regulation of mitochondrial genes and FGF21 signaling by hepatic LSD1.

Mitochondrial biogenesis and function are controlled by anterograde regulatory pathways involving more than 1000 nuclear-encoded proteins. Transcriptional networks controlling the nuclear-encoded mitochondrial genes remain to be fully elucidated. Here, we show that histone demethylase LSD1 KO from adult mouse liver (LSD1-LKO) reduces the expression of one-third of all nuclear-encoded mitochondrial genes and decreases mitochondrial biogenesis and function. LSD1-modulated histone methylation epigenetically regulates nuclear-encoded mitochondrial genes. Furthermore, LSD1 regulates gene expression and protein methylation of nicotinamide mononucleotide adenylyltransferase 1 (NMNAT1), which controls the final step of NAD+ synthesis and limits NAD+ availability in the nucleus. Lsd1 KO reduces NAD+-dependent SIRT1 and SIRT7 deacetylase activity, leading to hyperacetylation and hypofunctioning of GABPß and PGC-1α, the major transcriptional factor/cofactor for nuclear-encoded mitochondrial genes. Despite the reduced mitochondrial function in the liver, LSD1-LKO mice are protected from diet-induced hepatic steatosis and glucose intolerance, partially due to induction of hepatokine FGF21. Thus, LSD1 orchestrates a core regulatory network involving epigenetic modifications and NAD+ synthesis to control mitochondrial function and hepatokine production.

Sox17 and ß-catenin co-occupy Wnt-responsive enhancers to govern the endoderm gene regulatory network.

Lineage specification is governed by gene regulatory networks (GRNs) that integrate the activity of signaling effectors and transcription factors (TFs) on enhancers. Sox17 is a key transcriptional regulator of definitive endoderm development, and yet, its genomic targets remain largely uncharacterized. Here, using genomic approaches and epistasis experiments, we define the Sox17-governed endoderm GRN in Xenopus gastrulae. We show that Sox17 functionally interacts with the canonical Wnt pathway to specify and pattern the endoderm while repressing alternative mesectoderm fates. Sox17 and ß-catenin co-occupy hundreds of key enhancers. In some cases, Sox17 and ß-catenin synergistically activate transcription apparently independent of Tcfs, whereas on other enhancers, Sox17 represses ß-catenin/Tcf-mediated transcription to spatially restrict gene expression domains. Our findings establish Sox17 as a tissue-specific modifier of Wnt responses and point to a novel paradigm where genomic specificity of Wnt/ß-catenin transcription is determined through functional interactions between lineage-specific Sox TFs and ß-catenin/Tcf transcriptional complexes. Given the ubiquitous nature of Sox TFs and Wnt signaling, this mechanism has important implications across a diverse range of developmental and disease contexts.

Foxh1/Nodal Defines Context-Specific Direct Maternal Wnt/ß-Catenin Target Gene Regulation in Early Development.

Although Wnt/ß-catenin signaling is generally conserved and well understood, the regulatory mechanisms controlling context-specific direct Wnt target gene expression in development and disease are still unclear. The onset of zygotic gene transcription in early embryogenesis represents an ideal, accessible experimental system to investigate context-specific direct Wnt target gene regulation. We combine transcriptomics using RNA-seq with genome-wide ß-catenin association using ChIP-seq to identify stage-specific direct Wnt target genes. We propose coherent feedforward regulation involving two distinct classes of direct maternal Wnt target genes, which differ both in expression and persistence of ß-catenin association. We discover that genomic ß-catenin association overlaps with Foxh1-associated regulatory sequences and demonstrate that direct maternal Wnt target gene expression requires Foxh1 function and Nodal/Tgfß signaling. Our results support a new paradigm for direct Wnt target gene co-regulation with context-specific mechanisms that will inform future studies of embryonic development and more widely stem cell-mediated homeostasis and human disease.

Early Xenopus gene regulatory programs, chromatin states, and the role of maternal transcription factors.

For decades, the early development of the Xenopus embryo has been an essential model system to study the gene regulatory mechanisms that govern cellular specification. At the top of the hierarchy of gene regulatory networks, maternally deposited transcription factors initiate this process and regulate the expression of zygotic genes that give rise to three distinctive germ layer cell types (ectoderm, mesoderm, and endoderm), and subsequent generation of organ precursors. The onset of germ layer specification is also closely coupled with changes associated with chromatin modifications. This review will examine the timing of maternal transcription factors initiating the zygotic genome activation, the epigenetic landscape of embryonic chromatin, and the network structure that governs the process.

Endodermal Maternal Transcription Factors Establish Super-Enhancers during Zygotic Genome Activation.

Elucidation of the sequence of events underlying the dynamic interaction between transcription factors and chromatin states is essential. Maternal transcription factors function at the top of the regulatory hierarchy to specify the primary germ layers at the onset of zygotic genome activation (ZGA). We focus on the formation of endoderm progenitor cells and examine the interactions between maternal transcription factors and chromatin state changes underlying the cell specification process. Endoderm-specific factors Otx1 and Vegt together with Foxh1 orchestrate endoderm formation by coordinated binding to select regulatory regions. These interactions occur before the deposition of enhancer histone marks around the regulatory regions, and these TFs recruit RNA polymerase II, regulate enhancer activity, and establish super-enhancers associated with important endodermal genes. Therefore, maternal transcription factors Otx1, Vegt, and Foxh1 combinatorially regulate the activity of super-enhancers, which in turn activate key lineage-specifying genes during ZGA.

Morpholinos Do Not Elicit an Innate Immune Response during Early Xenopus Embryogenesis.

It has recently been reported that a common side effect of translation-blocking morpholino antisense oligonucleotides is the induction of a set of innate immune response genes in Xenopus embryos and that splicing-blocking morpholinos lead to unexpected off-target mis-splicing events. Here, we present an analysis of all publicly available Xenopus RNA sequencing (RNA-seq) data in a reexamination of the effects of translation-blocking morpholinos on the innate immune response. Our analysis does not support the authors' general conclusion, which was based on a limited number of RNA-seq datasets. Moreover, the strong induction of an immune response appears to be specific to the tbxt/tbxt2 morpholinos. The more comprehensive study presented here indicates that using morpholinos for targeted gene knockdowns remains of considerable value for the rapid identification of gene function.

Foxh1 Occupies cis-Regulatory Modules Prior to Dynamic Transcription Factor Interactions Controlling the Mesendoderm Gene Program.

The interplay between transcription factors and chromatin dictates gene regulatory network activity. Germ layer specification is tightly coupled with zygotic gene activation and, in most metazoans, is dependent upon maternal factors. We explore the dynamic genome-wide interactions of Foxh1, a maternal transcription factor that mediates Nodal/TGF-ß signaling, with cis-regulatory modules (CRMs) during mesendodermal specification. Foxh1 marks CRMs during cleavage stages and recruits the co-repressor Tle/Groucho in the early blastula. We highlight a population of CRMs that are continuously occupied by Foxh1 and show that they are marked by H3K4me1, Ep300, and Fox/Sox/Smad motifs, suggesting interplay between these factors in gene regulation. We also propose a molecular "hand-off" between maternal Foxh1 and zygotic Foxa at these CRMs to maintain enhancer activation. Our findings suggest that Foxh1 functions at the top of a hierarchy of interactions by marking developmental genes for activation, beginning with the onset of zygotic gene expression.

A gene regulatory program controlling early Xenopus mesendoderm formation: Network conservation and motifs.

Germ layer formation is among the earliest differentiation events in metazoan embryos. In triploblasts, three germ layers are formed, among which the endoderm gives rise to the epithelial lining of the gut tube and associated organs including the liver, pancreas and lungs. In frogs (Xenopus), where early germ layer formation has been studied extensively, the process of endoderm specification involves the interplay of dozens of transcription factors. Here, we review the interactions between these factors, summarized in a transcriptional gene regulatory network (GRN). We highlight regulatory connections conserved between frog, fish, mouse, and human endodermal lineages. Especially prominent is the conserved role and regulatory targets of the Nodal signaling pathway and the T-box transcription factors, Vegt and Eomes. Additionally, we highlight network topologies and motifs, and speculate on their possible roles in development.

Developmentally regulated long non-coding RNAs in Xenopus tropicalis.

Advances in RNA sequencing technologies have led to the surprising discovery that a vast number of transcripts emanate from regions of the genome that are not part of coding genes. Although some of the smaller ncRNAs such as microRNAs have well-characterized functions, the majority of long ncRNA (lncRNA) functions remain poorly understood. Understanding the significance of lncRNAs is an important challenge facing biology today. A powerful approach to uncovering the function of lncRNAs is to explore temporal and spatial expression profiling. This may be particularly useful for classes of lncRNAs that have developmentally important roles as the expression of such lncRNAs will be expected to be both spatially and temporally regulated during development. Here, we take advantage of our ultra-high frequency (temporal) sampling of Xenopus embryos to analyze gene expression trajectories of lncRNA transcripts over the first 3 days of development. We computationally identify 5689 potential single- and multi-exon lncRNAs. These lncRNAs demonstrate clear dynamic expression patterns. A subset of them displays highly correlative temporal expression profiles with respect to those of the neighboring genes. We also identified spatially localized lncRNAs in the gastrula stage embryo. These results suggest that lncRNAs have regulatory roles during early embryonic development.

A catalog of Xenopus tropicalis transcription factors and their regional expression in the early gastrula stage embryo.

Gene regulatory networks (GRNs) involve highly combinatorial interactions between transcription factors and short sequence motifs in cis-regulatory modules of target genes to control cellular phenotypes. The GRNs specifying most cell types are largely unknown and are the subject of wide interest. A catalog of transcription factors is a valuable tool toward obtaining a deeper understanding of the role of these critical effectors in any biological setting. Here we present a comprehensive catalog of the transcription factors for the diploid frog Xenopus tropicalis. We identify 1235 genes encoding DNA-binding transcription factors, comparable to the numbers found in typical mammalian species. In detail, the repertoire of X. tropicalis transcription factor genes is nearly identical to human and mouse, with the exception of zinc finger family members, and a small number of species/lineage-specific gene duplications and losses relative to the mammalian repertoires. We applied this resource to the identification of transcription factors differentially expressed in the early gastrula stage embryo. We find transcription factor enrichment in Spemann's organizer, the ventral mesoderm, ectoderm and endoderm, and report 218 TFs that show regionalized expression patterns at this stage. Many of these have not been previously reported as expressed in the early embryo, suggesting thus far unappreciated roles for many transcription factors in the GRNs regulating early development. We expect our transcription factor catalog will facilitate myriad studies using Xenopus as a model system to understand basic biology and human disease.