cis-Regulatory control of cyclophilin, a member of the ETS-DRI skeletogenic gene battery in the sea urchin embryo

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.

Architecture and evolution of the cis-regulatory system of the echinoderm kirrelL gene

The gene regulatory network (GRN) that underlies echinoderm skeletogenesis is a prominent model of GRN architecture and evolution. KirrelL is an essential downstream effector gene in this network and encodes an Ig-superfamily protein required for the fusion of skeletogenic cells and the formation of the skeleton. In this study, we dissected the transcriptional control region of the kirrelL gene of the purple sea urchin, Strongylocentrotus purpuratus. Using plasmid- and bacterial artificial chromosome-based transgenic reporter assays, we identified key cis-regulatory elements (CREs) and transcription factor inputs that regulate Sp-kirrelL, including direct, positive inputs from two key transcription factors in the skeletogenic GRN, Alx1 and Ets1. We next identified kirrelL cis-regulatory regions from seven other echinoderm species that together represent all classes within the phylum. By introducing these heterologous regulatory regions into developing sea urchin embryos we provide evidence of their remarkable conservation across ~500 million years of evolution. We dissected in detail the kirrelL regulatory region of the sea star, Patiria miniata, and demonstrated that it also receives direct inputs from Alx1 and Ets1. Our findings identify kirrelL as a component of the ancestral echinoderm skeletogenic GRN. They support the view that GRN subcircuits, including specific transcription factor–CRE interactions, can remain stable over vast periods of evolutionary history. Lastly, our analysis of kirrelL establishes direct linkages between a developmental GRN and an effector gene that controls a key morphogenetic cell behavior, cell–cell fusion, providing a paradigm for extending the explanatory power of GRNs.

Pattern of Repetitive Element Transcription Segregate Cell Lineages during the Embryogenesis of Sea Urchin Strongylocentrotus purpuratus

Repetitive elements (REs) occupy a significant part of eukaryotic genomes and are shown to play diverse roles in genome regulation. During embryogenesis of the sea urchin, a large number of REs are expressed, but the role of these elements in the regulation of biological processes remains unknown. The aim of this study was to identify the RE expression at different stages of embryogenesis. REs occupied 44% of genomic DNA of Strongylocentrotus purpuratus. The most prevalent among these elements were the unknown elements-in total, they contributed 78.5% of REs (35% in total genome occupancy). It was revealed that the transcription pattern of genes and REs changes significantly during gastrulation. Using the de novo transcriptome assembly, we showed that the expression of RE is independent of its copy number in the genome. We also identified copies that are expressed. Only active RE copies were used for mapping and quantification of RE expression in the single-cell RNA sequencing data. REs expression was observed in all cell lineages and they were detected as population markers. Moreover, the primary mesenchyme cell (PMC) line had the greatest diversity of REs among the markers. Our data suggest a role for RE in the organization of developmental domains during the sea urchin embryogenesis at the single-cell resolution level.

Pearl Sac Gene Expression Profiles Associated With Pearl Attributes in the Silver-Lip Pearl Oyster, Pinctada maxima

Pearls are highly prized biomineralized gemstones produced by molluscs. The appearance and mineralogy of cultured pearls can vary markedly, greatly affecting their commercial value. To begin to understand the role of pearl sacs—organs that form in host oysters from explanted mantle tissues that surround and synthesize pearls—we undertook transcriptomic analyses to identify genes that are differentially expressed in sacs producing pearls with different surface and structural characteristics. Our results indicate that gene expression profiles correlate with different pearl defects, suggesting that gene regulation in the pearl sac contributes to pearl appearance and quality. For instance, pearl sacs that produced pearls with surface non-lustrous calcification significantly down-regulate genes associated with cilia and microtubule function compared to pearl sacs giving rise to lustrous pearls. These results suggest that gene expression profiling can advance our understanding of processes that control biomineralization, which may be of direct value to the pearl industry, particularly in relation to defects that result in low value pearls.

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.

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.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

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

Large-scale gene expression study in the ophiuroid Amphiura filiformis provides insights into evolution of gene regulatory networks

Background

The evolutionary mechanisms involved in shaping complex gene regulatory networks (GRN) that encode for morphologically similar structures in distantly related animals remain elusive. In this context, echinoderm larval skeletons found in brittle stars and sea urchins provide an ideal system. Here, we characterize for the first time the development of the larval skeleton in the ophiuroid Amphiura filiformis and compare it systematically with its counterpart in sea urchin.

Results

We show that ophiuroids and euechinoids, that split at least 480 Million years ago (Mya), have remarkable similarities in tempo and mode of skeletal development. Despite morphological and ontological similarities, our high-resolution study of the dynamics of genetic regulatory states in A. filiformis highlights numerous differences in the architecture of their underlying GRNs. Importantly, the A.filiformispplx, the closest gene to the sea urchin double negative gate (DNG) repressor pmar1, fails to drive the skeletogenic program in sea urchin, showing important evolutionary differences in protein function. hesC, the second repressor of the DNG, is co-expressed with most of the genes that are repressed in sea urchin, indicating the absence of direct repression of tbr, ets1/2, and delta in A. filiformis. Furthermore, the absence of expression in later stages of brittle star skeleton development of key regulatory genes, such as foxb and dri, shows significantly different regulatory states.

