A genome-wide analysis of biomineralization-related proteins in the sea urchin Strongylocentrotus purpuratus

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.

The brittle star genome illuminates the genetic basis of animal appendage regeneration

Species within nearly all extant animal lineages are capable of regenerating body parts. However, it remains unclear whether the gene expression programme controlling regeneration is evolutionarily conserved. Brittle stars are a species-rich class of echinoderms with outstanding regenerative abilities, but investigations into the genetic bases of regeneration in this group have been hindered by the limited genomic resources. Here we report a chromosome-scale genome assembly for the brittle star Amphiura filiformis. We show that the brittle star genome is the most rearranged among echinoderms sequenced so far, featuring a reorganized Hox cluster reminiscent of the rearrangements observed in sea urchins. In addition, we performed an extensive profiling of gene expression during brittle star adult arm regeneration and identified sequential waves of gene expression governing wound healing, proliferation and differentiation. We conducted comparative transcriptomic analyses with other invertebrate and vertebrate models for appendage regeneration and uncovered hundreds of genes with conserved expression dynamics, particularly during the proliferative phase of regeneration. Our findings emphasize the crucial importance of echinoderms to detect long-range expression conservation between vertebrates and classical invertebrate regeneration model systems.

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.

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.

Spatial-temporal expression analysis of lineage-restricted shell matrix proteins reveals shell field regionalization and distinct cell populations in the slipper snail Crepidula atrasolea

Molluscs are one of the most morphologically diverse clades of metazoans, exhibiting an immense diversification of calcium carbonate structures, such as the shell. Biomineralization of the calcified shell is dependent on shell matrix proteins (SMPs). While SMP diversity is hypothesized to drive molluscan shell diversity, we are just starting to unravel SMP evolutionary history and biology. Here we leveraged two complementary model mollusc systems, Crepidula fornicata and Crepidula atrasolea , to determine the lineage-specificity of 185 Crepidula SMPs. We found that 95% of the adult C. fornicata shell proteome belongs to conserved metazoan and molluscan orthogroups, with molluscan-restricted orthogroups containing half of all SMPs in the shell proteome. The low number of C. fornicata -restricted SMPs contradicts the generally-held notion that an animal’s biomineralization toolkit is dominated by mostly novel genes. Next, we selected a subset of lineage-restricted SMPs for spatial-temporal analysis using in situ hybridization chain reaction (HCR) during larval stages in C. atrasolea . We found that 12 out of 18 SMPs analyzed are expressed in the shell field. Notably, these genes are present in 5 expression patterns, which define at least three distinct cell populations within the shell field. These results represent the most comprehensive analysis of gastropod SMP evolutionary age and shell field expression patterns to date. Collectively, these data lay the foundation for future work to interrogate the molecular mechanisms and cell fate decisions underlying molluscan mantle specification and diversification.

Worms and gills, plates and spines: the evolutionary origins and incredible disparity of deuterostomes revealed by fossils, genes, and development

Deuterostomes are the major division of animal life which includes sea stars, acorn worms, and humans, among a wide variety of ecologically and morphologically disparate taxa. However, their early evolution is poorly understood, due in part to their disparity, which makes identifying commonalities difficult, as well as their relatively poor early fossil record. Here, we review the available morphological, palaeontological, developmental, and molecular data to establish a framework for exploring the origins of this important and enigmatic group. Recent fossil discoveries strongly support a vermiform ancestor to the group Hemichordata, and a fusiform active swimmer as ancestor to Chordata. The diverse and anatomically bewildering variety of forms among the early echinoderms show evidence of both bilateral and radial symmetry. We consider four characteristics most critical for understanding the form and function of the last common ancestor to Deuterostomia: Hox gene expression patterns, larval morphology, the capacity for biomineralization, and the morphology of the pharyngeal region. We posit a deuterostome last common ancestor with a similar antero-posterior gene regulatory system to that found in modern acorn worms and cephalochordates, a simple planktonic larval form, which was later elaborated in the ambulacrarian lineage, the ability to secrete calcium minerals in a limited fashion, and a pharyngeal respiratory region composed of simple pores. This animal was likely to be motile in adult form, as opposed to the sessile origins that have been historically suggested. Recent debates regarding deuterostome monophyly as well as the wide array of deuterostome-affiliated problematica further suggest the possibility that those features were not only present in the last common ancestor of Deuterostomia, but potentially in the ur-bilaterian. The morphology and development of the early deuterostomes, therefore, underpin some of the most significant questions in the study of metazoan evolution.

EchinoDB: an update to the web-based application for genomic and transcriptomic data on echinoderms

BACKGROUND

