High regulatory gene use in sea urchin embryogenesis: Implications for bilaterian development and evolution

Widespread priming of transcriptional regulatory elements by incipient accessibility or RNA polymerase II pause in early embryos of the sea urchin Strongylocentrotus purpuratus

Transcriptional regulatory elements (TREs) are the primary nodes that control developmental gene regulatory networks. In embryo stages, larvae, and adult differentiated red spherule cells of the sea urchin Strongylocentrotus purpuratus, transcriptionally engaged TREs are detected by Precision Run-On Sequencing (PRO-seq), which maps genome-wide at base pair resolution the location of paused or elongating RNA polymerase II (Pol II). In parallel, TRE accessibility is estimated by the Assay for Transposase-Accessible Chromatin using Sequencing (ATAC-seq). Our analysis identifies surprisingly early and widespread TRE accessibility in 4-cell cleavage embryos that is not necessarily followed by concurrent or subsequent transcription. TRE transcriptional differences identified by PRO-seq provide more contrast among embryonic stages than ATAC-seq accessibility differences, in agreement with the apparent excess of accessible but inactive TREs during embryogenesis. Global TRE accessibility reaches a maximum around the 20-hour late blastula stage, which coincides with the consolidation of major embryo regionalizations and peak histone variant H2A.Z expression. A transcriptional potency model based on labile nucleosome TRE occupancy driven by DNA sequences and the prevalence of histone variants is proposed in order to explain the basal accessibility of transcriptionally inactive TREs during embryogenesis. However, our results would not reconcile well with labile nucleosome models based on simple A/T sequence enrichment. In addition, a large number of distal TREs become transcriptionally disengaged during developmental progression, in support of an early Pol II paused model for developmental gene regulation that eventually resolves in transcriptional activation or silencing. Thus, developmental potency in early embryos may be facilitated by incipient accessibility and transcriptional pause at TREs.

The temporal dynamics of the sea urchin regulome

In this work, we used Nanostring N-counter technology, to evaluate the mRNA expression level of more than 330 regulatory genes over 34 time points covering the first three days of development of the sea urchin larvae. The hierarchical clustering of the mRNAs expression levels has identified groups corresponding to the major developmental landmarks (e.g. maternal to zygotic transition and gastrulation). Furthermore, comparison with previous experiments indicates high reproducibility of mRNA level temporal dynamics across batches. Finally, we generated an online tool to visualise gene expression during sea urchin larval development. The site can be accessed at and https://www621.lamp.le.ac.uk/nanostring_app/nanostring/.

Summary: This study uses Nanostring N-counter technology to evaluate the mRNA expression level of more than 330 regulatory genes over 34 time points covering the first three days of development of the sea urchin larvae.

The developmental transcriptome for Lytechinus variegatus exhibits temporally punctuated gene expression changes

Embryonic development is arguably the most complex process an organism undergoes during its lifetime, and understanding this complexity is best approached with a systems-level perspective. The sea urchin has become a highly valuable model organism for understanding developmental specification, morphogenesis, and evolution. As a non-chordate deuterostome, the sea urchin occupies an important evolutionary niche between protostomes and vertebrates. Lytechinus variegatus (Lv) is an Atlantic species that has been well studied, and which has provided important insights into signal transduction, patterning, and morphogenetic changes during embryonic and larval development. The Pacific species, Strongylocentrotus purpuratus (Sp), is another well-studied sea urchin, particularly for gene regulatory networks (GRNs) and cis-regulatory analyses. A well-annotated genome and transcriptome for Sp are available, but similar resources have not been developed for Lv. Here, we provide an analysis of the Lv transcriptome at 11 timepoints during embryonic and larval development. Temporal analysis suggests that the gene regulatory networks that underlie specification are well-conserved among sea urchin species. We show that the major transitions in variation of embryonic transcription divide the developmental time series into four distinct, temporally sequential phases. Our work shows that sea urchin development occurs via sequential intervals of relatively stable gene expression states that are punctuated by abrupt transitions.

A Survey on Tubulin and Arginine Methyltransferase Families Sheds Light on P. lividus Embryo as Model System for Antiproliferative Drug Development

