Genome-wide identification of binding sites and gene targets of Alx1, a pivotal regulator of echinoderm skeletogenesis

Integrative multi-omics increase resolution of the sea urchin posterior gut gene regulatory network at single-cell level

Drafting gene regulatory networks (GRNs) requires embryological knowledge pertaining to the cell type families, information on the regulatory genes, causal data from gene knockdown experiments and validations of the identified interactions by cis-regulatory analysis. We use multi-omics involving next-generation sequencing to obtain the necessary information for drafting the Strongylocentrotus purpuratus (Sp) posterior gut GRN. Here, we present an update to the GRN using: (1) a single-cell RNA-sequencing-derived cell atlas highlighting the 2 day-post-fertilization (dpf) sea urchin gastrula cell type families, as well as the genes expressed at the single-cell level; (2) a set of putative cis-regulatory modules and transcription factor-binding sites obtained from chromatin accessibility ATAC-seq data; and (3) interactions directionality obtained from differential bulk RNA sequencing following knockdown of the transcription factor Sp-Pdx1, a key regulator of gut patterning in sea urchins. Combining these datasets, we draft the GRN for the hindgut Sp-Pdx1-positive cells in the 2 dpf gastrula embryo. Overall, our data suggest the complex connectivity of the posterior gut GRN and increase the resolution of gene regulatory cascades operating within it.

Molecular compartmentalization in a syncytium: restricted mobility of proteins within the sea urchin skeletogenic mesenchyme

Multinucleated cells, or syncytia, are found in diverse taxa. Their biological function is often associated with the compartmentalization of biochemical or cellular activities within the syncytium. How such compartments are generated and maintained is poorly understood. The sea urchin embryonic skeleton is secreted by a syncytium, and local patterns of skeletal growth are associated with distinct sub-domains of gene expression within the syncytium. For such molecular compartments to be maintained and to control local patterns of skeletal growth: (1) the mobility of TFs must be restricted to produce stable differences in the transcriptional states of nuclei within the syncytium; and (2) the mobility of biomineralization proteins must also be restricted to produce regional differences in skeletal growth. To test these predictions, we expressed fluorescently tagged forms of transcription factors and biomineralization proteins in sub-domains of the skeletogenic syncytium. We found that both classes of proteins have restricted mobility within the syncytium and identified motifs that limit their mobility. Our findings have general implications for understanding the functional and molecular compartmentalization of syncytia.

RAS-independent ERK activation by constitutively active KSR3 in non-chordate metazoa

During early development of the sea urchin embryo, activation of ERK signalling in mesodermal precursors is not triggered by extracellular RTK ligands but by a cell-autonomous, RAS-independent mechanism that was not understood. We discovered that in these cells, ERK signalling is activated through the transcriptional activation of a gene encoding a protein related to Kinase Suppressor of Ras, that we named KSR3. KSR3 belongs to a family of catalytically inactive allosteric activators of RAF. Phylogenetic analysis revealed that genes encoding kinase defective KSR3 proteins are present in most non-chordate metazoa but have been lost in flies and nematodes. We show that the structure of KSR3 factors resembles that of several oncogenic human RAF mutants and that KSR3 from echinoderms, cnidarians and hemichordates activate ERK signalling independently of RAS when overexpressed in cultured cells. Finally, we used the sequence of KSR3 factors to identify activating mutations of human B-RAF. These findings reveal key functions for this family of factors as activators of RAF in RAS-independent ERK signalling in invertebrates. They have implications on the evolution of the ERK signalling pathway and suggest a mechanism for its co-option in the course of evolution.

An optimized Tet-On system for conditional control of gene expression in sea urchins

Sea urchins and other echinoderms are important experimental models for studying developmental processes. The lack of approaches for conditional gene perturbation, however, has made it challenging to investigate the late developmental functions of genes that have essential roles during early embryogenesis and genes that have diverse functions in multiple tissues. The doxycycline-controlled Tet-On system is a widely used molecular tool for temporally and spatially regulated transgene expression. Here, we optimized the Tet-On system to conditionally induce gene expression in sea urchin embryos. Using this approach, we explored the roles the MAPK signaling plays in skeletogenesis by expressing genes that perturb the pathway specifically in primary mesenchyme cells during later stages of development. We demonstrated the wide utility of the Tet-On system by applying it to a second sea urchin species and in cell types other than the primary mesenchyme cells. Our work provides a robust and flexible platform for the spatiotemporal regulation of gene expression in sea urchins, which will considerably enhance the utility of this prominent model system.