Conclusion

Our data fill up an important gap in the picture of larval mesoderm in echinoderms and allows us to explore the evolutionary implications relative to the recently established phylogeny of echinoderm classes. In light of recent studies on other echinoderms, our data highlight a high evolutionary plasticity of the same nodes throughout evolution of echinoderm skeletogenesis. Finally, gene duplication, protein function diversification, and cis-regulatory element evolution all contributed to shape the regulatory program for larval skeletogenesis in different branches of echinoderms.

Electronic supplementary material

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

Expession patterns of mesenchyme specification genes in two distantly related echinoids, Glyptocidaris crenularis and Echinocardium cordatum

The molecular mechanism of the larval mesenchyme cell specification in echinoids has been well analyzed. However, most of the data have been provided by studies of a single group of echinoids, the order Camarodonta. Little is known about this mechanism in other echinoid orders. We examined the expression patterns of mesenchyme specification genes, micro1, hesC, alx1, tbr, ets1, cyp1, and gcm, in the two non-Camarodonta echinoids, Glyptocidaris crenularis and Echinocardium cordatum. We found that the expression patterns of some genes contained characteristics that were unique to one of the species; others were shared by the two species. Some of the shared characteristics of G. crenularis and E. cordatum are not found in the species belonging to Camarodonta, suggesting the derived status of this order. The expression of ets1 in E. cordatum aboral ectoderm is one of the molecular level modifications possibly related to an evolutionarily novel larval structure, the posterior process. Our results suggest that a considerable number of modifications in the mesenchyme specification mechanisms have been introduced during the echinoid evolution.

Larval mesenchyme cell specification in the primitive echinoid occurs independently of the double-negative gate

Echinoids (sea urchins) are divided into two major groups - cidaroids (a 'primitive' group) and euechinoids (a 'derived' group). The cidaroids are a promising model species for understanding the ancestral developmental mechanisms in echinoids, but little is known about the molecular mechanisms of cidaroid development. In euechinoids, skeletogenic mesenchyme cell specification is regulated by the double-negative gate (DNG), in which hesC represses the transcription of the downstream mesenchyme specification genes (alx1, tbr and ets1), thereby defining the prospective mesenchyme region. To estimate the ancestral mechanism of larval mesenchyme cell specification in echinoids, the expression patterns and roles of mesenchyme specification genes in the cidaroid Prionocidaris baculosa were examined. The present study reveals that the expression pattern and function of hesC in P. baculosa were inconsistent with the DNG model, suggesting that the euechinoid-type DNG is not utilized during cidaroid mesenchyme specification. In contrast with hesC, the expression patterns and functions of alx1, tbr and ets1 were similar between P. baculosa and euechinoids. Based on these results, we propose that the roles of alx1, tbr and ets1 in mesenchyme specification were established in the common ancestor of echinoids, and that the DNG system was acquired in the euechinoid lineage after divergence from the cidaroid ancestor. The evolutionary timing of the establishment of the DNG suggests that the DNG was originally related to micromere and/or primary mesenchyme cell formation but not to skeletogenic cell differentiation.

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.

High ordered biomineralization induced by carbon nanoparticles in the sea urchin Paracentrotus lividus

A surprising and unexpected biomineralization process was observed during toxicological assessment of carbon nanoparticles on Paracentrotus lividus (sea urchin) pluteus larvae. The larvae activate a process of defense against external material, by incorporating the nanoparticles into microstructures of aragonite similarly to pearl oysters. Aiming at a better understanding of this phenomenon, the larvae were exposed to increasing concentrations of carbon nanoparticles and the biomineralization products were analyzed by electron microscopy, x-ray diffraction and Raman spectroscopy. In order to evaluate the possible influence of Sp-CyP-1 expression on this biomineralization process by larvae, analyses of gene expression (Sp-CyP-1) and calcein labeling were performed. Overall, we report experimental evidence about the capability of carbon nanoparticles to induce an increment of Sp-CyP-1 expression with the consequent activation of a biomineralization process leading to the production of a new pearl-like biomaterial never previously observed in sea urchins.

The evolution of the neural crest: new perspectives from lamprey and invertebrate neural crest-like cells

The neural crest is an embryonic cell population that gives rise to an array of tissues and structures in adult vertebrates including most of the head skeleton. Because neural crest cells (NCCs), and many of their derivatives, are unique to vertebrates, the evolution of the neural crest is thought to have potentiated vertebrate origins and diversification. However, the lack of clear NCC homologs in invertebrate chordates has made it difficult to reconstruct the evolutionary history of modern NCCs. In this review, the development of NCCs in the basal jawless vertebrate, lamprey, is compared with the development of neural crest-like cells in a range of invertebrates to deduce features of the first NCCs and their evolutionary precursors. These comparisons demonstrate that most of the defining attributes of NCCs are widespread features of invertebrate embryonic ectoderm. In addition, they suggest ancient origins for the neural border domain and chondroid skeletal tissue in the first bilaterian, and show that NCCs must have evolved in a chordate with an unduplicated invertebrate-type genome. On the basis of these observations, a stepwise model for the evolution of NCCs involving heterotopic and heterochronic activation of ancient ectodermal gene programs and new responsiveness to preexisting inducing signals is proposed. In light of the phylogenetic distribution of neural crest-like cells, the deep homology of developmental gene networks, and the central role of evolutionary loss in deuterostome evolution, this article concludes with suggestions for future studies in a broad range of bilaterians to test key aspects of this model. WIREs Dev Biol 2013, 2:1-15. doi: 10.1002/wdev.85 For further resources related to this article, please visit the WIREs website.

