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

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.

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.

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.

Transcriptomic Responses of a Lightly Calcified Echinoderm to Experimental Seawater Acidification and Warming during Early Development

Ocean acidification (OA) and ocean warming (OW) are potential obstacles to the survival and growth of marine organisms, particularly those that rely on calcification. This study investigated the single and joint effects of OA and OW on sea cucumber Apostichopus japonicus larvae raised under combinations of two temperatures (19 °C or 22 °C) and two pCO2 levels (400 or 1000 μatm) that reflect the current and end-of-21st-century projected ocean scenarios. The investigation focused on assessing larval development and identifying differences in gene expression patterns at four crucial embryo-larval stages (blastula, gastrula, auricularia, and doliolaria) of sea cucumbers, using RNA-seq. Results showed the detrimental effect of OA on the early development and body growth of A. japonicus larvae and a reduction in the expression of genes associated with biomineralization, skeletogenesis, and ion homeostasis. This effect was particularly pronounced during the doliolaria stage, indicating the presence of bottlenecks in larval development at this transition phase between the larval and megalopa stages in response to OA. OW accelerated the larval development across four stages of A. japonicus, especially at the blastula and doliolaria stages, but resulted in a widespread upregulation of genes related to heat shock proteins, antioxidant defense, and immune response. Significantly, the negative effects of elevated pCO2 on the developmental process of larvae appeared to be mitigated when accompanied by increased temperatures at the expense of reduced immune resilience and increased system fragility. These findings suggest that alterations in gene expression within the larvae of A. japonicus provide a mechanism to adapt to stressors arising from a rapidly changing oceanic environment.

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

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

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

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

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.

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.

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.

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.

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.

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.

Extreme phenotypic divergence and the evolution of development

As analyses of developmental mechanisms extend to ever more species, it becomes important to understand not just what is conserved or altered during evolution, but why. Closely related species that exhibit extreme phenotypic divergence can be uniquely informative in this regard. A case in point is the sea urchin genus Heliocidaris, which contains species that recently evolved a life history involving nonfeeding larvae following nearly half a billion years of prior evolution with feeding larvae. The resulting shift in selective regimes produced rapid and surprisingly extensive changes in developmental mechanisms that are otherwise highly conserved among echinoderm species. The magnitude and extent of these changes challenges the notion that conservation of early development in echinoderms is largely due to internal constraints that prohibit modification and instead suggests that natural selection actively maintains stability of inherently malleable trait developmental mechanisms over immense time periods. Knowing how and why natural selection changed during the evolution of nonfeeding larvae can also reveal why developmental mechanisms do and do not change in particular ways.

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.

Pigment cells: Paragons of cellular development

Larvae of sea urchins have a population of conspicuous pigmented cells embedded in the outer surface epithelium. Pigment cells are a distinct mesodermal lineage that gives rise to a key component of the larval immune system. During cleavage, signaling from adjacent cells influences a small crescent of cells to initiate a network of genetic interactions that prepare the cells for morphogenesis and specializes them as immunocytes. The cells become active during gastrulation, detach from the epithelium, migrate through the blastocoel, and insert into the ectoderm where they complete their differentiation. Studies of pigment cell development have helped establish how cellular signaling controls networks of genetic interactions that bring about morphogenesis and differentiation. This review summarizes studies of pigment cell development and concludes that pigment cells are an excellent experimental model. Pigment cells provide several opportunities to further test and refine our understanding of the molecular basis of cellular development.

Perspectives on divergence of early developmental regulatory pathways: Insight from the evolution of echinoderm double negative gate

Evolution of gene regulatory networks (GRN) that orchestrate the highly coordinated course of development, is made possible by the network's robust nature for incorporating change without detrimental developmental outcome. It can be considered that the upstream network regulating early development, has immense influence over succeeding pathways thus may be less subjected to evolutionary modification. However, recent studies show incorporation of novel genes in such early developmental pathways such as the echinoderm pmar1 as evidence for drastic change occurring high in the GRN hierarchy. Here we discuss the mechanisms that underlie divergence of early developmental pathways utilizing promising insights from the evolution of echinoderm early mesoderm specification pathway of Pmar1-HesC double negative gate found solely in the euechinoid sea urchin lineage, as well as examples from other groups such as Spiralia and Drosophila.

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.

