ABSTRACT
The gastrointestinal (GI) tract is complex and consists of multiple organs with unique functions. Rare gene variants can cause congenital malformations of the human GI tract, although the molecular basis of these has been poorly studied. We identified a patient with compound-heterozygous variants in RFX6 presenting with duodenal malrotation and atresia, implicating RFX6 in development of the proximal intestine. To identify how mutations in RFX6 impact intestinal patterning and function, we derived induced pluripotent stem cells from this patient to generate human intestinal organoids (HIOs). We identified that the duodenal HIOs and human tissues had mixed regional identity, with gastric and ileal features. CRISPR-mediated correction of RFX6 restored duodenal identity. We then used gain- and loss-of-function and transcriptomic approaches in HIOs and Xenopus embryos to identify that PDX1 is a downstream transcriptional target of RFX6 required for duodenal development. However, RFX6 had additional PDX1-independent transcriptional targets involving multiple components of signaling pathways that are required for establishing early regional identity in the GI tract. In summary, we have identified RFX6 as a key regulator in intestinal patterning that acts by regulating transcriptional and signaling pathways.
Subject(s)
Gene Expression Regulation, Developmental , Homeodomain Proteins , Organoids , Regulatory Factor X Transcription Factors , Trans-Activators , Humans , Regulatory Factor X Transcription Factors/genetics , Regulatory Factor X Transcription Factors/metabolism , Animals , Homeodomain Proteins/metabolism , Homeodomain Proteins/genetics , Trans-Activators/metabolism , Trans-Activators/genetics , Organoids/metabolism , Organoids/embryology , Duodenum/metabolism , Duodenum/embryology , Intestines/embryology , Intestinal Atresia/genetics , Induced Pluripotent Stem Cells/metabolism , Body Patterning/genetics , Signal Transduction/genetics , Mutation/geneticsABSTRACT
To investigate the co-development of vasculature, mesenchyme, and epithelium crucial for organogenesis and the acquisition of organ-specific characteristics, we constructed a human pluripotent stem cell-derived organoid system comprising lung or intestinal epithelium surrounded by organotypic mesenchyme and vasculature. We demonstrated the pivotal role of co-differentiating mesoderm and endoderm via precise BMP regulation in generating multilineage organoids and gut tube patterning. Single-cell RNA-seq analysis revealed organ specificity in endothelium and mesenchyme, and uncovered key ligands driving endothelial specification in the lung (e.g., WNT2B and Semaphorins) or intestine (e.g., GDF15). Upon transplantation under the kidney capsule in mice, these organoids further matured and developed perfusable human-specific sub-epithelial capillaries. Additionally, our model recapitulated the abnormal endothelial-epithelial crosstalk in patients with FOXF1 deletion or mutations. Multilineage organoids provide a unique platform to study developmental cues guiding endothelial and mesenchymal cell fate determination, and investigate intricate cell-cell communications in human organogenesis and disease. Highlights: BMP signaling fine-tunes the co-differentiation of mesoderm and endoderm.The cellular composition in multilineage organoids resembles that of human fetal organs.Mesenchyme and endothelium co-developed within the organoids adopt organ-specific characteristics.Multilineage organoids recapitulate abnormal endothelial-epithelial crosstalk in FOXF1-associated disorders.
ABSTRACT
INTRODUCTION: Heterogeneity in reported outcomes of infants with oesophageal atresia (OA) with or without tracheo-oesophageal fistula (TOF) prevents effective data pooling. Core outcome sets (COS) have been developed for many conditions to standardise outcome reporting, facilitate meta-analysis and improve the relevance of research for patients and families. Our aim is to develop an internationally-agreed, comprehensive COS for OA-TOF, relevant from birth through to transition and adulthood. METHODS AND ANALYSIS: A long list of outcomes will be generated using (1) a systematic review of existing studies on OA-TOF and (2) qualitative research with children (patients), adults (patients) and families involving focus groups, semistructured interviews and self-reported outcome activity packs. A two-phase Delphi survey will then be completed by four key stakeholder groups: (1) patients (paediatric and adult); (2) families; (3) healthcare professionals; and (4) researchers. Phase I will include stakeholders individually rating the importance and relevance of each long-listed outcome using a 9-point Likert scale, with the option to suggest additional outcomes not already included. During phase II, stakeholders will review summarised results from phase I relative to their own initial score and then will be asked to rescore the outcome based on this information. Responses from phase II will be summarised using descriptive statistics and a predefined definition of consensus for inclusion or exclusion of outcomes. Following the Delphi process, stakeholder experts will be invited to review data at a consensus meeting and agree on a COS for OA-TOF. ETHICS AND DISSEMINATION: Ethical approval was sought through the Health Research Authority via the Integrated Research Application System, registration no. 297026. However, approval was deemed not to be required, so study sponsorship and oversight were provided by Alder Hey Children's NHS Foundation Trust. The study has been prospectively registered with the COMET Initiative. The study will be published in an open access forum.
Subject(s)
Delphi Technique , Esophageal Atresia , Tracheoesophageal Fistula , Humans , Child , Adult , Research Design , Outcome Assessment, Health Care , Qualitative Research , Systematic Reviews as Topic , InfantABSTRACT
Disruptions in foregut morphogenesis can result in life-threatening conditions where the trachea and esophagus fail to separate properly, such as esophageal atresia (EA) and tracheoesophageal fistulas (TEF). The developmental basis of these congenital anomalies is poorly understood, but recent genome sequencing reveals that de novo variants in intracellular trafficking genes are enriched in EA/TEF patients. Here we show that mutation of orthologous genes in Xenopus disrupts trachea-esophageal separation similar to EA/TEF patients. We show that the Rab11a recycling endosome pathway is required to localize Vangl-Celsr polarity complexes at the cell surface where opposite sides of the common foregut tube fuse. Partial loss of endosome trafficking or the Vangl/Celsr complex disrupts epithelial polarity and cell division orientation. Mutant cells accumulate at the fusion point, fail to downregulate cadherin, and do not separate into distinct trachea and esophagus. These data provide new insights into the mechanisms of congenital anomalies and general paradigms of tissue fusion during organogenesis.
ABSTRACT
BACKGROUND: Specialized valve endothelial cell (VEC) populations are localized oriented to blood flow in developing aortic and mitral valves, but their roles in valve development and disease are unknown. In the aortic valve (AoV), a population of VECs on the fibrosa side expresses the transcription factor Prox1 together with genes found in lymphatic ECs. In this study, we examine Prox1's role in regulating a lymphatic-like gene network and promoting VEC diversity required for the development of the stratified trilaminar extracellular matrix (ECM) of murine AoV leaflets. METHODS: To determine whether disruption of Prox1 localization affects heart valve development, we generated mice (NFATc1enCre Prox1 gain-of-function) in which Prox1 is overexpressed on the ventricularis side of the AoV beginning in embryonic development. To identify potential targets of Prox1, we performed cleavage under targets and release using nuclease on wild-type and NFATc1enCre Prox1 gain-of-function AoVs with validation by colocalization in vivo using RNA in situ hybridization in NFATc1enCre Prox1 gain-of-function AoVs. Natural induction of Prox1 and target gene expression was evaluated in myxomatous AoVs in a mouse model of Marfan syndrome (Fbn1C1039G/+). RESULTS: The overexpression of Prox1 is sufficient to cause enlargement of AoVs by postnatal day (P)0, as well as a decrease in ventricularis-specific gene expression and disorganized interstitial ECM layers at P7. We identified potential targets of Prox1 known to play roles in lymphatic ECs including Flt1, Efnb2, Egfl7, and Cx37. Ectopic Prox1 colocalized with induced Flt1, Efnb2, and Cx37 expression in NFATc1enCre Prox1 gain-of-function AoVs. Moreover, in Marfan syndrome myxomatous AoVs, endogenous Prox1, and its identified targets, were ectopically induced in ventricularis side VECs. CONCLUSIONS: Our results support a role for Prox1 in localized lymphatic-like gene expression on the fibrosa side of the AoV. Furthermore, localized VEC specialization is required for development of the stratified trilaminar ECM critical for AoV function and is dysregulated in congenitally malformed valves.
Subject(s)
Aortic Valve , Marfan Syndrome , Mice , Animals , Aortic Valve/metabolism , Marfan Syndrome/metabolism , Extracellular Matrix/metabolism , Transcription Factors/metabolism , Endothelial Cells/metabolismABSTRACT
BACKGROUND: Proteins of the TGFß family, which are largely studied as homodimers, are also known to form heterodimers with biological activity distinct from their component homodimers. For instance, heterodimers of bone morphogenetic proteins, including BMP2/BMP7, BMP2/BMP6, and BMP9/BMP10, among others, have illustrated the importance of these heterodimeric proteins within the context of TGFß signaling. RESULTS: In this study, we have determined that mature GDF5 can be combined with mature BMP2 or BMP4 to form BMP2/GDF5 and BMP4/GDF5 heterodimer. Intriguingly, this combination of a BMP2 or BMP4 monomer, which exhibit high affinity to heparan sulfate characteristic to the BMP class, with a GDF5 monomer with low heparan sulfate affinity produces a heterodimer with an intermediate affinity. Using heparin affinity chromatography to purify the heterodimeric proteins, we then determined that both the BMP2/GDF5 and BMP4/GDF5 heterodimers consistently signaled potently across an array of cellular and in vivo systems, while the activities of their homodimeric counterparts were more context dependent. These differences were likely driven by an increase in the combined affinities for the type 1 receptors, Alk3 and Alk6. Furthermore, the X-ray crystal structure of BMP2/GDF5 heterodimer was determined, highlighting the formation of two asymmetric type 1 receptor binding sites that are both unique relative to the homodimers. CONCLUSIONS: Ultimately, this method of heterodimer production yielded a signaling molecule with unique properties relative to the homodimeric ligands, including high affinity to multiple type 1 and moderate heparan binding affinity.
Subject(s)
Bone Morphogenetic Protein 2 , Bone Morphogenetic Proteins , Bone Morphogenetic Proteins/metabolism , Bone Morphogenetic Protein 2/genetics , Bone Morphogenetic Protein 2/metabolism , Transforming Growth Factor beta/metabolism , Protein Binding , Carrier Proteins/metabolism , Heparitin SulfateABSTRACT
Xenbase (https://www.xenbase.org/), the Xenopus model organism knowledgebase, is a web-accessible resource that integrates the diverse genomic and biological data from research on the laboratory frogs Xenopus laevis and Xenopus tropicalis. The goal of Xenbase is to accelerate discovery and empower Xenopus research, to enhance the impact of Xenopus research data, and to facilitate the dissemination of these data. Xenbase also enhances the value of Xenopus data through high-quality curation, data integration, providing bioinformatics tools optimized for Xenopus experiments, and linking Xenopus data to human data, and other model organisms. Xenbase also plays an indispensable role in making Xenopus data interoperable and accessible to the broader biomedical community in accordance with FAIR principles. Xenbase provides annotated data updates to organizations such as NCBI, UniProtKB, Ensembl, the Gene Ontology consortium, and most recently, the Alliance of Genomic Resources, a common clearing house for data from humans and model organisms. This article provides a brief overview of key and recently added features of Xenbase. New features include processing of Xenopus high-throughput sequencing data from the NCBI Gene Expression Omnibus; curation of anatomical, physiological, and expression phenotypes with the newly created Xenopus Phenotype Ontology; Xenopus Gene Ontology annotations; new anatomical drawings of the Normal Table of Xenopus development; and integration of the latest Xenopus laevis v10.1 genome annotations. Finally, we highlight areas for future development at Xenbase as we continue to support the Xenopus research community.
Subject(s)
Databases, Genetic , Genomics , Animals , Humans , Xenopus laevis/genetics , Xenopus/genetics , Computational BiologyABSTRACT
Human intestinal organoids (HIOs) derived from pluripotent stem cells provide a valuable model for investigating human intestinal organogenesis and physiology, but they lack the immune components required to fully recapitulate the complexity of human intestinal biology and diseases. To address this issue and to begin to decipher human intestinal-immune crosstalk during development, we generated HIOs containing immune cells by transplanting HIOs under the kidney capsule of mice with a humanized immune system. We found that human immune cells temporally migrate to the mucosa and form cellular aggregates that resemble human intestinal lymphoid follicles. Moreover, after microbial exposure, epithelial microfold cells are increased in number, leading to immune cell activation determined by the secretion of IgA antibodies in the HIO lumen. This in vivo HIO system with human immune cells provides a framework for future studies on infection- or allergen-driven intestinal diseases.
Subject(s)
Pluripotent Stem Cells , Transplants , Humans , Animals , Mice , Intestines , Intestinal Mucosa , OrganoidsABSTRACT
The larynx enables speech while regulating swallowing and respiration. Larynx function hinges on the laryngeal epithelium which originates as part of the anterior foregut and undergoes extensive remodeling to separate from the esophagus and form vocal folds that interface with the adjacent trachea. Here we find that sonic hedgehog (SHH) is essential for epithelial integrity in the mouse larynx as well as the anterior foregut. During larynx-esophageal separation, low Shh expression marks specific domains of actively remodeling epithelium that undergo an epithelial-to-mesenchymal transition (EMT) characterized by the induction of N-Cadherin and movement of cells out of the epithelial layer. Consistent with a role for SHH signaling in regulating this process, Shh mutants undergo an abnormal EMT throughout the anterior foregut and larynx, marked by a cadherin switch, movement out of the epithelial layer and cell death. Unexpectedly, Shh mutant epithelial cells are replaced by a new population of FOXA2-negative cells that likely derive from adjacent pouch tissues and form a rudimentary epithelium. These findings have important implications for interpreting the etiology of HH-dependent birth defects within the foregut. We propose that SHH signaling has a default role in maintaining epithelial identity throughout the anterior foregut and that regionalized reductions in SHH trigger epithelial remodeling.
Subject(s)
Hedgehog Proteins , Larynx , Animals , Mice , Morphogenesis , Epithelial Cells , Epithelium , CadherinsABSTRACT
Genotype-phenotype associations for common diseases are often compounded by pleiotropy and metabolic state. Here, we devised a pooled human organoid-panel of steatohepatitis to investigate the impact of metabolic status on genotype-phenotype association. En masse population-based phenotypic analysis under insulin insensitive conditions predicted key non-alcoholic steatohepatitis (NASH)-genetic factors including the glucokinase regulatory protein (GCKR)-rs1260326:C>T. Analysis of NASH clinical cohorts revealed that GCKR-rs1260326-T allele elevates disease severity only under diabetic state but protects from fibrosis under non-diabetic states. Transcriptomic, metabolomic, and pharmacological analyses indicate significant mitochondrial dysfunction incurred by GCKR-rs1260326, which was not reversed with metformin. Uncoupling oxidative mechanisms mitigated mitochondrial dysfunction and permitted adaptation to increased fatty acid supply while protecting against oxidant stress, forming a basis for future therapeutic approaches for diabetic NASH. Thus, "in-a-dish" genotype-phenotype association strategies disentangle the opposing roles of metabolic-associated gene variant functions and offer a rich mechanistic, diagnostic, and therapeutic inference toolbox toward precision hepatology. VIDEO ABSTRACT.
Subject(s)
Genetic Predisposition to Disease , Non-alcoholic Fatty Liver Disease , Humans , Non-alcoholic Fatty Liver Disease/genetics , Organoids , Genetic Association Studies , Alleles , LiverABSTRACT
Ventx2 is an Antennapedia superfamily/NK-like subclass homeodomain transcription factor best known for its roles in the regulation of early dorsoventral patterning during Xenopus gastrulation and in the maintenance of neural crest multipotency. In this work we characterize the spatiotemporal expression pattern of ventx2 in progenitor cells of the Xenopus respiratory system epithelium. We find that ventx2 is directly induced by BMP signaling in the ventral foregut prior to nkx2-1, the earliest epithelial marker of the respiratory lineage. Functional studies demonstrate that Ventx2 regulates the number of Nkx2-1/Sox9+ respiratory progenitor cells induced during foregut development, the timing and level of surfactant protein gene expression, and proper tracheal-esophageal separation. Our data suggest that Ventx2 regulates the balance of respiratory progenitor cell expansion and differentiation. While the ventx gene family has been lost from the mouse genome during evolution, humans have retained a ventx2-like gene (VENTX). Finally, we discuss how our findings might suggest a possible function of VENTX in human respiratory progenitor cells.
Subject(s)
DNA-Binding Proteins , Transcription Factors , Animals , Cell Differentiation/genetics , DNA-Binding Proteins/genetics , Gene Expression Regulation, Developmental , Homeodomain Proteins/genetics , Homeodomain Proteins/metabolism , Humans , Lung/metabolism , Mice , Stem Cells/metabolism , Surface-Active Agents , Trachea , Transcription Factors/genetics , Transcription Factors/metabolism , Xenopus Proteins/genetics , Xenopus Proteins/metabolism , Xenopus laevis/geneticsABSTRACT
WNT/ß-catenin signaling controls gene expression across biological contexts from development and stem cell homeostasis to diseases including cancer. How ß-catenin is recruited to distinct enhancers to activate context-specific transcription is unclear, given that most WNT/ß-catenin-responsive transcription is thought to be mediated by TCF/LEF transcription factors (TFs). With time-resolved multi-omic analyses, we show that SOX TFs can direct lineage-specific WNT-responsive transcription during the differentiation of human pluripotent stem cells (hPSCs) into definitive endoderm and neuromesodermal progenitors. We demonstrate that SOX17 and SOX2 are required to recruit ß-catenin to lineage-specific WNT-responsive enhancers, many of which are not occupied by TCFs. At TCF-independent enhancers, SOX TFs establish a permissive chromatin landscape and recruit a WNT-enhanceosome complex to activate SOX/ß-catenin-dependent transcription. Given that SOX TFs and the WNT pathway are critical for specification of most cell types, these results have broad mechanistic implications for the specificity of WNT responses across developmental and disease contexts.
Subject(s)
Pluripotent Stem Cells , beta Catenin , Humans , Pluripotent Stem Cells/metabolism , SOX Transcription Factors/genetics , SOX Transcription Factors/metabolism , TCF Transcription Factors/genetics , Wnt Proteins/metabolism , Wnt Signaling Pathway , beta Catenin/metabolismABSTRACT
Development of visceral organs such as the esophagus, lung, liver and stomach are coordinated by reciprocal signaling interactions between the endoderm and adjacent mesoderm cells in the fetal foregut. Although the recent successes in recapitulating developmental signaling in vitro has enabled the differentiation of human pluripotent stem cells (hPSCs) into various types of organ-specific endodermal epithelium, the generation of organ-specific mesenchyme has received much less attention. This is a major limitation in ongoing efforts to engineer complex human tissue. Here, we describe a protocol to differentiate hPSCs into different types of organ-specific mesoderm, leveraging signaling networks and molecular markers elucidated from single-cell transcriptomics of mouse foregut organogenesis. Building on established methods, hPSC-derived lateral plate mesoderm treated with either retinoic acid (RA) or RA together with a Hedgehog (HH) agonist generates posterior or anterior foregut splanchnic mesoderm, respectively, after 4-d cultures. These are directed into organ-specific mesenchyme lineages by the combinatorial activation or inhibition of WNT, BMP, RA or HH pathways from days 4 to 7 in cultures. By day 7, the cultures are enriched for different types of mesoderm with distinct molecular signatures: 60-90% pure liver septum transversum/mesothelium-like, 70-80% pure liver-like fibroblasts and populations of ~35% respiratory-like mesoderm, gastric-like mesoderm or esophageal-like mesoderm. This protocol can be performed by anyone with moderate experience differentiating hPSCs, provides a novel platform to study human mesoderm development and can be used to engineer more complex foregut tissue for disease modeling and regenerative medicine.
Subject(s)
Hedgehog Proteins , Pluripotent Stem Cells , Humans , Mice , Animals , Hedgehog Proteins/metabolism , Mesoderm , Endoderm , Cell Differentiation , Tretinoin/pharmacology , LungABSTRACT
Normal tables of development are essential for studies of embryogenesis, serving as an important resource for model organisms, including the frog Xenopus laevis. Xenopus has long been used to study developmental and cell biology, and is an increasingly important model for human birth defects and disease, genomics, proteomics and toxicology. Scientists utilize Nieuwkoop and Faber's classic 'Normal Table of Xenopus laevis (Daudin)' and accompanying illustrations to enable experimental reproducibility and reuse the illustrations in new publications and teaching. However, it is no longer possible to obtain permission for these copyrighted illustrations. We present 133 new, high-quality illustrations of X. laevis development from fertilization to metamorphosis, with additional views that were not available in the original collection. All the images are available on Xenbase, the Xenopus knowledgebase (http://www.xenbase.org/entry/zahn.do), for download and reuse under an attributable, non-commercial creative commons license. Additionally, we have compiled a 'Landmarks Table' of key morphological features and marker gene expression that can be used to distinguish stages quickly and reliably (https://www.xenbase.org/entry/landmarks-table.do). This new open-access resource will facilitate Xenopus research and teaching in the decades to come.
Subject(s)
Databases, Genetic , Genomics , Animals , Humans , Metamorphosis, Biological , Reproducibility of Results , Xenopus laevis/geneticsABSTRACT
[This corrects the article DOI: 10.1016/j.xhgg.2022.100107.].
ABSTRACT
Esophageal atresias/tracheoesophageal fistulas (EA/TEF) are rare congenital anomalies caused by aberrant development of the foregut. Previous studies indicate that rare or de novo genetic variants significantly contribute to EA/TEF risk, and most individuals with EA/TEF do not have pathogenic genetic variants in established risk genes. To identify the genetic contributions to EA/TEF, we performed whole genome sequencing of 185 trios (probands and parents) with EA/TEF, including 59 isolated and 126 complex cases with additional congenital anomalies and/or neurodevelopmental disorders. There was a significant burden of protein-altering de novo coding variants in complex cases (p = 3.3 × 10-4), especially in genes that are intolerant of loss-of-function variants in the population. We performed simulation analysis of pathway enrichment based on background mutation rate and identified a number of pathways related to endocytosis and intracellular trafficking that as a group have a significant burden of protein-altering de novo variants. We assessed 18 variants for disease causality using CRISPR-Cas9 mutagenesis in Xenopus and confirmed 13 with tracheoesophageal phenotypes. Our results implicate disruption of endosome-mediated epithelial remodeling as a potential mechanism of foregut developmental defects. Our results suggest significant genetic heterogeneity of EA/TEF and may have implications for the mechanisms of other rare congenital anomalies.
ABSTRACT
BACKGROUND: Ontologies of precisely defined, controlled vocabularies are essential to curate the results of biological experiments such that the data are machine searchable, can be computationally analyzed, and are interoperable across the biomedical research continuum. There is also an increasing need for methods to interrelate phenotypic data easily and accurately from experiments in animal models with human development and disease. RESULTS: Here we present the Xenopus phenotype ontology (XPO) to annotate phenotypic data from experiments in Xenopus, one of the major vertebrate model organisms used to study gene function in development and disease. The XPO implements design patterns from the Unified Phenotype Ontology (uPheno), and the principles outlined by the Open Biological and Biomedical Ontologies (OBO Foundry) to maximize interoperability with other species and facilitate ongoing ontology management. Constructed in Web Ontology Language (OWL) the XPO combines the existing uPheno library of ontology design patterns with additional terms from the Xenopus Anatomy Ontology (XAO), the Phenotype and Trait Ontology (PATO) and the Gene Ontology (GO). The integration of these different ontologies into the XPO enables rich phenotypic curation, whilst the uPheno bridging axioms allows phenotypic data from Xenopus experiments to be related to phenotype data from other model organisms and human disease. Moreover, the simple post-composed uPheno design patterns facilitate ongoing XPO development as the generation of new terms and classes of terms can be substantially automated. CONCLUSIONS: The XPO serves as an example of current best practices to help overcome many of the inherent challenges in harmonizing phenotype data between different species. The XPO currently consists of approximately 22,000 terms and is being used to curate phenotypes by Xenbase, the Xenopus Model Organism Knowledgebase, forming a standardized corpus of genotype-phenotype data that can be directly related to other uPheno compliant resources.
Subject(s)
Biological Ontologies , Animals , Gene Ontology , Humans , Phenotype , Xenopus laevisABSTRACT
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.