Synthetic in vivo validation of gene network circuitry

Embryonic development is controlled by networks of interacting regulatory genes. The individual linkages of gene regulatory networks (GRNs) are customarily validated by functional cis-regulatory analysis, but an additional approach to validation is to rewire GRN circuitry to test experimentally predictions derived from network structure. Here we use this synthetic method to challenge specific predictions of the sea urchin embryo endomesoderm GRN. Expression vectors generated by in vitro recombination of exogenous sequences into BACs were used to cause elements of a nonskeletogenic mesoderm GRN to be deployed in skeletogenic cells and to detect their effects. The result of reengineering the regulatory circuitry in this way was to divert the developmental program of these cells from skeletogenesis to pigment cell formation, confirming a direct prediction of the GRN. In addition, the experiment revealed previously undetected cryptic repression functions.

The genomic regulatory control of skeletal morphogenesis in the sea urchin

A central challenge of developmental and evolutionary biology is to understand how anatomy is encoded in the genome. Elucidating the genetic mechanisms that control the development of specific anatomical features will require the analysis of model morphogenetic processes and an integration of biological information at genomic, cellular and tissue levels. The formation of the endoskeleton of the sea urchin embryo is a powerful experimental system for developing such an integrated view of the genomic regulatory control of morphogenesis. The dynamic cellular behaviors that underlie skeletogenesis are well understood and a complex transcriptional gene regulatory network (GRN) that underlies the specification of embryonic skeletogenic cells (primary mesenchyme cells, PMCs) has recently been elucidated. Here, we link the PMC specification GRN to genes that directly control skeletal morphogenesis. We identify new gene products that play a proximate role in skeletal morphogenesis and uncover transcriptional regulatory inputs into many of these genes. Our work extends the importance of the PMC GRN as a model developmental GRN and establishes a unique picture of the genomic regulatory control of a major morphogenetic process. Furthermore, because echinoderms exhibit diverse programs of skeletal development, the newly expanded sea urchin skeletogenic GRN will provide a foundation for comparative studies that explore the relationship between GRN evolution and morphological evolution.

Expression pattern of polyketide synthase-2 during sea urchin development

Polyketide synthases (PKSs) are a large group of proteins responsible for the biosynthesis of polyketide compounds, which are mainly found in bacteria, fungi, and plants. Polyketides have a wide array of biological functions, including antibiotic, antifungal, predator defense, and light responses. In this study, we describe the developmental expression pattern of pks2, one of two pks found in the sea urchin genome. Throughout development, pks2 expression was restricted to skeletogenic cells and their precursors. Pks2 was first detected during the blastula stage. The transcript level peaked at hatched blastula, when all skeletogenic cell precursors expressed pks2. This was followed by a steady decline in expression in the skeletogenic cells on the aboral side of the embryo. By the prism stage, pks2 expression was limited to only 3-4 skeletogenic cells localized on the oral side.

The control of foxN2/3 expression in sea urchin embryos and its function in the skeletogenic gene regulatory network

Early development requires well-organized temporal and spatial regulation of transcription factors that are assembled into gene regulatory networks (GRNs). In the sea urchin, an endomesoderm GRN model explains much of the specification in the endoderm and mesoderm prior to gastrulation, yet some GRN connections remain incomplete. Here, we characterize FoxN2/3 in the primary mesenchyme cell (PMC) GRN state. Expression of foxN2/3 mRNA begins in micromeres at the hatched blastula stage and then is lost from micromeres at the mesenchyme blastula stage. foxN2/3 expression then shifts to the non-skeletogenic mesoderm and, later, to the endoderm. Here, we show that Pmar1, Ets1 and Tbr are necessary for activation of foxN2/3 in micromeres. The later endomesoderm expression of foxN2/3 is independent of the earlier expression of foxN2/3 in micromeres and is independent of signals from PMCs. FoxN2/3 is necessary for several steps in the formation of the larval skeleton. Early expression of genes for the skeletal matrix is dependent on FoxN2/3, but only until the mesenchyme blastula stage as foxN2/3 mRNA disappears from PMCs at that time and we assume that the protein is not abnormally long-lived. Knockdown of FoxN2/3 inhibits normal PMC ingression and foxN2/3 morphant PMCs do not organize in the blastocoel and fail to join the PMC syncytium. In addition, without FoxN2/3, the PMCs fail to repress the transfating of other mesodermal cells into the skeletogenic lineage. Thus, FoxN2/3 is necessary for normal ingression, for expression of several skeletal matrix genes, for preventing transfating and for fusion of the PMC syncytium.

Experimentally based sea urchin gene regulatory network and the causal explanation of developmental phenomenology

