Integrated microarray and ChIP analysis identifies multiple Foxa2 dependent target genes in the notochord

A single-cell atlas of spatial and temporal gene expression in the mouse cranial neural plate

The formation of the mammalian brain requires regionalization and morphogenesis of the cranial neural plate, which transforms from an epithelial sheet into a closed tube that provides the structural foundation for neural patterning and circuit formation. Sonic hedgehog (SHH) signaling is important for cranial neural plate patterning and closure, but the transcriptional changes that give rise to the spatially regulated cell fates and behaviors that build the cranial neural tube have not been systematically analyzed. Here, we used single-cell RNA sequencing to generate an atlas of gene expression at six consecutive stages of cranial neural tube closure in the mouse embryo. Ordering transcriptional profiles relative to the major axes of gene expression predicted spatially regulated expression of 870 genes along the anterior-posterior and mediolateral axes of the cranial neural plate and reproduced known expression patterns with over 85% accuracy. Single-cell RNA sequencing of embryos with activated SHH signaling revealed distinct SHH-regulated transcriptional programs in the developing forebrain, midbrain, and hindbrain, suggesting a complex interplay between anterior-posterior and mediolateral patterning systems. These results define a spatiotemporally resolved map of gene expression during cranial neural tube closure and provide a resource for investigating the transcriptional events that drive early mammalian brain development.

In vitro modelling of anterior primitive streak patterning with human pluripotent stem cells identifies the path to notochord progenitors

Notochord progenitors (NotoPs) represent a scarce yet crucial embryonic cell population, playing important roles in embryo patterning and eventually giving rise to the cells that form and maintain intervertebral discs. The mechanisms regulating NotoPs emergence are unclear. This knowledge gap persists due to the inherent complexity of cell fate patterning during gastrulation, particularly within the anterior primitive streak (APS), where NotoPs first arise alongside neuro-mesoderm and endoderm. To gain insights into this process, we use micropatterning together with FGF and the WNT pathway activator CHIR9901 to guide the development of human embryonic stem cells into reproducible patterns of APS cell fates. We show that CHIR9901 dosage dictates the downstream dynamics of endogenous TGFβ signalling, which in turn controls cell fate decisions. While sustained NODAL signalling defines endoderm and NODAL inhibition is imperative for neuro-mesoderm emergence, timely inhibition of NODAL signalling with spatial confinement potentiates WNT activity and enables us to generate NotoPs efficiently. Our work elucidates the signalling regimes underpinning NotoP emergence and provides insights into the regulatory mechanisms controlling the balance of APS cell fates during gastrulation.

Dedifferentiation-like reprogramming of degenerative nucleus pulposus cells into notochordal-like cells by defined factors

The extensive degeneration of functional somatic cells and the depletion of endogenous stem/progenitor populations present significant challenges to tissue regeneration in degenerative diseases. Currently, a cellular reprogramming approach enabling directly generating corresponding progenitor populations from degenerative somatic cells remains elusive. The present study focused on intervertebral disc degeneration (IVDD) and identified a three-factor combination (OCT4, FOXA2, TBXT [OFT]) that could induce the dedifferentiation-like reprogramming of degenerative nucleus pulposus cells (dNPCs) toward induced notochordal-like cells (iNCs). Single-cell transcriptomics dissected the transitions of cell identity during reprogramming. Further, OCT4 was found to directly interact with bromodomain PHD-finger transcription factor to remodel the chromatin during the early phases, which was crucial for initiating this dedifferentiation-like reprogramming. In rat models, intradiscal injection of adeno-associated virus carrying OFT generated iNCs from in situ dNPCs and reversed IVDD. These results collectively present a proof-of-concept for dedifferentiation-like reprogramming of degenerated somatic cells into corresponding progenitors through the development of a factor-based strategy, providing a promising approach for regeneration in degenerative disc diseases.

Cis-regulatory interfaces reveal the molecular mechanisms underlying the notochord gene regulatory network of Ciona

Tissue-specific gene expression is fundamental in development and evolution, and is mediated by transcription factors and by the cis-regulatory regions (enhancers) that they control. Transcription factors and their respective tissue-specific enhancers are essential components of gene regulatory networks responsible for the development of tissues and organs. Although numerous transcription factors have been characterized from different organisms, the knowledge of the enhancers responsible for their tissue-specific expression remains fragmentary. Here we use Ciona to study the enhancers associated with ten transcription factors expressed in the notochord, an evolutionary hallmark of the chordate phylum. Our results illustrate how two evolutionarily conserved transcription factors, Brachyury and Foxa2, coordinate the deployment of other notochord transcription factors. The results of these detailed cis-regulatory analyses delineate a high-resolution view of the essential notochord gene regulatory network of Ciona, and provide a reference for studies of transcription factors, enhancers, and their roles in development, disease, and evolution.

In vitro and in vivo models define a molecular signature reference for human embryonic notochordal cells

Understanding the emergence of human notochordal cells (NC) is essential for the development of regenerative approaches. We present a comprehensive investigation into the specification and generation of bona fide NC using a straightforward pluripotent stem cell (PSC)-based system benchmarked with human fetal notochord. By integrating in vitro and in vivo transcriptomic data at single-cell resolution, we establish an extended molecular signature and overcome the limitations associated with studying human notochordal lineage at early developmental stages. We show that TGF-β inhibition enhances the yield and homogeneity of notochordal lineage commitment in vitro. Furthermore, this study characterizes regulators of cell-fate decision and matrisome enriched in the notochordal niche. Importantly, we identify specific cell-surface markers opening avenues for differentiation refinement, NC purification, and functional studies. Altogether, this study provides a human notochord transcriptomic reference that will serve as a resource for notochord identification in human systems, diseased-tissues modeling, and facilitating future biomedical research.

Genome-wide identification of notochord enhancers comprising the regulatory landscape of the brachyury locus in mouse

