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

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

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

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

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

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

Sea Urchin Spicule Matrix Proteins Form Mesoscale "Smart" Hydrogels That Exhibit Selective Ion Interactions

In the sea urchin embryo spicule, there exists a proteome of >200 proteins that are responsible for controlling the mineralization of the spicule and the formation of a fracture-resistant composite. In this report, using recombinant proteins, we identify that two protein components of the spicule, SM30B/C and SM50, are hydrogelators. Because of the presence of intrinsic disorder and aggregation-prone regions, these proteins assemble to form porous mesoscale hydrogel particles in solution. These hydrogel particles change their size, organization, and internal structure in response to pH and ions, particularly Ca(II), which indicates that these behave as ion-responsive or "smart" hydrogels. Using diffusion-ordered spectroscopy NMR, we find that both hydrogels affect the diffusion of water, but only SM50 affects the diffusion of an anionic solute. Thus, the extracellular matrix of the spicule consists of several hydrogelator proteins which are responsive to solution conditions and can control the diffusion of water and solutes, and these proteins will serve as a model system for designing ion-responsive, composite, and smart hydrogels.

Endocytosis in primary mesenchyme cells during sea urchin larval skeletogenesis

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Ocean acidification impacts spine integrity but not regenerative capacity of spines and tube feet in adult sea urchins

Increasing atmospheric carbon dioxide (CO2) has resulted in a change in seawater chemistry and lowering of pH, referred to as ocean acidification. Understanding how different organisms and processes respond to ocean acidification is vital to predict how marine ecosystems will be altered under future scenarios of continued environmental change. Regenerative processes involving biomineralization in marine calcifiers such as sea urchins are predicted to be especially vulnerable. In this study, the effect of ocean acidification on regeneration of external appendages (spines and tube feet) was investigated in the sea urchin Lytechinus variegatus exposed to ambient (546 µatm), intermediate (1027 µatm) and high (1841 µatm) partial pressure of CO2 (pCO2) for eight weeks. The rate of regeneration was maintained in spines and tube feet throughout two periods of amputation and regrowth under conditions of elevated pCO2. Increased expression of several biomineralization-related genes indicated molecular compensatory mechanisms; however, the structural integrity of both regenerating and homeostatic spines was compromised in high pCO2 conditions. Indicators of physiological fitness (righting response, growth rate, coelomocyte concentration and composition) were not affected by increasing pCO2, but compromised spine integrity is likely to have negative consequences for defence capabilities and therefore survival of these ecologically and economically important organisms.

De novo transcriptome analysis of the excretory tubules of Carausius morosus (Phasmatodea) and possible functions of the midgut 'appendices'

The Malpighian tubules are the insect excretory organs, responsible for ion and water homeostasis and elimination of nitrogenous wastes. Post-genomic assays suggest they also metabolize and detoxify xenobiotic compounds and have antimicrobial properties. The Phasmatodea have an additional, unique set of excretory organs referred to predominantly as midgut appendices. Their function and how it compares to phasmid and other insect Malpighian tubules is unknown. Hypotheses include carbonic anhydrase activity, calcium and metal cation sequestration, and xenobiotic transport. This work presents the first comparative transcriptomic analysis of the Phasmatodean excretory organs, using the model insect Carausius morosus. I produced de novo transcriptomes of the midgut appendices, midgut wall, and Malpighian tubules, and looked for differentially expressed genes associated with putative organ functions. The appendices differentially and highly express lipid transport and metabolism proteins, and the biomineralization gene otopetrin. The Malpighian tubules differentially and highly express acid phosphatases and multiple transporter types, while appendices express fat-soluble vitamin and peptide transporters. Many defense proteins such as multidrug resistance proteins, ABC transporters, cytochrome P450's, and glutathione-S-transferases were differentially expressed in specific excretory organs. I hypothesize that the appendices and Malpighian tubules both have defensive / xenobiotic metabolism functions, but each likely target different substrates. Phasmid Malpighian tubules excrete as in other insects, while the appendices may predominantly regulate amino acids, fats, and fat-soluble compounds. Lipid metabolism in insects is poorly understood, and the Phasmatodea may thus serve as a model for studying this further.

Transcriptomic responses to seawater acidification among sea urchin populations inhabiting a natural pH mosaic

Increasing awareness of spatial and temporal variation in ocean pH suggests some marine populations may be adapted to local pH regimes and will therefore respond differently to present-day pH variation and to long-term ocean acidification. In the Northeast Pacific Ocean, differences in the strength of coastal upwelling cause latitudinal variation in prevailing pH regimes that are hypothesized to promote local adaptation and unequal pH tolerance among resident populations. In this study, responses to experimental seawater acidification were compared among embryos and larvae from six populations of purple sea urchins (Strongylocentrotus purpuratus) inhabiting areas that differ in their frequency of low pH exposure and that prior research suggests are locally adapted to seawater pH. Transcriptomic analyses demonstrate urchin populations most frequently exposed to low pH seawater responded to experimental acidification by expressing genes within major ATP-producing pathways at greater levels than populations encountering low pH less often. Multiple genes within the tricarboxylic acid cycle, electron transport chain and fatty acid beta oxidation pathways were upregulated in urchin populations experiencing low pH conditions most frequently. These same metabolic pathways were significantly over-represented among genes both expressed in a population-specific manner and putatively under selection to enhance low pH tolerance. Collectively, these data suggest natural selection is acting on metabolic gene networks to redirect ATP toward maintaining acid-base homeostasis and enhance tolerance of seawater acidification. As a trade-off, marine populations more tolerant of low pH may have less energy to put towards other aspects of fitness and to respond to additional ocean change.

A Model Sea Urchin Spicule Matrix Protein Self-Associates To Form Mineral-Modifying Protein Hydrogels

