Genome-wide analysis of the skeletogenic gene regulatory network of sea urchins

Single-Cell Transcriptomics Reveals Evolutionary Reconfiguration of Embryonic Cell Fate Specification in the Sea Urchin Heliocidaris erythrogramma

Altered regulatory interactions during development likely underlie a large fraction of phenotypic diversity within and between species, yet identifying specific evolutionary changes remains challenging. Analysis of single-cell developmental transcriptomes from multiple species provides a powerful framework for unbiased identification of evolutionary changes in developmental mechanisms. Here, we leverage a "natural experiment" in developmental evolution in sea urchins, where a major life history switch recently evolved in the lineage leading to Heliocidaris erythrogramma, precipitating extensive changes in early development. Comparative analyses of single-cell transcriptome analysis (scRNA-seq) developmental time courses from H. erythrogramma and Lytechinus variegatus (representing the derived and ancestral states, respectively) reveal numerous evolutionary changes in embryonic patterning. The earliest cell fate specification events and the primary signaling center are co-localized in the ancestral developmental gene regulatory network; remarkably, in H. erythrogramma, they are spatially and temporally separate. Fate specification and differentiation are delayed in most embryonic cell lineages, although in some cases, these processes are conserved or even accelerated. Comparative analysis of regulator-target gene co-expression is consistent with many specific interactions being preserved but delayed in H. erythrogramma, while some otherwise widely conserved interactions have likely been lost. Finally, specific patterning events are directly correlated with evolutionary changes in larval morphology, suggesting that they are directly tied to the life history shift. Together, these findings demonstrate that comparative scRNA-seq developmental time courses can reveal a diverse set of evolutionary changes in embryonic patterning and provide an efficient way to identify likely candidate regulatory interactions for subsequent experimental validation.

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.

Shared regulatory function of non-genomic thyroid hormone signaling in echinoderm skeletogenesis

Thyroid hormones are crucial regulators of metamorphosis and development in bilaterians, particularly in chordate deuterostomes. Recent evidence suggests a role for thyroid hormone signaling, principally via 3,5,3',5'-Tetraiodo-L-thyronine (T4), in the regulation of metamorphosis, programmed cell death and skeletogenesis in echinoids (sea urchins and sand dollars) and sea stars. Here, we test whether TH signaling in skeletogenesis is a shared trait of Echinozoa (Echinoida and Holothouroida) and Asterozoa (Ophiourida and Asteroida). We demonstrate dramatic acceleration of skeletogenesis after TH treatment in three classes of echinoderms: sea urchins, sea stars, and brittle stars (echinoids, asteroids, and ophiuroids). Fluorescently labeled thyroid hormone analogues reveal thyroid hormone binding to cells proximal to regions of skeletogenesis in the gut and juvenile rudiment. We also identify, for the first time, a potential source of thyroxine during gastrulation in sea urchin embryos. Thyroxine-positive cells are present in tip of the archenteron. In addition, we detect thyroid hormone binding to the cell membrane and nucleus during metamorphic development in echinoderms. Immunohistochemistry of phosphorylated MAPK in the presence and absence of TH-binding inhibitors suggests that THs may act via phosphorylation of MAPK (ERK1/2) to accelerate initiation of skeletogenesis in the three echinoderm groups. Together, these results indicate that TH regulation of mesenchyme cell activity via integrin-mediated MAPK signaling may be a conserved mechanism for the regulation of skeletogenesis in echinoderm development. In addition, TH action via a nuclear thyroid hormone receptor may regulate metamorphic development. Our findings shed light on potentially ancient pathways of thyroid hormone activity in echinoids, ophiuroids, and asteroids, or on a signaling system that has been repeatedly co-opted to coordinate metamorphic development in bilaterians.

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.

Semi-automated, high-content imaging of drug transporter knockout sea urchin (Lytechinus pictus) embryos