Gene regulatory networks for development underlie cell fate specification and differentiation. Network topology, logic and dynamics can be obtained by thorough experimental analysis. Our understanding of the gene regulatory network controlling endomesoderm specification in the sea urchin embryo has attained an advanced level such that it explains developmental phenomenology. Here we review how the network explains the mechanisms utilized in development to control the formation of dynamic expression patterns of transcription factors and signaling molecules. The network represents the genomic program controlling timely activation of specification and differentiation genes in the correct embryonic lineages. It can also be used to study evolution of body plans. We demonstrate how comparing the sea urchin gene regulatory network to that of the sea star and to that of later developmental stages in the sea urchin, reveals mechanisms underlying the origin of evolutionary novelty. The experimentally based gene regulatory network for endomesoderm specification in the sea urchin embryo provides unique insights into the system level properties of cell fate specification and its evolution.

Functional evolution of Ets in echinoderms with focus on the evolution of echinoderm larval skeletons

Convergent evolution of echinoderm pluteus larva was examined from the standpoint of functional evolution of a transcription factor Ets1/2. In sea urchins, Ets1/2 plays a central role in the differentiation of larval skeletogenic mesenchyme cells. In addition, Ets1/2 is suggested to be involved in adult skeletogenesis. Conversely, in starfish, although no skeletogenic cells differentiate during larval development, Ets1/2 is also expressed in the larval mesoderm. Here, we confirmed that the starfish Ets1/2 is indispensable for the differentiation of the larval mesoderm. This result led us to assume that, in the common ancestors of echinoderms, Ets1/2 activates the transcription of distinct gene sets, one for the differentiation of the larval mesoderm and the other for the development of the adult skeleton. Thus, the acquisition of the larval skeleton involved target switching of Ets1/2. Specifically, in the sea urchin lineage, Ets1/2 activated a downstream target gene set for skeletogenesis during larval development in addition to a mesoderm target set. We examined whether this heterochronic activation of the skeletogenic target set was achieved by the molecular evolution of the Ets1/2 transcription factor itself. We tested whether starfish Ets1/2 induced skeletogenesis when injected into sea urchin eggs. We found that, in addition to ectopic induction of mesenchyme cells, starfish Ets1/2 can activate some parts of the skeletogenic pathway in these mesenchyme cells. Thus, we suggest that the nature of the transcription factor Ets1/2 did not change, but rather that some unidentified co-factor(s) for Ets1/2 may distinguish between targets for the larval mesoderm and for skeletogenesis. Identification of the co-factor(s) will be key to understanding the molecular evolution underlying the evolution of the pluteus larvae.

Proteomic analysis of sea urchin (Strongylocentrotus purpuratus) spicule matrix

BACKGROUND

The sea urchin embryo has been an important model organism in developmental biology for more than a century. This is due to its relatively simple construction, translucent appearance, and the possibility to follow the fate of individual cells as development to the pluteus larva proceeds. Because the larvae contain tiny calcitic skeletal elements, the spicules, they are also important model organisms for biomineralization research. Similar to other biominerals the spicule contains an organic matrix, which is thought to play an important role in its formation. However, only few spicule matrix proteins were identified previously.

RESULTS

Using mass spectrometry-based methods we have identified 231 proteins in the matrix of the S. purpuratus spicule matrix. Approximately two thirds of the identified proteins are either known or predicted to be extracellular proteins or transmembrane proteins with large ectodomains. The ectodomains may have been solubilized by partial proteolysis and subsequently integrated into the growing spicule. The most abundant protein of the spicule matrix is SM50. SM50-related proteins, SM30-related proteins, MSP130 and related proteins, matrix metalloproteases and carbonic anhydrase are among the most abundant components.

CONCLUSIONS

The spicule matrix is a relatively complex mixture of proteins not only containing matrix-specific proteins with a function in matrix assembly or mineralization, but also: 1) proteins possibly important for the formation of the continuous membrane delineating the mineralization space; 2) proteins for secretory processes delivering proteinaceous or non-proteinaceous precursors; 3) or proteins reflecting signaling events at the cell/matrix interface. Comparison of the proteomes of different skeletal matrices allows prediction of proteins of general importance for mineralization in sea urchins, such as SM50, SM30-E, SM29 or MSP130. The comparisons also help point out putative tissue-specific proteins, such as tooth phosphodontin or specific spicule matrix metalloproteases of the MMP18/19 group. Furthermore, the direct sequence analysis of peptides by MS/MS validates many predicted genes and confirms the existence of the corresponding proteins.

Cis-regulatory analysis of the sea urchin pigment cell gene polyketide synthase

The Strongylocentrotus purpuratus polyketide synthase gene (SpPks) encodes an enzyme required for the biosynthesis of the larval pigment echinochrome. SpPks is expressed exclusively in pigment cells and their precursors starting at blastula stage. The 7th-9th cleavage Delta-Notch signaling, required for pigment cell development, positively regulates SpPks. In previous studies, the transcription factors glial cell missing (SpGcm), SpGatae and kruppel-like (SpKrl/z13) have been shown to positively regulate SpPks. To uncover the structure of the Gene Regulatory Network (GRN) regulating the specification and differentiation processes of pigment cells, we experimentally analyzed the putative SpPks cis-regulatory region. We established that the -1.5kb region is sufficient to recapitulate the correct spatial and temporal expression of SpPks. Predicted DNA-binding sites for SpGcm, SpGataE and SpKrl are located within this region. The mutagenesis of these DNA-binding sites indicated that SpGcm, SpGataE and SpKrl are direct positive regulators of SpPks. These results demonstrate that the sea urchin GRN for pigment cell development is quite shallow, which is typical of type I embryo development.