In the purple sea urchin Strongylocentrotus purpuratus, the formation and mineralization of fracture-resistant skeletal elements such as the embryonic spicule require the combinatorial participation of numerous spicule matrix proteins such as the SpSM30A-F isoforms. However, because of limited abundance, it has been difficult to pursue extensive biochemical studies of the SpSM30 proteins and deduce their role in spicule formation and mineralization. To circumvent these problems, we expressed a model recombinant spicule matrix protein, rSpSM30B/C, which possesses the key sequence attributes of isoforms "B" and "C". Our findings indicate that rSpSM30B/C is expressed in insect cells as a single polypeptide containing variations in glycosylation that create microheterogeneity in rSpSM30B/C molecular masses. These post-translational modifications incorporate O- and N-glycans and anionic mono- and bisialylated and mono- and bisulfated monosaccharides on the protein molecules and enhance its aggregation propensity. Bioinformatics and biophysical experiments confirm that rSpSM30B/C is an intrinsically disordered, aggregation-prone protein that forms porous protein hydrogels that control the in vitro mineralization process in three ways: (1) increase the time interval for prenucleation cluster formation and transiently stabilize an ACC polymorph, (2) promote and organize single-crystal calcite nanoparticles, and (3) promote faceted growth and create surface texturing of calcite crystals. These features are also common to mollusk shell nacre proteins, and we conclude that rSpSM30B/C is a spiculogenesis protein that exhibits traits found in other calcium carbonate mineral modification proteins.

The Widespread Prevalence and Functional Significance of Silk-Like Structural Proteins in Metazoan Biological Materials

In nature, numerous mechanisms have evolved by which organisms fabricate biological structures with an impressive array of physical characteristics. Some examples of metazoan biological materials include the highly elastic byssal threads by which bivalves attach themselves to rocks, biomineralized structures that form the skeletons of various animals, and spider silks that are renowned for their exceptional strength and elasticity. The remarkable properties of silks, which are perhaps the best studied biological materials, are the result of the highly repetitive, modular, and biased amino acid composition of the proteins that compose them. Interestingly, similar levels of modularity/repetitiveness and similar bias in amino acid compositions have been reported in proteins that are components of structural materials in other organisms, however the exact nature and extent of this similarity, and its functional and evolutionary relevance, is unknown. Here, we investigate this similarity and use sequence features common to silks and other known structural proteins to develop a bioinformatics-based method to identify similar proteins from large-scale transcriptome and whole-genome datasets. We show that a large number of proteins identified using this method have roles in biological material formation throughout the animal kingdom. Despite the similarity in sequence characteristics, most of the silk-like structural proteins (SLSPs) identified in this study appear to have evolved independently and are restricted to a particular animal lineage. Although the exact function of many of these SLSPs is unknown, the apparent independent evolution of proteins with similar sequence characteristics in divergent lineages suggests that these features are important for the assembly of biological materials. The identification of these characteristics enable the generation of testable hypotheses regarding the mechanisms by which these proteins assemble and direct the construction of biological materials with diverse morphologies. The SilkSlider predictor software developed here is available at https://github.com/wwood/SilkSlider.

Sea shell diversity and rapidly evolving secretomes: insights into the evolution of biomineralization

An external skeleton is an essential part of the body plan of many animals and is thought to be one of the key factors that enabled the great expansion in animal diversity and disparity during the Cambrian explosion. Molluscs are considered ideal to study the evolution of biomineralization because of their diversity of highly complex, robust and patterned shells. The molluscan shell forms externally at the interface of animal and environment, and involves controlled deposition of calcium carbonate within a framework of macromolecules that are secreted from the dorsal mantle epithelium. Despite its deep conservation within Mollusca, the mantle is capable of producing an incredible diversity of shell patterns, and macro- and micro-architectures. Here we review recent developments within the field of molluscan biomineralization, focusing on the genes expressed in the mantle that encode secreted proteins. The so-called mantle secretome appears to regulate shell deposition and patterning and in some cases becomes part of the shell matrix. Recent transcriptomic and proteomic studies have revealed marked differences in the mantle secretomes of even closely-related molluscs; these typically exceed expected differences based on characteristics of the external shell. All mantle secretomes surveyed to date include novel genes encoding lineage-restricted proteins and unique combinations of co-opted ancient genes. A surprisingly large proportion of both ancient and novel secreted proteins containing simple repetitive motifs or domains that are often modular in construction. These repetitive low complexity domains (RLCDs) appear to further promote the evolvability of the mantle secretome, resulting in domain shuffling, expansion and loss. RLCD families further evolve via slippage and other mechanisms associated with repetitive sequences. As analogous types of secreted proteins are expressed in biomineralizing tissues in other animals, insights into the evolution of the genes underlying molluscan shell formation may be applied more broadly to understanding the evolution of metazoan biomineralization.

Characterization of an Alpha Type Carbonic Anhydrase from Paracentrotus lividus Sea Urchin Embryos

Carbonic anhydrases (CA) are zinc metalloenzymes that catalyze the reversible hydration of carbon dioxide to bicarbonate. In the sea urchin, CA has a role in the formation of the calcitic skeleton during embryo development. Here, we report a newly identified mRNA sequence from embryos of the sea urchin Paracentrotus lividus, referred to as Pl-can. The complete coding sequence was identified with the aid of both EST databases and experimental procedures. Pl-CAN is a 447 aa-long protein, with an estimated molecular mass of 48.5 kDa and an isoelectric point of 6.83. The in silico study of functional domains showed, in addition to the alpha type CA-specific domain, the presence of an unexpected glycine-rich region at the N-terminal of the molecule. This is not found in any other species described so far, but probably it is restricted to the sea urchins. The phylogenetic analysis indicated that Pl-CAN is evolutionarily closer to human among chordates than to other species. The putative role(s) of the identified domains is discussed. The Pl-can temporal and spatial expression profiles, analyzed throughout embryo development by comparative qPCR and whole-mount in situ hybridization (WMISH), showed that Pl-can mRNA is specifically expressed in the primary mesenchyme cells (PMC) of the embryo and levels increase along with the growth of the embryonic skeleton, reaching a peak at the pluteus stage. A recombinant fusion protein was produced in E. coli and used to raise specific antibodies in mice recognized the endogenous Pl-CAN by Western blot in embryo extracts from gastrula and pluteus.

Comparative Developmental Transcriptomics Reveals Rewiring of a Highly Conserved Gene Regulatory Network during a Major Life History Switch in the Sea Urchin Genus Heliocidaris