A defining feature of sea urchins is their extreme fecundity. Urchins produce millions of transparent, synchronously developing embryos, ideal for spatial and temporal analysis of development. This biological feature has been effectively utilized for ensemble measurement of biochemical changes. However, it has been underutilized in imaging studies, where single embryo measurements are used. Here we present an example of how stable genetics and high content imaging, along with machine learning-based image analysis, can be used to exploit the fecundity and synchrony of sea urchins in imaging-based drug screens. Building upon our recently created sea urchin ABCB1 knockout line, we developed a high-throughput assay to probe the role of this drug transporter in embryos. We used high content imaging to compare accumulation and toxicity of canonical substrates and inhibitors of the transporter, including fluorescent molecules and antimitotic cancer drugs, in homozygous knockout and wildtype embryos. To measure responses from the resulting image data, we used a nested convolutional neural network, which rapidly classified embryos according to fluorescence or cell division. This approach identified sea urchin embryos with 99.8% accuracy and determined two-cell and aberrant embryos with 96.3% and 89.1% accuracy, respectively. The results revealed that ABCB1 knockout embryos accumulated the transporter substrate calcein 3.09 times faster than wildtypes. Similarly, knockouts were 4.71 and 3.07 times more sensitive to the mitotic poisons vinblastine and taxol. This study paves the way for large scale pharmacological screens in the sea urchin embryo.

Single-cell transcriptomics reveals evolutionary reconfiguration of embryonic cell fate specification in the sea urchin Heliocidaris erythrogramma

Altered regulatory interactions during development likely underlie a large fraction of phenotypic diversity within and between species, yet identifying specific evolutionary changes remains challenging. Analysis of single-cell developmental transcriptomes from multiple species provides a powerful framework for unbiased identification of evolutionary changes in developmental mechanisms. Here, we leverage a "natural experiment" in developmental evolution in sea urchins, where a major life history switch recently evolved in the lineage leading to Heliocidaris erythrogramma, precipitating extensive changes in early development. Comparative analyses of scRNA-seq developmental time courses from H. erythrogramma and Lytechinus variegatus (representing the derived and ancestral states respectively) reveals numerous evolutionary changes in embryonic patterning. The earliest cell fate specification events, and the primary signaling center are co-localized in the ancestral dGRN but remarkably, in H. erythrogramma they are spatially and temporally separate. Fate specification and differentiation are delayed in most embryonic cell lineages, although in some cases, these processes are conserved or even accelerated. Comparative analysis of regulator-target gene co-expression is consistent with many specific interactions being preserved but delayed in H. erythrogramma, while some otherwise widely conserved interactions have likely been lost. Finally, specific patterning events are directly correlated with evolutionary changes in larval morphology, suggesting that they are directly tied to the life history shift. Together, these findings demonstrate that comparative scRNA-seq developmental time courses can reveal a diverse set of evolutionary changes in embryonic patterning and provide an efficient way to identify likely candidate regulatory interactions for subsequent experimental validation.

microRNA-1 regulates sea urchin skeletogenesis by directly targeting skeletogenic genes and modulating components of signaling pathways

microRNAs are evolutionarily conserved non-coding RNAs that direct post-transcriptional regulation of target transcripts. In vertebrates, microRNA-1 (miR-1) is expressed in muscle and has been found to play critical regulatory roles in vertebrate angiogenesis, a process that has been proposed to be analogous to sea urchin skeletogenesis. Results indicate that both miR-1 inhibitor and miR-1 mimic-injected larvae have significantly less F-actin enriched circumpharyngeal muscle fibers and fewer gut contractions. In addition, miR-1 regulates the positioning of skeletogenic primary mesenchyme cells (PMCs) and skeletogenesis of the sea urchin embryo. Interestingly, the gain-of-function of miR-1 leads to more severe PMC patterning and skeletal branching defects than its loss-of-function. The results suggest that miR-1 directly suppresses Ets1/2, Tbr, and VegfR7 of the skeletogenic gene regulatory network, and Nodal, and Wnt1 signaling components. This study identifies potential targets of miR-1 that impacts skeletogenesis and muscle formation and contributes to a deeper understanding of miR-1's function during development.

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.

Solute carrier (SLC) expression reveals skeletogenic cell diversity