Tubulins and microtubules (MTs) represent targets for taxane-based chemotherapy. To date, several lines of evidence suggest that effectiveness of compounds binding tubulin often relies on different post-translational modifications on tubulins. Among them, methylation was recently associated to drug resistance mechanisms impairing taxanes binding. The sea urchin is recognized as a research model in several fields including fertilization, embryo development and toxicology. To date, some α- and β-tubulin genes have been identified in P. lividus, while no data are available in echinoderms for arginine methyl transferases (PRMT). To evaluate the exploiting of the sea urchin embryo in the field of antiproliferative drug development, we carried out a survey of the expressed α- and β-tubulin gene sets, together with a comprehensive analysis of the PRMT gene family and of the methylable arginine residues in P. lividus tubulins. Because of their specificities, the sea urchin embryo may represent an interesting tool for dissecting mechanisms of tubulin targeting drug action. Therefore, results herein reported provide evidences supporting the P. lividus embryo as animal system for testing antiproliferative drugs.

Methods for the experimental and computational analysis of gene regulatory networks in sea urchins

The discovery of gene regulatory networks (GRNs) has opened a gate to access the genomic mechanisms controlling development. GRNs are systems of transcriptional regulatory circuits that control the differential specification of cell fates during development by regulating gene expression. The experimental analysis of GRNs involves a collection of methods, each revealing aspects of the overall control process. This review provides an overview of experimental and computational methods that have been successfully applied for solving developmental GRNs in the sea urchin embryo. The key in this approach is to obtain experimental evidence for functional interactions between transcription factors and regulatory DNA. In the second part of this review, a more generally applicable strategy is discussed that shows a path from experimental evidence to annotation of regulatory linkages to the generation of GRN models.

Neurogenesis in the sea urchin embryo is initiated uniquely in three domains

Many marine larvae begin feeding within a day of fertilization, thus requiring rapid development of a nervous system to coordinate feeding activities. Here, we examine the patterning and specification of early neurogenesis in sea urchin embryos. Lineage analysis indicates that neurons arise locally in three regions of the embryo. Perturbation analyses showed that when patterning is disrupted, neurogenesis in the three regions is differentially affected, indicating distinct patterning requirements for each neural domain. Six transcription factors that function during proneural specification were identified and studied in detail. Perturbations of these proneural transcription factors showed that specification occurs differently in each neural domain prior to the Delta-Notch restriction signal. Though gene regulatory network state changes beyond the proneural restriction are largely unresolved, the data here show that the three neural regions already differ from each other significantly early in specification. Future studies that define the larval nervous system in the sea urchin must therefore separately characterize the three populations of neurons that enable the larva to feed, to navigate, and to move food particles through the gut.

Echinobase: an expanding resource for echinoderm genomic information

Echinobase, a web accessible information system of diverse genomics and biological data for the echinoderm clade, grew out of SpBase, the first echinoderm genome project for sea urchin, Strongylocentrotus purpuratus. Sea urchins and their relatives are utilitarian research models in fields ranging from marine biology to developmental biology and gene regulatory systems. Echinobase is a user-friendly web interface that links an array of biological data that would otherwise have been tedious and frustrating for researchers to extract and organize. The system hosts a powerful gene search engine, genomics browser and other bioinformatics tools to investigate genomics and high throughput data. The Echinobase information system now serves genomic information for eight echinoderm species: S. purpuratus, Strongylocentrotus fransciscanus, Allocentrotus fragilis, Lytechinus variegatus, Patiria miniata, Parastichopus parvimensis and Ophiothrix spiculata, Eucidaris tribuloides. Herein lies a description of the web information system, genomics data types and content hosted by Echinobase.org. The goal of Echinobase is to connect genomic information to various experimental data and accelerate the research in field of molecular biology, developmental process, gene regulatory networks and more recently engineering biological systems0.

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

Regulatory states in the developmental control of gene expression

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

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

Using Morpholinos to Probe Gene Networks in Sea Urchin

The control processes that underlie the progression of development can be summarized in maps of gene regulatory networks (GRNs). A critical step in their assembly is the systematic perturbation of network candidates. In sea urchins the most important method for interfering with expression in a gene-specific way is application of morpholino antisense oligonucleotides (MOs). MOs act by binding to their sequence complement in transcripts resulting in a block in translation or a change in splicing and thus result in a loss of function. Despite the tremendous success of this technology, recent comparisons to mutants generated by genome editing have led to renewed criticism and challenged its reliability. As with all methods based on sequence recognition, MOs are prone to off-target binding that may result in phenotypes that are erroneously ascribed to the loss of the intended target. However, the slow progression of development in sea urchins has enabled extremely detailed studies of gene activity in the embryo. This wealth of knowledge paired with the simplicity of the sea urchin embryo enables careful analysis of MO phenotypes through a variety of methods that do not rely on terminal phenotypes. This article summarizes the use of MOs in probing GRNs and the steps that should be taken to assure their specificity.

