Signal-dependent regulation of the sea urchin skeletogenic gene regulatory network

Localization and origins of juvenile skeletogenic cells in the sea urchin Lytechinuspictus

The development of the sea urchin larval body plan is well understood from extensive studies of embryonic patterning. However, fewer studies have investigated the late larval stages during which the unique pentaradial adult body plan develops. Previous work on late larval development highlights major tissue changes leading up to metamorphosis, but the location of specific cell types during juvenile development is less understood. Here, we improve on technical limitations by applying highly sensitive hybridization chain reaction fluorescent in situ hybridization (HCR-FISH) to the fast-developing and transparent sea urchin Lytechinus pictus, with a focus on skeletogenic cells. First, we show that HCR-FISH can be used in L. pictus to precisely localize skeletogenic cells in the rudiment. In doing so, we provide a detailed staging scheme for the appearance of skeletogenic cells around the rudiment prior to and during biomineralization and show that many skeletogenic cells unassociated with larval rods localize outside of the rudiment prior to localizing inside. Second, we show that downstream biomineralization genes have similar expression patterns during larval and juvenile skeletogenesis, suggesting some conservation of skeletogenic mechanisms during development between stages. Third, we find co-expression of blastocoelar and skeletogenic cell markers around juvenile skeleton located outside of the rudiment, which is consistent with data showing that cells from the non-skeletogenic mesoderm embryonic lineage contribute to the juvenile skeletogenic cell lineage. This work sets the foundation for subsequent studies of other cell types in the late larva of L. pictus to better understand juvenile body plan development, patterning, and evolution.

Expression and Transcriptional Targets of TGFβ-RII in Paracentrotus lividus Larval Skeletogenesis

Organisms from the five kingdoms of life use minerals to harden their tissues and make teeth, shells and skeletons, in the process of biomineralization. The sea urchin larval skeleton is an excellent system to study the biological regulation of biomineralization and its evolution. The gene regulatory network (GRN) that controls sea urchin skeletogenesis is known in great details and shows similarity to the GRN that controls vertebrates' vascularization while it is quite distinct from the GRN that drives vertebrates' bone formation. Yet, transforming growth factor beta (TGF-β) signaling regulates both sea urchin and vertebrates' skeletogenesis. Here, we study the upstream regulation and identify transcriptional targets of TGF-β in the Mediterranean Sea urchin species, Paracentrotus lividus. TGF-βRII is transiently active in the skeletogenic cells downstream of vascular endothelial growth factor (VEGF) signaling, in P. lividus. Continuous perturbation of TGF-βRII activity significantly impairs skeletal elongation and the expression of key skeletogenic genes. Perturbation of TGF-βRII after skeletal initiation leads to a delay in skeletal elongation and minor changes in gene expression. TGF-β targets are distinct from its transcriptional targets during vertebrates' bone formation, suggesting that the role of TGF-β in biomineralization in these two phyla results from convergent evolution.

Echinobase: a resource to support the echinoderm research community

Echinobase (www.echinobase.org) is a model organism knowledgebase serving as a resource for the community that studies echinoderms, a phylum of marine invertebrates that includes sea urchins and sea stars. Echinoderms have been important experimental models for over 100 years and continue to make important contributions to environmental, evolutionary, and developmental studies, including research on developmental gene regulatory networks. As a centralized resource, Echinobase hosts genomes and collects functional genomic data, reagents, literature, and other information for the community. This third-generation site is based on the Xenbase knowledgebase design and utilizes gene-centric pages to minimize the time and effort required to access genomic information. Summary gene pages display gene symbols and names, functional data, links to the JBrowse genome browser, and orthology to other organisms and reagents, and tabs from the Summary gene page contain more detailed information concerning mRNAs, proteins, diseases, and protein-protein interactions. The gene pages also display 1:1 orthologs between the fully supported species Strongylocentrotus purpuratus (purple sea urchin), Lytechinus variegatus (green sea urchin), Patiria miniata (bat star), and Acanthaster planci (crown-of-thorns sea star). JBrowse tracks are available for visualization of functional genomic data from both fully supported species and the partially supported species Anneissia japonica (feather star), Asterias rubens (sugar star), and L. pictus (painted sea urchin). Echinobase serves a vital role by providing researchers with annotated genomes including orthology, functional genomic data aligned to the genomes, and curated reagents and data. The Echinoderm Anatomical Ontology provides a framework for standardizing developmental data across the phylum, and knowledgebase content is formatted to be findable, accessible, interoperable, and reusable by the research community.

Effects of seawater acidification and warming on morphometrics and biomineralization-related gene expression during embryo-larval development of a lightly-calcified echinoderm

CO2-induced ocean acidification and warming pose ecological threats to marine life, especially calcifying species such as echinoderms, who rely on biomineralization for skeleton formation. However, previous studies on echinoderm calcification amid climate change had a strong bias towards heavily calcified echinoderms, with little research on lightly calcified ones, such as sea cucumbers. Here, we analyzed the embryo-larval development and their biomineralization-related gene expression of a lightly calcified echinoderm, the sea cucumber (Apostichopus japonicus), under experimental seawater acidification (OA) and/or warming (OW). Results showed that OA (- 0.37 units) delayed development and decreased body size (8.58-56.25 % and 0.36-19.66 % decreases in stage duration and body length, respectively), whereas OW (+3.1 °C) accelerated development and increased body size (33.99-55.28 % increase in stage duration and 2.44-14.41 % enlargement in body length). OW buffered the negative effects of OA on the development timing and body size of A. japonicus. Additionally, no target genes were expressed in the blastula stage, and only two biomineralization genes (colp3α, cyp2) and five TFs (erg, tgif, foxN2/3, gata1/2/3, and tbr) were expressed throughout the embryo-larval development. Our findings suggest that the low calcification in A. japonicus larvae may be caused by biomineralization genes contraction, and low expression of those genes. Furthermore, this study indicated that seawater acidification and warming affect expression of biomineralization-related genes, and had an effect on body size and development rate during the embryo-larval stage in sea cucumbers. Our study is a first step toward a better understanding of the complexity of high pCO2 on calcification and helpful for revealing the adaptive strategy of less-calcified echinoderms amid climate change.

ROCK and the actomyosin network control biomineral growth and morphology during sea urchin skeletogenesis

Biomineralization had apparently evolved independently in different phyla, using distinct minerals, organic scaffolds, and gene regulatory networks (GRNs). However, diverse eukaryotes from unicellular organisms, through echinoderms to vertebrates, use the actomyosin network during biomineralization. Specifically, the actomyosin remodeling protein, Rho-associated coiled-coil kinase (ROCK) regulates cell differentiation and gene expression in vertebrates' biomineralizing cells, yet, little is known on ROCK's role in invertebrates' biomineralization. Here, we reveal that ROCK controls the formation, growth, and morphology of the calcite spicules in the sea urchin larva. ROCK expression is elevated in the sea urchin skeletogenic cells downstream of the Vascular Endothelial Growth Factor (VEGF) signaling. ROCK inhibition leads to skeletal loss and disrupts skeletogenic gene expression. ROCK inhibition after spicule formation reduces the spicule elongation rate and induces ectopic spicule branching. Similar skeletogenic phenotypes are observed when ROCK is inhibited in a skeletogenic cell culture, indicating that these phenotypes are due to ROCK activity specifically in the skeletogenic cells. Reduced skeletal growth and enhanced branching are also observed under direct perturbations of the actomyosin network. We propose that ROCK and the actomyosin machinery were employed independently, downstream of distinct GRNs, to regulate biomineral growth and morphology in Eukaryotes.