Within the developing embryo is a microcosm of cell type diversity. Single cell RNA-sequencing (scRNA-seq) is used to reveal cell types, typically by grouping cells according to their gene regulatory states. However, both across and within these regulatory states are additional layers of cellular diversity represented by the differential expression of genes that govern cell function. Here, we analyzed scRNA-seq data representing the late gastrula stage of Strongylocentrotus purpuratus (purple sea urchin) to understand the patterning of transporters belonging to the ABC and SLC families. These transporters handle diverse substrates from amino acids to signaling molecules, nutrients and xenobiotics. Using transporter-based clustering, we identified unique transporter patterns that are both shared across cell lineages, as well as those that were unique to known cell types. We further explored three patterns of transporter expression in mesodermal cells including secondary mesenchyme cells (pigment cells and blastocoelar cells) and skeletogenic cells (primary mesenchyme cells). The results revealed the enrichment of SMTs potentially involved in nutrient absorption (SLC5A9, SLC7A11, SLC35F3, and SLC52A3) and skeletogenesis (SLC9A3, SLC13A2/3/5, and SLC39A13) in pigment cells and blastocoelar cells respectively. The results indicated that the strategy of clustering by cellular activity can be useful for discovering cellular populations that would otherwise remain obscured.

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.

Thyroid hormone membrane receptor binding and transcriptional regulation in the sea urchin Strongylocentrotus purpuratus

Thyroid hormones (THs) are small amino acid derived signaling molecules with broad physiological and developmental functions in animals. Specifically, their function in metamorphic development, ion regulation, angiogenesis and many others have been studied in detail in mammals and some other vertebrates. Despite extensive reports showing pharmacological responses of invertebrate species to THs, little is known about TH signaling mechanisms outside of vertebrates. Previous work in sea urchins suggests that non-genomic mechanisms are activated by TH ligands. Here we show that several THs bind to sea urchin (Strongylocentrotus purpuratus) cell membrane extracts and are displaced by ligands of RGD-binding integrins. A transcriptional analysis across sea urchin developmental stages shows activation of genomic and non-genomic pathways in response to TH exposure, suggesting that both pathways are activated by THs in sea urchin embryos and larvae. We also provide evidence associating TH regulation of gene expression with TH response elements in the genome. In ontogeny, we found more differentially expressed genes in older larvae compared to gastrula stages. In contrast to gastrula stages, the acceleration of skeletogenesis by thyroxine in older larvae is not fully inhibited by competitive ligands or inhibitors of the integrin membrane receptor pathway, suggesting that THs likely activate multiple pathways. Our data confirms a signaling function of THs in sea urchin development and suggests that both genomic and non-genomic mechanisms play a role, with genomic signaling being more prominent during later stages of larval development.

Early expression onset of tissue-specific effector genes during the specification process in sea urchin embryos

BACKGROUND

In the course of animal developmental processes, various tissues are differentiated through complex interactions within the gene regulatory network. As a general concept, differentiation has been considered to be the endpoint of specification processes. Previous works followed this view and provided a genetic control scheme of differentiation in sea urchin embryos: early specification genes generate distinct regulatory territories in an embryo to express a small set of differentiation driver genes; these genes eventually stimulate the expression of tissue-specific effector genes, which provide biological identity to differentiated cells, in each region. However, some tissue-specific effector genes begin to be expressed in parallel with the expression onset of early specification genes, raising questions about the simplistic regulatory scheme of tissue-specific effector gene expression and the current concept of differentiation itself.

RESULTS

Here, we examined the dynamics of effector gene expression patterns during sea urchin embryogenesis. Our transcriptome-based analysis indicated that many tissue-specific effector genes begin to be expressed and accumulated along with the advancing specification GRN in the distinct cell lineages of embryos. Moreover, we found that the expression of some of the tissue-specific effector genes commences before cell lineage segregation occurs.

CONCLUSIONS

Based on this finding, we propose that the expression onset of tissue-specific effector genes is controlled more dynamically than suggested in the previously proposed simplistic regulation scheme. Thus, we suggest that differentiation should be conceptualized as a seamless process of accumulation of effector expression along with the advancing specification GRN. This pattern of effector gene expression may have interesting implications for the evolution of novel cell types.

Soluble adenylyl cyclase coordinates intracellular pH homeostasis and biomineralization in calcifying cells of a marine animal