Embryonic neurogenesis in echinoderms

The phylogenetic position of echinoderms is well suited to revealing shared features of deuterostomes that distinguish them from other bilaterians. Although echinoderm neurobiology remains understudied, genomic resources, molecular methods, and systems approaches have enabled progress in understanding mechanisms of embryonic neurogenesis. Even though the morphology of echinoderm larvae is diverse, larval nervous systems, which arise during gastrulation, have numerous similarities in their organization. Diverse neural subtypes and specialized sensory neurons have been identified and details of neuroanatomy using neuron-specific labels provide hypotheses for neural function. The early patterning of ectoderm and specification of axes has been well studied in several species and underlying gene regulatory networks have been established. The cells giving rise to central and peripheral neural components have been identified in urchins and sea stars. Neurogenesis includes typical metazoan features of asymmetric division of neural progenitors and in some cases limited proliferation of neural precursors. Delta/Notch signaling has been identified as having critical roles in regulating neural patterning and differentiation. Several transcription factors functioning in pro-neural phases of specification, neural differentiation, and sub-type specification have been identified and structural or functional components of neurons are used as differentiation markers. Several methods for altering expression in embryos have revealed aspects of a regulatory hierarchy of transcription factors in neurogenesis. Interfacing neurogenic gene regulatory networks to the networks regulating ectodermal domains and identifying the spatial and temporal inputs that pattern the larval nervous system is a major challenge that will contribute substantially to our understanding of the evolution of metazoan nervous systems. This article is categorized under: Comparative Development and Evolution > Model Systems Comparative Development and Evolution > Body Plan Evolution Early Embryonic Development > Gastrulation and Neurulation.

EchinoBase: Tools for Echinoderm Genome Analyses

The echinoderms are a phylum of invertebrate deuterostome animals that constitute important research models for a number of biological disciplines. EchinoBase ( www.echinobase.org ) is an echinoderm-specific genome database and web information system that provides a platform for the interrogation and exploration of echinoderm genomic data. This chapter outlines the datasets available on EchinoBase; from assembled genomes and genome annotations, to spatial and quantitative expression data, as well as functional genomics datasets. We also highlight the bioinformatic tools available on the website to facilitate rapid inquiries using these data (genome browsers, precompiled BLAST databases, etc.), and suggest optimized strategies for performing these inquiries. We conclude with a perspective on how one could integrate various genomic resources to predict putative noncoding regulatory regions. The available datasets and analyses they permit provide the basic components required for developing an understanding of how echinoderm genomes are regulated, especially during early development, and provides a platform for comparative genomic inquiries among species in this phylum.

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.

Quantitative developmental transcriptomes of the sea urchin Strongylocentrotus purpuratus

Development depends on the precise control of gene expression in time and space. A critical step towards understanding the global gene regulatory networks underlying development is to obtain comprehensive information on gene expression. In this study, we measured expression profiles for the entire expressed gene set during sea urchin embryonic development. We confirmed the reliability of these profiles by comparison with NanoString measurements for a subset of genes and with literature values. The data show that ~16,500 genes have been activated by the end of embryogenesis, and for half of them the transcript abundance changes more than 10-fold during development. From this genome scale expression survey, we show that complex patterns of expression by many genes underlie embryonic development, particularly during the early stages before gastrulation. An intuitive web application for data query and visualization is presented to facilitate use of this large dataset.

Phylogenomics of strongylocentrotid sea urchins

Background

Strongylocentrotid sea urchins have a long tradition as model organisms for studying many fundamental processes in biology including fertilization, embryology, development and genome regulation but the phylogenetic relationships of the group remain largely unresolved. Although the differing isolating mechanisms of vicariance and rapidly evolving gamete recognition proteins have been proposed, a stable and robust phylogeny is unavailable.

Results