Modularity and design principles in the sea urchin embryo gene regulatory network

The gene regulatory network (GRN) established experimentally for the pre-gastrular sea urchin embryo provides causal explanations of the biological functions required for spatial specification of embryonic regulatory states. Here we focus on the structure of the GRN which controls the progressive increase in complexity of territorial regulatory states during embryogenesis; and on the types of modular subcircuits of which the GRN is composed. Each of these subcircuit topologies executes a particular operation of spatial information processing. The GRN architecture reflects the particular mode of embryogenesis represented by sea urchin development. Network structure not only specifies the linkages constituting the genomic regulatory code for development, but also indicates the various regulatory requirements of regional developmental processes.

Parallel evolution of nacre building gene sets in molluscs

The capacity to biomineralize is closely linked to the rapid expansion of animal life during the early Cambrian, with many skeletonized phyla first appearing in the fossil record at this time. The appearance of disparate molluscan forms during this period leaves open the possibility that shells evolved independently and in parallel in at least some groups. To test this proposition and gain insight into the evolution of structural genes that contribute to shell fabrication, we compared genes expressed in nacre (mother-of-pearl) forming cells in the mantle of the bivalve Pinctada maxima and the gastropod Haliotis asinina. Despite both species having highly lustrous nacre, we find extensive differences in these expressed gene sets. Following the removal of housekeeping genes, less than 10% of all gene clusters are shared between these molluscs, with some being conserved biomineralization genes that are also found in deuterostomes. These differences extend to secreted proteins that may localize to the organic shell matrix, with less than 15% of this secretome being shared. Despite these differences, H. asinina and P. maxima both secrete proteins with repetitive low-complexity domains (RLCDs). Pinctada maxima RLCD proteins-for example, the shematrins-are predominated by silk/fibroin-like domains, which are absent from the H. asinina data set. Comparisons of shematrin genes across three species of Pinctada indicate that this gene family has undergone extensive divergent evolution within pearl oysters. We also detect fundamental bivalve-gastropod differences in extracellular matrix proteins involved in mollusc-shell formation. Pinctada maxima expresses a chitin synthase at high levels and several chitin deacetylation genes, whereas only one protein involved in chitin interactions is present in the H. asinina data set, suggesting that the organic matrix on which calcification proceeds differs fundamentally between these species. Large-scale differences in genes expressed in nacre-forming cells of Pinctada and Haliotis are compatible with the hypothesis that gastropod and bivalve nacre is the result of convergent evolution. The expression of novel biomineralizing RLCD proteins in each of these two molluscs and, interestingly, sea urchins suggests that the evolution of such structural proteins has occurred independently multiple times in the Metazoa.

A perturbation model of the gene regulatory network for oral and aboral ectoderm specification in the sea urchin embryo

The current gene regulatory network (GRN) for the sea urchin embryo pertains to pregastrular specification functions in the endomesodermal territories. Here we extend gene regulatory network analysis to the adjacent oral and aboral ectoderm territories over the same period. A large fraction of the regulatory genes predicted by the sea urchin genome project and shown in ancillary studies to be expressed in either oral or aboral ectoderm by 24 h are included, though universally expressed and pan-ectodermal regulatory genes are in general not. The loci of expression of these genes have been determined by whole mount in situ hybridization. We have carried out a global perturbation analysis in which expression of each gene was interrupted by introduction of morpholino antisense oligonucleotide, and the effects on all other genes were measured quantitatively, both by QPCR and by a new instrumental technology (NanoString Technologies nCounter Analysis System). At its current stage the network model, built in BioTapestry, includes 22 genes encoding transcription factors, 4 genes encoding known signaling ligands, and 3 genes that are yet unknown but are predicted to perform specific roles. Evidence emerged from the analysis pointing to distinctive subcircuit features observed earlier in other parts of the GRN, including a double negative transcriptional regulatory gate, and dynamic state lockdowns by feedback interactions. While much of the regulatory apparatus is downstream of Nodal signaling, as expected from previous observations, there are also cohorts of independently activated oral and aboral ectoderm regulatory genes, and we predict yet unidentified signaling interactions between oral and aboral territories.

Lessons from a gene regulatory network: echinoderm skeletogenesis provides insights into evolution, plasticity and morphogenesis

Significant new insights have emerged from the analysis of a gene regulatory network (GRN) that underlies the development of the endoskeleton of the sea urchin embryo. Comparative studies have revealed ways in which this GRN has been modified (and conserved) during echinoderm evolution, and point to mechanisms associated with the evolution of a new cell lineage. The skeletogenic GRN has also recently been used to study the long-standing problem of developmental plasticity. Other recent findings have linked this transcriptional GRN to morphoregulatory proteins that control skeletal anatomy. These new studies highlight powerful new ways in which GRNs can be used to dissect development and the evolution of morphogenesis.

Functional studies of regulatory genes in the sea urchin embryo

Sea urchin embryos are characterized by an extremely simple mode of development, rapid cleavage, high transparency, and well-defined cell lineage. Although they are not suitable for genetic studies, other approaches are successfully used to unravel mechanisms and molecules involved in cell fate specification and morphogenesis. Microinjection is the elective method to study gene function in sea urchin embryos. It is used to deliver precise amounts of DNA, RNA, oligonucleotides, peptides, or antibodies into the eggs or even into blastomeres. Here we describe microinjection as it is currently applied in our laboratory and show how it has been used in gene perturbation analyses and dissection of cis-regulatory DNA elements.