Biomineralizing cells concentrate dissolved inorganic carbon (DIC) and remove protons from the site of mineral precipitation. However, the molecular regulatory mechanisms that orchestrate pH homeostasis and biomineralization of calcifying cells are poorly understood. Here, we report that the acid-base sensing enzyme soluble adenylyl cyclase (sAC) coordinates intracellular pH (pH_i) regulation in the calcifying primary mesenchyme cells (PMCs) of sea urchin larvae. Single-cell transcriptomics, in situ hybridization, and immunocytochemistry elucidated the spatiotemporal expression of sAC during skeletogenesis. Live pH_i imaging of PMCs revealed that the downregulation of sAC activity with two structurally unrelated small molecules inhibited pH_i regulation of PMCs, an effect that was rescued by the addition of cell-permeable cAMP. Pharmacological sAC inhibition also significantly reduced normal spicule growth and spicule regeneration, establishing a link between PMC pH_i regulation and biomineralization. Finally, increased expression of sAC mRNA was detected during skeleton remineralization and exposure to CO₂-induced acidification. These findings suggest that transcriptional regulation of sAC is required to promote remineralization and to compensate for acidic stress. This work highlights the central role of sAC in coordinating acid-base regulation and biomineralization in calcifying cells of a marine animal.

Deep resilience: An evolutionary perspective on calcification in an age of ocean acidification

The success of today's calcifying organisms in tomorrow's oceans depends, in part, on the resilience of their skeletons to ocean acidification. To the extent this statement is true there is reason to have hope. Many marine calcifiers demonstrate resilience when exposed to environments that mimic near-term ocean acidification. The fossil record similarly suggests that resilience in skeletons has increased dramatically over geologic time. This "deep resilience" is seen in the long-term stability of skeletal chemistry, as well as a decreasing correlation between skeletal mineralogy and extinction risk over time. Such resilience over geologic timescales is often attributed to genetic canalization-the hardening of genetic pathways due to the evolution of increasingly complex regulatory systems. But paradoxically, our current knowledge on biomineralization genetics suggests an opposing trend, where genes are co-opted and shuffled at an evolutionarily rapid pace. In this paper we consider two possible mechanisms driving deep resilience in skeletons that fall outside of genetic canalization: microbial co-regulation and macroevolutionary trends in skeleton structure. The mechanisms driving deep resilience should be considered when creating risk assessments for marine organisms facing ocean acidification and provide a wealth of research avenues to explore.

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-like hormone signaling in invertebrates and its potential role in initial screening of thyroid hormone system disrupting chemicals

This review examines the presence and evolution of thyroid-like systems in selected aquatic invertebrates to determine the potential use of these organisms in screens for vertebrate thyroid hormone axis disrupting chemicals (THADCs). Such a screen might support the phasing out of some vertebrate testing. Although arthropods including crustaceans do not contain a functional thyroid signaling system, elements of such a system exist in the aquatic phyla mollusks, echinoderms, tunicates, and cephalochordates. These phyla can synthesize thyroid hormone, which has been demonstrated in some groups to induce the nuclear thyroid hormone receptor (THR). Thyroid hormone may act in these phyla through interaction with a membrane integrin receptor. Thyroid hormone regulates inter alia metamorphosis but, unlike in vertebrates, this does not occur via receptor activation by the ligands triiodothyronine (T3) and thyroxine (T4). Instead, the unliganded nuclear receptor itself controls metamorphosis in mollusks, echinoderms, and tunicates, whereas the T3 derivative tri-iodothyroacetic acid (TRIAC) acts as a THR ligand in cephalochordates. In view of this, it may be possible to develop an invertebrate-based screen that is sensitive to vertebrate THADCs that interfere with thyroid hormone synthesis or metabolism along with interaction with membrane receptors. The review makes some recommendations for the need to develop an appropriate test method. Integr Environ Assess Manag 2023;19:63-82. © 2022 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

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.

The evolution of neurosensation provides opportunities and constraints for phenotypic plasticity

Phenotypic plasticity is widely regarded as important for enabling species resilience to environmental change and for species evolution. However, insight into the complex mechanisms by which phenotypic plasticity evolves in nature is limited by our ability to reconstruct evolutionary histories of plasticity. By using part of the molecular mechanism, we were able to trace the evolution of pre-feeding phenotypic plasticity across the class Echinoidea and identify the origin of plasticity at the base of the regular urchins. The neurosensory foundation for plasticity was ancestral within the echinoids. However, coincident development of the plastic trait and the neurosensory system was not achieved until the regular urchins, likely due to pleiotropic effects and linkages between the two colocalized systems. Plasticity continues to evolve within the urchins with numerous instances of losses associated with loss of sensory abilities and neurons, consistent with a cost of maintaining these capabilities. Thus, evidence was found for the neurosensory system providing opportunities and constraints to the evolution of phenotypic plasticity.

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.

Biomineralization: Integrating mechanism and evolutionary history