Here we release a new version of EchinoDB, EchinoDB v2.0 ( https://echinodb.uncc.edu ). EchinoDB is a database of genomic and transcriptomic data on echinoderms. The initial database consisted of groups of 749,397 orthologous and paralogous transcripts arranged in orthoclusters by sequence similarity.

RESULTS

The updated version of EchinoDB includes two new major datasets: the RNA-Seq data of the brittle star Ophioderma brevispinum and the high-quality genomic assembly data of the green sea urchin Lytechinus variegatus. In addition, we enabled keyword searches for annotated data and installed an updated version of Sequenceserver to allow Basic Local Alignment Search Tool (BLAST) searches. The data are downloadable in FASTA format. The first version of EchinoDB appeared in 2016 and was implemented in GO on a local server. The new version has been updated using R Shiny to include new features and improvements in the application. Furthermore, EchinoDB now runs entirely in the cloud for increased reliability and scaling.

CONCLUSION

EchinoDB serves a user base drawn from the fields of phylogenetics, developmental biology, genomics, physiology, neurobiology, and regeneration. As use cases, we illustrate the function of EchinoDB in retrieving components of signaling pathways involved in the tissue regeneration process of different echinoderms, including the emerging model species Ophioderma brevispinum. Moreover, we use EchinoDB to shed light on the conservation of the molecular components involved in two echinoderm-specific phenomena: spicule matrix proteins involved in the formation of stereom endoskeleton and the tensilin protein that contributes to the capacity of the connective tissues to quickly change its mechanical properties. The genes involved in the former had been previously studied in echinoids, while gene sequences involved in the latter had been previously described in holothuroids. Specifically, we ask (a) if the biomineralization-related proteins previously reported only in sea urchins are also present in other, non-echinoid, echinoderms and (b) if tensilin, the protein responsible for the control of stiffness of the mutable collagenous tissue, previously described in sea cucumbers, is conserved across the phylum.

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.

Brachiopod and mollusc biomineralisation is a conserved process that was lost in the phoronid-bryozoan stem lineage

BACKGROUND

Brachiopods and molluscs are lophotrochozoans with hard external shells which are often believed to have evolved convergently. While palaeontological data indicate that both groups are descended from biomineralising Cambrian ancestors, the closest relatives of brachiopods, phoronids and bryozoans, are mineralised to a much lower extent and are comparatively poorly represented in the Palaeozoic fossil record. Although brachiopod and mollusc shells are structurally analogous, genomic and proteomic evidence indicates that their formation involves a complement of conserved, orthologous genes. Here, we study a set of genes comprised of 3 homeodomain transcription factors, one signalling molecule and 6 structural proteins which are implicated in mollusc and brachiopod shell formation, search for their orthologs in transcriptomes or genomes of brachiopods, phoronids and bryozoans, and present expression patterns of 8 of the genes in postmetamorphic juveniles of the rhynchonelliform brachiopod T. transversa.

RESULTS

Transcriptome and genome searches for the 10 target genes in the brachiopods Terebratalia transversa, Lingula anatina, Novocrania anomala, the bryozoans Bugula neritina and Membranipora membranacea, and the phoronids Phoronis australis and Phoronopsis harmeri resulted in the recovery of orthologs of the majority of the genes in all taxa. While the full complement of genes was present in all brachiopods with a single exception in L. anatina, a bloc of four genes could consistently not be retrieved from bryozoans and phoronids. The genes engrailed, distal-less, ferritin, perlucin, sp1 and sp2 were shown to be expressed in the biomineralising mantle margin of T. transversa juveniles.

CONCLUSIONS

The gene expression patterns we recovered indicate that while mineralised shells in brachiopods and molluscs are structurally analogous, their formation builds on a homologous process that involves a conserved complement of orthologous genes. Losses of some of the genes related to biomineralisation in bryozoans and phoronids indicate that loss of the capacity to form mineralised structures occurred already in the phoronid-bryozoan stem group and supports the idea that mineralised skeletons evolved secondarily in some of the bryozoan subclades.

Investigating calcification-related candidates in a non-symbiotic scleractinian coral, Tubastraea spp

In hermatypic scleractinian corals, photosynthetic fixation of CO2 and the production of CaCO3 are intimately linked due to their symbiotic relationship with dinoflagellates of the Symbiodiniaceae family. This makes it difficult to study ion transport mechanisms involved in the different pathways. In contrast, most ahermatypic scleractinian corals do not share this symbiotic relationship and thus offer an advantage when studying the ion transport mechanisms involved in the calcification process. Despite this advantage, non-symbiotic scleractinian corals have been systematically neglected in calcification studies, resulting in a lack of data especially at the molecular level. Here, we combined a tissue micro-dissection technique and RNA-sequencing to identify calcification-related ion transporters, and other candidates, in the ahermatypic non-symbiotic scleractinian coral Tubastraea spp. Our results show that Tubastraea spp. possesses several calcification-related candidates previously identified in symbiotic scleractinian corals (such as SLC4-γ, AMT-1like, CARP, etc.). Furthermore, we identify and describe a role in scleractinian calcification for several ion transporter candidates (such as SLC13, -16, -23, etc.) identified for the first time in this study. Taken together, our results provide not only insights about the molecular mechanisms underlying non-symbiotic scleractinian calcification, but also valuable tools for the development of biotechnological solutions to better control the extreme invasiveness of corals belonging to this particular genus.

Comparative Single-Cell Transcriptomics Reveals Novel Genes Involved in Bivalve Embryonic Shell Formation and Questions Ontogenetic Homology of Molluscan Shell Types

Mollusks are known for their highly diverse repertoire of body plans that often includes external armor in form of mineralized hardparts. Representatives of the Conchifera, one of the two major lineages that comprises taxa which originated from a uni-shelled ancestor (Monoplacophora, Gastropoda, Cephalopoda, Scaphopoda, Bivalvia), are particularly relevant regarding the evolution of mollusk shells. Previous studies have found that the shell matrix of the adult shell (teleoconch) is rapidly evolving and that the gene set involved in shell formation is highly taxon-specific. However, detailed annotation of genes expressed in tissues involved in the formation of the embryonic shell (protoconch I) or the larval shell (protoconch II) are currently lacking. Here, we analyzed the genetic toolbox involved in embryonic and larval shell formation in the quagga mussel Dreissena rostriformis using single cell RNA sequencing. We found significant differences in genes expressed during embryonic and larval shell secretion, calling into question ontogenetic homology of these transitory bivalve shell types. Further ortholog comparisons throughout Metazoa indicates that a common genetic biomineralization toolbox, that was secondarily co-opted into molluscan shell formation, was already present in the last common metazoan ancestor. Genes included are engrailed, carbonic anhydrase, and tyrosinase homologs. However, we found that 25% of the genes expressed in the embryonic shell field of D. rostriformis lack an ortholog match with any other metazoan. This indicates that not only adult but also embryonic mollusk shells may be fast-evolving structures. We raise the question as to what degree, and on which taxonomic level, the gene complement involved in conchiferan protoconch formation may be lineage-specific or conserved across taxa.

Toxicity of Vanadium during Development of Sea Urchin Embryos: Bioaccumulation, Calcium Depletion, ERK Modulation and Cell-Selective Apoptosis

Vanadium toxicology is a topic of considerable importance as this metal is widely used in industrial and biomedical fields. However, it represents a potential emerging environmental pollutant because wastewater treatment plants do not adequately remove metal compounds that are subsequently released into the environment. Vanadium applications are limited due to its toxicity, so it is urgent to define this aspect. This metal is associated with sea urchin embryo toxicity as it perturbs embryogenesis and skeletogenesis, triggering several stress responses. Here we investigated its bioaccumulation and the correlation with cellular and molecular developmental pathways. We used cytotoxic concentrations of 1 mM and 500 μM to perform quantitative analyses, showing that vanadium accumulation interferes with calcium uptake during sea urchin development and provokes a disruption in the biomineralization process. At the end of the whole treatment, the accumulation of vanadium was about 14 and 8 μg for embryos treated respectively with 1 mM and 500 μM, showing a dose-dependent response. Then, we monitored the cell signaling perturbation, analyzing key molecular markers of cell survival/cell death mechanisms and the DNA fragmentation associated with apoptosis. This paper clarifies vanadium’s trend to accumulate directly into embryonic cells, interfering with calcium uptake. In addition, our results indicate that vanadium can modulate the ERK pathway and activate a cell-selective apoptosis. These results endorse the sea urchin embryo as an adequate experimental model to study metal-related cellular/molecular responses.

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.

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.

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.

Enzymatic Approach in Calcium Phosphate Biomineralization: A Contribution to Reconcile the Physicochemical with the Physiological View

Biomineralization is the process by which organisms produce hard inorganic matter from soft tissues with outstanding control of mineral deposition in time and space. For this purpose, organisms deploy a sophisticated "toolkit" that has resulted in significant evolutionary innovations, for which calcium phosphate (CaP) is the biomineral selected for the skeleton of vertebrates. While CaP mineral formation in aqueous media can be investigated by studying thermodynamics and kinetics of phase transitions in supersaturated solutions, biogenic mineralization requires coping with the inherent complexity of biological systems. This mainly includes compartmentalization and homeostatic processes used by organisms to regulate key physiological factors, including temperature, pH and ion concentration. A detailed analysis of the literature shows the emergence of two main views describing the mechanism of CaP biomineralization. The first one, more dedicated to the study of in vivo systems and supported by researchers in physiology, often involves matrix vesicles (MVs). The second one, more investigated by the physicochemistry community, involves collagen intrafibrillar mineralization particularly through in vitro acellular models. Herein, we show that there is an obvious need in the biological systems to control both where and when the mineral forms through an in-depth survey of the mechanism of CaP mineralization. This necessity could gather both communities of physiologists and physicochemists under a common interest for an enzymatic approach to better describe CaP biomineralization. Both homogeneous and heterogeneous enzymatic catalyses are conceivable for these systems, and a few preliminary promising results on CaP mineralization for both types of enzymatic catalysis are reported in this work. Through them, we aim to describe the relevance of our point of view and the likely findings that could be obtained when adding an enzymatic approach to the already rich and creative research field dealing with CaP mineralization. This complementary approach could lead to a better understanding of the biomineralization mechanism and inspire the biomimetic design of new materials.

Forced Biomineralization: A Review

Biologically induced and controlled mineralization of metals promotes the development of protective structures to shield cells from thermal, chemical, and ultraviolet stresses. Metal biomineralization is widely considered to have been relevant for the survival of life in the environmental conditions of ancient terrestrial oceans. Similar behavior is seen among extremophilic biomineralizers today, which have evolved to inhabit a variety of industrial aqueous environments with elevated metal concentrations. As an example of extreme biomineralization, we introduce the category of "forced biomineralization", which we use to refer to the biologically mediated sequestration of dissolved metals and metalloids into minerals. We discuss forced mineralization as it is known to be carried out by a variety of organisms, including polyextremophiles in a range of psychrophilic, thermophilic, anaerobic, alkaliphilic, acidophilic, and halophilic conditions, as well as in environments with very high or toxic metal ion concentrations. While much additional work lies ahead to characterize the various pathways by which these biominerals form, forced biomineralization has been shown to provide insights for the progression of extreme biomimetics, allowing for promising new forays into creating the next generation of composites using organic-templating approaches under biologically extreme laboratory conditions relevant to a wide range of industrial conditions.

FGF signalling plays similar roles in development and regeneration of the skeleton in the brittle star Amphiura filiformis

Regeneration as an adult developmental process is in many aspects similar to embryonic development. Although many studies point out similarities and differences, no large-scale, direct and functional comparative analyses between development and regeneration of a specific cell type or structure in one animal exist. Here, we use the brittle star Amphiura filiformis to characterise the role of the FGF signalling pathway during skeletal development in embryos and arm regeneration. In both processes, we find ligands expressed in ectodermal cells that flank underlying skeletal mesenchymal cells, which express the receptors. Perturbation of FGF signalling showed inhibited skeleton formation in both embryogenesis and regeneration, without affecting other key developmental processes. Differential transcriptome analysis finds mostly differentiation genes rather than transcription factors to be downregulated in both contexts. Moreover, comparative gene analysis allowed us to discover brittle star-specific differentiation genes. In conclusion, our results show that the FGF pathway is crucial for skeletogenesis in the brittle star, as in other deuterostomes, and provide evidence for the re-deployment of a developmental gene regulatory module during regeneration.

Somite Compartments in Amphioxus and Its Implications on the Evolution of the Vertebrate Skeletal Tissues

Mineralized skeletal tissues of vertebrates are an evolutionary novelty within the chordate lineage. While the progenitor cells that contribute to vertebrate skeletal tissues are known to have two embryonic origins, the mesoderm and neural crest, the evolutionary origin of their developmental process remains unclear. Using cephalochordate amphioxus as our model, we found that cells at the lateral wall of the amphioxus somite express SPARC (a crucial gene for tissue mineralization) and various collagen genes. During development, some of these cells expand medially to surround the axial structures, including the neural tube, notochord and gut, while others expand laterally and ventrally to underlie the epidermis. Eventually these cell populations are found closely associated with the collagenous matrix around the neural tube, notochord, and dorsal aorta, and also with the dense collagen sheets underneath the epidermis. Using known genetic markers for distinct vertebrate somite compartments, we showed that the lateral wall of amphioxus somite likely corresponds to the vertebrate dermomyotome and lateral plate mesoderm. Furthermore, we demonstrated a conserved role for BMP signaling pathway in somite patterning of both amphioxus and vertebrates. These results suggest that compartmentalized somites and their contribution to primitive skeletal tissues are ancient traits that date back to the chordate common ancestor. The finding of SPARC-expressing skeletal scaffold in amphioxus further supports previous hypothesis regarding SPARC gene family expansion in the elaboration of the vertebrate mineralized skeleton.

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.

Evolution of Protein-Mediated Biomineralization in Scleractinian Corals

While recent strides have been made in understanding the biological process by which stony corals calcify, much remains to be revealed, including the ubiquity across taxa of specific biomolecules involved. Several proteins associated with this process have been identified through proteomic profiling of the skeletal organic matrix (SOM) extracted from three scleractinian species. However, the evolutionary history of this putative “biomineralization toolkit,” including the appearance of these proteins’ throughout metazoan evolution, remains to be resolved. Here we used a phylogenetic approach to examine the evolution of the known scleractinians’ SOM proteins across the Metazoa. Our analysis reveals an evolutionary process dominated by the co-option of genes that originated before the cnidarian diversification. Each one of the three species appears to express a unique set of the more ancient genes, representing the independent co-option of SOM proteins, as well as a substantial proportion of proteins that evolved independently. In addition, in some instances, the different species expressed multiple orthologous proteins sharing the same evolutionary history. Furthermore, the non-random clustering of multiple SOM proteins within scleractinian-specific branches suggests the conservation of protein function between distinct species for what we posit is part of the scleractinian “core biomineralization toolkit.” This “core set” contains proteins that are likely fundamental to the scleractinian biomineralization mechanism. From this analysis, we infer that the scleractinians’ ability to calcify was achieved primarily through multiple lineage-specific protein expansions, which resulted in a new functional role that was not present in the parent gene.

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.

Acidic Monosaccharides become Incorporated into Calcite Single Crystals*

Carbohydrates, along with proteins and peptides, are known to represent a major class of biomacromolecules involved in calcium carbonate biomineralization. However, in spite of multiple physical and biochemical characterizations, the explicit role of saccharide macromolecules (long chains of carbohydrate molecules) in mineral deposition is not yet understood. In this study, we investigated the influence of two common acidic monosaccharides (MSs), the two simplest forms of acidic carbohydrates, namely glucuronic and galacturonic acids, on the formation of calcite crystals in vitro. We show here that the size, morphology, and microstructure of calcite crystals are altered when they are grown in the presence of these MSs. More importantly, these MSs were found to become incorporated into the calcite crystalline lattice and induce anisotropic lattice distortions, a phenomenon widely studied for other biomolecules related to CaCO₃ biomineralization, but never before reported in the case of single MSs. Changes in the calcite lattice induced by MSs incorporation were precisely determined by high-resolution synchrotron powder X-ray diffraction. We believe that the results of this research may deepen our understanding of the interaction of saccharide polymers with an inorganic host and shed light on the implications of carbohydrates for biomineralization processes.

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.

Integrated proteomic and transcriptomic analysis reveals that polymorphic shell colors vary with melanin synthesis in Bellamya purificata snail

The snail Bellamya purificata is an ecologically and economically important freshwater gastropod species. However, limited genomic resources are available for this snail. In this study, the transcriptome of mantle tissues and proteome of shells of B. purificata with two shell colors (namely light-cyan line (LC) and light-purple line (LP)) were deeply sequenced and characterized. A total of 5.72 million contigs were assembled into 157,015 unigenes, 21,455 (13.66%) of these unigenes were significantly matched to NR, Swiss-Prot, KOG, GO and KEGG database. 1807 differentially expressed genes (DEGs) were identified between the two different shell color lines. These DEGs were significantly enriched in five KEGG pathways including tyrosine metabolism, tryptophan metabolism, phenylalanine, tyrosine and tryptophan biosynthesis, phenylalanine metabolism, and histidine metabolism, which suggested that the shell color polymorphism in B. purificata was a result of melanin synthesis variation. A total of 1521 proteins were identified in B. purificata here as well. The differentially expressed protein analysis showed that the tyrosinase content in LP was significantly decreased in comparison to LC, which agreed with the transcriptome analysis results. This study provides valuable genomic resources of B. purificata and improves our understanding of molecular mechanisms of biomineralization and shell color polymorphism in snail.

Genomic insights of body plan transitions from bilateral to pentameral symmetry in Echinoderms

Echinoderms are an exceptional group of bilaterians that develop pentameral adult symmetry from a bilaterally symmetric larva. However, the genetic basis in evolution and development of this unique transformation remains to be clarified. Here we report newly sequenced genomes, developmental transcriptomes, and proteomes of diverse echinoderms including the green sea urchin (L. variegatus), a sea cucumber (A. japonicus), and with particular emphasis on a sister group of the earliest-diverged echinoderms, the feather star (A. japonica). We learned that the last common ancestor of echinoderms retained a well-organized Hox cluster reminiscent of the hemichordate, and had gene sets involved in endoskeleton development. Further, unlike in other animal groups, the most conserved developmental stages were not at the body plan establishing phase, and genes normally involved in bilaterality appear to function in pentameric axis development. These results enhance our understanding of the divergence of protostomes and deuterostomes almost 500 Mya.

Li et al. investigate the evolution and genetic basis of the adult pentameral body plan in echinoderms using genomic, transcriptomic, and proteomic data. They determine that the last common ancestor of echinoderms contained an organized Hox cluster and endoskeleton genes, and suggest that cooption of bilateral development genes was involved in evolution of the pentameric body plan.

Molecular mechanisms of biomineralization in marine invertebrates

Much recent marine research has been directed towards understanding the effects of anthropogenic-induced environmental change on marine biodiversity, particularly for those animals with heavily calcified exoskeletons, such as corals, molluscs and urchins. This is because life in our oceans is becoming more challenging for these animals with changes in temperature, pH and salinity. In the future, it will be more energetically expensive to make marine skeletons and the increasingly corrosive conditions in seawater are expected to result in the dissolution of these external skeletons. However, initial predictions of wide-scale sensitivity are changing as we understand more about the mechanisms underpinning skeletal production (biomineralization). These studies demonstrate the complexity of calcification pathways and the cellular responses of animals to these altered conditions. Factors including parental conditioning, phenotypic plasticity and epigenetics can significantly impact the production of skeletons and thus future population success. This understanding is paralleled by an increase in our knowledge of the genes and proteins involved in biomineralization, particularly in some phyla, such as urchins, molluscs and corals. This Review will provide a broad overview of our current understanding of the factors affecting skeletal production in marine invertebrates. It will focus on the molecular mechanisms underpinning biomineralization and how knowledge of these processes affects experimental design and our ability to predict responses to climate change. Understanding marine biomineralization has many tangible benefits in our changing world, including improvements in conservation and aquaculture and exploitation of natural calcified structure design using biomimicry approaches that are aimed at producing novel biocomposites.

Skeletal development in the sea urchin relies upon protein families that contain intrinsic disorder, aggregation-prone, and conserved globular interactive domains

The formation of the sea urchin spicule skeleton requires the participation of hydrogel-forming protein families that regulate mineral nucleation and nanoparticle assembly processes that give rise to the spicule. However, the structure and molecular behavior of these proteins is not well established, and thus our ability to understand this process is hampered. We embarked on a study of sea urchin spicule proteins using a combination of biophysical and bioinformatics techniques. Our biophysical findings indicate that recombinant variants of the two most studied spicule matrix proteins, SpSM50 and SpSM30B/C (S. purpuratus) have a conformational landscape that include a C-terminal random coil/intrinsically disordered MAPQG sequence coupled to a conserved, folded N-terminal C-type lectin-like (CTLL) domain, with SpSM50 > SpSM30B/C with regard to intrinsic disorder. Both proteins possess solvent-accessible unfolded MAQPG sequence regions where Asn, Gln, and Arg residues may be accessible for protein hydrogel interactions with water molecules. Our bioinformatics study included seven other spicule matrix proteins where we note similarities between these proteins and rare, unusual proteins that possess folded and unfolded traits. Moreover, spicule matrix proteins possess three types of sequences: intrinsically disordered, amyloid-like, and folded protein-protein interactive. Collectively these reactive domains would be capable of driving protein assembly and hydrogel formation. Interestingly, three types of global conformations are predicted for the nine member protein set, wherein we note variations in the arrangement of intrinsically disordered and interactive globular domains. These variations may reflect species-specific requirements for spiculogenesis. We conclude that the molecular landscape of spicule matrix protein families enables them to function as hydrogelators, nucleators, and assemblers of mineral nanoparticles.

Glycan Binding Profiling of Jacalin-Related Lectins from the Pteria Penguin Pearl Shell

We determined the primary structures of jacalin-related lectins termed PPL3s (PPL3A, 3B, and 3C, which are dimers consisting of sequence variants α + α, α + β, β + β, respectively) and PPL4, which is heterodimer consisting of α + β subunits, isolated from mantle secretory fluid of Pteria penguin (Mabe) pearl shell. Their carbohydrate-binding properties were analyzed, in addition to that of PPL2A, which was previously reported as a matrix protein. PPL3s and PPL4 shared only 35–50% homology to PPL2A, respectively; they exhibited significantly different carbohydrate-binding specificities based on the multiple glycan binding profiling data sets from frontal affinity chromatography analysis. The carbohydrate-binding specificity of PPL3s was similar to that of PPL2A, except only for Man₃Fuc₁Xyl₁GlcNAc₂ oligosaccharide, while PPL4 showed different carbohydrate-binding specificity compared with PPL2A and PPL3s. PPL2A and PPL3s mainly recognize agalactosylated- and galactosylated-type glycans. On the other hand, PPL4 binds to high-mannose-and hybrid-type N-linked glycans but not agalactosylated- and galactosylated-type glycans.

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.

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.

Whole transcriptome sequencing and biomineralization gene architecture associated with cultured pearl quality traits in the pearl oyster, Pinctada margaritifera

Background

Cultured pearls are unique gems produced by living organisms, mainly molluscs of the Pinctada genus, through the biomineralization properties of pearl sac tissue. Improvement of P. margaritifera pearl quality is one of the biggest challenges that Polynesian research has faced to date. To achieve this goal, a better understanding of the complex mechanisms related to nacre and pearl formation is essential and can now be approached through the use of massive parallel sequencing technologies. The aim of this study was to use RNA-seq to compare whole transcriptome expression of pearl sacs that had producing pearls with high and low quality. For this purpose, a comprehensive reference transcriptome of P. margaritifera was built based on multi-tissue sampling (mantle, gonad, whole animal), including different living stages (juvenile, adults) and phenotypes (colour morphotypes, sex).

Results

Strikingly, few genes were found to be up-regulated for high quality pearls (n = 16) compared to the up-regulated genes in low quality pearls (n = 246). Biomineralization genes up-regulated in low quality pearls were specific to prismatic and prism-nacre layers. Alternative splicing was further identified in several key biomineralization genes based on a recent P. margaritifera draft genome.

Conclusion

This study lifts the veil on the multi-level regulation of biomineralization genes associated with pearl quality determination.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5443-5) contains supplementary material, which is available to authorized users.