The node and notochord are important signaling centers organizing the dorso-ventral patterning of cells arising from neuro-mesodermal progenitors forming the embryonic body anlage. Owing to the scarcity of notochord progenitors and notochord cells, a comprehensive identification of regulatory elements driving notochord-specific gene expression has been lacking. Here, we have used ATAC-seq analysis of FACS-purified notochord cells from Theiler stage 12-13 mouse embryos to identify 8921 putative notochord enhancers. In addition, we established a new model for generating notochord-like cells in culture, and found 3728 of these enhancers occupied by the essential notochord control factors brachyury (T) and/or Foxa2. We describe the regulatory landscape of the T locus, comprising ten putative enhancers occupied by these factors, and confirmed the regulatory activity of three of these elements. Moreover, we characterized seven new elements by knockout analysis in embryos and identified one new notochord enhancer, termed TNE2. TNE2 cooperates with TNE in the trunk notochord, and is essential for notochord differentiation in the tail. Our data reveal an essential role of Foxa2 in directing T-expressing cells towards the notochord lineage.

Nodal flow transfers polycystin to determine mouse left-right asymmetry

Left-dominant [Ca2+]i elevation on the left margin of the ventral node furnishes the initial laterality of mouse embryos. It depends on extracellular leftward fluid flow (nodal flow), fibroblast growth factor receptor (FGFR)/sonic hedgehog (Shh) signaling, and the PKD1L1 polycystin subunit, of which interrelationship is still elusive. Here, we show that leftward nodal flow directs PKD1L1-containing fibrous strands and facilitates Nodal-mediated [Ca2+]i elevation on the left margin. We generate KikGR-PKD1L1 knockin mice in order to monitor protein dynamics with a photoconvertible fluorescence protein tag. By imaging those embryos, we have identified fragile meshwork being gradually transferred leftward involving pleiomorphic extracellular events. A portion of the meshwork finally bridges over the left nodal crown cells in an FGFR/Shh-dependent manner. As PKD1L1 N-term is predominantly associated with Nodal on the left margin and that PKD1L1/PKD2 overexpression significantly augments cellular Nodal sensitivity, we propose that leftward transfer of polycystin-containing fibrous strands determines left-right asymmetry in developing embryos.

Next-generation plasmids for transgenesis in zebrafish and beyond

Transgenesis is an essential technique for any genetic model. Tol2-based transgenesis paired with Gateway-compatible vector collections has transformed zebrafish transgenesis with an accessible modular system. Here, we establish several next-generation transgenesis tools for zebrafish and other species to expand and enhance transgenic applications. To facilitate gene regulatory element testing, we generated Gateway middle entry vectors harboring the small mouse beta-globin minimal promoter coupled to several fluorophores, CreERT2 and Gal4. To extend the color spectrum for transgenic applications, we established middle entry vectors encoding the bright, blue-fluorescent protein mCerulean and mApple as an alternative red fluorophore. We present a series of p2A peptide-based 3' vectors with different fluorophores and subcellular localizations to co-label cells expressing proteins of interest. Finally, we established Tol2 destination vectors carrying the zebrafish exorh promoter driving different fluorophores as a pineal gland-specific transgenesis marker that is active before hatching and through adulthood. exorh-based reporters and transgenesis markers also drive specific pineal gland expression in the eye-less cavefish (Astyanax). Together, our vectors provide versatile reagents for transgenesis applications in zebrafish, cavefish and other models.

Xbp1 and Brachyury establish an evolutionarily conserved subcircuit of the notochord gene regulatory network

Gene regulatory networks coordinate the formation of organs and structures that compose the evolving body plans of different organisms. We are using a simple chordate model, the Ciona embryo, to investigate the essential gene regulatory network that orchestrates morphogenesis of the notochord, a structure necessary for the proper development of all chordate embryos. Although numerous transcription factors expressed in the notochord have been identified in different chordates, several of them remain to be positioned within a regulatory framework. Here, we focus on Xbp1, a transcription factor expressed during notochord formation in Ciona and other chordates. Through the identification of Xbp1-downstream notochord genes in Ciona, we found evidence of the early co-option of genes involved in the unfolded protein response to the notochord developmental program. We report the regulatory interplay between Xbp1 and Brachyury, and by extending these results to Xenopus, we show that Brachyury and Xbp1 form a cross-regulatory subcircuit of the notochord gene regulatory network that has been consolidated during chordate evolution.

Functional in vivo characterization of sox10 enhancers in neural crest and melanoma development

The role of a neural crest developmental transcriptional program, which critically involves Sox10 upregulation, is a key conserved aspect of melanoma initiation in both humans and zebrafish, yet transcriptional regulation of sox10 expression is incompletely understood. Here we used ATAC-Seq analysis of multiple zebrafish melanoma tumors to identify recurrently open chromatin domains as putative melanoma-specific sox10 enhancers. Screening in vivo with EGFP reporter constructs revealed 9 of 11 putative sox10 enhancers with embryonic activity in zebrafish. Focusing on the most active enhancer region in melanoma, we identified a region 23 kilobases upstream of sox10, termed peak5, that drives EGFP reporter expression in a subset of neural crest cells, Kolmer-Agduhr neurons, and early melanoma patches and tumors with high specificity. A ~200 base pair region, conserved in Cyprinidae, within peak5 is required for transgenic reporter activity in neural crest and melanoma. This region contains dimeric SoxE/Sox10 dimeric binding sites essential for peak5 neural crest and melanoma activity. We show that deletion of the endogenous peak5 conserved genomic locus decreases embryonic sox10 expression and disrupts adult stripe patterning in our melanoma model background. Our work demonstrates the power of linking developmental and cancer models to better understand neural crest identity in melanoma.

Long-Range Optogenetic Control of Axon Guidance Overcomes Developmental Boundaries and Defects