We used a phylogenomic approach with mitochondrial and nuclear genes taking advantage of the whole-genome sequencing of nine species in the group to establish a stable (i.e. concordance in tree topology among multiple lies of evidence) and robust (i.e. high nodal support) phylogenetic hypothesis for the family Strongylocentrotidae. We generated eight draft mitochondrial genome assemblies and obtained 13 complete mitochondrial genes for each species. Consistent with previous studies, mitochondrial sequences failed to provide a reliable phylogeny. In contrast, we obtained a very well-supported phylogeny from 2301 nuclear genes without evidence of positive Darwinian selection both from the majority of most-likely gene trees and the concatenated fourfold degenerate sites: ((P. depressus, (M. nudus, M. franciscanus), (H. pulcherrimus, (S. purpuratus, (S. fragilis, (S. pallidus, (S. droebachiensis, S. intermedius)). This phylogeny was consistent with a single invasion of deep-water environments followed by a holarctic expansion by Strongylocentrotus. Divergence times for each species estimated with reference to the divergence times between the two major clades of the group suggest a correspondence in the timing with the opening of the Bering Strait and the invasion of the holarctic regions.

Conclusions

Nuclear genome data contains phylogenetic signal informative for understanding the evolutionary history of this group. However, mitochondrial genome data does not. Vicariance can explain major patterns observed in the phylogeny. Other isolating mechanisms are appropriate to explore in this system to help explain divergence patterns not well supported by vicariance, such as the effects of rapidly evolving gamete recognition proteins on isolating populations. Our findings of a stable and robust phylogeny, with the increase in mitochondrial and nuclear comparative genomic data, provide a system in which we can enhance our understanding of molecular evolution and adaptation in this group of sea urchins.

Gene structure in the sea urchin Strongylocentrotus purpuratus based on transcriptome analysis

A comprehensive transcriptome analysis has been performed on protein-coding RNAs of Strongylocentrotus purpuratus, including 10 different embryonic stages, six feeding larval and metamorphosed juvenile stages, and six adult tissues. In this study, we pooled the transcriptomes from all of these sources and focused on the insights they provide for gene structure in the genome of this recently sequenced model system. The genome had initially been annotated by use of computational gene model prediction algorithms. A large fraction of these predicted genes were recovered in the transcriptome when the reads were mapped to the genome and appropriately filtered and analyzed. However, in a manually curated subset, we discovered that more than half the computational gene model predictions were imperfect, containing errors such as missing exons, prediction of nonexistent exons, erroneous intron/exon boundaries, fusion of adjacent genes, and prediction of multiple genes from single genes. The transcriptome data have been used to provide a systematic upgrade of the gene model predictions throughout the genome, very greatly improving the research usability of the genomic sequence. We have constructed new public databases that incorporate information from the transcriptome analyses. The transcript-based gene model data were used to define average structural parameters for S. purpuratus protein-coding genes. In addition, we constructed a custom sea urchin gene ontology, and assigned about 7000 different annotated transcripts to 24 functional classes. Strong correlations became evident between given functional ontology classes and structural properties, including gene size, exon number, and exon and intron size.

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 cis-regulatory genomics for systems biology

Gene expression is controlled by interactions between trans-regulatory factors and cis-regulatory DNA sequences, and these interactions constitute the essential functional linkages of gene regulatory networks (GRNs). Validation of GRN models requires experimental cis-regulatory tests of predicted linkages to authenticate their identities and proposed functions. However, cis-regulatory analysis is, at present, at a severe bottleneck in genomic system biology because of the demanding experimental methodologies currently in use for discovering cis-regulatory modules (CRMs), in the genome, and for measuring their activities. Here we demonstrate a high-throughput approach to both discovery and quantitative characterization of CRMs. The unique aspect is use of DNA sequence tags to "barcode" CRM expression constructs, which can then be mixed, injected together into sea urchin eggs, and subsequently deconvolved. This method has increased the rate of cis-regulatory analysis by >100-fold compared with conventional one-by-one reporter assays. The utility of the DNA-tag reporters was demonstrated by the rapid discovery of 81 active CRMs from 37 previously unexplored sea urchin genes. We then obtained simultaneous high-resolution temporal characterization of the regulatory activities of more than 80 CRMs. On average 2-3 CRMs were discovered per gene. Comparison of endogenous gene expression profiles with those of the CRMs recovered from each gene showed that, for most cases, at least one CRM is active in each phase of endogenous expression, suggesting that CRM recovery was comprehensive. This approach will qualitatively alter the practice of GRN construction as well as validation, and will impact many additional areas of regulatory system biology.

Echinoderms: potential model systems for studies on muscle regeneration

Organisms of the phylum Echinodermata show some of the most impressive regenerative feats within the animal kingdom. Following injury or self-induced autotomy, species in this phylum can regenerate most tissues and organs, being the regeneration of the muscular systems one of the best studied. Even though echinoderms are closely related to chordates, they are little known in the biomedical field, and therefore their uses to study pharmacological effects on muscle formation and/or regeneration have been extremely limited. In order to rectify this lack of knowledge, we describe here the echinoderm muscular systems, particularly the somatic and visceral muscle components. In addition, we provide details of the processes that are known to take place during muscle regeneration, namely dedifferentiation, myogenesis and new muscle formation. Finally, we provide the available information on molecular and pharmacological studies that involve echinoderm muscle regeneration. We expect that by making this information accessible, researchers consider the use of echinoderms as model systems for pharmacological studies in muscle development and regeneration.

The endoderm gene regulatory network in sea urchin embryos up to mid-blastula stage

As the result of early specification processes, sea urchin embryos eventually form various mesodermal cell lineages and a gut consisting of fore-, mid- and hindgut. The progression of specification as well as the overall spatial organization of the organism is encoded in its gene regulatory networks (GRNs). We have analyzed the GRN driving endoderm specification up to the onset of gastrulation and present in this paper the mechanisms which determine this process up to mid-blastula stage. At this stage, the embryo consists of two separate lineages of endoderm precursor cells with distinct regulatory states. One of these lineages, the veg2 cell lineage, gives rise to endoderm and mesoderm cell types. The separation of these cell fates is initiated by the spatially confined activation of the mesoderm GRN superimposed on a generally activated endoderm GRN within veg2 descendants. Here we integrate the architecture of regulatory interactions with the spatial restriction of regulatory gene expression to model the logic control of endoderm development.

Experimental advances in the characterization of metazoan gene regulatory networks

Gene regulatory networks (GRNs) play a vital role in metazoan development and function, and deregulation of these networks is often implicated in disease. GRNs depict the dynamic interactions between genomic and regulatory state components. The genomic components comprise genes and their associated cis-regulatory elements. The regulatory state components consist primarily of transcriptional complexes that bind the latter elements. With the availability of complete genome sequences, several approaches have recently been developed which promise to significantly enhance our ability to identify either the genomic or regulatory state components, or the interactions between these two. In this review, I provide an in-depth overview of these approaches and detail how each contributes to a more comprehensive understanding of GRN composition and function.

The gene regulatory network basis of the "community effect," and analysis of a sea urchin embryo example

The "Community Effect" denotes intra-territorial signaling amongst cells which constitute a particular tissue or embryonic progenitor field. The cells of the territory express the same transcriptional regulatory state, and the intra-territorial signaling is essential to maintenance of this specific regulatory state. The structure of the underlying gene regulatory network (GRN) subcircuitry explains the genomically wired mechanism by which community effect signaling is linked to the continuing transcriptional generation of the territorial regulatory state. A clear example is afforded by the oral ectoderm GRN of the sea urchin embryo where cis-regulatory evidence, experimental embryology, and network analysis combine to provide a complete picture. We review this example and consider less well known but similar cases in other developing systems where the same subcircuit GRN topology is present. To resolve mechanistic issues that arise in considering how community effect signaling could operate to produce its observed effects, we construct and analyze the behavior of a quantitative model of community effect signaling in the sea urchin embryo oral ectoderm. Community effect network topology could constitute part of the genomic regulatory code that defines transcriptional function in multicellular tissues composed of cells in contact, and hence may have arisen as a metazoan developmental strategy.

SpBase: the sea urchin genome database and web site

SpBase is a system of databases focused on the genomic information from sea urchins and related echinoderms. It is exposed to the public through a web site served with open source software (http://spbase.org/). The enterprise was undertaken to provide an easily used collection of information to directly support experimental work on these useful research models in cell and developmental biology. The information served from the databases emerges from the draft genomic sequence of the purple sea urchin, Strongylocentrotus purpuratus and includes sequence data and genomic resource descriptions for other members of the echinoderm clade which in total span 540 million years of evolutionary time. This version of the system contains two assemblies of the purple sea urchin genome, associated expressed sequences, gene annotations and accessory resources. Search mechanisms for the sequences and the gene annotations are provided. Because the system is maintained along with the Sea Urchin Genome resource, a database of sequenced clones is also provided.

The sea urchin genome as a window on function

The emphasis on the sequencing of genomes seems to make this task an end in itself. However, genome sequences and the genes that are predicted from them are really an opportunity to examine the biological function of the organism constructed by that genome. This point is illustrated here by examples in which the newly annotated gene complement reveals surprises about the way Strongylocentrotus purpuratus, the purple sea urchin, goes about its business. The three topics considered here are the nature of the innate immune system; the unexpected complexity of sensory function implied by genes encoding sensory proteins; and the remarkable intricacy of the regulatory gene complement in embryogenesis.

A new method, using cis-regulatory control, for blocking embryonic gene expression

Many genes, and particularly regulatory genes, are utilized multiple times in unrelated phases of development. For studies of gene function during embryogenesis, there is often need of a method for interfering with expression only at a specific developmental time or place. Here we show that in sea urchin embryos cis-regulatory control systems which operate only at specific times and places can be used to drive expression of short designed sequences targeting given primary transcripts, thereby effectively taking out the function of the target genes. The active sequences are designed to be complementary to intronic sequences of the primary transcript of the target genes. In this work, the target genes were the transcription factors alx1 and ets1, both required for skeletogenesis, and the regulatory drivers were from the sm30 and tbr genes. The sm30 gene is expressed only after skeletogenic cell ingression. When its regulatory apparatus was used as driver, the alx1 and ets1 repression constructs had the effect of preventing postgastrular skeletogenesis, while not interfering with earlier alx1 and ets1 function in promoting skeletogenic mesenchyme ingression. In contrast, repression constructs using the tbr driver, which is active in blastula stage, block ingression. This method thus provides the opportunity to study regulatory requirements of skeletogenesis after ingression, and may be similarly useful in many other developmental contexts.

The roles of Ca2+, downstream protein kinases, and oscillatory signaling in regulating fertilization and the activation of development

Reviews in Developmental Biology have covered the pathways that generate the all-important intracellular calcium (Ca(2+)) signal at fertilization [Miyazaki, S., Shirakawa, H., Nakada, K., Honda, Y., 1993a. Essential role of the inositol 1,4,5-trisphosphate receptor/Ca(2+) release channel in Ca(2+) waves and Ca(2+) oscillations at fertilization of mammalian eggs. Dev. Biol. 158, 62-78; Runft, L., Jaffe, L., Mehlmann, L., 2002. Egg activation at fertilization: where it all begins. Dev. Biol. 245, 237-254] and the different temporal responses of Ca(2+) in many organisms [Stricker, S., 1999. Comparative biology of calcium signaling during fertilization and egg activation in animals. Dev. Biol. 211, 157-176]. Those reviews raise the importance of identifying how Ca(2+) causes the events of egg activation (EEA) and to what extent these temporal Ca(2+) responses encode developmental information. This review covers recent studies that have analyzed how these Ca(2+) signals are interpreted by specific proteins, and how these proteins regulate various EEA responsible for the onset of development. Many of these proteins are protein kinases (CaMKII, PKC, MPF, MAPK, MLCK) whose activity is directly or indirectly regulated by Ca(2+), and whose amount increases during late oocyte maturation. We cover biochemical progress in defining the signaling pathways between Ca(2+) and the EEA, as well as discuss how oscillatory or multiple Ca(2+) signals are likely to have specific advantages biochemically and/or developmentally. These emerging concepts are put into historical context, emphasizing that key contributions have come from many organisms. The intricate interdependence of Ca(2+), Ca(2+)-dependent proteins, and the EEA raise many new questions for future investigations that will provide insight into the extent to which fertilization-associated signaling has long-range implications for development. In addition, answers to these questions should be beneficial to establishing parameters of egg quality for human and animal IVF, as well as improving egg activation protocols for somatic cell nuclear transfer to generate stem cells and save endangered species.

Logic of gene regulatory networks

Regulatory networks of transcription factors and signaling molecules lie at the heart of development. Their architecture implements logic functions whose execution propels cells from one regulatory state to the next, thus driving development forward. As an example of a subcircuit that translates transcriptional input into developmental output, we consider a particularly simple case, the regulatory processes underlying pigment cell formation in sea urchin embryos. The regulatory events in this process can be represented as elementary logic functions.

A missing link in the sea urchin embryo gene regulatory network: hesC and the double-negative specification of micromeres

Specification of sea urchin embryo micromeres occurs early in cleavage, with the establishment of a well defined regulatory state. The architecture of the gene regulatory network controlling the specification process indicates that transcription of the initial tier of control genes depends on a double-negative gate. A gene encoding a transcriptional repressor, pmar1, is activated specifically in micromeres, where it represses transcription of a second repressor that is otherwise active globally. Thus, the micromere-specific control genes, which are the target of the second repressor, are expressed exclusively in this lineage. The double-negative specification gate was logically required from the results of numerous prior experiments, but the identity of the gene encoding the second repressor remained elusive. Here we show that hesC is this gene, and we demonstrate experimentally all of its predicted functions, including global repression of micromere-specific regulatory genes. As logically required, blockade of hesC mRNA translation and global overexpression of pmar1 mRNA have the same effect, which is to cause all of the cells of the embryo to express micromere-specific genes.

Real-time monitoring of functional interactions between upstream and core promoter sequences in living cells of sea urchin embryos

There are some functional compatibilities between upstream and core promoter sequences for transcriptional activation in yeast, Drosophila and mammalian cells. Here we examined whether similar compatibilities exist in sea urchin embryos, and if so, whether they are dynamically regulated during early development. Two reporter plasmids, each containing a test promoter conjugated to either CFP or YFP, were concurrently introduced into embryos, and their expression patterns were studied by fluorescence microscopy. The upstream sequence of the Hemicentrotus pulcherrimus (Hp) OtxE promoter drives the expression of its own core promoter and that of Strongylocentrotus purpuratus (Sp) Spec2a in different embryonic regions, especially at the late gastrula stage. Interestingly, when the four putative transcription factor binding sites of this upstream sequence were individually mutated, the resulting sequences directed different spatiotemporal expression from the same set of two core promoters, indicating that combinations of upstream factors may determine core promoter usage in sea urchin embryos. In addition, the insertion or deletion of consensus or nonconsensus TATA sequences changed the expression profile significantly, irrespective of whether the upstream sequence was intact or mutated. Thus, the TATA sequence may serve as a primary determinant for core promoter selection in these cells.

Sea urchin Forkhead gene family: Phylogeny and embryonic expression

Transcription factors of the Forkhead (Fox) family have been identified in many metazoans, and play important roles in diverse biological processes. Here we define the set of fox genes present in the sea urchin genome, and survey their usage during development. This genome includes 22 fox genes, only three of which were previously known. Of the 23 fox gene subclasses identified in vertebrate genomes, the Strongylocentrotus purpuratus genome has orthologues of all but four (E, H, R and S). Phylogenetic analysis suggests that one S. purpuratus fox gene is equally related to foxA and foxB of vertebrates; this gene defines a new class. Two other genes appear to be specific to the sea urchin, with respect to the genomes so far sequenced. Fox genes orthologous with those of vertebrates but lacking in arthropod or nematode genomes may be deuterostome-specific (subclasses I, J1, J2, L1, M and Q1), while the majority are pan-bilaterian. All but one of the S. purpuratus fox genes (SpfoxQ1) are expressed during embryogenesis, most in a very specific temporal and spatial manner. The sea urchin fox genes clearly execute many different regulatory functions, and almost all of them participate in the process of embryonic development.

The Transcriptome of the Sea Urchin Embryo

The sea urchin Strongylocentrotus purpuratus is a model organism for study of the genomic control circuitry underlying embryonic development. We examined the complete repertoire of genes expressed in the S. purpuratus embryo, up to late gastrula stage, by means of high-resolution custom tiling arrays covering the whole genome. We detected complete spliced structures even for genes known to be expressed at low levels in only a few cells. At least 11,000 to 12,000 genes are used in embryogenesis. These include most of the genes encoding transcription factors and signaling proteins, as well as some classes of general cytoskeletal and metabolic proteins, but only a minor fraction of genes encoding immune functions and sensory receptors. Thousands of small asymmetric transcripts of unknown function were also detected in intergenic regions throughout the genome. The tiling array data were used to correct and authenticate several thousand gene models during the genome annotation process.