Insights into the Evolution of Shells and Love Darts of Land Snails Revealed from Their Matrix Proteins

Over the past decade, many skeletal matrix proteins that are possibly related to calcification have been reported in various calcifying animals. Molluscs are among the most diverse calcifying animals and some gastropods have adapted to terrestrial ecological niches. Although many shell matrix proteins (SMPs) have already been reported in molluscs, most reports have focused on marine molluscs, and the SMPs of terrestrial snails remain unclear. In addition, some terrestrial stylommatophoran snails have evolved an additional unique calcified character, called a "love dart," used for mating behavior. We identified 54 SMPs in the terrestrial snail Euhadra quaesita, and found that they contain specific domains that are widely conserved in molluscan SMPs. However, our results also suggest that some of them possibly have evolved independently by domain shuffling, domain recruitment, or gene co-option. We then identified four dart matrix proteins, and found that two of them are the same proteins as those identified as SMPs. Our results suggest that some dart matrix proteins possibly have evolved by independent gene co-option from SMPs during dart evolution events. These results provide a new perspective on the evolution of SMPs and "love darts" in land snails.

Developmental effector gene regulation: Multiplexed strategies for functional analysis

The staggering complexity of the genome controls for developmental processes is revealed through massively parallel cis-regulatory analysis using new methods of perturbation and readout. The choice of combinations of these new methods is tailored to the system, question and resources at hand. Our focus is on issues that include the necessity or sufficiency of given cis-regulatory modules, cis-regulatory function in the normal spatial genomic context, and easily accessible high throughput and multiplexed analysis methods. In the sea urchin embryonic model, recombineered BACs offer new opportunities for consecutive modes of cis-regulatory analyses that answer these requirements, as we here demonstrate on a diverse suite of previously unstudied sea urchin effector genes expressed in skeletogenic cells. Positively active cis-regulatory modules were located in single Nanostring experiments per BAC containing the gene of interest, by application of our previously reported "barcode" tag vectors of which> 100 can be analyzed at one time. Computational analysis of DNA sequences that drive expression, based on the known skeletogenic regulatory state, then permitted effective identification of functional target site clusters. Deletion of these sub-regions from the parent BACs revealed module necessity, as simultaneous tests of the same regions in short constructs revealed sufficiency. Predicted functional inputs were then confirmed by site mutations, all generated and tested in multiplex formats. There emerged the simple conclusion that each effector gene utilizes a small subset of inputs from the skeletogenic GRN. These inputs may function to only adjust expression levels or in some cases necessary for expression. Since we know the GRN architecture upstream of the effector genes, we could then conceptually isolate and compare the wiring of the effector gene driver sub-circuits and identify the inputs whose removal abolish expression.

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.