Axons connect neurons together, establishing the wiring architecture of neuronal networks. Axonal connectivity is largely built during embryonic development through highly constrained processes of axon guidance, which have been extensively studied. However, the inability to control axon guidance, and thus neuronal network architecture, has limited investigation of how axonal connections influence subsequent development and function of neuronal networks. Here, we use zebrafish motor neurons expressing a photoactivatable Rac1 to co-opt endogenous growth cone guidance machinery to precisely and non-invasively direct axon growth using light. Axons can be guided over large distances, within complex environments of living organisms, overriding competing endogenous signals and redirecting axons across potent repulsive barriers to construct novel circuitry. Notably, genetic axon guidance defects can be rescued, restoring functional connectivity. These data demonstrate that intrinsic growth cone guidance machinery can be co-opted to non-invasively build new connectivity, allowing investigation of neural network dynamics in intact living organisms.

Embryonic ethanol exposure alters expression of sox2 and other early transcripts in zebrafish, producing gastrulation defects

Ethanol exposure during prenatal development causes fetal alcohol spectrum disorder (FASD), the most frequent preventable birth defect and neurodevelopmental disability syndrome. The molecular targets of ethanol toxicity during development are poorly understood. Developmental stages surrounding gastrulation are very sensitive to ethanol exposure. To understand the effects of ethanol on early transcripts during embryogenesis, we treated zebrafish embryos with ethanol during pre-gastrulation period and examined the transcripts by Affymetrix GeneChip microarray before gastrulation. We identified 521 significantly dysregulated genes, including 61 transcription factors in ethanol-exposed embryos. Sox2, the key regulator of pluripotency and early development was significantly reduced. Functional annotation analysis showed enrichment in transcription regulation, embryonic axes patterning, and signaling pathways, including Wnt, Notch and retinoic acid. We identified all potential genomic targets of 25 dysregulated transcription factors and compared their interactions with the ethanol-dysregulated genes. This analysis predicted that Sox2 targeted a large number of ethanol-dysregulated genes. A gene regulatory network analysis showed that many of the dysregulated genes are targeted by multiple transcription factors. Injection of sox2 mRNA partially rescued ethanol-induced gene expression, epiboly and gastrulation defects. Additional studies of this ethanol dysregulated network may identify therapeutic targets that coordinately regulate early development.

NOTO Transcription Factor Directs Human Induced Pluripotent Stem Cell-Derived Mesendoderm Progenitors to a Notochordal Fate

The founder cells of the Nucleus pulposus, the centre of the intervertebral disc, originate in the embryonic notochord. After birth, mature notochordal cells (NC) are identified as key regulators of disc homeostasis. Better understanding of their biology has great potential in delaying the onset of disc degeneration or as a regenerative-cell source for disc repair. Using human pluripotent stem cells, we developed a two-step method to generate a stable NC-like population with a distinct molecular signature. Time-course analysis of lineage-specific markers shows that WNT pathway activation and transfection of the notochord-related transcription factor NOTO are sufficient to induce high levels of mesendoderm progenitors and favour their commitment toward the notochordal lineage instead of paraxial and lateral mesodermal or endodermal lineages. This study results in the identification of NOTO-regulated genes including some that are found expressed in human healthy disc tissue and highlights NOTO function in coordinating the gene network to human notochord differentiation.

The notochord gene regulatory network in chordate evolution: Conservation and divergence from Ciona to vertebrates

The notochord is a structure required for support and patterning of all chordate embryos, from sea squirts to humans. An increasing amount of information on notochord development and on the molecular strategies that ensure its proper morphogenesis has been gleaned through studies in the sea squirt Ciona. This invertebrate chordate offers a fortunate combination of experimental advantages, ranging from translucent, fast-developing embryos to a compact genome and impressive biomolecular resources. These assets have enabled the rapid identification of numerous notochord genes and cis-regulatory regions, and provide a rather unique opportunity to reconstruct the gene regulatory network that controls the formation of this developmental and evolutionary chordate landmark. This chapter summarizes the morphogenetic milestones that punctuate notochord formation in Ciona, their molecular effectors, and the current knowledge of the gene regulatory network that ensures the accurate spatial and temporal orchestration of these processes.

Active receptor tyrosine kinases, but not Brachyury, are sufficient to trigger chordoma in zebrafish

The aberrant activation of developmental processes triggers diverse cancer types. Chordoma is a rare, aggressive tumor arising from transformed notochord remnants. Several potentially oncogenic factors have been found to be deregulated in chordoma, yet causation remains uncertain. In particular, sustained expression of TBXT – encoding the notochord regulator protein brachyury – is hypothesized as a key driver of chordoma, yet experimental evidence is absent. Here, we employ a zebrafish chordoma model to identify the notochord-transforming potential of implicated genes in vivo. We find that Brachyury, including a form with augmented transcriptional activity, is insufficient to initiate notochord hyperplasia. In contrast, the chordoma-implicated receptor tyrosine kinases (RTKs) EGFR and Kdr/VEGFR2 are sufficient to transform notochord cells. Aberrant activation of RTK/Ras signaling attenuates processes required for notochord differentiation, including the unfolded protein response and endoplasmic reticulum stress pathways. Our results provide the first in vivo evidence against a tumor-initiating potential of Brachyury in the notochord, and imply activated RTK signaling as a possible initiating event in chordoma. Furthermore, our work points at modulating endoplasmic reticulum and protein stress pathways as possible therapeutic avenues against chordoma.

Summary: An injection-based chordoma model in zebrafish shows that the hypothesized chordoma oncogene brachyury is insufficient, whereas EGFR and VEGFR2 are sufficient, to trigger notochord hyperplasia in our model.

Bioinformatic Analysis of Gene Variants from Gastroschisis Recurrence Identifies Multiple Novel Pathogenetic Pathways: Implication for the Closure of the Ventral Body Wall