Evolution of gene regulatory network architectures: examples of subcircuit conservation and plasticity between classes of echinoderms

Developmental gene regulatory networks (GRNs) explain how regulatory states are established in particular cells during development and how these states then determine the final form of the embryo. Evolutionary changes to the sequence of the genome will direct reorganization of GRN architectures, which in turn will lead to the alteration of developmental programs. A comparison of GRN architectures must consequently reveal the molecular basis for the evolution of developmental programs among different organisms. This review highlights some of the important findings that have emerged from the most extensive direct comparison of GRN architectures to date. Comparison of the orthologous GRNs for endomesodermal specification in the sea urchin and sea star, provides examples of several discrete, functional GRN subcircuits and shows that they are subject to diverse selective pressures. This demonstrates that different regulatory linkages may be more or less amenable to evolutionary change. One of the more surprising findings from this comparison is that GRN-level functions may be maintained while the factors performing the functions have changed, suggesting that GRNs have a high capacity for compensatory changes involving transcription factor binding to cis regulatory modules.

Modeling the dynamics of transcriptional gene regulatory networks for animal development

The dynamic process of cell fate specification is regulated by networks of regulatory genes. The architecture of the network defines the temporal order of specification events. To understand the dynamic control of the developmental process, the kinetics of mRNA and protein synthesis and the response of the cis-regulatory modules to transcription factor concentration must be considered. Here we review mathematical models for mRNA and protein synthesis kinetics which are based on experimental measurements of the rates of the relevant processes. The model comprises the response functions of cis-regulatory modules to their transcription factor inputs, by incorporating binding site occupancy and its dependence on biologically measurable quantities. We use this model to simulate gene expression, to distinguish between cis-regulatory execution of "AND" and "OR" logic functions, rationalize the oscillatory behavior of certain transcriptional auto-repressors and to show how linked subcircuits can be dealt with. Model simulations display the effects of mutation of binding sites, or perturbation of upstream gene expression. The model is a generally useful tool for understanding gene regulation and the dynamics of cell fate specification.

Compositional genome contexts affect gene expression control in sea urchin embryo

Gene expression is widely perceived as exclusively controlled by the information contained in cis-regulatory regions. These are built in a modular way, each module being a cluster of binding sites for the transcription factors that control the level, the location and the time at which gene transcription takes place. On the other hand, results from our laboratory have shown that gene expression is affected by the compositional properties (GC levels) of the isochores in which genes are embedded, i.e. the genome context. To clarify how compositional genomic properties affect the way cis-regulatory information is utilized, we have changed the genome context of a GFP-reporter gene containing the complete cis-regulatory region of the gene spdeadringer (spdri), expressed during sea urchin embryogenesis. We have observed that GC levels higher or lower than those found in the natural genome context can alter the reporter expression pattern. We explain this as the result of an interference with the functionality of specific modules in the gene's cis-regulatory region. From these observations we derive the notion that the compositional properties of the genome context can affect cis-regulatory control of gene expression. Therefore although the way a gene works depends on the information contained in its cis-regulatory region, availability of such information depends on the compositional properties of the genomic context.

In-depth, high-accuracy proteomics of sea urchin tooth organic matrix

BACKGROUND

The organic matrix contained in biominerals plays an important role in regulating mineralization and in determining biomineral properties. However, most components of biomineral matrices remain unknown at present. In sea urchin tooth, which is an important model for developmental biology and biomineralization, only few matrix components have been identified. The recent publication of the Strongylocentrotus purpuratus genome sequence rendered possible not only the identification of genes potentially coding for matrix proteins, but also the direct identification of proteins contained in matrices of skeletal elements by in-depth, high-accuracy proteomic analysis.

RESULTS

We identified 138 proteins in the matrix of tooth powder. Only 56 of these proteins were previously identified in the matrices of test (shell) and spine. Among the novel components was an interesting group of five proteins containing alanine- and proline-rich neutral or basic motifs separated by acidic glycine-rich motifs. In addition, four of the five proteins contained either one or two predicted Kazal protease inhibitor domains. The major components of tooth matrix were however largely identical to the set of spicule matrix proteins and MSP130-related proteins identified in test (shell) and spine matrix. Comparison of the matrices of crushed teeth to intact teeth revealed a marked dilution of known intracrystalline matrix proteins and a concomitant increase in some intracellular proteins.

CONCLUSION

This report presents the most comprehensive list of sea urchin tooth matrix proteins available at present. The complex mixture of proteins identified may reflect many different aspects of the mineralization process. A comparison between intact tooth matrix, presumably containing odontoblast remnants, and crushed tooth matrix served to differentiate between matrix components and possible contributions of cellular remnants. Because LC-MS/MS-based methods directly measures peptides our results validate many predicted genes and confirm the existence of the corresponding proteins. Knowledge of the components of this model system may stimulate further experiments aiming at the elucidation of structure, function, and interaction of biomineral matrix components.

A protocol describing the principles of cis-regulatory analysis in the sea urchin