Secrets of the Sea Urchin Spicule Revealed: Protein Cooperativity Is Responsible for ACC Transformation, Intracrystalline Incorporation, and Guided Mineral Particle Assembly in Biocomposite Material Formation

The formation of the sea urchin spicule involves the stabilization and transformation of amorphous calcium carbonate (ACC) and assembly of ACC nanoparticle precursors into a mesoscale single crystal of fracture-resistant calcite. This process of particle assembly or attachment is under the control of a family of proteins known as the spicule matrix [Strongylocentrotus purpuratus (SpSM)] proteome. Recently, two members of this proteome, SpSM50 and the glycoprotein SpSM30B/C-G (in recombinant forms), were found to interact together via SpSM30B/C-G oligosaccharide-SpSM50 protein interactions to form hybrid protein hydrogels with unique physical properties. In this study, we investigate the mineralization properties of this hybrid hydrogel alongside the hydrogels formed by SpSM50 and SpSM30B/C-G individually. We find that the SpSM50 + SpSM30B/C-G hybrid hydrogel is synergistic with regard to surface modifications and intracrystalline inclusions of existing calcite crystals, the inhibition of ACC formation, and the kinetic destabilization of ACC to form a crystalline phase. Most importantly, the hybrid hydrogel phase assembles and organizes mineral particles into discrete clusters or domains within in vitro mineralization environments. Thus, the interactions of SpSM50 and SpSM30B/C-G, mediated by carbohydrate-protein binding, reflect the need for protein cooperativity for the ACC-to-crystalline transformation, intracrystalline void formation, and guided mineral particle assembly processes that are instrumental in spicule formation.