We investigated whether likely pathogenic variants co-segregating with gastroschisis through a family-based approach using bioinformatic analyses were implicated in body wall closure. Gene Ontology (GO)/Panther functional enrichment and protein-protein interaction analysis by String identified several biological networks of highly connected genes in UGT1A3, UGT1A4, UGT1A5, UGT1A6, UGT1A7, UGT1A8, UGT1A9, UGT1A10, AOX1, NOTCH1, HIST1H2BB, RPS3, THBS1, ADCY9, and FGFR4. SVS–PhoRank identified a dominant model in OR10G4 (also as heterozygous de novo), ITIH3, PLEKHG4B, SLC9A3, ITGA2, AOX1, and ALPP, including a recessive model in UGT1A7, UGT1A6, PER2, PTPRD, and UGT1A3. A heterozygous compound model was observed in CDYL, KDM5A, RASGRP1, MYBPC2, PDE4DIP, F5, OBSCN, and UGT1A. These genes were implicated in pathogenetic pathways involving the following GO related categories: xenobiotic, regulation of metabolic process, regulation of cell adhesion, regulation of gene expression, inflammatory response, regulation of vascular development, keratinization, left-right symmetry, epigenetic, ubiquitination, and regulation of protein synthesis. Multiple background modifiers interacting with disease-relevant pathways may regulate gastroschisis susceptibility. Based in our findings and considering the plausibility of the biological pattern of mechanisms and gene network modeling, we suggest that the gastroschisis developmental process may be the consequence of several well-orchestrated biological and molecular mechanisms which could be interacting with gastroschisis predispositions within the first ten weeks of development.

The dyslexia susceptibility KIAA0319 gene shows a specific expression pattern during zebrafish development supporting a role beyond neuronal migration

Dyslexia is a common neurodevelopmental disorder caused by a significant genetic component. The KIAA0319 gene is one of the most robust dyslexia susceptibility factors but its function remains poorly understood. Initial RNA-interference studies in rats suggested a role in neuronal migration whereas subsequent work with double knock-out mouse models for both Kiaa0319 and its paralogue Kiaa0319-like reported effects in the auditory system but not in neuronal migration. To further understand the role of KIAA0319 during neurodevelopment, we carried out an expression study of its zebrafish orthologue at different embryonic stages. We used different approaches including RNAscope in situ hybridization combined with light-sheet microscopy. The results show particularly high expression during the first few hours of development. Later, expression becomes localized in well-defined structures. In addition to high expression in the brain, we report for the first time expression in the eyes and the notochord. Surprisingly, kiaa0319-like, which generally shows a similar expression pattern to kiaa0319, was not expressed in the notochord suggesting a distinct role for kiaa0319 in this structure. This observation was supported by the identification of notochord enhancers enriched upstream of the KIAA0319 transcription start site, in both zebrafish and humans. This study supports a developmental role for KIAA0319 in the brain as well as in other developing structures, particularly in the notochord which, is key for establishing body patterning in vertebrates.

Transcriptionally dynamic progenitor populations organised around a stable niche drive axial patterning

The elongating mouse anteroposterior axis is supplied by progenitors with distinct tissue fates. It is not known whether these progenitors confer anteroposterior pattern to the embryo. We have analysed the progenitor population transcriptomes in the mouse primitive streak and tail bud throughout axial elongation. Transcriptomic signatures distinguish three known progenitor types (neuromesodermal, lateral/paraxial mesoderm and notochord progenitors; NMPs, LPMPs and NotoPs). Both NMP and LPMP transcriptomes change extensively over time. In particular, NMPs upregulate Wnt, Fgf and Notch signalling components, and many Hox genes as progenitors transit from production of the trunk to the tail and expand in number. In contrast, the transcriptome of NotoPs is stable throughout axial elongation and they are required for normal axis elongation. These results suggest that NotoPs act as a progenitor niche whereas anteroposterior patterning originates within NMPs and LPMPs.

Development of the axial skeleton and intervertebral disc

Development of the axial skeleton is a complex, stepwise process that relies on intricate signaling and coordinated cellular differentiation. Disruptions to this process can result in a myriad of skeletal malformations that range in severity. The notochord and the sclerotome are embryonic tissues that give rise to the major components of the intervertebral discs and the vertebral bodies of the spinal column. Through a number of mouse models and characterization of congenital abnormalities in human patients, various growth factors, transcription factors, and other signaling proteins have been demonstrated to have critical roles in the development of the axial skeleton. Balance between opposing growth factors as well as other environmental cues allows for cell fate specification and divergence of tissue types during development. Furthermore, characterization of progenitor cells for specific cell lineages has furthered the understanding of specific spatiotemporal cues that cells need in order to initiate and complete development of distinct tissues. Identifying specific marker genes that can distinguish between the various embryonic and mature cell types is also of importance. Clinically, understanding developmental clues can aid in the generation of therapeutics for musculoskeletal disease through the process of developmental engineering. Studies into potential stem cell therapies are based on knowledge of the normal processes that occur in the embryo, which can then be applied to stepwise tissue engineering strategies.

Latest advances in intervertebral disc development and progenitor cells

This paper is a concise review aiming to assemble the most relevant topics presented by the authors at ORS-Philadelphia Spine Research Society Fourth International Spine Research Symposium. It centers on the latest advances in disc development, its main structural entities, and the populating cells, with emphasis on the advances in pivotal molecular pathways responsible for forming the intervertebral discs (IVD). The objective of finding and emphasizing pathways and mechanisms that function to control tissue formation is to identify and to explore modifications occurring during normal aging, disease, and tissue repair. Thus, to comprehend that the cellular and molecular basis of tissue degeneration are crucial in the study of the dynamic interplay that includes cell-cell communication, gene regulation, and growth factors required to form a healthy and functional tissue during normal development.

FACS-Seq analysis of Pax3-derived cells identifies non-myogenic lineages in the embryonic forelimb