Composite material in the sea urchin Cidaris rugosa: ordered and disordered micrometre-scale bicontinuous geometries

The sponge-like biomineralized calcite materials found in echinoderm skeletons are of interest in terms of both structure formation and biological function. Despite their crystalline atomic structure, they exhibit curved interfaces that have been related to known triply periodic minimal surfaces. Here, we investigate the endoskeleton of the sea urchin Cidaris rugosa that has long been known to form a microstructure related to the Primitive surface. Using X-ray tomography, we find that the endoskeleton is organized as a composite material consisting of domains of bicontinuous microstructures with different structural properties. We describe, for the first time, the co-occurrence of ordered single Primitive and single Diamond structures and of a disordered structure within a single skeletal plate. We show that these structures can be distinguished by structural properties including solid volume fraction, trabeculae width and, to a lesser extent, interface area and mean curvature. In doing so, we present a robust method that extracts interface areas and curvature integrals from voxelized datasets using the Steiner polynomial for parallel body volumes. We discuss these very large-scale bicontinuous structures in the context of their function, formation and evolution.

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.

Polychrome labeling reveals skeletal triradiate and elongation dynamics and abnormalities in patterning cue-perturbed embryos

The larval skeleton of the sea urchin Lytechinus variegatus is an ideal model system for studying skeletal patterning; however, our understanding of the etiology of skeletal patterning in sea urchin larvae is limited due to the lack of approaches to live-image skeleton formation. Calcium-binding fluorochromes have been used to study the temporal dynamics of bone growth and healing. To date, only calcein green has been used in sea urchin larvae to fluorescently label the larval skeleton. Here, we optimize labeling protocols for two additional calcium-binding fluorochromes: xylenol orange and calcein blue- and demonstrate that these fluorochromes can be used individually or in nested pulse-chase experiments to understand the temporal dynamics of skeletogenesis and patterning. Using a pulse-chase approach, we show that the initiation of skeletogenesis begins around 15 h post fertilization. We also assess the timing of triradiate formation in embryos treated with a range of patterning perturbagens and demonstrate that triradiate formation is delayed and asynchronous in embryos ventralized via treatment with either nickel or chlorate. Finally, we measure the extent of fluorochrome incorporation in triple-labeled embryos to determine the elongation rate of numerous skeletal elements throughout early skeletal patterning and compare this to the rate of skeletal growth in embryos treated with axitinib to inhibit VEGFR. We find that skeletal elements elongate much more slowly in axitinib-treated embryos, and that axitinib treatment is sufficient to induce abnormal orientation of the triradiates.

The dopamine effect on sea urchin larvae depends on their age

Activation of the dopamine type-D₂ receptor in late gastrula of sea urchins is known to decrease the growth rate of post-oral arms of larvae, and, as a result, the phenotype of these larvae mimics that of larvae developing in the abundance of food. Our data indicate that the effect of dopamine on sea urchin larvae is stage-dependent. In our experiment, the early four-armed plutei (96 hours post fertilization, hpf) of Strongylocentrotus intermedius had substantially shorter post-oral arms if they developed from the larvae treated with dopamine at the early pluteus stage (48 hpf), when they had already formed the first dopaminergic neurons, as compared to the plutei from the larvae treated with dopamine at the mid to late gastrula stage (24 hpf), when they did not have any neurons yet. The pre-treatment of larvae in 6-hydroxydopamine, a neurotoxic analog of dopamine that specifically disrupts activity of dopaminergic neurons, prevented the development of the short post-oral arms phenotype in larvae. These results confirm the assumption that dopaminergic neurons play an important role in the development of the short post-oral arms phenotype in sea urchin larvae. Another finding of our study is that the dopamine treatment also affects the growth of the body rods and the overall larval body growth. Based on these observations, we suggest researchers to carefully select the developmental stage, pharmacological agents, and incubation time for experimental manipulation of sea urchin larvae phenotypes through dopaminergic nervous system.

Ethanol exposure perturbs sea urchin development and disrupts developmental timing

Ethanol is a known vertebrate teratogen that causes craniofacial defects as a component of fetal alcohol syndrome (FAS). Our results show that sea urchin embryos treated with ethanol similarly show broad skeletal patterning defects, potentially analogous to the defects associated with FAS. The sea urchin larval skeleton is a simple patterning system that involves only two cell types: the primary mesenchymal cells (PMCs) that secrete the calcium carbonate skeleton and the ectodermal cells that provide migratory, positional, and differentiation cues for the PMCs. Perturbations in RA biosynthesis and Hh signaling pathways are thought to be causal for the FAS phenotype in vertebrates. Surprisingly, our results indicate that these pathways are not functionally relevant for the teratogenic effects of ethanol in developing sea urchins. We found that developmental morphology as well as the expression of some ectodermal and PMC genes was delayed by ethanol exposure. Temporal transcriptome analysis revealed significant impacts of ethanol on signaling and metabolic gene expression, and a disruption in the timing of GRN gene expression that includes both delayed and precocious gene expression throughout the specification network. We conclude that the skeletal patterning perturbations in ethanol-treated embryos likely arise from a loss of temporal synchrony within and between the instructive and responsive tissues.

Carbonic anhydrases in development: morphological observations and gene expression profiling in sea urchin embryos exposed to acetazolamide

Carbonic anhydrases (CANs) are conserved metalloenzymes catalysing the reversible hydration of carbon dioxide into protons and bicarbonate, with important roles in cells physiology. Some CAN-coding genes were found in sea urchin genome, although only one involved in embryonic skeletogenesis was described in Paracentrotus lividus. Here, we investigated gene expression patterns of P. lividus embryos cultured in the presence of acetazolamide (AZ), a CAN inhibitor, to combine morphological defects with their molecular underpinning. CAN inhibition blocked skeletogenesis, affected the spatial/temporal expression of some biomineralization-related genes, inhibited embryos swimming. A comparative analysis on the expression of 127 genes in control and 3 h/24 h AZ-treated embryos, using NanoString technology, showed the differential expression of genes encoding for structural/regulatory proteins, with different embryonic roles: biomineralization, transcriptional regulation, signalling, development and defence response. The study of the differentially expressed genes and the signalling pathways affected, besides in silico analyses and a speculative 'interactomic model', leads to predicting the presence of various CAN isoforms, possibly involved in different physiological processes/activities in sea urchin embryo, and their potential target genes/proteins. Our findings provide new valuable molecular data for further studies in several biological fields: developmental biology (biomineralization, axes patterning), cell differentiation (neural development) and drug toxicology (AZ effects on embryos/tissues).

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.

Thyroid hormone-induced cell death in sea urchin metamorphic development

Thyroid hormones (THs) are important regulators of development, metabolism and homeostasis in metazoans. Specifically, they have been shown to regulate the metamorphic transitions of vertebrates and invertebrates alike. Indirectly developing sea urchin larvae accelerate the formation of juvenile structures in response to thyroxine (T4) treatment, while reducing their larval arm length. The mechanisms underlying larval arm reduction are unknown and we hypothesized that programmed cell death (PCD) is linked to this process. To test this hypothesis, we measured larval arm retraction in response to different THs (T4, T3, rT3, Tetrac) and assessed cell death in larvae using three different methods (TUNEL, YO-PRO-1 and caspase-3 activity) in the sea urchin Strongylocentrotus purpuratus. We also compared the extent of PCD in response to TH treatment before and after the invagination of the larval ectoderm, which marks the initiation of juvenile development in larval sea urchin species. We found that T4 treatment results in the strongest reduction of larval arms but detected a significant increase of PCD in response to T4, T3 and Tetrac in post-ingression but not pre-ingression larvae. As post-ingression larvae have initiated metamorphic development and therefore allocate resources to both larval and the juvenile structures, these results provide evidence that THs regulate larval development differentially via PCD. PCD in combination with cell proliferation likely has a key function in sea urchin development.