Exploring the drivers of early biomineralization

The first biomineralized hard parts are known from ∼810 Million years ago (Ma), consisting of phosphatic plates of probable protists formed under active biological control. Large skeletons in diverse taxa, probably including total-group poriferans and total-group cnidarians, first appear in the terminal Ediacaran, ∼550 Ma. This is followed by a substantial increase in abundance, diversity and mineralogy during the early Cambrian. The biological relationship of Ediacaran to early Cambrian skeletal biota is unclear, but tubular skeletal fossils such as Cloudina and Anabarites straddle the transition. Many Ediacaran skeletal biota are found exclusively in carbonate settings, and present skeletons whose form infers an organic scaffold which provided the framework for interactions between extracellular matrix and mineral ions. Several taxa have close soft-bodied counterparts hosted in contemporary clastic rocks. This supports the assertion that the calcification was an independent and derived feature that appeared in diverse groups, which was initially acquired with minimal biological control in the highly saturated, high-alkalinity carbonate settings of the Ediacaran, where the carbonate polymorph was further controlled by seawater chemistry. The trigger for Ediacaran-Cambrian biomineralization is far from clear, but may have been either changing seawater Mg/Ca ratios that facilitated widespread aragonite and high-Mg calcite precipitation, and/or increasing or stabilizing oxygen levels. By the Early Cambrian, the diversity of biomineralization styles may have been an escalating defensive response to increasing predation pressure, with skeletal hard parts first appearing in abundance in clastic settings by the Fortunian. This marks full independence from ambient seawater chemistry and significant biological control of biomineralization.