Acorn worm ossicle ultrastructure and composition and the origin of the echinoderm skeleton

Here, we describe the shape and mineral composition of ossicles from eight acorn worm species, bringing the total known biomineralizing enteropneusts to 10 and confirming that ossicles are widespread in Enteropneusta. Three general forms were identified including a globular form that occurs in all three major enteropneust families. The biomineral compositions included all three polymorphs of calcium carbonate; calcite, aragonite and vaterite, and low to high magnesium concentrations. Calcite was the most common and characteristic of echinoderm ossicles. Based on these findings we hypothesize that an enteropneust-like ancestor to the Ambulacraria had ectodermal ossicles, formed in an extracellular occluded space bordered by a sheath of sclerocyte cells. The ossicles were microscopic, monotypic globular shaped, calcite ossicles with low to high Mg content and MSP130 proteins. The ossicles lacked intercalation with other ossicles. The function of acorn worm ossicles is unknown, but the position of ossicles in the trunk epithelia and near to the surface suggests predator deterrence, to provide grip on the walls of a burrow or tube, as storage of metabolic waste, or to regulate blood pH, rather than as an endoskeleton function seen in fossil and crown group Echinodermata.

Acorn worm ossicle ultrastructure and composition and the origin of the echinoderm skeleton

Here, we describe the shape and mineral composition of ossicles from eight acorn worm species, bringing the total known biomineralizing enteropneusts to 10 and confirming that ossicles are widespread in Enteropneusta. Three general forms were identified including a globular form that occurs in all three major enteropneust families. The biomineral compositions included all three polymorphs of calcium carbonate; calcite, aragonite and vaterite, and low to high magnesium concentrations. Calcite was the most common and characteristic of echinoderm ossicles. Based on these findings we hypothesize that an enteropneust-like ancestor to the Ambulacraria had ectodermal ossicles, formed in an extracellular occluded space bordered by a sheath of sclerocyte cells. The ossicles were microscopic, monotypic globular shaped, calcite ossicles with low to high Mg content and MSP130 proteins. The ossicles lacked intercalation with other ossicles. The function of acorn worm ossicles is unknown, but the position of ossicles in the trunk epithelia and near to the surface suggests predator deterrence, to provide grip on the walls of a burrow or tube, as storage of metabolic waste, or to regulate blood pH, rather than as an endoskeleton function seen in fossil and crown group Echinodermata.

Architecture and evolution of the cis-regulatory system of the echinoderm kirrelL gene

The gene regulatory network (GRN) that underlies echinoderm skeletogenesis is a prominent model of GRN architecture and evolution. KirrelL is an essential downstream effector gene in this network and encodes an Ig-superfamily protein required for the fusion of skeletogenic cells and the formation of the skeleton. In this study, we dissected the transcriptional control region of the kirrelL gene of the purple sea urchin, Strongylocentrotus purpuratus. Using plasmid- and bacterial artificial chromosome-based transgenic reporter assays, we identified key cis-regulatory elements (CREs) and transcription factor inputs that regulate Sp-kirrelL, including direct, positive inputs from two key transcription factors in the skeletogenic GRN, Alx1 and Ets1. We next identified kirrelL cis-regulatory regions from seven other echinoderm species that together represent all classes within the phylum. By introducing these heterologous regulatory regions into developing sea urchin embryos we provide evidence of their remarkable conservation across ~500 million years of evolution. We dissected in detail the kirrelL regulatory region of the sea star, Patiria miniata, and demonstrated that it also receives direct inputs from Alx1 and Ets1. Our findings identify kirrelL as a component of the ancestral echinoderm skeletogenic GRN. They support the view that GRN subcircuits, including specific transcription factor–CRE interactions, can remain stable over vast periods of evolutionary history. Lastly, our analysis of kirrelL establishes direct linkages between a developmental GRN and an effector gene that controls a key morphogenetic cell behavior, cell–cell fusion, providing a paradigm for extending the explanatory power of GRNs.