Calcium carbonate (CaCO₃) biomineralizing organisms have played major roles in the history of life and the global carbon cycle during the past 541 Ma. Both marine diversification and mass extinctions reflect physiological responses to environmental changes through time. An integrated understanding of carbonate biomineralization is necessary to illuminate this evolutionary record and to understand how modern organisms will respond to 21st century global change. Biomineralization evolved independently but convergently across phyla, suggesting a unity of mechanism that transcends biological differences. In this review, we combine CaCO₃ skeleton formation mechanisms with constraints from evolutionary history, omics, and a meta-analysis of isotopic data to develop a plausible model for CaCO₃ biomineralization applicable to all phyla. The model provides a framework for understanding the environmental sensitivity of marine calcifiers, past mass extinctions, and resilience in 21st century acidifying oceans. Thus, it frames questions about the past, present, and future of CaCO₃ biomineralizing organisms.

Toxicological Impact of Rare Earth Elements (REEs) on the Reproduction and Development of Aquatic Organisms Using Sea Urchins as Biological Models

The growing presence of lanthanides in the environment has drawn the attention of the scientific community on their safety and toxicity. The sources of lanthanides in the environment include diagnostic medicine, electronic devices, permanent magnets, etc. Their exponential use and the poor management of waste disposal raise serious concerns about the quality and safety of the ecosystems at a global level. This review focused on the impact of lanthanides in marine organisms on reproductive fitness, fertilization and embryonic development, using the sea urchin as a biological model system. Scientific evidence shows that exposure to lanthanides triggers a wide variety of toxic insults, including reproductive performance, fertilization, redox metabolism, embryogenesis, and regulation of embryonic gene expression. This was thoroughly demonstrated for gadolinium, the most widely used lanthanide in diagnostic medicine, whose uptake in sea urchin embryos occurs in a time- and concentration-dependent manner, correlates with decreased calcium absorption and primarily affects skeletal growth, with incorrect regulation of the skeletal gene regulatory network. The results collected on sea urchin embryos demonstrate a variable sensitivity of the early life stages of different species, highlighting the importance of testing the effects of pollution in different species. The accumulation of lanthanides and their emerging negative effects make risk assessment and consequent legislative intervention on their disposal mandatory.

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.

Ion Pathways in Biomineralization: Perspectives on Uptake, Transport, and Deposition of Calcium, Carbonate, and Phosphate

Minerals are formed by organisms in all of the kingdoms of life. Mineral formation pathways all involve uptake of ions from the environment, transport of ions by cells, sometimes temporary storage, and ultimately deposition in or outside of the cells. Even though the details of how all this is achieved vary enormously, all pathways need to respect both the chemical limitations of ion manipulation, as well as the many "housekeeping" roles of ions in cell functioning. Here we provide a chemical perspective on the biological pathways of biomineralization. Our approach is to compare and contrast the ion pathways involving calcium, phosphate, and carbonate in three very different organisms: the enormously abundant unicellular marine coccolithophores, the well investigated sea urchin larval model for single crystal formation, and the complex pathways used by vertebrates to form their bones. The comparison highlights both common and unique processes. Significantly, phosphate is involved in regulating calcium carbonate deposition and carbonate is involved in regulating calcium phosphate deposition. One often overlooked commonality is that, from uptake to deposition, the solutions involved are usually supersaturated. This therefore requires not only avoiding mineral deposition where it is not needed but also exploiting this saturated state to produce unstable mineral precursors that can be conveniently stored, redissolved, and manipulated into diverse shapes and upon deposition transformed into more ordered and hence often functional final deposits.

Developmental single-cell transcriptomics in the Lytechinus variegatus sea urchin embryo

Using scRNA-seq coupled with computational approaches, we studied transcriptional changes in cell states of sea urchin embryos during development to the larval stage. Eighteen closely spaced time points were taken during the first 24 h of development of Lytechinus variegatus (Lv). Developmental trajectories were constructed using Waddington-OT, a computational approach to 'stitch' together developmental time points. Skeletogenic and primordial germ cell trajectories diverged early in cleavage. Ectodermal progenitors were distinct from other lineages by the 6th cleavage, although a small percentage of ectoderm cells briefly co-expressed endoderm markers that indicated an early ecto-endoderm cell state, likely in cells originating from the equatorial region of the egg. Endomesoderm cells also originated at the 6th cleavage and this state persisted for more than two cleavages, then diverged into distinct endoderm and mesoderm fates asynchronously, with some cells retaining an intermediate specification status until gastrulation. Seventy-nine out of 80 genes (99%) examined, and included in published developmental gene regulatory networks (dGRNs), are present in the Lv-scRNA-seq dataset and are expressed in the correct lineages in which the dGRN circuits operate.

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.

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.

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.

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.