Skeletal muscle in the forelimb develops during embryonic and fetal development and perinatally. While much is known regarding the molecules involved in forelimb myogenesis, little is known about the specific mechanisms and interactions. Migrating skeletal muscle precursor cells express Pax3 as they migrate into the forelimb from the dermomyotome. To compare gene expression profiles of the same cell population over time, we isolated lineage-traced Pax3⁺ cells (Pax3^EGFP) from forelimbs at different embryonic days. We performed whole transcriptome profiling via RNA-Seq of Pax3⁺ cells to construct gene networks involved in different stages of embryonic and fetal development. With this, we identified genes involved in the skeletal, muscular, vascular, nervous and immune systems. Expression of genes related to the immune, skeletal and vascular systems showed prominent increases over time, suggesting a non-skeletal myogenic context of Pax3-derived cells. Using co-expression analysis, we observed an immune-related gene subnetwork active during fetal myogenesis, further implying that Pax3-derived cells are not a strictly myogenic lineage, and are involved in patterning and three-dimensional formation of the forelimb through multiple systems.

Synergistic co-regulation and competition by a SOX9-GLI-FOXA phasic transcriptional network coordinate chondrocyte differentiation transitions

The growth plate mediates bone growth where SOX9 and GLI factors control chondrocyte proliferation, differentiation and entry into hypertrophy. FOXA factors regulate hypertrophic chondrocyte maturation. How these factors integrate into a Gene Regulatory Network (GRN) controlling these differentiation transitions is incompletely understood. We adopted a genome-wide whole tissue approach to establish a Growth Plate Differential Gene Expression Library (GP-DGEL) for fractionated proliferating, pre-hypertrophic, early and late hypertrophic chondrocytes, as an overarching resource for discovery of pathways and disease candidates. De novo motif discovery revealed the enrichment of SOX9 and GLI binding sites in the genes preferentially expressed in proliferating and prehypertrophic chondrocytes, suggesting the potential cooperation between SOX9 and GLI proteins. We integrated the analyses of the transcriptome, SOX9, GLI1 and GLI3 ChIP-seq datasets, with functional validation by transactivation assays and mouse mutants. We identified new SOX9 targets and showed SOX9-GLI directly and cooperatively regulate many genes such as Trps1, Sox9, Sox5, Sox6, Col2a1, Ptch1, Gli1 and Gli2. Further, FOXA2 competes with SOX9 for the transactivation of target genes. The data support a model of SOX9-GLI-FOXA phasic GRN in chondrocyte development. Together, SOX9-GLI auto-regulate and cooperate to activate and repress genes in proliferating chondrocytes. Upon hypertrophy, FOXA competes with SOX9, and control toward terminal differentiation passes to FOXA, RUNX, AP1 and MEF2 factors.

In the development of the mammalian growth plate, while several transcription factors are individually well known for their key roles in regulating phases of chondrocyte differentiation, there is little information on how they interact and cooperate with each other. We took an unbiased genome wide approach to identify the transcription factors and signaling pathways that play dominant roles in the chondrocyte differentiation cascade. We developed a searchable library of differentially expressed genes, GP-DGEL, which has fine spatial resolution and global transcriptomic coverage for discovery of processes, pathways and disease candidates. Our work identifies a novel regulatory mechanism that integrates the action of three transcription factors, SOX9, GLI and FOXA. SOX9-GLI auto-regulate and cooperate to activate and repress genes in proliferating chondrocytes. Upon entry into prehypertrophy, FOXA competes with SOX9, and control of hypertrophy passes to FOXA, RUNX, AP1 and MEF2 factors.

Identification and fate mapping of the pancreatic mesenchyme

The murine pancreas buds from the ventral embryonic endoderm at approximately 8.75 dpc and a second pancreas bud emerges from the dorsal endoderm by 9.0 dpc. Although it is clear that secreted signals from adjacent mesoderm-derived sources are required for both the appropriate emergence and further refinement of the pancreatic endoderm, neither the exact signals nor the requisite tissue sources have been defined in mammalian systems. Herein we use DiI fate mapping of cultured murine embryos to identify the embryonic sources of both the early inductive and later condensed pancreatic mesenchyme. Despite being capable of supporting pancreas induction from dorsal endoderm in co-culture experiments, we find that in the context of the developing embryo, the dorsal aortae as well as the paraxial, intermediate, and lateral mesoderm derivatives only transiently associate with the dorsal pancreas bud, producing descendants that are decidedly anterior to the pancreas bud. Unlike these other mesoderm derivatives, the axial (notochord) descendants maintain association with the dorsal pre-pancreatic endoderm and early pancreas bud. This fate mapping data points to the notochord as the likely inductive source in vivo while also revealing dynamic morphogenetic movements displayed by individual mesodermal subtypes. Because none of the mesoderm examined above produced the pancreatic mesenchyme that condenses around the induced bud to support exocrine and endocrine differentiation, we also sought to identify the mesodermal origins of this mesenchyme. We identify a portion of the coelomic mesoderm that contributes to the condensed pancreatic mesenchyme. In conclusion, we identify a portion of the notochord as a likely source of the signals required to induce and maintain the early dorsal pancreas bud, demonstrate that the coelomic mesothelium contributes to the dorsal and ventral pancreatic mesenchyme, and provide insight into the dynamic morphological rearrangements of mesoderm-derived tissues during early organogenesis stages of mammalian development.

Functional and evolutionary insights from the Ciona notochord transcriptome

The notochord of the ascidian Ciona consists of only 40 cells, and is a longstanding model for studying organogenesis in a small, simple embryo. Here, we perform RNAseq on flow-sorted notochord cells from multiple stages to define a comprehensive Ciona notochord transcriptome. We identify 1364 genes with enriched expression and extensively validate the results by in situ hybridization. These genes are highly enriched for Gene Ontology terms related to the extracellular matrix, cell adhesion and cytoskeleton. Orthologs of 112 of the Ciona notochord genes have known notochord expression in vertebrates, more than twice as many as predicted by chance alone. This set of putative effector genes with notochord expression conserved from tunicates to vertebrates will be invaluable for testing hypotheses about notochord evolution. The full set of Ciona notochord genes provides a foundation for systems-level studies of notochord gene regulation and morphogenesis. We find only modest overlap between this set of notochord-enriched transcripts and the genes upregulated by ectopic expression of the key notochord transcription factor Brachyury, indicating that Brachyury is not a notochord master regulator gene as strictly defined.