Gene regulation of adult skeletogenesis in starfish and modifications during gene network co-option

The larval skeleton of the echinoderm is believed to have been acquired through co-option of a pre-existing gene regulatory network (GRN); that is, the mechanism for adult skeleton formation in the echinoderm was deployed in early embryogenesis during echinoderm diversification. To explore the evolutionary changes that occurred during co-option, we examined the mechanism for adult skeletogenesis using the starfish Patiria pectinifera. Expression patterns of skeletogenesis-related genes (vegf, vegfr, ets1/2, erg, alx1, ca1, and clect) suggest that adult skeletogenic cells develop from the posterior coelom after the start of feeding. Treatment with inhibitors and gene knockout using transcription activator-like effector nucleases (TALENs) suggest that the feeding-nutrient sensing pathway activates Vegf signaling via target of rapamycin (TOR) activity, leading to the activation of skeletogenic regulatory genes in starfish. In the larval skeletogenesis of sea urchins, the homeobox gene pmar1 activates skeletogenic regulatory genes, but in starfish, localized expression of the pmar1-related genes phbA and phbB was not detected during the adult skeleton formation stage. Based on these data, we provide a model for the adult skeletogenic GRN in the echinoderm and propose that the upstream regulatory system changed from the feeding-TOR-Vegf pathway to a homeobox gene-system during co-option of the skeletogenic GRN.

A nomenclature for echinoderm genes

Echinoderm embryos and larvae are prominent experimental model systems for studying developmental mechanisms. High-quality, assembled, annotated genome sequences are now available for several echinoderm species, including representatives from most classes. The increased availability of these data necessitates the development of a nomenclature that assigns universally interpretable gene symbols to echinoderm genes to facilitate cross-species comparisons of gene functions, both within echinoderms and across other phyla. This paper describes the implementation of an improved set of echinoderm gene nomenclature guidelines that both communicates meaningful orthology information in protein-coding gene symbols and names and establishes continuity with nomenclatures developed for major vertebrate model organisms, including humans. Differences between the echinoderm gene nomenclature guidelines and vertebrate guidelines are examined and explained. This nomenclature incorporates novel solutions to allow for several types of orthologous relationships, including the single echinoderm genes with multiple vertebrate co-orthologs that result from whole-genome-duplication events. The current version of the Echinoderm Gene Nomenclature Guidelines can be found at https://www.echinobase.org/gene/static/geneNomenclature.jsp

Database URL https://www.echinobase.org/

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.

Conditional gene knockdowns in sea urchins using caged morpholinos

Echinoderms are important experimental models for analyzing embryonic development, but a lack of spatial and temporal control over gene perturbations has hindered developmental studies using these animals. Morpholino antisense oligonucleotides (MOs) have been used successfully by the echinoderm research community for almost two decades, and MOs remain the most widely used tool for acute gene knockdowns in these organisms. Echinoderm embryos develop externally and are optically transparent, making them ideally-suited to many light-based approaches for analyzing and manipulating development. Studies using zebrafish embryos have demonstrated the effectiveness of photoactivatable (caged) MOs for conditional gene knockdowns. Here we show that caged MOs, synthesized using nucleobase-caged monomers, provide light-regulated control over gene expression in sea urchin embryos. Our work provides the first robust approach for conditional gene silencing in this prominent model system.

Coup-TF: A maternal factor essential for differentiation along the embryonic axes in the sea urchin Paracentrotus lividus