Sea Urchin as a Universal Model for Studies of Gene Networks

The purple sea urchin Strongylocentrotus purpuratus has been used for over 150 years as a model organism in developmental biology. Using this model species, scientists have been able to describe, in detail, the mechanisms of cell cycle control and cell adhesion, fertilization, calcium signaling, cell differentiation, and death. Massive parallel sequencing of the sea urchin genome enabled the deciphering of the main components of gene regulatory networks during the activation of embryonic signaling pathways. This knowledge helped to extrapolate aberrations in somatic cells that may lead to diseases, including cancer in humans. Furthermore, since many, if not all, developmental signaling pathways were shown to be controlled by non-coding RNAs (ncRNAs), the sea urchin organism represents an attractive experimental model. In this review, we discuss the main discoveries in the genetics, genomics, and transcriptomics of sea urchins during embryogenesis with the main focus on the role of ncRNAs. This information may be useful for comparative studies between different organisms, and may help identify new regulatory networks controlled by ncRNAs.

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.

Ocean acidification induces distinct transcriptomic responses across life history stages of the sea urchin Heliocidaris erythrogramma

Ocean acidification (OA) from seawater uptake of rising carbon dioxide emissions impairs development in marine invertebrates, particularly in calcifying species. Plasticity in gene expression is thought to mediate many of these physiological effects, but how these responses change across life history stages remains unclear. The abbreviated lecithotrophic development of the sea urchin Heliocidaris erythrogramma provides a valuable opportunity to analyse gene expression responses across a wide range of life history stages, including the benthic, post-metamorphic juvenile. We measured the transcriptional response to OA in H. erythrogramma at three stages of the life cycle (embryo, larva, and juvenile) in a controlled breeding design. The results reveal a broad range of strikingly stage-specific impacts of OA on transcription, including changes in the number and identity of affected genes; the magnitude, sign, and variance of their expression response; and the developmental trajectory of expression. The impact of OA on transcription was notably modest in relation to gene expression changes during unperturbed development and much smaller than genetic contributions from parentage. The latter result suggests that natural populations may provide an extensive genetic reservoir of resilience to OA. Taken together, these results highlight the complexity of the molecular response to OA, its substantial life history stage specificity, and the importance of contextualizing the transcriptional response to pH stress in light of normal development and standing genetic variation to better understand the capacity for marine invertebrates to adapt to OA.

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.

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.

Genetic basis for divergence in developmental gene expression in two closely related sea urchins

The genetic basis for divergence in developmental gene expression among species is poorly understood, despite growing evidence that such changes underlie many interesting traits. Here we quantify transcription in hybrids of Heliocidaris tuberculata and Heliocidaris erythrogramma, two closely related sea urchins with highly divergent developmental gene expression and life histories. We find that most expression differences between species result from genetic influences that affect one stage of development, indicating limited pleiotropic consequences for most mutations that contribute to divergence in gene expression. Activation of zygotic transcription is broadly delayed in H. erythrogramma, the species with the derived life history, despite its overall faster premetamorphic development. Altered expression of several terminal differentiation genes associated with the derived larval morphology of H. erythrogramma is based largely on differences in the expression or function of their upstream regulators, providing insights into the genetic basis for the evolution of key life history traits.

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.

Gastrulation in the sea urchin

Gastrulation is arguably the most important evolutionary innovation in the animal kingdom. This process provides the basic embryonic architecture, an inner layer separated from an outer layer, from which all animal forms arise. An extraordinarily simple and elegant process of gastrulation is observed in the sea urchin embryo. The cells participating in sea urchin gastrulation are specified early during cleavage. One outcome of that specification is the expression of transcription factors that control each of the many subsequent morphogenetic changes. The first of these movements is an epithelial-mesenchymal transition (EMT) of skeletogenic mesenchyme cells, then EMT of pigment cell progenitors. Shortly thereafter, invagination of the archenteron occurs. At the end of archenteron extension, a second wave of EMT occurs to release immune cells into the blastocoel and primordial germ cells that will home to the coelomic pouches. The archenteron then remodels to establish the three parts of the gut, and at the anterior end, the gut fuses with the stomodaeum to form the through-gut. As part of the anterior remodeling, mesodermal coelomic pouches bud off the lateral sides of the archenteron tip. Multiple cell biological processes conduct each of these movements and in some cases the upstream transcription factors controlling this process have been identified. Remarkably, each event seamlessly occurs at the right time to orchestrate formation of the primitive body plan. This review covers progress toward understanding many of the molecular mechanisms underlying this sequence of morphogenetic events.

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.