Glycosylation Fosters Interactions between Model Sea Urchin Spicule Matrix Proteins. Implications for Embryonic Spiculogenesis and Biomineralization

The formation of embryonic mineralized skeletal elements (spicules) in the sea urchin requires the participation of proteins, many of which may interact with one another and assist in the creation of an extracellular matrix wherein mineral formation takes place. To probe this, we created a sea urchin spicule recombinant model protein pair system wherein we tested the interactions between two major spicule proteins, SpSM50 and the glycoprotein, SpSM30B/C. Both proteins are strong hydrogelators that manipulate early and later events in mineral formation. We discovered that the anionic glycan moieties of SpSM30B/C are required for interaction with the SpSM50 protein and that these interactions are Ca(II)-independent. In addition, when these proteins form a complex, they create hybrid hydrogel particles that are physically distinct from their individual counterparts. Thus, glycan-mediated interactions play an important role in in vitro spicule protein assembly and most likely within the spicule itself.

Global analysis of primary mesenchyme cell cis-regulatory modules by chromatin accessibility profiling

Background

The developmental gene regulatory network (GRN) that underlies skeletogenesis in sea urchins and other echinoderms is a paradigm of GRN structure, function, and evolution. This transcriptional network is deployed selectively in skeleton-forming primary mesenchyme cells (PMCs) of the early embryo. To advance our understanding of this model developmental GRN, we used genome-wide chromatin accessibility profiling to identify and characterize PMC cis-regulatory modules (CRMs).