The ecologically significant shift in developmental strategy from planktotrophic (feeding) to lecithotrophic (nonfeeding) development in the sea urchin genus Heliocidaris is one of the most comprehensively studied life history transitions in any animal. Although the evolution of lecithotrophy involved substantial changes to larval development and morphology, it is not known to what extent changes in gene expression underlie the developmental differences between species, nor do we understand how these changes evolved within the context of the well-defined gene regulatory network (GRN) underlying sea urchin development. To address these questions, we used RNA-seq to measure expression dynamics across development in three species: the lecithotroph Heliocidaris erythrogramma, the closely related planktotroph H. tuberculata, and an outgroup planktotroph Lytechinus variegatus. Using well-established statistical methods, we developed a novel framework for identifying, quantifying, and polarizing evolutionary changes in gene expression profiles across the transcriptome and within the GRN. We found that major changes in gene expression profiles were more numerous during the evolution of lecithotrophy than during the persistence of planktotrophy, and that genes with derived expression profiles in the lecithotroph displayed specific characteristics as a group that are consistent with the dramatically altered developmental program in this species. Compared to the transcriptome, changes in gene expression profiles within the GRN were even more pronounced in the lecithotroph. We found evidence for conservation and likely divergence of particular GRN regulatory interactions in the lecithotroph, as well as significant changes in the expression of genes with known roles in larval skeletogenesis. We further use coexpression analysis to identify genes of unknown function that may contribute to both conserved and derived developmental traits between species. Collectively, our results indicate that distinct evolutionary processes operate on gene expression during periods of life history conservation and periods of life history divergence, and that this contrast is even more pronounced within the GRN than across the transcriptome as a whole.

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

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

Reorganization of sea urchin gene regulatory networks at least 268 million years ago as revealed by oldest fossil cidaroid echinoid

Echinoids, or sea urchins, are rare in the Palaeozoic fossil record, and thus the details regarding the early diversification of crown group echinoids are unclear. Here we report on the earliest probable crown group echinoid from the fossil record, recovered from Permian (Roadian-Capitanian) rocks of west Texas, which has important implications for the timing of the divergence of crown group echinoids. The presence of apophyses and rigidly sutured interambulacral areas with two columns of plates indicates this species is a cidaroid echinoid. The species, Eotiaris guadalupensis, n. sp. is therefore the earliest stem group cidaroid. The occurrence of this species in Roadian strata pushes back the divergence of cidaroids and euechinoids, the clades that comprise all living echinoids, to at least 268.8 Ma, ten million years older than the previously oldest known cidaroid. Furthermore, the genomic regulation of development in echinoids is amongst the best known, and this new species informs the timing of large-scale reorganization in echinoid gene regulatory networks that occurred at the cidaroid-euechinoid divergence, indicating that these changes took place by the Roadian stage of the Permian.

Genome-wide assessment of differential effector gene use in embryogenesis

Six different populations of cells were isolated by fluorescence-activated cell sorting from disaggregated late blastula- and gastrula-stage sea urchin embryos according to the regulatory states expressed in these cells, as reported by recombineered bacterial artificial chromosomes producing fluorochromes. Transcriptomes recovered from these embryonic cell populations revealed striking, early differential expression of large cohorts of effector genes. The six cell populations were presumptive pigment cells, presumptive neurogenic cells, presumptive skeletogenic cells, cells from the stomodeal region of the oral ectoderm, ciliated band cells and cells from the endoderm/ectoderm boundary that will give rise both to hindgut and to border ectoderm. Transcriptome analysis revealed that each of these domains specifically expressed several hundred effector genes at significant levels. Annotation indicated the qualitative individuality of the functional nature of each cell population, even though they were isolated from embryos only 1-2 days old. In no case was more than a tiny fraction of the transcripts enriched in one population also enriched in any other of the six populations studied. As was particularly clear in the cases of the presumptive pigment, neurogenic and skeletogenic cells, all three of which represent precociously differentiating cell types of this embryo, most specifically expressed genes of given cell types are not significantly expressed at all in the other cell types. Thus, at the effector gene level, a dramatic, cell type-specific pattern of differential gene regulation is established well before any significant embryonic morphogenesis has occurred.

The Lingula genome provides insights into brachiopod evolution and the origin of phosphate biomineralization

The evolutionary origins of lingulid brachiopods and their calcium phosphate shells have been obscure. Here we decode the 425-Mb genome of Lingula anatina to gain insights into brachiopod evolution. Comprehensive phylogenomic analyses place Lingula close to molluscs, but distant from annelids. The Lingula gene number has increased to ∼34,000 by extensive expansion of gene families. Although Lingula and vertebrates have superficially similar hard tissue components, our genomic, transcriptomic and proteomic analyses show that Lingula lacks genes involved in bone formation, indicating an independent origin of their phosphate biominerals. Several genes involved in Lingula shell formation are shared by molluscs. However, Lingula has independently undergone domain combinations to produce shell matrix collagens with EGF domains and carries lineage-specific shell matrix proteins. Gene family expansion, domain shuffling and co-option of genes appear to be the genomic background of Lingula's unique biomineralization. This Lingula genome provides resources for further studies of lophotrochozoan evolution.

Lingulid brachiopods possess calcium phosphate shells. Here, the authors sequence the genome of Lingula anatine to show that Lingula is evolutionary close to molluscs, but distant from annelids, and identify the genomic background of Lingula's unique biomineralization mechanism.

H(+)/K(+) ATPase activity is required for biomineralization in sea urchin embryos

The bioelectrical signatures associated with regeneration, wound healing, development, and cancer are changes in the polarization state of the cell that persist over long durations, and are mediated by ion channel activity. To identify physiologically relevant bioelectrical changes that occur during normal development of the sea urchin Lytechinus variegatus, we tested a range of ion channel inhibitors, and thereby identified SCH28080, a chemical inhibitor of the H(+)/K(+) ATPase (HKA), as an inhibitor of skeletogenesis. In sea urchin embryos, the primary mesodermal lineage, the PMCs, produce biomineral in response to signals from the ectoderm. However, in SCH28080-treated embryos, aside from randomization of the left-right axis, the ectoderm is normally specified and differentiated, indicating that the block to skeletogenesis observed in SCH28080-treated embryos is PMC-specific. HKA inhibition did not interfere with PMC specification, and was sufficient to block continuing biomineralization when embryos were treated with SCH28080 after the initiation of skeletogenesis, indicating that HKA activity is continuously required during biomineralization. Ion concentrations and voltage potential were abnormal in the PMCs in SCH28080-treated embryos, suggesting that these bioelectrical abnormalities prevent biomineralization. Our results indicate that this effect is due to the inhibition of amorphous calcium carbonate precipitation within PMC vesicles.