Novel Genes, Ancient Genes, and Gene Co-Option Contributed to the Genetic Basis of the Radula, a Molluscan Innovation

The radula is the central foraging organ and apomorphy of the Mollusca. However, in contrast to other innovations, including the mollusk shell, genetic underpinnings of radula formation remain virtually unknown. Here, we present the first radula formative tissue transcriptome using the viviparous freshwater snail Tylomelania sarasinorum and compare it to foot tissue and the shell-building mantle of the same species. We combine differential expression, functional enrichment, and phylostratigraphic analyses to identify both specific and shared genetic underpinnings of the three tissues as well as their dominant functions and evolutionary origins. Gene expression of radula formative tissue is very distinct, but nevertheless more similar to mantle than to foot. Generally, the genetic bases of both radula and shell formation were shaped by novel orchestration of preexisting genes and continuous evolution of novel genes. A significantly increased proportion of radula-specific genes originated since the origin of stem-mollusks, indicating that novel genes were especially important for radula evolution. Genes with radula-specific expression in our study are frequently also expressed during the formation of other lophotrochozoan hard structures, like chaetae (hes1, arx), spicules (gbx), and shells of mollusks (gbx, heph) and brachiopods (heph), suggesting gene co-option for hard structure formation. Finally, a Lophotrochozoa-specific chitin synthase with a myosin motor domain (CS-MD), which is expressed during mollusk and brachiopod shell formation, had radula-specific expression in our study. CS-MD potentially facilitated the construction of complex chitinous structures and points at the potential of molecular novelties to promote the evolution of different morphological innovations.

How Does the Regulatory Genome Work?

The regulatory genome controls genome activity throughout the life of an organism. This requires that complex information processing functions are encoded in, and operated by, the regulatory genome. Although much remains to be learned about how the regulatory genome works, we here discuss two cases where regulatory functions have been experimentally dissected in great detail and at the systems level, and formalized by computational logic models. Both examples derive from the sea urchin embryo, but assess two distinct organizational levels of genomic information processing. The first example shows how the regulatory system of a single gene, endo16, executes logic operations through individual transcription factor binding sites and cis-regulatory modules that control the expression of this gene. The second example shows information processing at the gene regulatory network (GRN) level. The GRN controlling development of the sea urchin endomesoderm has been experimentally explored at an almost complete level. A Boolean logic model of this GRN suggests that the modular logic functions encoded at the single-gene level show compositionality and suffice to account for integrated function at the network level. We discuss these examples both from a biological-experimental point of view and from a computer science-informational point of view, as both illuminate principles of how the regulatory genome works.

Inference of Developmental Gene Regulatory Networks Beyond Classical Model Systems: New Approaches in the Post-genomic Era

The advent of high-throughput sequencing (HTS) technologies has revolutionized the way we understand the transformation of genetic information into morphological traits. Elucidating the network of interactions between genes that govern cell differentiation through development is one of the core challenges in genome research. These networks are known as developmental gene regulatory networks (dGRNs) and consist largely of the functional linkage between developmental control genes, cis-regulatory modules, and differentiation genes, which generate spatially and temporally refined patterns of gene expression. Over the last 20 years, great advances have been made in determining these gene interactions mainly in classical model systems, including human, mouse, sea urchin, fruit fly, and worm. This has brought about a radical transformation in the fields of developmental biology and evolutionary biology, allowing the generation of high-resolution gene regulatory maps to analyze cell differentiation during animal development. Such maps have enabled the identification of gene regulatory circuits and have led to the development of network inference methods that can recapitulate the differentiation of specific cell-types or developmental stages. In contrast, dGRN research in non-classical model systems has been limited to the identification of developmental control genes via the candidate gene approach and the characterization of their spatiotemporal expression patterns, as well as to the discovery of cis-regulatory modules via patterns of sequence conservation and/or predicted transcription-factor binding sites. However, thanks to the continuous advances in HTS technologies, this scenario is rapidly changing. Here, we give a historical overview on the architecture and elucidation of the dGRNs. Subsequently, we summarize the approaches available to unravel these regulatory networks, highlighting the vast range of possibilities of integrating multiple technical advances and theoretical approaches to expand our understanding on the global gene regulation during animal development in non-classical model systems. Such new knowledge will not only lead to greater insights into the evolution of molecular mechanisms underlying cell identity and animal body plans, but also into the evolution of morphological key innovations in animals.

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.