Single-cell RNA-sequencing analysis of early sea star development

Echinoderms represent a broad phylum with many tractable features to test evolutionary changes and constraints. Here, we present a single-cell RNA-sequencing analysis of early development in the sea star Patiria miniata, to complement the recent analysis of two sea urchin species. We identified 20 cell states across six developmental stages from 8 hpf to mid-gastrula stage, using the analysis of 25,703 cells. The clusters were assigned cell states based on known marker gene expression and by in situ RNA hybridization. We found that early (morula, 8-14 hpf) and late (blastula-to-mid-gastrula) cell states are transcriptionally distinct. Cells surrounding the blastopore undergo rapid cell state changes that include endomesoderm diversification. Of particular import to understanding germ cell specification is that we never see Nodal pathway members within Nanos/Vasa-positive cells in the region known to give rise to the primordial germ cells (PGCs). The results from this work contrast the results of PGC specification in the sea urchin, and the dataset presented here enables deeper comparative studies in tractable developmental models for testing a variety of developmental mechanisms.

Extracellular carbonic anhydrase activity promotes a carbon concentration mechanism in metazoan calcifying cells

Many calcifying organisms utilize metabolic CO₂ to generate CaCO₃ minerals to harden their shells and skeletons. Carbonic anhydrases are evolutionary ancient enzymes that have been proposed to play a key role in the calcification process, with the underlying mechanisms being little understood. Here, we used the calcifying primary mesenchyme cells (PMCs) of sea urchin larva to study the role of cytosolic (iCAs) and extracellular carbonic anhydrases (eCAs) in the cellular carbon concentration mechanism (CCM). Molecular analyses identified iCAs and eCAs in PMCs and highlight the prominent expression of a glycosylphosphatidylinositol-anchored membrane-bound CA (Cara7). Intracellular pH recordings in combination with CO₂ pulse experiments demonstrated iCA activity in PMCs. iCA activity measurements, together with pharmacological approaches, revealed an opposing contribution of iCAs and eCAs on the CCM. H⁺-selective electrodes were used to demonstrate eCA-catalyzed CO₂ hydration rates at the cell surface. Knockdown of Cara7 reduced extracellular CO₂ hydration rates accompanied by impaired formation of specific skeletal segments. Finally, reduced pH_i regulatory capacities during inhibition and knockdown of Cara7 underscore a role of this eCA in cellular HCO₃^- uptake. This work reveals the function of CAs in the cellular CCM of a marine calcifying animal. Extracellular hydration of metabolic CO₂ by Cara7 coupled to HCO₃^- uptake mechanisms mitigates the loss of carbon and reduces the cellular proton load during the mineralization process. The findings of this work provide insights into the cellular mechanisms of an ancient biological process that is capable of utilizing CO₂ to generate a versatile construction material.

Evolutionary analyses of genes in Echinodermata offer insights towards the origin of metazoan phyla

Despite recent studies discussing the evolutionary impacts of gene duplications and losses among metazoans, the genomic basis for the evolution of phyla remains enigmatic. Here, we employ phylogenomic approaches to search for orthologous genes without known functions among echinoderms, and subsequently use them to guide the identification of their homologs across other metazoans. Our final set of 14 genes was obtained via a suite of homology prediction tools, gene expression data, gene ontology, and generating the Strongylocentrotus purpuratus phylome. The gene set was subjected to selection pressure analyses, which indicated that they are highly conserved and under negative selection. Their presence across broad taxonomic depths suggests that genes required to form a phylum are ancestral to that phylum. Therefore, rather than de novo gene genesis, we posit that evolutionary forces such as selection on existing genomic elements over large timescales may drive divergence and contribute to the emergence of phyla.

Distinct regulatory states control the elongation of individual skeletal rods in the sea urchin embryo

Background

Understanding how gene regulatory networks (GRNs) control developmental progression is a key to the mechanistic understanding of morphogenesis. The sea urchin larval skeletogenesis provides an excellent platform to tackle this question. In the early stages of sea urchin skeletogenesis, skeletogenic genes are uniformly expressed in the skeletogenic lineage. Yet, during skeletal elongation, skeletogenic genes are expressed in distinct spatial sub‐domains. The regulation of differential gene expression during late skeletogenesis is not well understood.

Results

Here we reveal the dynamic expression of the skeletogenic regulatory genes that define a specific regulatory state for each pair of skeletal rods, in the sea urchin Paracentrotus lividus. The vascular endothelial growth factor (VEGF) signaling, essential for skeleton formation, specifically controls the migration of cells that form the postoral and distal anterolateral skeletogenic rods. VEGF signaling also controls the expression of regulatory genes in cells at the tips of the postoral rods, including the transcription factors Pitx1 and MyoD1. Pitx1 activity is required for normal skeletal elongation and for the expression of some of VEGF target genes.

Conclusions

Our study illuminates the fine‐tuning of the regulatory system during the transition from early to late skeletogenesis that gives rise to rod‐specific regulatory states.

The skeletogenic transcription factors form specific regulatory states in various skeletogenic sub‐populations.

Late VEGF signaling controls the regulatory states at the tips of the post‐oral and anterolateral skeletal rods.

VEGF signaling controls the expression of the transcription factors, MyoD1 and Pitx1.

Pitx1 activity is required for normal skeletal elongation and for the expression of some of VEGF target genes.

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.

The Evolution of Biomineralization through the Co-Option of Organic Scaffold Forming Networks

Biomineralization is the process in which organisms use minerals to generate hard structures like teeth, skeletons and shells. Biomineralization is proposed to have evolved independently in different phyla through the co-option of pre-existing developmental programs. Comparing the gene regulatory networks (GRNs) that drive biomineralization in different species could illuminate the molecular evolution of biomineralization. Skeletogenesis in the sea urchin embryo was extensively studied and the underlying GRN shows high conservation within echinoderms, larval and adult skeletogenesis. The organic scaffold in which the calcite skeletal elements form in echinoderms is a tubular compartment generated by the syncytial skeletogenic cells. This is strictly different than the organic cartilaginous scaffold that vertebrates mineralize with hydroxyapatite to make their bones. Here I compare the GRNs that drive biomineralization and tubulogenesis in echinoderms and in vertebrates. The GRN that drives skeletogenesis in the sea urchin embryo shows little similarity to the GRN that drives bone formation and high resemblance to the GRN that drives vertebrates’ vascular tubulogenesis. On the other hand, vertebrates’ bone-GRNs show high similarity to the GRNs that operate in the cells that generate the cartilage-like tissues of basal chordate and invertebrates that do not produce mineralized tissue. These comparisons suggest that biomineralization in deuterostomes evolved through the phylum specific co-option of GRNs that control distinct organic scaffolds to mineralization.

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.

The biological regulation of sea urchin larval skeletogenesis - From genes to biomineralized tissue