SM50 repeat-polypeptides self-assemble into discrete matrix subunits and promote appositional calcium carbonate crystal growth during sea urchin tooth biomineralization

The two major proteins involved in vertebrate enamel formation and echinoderm sea urchin tooth biomineralization, amelogenin and SM50, are both characterized by elongated polyproline repeat domains in the center of the macromolecule. To determine the role of polyproline repeat polypeptides in basal deuterostome biomineralization, we have mapped the localization of SM50 as it relates to crystal growth, conducted self-assembly studies of SM50 repeat polypeptides, and examined their effect on calcium carbonate and apatite crystal growth. Electron micrographs of the growth zone of Strongylocentrotus purpuratus sea urchin teeth documented a series of successive events from intravesicular mineral nucleation to mineral deposition at the interface between tooth surface and odontoblast syncytium. Using immunohistochemistry, SM50 was detected within the cytoplasm of cells associated with the developing tooth mineral, at the mineral secreting front, and adjacent to initial mineral deposits, but not in muscles and ligaments. Polypeptides derived from the SM50 polyproline alternating hexa- and hepta-peptide repeat region (SM50P6P7) formed highly discrete, donut-shaped self-assembly patterns. In calcium carbonate crystal growth studies, SM50P6P7 repeat peptides triggered the growth of expansive networks of fused calcium carbonate crystals while in apatite growth studies, SM50P6P7 peptides facilitated the growth of needle-shaped and parallel arranged crystals resembling those found in developing vertebrate enamel. In comparison, SM50P6P7 surpassed the PXX24 polypeptide repeat region derived from the vertebrate enamel protein amelogenin in its ability to promote crystal nucleation and appositional crystal growth. Together, these studies establish the SM50P6P7 polyproline repeat region as a potent regulator in the protein-guided appositional crystal growth that occurs during continuous tooth mineralization and eruption. In addition, our studies highlight the role of species-specific polyproline repeat motifs in the formation of discrete self-assembled matrices and the resulting control of mineral growth.

Carbonic anhydrase inhibition blocks skeletogenesis and echinochrome production in Paracentrotus lividus and Heliocidaris tuberculata embryos and larvae

Carbonic anhydrases (CAs) are a family of widely distributed metalloenzymes, involved in diverse physiological processes. These enzymes catalyse the reversible conversion of carbon dioxide to protons and bicarbonate. At least 19 genes encoding for CAs have been identified in the sea urchin genome, with one of these localized to the skeletogenic mesoderm (primary mesenchyme cells, PMCs). We investigated the effects of a specific inhibitor of CA, acetazolamide (AZ), on development of two sea urchin species with contrasting investment in skeleton production, Paracentrotus lividus and Heliocidaris tuberculata, to determine the role of CA on PMC differentiation, skeletogenesis and on non-skeletogenic mesodermal (NSM) cells. Embryos were cultured in the presence of AZ from the blastula stage prior to skeleton formation and development to the larval stage was monitored. At the dose of 8 mmol/L AZ, 98% and 90% of P. lividus and H. tuberculata embryos lacked skeleton, respectively. Nevertheless, an almost normal PMC differentiation was indicated by the expression of msp130, a PMC-specific marker. Strikingly, the AZ-treated embryos also lacked the echinochrome pigment produced by the pigment cells, a subpopulation of NSM cells with immune activities within the larva. Conversely, all ectoderm and endoderm derivatives and other subpopulations of mesoderm developed normally. The inhibitory effects of AZ were completely reversed after removal of the inhibitor from the medium. Our data, together with new information concerning the involvement of CA on skeleton formation, provide evidence for the first time of a possible role of the CAs in larval immune pigment cells.

Spine and test skeletal matrices of the Mediterranean sea urchin Arbacia lixula--a comparative characterization of their sugar signature

Calcified structures of sea urchins are biocomposite materials that comprise a minor fraction of organic macromolecules, such as proteins, glycoproteins and polysaccharides. These macromolecules are thought to collectively regulate mineral deposition during the process of calcification. When occluded, they modify the properties of the mineral. In the present study, the organic matrices (both soluble and insoluble in acetic acid) of spines and tests from the Mediterranean black sea urchin Arbacia lixula were extracted and characterized, in order to determine whether they exhibit similar biochemical signatures. Bulk characterizations were performed by mono-dimensional SDS/PAGE, FT-IR spectroscopy, and an in vitro crystallization assay. We concentrated our efforts on characterization of the sugar moieties. To this end, we determined the monosaccharide content of the soluble and insoluble organic matrices of A. lixula spines and tests by HPAE-PAD, together with their respective lectin-binding profiles via enzyme-linked lectin assay. Finally, we performed in situ localization of N-acetyl glucosamine-containing saccharides on spines and tests using gold-conjugated wheatgerm agglutinin. Our data show that the test and spine matrices exhibit different biochemical signatures with regard to their saccharidic fraction, suggesting that future studies should analyse the regulation of mineral deposition by the matrix in these two mineralized structures in detail. This study re-emphasizes the importance of non-protein moieties, i.e. sugars, in calcium carbonate systems, and highlights the need to clearly identify their function in the biomineralization process.

Examination of the skeletal proteome of the brittle star Ophiocoma wendtii reveals overall conservation of proteins but variation in spicule matrix proteins

Background

While formation of mineralized tissue is characteristic of many animal taxa, the proteins that interact with mineral are diverse and appear in many cases to be of independent origin. Extracellular matrix proteins involved in mineralization do share some common features. They tend to be disordered, secreted proteins with repetitive, low complexity. The genes encoding these proteins are often duplicated and undergo concerted evolution, further diversifying the repetitive domains. This makes it difficult to identify mineralization genes and the proteins they encode using bioinformatics techniques. Here we describe the use of proteomics to identify mineralization genes in an ophiuroid echinoderm, Ophiocoma wendtii (O. wendtii).

Results