cis-Regulatory analysis (CRA) is the precise identification of the cis-acting genomic sequences regulating gene transcription. As such, CRA provides essential mechanistic insight into key biological processes such as development. The first phase of this protocol involves identification of a large (approximately 100 kb) clone of genomic material surrounding the gene of interest and use of this clone to establish a reliable and unambiguous reporter assay. In the second phase, phylogenetic footprinting is used to identify candidate regulatory modules; these genomic sequences are then recursively tested for reporter activity. In the final phase, potential transcription factor binding sites are identified and disrupted in reporter constructs for individual testing. The strengths of this method reflect the use of large clones containing all relevant genomic regulatory sequences to establish a reporter assay with high fidelity. Given these foundational elements, predicted or suspected regulatory inputs can be rigorously tested and novel regulatory inputs identified. Although the expected time line varies greatly with the depth of information required, results may be obtained in as little as 4-6 months, but more detailed analyses will require several years to complete.

The surprising complexity of the transcriptional regulation of the spdri gene reveals the existence of new linkages inside sea urchin's PMC and Oral Ectoderm Gene Regulatory Networks

During sea urchin embryogenesis the spdri gene participates in two separate Gene Regulatory Networks (GRNs): the Primary Mesenchyme Cells' (PMCs) and the Oral Ectoderm's one. In both cases, activation of the gene follows initial specification events [Amore, G., Yavrouian, R., Peterson, K., Ransick, A., McClay, D., Davidson, E., 2003. Spdeadringer, a sea urchin embryo gene required separately in skeletogenic and oral ectoderm gene regulatory networks. Dev. Biol. 261, 55-81.]. We identified a portion of genomic DNA ("4.7IL" -3456;+389) which is sufficient to replicate sdpri's expression pattern in experiments of transgenesis, using a GFP reporter. In our experiments, the activation kinetic of 4.7IL-GFP was similar to that of the endogenous gene and the reporter responded to known spdri's transcriptional regulators (Ets1, Alx1, Gsc and Dri). Here we present a dissection of this regulatory region and a description of the modules involved in spdri's transcriptional regulation. Both in the PMCs' and Oral Ectoderm's expression phases, activation of spdri is obtained through the integration of three kinds of inputs: positive and globally distributed ones; negative ones (that prevent ectopic expression); positive and tissue-specific ones. Our results allow to expand the map of the regulatory connections at the spdri node, both in the PMCs and in the Oral Ectoderm Gene Regulatory Networks (GRNs).

cis-Regulatory sequences driving the expression of the Hbox12 homeobox-containing gene in the presumptive aboral ectoderm territory of the Paracentrotus lividus sea urchin embryo

Embryonic development is coordinated by networks of evolutionary conserved regulatory genes encoding transcription factors and components of cell signalling pathways. In the sea urchin embryo, a number of genes encoding transcription factors display territorial restricted expression. Among these, the zygotic Hbox12 homeobox gene is transiently transcribed in a limited number of cells of the animal-lateral half of the early Paracentrotus lividus embryo, whose descendants will constitute part of the ectoderm territory. To obtain insights on the regulation of Hbox12 expression, we have explored the cis-regulatory apparatus of the gene. In this paper, we show that the intergenic region of the tandem Hbox12 repeats drives GFP expression in the presumptive aboral ectoderm and that a 234 bp fragment, defined aboral ectoderm (AE) module, accounts for the restricted expression of the transgene. Within this module, a consensus sequence for a Sox factor and the binding of the Otx activator are both required for correct Hbox12 gene expression. Spatial restriction to the aboral ectoderm is achieved by a combination of different repressive sequence elements. Negative sequence elements necessary for repression in the endomesoderm map within the most upstream 60 bp region and nearby the Sox binding site. Strikingly, a Myb-like consensus is necessary for repression in the oral ectoderm, while down-regulation at the gastrula stage depends on a GA-rich region. These results suggest a role for Hbox12 in aboral ectoderm specification and represent our first attempt in the identification of the gene regulatory circuits involved in this process.

Global regulatory logic for specification of an embryonic cell lineage

Explanation of a process of development must ultimately be couched in the terms of the genomic regulatory code. Specification of an embryonic cell lineage is driven by a network of interactions among genes encoding transcription factors. Here, we present the gene regulatory network (GRN) that directs the specification of the skeletogenic micromere lineage of the sea urchin embryo. The GRN now includes all regulatory genes expressed in this lineage up to late blastula stage, as identified in a genomewide survey. The architecture of the GRN was established by a large-scale perturbation analysis in which the expression of each gene in the GRN was cut off by use of morpholinos, and the effects on all other genes were measured quantitatively. Several cis-regulatory analyses provided additional evidence. The explanatory power of the GRN suffices to provide a causal explanation for all observable developmental functions of the micromere lineage during the specification period. These functions are: (i) initial acquisition of identity through transcriptional interpretation of localized maternal cues; (ii) activation of specific regulatory genes by use of a double negative gate; (iii) dynamic stabilization of the regulatory state by activation of a feedback subcircuit; (iv) exclusion of alternative regulatory states; (v) presentation of a signal required by the micromeres themselves and of two different signals required for development of adjacent endomesodermal lineages; and (vi) lineage-specific activation of batteries of skeletogenic genes. The GRN precisely predicts gene expression responses and provides a coherent explanation of the biology of specification.