Coup-TF, a member of the nuclear receptor super-family, is present in the pool of maternal mRNAs and proteins in the sea urchin egg. The presence of this protein seems to be essential for the execution of the early developmental program, leading to all three embryonic layers. Our results demonstrate that Pl-Coup-TF morphants, i.e. Pl-Coup-TF morpholino knockdown embryos, resemble blastulae that lack archenteron at 24 hpf (hours post fertilization), a stage at which normal embryos reach the end of gastrulation in Paracentrotus lividus. At 48 hpf, when normal embryos reach the pluteus larva stage, the morphants are seemingly underdeveloped and lack the characteristic skeletal rods. Nevertheless, the morphant embryos express vegetal endomesodermal marker genes, such as Pl-Blimp1, Pl-Endo16, Pl-Alx1 and Pl-Tbr as judged by in situ hybridization experiments. The anterior neuroectoderm genes, Pl-FoxQ2, Pl-Six3 and Pl-Pax6, are also expressed in the morphant embryos, but Pl-Hbn and Pl-Fez mRNAs, which encode proteins significant for the differentiation of serotonergic neurons, are not detected. Consequently, Pl-Coup-TF morphants at 48 hpf lack serotonergic neurons, whereas normal 48 hpf plutei exhibit the formation of two bilateral pairs of such neurons in the apical organ. Furthermore, genes indicative of the ciliary band formation, Pl-Hnf6, Pl-Dri, Pl-FoxG and Pl-Otx, are not expressed in Pl-Coup-TF morphants, suggesting the disruption of this neurogenic territory as well. In addition, the Pl-SynB gene, a marker of differentiated neurons, is silent leading to the hypothesis that Pl-Coup-TF morphants might lack all types of neurons. On the contrary, the genes expressing signaling molecules, which establish the ventral/dorsal axis, Pl-Nodal and Pl-Lefty show the characteristic ventral lateral expression pattern, Pl-Bmp2/4, which activates the dorsal ectoderm GRN is down-regulated and Pl-Chordin is aberrantly over-expressed in the entire ectoderm. The identity of ectodermal cells in Pl-Coup-TF morphant embryos, was probed for expression of the ventral marker Pl-Gsc which was over-expressed and dorsal markers, Pl-IrxA and Pl-Hox7, which were silent. Therefore, we propose that maternal Pl-Coup-TF is essential for correct dissemination of the early embryonic signaling along both animal/vegetal and ventral/dorsal axes. Limiting Pl-Coup-TF's quantity, results in an embryo without digestive and nervous systems, skeleton and ciliary band that cannot survive past the initial 48 h of development.

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.

Multiple origins of feeding head larvae by the Early Cambrian

In many animals the head develops early, most of the body axis later. A larva composed mostly of the developing front end therefore can attain mobility and feeding earlier in development. Fossils, functional morphology, and inferred homologies indicate that feeding head larvae existed by the Early Cambrian in members of three major clades of animals: ecdysozoans, lophotrochozoans, and deuterostomes. Some of these early larval feeding mechanisms were also those of juveniles and adults (the lophophore of brachiopod larvae and possibly the ciliary band of the dipleurula of hemichordates and echinoderms); some were derived from structures that previously had other functions (appendages of the nauplius). Trochophores that swim with a preoral band of cilia, the prototroch, originated before divergence of annelids and molluscs, but evidence of larval growth and thus a prototrochal role in feeding is lacking for molluscs until the Ordovician. Feeding larvae that definitely originated much later, as in insects, teleost fish, and amphibians, develop all or nearly all of what will become the adult body axis before they begin feeding. On present evidence, head larvae, including feeding head larvae, evolved multiple times early in the evolution of bilaterian animals and never since.

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.

Optimization of skeletal protein preparation for LC-MS/MS sequencing yields additional coral skeletal proteins in Stylophora pistillata

Stony corals generate their calcium carbonate exoskeleton in a highly controlled biomineralization process mediated by a variety of macromolecules including proteins. Fully identifying and classifying these proteins is crucial to understanding their role in exoskeleton formation, yet no optimal method to purify and characterize the full suite of extracted coral skeletal proteins has been established and hence their complete composition remains obscure. Here, we tested four skeletal protein purification protocols using acetone precipitation and ultrafiltration dialysis filters to present a comprehensive scleractinian coral skeletal proteome. We identified a total of 60 proteins in the coral skeleton, 44 of which were not present in previously published stony coral skeletal proteomes. Extracted protein purification protocols carried out in this study revealed that no one method captures all proteins and each protocol revealed a unique set of method-exclusive proteins. To better understand the general mechanism of skeletal protein transportation, we further examined the proteins’ gene ontology, transmembrane domains, and signal peptides. We found that transmembrane domain proteins and signal peptide secretion pathways, by themselves, could not explain the transportation of proteins to the skeleton. We therefore propose that some proteins are transported to the skeleton via non-traditional secretion pathways.

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.