Results

ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) analysis of purified PMCs provided a global picture of chromatin accessibility in these cells. We used both ATAC-seq and DNase-seq (DNase I hypersensitive site sequencing) to identify > 3000 sites that exhibited increased accessibility in PMCs relative to other embryonic cell lineages, and provide both computational and experimental evidence that a large fraction of these sites represent bona fide skeletogenic CRMs. Putative PMC CRMs were preferentially located near genes differentially expressed by PMCs and consensus binding sites for two key transcription factors in the PMC GRN, Alx1 and Ets1, were enriched in these CRMs. Moreover, a high proportion of candidate CRMs drove reporter gene expression specifically in PMCs in transgenic embryos. Surprisingly, we found that PMC CRMs were partially open in other embryonic lineages and exhibited hyperaccessibility as early as the 128-cell stage.

Conclusions

Our work provides a comprehensive picture of chromatin accessibility in an early embryonic cell lineage. By identifying thousands of candidate PMC CRMs, we significantly enhance the utility of the sea urchin skeletogenic network as a general model of GRN architecture and evolution. Our work also shows that differential chromatin accessibility, which has been used for the high-throughput identification of enhancers in differentiated cell types, is a powerful approach for the identification of CRMs in early embryonic cells. Lastly, we conclude that in the sea urchin embryo, CRMs that control the cell type-specific expression of effector genes are hyperaccessible several hours in advance of gene activation.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4542-z) contains supplementary material, which is available to authorized users.