We have isolated the occluded proteins within the mineralized tissue of the brittle star Ophiocoma wendtii. The proteins were analyzed both unfractionated and separated on SDS-PAGE gels. Each slice was analyzed using mass spectroscopy and the amino acid sequence of the most prevalent peptides was obtained. This was compared to both an embryonic transcriptome from the gastrula stage when skeleton is being formed and a tube foot (an adult mineralized tissue) transcriptome. Thirty eight proteins were identified which matched known proteins or protein domains in the NCBI databases. These include C-type lectins, ECM proteins, Kazal-type protease inhibitors, matrix metalloproteases as well as more common cellular proteins. Many of these are similar to those found in the sea urchin Strongylocentrotus purpuratus (S. purpuratus) skeleton. We did not, however, identify clear homologs to the sea urchin spicule matrix proteins, and the number of C-type lectin containing genes was much reduced compared to sea urchins. Also notably absent was MSP-130.

Conclusions

Our results show an overall conservation of the types of proteins found in the mineralized tissues of two divergent groups of echinoderms, as well as in mineralized tissues in general. However, the extensive gene duplication and concerted evolution seen in the spicule matrix proteins found in the sea urchin skeleton was not observed in the brittle star.

Electronic supplementary material

The online version of this article (doi:10.1186/s12953-015-0064-7) contains supplementary material, which is available to authorized users.

Juvenile skeletogenesis in anciently diverged sea urchin clades

Mechanistic understanding of evolutionary divergence in animal body plans devolves from analysis of those developmental processes that, in forms descendant from a common ancestor, are responsible for their morphological differences. The last common ancestor of the two extant subclasses of sea urchins, i.e., euechinoids and cidaroids, existed well before the Permian/Triassic extinction (252 mya). Subsequent evolutionary divergence of these clades offers in principle a rare opportunity to solve the developmental regulatory events underlying a defined evolutionary divergence process. Thus (i) there is an excellent and fairly dense (if yet incompletely analyzed) fossil record; (ii) cladistically confined features of the skeletal structures of modern euechinoid and cidaroid sea urchins are preserved in fossils of ancestral forms; (iii) euechinoids and cidaroids are among current laboratory model systems in molecular developmental biology (here Strongylocentrotus purpuratus [Sp] and Eucidaris tribuloides [Et]); (iv) skeletogenic specification in sea urchins is uncommonly well understood at the causal level of interactions of regulatory genes with one another, and with known skeletogenic effector genes, providing a ready arsenal of available molecular tools. Here we focus on differences in test and perignathic girdle skeletal morphology that distinguish all modern euechinoid from all modern cidaroid sea urchins. We demonstrate distinct canonical test and girdle morphologies in juveniles of both species by use of SEM and X-ray microtomography. Among the sharply distinct morphological features of these clades are the internal skeletal structures of the perignathic girdle to which attach homologous muscles utilized for retraction and protraction of Aristotles׳ lantern and its teeth. We demonstrate that these structures develop de novo between one and four weeks after metamorphosis. In order to study the underlying developmental processes, a method of section whole mount in situ hybridization was adapted. This method displays current gene expression in the developing test and perignathic girdle skeletal elements of both Sp and Et juveniles. Active, specific expression of the sm37 biomineralization gene in these muscle attachment structures accompanies morphogenetic development of these clade-specific features in juveniles of both species. Skeletogenesis at these clade-specific muscle attachment structures displays molecular earmarks of the well understood embryonic skeletogenic GRN: thus the upstream regulatory gene alx1 and the gene encoding the vegfR signaling receptor are both expressed at the sites where they are formed. This work opens the way to analysis of the alternative spatial specification processes that were installed at the evolutionary divergence of the two extant subclasses of sea urchins.

Calcareous sponge genomes reveal complex evolution of α-carbonic anhydrases and two key biomineralization enzymes

BACKGROUND

Calcium carbonate biominerals form often complex and beautiful skeletal elements, including coral exoskeletons and mollusc shells. Although the ability to generate these carbonate structures was apparently gained independently during animal evolution, it sometimes involves the same gene families. One of the best-studied of these gene families comprises the α- carbonic anhydrases (CAs), which catalyse the reversible transformation of CO2 to HCO3 - and fulfill many physiological functions. Among Porifera -the oldest animal phylum with the ability to produce skeletal elements- only the class of calcareous sponges can build calcitic spicules, which are the extracellular products of specialized cells, the sclerocytes. Little is known about the molecular mechanisms of their synthesis, but inhibition studies suggest an essential role of CAs. In order to gain insight into the evolution and function of CAs in biomineralization of a basal metazoan species, we determined the diversity and expression of CAs in the calcareous sponges Sycon ciliatum and Leucosolenia complicata by means of genomic screening, RNA-Seq and RNA in situ hybridization expression analysis. Active biomineralization was located with calcein-staining.

RESULTS

We found that the CA repertoires of two calcareous sponge species are strikingly more complex than those of other sponges. By characterizing their expression patterns, we could link two CAs (one intracellular and one extracellular) to the process of calcite spicule formation in both studied species. The extracellular biomineralizing CAs seem to be of paralogous origin, a finding that advises caution against assuming functional conservation of biomineralizing genes based upon orthology assessment alone. Additionally, calcareous sponges possess acatalytic CAs related to human CAs X and XI, suggesting an ancient origin of these proteins. Phylogenetic analyses including CAs from genomes of all non-bilaterian phyla suggest multiple gene losses and duplications and presence of several CAs in the last common ancestor of metazoans.

CONCLUSIONS

We identified two key biomineralization enzymes from the CA-family in calcareous sponges and propose their possible interaction in spicule formation. The complex evolutionary history of the CA family is driven by frequent gene diversification and losses. These evolutionary patterns likely facilitated the numerous events of independent recruitment of CAs into biomineralization within Metazoa.

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