Simulations of sea urchin early development delineate the role of oriented cell division in the morula-to-blastula transition

The sea urchin morula to blastula transition has long been thought to require oriented cell divisions and blastomere adherence to the enveloping hyaline layer. In a computer simulation model, cell divisions constrained by a surface plane division rule are adequate to effect morphological transition. The hyaline membrane acts as an enhancer but is not essential. The model is consistent with the orientation of micromere divisions and the open blastulae of direct developing species. The surface plane division rule precedes overt epithelization of surface cells and acts to organize the developing epithelium. It is a universal feature of early metazoan development and simulations of non-echinoid cleavage patterns support its role throughout Metazoa. The surface plane division rule requires only local cues and cells need not reference global positional information or embryonic axes.

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.

pmar1/phb homeobox genes and the evolution of the double-negative gate for endomesoderm specification in echinoderms

In several model animals, the earliest phases of embryogenesis are regulated by lineage-specific genes, such as Drosophila bicoid Sea urchin (echinoid) embryogenesis is initiated by zygotic expression of pmar1, a paired-class homeobox gene that has been considered to be present only in the lineage of modern urchins (euechinoids). In euechinoids, Pmar1 promotes endomesoderm specification by repressing the hairy and enhancer of split C (hesC) gene. Here, we have identified the basal echinoid (cidaroid) pmar1 gene, which also promotes endomesoderm specification but not by repressing hesC A further search for related genes demonstrated that other echinoderms have pmar1-related genes named phb Functional analyses of starfish Phb proteins indicated that, similar to cidaroid Pmar1, they promote activation of endomesoderm regulatory gene orthologs via an unknown repressor that is not HesC. Based on these results, we propose that Pmar1 may have recapitulated the regulatory function of Phb during the early diversification of echinoids and that the additional repressor HesC was placed under the control of Pmar1 in the euechinoid lineage. This case provides an exceptional model for understanding how early developmental processes diverge.

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.

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

Alx1 is a conserved regulator of skeletogenesis in echinoderms and evolutionary changes in Alx1 sequence and expression have played a pivotal role in modifying programs of skeletogenesis within the phylum. Alx1 regulates a large suite of effector genes that control the morphogenetic behaviors and biomineral-forming activities of skeletogenic cells. To better understand the gene regulatory control of skeletogenesis by Alx1, we used genome-wide ChIP-seq to identify Alx1-binding sites and direct gene targets. Our analysis revealed that many terminal differentiation genes receive direct transcriptional inputs from Alx1. In addition, we found that intermediate transcription factors previously shown to be downstream of Alx1 all receive direct inputs from Alx1. Thus, Alx1 appears to regulate effector genes by indirect, as well as direct, mechanisms. We tested 23 high-confidence ChIP-seq peaks using GFP reporters and identified 18 active cis-regulatory modules (CRMs); this represents a high success rate for CRM discovery. Detailed analysis of a representative CRM confirmed that a conserved, palindromic Alx1-binding site was essential for expression. Our work significantly advances our understanding of the gene regulatory circuitry that controls skeletogenesis in sea urchins and provides a framework for evolutionary studies.

Cell type phylogenetics informs the evolutionary origin of echinoderm larval skeletogenic cell identity

The multiplicity of cell types comprising multicellular organisms begs the question as to how cell type identities evolve over time. Cell type phylogenetics informs this question by comparing gene expression of homologous cell types in distantly related taxa. We employ this approach to inform the identity of larval skeletogenic cells of echinoderms, a clade for which there are phylogenetically diverse datasets of spatial gene expression patterns. We determined ancestral spatial expression patterns of alx1, ets1, tbr, erg, and vegfr, key components of the skeletogenic gene regulatory network driving identity of the larval skeletogenic cell. Here we show ancestral state reconstructions of spatial gene expression of extant eleutherozoan echinoderms support homology and common ancestry of echinoderm larval skeletogenic cells. We propose larval skeletogenic cells arose in the stem lineage of eleutherozoans during a cell type duplication event that heterochronically activated adult skeletogenic cells in a topographically distinct tissue in early development.

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.