An Integrative Developmental Genomics and Systems Biology Approach to Identify an In Vivo Sox Trio-Mediated Gene Regulatory Network in Murine Embryos

Embryogenesis is an intricate process involving multiple genes and pathways. Some of the key transcription factors controlling specific cell types are the Sox trio, namely, Sox5, Sox6, and Sox9, which play crucial roles in organogenesis working in a concerted manner. Much however still needs to be learned about their combinatorial roles during this process. A developmental genomics and systems biology approach offers to complement the reductionist methodology of current developmental biology and provide a more comprehensive and integrated view of the interrelationships of complex regulatory networks that occur during organogenesis. By combining cell type-specific transcriptome analysis and in vivo ChIP-Seq of the Sox trio using mouse embryos, we provide evidence for the direct control of Sox5 and Sox6 by the transcriptional trio in the murine model and by Morpholino knockdown in zebrafish and demonstrate the novel role of Tgfb2, Fbxl18, and Tle3 in formation of Sox5, Sox6, and Sox9 dependent tissues. Concurrently, a complete embryonic gene regulatory network has been generated, identifying a wide repertoire of genes involved and controlled by the Sox trio in the intricate process of normal embryogenesis.

Making a Kidney Organoid Using the Directed Differentiation of Human Pluripotent Stem Cells

An organoid can be defined as a three-dimensional organ-like structure formed from organ-specific progenitor cells. Organ progenitor cells were empirically found to self-organize three-dimensional tissues when they were aggregated and cultivated in vitro. While this nature power of progenitor cells has an amazing potential to recreate artificial organs in vitro, there had been difficulty to apply this technology to human organs due to the inaccessibility to human progenitor cells until human-induced pluripotent stem cell (hiPSC) was invented by Takahashi and Yamanaka in 2007. As embryonic stem cells do, hiPSCs also have pluripotency to give rise to any organs/tissues cell types, including the kidney, via directed differentiation. Here, we provide a detailed protocol for generating kidney organoids using human pluripotent stem cells. The protocol differentiates human pluripotent stem cells into the posterior primitive streak. This is followed by the simultaneous induction of posterior and anterior intermediate mesoderm that are subsequently aggregated and undergo self-organization into the kidney organoid. Such kidney organoids are comprised of all anticipated kidney cell types including nephrons segmented into the glomerulus, proximal tubule, loop of Henle, and distal tubule as well as the collecting duct, endothelial network, and renal interstitium.

Lhx1 functions together with Otx2, Foxa2, and Ldb1 to govern anterior mesendoderm, node, and midline development

Costello et al. demonstrate that Smad4/Eomes-dependent Lhx1 expression in the epiblast marks the entire definitive endoderm lineage, the anterior mesendoderm, and midline progenitors. In proteomic experiments, they characterize a complex comprised of Lhx1, Otx2, and Foxa2 as well as the chromatin-looping protein Ldb1.

Gene regulatory networks controlling functional activities of spatially and temporally distinct endodermal cell populations in the early mouse embryo remain ill defined. The T-box transcription factor Eomes, acting downstream from Nodal/Smad signals, directly activates the LIM domain homeobox transcription factor Lhx1 in the visceral endoderm. Here we demonstrate Smad4/Eomes-dependent Lhx1 expression in the epiblast marks the entire definitive endoderm lineage, the anterior mesendoderm, and midline progenitors. Conditional inactivation of Lhx1 disrupts anterior definitive endoderm development and impedes node and midline morphogenesis in part due to severe disturbances in visceral endoderm displacement. Transcriptional profiling and ChIP-seq (chromatin immunoprecipitation [ChIP] followed by high-throughput sequencing) experiments identified Lhx1 target genes, including numerous anterior definitive endoderm markers and components of the Wnt signaling pathway. Interestingly, Lhx1-binding sites were enriched at enhancers, including the Nodal-proximal epiblast enhancer element and enhancer regions controlling Otx2 and Foxa2 expression. Moreover, in proteomic experiments, we characterized a complex comprised of Lhx1, Otx2, and Foxa2 as well as the chromatin-looping protein Ldb1. These partnerships cooperatively regulate development of the anterior mesendoderm, node, and midline cell populations responsible for establishment of the left–right body axis and head formation.

A zebrafish melanoma model reveals emergence of neural crest identity during melanoma initiation

The "cancerized field" concept posits that cancer-prone cells in a given tissue share an oncogenic mutation, but only discreet clones within the field initiate tumors. Most benign nevi carry oncogenic BRAF(V600E) mutations but rarely become melanoma. The zebrafish crestin gene is expressed embryonically in neural crest progenitors (NCPs) and specifically reexpressed in melanoma. Live imaging of transgenic zebrafish crestin reporters shows that within a cancerized field (BRAF(V600E)-mutant; p53-deficient), a single melanocyte reactivates the NCP state, revealing a fate change at melanoma initiation in this model. NCP transcription factors, including sox10, regulate crestin expression. Forced sox10 overexpression in melanocytes accelerated melanoma formation, which is consistent with activation of NCP genes and super-enhancers leading to melanoma. Our work highlights NCP state reemergence as a key event in melanoma initiation.

Brachyury, Foxa2 and the cis-Regulatory Origins of the Notochord