The endoskeleton of the sea urchin embryo is produced by primary mesenchyme cells (PMCs). Maternal inputs activate a complex gene regulatory network (GRN) in the PMC lineage in a cell-autonomous fashion during early development, initially creating a uniform population of prospective skeleton-forming cells. Previous studies showed that at post-blastula stages of development, several effector genes in the network exhibit non-uniform patterns of expression, suggesting that their regulation becomes subject to local, extrinsic cues. Other studies have identified the VEGF and MAPK pathways as regulators of PMC migration, gene expression, and biomineralization. In this study, we used whole mount in situ hybridization (WMISH) to examine the spatial expression patterns of 39 PMC-specific/enriched mRNAs in Strongylocentrotus purpuratus embryos at the late gastrula, early prism and pluteus stages. We found that all 39 mRNAs (including several regulatory genes) showed non-uniform patterns of expression within the PMC syncytium, revealing a global shift in the regulation of the skeletogenic GRN from a cell-autonomous to a signal-dependent mode. In general, localized regions of elevated gene expression corresponded to sites of rapid biomineral deposition. We used a VEGFR inhibitor (axitinib) and a MEK inhibitor (U0126) to show that VEGF signaling and the MAPK pathway are essential for maintaining high levels of gene expression in PMCs at the tips of rods that extend from the ventral region of the embryo. These inhibitors affected gene expression in the PMCs in similar ways, suggesting that VEGF acts via the MAPK pathway. In contrast, axitinib and U0126 did not affect the localized expression of genes in PMCs at the tips of the body rods, which form on the dorsal side of the embryo. Our results therefore indicate that multiple signaling pathways regulate the skeletogenic GRN during late stages of embryogenesis-VEGF/MAPK signaling on the ventral side and a separate, unidentified pathway on the dorsal side. These two signaling pathways appear to be activated sequentially (ventral followed by dorsal) and many effector genes are subject to regulation by both pathways.

Specification to biomineralization: following a single cell type as it constructs a skeleton

The sea urchin larva is shaped by a calcite endoskeleton. That skeleton is built by 64 primary mesenchyme cells (PMCs) in Lytechinus variegatus. The PMCs originate as micromeres due to an unequal fourth cleavage in the embryo. Micromeres are specified in a well-described molecular sequence and enter the blastocoel at a precise time using a classic epithelial-mesenchymal transition. To make the skeleton, the PMCs receive signaling inputs from the overlying ectoderm, which provides positional information as well as control of the growth of initial skeletal tri-radiates. The patterning of the skeleton is the result both of autonomous inputs from PMCs, including production of proteins that are included in the skeletal matrix, and of non-autonomous dynamic information from the ectoderm. Here, we summarize the wealth of information known about how a PMC contributes to the skeletal structure. The larval skeleton is a model for understanding how information encoded in DNA is translated into a three-dimensional crystalline structure.

Investigating protein function in biomineralized tissues using molecular biology techniques

We describe modern molecular biology methods currently used in the study of biomineralization. We focus our descriptions on two areas of biomineralization research in which these methods have been particularly powerful. The first area is the use of modern molecular methods to identify and characterize the so-called occluded matrix proteins present in mineralized tissues. More specifically, we describe the use of RNA-seq and the next generation of DNA sequencers and the use of direct protein sequencing and mass spectrometers as ways of identifying proteins present in mineralized tissues. The second area is the use of molecular methods to examine the function of proteins in biomineralization. RNA interference (RNAi), morpholino antisense, and other methods are described and discussed as ways of elucidating protein function.

The skeleton of postmetamorphic echinoderms in a changing world

Available evidence on the impact of acidification and its interaction with warming on the skeleton of postmetamorphic (juvenile and adult) echinoderms is reviewed. Data are available on sea urchins, starfish, and brittle stars in 33 studies. Skeleton growth of juveniles of all sea urchin species studied so far is affected from pH 7.8 to 7.6 in seawater, values that are expected to be reached during the 21st century. Growth in adult sea urchins (six species studied) is apparently only marginally affected at seawater pH relevant to this century. The interacting effect of temperature differed according to studies. Juvenile starfish as well as adults seem to be either not impacted or even boosted by acidification. Brittle stars show moderate effects at pH below or equal to 7.4. Dissolution of the body wall skeleton is unlikely to be a major threat to sea urchins. Spines, however, due to their exposed position, are more prone to this threat, but their regeneration abilities can probably ensure their maintenance, although this could have an energetic cost and induce changes in resource allocation. No information is available on skeleton dissolution in starfish, and the situation in brittle stars needs further assessment. Very preliminary evidence indicates that mechanical properties in sea urchins could be affected. So, although the impact of ocean acidification on the skeleton of echinoderms has been considered as a major threat from the first studies, we need a better understanding of the induced changes, in particular the functional consequences of growth modifications and dissolution related to mechanical properties. It is suggested to focus studies on these aspects.

Effects of seawater acidification on gene expression: resolving broader-scale trends in sea urchins

Sea urchins are ecologically and economically important calcifying organisms threatened by acidification of the global ocean caused by anthropogenic CO2 emissions. Propelled by the sequencing of the purple sea urchin (Strongylocentrotus purpuratus) genome, profiling changes in gene expression during exposure to high pCO2 seawater has emerged as a powerful and increasingly common method to infer the response of urchins to ocean change. However, analyses of gene expression are sensitive to experimental methodology, and comparisons between studies of genes regulated by ocean acidification are most often made in the context of major caveats. Here we perform meta-analyses as a means of minimizing experimental discrepancies and resolving broader-scale trends regarding the effects of ocean acidification on gene expression in urchins. Analyses across eight studies and four urchin species largely support prevailing hypotheses about the impact of ocean acidification on marine calcifiers. The predominant expression pattern involved the down-regulation of genes within energy-producing pathways, a clear indication of metabolic depression. Genes with functions in ion transport were significantly over-represented and are most plausibly contributing to intracellular pH regulation. Expression profiles provided extensive evidence for an impact on biomineralization, epitomized by the down-regulation of seven spicule matrix proteins. In contrast, expression profiles provided limited evidence for CO2-mediated developmental delay or induction of a cellular stress response. Congruence between studies of gene expression and the ocean acidification literature in general validates the accuracy of gene expression in predicting the consequences of ocean change and justifies its continued use in future studies.

Lectin uptake and incorporation into the calcitic spicule of sea urchin embryos

Primary mesenchyme cells (PMCs) are skeletogenenic cells that produce a calcareous endoskeleton in developing sea urchin larvae. The PMCs fuse to form a cavity in which spicule matrix proteins and calcium are secreted forming the mineralized spicule. In this study, living sea urchin embryos were stained with fluorescently conjugated wheat germ agglutinin, a lectin that preferentially binds to PMCs, and the redistribution of this fluorescent tag was examined during sea urchin development. Initially, fluorescence was associated primarily with the surface of PMCs. Subsequently, the fluorescent label redistributed to intracellular vesicles in the PMCs. As the larval skeleton developed, intracellular granular staining diminished and fluorescence appeared in the spicules. Spicules that were cleaned to remove membranous material associated with the surface exhibited bright fluorescence, which indicated that fluorescently labelled lectin had been incorporated into the spicule matrix. The results provide evidence for a cellular pathway in which material is taken up at the cell surface, sequestered in intracellular vesicles and then incorporated into the developing spicule.

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