Expression of exogenous mRNAs to study gene function in echinoderm embryos

With the completion of the genome sequencing projects, a new challenge for developmental biologists is to assign a function to the thousands of genes identified. Expression of exogenous mRNAs is a powerful, versatile and rapid technique that can be used to study gene function during development of the sea urchin. This chapter describes how this technique can be used to analyze gene function in echinoderm embryos, how it can be combined with cell transplantation to perform mosaic analysis and how it can be applied to identify downstream targets genes of transcription factors and signaling pathways. We describe specific examples of the use of overexpression of mRNA to analyze gene function, mention the benefits and current limitations of the technique and emphasize the importance of using different controls to assess the specificity of the effects observed. Finally, this chapter details the different steps, vectors and protocols for in vitro production of mRNA and phenotypic analysis.

Genome-wide analysis of chromatin accessibility using ATAC-seq

Programs of gene transcription are controlled by cis-acting DNA elements, including enhancers, silencers, and promoters. Local accessibility of chromatin has proven to be a highly informative structural feature for identifying such regulatory elements, which tend to be relatively open due to their interactions with proteins. Recently, ATAC-seq (assay for transposase-accessible chromatin using sequencing) has emerged as one of the most powerful approaches for genome-wide chromatin accessibility profiling. This method assesses DNA accessibility using hyperactive Tn5 transposase, which simultaneously cuts DNA and inserts sequencing adaptors, preferentially in regions of open chromatin. ATAC-seq is a relatively simple procedure which can be applied to only a few thousand cells. It is well-suited to developing embryos of sea urchins and other echinoderms, which are a prominent experimental model for understanding the genomic control of animal development. In this chapter, we present a protocol for applying ATAC-seq to embryonic cells of sea urchins.

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.

Regulatory states in the developmental control of gene expression

A growing body of evidence shows that gene expression in multicellular organisms is controlled by the combinatorial function of multiple transcription factors. This indicates that not the individual transcription factors or signaling molecules, but the combination of expressed regulatory molecules, the regulatory state, should be viewed as the functional unit in gene regulation. Here, I discuss the concept of the regulatory state and its proposed role in the genome-wide control of gene expression. Recent analyses of regulatory gene expression in sea urchin embryos have been instrumental for solving the genomic control of cell fate specification in this system. Some of the approaches that were used to determine the expression of regulatory states during sea urchin embryogenesis are reviewed. Significant developmental changes in regulatory state expression leading to the distinct specification of cell fates are regulated by gene regulatory network circuits. How these regulatory state transitions are encoded in the genome is illuminated using the sea urchin endoderm-mesoderms cell fate decision circuit as an example. These observations highlight the importance of considering developmental gene regulation, and the function of individual transcription factors, in the context of regulatory states.

Assessing regulatory information in developmental gene regulatory networks

Gene regulatory networks (GRNs) provide a transformation function between the static genomic sequence and the primary spatial specification processes operating development. The regulatory information encompassed in developmental GRNs thus goes far beyond the control of individual genes. We here address regulatory information at different levels of network organization, from single node to subcircuit to large-scale GRNs and discuss how regulatory design features such as network architecture, hierarchical organization, and cis-regulatory logic contribute to the developmental function of network circuits. Using specific subcircuits from the sea urchin endomesoderm GRN, for which both circuit design and biological function have been described, we evaluate by Boolean modeling and in silico perturbations the import of given circuit features on developmental function. The examples include subcircuits encoding positive feedback, mutual repression, and coherent feedforward, as well as signaling interaction circuitry. Within the hierarchy of the endomesoderm GRN, these subcircuits are organized in an intertwined and overlapping manner. Thus, we begin to see how regulatory information encoded at individual nodes is integrated at all levels of network organization to control developmental process.