Lessons from a transcription factor: Alx1 provides insights into gene regulatory networks, cellular reprogramming, and cell type evolution

The skeleton-forming cells of sea urchins and other echinoderms have been studied by developmental biologists as models of cell specification and morphogenesis for many decades. The gene regulatory network (GRN) deployed in the embryonic skeletogenic cells of euechinoid sea urchins is one of the best understood in any developing animal. Recent comparative studies have leveraged the information contained in this GRN, bringing renewed attention to the diverse patterns of skeletogenesis within the phylum and the evolutionary basis for this diversity. The homeodomain-containing transcription factor, Alx1, was originally shown to be a core component of the skeletogenic GRN of the sea urchin embryo. Alx1 has since been found to be key regulator of skeletal cell identity throughout the phylum. As such, Alx1 is currently serving as a lens through which multiple developmental processes are being investigated. These include not only GRN organization and evolution, but also cell reprogramming, cell type evolution, and the gene regulatory control of morphogenesis. This review summarizes our current state of knowledge concerning Alx1 and highlights the insights it is yielding into these important developmental and evolutionary processes.

Developmental gene regulatory network connections predicted by machine learning from gene expression data alone

Gene regulatory network (GRN) inference can now take advantage of powerful machine learning algorithms to complement traditional experimental methods in building gene networks. However, the dynamical nature of embryonic development-representing the time-dependent interactions between thousands of transcription factors, signaling molecules, and effector genes-is one of the most challenging arenas for GRN prediction. In this work, we show that successful GRN predictions for a developmental network from gene expression data alone can be obtained with the Priors Enriched Absent Knowledge (PEAK) network inference algorithm. PEAK is a noise-robust method that models gene expression dynamics via ordinary differential equations and selects the best network based on information-theoretic criteria coupled with the machine learning algorithm Elastic Net. We test our GRN prediction methodology using two gene expression datasets for the purple sea urchin, Stronglyocentrotus purpuratus, and cross-check our results against existing GRN models that have been constructed and validated by over 30 years of experimental results. Our results find a remarkably high degree of sensitivity in identifying known gene interactions in the network (maximum 81.58%). We also generate novel predictions for interactions that have not yet been described, which provide a resource for researchers to use to further complete the sea urchin GRN. Published ChIPseq data and spatial co-expression analysis further support a subset of the top novel predictions. We conclude that GRN predictions that match known gene interactions can be produced using gene expression data alone from developmental time series experiments.

Single-cell RNA sequencing of the Strongylocentrotus purpuratus larva reveals the blueprint of major cell types and nervous system of a non-chordate deuterostome

Identifying the molecular fingerprint of organismal cell types is key for understanding their function and evolution. Here, we use single-cell RNA sequencing (scRNA-seq) to survey the cell types of the sea urchin early pluteus larva, representing an important developmental transition from non-feeding to feeding larva. We identify 21 distinct cell clusters, representing cells of the digestive, skeletal, immune, and nervous systems. Further subclustering of these reveal a highly detailed portrait of cell diversity across the larva, including the identification of neuronal cell types. We then validate important gene regulatory networks driving sea urchin development and reveal new domains of activity within the larval body. Focusing on neurons that co-express Pdx-1 and Brn1/2/4, we identify an unprecedented number of genes shared by this population of neurons in sea urchin and vertebrate endocrine pancreatic cells. Using differential expression results from Pdx-1 knockdown experiments, we show that Pdx1 is necessary for the acquisition of the neuronal identity of these cells. We hypothesize that a network similar to the one orchestrated by Pdx1 in the sea urchin neurons was active in an ancestral cell type and then inherited by neuronal and pancreatic developmental lineages in sea urchins and vertebrates.

Global patterns of enhancer activity during sea urchin embryogenesis assessed by eRNA profiling