A main challenge of modern biology is to understand how specific constellations of genes are activated to differentiate cells and give rise to distinct tissues. This study focuses on elucidating how gene expression is initiated in the notochord, an axial structure that provides support and patterning signals to embryos of humans and all other chordates. Although numerous notochord genes have been identified, the regulatory DNAs that orchestrate development and propel evolution of this structure by eliciting notochord gene expression remain mostly uncharted, and the information on their configuration and recurrence is still quite fragmentary. Here we used the simple chordate Ciona for a systematic analysis of notochord cis-regulatory modules (CRMs), and investigated their composition, architectural constraints, predictive ability and evolutionary conservation. We found that most Ciona notochord CRMs relied upon variable combinations of binding sites for the transcription factors Brachyury and/or Foxa2, which can act either synergistically or independently from one another. Notably, one of these CRMs contains a Brachyury binding site juxtaposed to an (AC) microsatellite, an unusual arrangement also found in Brachyury-bound regulatory regions in mouse. In contrast, different subsets of CRMs relied upon binding sites for transcription factors of widely diverse families. Surprisingly, we found that neither intra-genomic nor interspecific conservation of binding sites were reliably predictive hallmarks of notochord CRMs. We propose that rather than obeying a rigid sequence-based cis-regulatory code, most notochord CRMs are rather unique. Yet, this study uncovered essential elements recurrently used by divergent chordates as basic building blocks for notochord CRMs.

Hedgehog/Patched-associated rhabdomyosarcoma formation from delta1-expressing mesodermal cells

Rhabdomyosarcoma (RMS) is the most common pediatric soft tissue sarcoma. In children, the 2 major RMS subtypes are alveolar and embryonal RMS. Aberrant Hedgehog/Patched1 (Hh/Ptch) signaling is a hallmark of embryonal RMS. We demonstrate that mice carrying a Ptch mutation in mesodermal Delta1-expressing cells develop embryonal-like RMS at a similar rate as mice harboring a Ptch mutation in the germline or the brachury-expressing mesoderm. The tumor incidence decreases dramatically when Ptch is mutated in Myf5- or Pax3-expressing cells. No RMS develop from Myogenin/Mef2c-expressing cells. This suggests that Hh/Ptch-associated RMS are derived from Delta1-positive, Myf5-negative, Myogenin-negative and Pax3-negative mesodermal progenitors that can undergo myogenic differentiation but lack stable lineage commitment. Additional preliminary genetic data and data on mesodermal progenitors further imply an interplay of Hh/Ptch and Delta/Notch signaling activity during RMS initiation. In contrast, Wnt signals supposedly suppress RMS formation because RMS multiplicity decreases after inactivation of the Wnt-inhibitor Wif1. Finally, our results strongly suggest that the tumor-initiating event determines the lineage of RMS origin.

FOXA2 attenuates the epithelial to mesenchymal transition by regulating the transcription of E-cadherin and ZEB2 in human breast cancer

The Forkhead Box A2 (FOXA2) transcription factor is required for embryonic development and for normal functions of multiple adult tissues, in which the maintained expression of FOXA2 is usually related to preventing the progression of malignant transformation. In this study, we found that FOXA2 prevented the epithelial to mesenchymal transition (EMT) in human breast cancer. We observed a strong correlation between the expression levels of FOXA2 and the epithelial phenotype. Knockdown of FOXA2 promoted the mesenchymal phenotype, whereas stable overexpression of FOXA2 attenuated EMT in breast cancer cells. FOXA2 was found to endogenously bind to and stimulate the promoter of E-cadherin that is crucial for epithelial phenotype of the tumor cells. Meanwhile, FOXA2 prevented EMT of breast cancer cells by repressing the expression of EMT-related transcription factor ZEB2 through recruiting a transcriptional corepressor TLE3 to the ZEB2 promoter. The stable overexpression of FOXA2 abolished metastasis of breast cancer cells in vivo. This study confirmed that FOXA2 inhibited EMT in breast cancer cells by regulating the transcription of EMT-related genes such as E-cadherin and ZEB2.

Proteomic signature of the murine intervertebral disc

Low back pain is the most common musculoskeletal problem and the single most common cause of disability, often attributed to degeneration of the intervertebral disc. Lack of effective treatment is directly related to our limited understanding of the pathways responsible for maintaining disc health. While transcriptional analysis has permitted initial insights into the biology of the intervertebral disc, complete proteomic characterization is required. We therefore employed liquid chromatography electrospray ionization tandem mass spectrometry (LC-ESI-MS/MS) protein/peptide separation and mass spectrometric analyses to characterize the protein content of intervertebral discs from skeletally mature wild-type mice. A total of 1360 proteins were identified and categorized using PANTHER. Identified proteins were primarily intracellular/plasma membrane (35%), organelle (30%), macromolecular complex (10%), extracellular region (9%). Molecular function categorization resulted in three distinct categories: catalytic activity (33%), binding (molecule interactions) (29%), and structural activity (13%). To validate our list, we confirmed the presence of 14 of 20 previously identified IVD-associated markers, including matrix proteins, transcriptional regulators, and secreted proteins. Immunohistochemical analysis confirmed distinct localization patterns of select protein with the intervertebral disc. Characterization of the protein composition of healthy intervertebral disc tissue is an important first step in identifying cellular processes and pathways disrupted during aging or disease progression.

Investigating microenvironmental regulation of human chordoma cell behaviour

The tumour microenvironment is complex and composed of many different constituents, including matricellular proteins such as connective tissue growth factor (CCN2), and is characterized by gradients in oxygen levels. In various cancers, hypoxia and CCN2 promote stem and progenitor cell properties, and regulate the proliferation, migration and phenotype of cancer cells. Our study was aimed at investigating the effects of hypoxia and CCN2 on chordoma cells, using the human U-CH1 cell line. We demonstrate that under basal conditions, U-CH1 cells express multiple CCN family members including CCN1, CCN2, CCN3 and CCN5. Culture of U-CH1 cells in either hypoxia or in the presence of recombinant CCN2 peptide promoted progenitor cell-like characteristics specific to the notochordal tissue of origin. Specifically, hypoxia induced the most robust increase in progenitor-like characteristics in U-CH1 cells, including increased expression of the notochord-associated markers T, CD24, FOXA1, ACAN and CA12, increased cell growth and tumour-sphere formation, and a decrease in the percentage of vacuolated cells present in the heterogeneous population. Interestingly, the effects of recombinant CCN2 peptide on U-CH1 cells were more pronounced under normoxia than hypoxia, promoting increased expression of CCN1, CCN2, CCN3 and CCN5, the notochord-associated markers SOX5, SOX6, T, CD24, and FOXA1 as well as increased tumour-sphere formation. Overall, this study highlights the importance of multiple factors within the tumour microenvironment and how hypoxia and CCN2 may regulate human chordoma cell behaviour.