Biomineralization is the process in which soft organic tissues use minerals to produce shells, skeletons and teeth for various functions such as protection and physical support. The ability of the cells to control the time and place of crystal nucleation as well as crystal orientation and stiffness is far beyond the state-of-the art of human technologies. Thus, understanding the biological control of biomineralization will promote our understanding of embryo development as well as provide novel approaches for material engineering. Sea urchin larval skeletogenesis offers an excellent platform for functional analyses of both the molecular control system and mineral uptake and deposition. Here we describe the current understanding of the genetic, molecular and cellular processes that underlie sea urchin larval skeletogenesis. We portray the regulatory genes that define the specification of the skeletogenic cells and drive the various morphogenetic processes that occur in the skeletogenic lineage, including: epithelial to mesenchymal transition, cell migration, spicule cavity formation and mineral deposition into the spicule cavity. We describe recent characterizations of the size, motion and mineral concentration of the calcium-bearing vesicles in the skeletogenic cells. We review the distinct specification states within the skeletogenic lineage that drive localized skeletal growth at the tips of the spicules. Finally, we discuss the surprising similarity between the regulatory network and cellular processes that drive sea urchin skeletogenesis and those that control vertebrate vascularization. Overall, we illustrate the novel insights on the biological regulation and evolution of biomineralization, gained from studies of the sea urchin larval skeletogenesis.

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.

An otopetrin family proton channel promotes cellular acid efflux critical for biomineralization in a marine calcifier

Otopetrins comprise a family of proton-selective channels that are critically important for the mineralization of otoliths and statoconia in vertebrates but whose underlying cellular mechanisms remain largely unknown. Here, we demonstrate that otopetrins are critically involved in the calcification process by providing an exit route for protons liberated by the formation of CaCO₃ Using the sea urchin larva, we examined the otopetrin ortholog otop2l, which is exclusively expressed in the calcifying primary mesenchymal cells (PMCs) that generate the calcitic larval skeleton. otop2l expression is stimulated during skeletogenesis, and knockdown of otop2l impairs spicule formation. Intracellular pH measurements demonstrated Zn²⁺-sensitive H⁺ fluxes in PMCs that regulate intracellular pH in a Na⁺/HCO₃ ^--independent manner, while Otop2l knockdown reduced membrane proton permeability. Furthermore, Otop2l displays unique features, including strong activation by high extracellular pH (>8.0) and check-valve-like outwardly rectifying H⁺ flux properties, making it into a cellular proton extrusion machine adapted to oceanic living conditions. Our results provide evidence that otopetrin family proton channels are a central component of the cellular pH regulatory machinery in biomineralizing cells. Their ubiquitous occurrence in calcifying systems across the animal kingdom suggest a conserved physiological function by mediating pH at the site of mineralization. This important role of otopetrin family proton channels has strong implications for our view on the cellular mechanisms of biomineralization and their response to changes in oceanic pH.

Integration of 1:1 orthology maps and updated datasets into Echinobase