We used capped analysis of gene expression with sequencing (CAGE-seq) to profile eRNA expression and enhancer activity during embryogenesis of a model echinoderm: the sea urchin, Strongylocentrotus purpuratus. We identified more than 18,000 enhancers that were active in mature oocytes and developing embryos and documented a burst of enhancer activation during cleavage and early blastula stages. We found that a large fraction (73.8%) of all enhancers active during the first 48 h of embryogenesis were hyperaccessible no later than the 128-cell stage and possibly even earlier. Most enhancers were located near gene bodies, and temporal patterns of eRNA expression tended to parallel those of nearby genes. Furthermore, enhancers near lineage-specific genes contained signatures of inputs from developmental gene regulatory networks deployed in those lineages. A large fraction (60%) of sea urchin enhancers previously shown to be active in transgenic reporter assays was associated with eRNA expression. Moreover, a large fraction (50%) of a representative subset of enhancers identified by eRNA profiling drove tissue-specific gene expression in isolation when tested by reporter assays. Our findings provide an atlas of developmental enhancers in a model sea urchin and support the utility of eRNA profiling as a tool for enhancer discovery and regulatory biology. The data generated in this study are available at Echinobase, the public database of information related to echinoderm genomics.

Analysis of the DNA-binding properties of Alx1, an evolutionarily conserved regulator of skeletogenesis in echinoderms

Alx1, a homeodomain-containing transcription factor, is a highly conserved regulator of skeletogenesis in echinoderms. In sea urchins, Alx1 plays a central role in the differentiation of embryonic primary mesenchyme cells (PMCs) and positively regulates the transcription of most biomineralization genes expressed by these cells. The alx1 gene arose via duplication and acquired a skeletogenic function distinct from its paralog (alx4) through the exonization of a 41–amino acid motif (the D2 domain). Alx1 and Alx4 contain glutamine-50 paired-type homeodomains, which interact preferentially with palindromic binding sites in vitro. Chromatin immunoprecipitation sequencing (ChIP-seq) studies have shown, however, that Alx1 binds both to palindromic and half sites in vivo. To address this apparent discrepancy and explore the function of the D2 domain, we used an endogenous cis-regulatory module associated with Sp-mtmmpb, a gene that encodes a PMC-specific metalloprotease, to analyze the DNA-binding properties of Alx1. We find that Alx1 forms dimeric complexes on TAAT-containing half sites by a mechanism distinct from the well-known mechanism of dimerization on palindromic sites. We used transgenic reporter assays to analyze the functional roles of half sites in vivo and demonstrate that two sites with partially redundant functions are essential for the PMC-specific activity of the Sp-mtmmpb cis-regulatory module. Finally, we show that the D2 domain influences the DNA-binding properties of Alx1 in vitro, suggesting that the exonization of this motif may have facilitated the acquisition of new transcriptional targets and consequently a novel developmental function.

hReg-CNCC reconstructs a regulatory network in human cranial neural crest cells and annotates variants in a developmental context

Cranial Neural Crest Cells (CNCC) originate at the cephalic region from forebrain, midbrain and hindbrain, migrate into the developing craniofacial region, and subsequently differentiate into multiple cell types. The entire specification, delamination, migration, and differentiation process is highly regulated and abnormalities during this craniofacial development cause birth defects. To better understand the molecular networks underlying CNCC, we integrate paired gene expression & chromatin accessibility data and reconstruct the genome-wide human Regulatory network of CNCC (hReg-CNCC). Consensus optimization predicts high-quality regulations and reveals the architecture of upstream, core, and downstream transcription factors that are associated with functions of neural plate border, specification, and migration. hReg-CNCC allows us to annotate genetic variants of human facial GWAS and disease traits with associated cis-regulatory modules, transcription factors, and target genes. For example, we reveal the distal and combinatorial regulation of multiple SNPs to core TF ALX1 and associations to facial distances and cranial rare disease. In addition, hReg-CNCC connects the DNA sequence differences in evolution, such as ultra-conserved elements and human accelerated regions, with gene expression and phenotype. hReg-CNCC provides a valuable resource to interpret genetic variants as early as gastrulation during embryonic development. The network resources are available at https://github.com/AMSSwanglab/hReg-CNCC .