A spatially dynamic cohort of regulatory genes in the endomesodermal gene network of the sea urchin embryo

A gene regulatory network subcircuit comprising the otx, wnt8, and blimp1 genes accounts for a moving torus of gene expression that sweeps concentrically across the vegetal domain of the sea urchin embryo. Here we confirm by mutation the inputs into the blimp1cis-regulatory module predicted by network analysis. Its essential design feature is that it includes both activation and autorepression sites. The wnt8 gene is functionally linked into the subcircuit in that cells receiving this ligand generate a beta-catenin/Tcf input required for blimp1 expression, while the wnt8 gene in turn requires a Blimp1 input. Their torus-like spatial expression patterns and gene regulatory analysis indicate that the genes even-skipped and hox11/13b are also entrained by this subcircuit. We verify the cis-regulatory inputs of even-skipped predicted by network analysis. These include activation by beta-catenin/Tcf and Blimp1, repression within the torus by Hox11/13b, and repression outside the torus by Tcf in the absence of Wnt8 signal input. Thus even-skipped and hox11/13b, along with blimp1 and wnt8, are members of a cohort of torus genes with similar regulatory inputs and similar, though slightly out-of-phase, expression patterns, which reflect differences in cis-regulatory design.

Insights from amphioxus into the evolution of vertebrate cartilage

Central to the story of vertebrate evolution is the origin of the vertebrate head, a problem difficult to approach using paleontology and comparative morphology due to a lack of unambiguous intermediate forms. Embryologically, much of the vertebrate head is derived from two ectodermal tissues, the neural crest and cranial placodes. Recent work in protochordates suggests the first chordates possessed migratory neural tube cells with some features of neural crest cells. However, it is unclear how and when these cells acquired the ability to form cellular cartilage, a cell type unique to vertebrates. It has been variously proposed that the neural crest acquired chondrogenic ability by recruiting proto-chondrogenic gene programs deployed in the neural tube, pharynx, and notochord. To test these hypotheses we examined the expression of 11 amphioxus orthologs of genes involved in neural crest chondrogenesis. Consistent with cellular cartilage as a vertebrate novelty, we find that no single amphioxus tissue co-expresses all or most of these genes. However, most are variously co-expressed in mesodermal derivatives. Our results suggest that neural crest-derived cartilage evolved by serial cooption of genes which functioned primitively in mesoderm.

Paleogenomics of echinoderms

Paleogenomics propels the meaning of genomic studies back through hundreds of millions of years of deep time. Now that the genome of the echinoid Strongylocentrotus purpuratus is sequenced, the operation of its genes can be interpreted in light of the well-understood echinoderm fossil record. Characters that first appear in Early Cambrian forms are still characteristic of echinoderms today. Key genes for one of these characters, the biomineralized tissue stereom, can be identified in the S. purpuratus genome and are likely to be the same genes that were involved with stereom formation in the earliest echinoderms some 520 million years ago.

Deciphering the Underlying Mechanism of Specification and Differentiation: The Sea Urchin Gene Regulatory Network

The regulatory genome functions as a vast information processor through development. It processes the initial conditions that are set by asymmetric distributions of cellular components in the egg and translates them into the onset of spatially localized specification states. It regulates the timely differential activation of signaling molecules and transcription factors that divide the emerging domains into subdomains. It also governs the activation of groups of differentiation genes, the genes that encode, at the protein level, the functional and the structural properties of a cell type. The sea urchin endomesoderm gene regulatory network provides a window into the different levels of the regulatory apparatus. It demonstrates how the static physical genomic components define functional connections between the various regulatory genes that act together to conduct the dynamical developmental program.

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

Biomineralization, the biologically controlled formation of mineral deposits, is of widespread importance in biology, medicine, and engineering. Mineralized structures are found in most metazoan phyla and often have supportive, protective, or feeding functions. Among deuterostomes, only echinoderms and vertebrates produce extensive biomineralized structures. Although skeletons appeared independently in these two groups, ancestors of the vertebrates and echinoderms may have utilized similar components of a shared genetic "toolkit" to carry out biomineralization. The present study had two goals. First, we sought to expand our understanding of the proteins involved in biomineralization in the sea urchin, a powerful model system for analyzing the basic cellular and molecular mechanisms that underlie this process. Second, we sought to shed light on the possible evolutionary relationships between biomineralization in echinoderms and vertebrates. We used several computational methods to survey the genome of the purple sea urchin Strongylocentrotus purpuratus for gene products involved in biomineralization. Our analysis has greatly expanded the collection of biomineralization-related proteins. We have found that these proteins are often members of small families encoded by genes that are clustered in the genome. Most of the proteins are sea urchin-specific; that is, they have no apparent homologues in other invertebrate deuterostomes or vertebrates. Similarly, many of the vertebrate proteins that mediate mineral deposition do not have counterparts in the S. purpuratus genome. Our findings therefore reveal substantial differences in the primary sequences of proteins that mediate biomineral formation in echinoderms and vertebrates, possibly reflecting loose constraints on the primary structures of the proteins involved. On the other hand, certain cellular and molecular processes associated with earlier events in skeletogenesis appear similar in echinoderms and vertebrates, leaving open the possibility of deeper evolutionary relationships.