A central challenge of developmental and evolutionary biology is to understand the transformation of genetic information into morphology. Elucidating the connections between genes and anatomy will require model morphogenetic processes that are amenable to detailed analysis of cell/tissue behaviors and to systems-level approaches to gene regulation. The formation of the calcified endoskeleton of the sea urchin embryo is a valuable experimental system for developing such an integrated view of the genomic regulatory control of morphogenesis. A transcriptional gene regulatory network (GRN) that underlies the specification of skeletogenic cells (primary mesenchyme cells, or PMCs) has recently been elucidated. In this study, we carried out a genome-wide analysis of mRNAs encoded by effector genes in the network and uncovered transcriptional inputs into many of these genes. We used RNA-seq to identify >400 transcripts differentially expressed by PMCs during gastrulation, when these cells undergo a striking sequence of behaviors that drives skeletal morphogenesis. Our analysis expanded by almost an order of magnitude the number of known (and candidate) downstream effectors that directly mediate skeletal morphogenesis. We carried out genome-wide analysis of (1) functional targets of Ets1 and Alx1, two pivotal, early transcription factors in the PMC GRN, and (2) functional targets of MAPK signaling, a pathway that plays an essential role in PMC specification. These studies identified transcriptional inputs into >200 PMC effector genes. Our work establishes a framework for understanding the genomic regulatory control of a major morphogenetic process and has important implications for reconstructing the evolution of biomineralization in metazoans.

Metazoan calcium carbonate biomineralizations: macroevolutionary trends – challenges for the coming decade

Calcium carbonate-based biominerals, also referred as biocalcifications, are the most abundant biogenic mineralized products at the surface of the Earth. In this paper, we summarize general concepts on biocalcifications and we sketch macro-evolutionary trends throughout the history of the Earth, from Archean to Phanerozoic times. Then, we expose five fundamental issues that represent key-challenges in biocalcification researches for the coming decade: the first one concerns the comprehension of the micro- and nano-structure of calcium carbonate biominerals from a mineral viewpoint, while the second one deals with the understanding of the dynamic process of their fabrication. The third one treats the subtle interplay between organics and the mineral phase. The fourth issue focuses on an environmental challenge related to ocean acidification (OA); at last, the diagenetic processes that affect biogenic calcium carbonate mineral constitute the fifth issue.

General approach for in vivo recovery of cell type-specific effector gene sets

Differentially expressed, cell type-specific effector gene sets hold the key to multiple important problems in biology, from theoretical aspects of developmental gene regulatory networks (GRNs) to various practical applications. Although individual cell types of interest have been recovered by various methods and analyzed, systematic recovery of multiple cell type-specific gene sets from whole developing organisms has remained problematic. Here we describe a general methodology using the sea urchin embryo, a material of choice because of the large-scale GRNs already solved for this model system. This method utilizes the regulatory states expressed by given cells of the embryo to define cell type and includes a fluorescence activated cell sorting (FACS) procedure that results in no perturbation of transcript representation. We have extensively validated the method by spatial and qualitative analyses of the transcriptome expressed in isolated embryonic skeletogenic cells and as a consequence, generated a prototypical cell type-specific transcriptome database.

The echinoderm larval skeleton as a possible model system for experimental evolutionary biology

The evolution of various body plans results from the acquisition of novel structures as well as the loss of existing structures. Some novel structures necessitate multiple evolutionary steps, requiring organisms to overcome the intermediate steps, which might be less adaptive or neutral. To examine this issue, echinoderms might provide an ideal experimental system. A larval skeleton is acquired in some echinoderm lineages, such as sea urchins, probably via the co-option of the skeletogenic machinery that was already established to produce the adult skeleton. The acquisition of a larval skeleton was found to require multiple steps and so provides a model experimental system for reproducing intermediate evolutionary stages. The fact that echinoderm embryology has been studied with various natural populations also presents an advantage.

Branching out: origins of the sea urchin larval skeleton in development and evolution

It is a challenge to understand how the information encoded in DNA is used to build a three-dimensional structure. To explore how this works the assembly of a relatively simple skeleton has been examined at multiple control levels. The skeleton of the sea urchin embryo consists of a number of calcite rods produced by 64 skeletogenic cells. The ectoderm supplies spatial cues for patterning, essentially telling the skeletogenic cells where to position themselves and providing the factors for skeletal growth. Here, we describe the information known about how this works. First the ectoderm must be patterned so that the signaling cues are released from precise positions. The skeletogenic cells respond by initiating skeletogenesis immediately beneath two regions (one on the right and the other on the left side). Growth of the skeletal rods requires additional signaling from defined ectodermal locations, and the skeletogenic cells respond to produce a membrane-bound template in which the calcite crystal grows. Important in this process are three signals, fibroblast growth factor, vascular endothelial growth factor, and Wnt5. Each is necessary for explicit tasks in skeleton production.

Crystallization and preliminary X-ray analysis of the C-type lectin domain of the spicule matrix protein SM50 from Strongylocentrotus purpuratus

Sea urchin spicules have a calcitic mesocrystalline architecture that is closely associated with a matrix of proteins and amorphous minerals. The mechanism underlying spicule formation involves complex processes encompassing spatio-temporally regulated organic-inorganic interactions. C-type lectin domains are present in several spicule matrix proteins in Strongylocentrotus purpuratus, implying their role in spiculogenesis. In this study, the C-type lectin domain of SM50 was overexpressed, purified and crystallized using a vapour-diffusion method. The crystal diffracted to a resolution of 2.85 Å and belonged to space group P212121, with unit-cell parameters a = 100.6, b = 115.4, c = 130.6 Å, α = β = γ = 90°. Assuming 50% solvent content, six chains are expected to be present in the asymmetric unit.

The impact of gene expression variation on the robustness and evolvability of a developmental gene regulatory network