In vivo genome-wide analysis of multiple tissues identifies gene regulatory networks, novel functions and downstream regulatory genes for Bapx1 and its co-regulation with Sox9 in the mammalian vertebral column

Background

Vertebrate organogenesis is a highly complex process involving sequential cascades of transcription factor activation or repression. Interestingly a single developmental control gene can occasionally be essential for the morphogenesis and differentiation of tissues and organs arising from vastly disparate embryological lineages.

Results

Here we elucidated the role of the mammalian homeobox gene Bapx1 during the embryogenesis of five distinct organs at E12.5 - vertebral column, spleen, gut, forelimb and hindlimb - using expression profiling of sorted wildtype and mutant cells combined with genome wide binding site analysis. Furthermore we analyzed the development of the vertebral column at the molecular level by combining transcriptional profiling and genome wide binding data for Bapx1 with similarly generated data sets for Sox9 to assemble a detailed gene regulatory network revealing genes previously not reported to be controlled by either of these two transcription factors.

Conclusions

The gene regulatory network appears to control cell fate decisions and morphogenesis in the vertebral column along with the prevention of premature chondrocyte differentiation thus providing a detailed molecular view of vertebral column development.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1072) contains supplementary material, which is available to authorized users.

Development of the annelid axochord: insights into notochord evolution

The origin of chordates has been debated for more than a century, with one key issue being the emergence of the notochord. In vertebrates, the notochord develops by convergence and extension of the chordamesoderm, a population of midline cells of unique molecular identity. We identify a population of mesodermal cells in a developing invertebrate, the marine annelid Platynereis dumerilii, that converges and extends toward the midline and expresses a notochord-specific combination of genes. These cells differentiate into a longitudinal muscle, the axochord, that is positioned between central nervous system and axial blood vessel and secretes a strong collagenous extracellular matrix. Ancestral state reconstruction suggests that contractile mesodermal midline cells existed in bilaterian ancestors. We propose that these cells, via vacuolization and stiffening, gave rise to the chordate notochord.

Gene transcription in the zebrafish embryo: regulators and networks

The precise spatial and temporal control of gene expression is a key process in the development, maintenance and regeneration of the vertebrate body. A substantial proportion of vertebrate genomes encode genes that control the transcription of the genetic information into mRNA. The zebrafish is particularly well suited to investigate gene regulatory networks underlying the control of gene expression during development due to the external development of its transparent embryos and the increasingly sophisticated tools for genetic manipulation available for this model system. We review here recent data on the analysis of cis-regulatory modules, transcriptional regulators and their integration into gene regulatory networks in the zebrafish, using the developing spinal cord as example.

Foxa1 and Foxa2 are required for formation of the intervertebral discs

The intervertebral disc (IVD) is composed of 3 main structures, the collagenous annulus fibrosus (AF), which surrounds the gel-like nucleus pulposus (NP), and hyaline cartilage endplates, which are attached to the vertebral bodies. An IVD is located between each vertebral body. Degeneration of the IVD is thought to be a major cause of back pain, a potentially chronic condition for which there exist few effective treatments. The NP forms from the embryonic notochord. Foxa1 and Foxa2, transcription factors in the forkhead box family, are expressed early during notochord development. However, embryonic lethality and the absence of the notochord in Foxa2 null mice have precluded the study of potential roles these genes may play during IVD formation. Using a conditional Foxa2 allele in conjunction with a tamoxifen-inducible Cre allele (ShhcreER^T2), we removed Foxa2 from the notochord of E7.5 mice null for Foxa1. Foxa1^−/−;Foxa2^c/c;ShhcreER^T2 double mutant animals had a severely deformed nucleus pulposus, an increase in cell death in the tail, decreased hedgehog signaling, defects in the notochord sheath, and aberrant dorsal-ventral patterning of the neural tube. Embryos lacking only Foxa1 or Foxa2 from the notochord were indistinguishable from control animals, demonstrating a functional redundancy for these genes in IVD formation. In addition, we provide in vivo genetic evidence that Foxa genes are required for activation of Shh in the notochord.

A novel mammal-specific three partite enhancer element regulates node and notochord-specific Noto expression

The vertebrate organizer and notochord have conserved, essential functions for embryonic development and patterning. The restricted expression of developmental regulators in these tissues is directed by specific cis-regulatory modules (CRMs) whose sequence conservation varies considerably. Some CRMs have been conserved throughout vertebrates and likely represent ancestral regulatory networks, while others have diverged beyond recognition but still function over a wide evolutionary range. Here we identify and characterize a mammalian-specific CRM required for node and notochord specific (NNC) expression of NOTO, a transcription factor essential for node morphogenesis, nodal cilia movement and establishment of laterality in mouse. A 523 bp enhancer region (NOCE) upstream the Noto promoter was necessary and sufficient for NNC expression from the endogenous Noto locus. Three subregions in NOCE together mediated full activity in vivo. Binding sites for known transcription factors in NOCE were functional in vitro but dispensable for NOCE activity in vivo. A FOXA2 site in combination with a novel motif was necessary for NOCE activity in vivo. Strikingly, syntenic regions in non-mammalian vertebrates showed no recognizable sequence similarities. In contrast to its activity in mouse NOCE did not drive NNC expression in transgenic fish. NOCE represents a novel, mammal-specific CRM required for the highly restricted Noto expression in the node and nascent notochord and thus regulates normal node development and function.