microRNA-31 regulates skeletogenesis by direct suppression of Eve and Wnt1

microRNAs (miRNAs) play a critical role in a variety of biological processes, including embryogenesis and the physiological functions of cells. Evolutionarily conserved microRNA-31 (miR-31) has been found to be involved in cancer, bone formation, and lymphatic development. We previously discovered that, in the sea urchin, miR-31 knockdown (KD) embryos have shortened dorsoventral connecting rods, mispatterned skeletogenic primary mesenchyme cells (PMCs) and shifted and expanded Vegf3 expression domain. Vegf3 itself does not contain miR-31 binding sites; however, we identified its upstream regulators Eve and Wnt1 to be directly suppressed by miR-31. Removal of miR-31's suppression of Eve and Wnt1 resulted in skeletal and PMC patterning defects, similar to miR-31 KD phenotypes. Additionally, removal of miR-31's suppression of Eve and Wnt1 results in an expansion and anterior shift in expression of Veg1 ectodermal genes, including Vegf3 in the blastulae. This indicates that miR-31 indirectly regulates Vegf3 expression through directly suppressing Eve and Wnt1. Furthermore, removing miR-31 suppression of Eve is sufficient to cause skeletogenic defects, revealing a novel regulatory role of Eve in skeletogenesis and PMC patterning. Overall, this study provides a proposed molecular mechanism of miR-31's regulation of skeletogenesis and PMC patterning through its cross-regulation of a Wnt signaling ligand and a transcription factor of the endodermal and ectodermal gene regulatory network.

Ultrastructural and molecular analysis of the origin and differentiation of cells mediating brittle star skeletal regeneration

BACKGROUND

Regeneration is the ability to re-grow body parts or tissues after trauma, and it is widespread across metazoans. Cells involved in regeneration can arise from a pool of undifferentiated proliferative cells or be recruited from pre-existing differentiated tissues. Both mechanisms have been described in different phyla; however, the cellular and molecular mechanisms employed by different animals to restore lost tissues as well as the source of cells involved in regeneration remain largely unknown. Echinoderms are a clade of deuterostome invertebrates that show striking larval and adult regenerative abilities in all extant classes. Here, we use the brittle star Amphiura filiformis to investigate the origin and differentiation of cells involved in skeletal regeneration using a combination of microscopy techniques and molecular markers.

RESULTS

Our ultrastructural analyses at different regenerative stages identify a population of morphologically undifferentiated cells which appear in close contact with the proliferating epithelium of the regenerating aboral coelomic cavity. These cells express skeletogenic marker genes, such as the transcription factor alx1 and the differentiation genes c-lectin and msp130L, and display a gradient of morphological differentiation from the aboral coelomic cavity towards the epidermis. Cells closer to the epidermis, which are in contact with developing spicules, have the morphology of mature skeletal cells (sclerocytes), and express several skeletogenic transcription factors and differentiation genes. Moreover, as regeneration progresses, sclerocytes show a different combinatorial expression of genes in various skeletal elements.

CONCLUSIONS

We hypothesize that sclerocyte precursors originate from the epithelium of the proliferating aboral coelomic cavity. As these cells migrate towards the epidermis, they differentiate and start secreting spicules. Moreover, our study shows that molecular and cellular processes involved in skeletal regeneration resemble those used during skeletal development, hinting at a possible conservation of developmental programmes during adult regeneration. Finally, we highlight that many genes involved in echinoderm skeletogenesis also play a role in vertebrate skeleton formation, suggesting a possible common origin of the deuterostome endoskeleton pathway.

Transcription Factors of the Alx Family: Evolutionarily Conserved Regulators of Deuterostome Skeletogenesis

Members of the alx gene family encode transcription factors that contain a highly conserved Paired-class, DNA-binding homeodomain, and a C-terminal OAR/Aristaless domain. Phylogenetic and comparative genomic studies have revealed complex patterns of alx gene duplications during deuterostome evolution. Remarkably, alx genes have been implicated in skeletogenesis in both echinoderms and vertebrates. In this review, we provide an overview of current knowledge concerning alx genes in deuterostomes. We highlight their evolutionarily conserved role in skeletogenesis and draw parallels and distinctions between the skeletogenic gene regulatory circuitries of diverse groups within the superphylum.