Echinobase (https://echinobase.org) is a central online platform that generates, manages and hosts genomic data relevant to echinoderm research. While the resource primarily serves the echinoderm research community, the recent release of an excellent quality genome for the frequently studied purple sea urchin (Strongylocentrotus purpuratus genome, v5.0) has provided an opportunity to adapt to the needs of a broader research community across other model systems. To this end, establishing pipelines to identify orthologous genes between echinoderms and other species has become a priority in many contexts including nomenclature, linking to data in other model organisms, and in internal functionality where data gathered in one hosted species can be associated with genes in other hosted echinoderms. This paper describes the orthology pipelines currently employed by Echinobase and how orthology data are processed to yield 1:1 ortholog mappings between a variety of echinoderms and other model taxa. We also describe functions of interest that have recently been included on the resource, including an updated developmental time course for S.purpuratus, and additional tracks for genome browsing. These data enhancements will increase the accessibility of the resource to non-echinoderm researchers and simultaneously expand the data quality and quantity available to core Echinobase users. Database URL: https://echinobase.org.

The tolerance to hypoxia is defined by a time-sensitive response of the gene regulatory network in sea urchin embryos

Deoxygenation, the reduction of oxygen level in the oceans induced by global warming and anthropogenic disturbances, is a major threat to marine life. This change in oxygen level could be especially harmful to marine embryos that use endogenous hypoxia and redox gradients as morphogens during normal development. Here, we show that the tolerance to hypoxic conditions changes between different developmental stages of the sea urchin embryo, possibly due to the structure of the gene regulatory networks (GRNs). We demonstrate that during normal development, the bone morphogenetic protein (BMP) pathway restricts the activity of the vascular endothelial growth factor (VEGF) pathway to two lateral domains and this restriction controls proper skeletal patterning. Hypoxia applied during early development strongly perturbs the activity of Nodal and BMP pathways that affect the VEGF pathway, dorsal-ventral (DV) and skeletogenic patterning. These pathways are largely unaffected by hypoxia applied after DV-axis formation. We propose that the use of redox and hypoxia as morphogens makes the sea urchin embryo highly sensitive to environmental hypoxia during early development, but the GRN structure provides higher tolerance to hypoxia at later stages.

Summary: The use of hypoxia and redox gradients as morphogens makes sea urchin early development sensitive to environmental hypoxia. This sensitivity decreases later, possibly due to the gene regulatory network structure.

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.

Post-metamorphic skeletal growth in the sea urchin Paracentrotus lividus and implications for body plan evolution

Background

Understanding the molecular and cellular processes that underpin animal development are crucial for understanding the diversity of body plans found on the planet today. Because of their abundance in the fossil record, and tractability as a model system in the lab, skeletons provide an ideal experimental model to understand the origins of animal diversity. We herein use molecular and cellular markers to understand the growth and development of the juvenile sea urchin (echinoid) skeleton.

Results

We developed a detailed staging scheme based off of the first ~ 4 weeks of post-metamorphic life of the regular echinoid Paracentrotus lividus. We paired this scheme with immunohistochemical staining for neuronal, muscular, and skeletal tissues, and fluorescent assays of skeletal growth and cell proliferation to understand the molecular and cellular mechanisms underlying skeletal growth and development of the sea urchin body plan.

Conclusions

Our experiments highlight the role of skeletogenic proteins in accretionary skeletal growth and cell proliferation in the addition of new metameric tissues. Furthermore, this work provides a framework for understanding the developmental evolution of sea urchin body plans on macroevolutionary timescales.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13227-021-00174-1.

Calcium-vesicles perform active diffusion in the sea urchin embryo during larval biomineralization

Biomineralization is the process by which organisms use minerals to harden their tissues and provide them with physical support. Biomineralizing cells concentrate the mineral in vesicles that they secret into a dedicated compartment where crystallization occurs. The dynamics of vesicle motion and the molecular mechanisms that control it, are not well understood. Sea urchin larval skeletogenesis provides an excellent platform for investigating the kinetics of mineral-bearing vesicles. Here we used lattice light-sheet microscopy to study the three-dimensional (3D) dynamics of calcium-bearing vesicles in the cells of normal sea urchin embryos and of embryos where skeletogenesis is blocked through the inhibition of Vascular Endothelial Growth Factor Receptor (VEGFR). We developed computational tools for displaying 3D-volumetric movies and for automatically quantifying vesicle dynamics. Our findings imply that calcium vesicles perform an active diffusion motion in both, calcifying (skeletogenic) and non-calcifying (ectodermal) cells of the embryo. The diffusion coefficient and vesicle speed are larger in the mesenchymal skeletogenic cells compared to the epithelial ectodermal cells. These differences are possibly due to the distinct mechanical properties of the two tissues, demonstrated by the enhanced f-actin accumulation and myosinII activity in the ectodermal cells compared to the skeletogenic cells. Vesicle motion is not directed toward the biomineralization compartment, but the vesicles slow down when they approach it, and probably bind for mineral deposition. VEGFR inhibition leads to an increase of vesicle volume but hardly changes vesicle kinetics and doesn’t affect f-actin accumulation and myosinII activity. Thus, calcium vesicles perform an active diffusion motion in the cells of the sea urchin embryo, with diffusion length and speed that inversely correlate with the strength of the actomyosin network. Overall, our studies provide an unprecedented view of calcium vesicle 3D-dynamics and point toward cytoskeleton remodeling as an important effector of the motion of mineral-bearing vesicles.

Biomineralization is a widespread, fundamental process by which organisms use minerals to harden their tissues. Mineral-bearing vesicles were observed in biomineralizing cells and play an essential role in biomineralization, yet little is known about their three-dimensional (3D) dynamics. Here we quantify 3D-vesicle-dynamics during calcite skeleton formation in sea urchin larvae, using lattice-light-sheet microscopy. We discover that calcium vesicles perform a diffusive motion in both calcifying and non-calcifying cells of the embryo. The diffusion coefficient and vesicle speed are higher in the mesenchymal skeletogenic cells compared to the epithelial ectodermal cells. This difference is possibly due to the higher rigidity of the ectodermal cells as demonstrated by the enhanced signal of f-actin and myosinII activity in these cells compared to the skeletogenic cells. The motion of the vesicles in the skeletogenic cells, is not directed toward the biomineralization compartment but the vesicles slow down near it, possibly to deposit their content. Blocking skeletogenesis through the inhibition of Vascular Endothelial Growth Factor Receptor (VEGFR), increases vesicle volume but doesn’t change the diffusion mode and the cytoskeleton markers in the cells. Our studies reveal the active diffusive motion of mineral bearing vesicles that is apparently defined by the mechanical properties of the cells.

VEGF signaling activates the matrix metalloproteinases, MmpL7 and MmpL5 at the sites of active skeletal growth and MmpL7 regulates skeletal elongation

Organisms can uptake minerals, shape them in different forms and generate teeth, skeletons or shells that support and protect them. Mineral uptake, trafficking and nucleation are tightly regulated by the biomineralizing cells through networks of specialized proteins. Specifically, matrix metalloproteases (MMPs) digest various extracellular substrates and allow for mineralization in the vertebrates' teeth and bones, but little is known about their role in invertebrates' systems. The sea urchin embryo provides an excellent invertebrate model for genetic and molecular studies of biomineralization. MMP inhibition prevents the growth of the calcite spicules of the sea urchin larval skeleton, however, the molecular mechanisms and genes that underlie this response are not well understood. Here we study the spatial expression and regulation of two membrane type MMPs that were found to be occluded in the sea urchin spicules, Pl-MmpL7 and Pl-MmpL5, and investigate the function of Pl-MmpL7 in skeletogenesis. The inhibition of MMPs does not change the volume of the calcium vesicles in the skeletogenic cells. The expression of Pl-MmpL7 and Pl-MmpL5 is regulated by the Vascular Endothelial Growth Factor (VEGF) signaling, from the time of skeleton initiation and on. The expression of these genes is localized to the subsets of skeletogenic cells where active spicule growth occurs throughout skeletogenesis. Downregulation of Pl-MmpL7 expression delays the growth of the skeletal rods and in some cases, strongly perturbs skeletal shape. The localized expression of Pl-MmpL7 and Pl-MmpL5 to the active growth zone and the effect of Pl-MmpL7 perturbations on skeletal growth, suggest that these genes are essential for normal spicule elongation in the sea urchin embryo.

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.

Cellular bicarbonate accumulation and vesicular proton transport promote calcification in the sea urchin larva

The sea urchin embryo develops a calcitic endoskeleton through intracellular formation of amorphous calcium carbonate (ACC). Intracellular precipitation of ACC, requires [Formula: see text] concentrating as well as proton export mechanisms to promote calcification. These processes are of fundamental importance in biological mineralization, but remain largely unexplored. Here, we demonstrate that the calcifying primary mesenchyme cells (PMCs) use Na+/H+-exchange (NHE) mechanisms to control cellular pH homeostasis during maintenance of the skeleton. During skeleton re-calcification, pHi of PMCs is increased accompanied by substantial elevation in intracellular [Formula: see text] mediated by the [Formula: see text] cotransporter Sp_Slc4a10. However, PMCs lower their pHi regulatory capacities associated with a reduction in NHE activity. Live-cell imaging using green fluorescent protein reporter constructs in combination with intravesicular pH measurements demonstrated alkaline and acidic populations of vesicles in PMCs and extensive trafficking of large V-type H+-ATPase (VHA)-rich acidic vesicles in blastocoelar filopodial cells. Pharmacological and gene expression analyses underline a central role of the VHA isoforms Sp_ATP6V0a1, Sp_ATP6V01_1 and Sp_ATPa1-4 for the process of skeleton re-calcification. These results highlight novel pH regulatory strategies in calcifying cells of a marine species with important implications for our understanding of the mineralization process in times of rapid changes in oceanic pH.

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.

PI3K inhibition highlights new molecular interactions involved in the skeletogenesis of Paracentrotus lividus embryos

The sea urchin embryo develops a well-defined biomineralized endoskeleton, synthesized exclusively by the skeletogenic cells, supported by ectodermal cues for the correct skeleton patterning. The biomineralization process is tightly regulated via a hierarchical order of gene expression, including transcription and growth factors, biomineralization proteins. Recently, the role of kinases and intracellular signaling pathways in sea urchin skeletogenesis has been addressed, although the downstream components still remain unknown. In this study, we investigated the role of phosphatidylinositide 3-kinase (PI3K)-mediated signaling pathway in Paracentrotus lividus, to identify its genes/proteins targets. The effects of LY294002 (LY), a PI3K-specific inhibitor, were evaluated at morphological and molecular levels. Treatment with 40 μM LY from the blastula stage completely blocked skeleton deposition, which was reversed by wash out experiments. Besides, LY caused a slight delay in the tripartite gut development. Despite the skeleton absence, a few skeleton-specific proteins/mRNAs were regularly expressed and localized in LY-treated embryos, as shown for MSP130 and SM50 by immunofluorescence and in situ hybridization experiments. QPCR analyses showed that LY differently affected the expression of genes coding for other biomineralization proteins, transcription and growth factors. SM30 and carbonic anhydrase expression was severely downregulated, while almost all the transcription factors analyzed were upregulated. Based on the present results and in silico analyses, we propose an "interactomic" model simulating PI3K connections in P. lividus embryos. Our findings define a novel regulatory step in the embryonic skeletogenesis, and provide valuable molecular data for further studies on the role of PI3K signaling in invertebrate biomineralization.

Spatially mapping gene expression in sea urchin primary mesenchyme cells

During sea urchin embryogenesis, primary mesenchyme cells (PMCs) follow a stereotypical migratory program, arrange into a primary pattern, then begin to secrete a bilaterally symmetric calcium carbonate skeleton. Recently identified genes are expressed in spatially-restricted domains within the PMC population (Sun & Ettensohn, 2014). To better understand the molecular mechanisms orchestrating PMC positioning, we are characterizing the expression profiles of PMC subset-specific genes. To deconvolve the spatiotemporal expression patterns within PMCs, we detect cell-specific mRNA expression with combined RNA fluorescence in situ hybridization and immunolabeling of PMCs. Subsequent confocal microscopy provides 3D position and expression information for individual PMCs. We extract PMC positions and relative gene expression levels, then model these results using open-source 3D modeling software. This versatile protocol can be extended to other models and systems.

Genomic resources for the study of echinoderm development and evolution

Echinoderms are important research models for a wide range of biological questions. In particular, echinoderm embryos are exemplary models for dissecting the molecular and cellular processes that drive development and testing how these processes can be modified through evolution to produce the extensive morphological diversity observed in the phylum. Modern attempts to characterize these processes depend on some level of genomic analysis; from querying annotated gene sets to functional genomics experiments to identify candidate cis-regulatory sequences. Given how essential these data have become, it is important that researchers using available datasets or performing their own genome-scale experiments understand the nature and limitations of echinoderm genomic analyses. In this chapter we highlight the current state of echinoderm genomic data and provide methodological considerations for common approaches, including analysis of transcriptome and functional genomics datasets.

Ectomesoderm and epithelial-mesenchymal transition-related genes in spiralian development

BACKGROUND

Spiralians (e.g., annelids, molluscs, and flatworms) possess two sources of mesoderm. One is from endodermal precursors (endomesoderm), which is considered to be the ancestral source in metazoans. The second is from ectoderm (ectomesoderm) and may represent a novel cell type in the Spiralia. In the mollusc Crepidula fornicata, ectomesoderm is derived from micromere daughters within the A and B cell quadrants. Their progeny lie along the anterolateral edges of the blastopore. There they undergo epithelial-mesenchymal transition (EMT), become rounded and undergo delamination/ingression. Subsequently, they assume the mesenchymal phenotype, and migrate beneath the surface ectoderm to differentiate various cell types, including muscles and pigment cells.

RESULTS

We examined expression of several genes whose homologs are known to regulate Type 1 EMT in other metazoans. Most of these genes were expressed within spiralian ectomesoderm during EMT.

CONCLUSIONS

We propose that spiralian ectomesoderm, which exhibits analogous cellular behaviors to other populations of mesenchymal cells, may be controlled by the same genes that drive EMT in other metazoans. Perhaps these genes comprise a conserved metazoan EMT gene regulatory network (GRN). This study represents the first step in elucidating the GRN controlling the development of a novel spiralian cell type (ectomesoderm). Developmental Dynamics 247:1097-1120, 2018. © 2018 Wiley Periodicals, Inc.

From genome to anatomy: The architecture and evolution of the skeletogenic gene regulatory network of sea urchins and other echinoderms

The skeletogenic gene regulatory network (GRN) of sea urchins and other echinoderms is one of the most intensively studied transcriptional networks in any developing organism. As such, it serves as a preeminent model of GRN architecture and evolution. This review summarizes our current understanding of this developmental network. We describe in detail the most comprehensive model of the skeletogenic GRN, one developed for the euechinoid sea urchin Strongylocentrotus purpuratus, including its initial deployment by maternal inputs, its elaboration and stabilization through regulatory gene interactions, and its control of downstream effector genes that directly drive skeletal morphogenesis. We highlight recent comparative studies that have leveraged the euechinoid GRN model to examine the evolution of skeletogenic programs in diverse echinoderms, studies that have revealed both conserved and divergent features of skeletogenesis within the phylum. Last, we summarize the major insights that have emerged from analysis of the structure and evolution of the echinoderm skeletogenic GRN and identify key, unresolved questions as a guide for future work.

Genome-wide use of high- and low-affinity Tbrain transcription factor binding sites during echinoderm development

Sea stars and sea urchins are model systems for interrogating the types of deep evolutionary changes that have restructured developmental gene regulatory networks (GRNs). Although cis-regulatory DNA evolution is likely the predominant mechanism of change, it was recently shown that Tbrain, a Tbox transcription factor protein, has evolved a changed preference for a low-affinity, secondary binding motif. The primary, high-affinity motif is conserved. To date, however, no genome-wide comparisons have been performed to provide an unbiased assessment of the evolution of GRNs between these taxa, and no study has attempted to determine the interplay between transcription factor binding motif evolution and GRN topology. The study here measures genome-wide binding of Tbrain orthologs by using ChIP-sequencing and associates these orthologs with putative target genes to assess global function. Targets of both factors are enriched for other regulatory genes, although nonoverlapping sets of functional enrichments in the two datasets suggest a much diverged function. The number of low-affinity binding motifs is significantly depressed in sea urchins compared with sea star, but both motif types are associated with genes from a range of functional categories. Only a small fraction (∼10%) of genes are predicted to be orthologous targets. Collectively, these data indicate that Tbr has evolved significantly different developmental roles in these echinoderms and that the targets and the binding motifs in associated cis-regulatory sequences are dispersed throughout the hierarchy of the GRN, rather than being biased toward terminal process or discrete functional blocks, which suggests extensive evolutionary tinkering.

SoxB2 in sea urchin development: implications in neurogenesis, ciliogenesis and skeletal patterning

Background

Current studies in evolutionary developmental biology are focused on the reconstruction of gene regulatory networks in target animal species. From decades, the scientific interest on genetic mechanisms orchestrating embryos development has been increasing in consequence to the fact that common features shared by evolutionarily distant phyla are being clarified. In 2011, a study across eumetazoan species showed for the first time the existence of a highly conserved non-coding element controlling the SoxB2 gene, which is involved in the early specification of the nervous system. This discovery raised several questions about SoxB2 function and regulation in deuterostomes from an evolutionary point of view.

Results

Due to the relevant phylogenetic position within deuterostomes, the sea urchin Strongylocentrotus purpuratus represents an advantageous animal model in the field of evolutionary developmental biology. Herein, we show a comprehensive study of SoxB2 functions in sea urchins, in particular its expression pattern in a wide range of developmental stages, and its co-localization with other neurogenic markers, as SoxB1, SoxC and Elav. Moreover, this work provides a detailed description of the phenotype of sea urchin SoxB2 knocked-down embryos, confirming its key function in neurogenesis and revealing, for the first time, its additional roles in oral and aboral ectoderm cilia and skeletal rod morphology.

Conclusions

We concluded that SoxB2 in sea urchins has a neurogenic function; however, this gene could have multiple roles in sea urchin embryogenesis, expanding its expression in non-neurogenic cells. We showed that SoxB2 is functionally conserved among deuterostomes and suggested that in S. purpuratus this gene acquired additional functions, being involved in ciliogenesis and skeletal patterning.

Electronic supplementary material

The online version of this article (10.1186/s13227-018-0094-1) contains supplementary material, which is available to authorized users.

Thyroid Hormones Accelerate Initiation of Skeletogenesis via MAPK (ERK1/2) in Larval Sea Urchins (Strongylocentrotus purpuratus)

Thyroid hormones are important regulators of development and metabolism in animals. Their function via genomic and non-genomic actions is well-established in vertebrate species but remains largely elusive among invertebrates. Previous work suggests that thyroid hormones, principally 3,5,3',5'-Tetraiodo-L-thyronine (T4), regulate development to metamorphosis in sea urchins. Here we show that thyroid hormones, including T4, 3,5,3'-triiodo-l-thyronine (T3), and 3,5-Diiodothyronine (T2) accelerate initiation of skeletogenesis in sea urchin gastrulae and pluteus larvae of the sea urchin Strongylocentrotus purpuratus, as measured by skeletal spicule formation. Fluorescently conjugated hormones show T4 binding to primary mesenchyme cells in sea urchin gastrulae. Furthermore, our investigation of TH mediated skeletogenesis shows that Ets1, a transcription factor controlling initiation of skeletogenesis, is a target of activated (phosphorylated) mitogen-activated protein kinase [MAPK; extracellular signal-regulated kinase 1/2 (ERK1/2)]. As well, we show that PD98059, an inhibitor of ERK1/2 MAPK signaling, prevents the T4 mediated acceleration of skeletogenesis and upregulation of Ets1. In contrast, SB203580, an inhibitor of p38 MAPK signaling, did not inhibit the effect of T4. Immunohistochemistry revealed that T4 causes phosphorylation of ERK1/2 in presumptive primary mesenchyme cells and the basal membrane of epithelial cells in the gastrula. Pre-incubation of sea urchin gastrulae with RGD peptide, a competitive inhibitor of TH binding to integrins, inhibited the effect of T4 on skeletogenesis. Together, these experiments provide evidence that T4 acts via a MAPK- (ERK1/2) mediated integrin membrane receptor to accelerate skeletogenesis in sea urchin mesenchyme cells. These findings shed light, for the first time, on a putative non-genomic pathway of TH action in a non-chordate deuterostome and help elucidate the evolutionary history of TH signaling in animals.

The sea cucumber genome provides insights into morphological evolution and visceral regeneration

Apart from sharing common ancestry with chordates, sea cucumbers exhibit a unique morphology and exceptional regenerative capacity. Here we present the complete genome sequence of an economically important sea cucumber, A. japonicus, generated using Illumina and PacBio platforms, to achieve an assembly of approximately 805 Mb (contig N50 of 190 Kb and scaffold N50 of 486 Kb), with 30,350 protein-coding genes and high continuity. We used this resource to explore key genetic mechanisms behind the unique biological characters of sea cucumbers. Phylogenetic and comparative genomic analyses revealed the presence of marker genes associated with notochord and gill slits, suggesting that these chordate features were present in ancestral echinoderms. The unique shape and weak mineralization of the sea cucumber adult body were also preliminarily explained by the contraction of biomineralization genes. Genome, transcriptome, and proteome analyses of organ regrowth after induced evisceration provided insight into the molecular underpinnings of visceral regeneration, including a specific tandem-duplicated prostatic secretory protein of 94 amino acids (PSP94)-like gene family and a significantly expanded fibrinogen-related protein (FREP) gene family. This high-quality genome resource will provide a useful framework for future research into biological processes and evolution in deuterostomes, including remarkable regenerative abilities that could have medical applications. Moreover, the multiomics data will be of prime value for commercial sea cucumber breeding programs.

Echinoderms, ubiquitous in the marine environment, are important from evolutionary, ecological, and socioeconomic perspectives. Together with chordates and hemichordates, they form the deuterostome clade, making them a crucial node in the study of chordate ancestry. Within echinoderms, class Holothuroidea is unique; its members (the sea cucumbers) display remarkable regenerative abilities and play key roles as sediment bioturbators and symbiotic hosts, and many are prized in the seafood and pharmaceutical industries. The sea cucumber genome therefore has the potential to significantly contribute to our understanding of important evolutionary and biological processes and help enhance aquaculture programs. Here we present a high-quality genome sequence for the economically important species Apostichopus japonicus. Through comparative analyses, we identified 763 echinoderm-specific gene families enriched in genes encoding membrane proteins, ion channels, and signal transduction proteins. Marker genes associated with the notochord and gill slits were also found, providing valuable insight into the origin of chordates. The reduced number and low expression levels of biomineralization genes reflect the skeletal degeneration seen in sea cucumbers. Importantly, 2 gene families appeared to be expanded in A. japonicus and may play crucial roles in its heightened regenerative potential. Together, findings from the sea cucumber genome provide important and novel insights into echinoderm and deuterostome biology.

Endocytosis in primary mesenchyme cells during sea urchin larval skeletogenesis

The sea urchin larval embryo elaborates two calcitic endoskeletal elements called spicules. Spicules are synthesized by the primary mesenchyme cells (PMCs) and begin to form at early gastrula stage. It is known that the calcium comprising the spicules comes from the seawater and we wish to further consider the mode of calcium transport from the extracellular seawater to the PMCs and then onto the forming spicules. We used PMC in vitro cultures, calcein, fluorescently labeled dextran, and fluorescently labeled Wheat Germ Agglutinin (WGA) to track calcium transport from the seawater into PMCs and spicules and to determine how molecules from the surface of PMCs interact with the incoming calcium. Labeling of PMC endocytic vesicles and forming spicules by both calcein and fluorescently tagged dextran indicate that calcium is taken up from the seawater by endocytosis and directly incorporated into spicules. Calcein labeling studies also indicate that calcium from the extracellular seawater begins to be incorporated into spicules within 30min of uptake. In addition, we demonstrate that fluorescently labeled WGA and calcein are taken up by many of the same endocytic vesicles and are incorporated into growing spicules. These findings suggest that PMC specific surface molecules accompany calcium ions as they enter PMCs via endocytosis and are incorporated together in the growing spicule. Using anti-spicule matrix protein antibodies, we pinpoint a subset of spicule matrix proteins that may accompany calcium ions from the surface of the PMCs until they are incorporated into spicules. Msp130 is identified as one of these spicule matrix proteins.

Biomineralization-related specialization of hemocytes and mantle tissues of the Pacific oyster Crassostrea gigas

The molluscan exoskeleton (shell) plays multiple important roles including structural support, protection from predators and stressors, and physiological homeostasis. Shell formation is a tightly regulated biological process that allows molluscs to build their shells even in environments unfavorable for mineral precipitation. Outer mantle edge epithelial cells (OME) and hemocytes were implicated in this process; however, the exact functions of these cell types in biomineralization are not clear. Pacific oysters (Crassostrea gigas) were used to study differences in the expression profiles of selected biomineralization-related genes in hemocytes and mantle cells, and the functional characteristics of hemocytes such as adhesion, motility and phagocytosis. The specialized role of OME in shell formation was supported by high expression levels of the extracellular matrix (ECM) related and cell-cell interaction genes. Density gradient separation of hemocytes revealed distinct phenotypes based on the cell morphology, gene expression patterns, motility and adhesion characteristics. These hemocyte fractions can be categorized into two functional groups, i.e. biomineralization and immune response cells. Gene expression profiles of the putative biomineralizing hemocytes indicate that in addition to their proposed role in mineral transport, hemocytes also contribute to the formation of the ECM, thus challenging the current paradigm of the mantle as the sole source of the ECM for shell formation. Our findings corroborate the specialized roles of hemocytes and the OME in biomineralization and emphasize complexity of the biological controls over shell formation in bivalves.

The skeletal proteome of the sea star Patiria miniata and evolution of biomineralization in echinoderms

BACKGROUND

Proteomic studies of skeletal proteins have revealed large, complex mixtures of proteins occluded within the mineral. Many skeletal proteomes contain rapidly evolving proteins with repetitive domains, further complicating our understanding. In echinoderms, proteomic analysis of the skeletal proteomes of mineralized tissues of the sea urchin Strongylocentrotus purpuratus prominently featured spicule matrix proteins with repetitive sequences linked to a C-type lectin domain. A comparative study of the brittle star Ophiocoma wendtii skeletal proteome revealed an order of magnitude fewer proteins containing C-type lectin domains. A number of other proteins conserved in the skeletons of the two groups were identified. Here we report the complete skeletal proteome of the sea star Patiria miniata and compare it to that of the other echinoderm groups.

RESULTS

We have identified eighty-five proteins in the P. miniata skeletal proteome. Forty-two percent of the proteins were determined to be homologous to proteins found in the S. purpuratus skeletal proteomes. An additional 34 % were from similar functional classes as proteins in the urchin proteomes. Thirteen percent of the P. miniata proteins had homologues in the O. wendtii skeletal proteome with an additional 29% showing similarity to brittle star skeletal proteins. The P. miniata skeletal proteome did not contain any proteins with C-lectin domains or with acidic repetitive regions similar to the sea urchin or brittle star spicule matrix proteins. MSP130 proteins were also not found. We did identify a number of proteins homologous between the three groups. Some of the highly conserved proteins found in echinoderm skeletons have also been identified in vertebrate skeletons.

CONCLUSIONS

The presence of proteins conserved in the skeleton in three different echinoderm groups indicates these proteins are important in skeleton formation. That a number of these proteins are involved in skeleton formation in vertebrates suggests a common origin for some of the fundamental processes co-opted for skeleton formation in deuterostomes. The proteins we identify suggest transport of proteins and calcium via endosomes was co-opted to this function in a convergent fashion. Our data also indicate that modifications to the process of skeleton formation can occur through independent co-option of proteins following species divergence as well as through domain shuffling.

A Conserved Role for VEGF Signaling in Specification of Homologous Mesenchymal Cell Types Positioned at Spatially Distinct Developmental Addresses in Early Development of Sea Urchins

Comparative studies of early development in echinoderms are revealing the tempo and mode of alterations to developmental gene regulatory networks and to the cell types they specify. In euechinoid sea urchins, skeletogenic mesenchyme (SM) ingresses prior to gastrulation at the vegetal pole and aligns into a ring-like array with two bilateral pockets of cells, the sites where spiculogenesis will later occur. In cidaroid sea urchins, the anciently diverged sister clade to euechinoid sea urchins, a homologous SM cell type ingresses later in development, after gastrulation has commenced, and consequently at a distinct developmental address. Thus, a heterochronic shift of ingression of the SM cell type occurred in one of the echinoid lineages. In euechinoids, specification and migration of SM are facilitated by vascular endothelial growth factor (VEGF) signaling. We describe spatiotemporal expression of vegf and vegfr and experimental manipulations targeting VEGF signaling in the cidaroid Eucidaris tribuloides. Spatially, vegf and vegfr mRNA localizes similarly as in euechinoids, suggesting conserved deployment in echinoids despite their spatially distinct development addresses of ingression. Inhibition of VEGF signaling in E. tribuloides suggests its role in SM specification is conserved in echinoids. Temporal discrepancies between the onset of vegf expression and SM ingression likely result in previous observations of SM "random wandering" behavior. Our results indicate that, although the SM cell type in echinoids ingresses into distinct developmental landscapes, it retains a signaling mechanism that restricts their spatial localization to a conserved developmental address where spiculogenesis later occurs.

TGF-β sensu stricto signaling regulates skeletal morphogenesis in the sea urchin embryo

Cell-cell signaling plays a prominent role in the formation of the embryonic skeleton of sea urchins, but the mechanisms are poorly understood. In the present study, we uncover an essential role for TGF-β sensu stricto signaling in this process. We show that TgfbrtII, a type II receptor dedicated to signaling through TGF-β sensu stricto, is expressed selectively in skeletogenic primary mesenchyme cells (PMCs) during skeleton formation. Morpholino (MO) knockdowns and studies with a specific TgfbrtII inhibitor (ITD-1) in both S. purpuratus and Lytechinus variegatus embryos show that this receptor is required for biomineral deposition. We provide pharmacological evidence that Alk4/5/7 is the cognate TGF-β type I receptor that pairs with TgfbrtII and show by inhibitor treatments of isolated micromeres cultured in vitro that both Alk4/5/7 and TgfbrtII function cell-autonomously in PMCs. Gene expression and gene knockdown studies suggest that TGF-β sensu stricto may be the ligand that interacts with TgfbrtII and support the view that this TGF-β superfamily ligand provides an essential, permissive cue for skeletogenesis, although it is unlikely to provide spatial patterning information. Taken together, our findings reveal that this model morphogenetic process involves an even more diverse suite of cell signaling pathways than previously appreciated and show that PMCs integrate a complex set of both generalized and spatially localized cues in assembling the endoskeleton.

Experimental Approach Reveals the Role of alx1 in the Evolution of the Echinoderm Larval Skeleton

Over the course of evolution, the acquisition of novel structures has ultimately led to wide variation in morphology among extant multicellular organisms. Thus, the origins of genetic systems for new morphological structures are a subject of great interest in evolutionary biology. The larval skeleton is a novel structure acquired in some echinoderm lineages via the activation of the adult skeletogenic machinery. Previously, VEGF signaling was suggested to have played an important role in the acquisition of the larval skeleton. In the present study, we compared expression patterns of Alx genes among echinoderm classes to further explore the factors involved in the acquisition of a larval skeleton. We found that the alx1 gene, originally described as crucial for sea urchin skeletogenesis, may have also played an essential role in the evolution of the larval skeleton. Unlike those echinoderms that have a larval skeleton, we found that alx1 of starfish was barely expressed in early larvae that have no skeleton. When alx1 overexpression was induced via injection of alx1 mRNA into starfish eggs, the expression patterns of certain genes, including those possibly involved in skeletogenesis, were altered. This suggested that a portion of the skeletogenic program was induced solely by alx1. However, we observed no obvious external phenotype or skeleton. We concluded that alx1 was necessary but not sufficient for the acquisition of the larval skeleton, which, in fact, requires several genetic events. Based on these results, we discuss how the larval expression of alx1 contributed to the acquisition of the larval skeleton in the putative ancestral lineage of echinoderms.

microRNA-31 modulates skeletal patterning in the sea urchin embryo

MicroRNAs (miRNAs) are small non-coding RNAs that repress the translation and reduce the stability of target mRNAs in animal cells. microRNA-31 (miR-31) is known to play a role in cancer, bone formation and lymphatic development. However, studies to understand the function of miR-31 in embryogenesis have been limited. We examined the regulatory role of miR-31 in early development using the sea urchin as a model. miR-31 is expressed at all stages of development and its knockdown (KD) disrupts the patterning and function of primary mesenchyme cells (PMCs), which form the embryonic skeleton spicules. We identified that miR-31 directly represses Pmar1, Alx1, Snail and VegfR7 within the PMC gene regulatory network using reporter constructs. Further, blocking the miR-31-mediated repression of Alx1 and/or VegfR7 in the developing embryo resulted in defects in PMC patterning and skeletogenesis. The majority of the mislocalized PMCs in miR-31 KD embryos did not express VegfR10, indicating that miR-31 regulates VegfR gene expression within PMCs. In addition, miR-31 indirectly suppresses Vegf3 expression in the ectoderm. These results indicate that miR-31 coordinately suppresses genes within the PMCs and in the ectoderm to impact PMC patterning and skeletogenesis. This study identifies the novel function and molecular mechanism of miR-31-mediated regulation in the developing embryo.