Regulatory interactions buffer development against genetic and environmental perturbations, but adaptation requires phenotypes to change. We investigated the relationship between robustness and evolvability within the gene regulatory network underlying development of the larval skeleton in the sea urchin Strongylocentrotus purpuratus. We find extensive variation in gene expression in this network throughout development in a natural population, some of which has a heritable genetic basis. Switch-like regulatory interactions predominate during early development, buffer expression variation, and may promote the accumulation of cryptic genetic variation affecting early stages. Regulatory interactions during later development are typically more sensitive (linear), allowing variation in expression to affect downstream target genes. Variation in skeletal morphology is associated primarily with expression variation of a few, primarily structural, genes at terminal positions within the network. These results indicate that the position and properties of gene interactions within a network can have important evolutionary consequences independent of their immediate regulatory role.

Expression of skeletogenic genes during arm regeneration in the brittle star Amphiura filiformis

•

Analysis of brittle star regenerating arms using differentiation markers.

•

Identification of the early segregation of skeletal and muscle progenitor cells.

•

Expression of skeletal and non-skeletal genes at different stages of regeneration.

•

Combinatorial role of TF genes in early specification of skeletal cells.

•

Same TF genes identify different skeletal structures later in regeneration.

The brittle star Amphiura filiformis, which regenerates its arms post autotomy, is emerging as a useful model for studying the molecular underpinnings of regeneration, aided by the recent availability of some molecular resources. During regeneration a blastema initially is formed distally to the amputation site, and then a rapid rebuild is obtained by adding metameric units, which will eventually differentiate and become fully functional. In this work we first characterize the developmental process of the regenerating arms using two differentiation markers for muscle and skeletal structures – Afi-trop-1 and Afi-αcoll. Both genes are not expressed in the blastema and newly added undifferentiated metameric units. Their expression at different regenerating stages shows an early segregation of muscle and skeletal cells during the regenerating process, long before the metameric units become functional. We then studied the expression of a set of genes orthologous of the sea urchin transcription factors involved in the development of skeletal and non-skeletal mesoderm: Afi-ets1/2, Afi-alx1, Afi-tbr, Afi-foxB and Afi-gataC. We found that Afi-ets1/2, Afi-alx1, Afi-foxB and Afi-gataC are all expressed at the blastemal stage. As regeneration progresses those genes are expressed in a similar small undifferentiated domain beneath the distal growth cap, while in more advanced metameric units they become restricted to different skeletal domains. Afi-foxB becomes expressed in non-skeletal structures. This suggests that they might play a combinatorial role only in the early cell specification process and that subsequently they function independently in the differentiation of different structures. Afi-tbr is not present in the adult arm tissue at any stage of regeneration. In situ hybridization results have been confirmed with a new strategy for quantitative PCR (QPCR), using a subdivision of the three stages of regeneration into proximal (differentiated) and distal (undifferentiated) arm segments.

Roles of larval sea urchin spicule SM50 domains in organic matrix self-assembly and calcium carbonate mineralization

The larval spicule matrix protein SM50 is the most abundant occluded matrix protein present in the mineralized larval sea urchin spicule. Recent evidence implicates SM50 in the stabilization of amorphous calcium carbonate (ACC). Here, we investigate the molecular interactions of SM50 and CaCO3 by investigating the function of three major domains of SM50 as small ubiquitin-like modifier (SUMO) fusion proteins - a C-type lectin domain (CTL), a glycine rich region (GRR) and a proline rich region (PRR). Under various mineralization conditions, we find that SUMO-CTL is monomeric and influences CaCO3 mineralization, SUMO-GRR aggregates into large protein superstructures and SUMO-PRR modifies the early CaCO3 mineralization stages as well as growth. The combination of these mineralization and self-assembly properties of the major domains synergistically enable the full-length SM50 to fulfill functions of constructing the organic spicule matrix as well as performing necessary mineralization activities such as Ca(2+) ion recruitment and organization to allow for proper growth and development of the mineralized larval sea urchin spicule.

SM30 protein function during sea urchin larval spicule formation

A central issue in better understanding the process of biomineralization is to elucidate the function of occluded matrix proteins present in mineralized tissues. A potent approach to addressing this issue utilizes specific inhibitors of expression of known genes. Application of antisense oligonucleotides that specifically suppress translation of a given mRNA are capable of causing aberrant biomineralization, thereby revealing, at least in part, a likely function of the protein and gene under investigation. We have applied this approach to study the possible function(s) of the SM30 family of proteins, which are found in spicules, teeth, spines, and tests of Strongylocentrotus purpuratus as well as other euechinoid sea urchins. It is possible using the anti-SM30 morpholino-oligonucleotides (MO's) to reduce the level of these proteins to very low levels, yet the development of skeletal spicules in the embryo shows little or no aberration. This surprising result requires re-thinking about the role of these, and possibly other occluded matrix proteins.

Evolutionary change during experimental ocean acidification

Rising atmospheric carbon dioxide (CO2) conditions are driving unprecedented changes in seawater chemistry, resulting in reduced pH and carbonate ion concentrations in the Earth's oceans. This ocean acidification has negative but variable impacts on individual performance in many marine species. However, little is known about the adaptive capacity of species to respond to an acidified ocean, and, as a result, predictions regarding future ecosystem responses remain incomplete. Here we demonstrate that ocean acidification generates striking patterns of genome-wide selection in purple sea urchins (Strongylocentrotus purpuratus) cultured under different CO2 levels. We examined genetic change at 19,493 loci in larvae from seven adult populations cultured under realistic future CO2 levels. Although larval development and morphology showed little response to elevated CO2, we found substantial allelic change in 40 functional classes of proteins involving hundreds of loci. Pronounced genetic changes, including excess amino acid replacements, were detected in all populations and occurred in genes for biomineralization, lipid metabolism, and ion homeostasis--gene classes that build skeletons and interact in pH regulation. Such genetic change represents a neglected and important impact of ocean acidification that may influence populations that show few outward signs of response to acidification. Our results demonstrate the capacity for rapid evolution in the face of ocean acidification and show that standing genetic variation could be a reservoir of resilience to climate change in this coastal upwelling ecosystem. However, effective response to strong natural selection demands large population sizes and may be limited in species impacted by other environmental stressors.