ALX1-related frontonasal dysplasia results from defective neural crest cell development and migration

A pedigree of subjects presented with frontonasal dysplasia (FND). Genome sequencing and analysis identified a p.L165F missense variant in the homeodomain of the transcription factor ALX1 which was imputed to be pathogenic. Induced pluripotent stem cells (iPSC) were derived from the subjects and differentiated to neural crest cells (NCC). NCC derived from ALX1L165F/L165F iPSC were more sensitive to apoptosis, showed an elevated expression of several neural crest progenitor state markers, and exhibited impaired migration compared to wild-type controls. NCC migration was evaluated in vivo using lineage tracing in a zebrafish model, which revealed defective migration of the anterior NCC stream that contributes to the median portion of the anterior neurocranium, phenocopying the clinical presentation. Analysis of human NCC culture media revealed a change in the level of bone morphogenic proteins (BMP), with a low level of BMP2 and a high level of BMP9. Soluble BMP2 and BMP9 antagonist treatments were able to rescue the defective migration phenotype. Taken together, these results demonstrate a mechanistic requirement of ALX1 in NCC development and migration.

The gene regulatory control of sea urchin gastrulation

The cell behaviors associated with gastrulation in sea urchins have been well described. More recently, considerable progress has been made in elucidating gene regulatory networks (GRNs) that underlie the specification of early embryonic territories in this experimental model. This review integrates information from these two avenues of work. I discuss the principal cell movements that take place during sea urchin gastrulation, with an emphasis on molecular effectors of the movements, and summarize our current understanding of the gene regulatory circuitry upstream of those effectors. A case is made that GRN biology can provide a causal explanation of gastrulation, although additional analysis is needed at several levels of biological organization in order to provide a deeper understanding of this complex morphogenetic process.

The evolution of a new cell type was associated with competition for a signaling ligand

There is presently a very limited understanding of the mechanisms that underlie the evolution of new cell types. The skeleton-forming primary mesenchyme cells (PMCs) of euechinoid sea urchins, derived from the micromeres of the 16-cell embryo, are an example of a recently evolved cell type. All adult echinoderms have a calcite-based endoskeleton, a synapomorphy of the Ambulacraria. Only euechinoids have a micromere-PMC lineage, however, which evolved through the co-option of the adult skeletogenic program into the embryo. During normal development, PMCs alone secrete the embryonic skeleton. Other mesoderm cells, known as blastocoelar cells (BCs), have the potential to produce a skeleton, but a PMC-derived signal ordinarily prevents these cells from expressing a skeletogenic fate and directs them into an alternative developmental pathway. Recently, it was shown that vascular endothelial growth factor (VEGF) signaling plays an important role in PMC differentiation and is part of a conserved program of skeletogenesis among echinoderms. Here, we report that VEGF signaling, acting through ectoderm-derived VEGF3 and its cognate receptor, VEGF receptor (VEGFR)-10-Ig, is also essential for the deployment of the skeletogenic program in BCs. This VEGF-dependent program includes the activation of aristaless-like homeobox 1 (alx1), a conserved transcriptional regulator of skeletogenic specification across echinoderms and an example of a “terminal selector” gene that controls cell identity. We show that PMCs control BC fate by sequestering VEGF3, thereby preventing activation of alx1 and the downstream skeletogenic network in BCs. Our findings provide an example of the regulation of early embryonic cell fates by direct competition for a secreted signaling ligand, a developmental mechanism that has not been widely recognized. Moreover, they reveal that a novel cell type evolved by outcompeting other embryonic cell lineages for an essential signaling ligand that regulates the expression of a gene controlling cell identity.

How do new cell types evolve? This study shows that mesoderm cells in sea urchin embryos diversified, at least in part, through a heterochronic shift in the expression of a key transcription factor, which led to competition for a signaling ligand and subsequent gene regulatory independence of the two cell types.