Developmental transcriptomics of the brittle star Amphiura filiformis reveals gene regulatory network rewiring in echinoderm larval skeleton evolution

Background

Amongst the echinoderms the class Ophiuroidea is of particular interest for its phylogenetic position, ecological importance and developmental and regenerative biology. However, compared to other echinoderms, notably echinoids (sea urchins), relatively little is known about developmental changes in gene expression in ophiuroids. To address this issue, we have generated and assembled a large RNAseq data set of four key stages of development in the brittle star Amphiura filiformis and a de novo reference transcriptome of comparable quality to that of a model echinoderm—the sea urchin Strongylocentrotus purpuratus. Furthermore, we provide access to the new data via a web interface: http://www.echinonet.eu/shiny/Amphiura_filiformis/.

Results

We have identified highly conserved genes associated with the development of a biomineralised skeleton. We also identify important class-specific characters, including the independent duplication of the msp130 class of genes in different echinoderm classes and the unique occurrence of spicule matrix (sm) genes in echinoids. Using a new quantification pipeline for our de novo transcriptome, validated with other methodologies, we find major differences between brittle stars and sea urchins in the temporal expression of many transcription factor genes. This divergence in developmental regulatory states is more evident in early stages of development when cell specification begins, rather than when cells initiate differentiation.

Conclusions

Our findings indicate that there has been a high degree of gene regulatory network rewiring and clade-specific gene duplication, supporting the hypothesis of a convergent evolution of larval skeleton development in echinoderms.

Electronic supplementary material

The online version of this article (10.1186/s13059-018-1402-8) contains supplementary material, which is available to authorized users.

First macrobiota biomineralization was environmentally triggered

Why large and diverse skeletons first appeared ca 550 Ma is not well understood. Many Ediacaran skeletal biota show evidence of flexibility, and bear notably thin skeletal walls with simple, non-hierarchical microstructures of either aragonite or high-Mg calcite. We present evidence that the earliest skeletal macrobiota, found only in carbonate rocks, had close soft-bodied counterparts hosted in contemporary clastic rocks. This includes the calcareous discoidal fossil Suvorovella, similar to holdfasts of Ediacaran biota taxa previously known only as casts and moulds, as well as tubular and vase-shaped fossils. In sum, these probably represent taxa of diverse affinity including unicellular eukaryotes, total group cnidarians and problematica. Our findings support the assertion that the calcification was an independent and derived feature that appeared in diverse groups where an organic scaffold was the primitive character, which provided the framework for interactions between the extracellular matrix and mineral ions. We conclude that such skeletons may have been acquired with relative ease in the highly saturated, high alkalinity carbonate settings of the Ediacaran, where carbonate polymorph was further controlled by seawater chemistry. The trigger for Ediacaran biomineralization may have been either changing seawater Mg/Ca and/or increasing oxygen levels. By the Early Cambrian, however, biomineralization styles and the range of biominerals had significantly diversified, perhaps as an escalating defensive response to increasing predation pressure. Indeed skeletal hardparts had appeared in clastic settings by Cambrian Stage 1, suggesting independence from ambient seawater chemistry where genetic and molecular mechanisms controlled biomineralization and mineralogy had become evolutionarily constrained.