Analysis of the P. lividus sea urchin genome highlights contrasting trends of genomic and regulatory evolution in deuterostomes

Sea urchins are emblematic models in developmental biology and display several characteristics that set them apart from other deuterostomes. To uncover the genomic cues that may underlie these specificities, we generated a chromosome-scale genome assembly for the sea urchin Paracentrotus lividus and an extensive gene expression and epigenetic profiles of its embryonic development. We found that, unlike vertebrates, sea urchins retained ancestral chromosomal linkages but underwent very fast intrachromosomal gene order mixing. We identified a burst of gene duplication in the echinoid lineage and showed that some of these expanded genes have been recruited in novel structures (water vascular system, Aristotle's lantern, and skeletogenic micromere lineage). Finally, we identified gene-regulatory modules conserved between sea urchins and chordates. Our results suggest that gene-regulatory networks controlling development can be conserved despite extensive gene order rearrangement.

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.

Monoplacophoran mitochondrial genomes: convergent gene arrangements and little phylogenetic signal

BACKGROUND

Although recent studies have greatly advanced understanding of deep molluscan phylogeny, placement of some taxa remains uncertain as different datasets support competing class-relationships. Traditionally, morphologists have placed Monoplacophora, a group of morphologically simple, limpet-like molluscs as sister group to all other conchiferans (shelled molluscs other than Polyplacophora), a grouping that is supported by the latest large-scale phylogenomic study that includes Laevipilina. However, molecular datasets dominated by nuclear ribosomal genes support Monoplacophora + Polyplacophora (Serialia). Here, we evaluate the potential of mitochondrial genome data for resolving placement of Monoplacophora.

RESULTS

Two complete (Laevipilina antarctica and Vema ewingi) and one partial (Laevipilina hyalina) mitochondrial genomes were sequenced, assembled, and compared. All three genomes show a highly similar architecture including an unusually high number of non-coding regions. Comparison of monoplacophoran gene order shows a gene arrangement pattern not previously reported; there is an inversion of one large gene cluster. Our reanalyses of recently published polyplacophoran mitogenomes show, however, that this feature is also present in some chiton species. Maximum Likelihood and Bayesian Inference analyses of 13 mitochondrial protein-coding genes failed to robustly place Monoplacophora and hypothesis testing could not reject any of the evaluated placements of Monoplacophora.

CONCLUSIONS

Under both serialian or aculiferan-conchiferan scenarios, the observed gene cluster inversion appears to be a convergent evolution of gene arrangements in molluscs. Our phylogenetic results are inconclusive and sensitive to taxon sampling. Aculifera (Polyplacophora + Aplacophora) and Conchifera were never recovered. However, some analyses recovered Serialia (Monoplacophora + Polyplacophora), Diasoma (Bivalvia + Scaphopoda) or Pleistomollusca (Bivalvia + Gastropoda). Although we could not shed light on deep evolutionary traits of Mollusca we found unique patterns of gene arrangements that are common to monoplacophoran and chitonine polyplacophoran species but not to acanthochitonine Polyplacophora. Complete mitochondrial genome of Laevipilina antarctica.

RNA-Seq identifies SPGs as a ventral skeletal patterning cue in sea urchins

The sea urchin larval skeleton offers a simple model for formation of developmental patterns. The calcium carbonate skeleton is secreted by primary mesenchyme cells (PMCs) in response to largely unknown patterning cues expressed by the ectoderm. To discover novel ectodermal cues, we performed an unbiased RNA-Seq-based screen and functionally tested candidates; we thereby identified several novel skeletal patterning cues. Among these, we show that SLC26a2/7 is a ventrally expressed sulfate transporter that promotes a ventral accumulation of sulfated proteoglycans, which is required for ventral PMC positioning and skeletal patterning. We show that the effects of SLC perturbation are mimicked by manipulation of either external sulfate levels or proteoglycan sulfation. These results identify novel skeletal patterning genes and demonstrate that ventral proteoglycan sulfation serves as a positional cue for sea urchin skeletal patterning.

The nervous system of Xenacoelomorpha: a genomic perspective

Xenacoelomorpha is, most probably, a monophyletic group that includes three clades: Acoela, Nemertodermatida and Xenoturbellida. The group still has contentious phylogenetic affinities; though most authors place it as the sister group of the remaining bilaterians, some would include it as a fourth phylum within the Deuterostomia. Over the past few years, our group, along with others, has undertaken a systematic study of the microscopic anatomy of these worms; our main aim is to understand the structure and development of the nervous system. This research plan has been aided by the use of molecular/developmental tools, the most important of which has been the sequencing of the complete genomes and transcriptomes of different members of the three clades. The data obtained has been used to analyse the evolutionary history of gene families and to study their expression patterns during development, in both space and time. A major focus of our research is the origin of 'cephalized' (centralized) nervous systems. How complex brains are assembled from simpler neuronal arrays has been a matter of intense debate for at least 100 years. We are now tackling this issue using Xenacoelomorpha models. These represent an ideal system for this work because the members of the three clades have nervous systems with different degrees of cephalization; from the relatively simple sub-epithelial net of Xenoturbella to the compact brain of acoels. How this process of 'progressive' cephalization is reflected in the genomes or transcriptomes of these three groups of animals is the subject of this paper.

The complete mitochondrial genome of Telmatobufo australis (Amphibia: Anura: Calyptocephalellidae)

The mitochondrial (mt) genome of Telmatobufo australis is a circular molecule of 17,989 bp in length, comprising 13 protein-coding genes, 22 transfer RNA genes, and two ribosomal RNA genes. Gene order and content are identical to those previously reported from other neobatrachian mt genomes. Two protein-coding genes (COI and ATP6) presumably used GTG as start codons while COIII possessed an incomplete stop codon.

LTR retroelements are intrinsic components of transcriptional networks in frogs

Background

LTR retroelements (LTR REs) constitute a major group of transposable elements widely distributed in eukaryotic genomes. Through their own mechanism of retrotranscription LTR REs enrich the genomic landscape by providing genetic variability, thus contributing to genome structure and organization. Nonetheless, transcriptomic activity of LTR REs still remains an obscure domain within cell, developmental, and organism biology.

Results

Here we present a first comparative analysis of LTR REs for anuran amphibians based on a full depth coverage transcriptome of the European pool frog, Pelophylax lessonae, the genome of the African clawed frog, Silurana tropicalis (release v7.1), and additional transcriptomes of S. tropicalis and Cyclorana alboguttata. We identified over 1000 copies of LTR REs from all four families (Bel/Pao, Ty1/Copia, Ty3/Gypsy, Retroviridae) in the genome of S. tropicalis and discovered transcripts of several of these elements in all RNA-seq datasets analyzed. Elements of the Ty3/Gypsy family were most active, especially Amn-san elements, which accounted for approximately 0.27% of the genome in Silurana. Some elements exhibited tissue specific expression patterns, for example Hydra1.1 and MuERV-like elements in Pelophylax. In S. tropicalis considerable transcription of LTR REs was observed during embryogenesis as soon as the embryonic genome became activated, i.e. at midblastula transition. In the course of embryonic development the spectrum of transcribed LTR REs changed; during gastrulation and neurulation MuERV-like and SnRV like retroviruses were abundantly transcribed while during organogenesis transcripts of the XEN1 retroviruses became much more active.

Conclusions

The differential expression of LTR REs during embryogenesis in concert with their tissue-specificity and the protein domains they encode are evidence for the functional roles these elements play as integrative parts of complex regulatory networks. Our results support the meanwhile widely accepted concept that retroelements are not simple “junk DNA” or “harmful genomic parasites” but essential components of the transcriptomic machinery in vertebrates.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-626) contains supplementary material, which is available to authorized users.

Massively parallel RNA sequencing identifies a complex immune gene repertoire in the lophotrochozoan Mytilus edulis

The marine mussel Mytilus edulis and its closely related sister species are distributed world-wide and play an important role in coastal ecology and economy. The diversification in different species and their hybrids, broad ecological distribution, as well as the filter feeding mode of life has made this genus an attractive model to investigate physiological and molecular adaptations and responses to various biotic and abiotic environmental factors. In the present study we investigated the immune system of Mytilus, which may contribute to the ecological plasticity of this species. We generated a large Mytilus transcriptome database from different tissues of immune challenged and stress treated individuals from the Baltic Sea using 454 pyrosequencing. Phylogenetic comparison of orthologous groups of 23 species demonstrated the basal position of lophotrochozoans within protostomes. The investigation of immune related transcripts revealed a complex repertoire of innate recognition receptors and downstream pathway members including transcripts for 27 toll-like receptors and 524 C1q domain containing transcripts. NOD-like receptors on the other hand were absent. We also found evidence for sophisticated TNF, autophagy and apoptosis systems as well as for cytokines. Gill tissue and hemocytes showed highest expression of putative immune related contigs and are promising tissues for further functional studies. Our results partly contrast with findings of a less complex immune repertoire in ecdysozoan and other lophotrochozoan protostomes. We show that bivalves are interesting candidates to investigate the evolution of the immune system from basal metazoans to deuterostomes and protostomes and provide a basis for future molecular work directed to immune system functioning in Mytilus.

Data integration using scanners with SQL output - The Bioscanners project at sourceforge

There is currently no standardized approach for parsing output that the numerous bioinformatical tools generate. Because the framework approach of the Bio-toolkits has some shortcomings, we searched for alternative approaches. To this end, we evaluated scanner generators for various programming languages with respect to their potential of creating standalone, small, and fast applications that can easily be delivered on any modern and many ancient operating systems. We developed sample applications that generate standard SQL database code and thereby greatly simplify the parsing work of data integration and data analysis. At the sourceforge project page the source code and some binaries for a selection of our applications are freely available at http://bioscanners.sourceforge.net.

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.

Data integration using scanners with SQL output--the bioscanners project at sourceforge

There is currently no standardized approach for parsing output that the numerous bioinformatical tools generate. Because the framework approach of the Bio-toolkits has some shortcomings, we searched for alternative approaches. To this end, we evaluated scanner generators for various programming languages with respect to their potential of standalone, small and fast applications that can easily delivered on any modern and many ancient operating system. We developed sample applications that generate standard SQL database code and thereby greatly simplify the parsing work of data integration and data analysis. At the sourceforge project page the source code and some binaries for a selection of our applications are freely available at http://bioscanners.sourceforge.net.

Phosphoproteomes of Strongylocentrotus purpuratus shell and tooth matrix: identification of a major acidic sea urchin tooth phosphoprotein, phosphodontin

BACKGROUND

Sea urchin is a major model organism for developmental biology and biomineralization research. However, identification of proteins involved in larval skeleton formation and mineralization processes in the embryo and adult, and the molecular characterization of such proteins, has just gained momentum with the sequencing of the Strongylocentrotus purpuratus genome and the introduction of high-throughput proteomics into the field.

RESULTS

The present report contains the determination of test (shell) and tooth organic matrix phosphoproteomes. Altogether 34 phosphoproteins were identified in the biomineral organic matrices. Most phosphoproteins were specific for one compartment, only two were identified in both matrices. The sea urchin phosphoproteomes contained several obvious orthologs of mammalian proteins, such as a Src family tyrosine kinase, protein kinase C-delta 1, Dickkopf-1 and other signal transduction components, or nucleobindin. In most cases phosphorylation sites were conserved between sea urchin and mammalian proteins. However, the majority of phosphoproteins had no mammalian counterpart. The most interesting of the sea urchin-specific phosphoproteins, from the perspective of biomineralization research, was an abundant highly phosphorylated and very acidic tooth matrix protein composed of 35 very similar short sequence repeats, a predicted N-terminal secretion signal sequence, and an Asp-rich C-terminal motif, contained in [Glean3:18919].

CONCLUSIONS

The 64 phosphorylation sites determined represent the most comprehensive list of experimentally identified sea urchin protein phosphorylation sites at present and are an important addition to the recently analyzed Strongylocentrotus purpuratus shell and tooth proteomes. The identified phosphoproteins included a major, highly phosphorylated protein, [Glean3:18919], for which we suggest the name phosphodontin. Although not sequence-related to such highly phosphorylated acidic mammalian dental phosphoproteins as phosphoryn or dentin matrix protein-1, phosphodontin may perform similar functions in the sea urchin tooth. More than half of the detected proteins were not previously identified at the protein level, thus confirming the existence of proteins only known as genomic sequences previously.

Deeply conserved chordate noncoding sequences preserve genome synteny but do not drive gene duplicate retention

Animal genomes possess highly conserved cis-regulatory sequences that are often found near genes that regulate transcription and development. Researchers have proposed that the strong conservation of these sequences may affect the evolution of the surrounding genome, both by repressing rearrangement, and possibly by promoting duplicate gene retention. Conflicting data, however, have made the validity of these propositions unclear. Here, we use a new computational method to identify phylogenetically conserved noncoding elements (PCNEs) in a manner that is not biased by rearrangement and duplication. This method is powerful enough to identify more than a thousand PCNEs that have been conserved between vertebrates and the basal chordate amphioxus. We test 42 of our PCNEs in transgenic zebrafish assays--including examples from vertebrates and amphioxus--and find that the majority are functional enhancers. We find that PCNEs are enriched around genes with ancient synteny conservation, and that this association is strongest for extragenic PCNEs, suggesting that cis-regulatory interdigitation plays a key role in repressing genome rearrangement. Next, we classify mouse and zebrafish genes according to association with PCNEs, synteny conservation, duplication history, and presence in bidirectional promoter pairs, and use these data to cluster gene functions into a series of distinct evolutionary patterns. These results demonstrate that subfunctionalization of conserved cis-regulation has not been the primary determinate of gene duplicate retention in vertebrates. Instead, the data support the gene balance hypothesis, which proposes that duplicate retention has been driven by selection against dosage imbalances in genes with many protein connections.

Monte Carlo analysis of an ODE Model of the Sea Urchin Endomesoderm Network

BACKGROUND

Gene Regulatory Networks (GRNs) control the differentiation, specification and function of cells at the genomic level. The levels of interactions within large GRNs are of enormous depth and complexity. Details about many GRNs are emerging, but in most cases it is unknown to what extent they control a given process, i.e. the grade of completeness is uncertain. This uncertainty stems from limited experimental data, which is the main bottleneck for creating detailed dynamical models of cellular processes. Parameter estimation for each node is often infeasible for very large GRNs. We propose a method, based on random parameter estimations through Monte-Carlo simulations to measure completeness grades of GRNs.

RESULTS

We developed a heuristic to assess the completeness of large GRNs, using ODE simulations under different conditions and randomly sampled parameter sets to detect parameter-invariant effects of perturbations. To test this heuristic, we constructed the first ODE model of the whole sea urchin endomesoderm GRN, one of the best studied large GRNs. We find that nearly 48% of the parameter-invariant effects correspond with experimental data, which is 65% of the expected optimal agreement obtained from a submodel for which kinetic parameters were estimated and used for simulations. Randomized versions of the model reproduce only 23.5% of the experimental data.

CONCLUSION

The method described in this paper enables an evaluation of network topologies of GRNs without requiring any parameter values. The benefit of this method is exemplified in the first mathematical analysis of the complete Endomesoderm Network Model. The predictions we provide deliver candidate nodes in the network that are likely to be erroneous or miss unknown connections, which may need additional experiments to improve the network topology. This mathematical model can serve as a scaffold for detailed and more realistic models. We propose that our method can be used to assess a completeness grade of any GRN. This could be especially useful for GRNs involved in human diseases, where often the amount of connectivity is unknown and/or many genes/interactions are missing.

Chordin is required for neural but not axial development in sea urchin embryos

The oral-aboral (OA) axis in the sea urchin is specified by the TGFbeta family members Nodal and BMP2/4. Nodal promotes oral specification, whereas BMP2/4, despite being expressed in the oral territory, is required for aboral specification. This study explores the role of Chordin (Chd) during sea urchin embryogenesis. Chd is a secreted BMP inhibitor that plays an important role in axial and neural specification and patterning in Drosophila and vertebrate embryos. In Lytechinus variegatus embryos, Chd and BMP2/4 are functionally antagonistic. Both are expressed in overlapping domains in the oral territory prior to and during gastrulation. Perturbation shows that, surprisingly, Chd is not involved in OA axis specification. Instead, Chd is required both for normal patterning of the ciliary band at the OA boundary and for development of synaptotagmin B-positive (synB) neurons in a manner that is reciprocal with BMP2/4. Chd expression and synB-positive neural development are both downstream from p38 MAPK and Nodal, but not Goosecoid. These data are summarized in a model for synB neural development.

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.

The sea urchin (Strongylocentrotus purpuratus) test and spine proteomes

BACKGROUND

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

RESULTS

We identified 110 proteins as components of sea urchin test and spine organic matrix. Fourty of these proteins occurred in both compartments while others were unique to their respective compartment. More than 95% of the proteins were detected in sea urchin skeletal matrices for the first time. The most abundant protein in both matrices was the previously characterized spicule matrix protein SM50, but at least eight other members of this group, many of them only known as conceptual translation products previously, were identified by mass spectrometric sequence analysis of peptides derived from in vitro matrix degradation. The matrices also contained proteins implicated in biomineralization processes previously by inhibition studies using antibodies or specific enzyme inhibitors, such as matrix metalloproteases and members of the mesenchyme-specific MSP130 family. Other components were carbonic anhydrase, collagens, echinonectin, a alpha2-macroglobulin-like protein and several proteins containing scavenger receptor cysteine-rich domains. A few possible signal transduction pathway components, such as GTP-binding proteins, a semaphorin and a possible tyrosine kinase were also identified.

CONCLUSION

This report presents the most comprehensive list of sea urchin skeletal matrix proteins available at present. The complex mixture of proteins identified in matrices of the sea urchin skeleton may reflect many different aspects of the mineralization process. Because LC-MS/MS-based methods directly measures peptides our results validate many predicted genes and confirm the existence of the corresponding proteins. Considering the many newly identified matrix proteins, this proteomic study may serve as a road map for the further exploration of biomineralization processes in an important model organism.

Early vertebrate whole genome duplications were predated by a period of intense genome rearrangement

Researchers, supported by data from polyploid plants, have suggested that whole genome duplication (WGD) may induce genomic instability and rearrangement, an idea which could have important implications for vertebrate evolution. Benefiting from the newly released amphioxus genome sequence (Branchiostoma floridae), an invertebrate that researchers have hoped is representative of the ancestral chordate genome, we have used gene proximity conservation to estimate rates of genome rearrangement throughout vertebrates and some of their invertebrate ancestors. We find that, while amphioxus remains the best single source of invertebrate information about the early chordate genome, its genome structure is not particularly well conserved and it cannot be considered a fossilization of the vertebrate preduplication genome. In agreement with previous reports, we identify two WGD events in early vertebrates and another in teleost fish. However, we find that the early vertebrate WGD events were not followed by increased rates of genome rearrangement. Indeed, we measure massive genome rearrangement prior to these WGD events. We propose that the vertebrate WGD events may have been symptoms of a preexisting predisposition toward genomic structural change.

GOblet: Annotation of anonymous sequence data with Gene Ontology and Pathway terms

The functional annotation of genomic data has become a major task for the ever-growing number of sequencing projects. In order to address this challenge, we recently developed GOblet, a free web service for the annotation of anonymous sequences with Gene Ontology (GO) terms. However, to overcome limitations of the GO terminology, and to aid in understanding not only single components but as well systemic interactions between the individual components, we have now extended the GOblet web service to integrate also pathway annotations. Furthermore, we extended and upgraded the data analysis pipeline with improved summaries, and added term enrichment and clustering algorithms. Finally, we are now making GOblet available as a stand-alone application for high-throughput processing on local machines. The advantages of this frequently requested feature is that a) the user can avoid restrictions of our web service for uploading and processing large amounts of data, and that b) confidential data can be analysed without insecure transfer to a public web server. The stand-alone version of the web service has been implemented using platform independent Tcl-scripts, which can be run with just a single runtime file utilizing the Starkit technology. The GOblet web service and the stand-alone application are freely available at http://goblet.molgen.mpg.de.

Evolution of a core gene network for skeletogenesis in chordates

The skeleton is one of the most important features for the reconstruction of vertebrate phylogeny but few data are available to understand its molecular origin. In mammals the Runt genes are central regulators of skeletogenesis. Runx2 was shown to be essential for osteoblast differentiation, tooth development, and bone formation. Both Runx2 and Runx3 are essential for chondrocyte maturation. Furthermore, Runx2 directly regulates Indian hedgehog expression, a master coordinator of skeletal development. To clarify the correlation of Runt gene evolution and the emergence of cartilage and bone in vertebrates, we cloned the Runt genes from hagfish as representative of jawless fish (MgRunxA, MgRunxB) and from dogfish as representative of jawed cartilaginous fish (ScRunx1-3). According to our phylogenetic reconstruction the stem species of chordates harboured a single Runt gene and thereafter Runt locus duplications occurred during early vertebrate evolution. All newly isolated Runt genes were expressed in cartilage according to quantitative PCR. In situ hybridisation confirmed high MgRunxA expression in hard cartilage of hagfish. In dogfish ScRunx2 and ScRunx3 were expressed in embryonal cartilage whereas all three Runt genes were detected in teeth and placoid scales. In cephalochordates (lancelets) Runt, Hedgehog and SoxE were strongly expressed in the gill bars and expression of Runt and Hedgehog was found in endo- as well as ectodermal cells. Furthermore we demonstrate that the lancelet Runt protein binds to Runt binding sites in the lancelet Hedgehog promoter and regulates its activity. Together, these results suggest that Runt and Hedgehog were part of a core gene network for cartilage formation, which was already active in the gill bars of the common ancestor of cephalochordates and vertebrates and diversified after Runt duplications had occurred during vertebrate evolution. The similarities in expression patterns of Runt genes support the view that teeth and placoid scales evolved from a homologous developmental module.

A global view of gene expression in lithium and zinc treated sea urchin embryos: new components of gene regulatory networks

BACKGROUND

The genome of the sea urchin Strongylocentrotus purpuratus has recently been sequenced because it is a major model system for the study of gene regulatory networks. Embryonic expression patterns for most genes are unknown, however.

RESULTS

Using large-scale screens on arrays carrying 50% to 70% of all genes, we identified novel territory-specific markers. Our strategy was based on computational selection of genes that are differentially expressed in lithium-treated embryos, which form excess endomesoderm, and in zinc-treated embryos, in which endomesoderm specification is blocked. Whole-mount in situ hybridization (WISH) analysis of 700 genes indicates that the apical organ region is eliminated in lithium-treated embryos. Conversely, apical and specifically neural markers are expressed more broadly in zinc-treated embryos, whereas endomesoderm signaling is severely reduced. Strikingly, the number of serotonergic neurons is amplified by at least tenfold in zinc-treated embryos. WISH analysis further indicates that there is crosstalk between the Wnt (wingless int), Notch, and fibroblast growth factor signaling pathways in secondary mesoderm cell specification and differentiation, similar to signaling cascades that function during development of presomitic mesoderm in mouse embryogenesis. We provide differential expression data for more than 4,000 genes and WISH patterns of more than 250 genes, and more than 2,400 annotated WISH images.

CONCLUSION

Our work provides tissue-specific expression patterns for a large fraction of the sea urchin genes that have not yet been included in existing regulatory networks and await functional integration. Furthermore, we noted neuron-inducing activity of zinc on embryonic development; this is the first observation of such activity in any organism.

Modeling development: spikes of the sea urchin

Modeling of specification events during development poses new challenges to biochemical modeling. These include data limitations and a notorious absence of homeostasis in developing systems. The sea urchin is one of the best studied model organisms concerning development and a network, the Endomesoderm Network, has been proposed that is presumed to control endoderm and mesoderm specification in the embryo of Strongy-locentrotus purpuratus. We have constructed a dynamic model of a subnetwork of the Endomesoderm Network. In constructing the model, we had to resolve the following issues: choice of appropriate subsystem, assignment of embryonic data to cellular model, choice of appropriate kinetics. Although the resulting model is capable of reproducing fractions of the experimental data, it falls short of reproducing specification of cell types. These findings can facilitate the refinement of the Endomesoderm Network.

The genome of the sea urchin Strongylocentrotus purpuratus

We report the sequence and analysis of the 814-megabase genome of the sea urchin Strongylocentrotus purpuratus, a model for developmental and systems biology. The sequencing strategy combined whole-genome shotgun and bacterial artificial chromosome (BAC) sequences. This use of BAC clones, aided by a pooling strategy, overcame difficulties associated with high heterozygosity of the genome. The genome encodes about 23,300 genes, including many previously thought to be vertebrate innovations or known only outside the deuterostomes. This echinoderm genome provides an evolutionary outgroup for the chordates and yields insights into the evolution of deuterostomes.

Gene identification and analysis of transcripts differentially regulated in fracture healing by EST sequencing in the domestic sheep

BACKGROUND

The sheep is an important model animal for testing novel fracture treatments and other medical applications. Despite these medical uses and the well known economic and cultural importance of the sheep, relatively little research has been performed into sheep genetics, and DNA sequences are available for only a small number of sheep genes.

RESULTS

In this work we have sequenced over 47 thousand expressed sequence tags (ESTs) from libraries developed from healing bone in a sheep model of fracture healing. These ESTs were clustered with the previously available 10 thousand sheep ESTs to a total of 19087 contigs with an average length of 603 nucleotides. We used the newly identified sequences to develop RT-PCR assays for 78 sheep genes and measured differential expression during the course of fracture healing between days 7 and 42 postfracture. All genes showed significant shifts at one or more time points. 23 of the genes were differentially expressed between postfracture days 7 and 10, which could reflect an important role for these genes for the initiation of osteogenesis.

CONCLUSION

The sequences we have identified in this work are a valuable resource for future studies on musculoskeletal healing and regeneration using sheep and represent an important head-start for genomic sequencing projects for Ovis aries, with partial or complete sequences being made available for over 5,800 previously unsequenced sheep genes.

Timing and mechanism of ancient vertebrate genome duplications – the adventure of a hypothesis

Complete genome doubling has long-term consequences for the genome structure and the subsequent evolution of an organism. It has been suggested that two genome duplications occurred at the origin of vertebrates (known as the 2R hypothesis). However, there has been considerable debate as to whether these were two successive duplications, or whether a single duplication occurred, followed by large-scale segmental duplications. In this article, we review and compare the evidence for the 2R duplications from vertebrate genomes with similar data from other more recent polyploids.

Nodal/activin signaling establishes oral-aboral polarity in the early sea urchin embryo

Components of the Wnt signaling pathway are involved in patterning the sea urchin primary or animal-vegetal (AV) axis, but the molecular cues that pattern the secondary embryonic axis, the aboral/oral (AO) axis, are not known. In an analysis of signaling molecules that influence patterning along the sea urchin embryonic axes, we found that members of the activin subfamily of transforming growth factor-beta (TGF-beta) signaling molecules influence the establishment of AO polarities in the early embryo. Injection of activin mRNAs into fertilized eggs or treatment with exogenously applied recombinant activin altered the allocation of ectodermal fates and ventralized the embryo. The phenotypes observed resemble the ventralized phenotype previously reported for NiCl2, a known disrupter of AO patterning. Sensitivity to exogenous activin occurs between fertilization and the late blastula stage, which is also the time of highest NiCl2 sensitivity. These results argue that specification of fates along the embryonic AO axis involves TGF-beta signaling. To further examine TGF-beta signaling in these embryos, we cloned an endogenous TGF-beta from sea urchin embryos that is a member of the activin subfamily, SpNodal, and show through gain of function analysis that it recapitulates results obtained with exogenous activins and NiCl2. The expression pattern of SpNodal is consistent with a role for nodal signaling in the establishment of fates along the AO axis. Loss of function experiments using SpNodal antisense morpholinos also support a role for SpNodal in the establishment of the AO axis.

On the origin of the chordate central nervous system: expression of onecut in the sea urchin embryo

We identified a transcription factor of the onecut class in the sea urchin Strongylocentrotus purpuratus that represents an ortholog of the mammalian gene HNF6, the founding member of the onecut class of proteins. The isolated sea urchin gene, named SpOnecut, encodes a protein of 483 amino acids with one cut domain and a homeodomain. Phylogenetic analysis clearly places the sea urchin gene into this family, most closely related to the ascidian onecut gene HNF-6. Nevertheless, phylogenetic analysis reveals a difficult phylogeny indicating that certain members of the family evolve more rapidly than others and also that the cut domain and homeodomain evolve at a different pace. In fly, worm, ascidian, and teleost fish, the onecut genes isolated so far are exclusively expressed in cells of the central nervous system (CNS), whereas in mammals the two copies of the gene have acquired additional functions in liver and pancreas development. In the sea urchin embryo, expression is first detected in the emerging ciliary band at the late blastula stage. During the gastrula stage, expression is limited to the ciliary band. In the early pluteus stage, SpOnecut is expressed at the apical organ and the elongating arms but continues most prominently in the ciliary band. This is the first gene known that exclusively marks the ciliary band and therein the apical organ in a pluteus larva, whereas chordate orthologs execute essential functions in dorsal CNS development. The significance of this finding for the hypothesis that the ciliary bands and apical organs of the hypothetical "dipleurula"-like chordate ancestor and the chordate/vertebrate CNS are of common origin is discussed.

Expression of an NK2 homeodomain gene in the apical ectoderm defines a new territory in the early sea urchin embryo

We have identified an NK2 family homeodomain transcription factor, SpNK2.1, in the sea urchin Strongylocentrotus purpuratus whose transcripts are initially detected within the apical plate ectoderm of the hatching blastula and are confined to the apical organ at least through 2 weeks of development. Protein localization studies demonstrate that SpNK2.1 is restricted to the apical plate epithelium, but is excluded from the nucleus of serotonergic neurons. The expression profile of SpNK2.1 is dictated via two separate regulatory systems. Initially, SpNK2.1 is restricted to the apical pole domain by beta-catenin-dependent processes operating along the animal-vegetal axis, as evidenced by an expansion of SpNK2.1 expression upon cadherin overexpression. Starting at gastrulation, expression in the apical plate is maintained by SpDri, the sea urchin orthologue of dead ringer. Abrogation of SpDri results in the downregulation of SpNK2.1 after gastrulation, but SpDri is not necessary for the initial activation of SpNK2.1. Loss of function experiments using SpNK2.1-specific morpholino antisense oligonucleotides and SpNK2.1 overexpression experiments do not disrupt embryonic development and have no effect upon the development of neuronal components of the apical organ. Nonetheless, SpNK2.1 defines a new early territory of the sea urchin embryo.

Cloning of a novel phospholipase C-delta isoform from pacific purple sea urchin (Strongylocentrotus purpuratus) gametes and its expression during early embryonic development

Calcium (Ca(2+)) is a ubiquitous intracellular messenger, controlling a diverse range of cellular processes, including fertilization and development of the embryo. One of the key mechanisms involved in triggering intracellular calcium release is the generation of the second messenger inositol-1,4,5-phosphate (IP(3)) by the phospholipase C (PLC) class of enzymes. Although five distinct forms of PLC have been identified in mammals (beta, gamma, delta, epsilon, and zeta), only one, PLCgamma, has thus far been detected in echinoderms. In the present study, we describe the isolation of a cDNA encoding a novel PLC isoform of the delta (delta) subclass, PLC-deltasu, from the egg of the Pacific purple sea urchin Strongylocentrotus purpuratus. We also demonstrate the presence of this PLC within the sperm and in the early embryo. The PLC-deltasu cDNA (2.44kb) encodes a 742 amino acid polypeptide with an open reading frame of 84.6kDa and a pI of 6.04. All of the characteristic domains found in mammalian PLCdelta isoforms (PH domain, EF hands, an X-Y catalytic region, and a C2 domain) are present in PLC-deltasu. A homology search revealed that PLC-deltasu shares most sequence identity with bovine PLCdelta2 (39%). We present evidence that PLC-deltasu is expressed in unfertilized eggs, fertilized eggs, and in the early embryo. In addition to Northern and polymerase chain reaction (PCR) analyses, in situ hybridization experiments further demonstrated that the embryonic regions within which the PLC-deltasu transcript can be detected during early embryonic development are associated with the highest levels of proliferative activity, suggesting a possible involvement with metabolism or cell cycle regulation.

Generation, annotation, evolutionary analysis, and database integration of 20,000 unique sea urchin EST clusters

Together with the hemichordates, sea urchins represent basal groups of nonchordate invertebrate deuterostomes that occupy a key position in bilaterian evolution. Because sea urchin embryos are also amenable to functional studies, the sea urchin system has emerged as one of the leading models for the analysis of the function of genomic regulatory networks that control development. We have analyzed a total of 107,283 cDNA clones of libraries that span the development of the sea urchin Strongylocentrotus purpuratus. Normalization by oligonucleotide fingerprinting, EST sequencing and sequence clustering resulted in an EST catalog comprised of 20,000 unique genes or gene fragments. Around 7000 of the unique EST consensus sequences were associated with molecular and developmental functions. Phylogenetic comparison of the identified genes to the genome of the urochordate Ciona intestinalis indicate that at least one quarter of the genes thought to be chordate specific were already present at the base of deuterostome evolution. Comparison of the number of gene copies in sea urchins to those in chordates and vertebrates indicates that the sea urchin genome has not undergone extensive gene or complete genome duplications. The established unique gene set represents an essential tool for the annotation and assembly of the forthcoming sea urchin genome sequence. All cDNA clones and filters of all analyzed libraries are available from the resource center of the German genome project at http://www.rzpd.de.

A single amphioxus and sea urchin runt-gene suggests that runt-gene duplications occurred in early chordate evolution

Runt-homologous molecules are characterized by their DNA binding runt-domain which is highly conserved within bilaterians. The three mammalian runt-genes are master regulators in cartilage/bone formation and hematopoiesis. Historically these features evolved in Craniota and might have been promoted by runt-gene duplication events. The purpose of this study was therefore to investigate how many runt-genes exist in the stem species of chordates, by analyzing the number of runt-genes in what is likely to be the closest living relative of Craniota-amphioxus. To acquire further insight into the possible role of runt-genes in early chordate evolution we have determined the number of runt-genes in sea urchins and have analyzed the runt-expression pattern in this species. Our findings demonstrate the presence of a single runt-gene in amphioxus and sea urchin, which makes it highly likely that the stem species of chordates harbored only a single runt-gene. This suggests that runt-gene duplications occurred later in chordate phylogeny, and are possibly also associated with the evolution of features such as hematopoiesis, cartilage and bone development. In sea urchin embryos runt-expression involves cells of endodermal, mesodermal and ectodermal origin. This complex pattern of expression might reflect the multiple roles played by runt-genes in mammals. A strong runt-signal in the gastrointestinal tract of the sea urchin is in line with runt-expression in the intestine of nematodes and in the murine gastrointestinal tract, and seems to be one of the phylogenetically ancient runt-expression domains.

New evidence for genome-wide duplications at the origin of vertebrates using an amphioxus gene set and completed animal genomes

The 2R hypothesis predicting two genome duplications at the origin of vertebrates is highly controversial. Studies published so far include limited sequence data from organisms close to the hypothesized genome duplications. Through the comparison of a gene catalog from amphioxus, the closest living invertebrate relative of vertebrates, to 3453 single-copy genes orthologous between Caenorhabditis elegans (C), Drosophila melanogaster (D), and Saccharomyces cerevisiae (Y), and to Ciona intestinalis ESTs, mouse, and human genes, we show with a large number of genes that the gene duplication activity is significantly higher after the separation of amphioxus and the vertebrate lineages, which we estimate at 650 million years (Myr). The majority of human orthologs of 195 CDY groups that could be dated by the molecular clock appear to be duplicated between 300 and 680 Myr with a mean at 488 million years ago (Mya). We detected 485 duplicated chromosomal segments in the human genome containing CDY orthologs, 331 of which are found duplicated in the mouse genome and within regions syntenic between human and mouse, indicating that these were generated earlier than the human-mouse split. Model based calculations of the codon substitution rate of the human genes included in these segments agree with the molecular clock duplication time-scale prediction. Our results favor at least one large duplication event at the origin of vertebrates, followed by smaller scale duplication closer to the bird-mammalian split.

brachyury Target genes in the early sea urchin embryo isolated by differential macroarray screening

Brachyury is a transcription factor that functions in gastrulation and endoderm development throughout the Bilateria. Here, we identify genes that are expressed downstream of brachyury during gastrulation of the sea urchin embryo. Screens with two different complex probes generated by subtractive hybridization were carried out on high-density arrays of embryonic cDNA libraries. An mRNA sequence population from embryos expressing brachyury at its peak stage of expression was subtracted with message sequence from embryos in which Brachyury function had been "knocked-out" by injection of a morpholine-substituted antisense oligonucleotide to generate a differential probe for brachyury target genes. Another probe was made by using an mRNA population from embryos that mis-express brachyury at a stage just prior to the normal onset of expression, subtracted with message sequence taken from normal embryos at this stage. Screens carried out with these probes target overlapping but distinct sets of downstream genes. After partial sequence characterization, promising genes were independently analyzed by quantitative real-time PCR and by in situ hybridization. Two major classes of genes emerge in this study: genes expressed in the subset of the secondary mesenchyme cells (SMC) that will become pigment cells, and genes that are expressed in portions of the endoderm coincident with brachyury expression. The latter genes are candidates for direct transcriptional targets of Brachyury. Some of the endodermal genes that respond to Brachyury are cytoskeletal modulators that may play a role in gut morphogenesis. This finding is consistent with the block in gastrulation induced by interfering with Brachyury function in sea urchins, and with known or suggested Brachyury function in other species. Other endodermal target genes are expressed in the archenteron and might be terminal differentiation enzymes of the gut. Brachyury expression patterns for Strongylocentrotus purpuratus reported in this paper are entirely consistent with data from other echinoderm species. Brachyury expression in the vegetal plate is confined to the presumptive endodermal cells. Therefore, the SMC genes are likely to be indirect targets of Brachyury-induced signaling from the surrounding endoderm to the central mesoderm, or the effects on these genes may be indirect consequences of gross disruption of the vegetal plate. These results and other data suggest that the brachyury gene transduces information about the state of endodermal specification to genes that modulate morphogenesis and genes that perform terminal functions in the gut.

A sea urchin genome project: sequence scan, virtual map, and additional resources

Results of a first-stage Sea Urchin Genome Project are summarized here. The species chosen was Strongylocentrotus purpuratus, a research model of major importance in developmental and molecular biology. A virtual map of the genome was constructed by sequencing the ends of 76,020 bacterial artificial chromosome (BAC) recombinants (average length, 125 kb). The BAC-end sequence tag connectors (STCs) occur an average of 10 kb apart, and, together with restriction digest patterns recorded for the same BAC clones, they provide immediate access to contigs of several hundred kilobases surrounding any gene of interest. The STCs survey >5% of the genome and provide the estimate that this genome contains approximately 27,350 protein-coding genes. The frequency distribution and canonical sequences of all middle and highly repetitive sequence families in the genome were obtained from the STCs as well. The 500-kb Hox gene complex of this species is being sequenced in its entirety. In addition, arrayed cDNA libraries of >10(5) clones each were constructed from every major stage of embryogenesis, several individual cell types, and adult tissues and are available to the community. The accumulated STC data and an expanding expressed sequence tag database (at present including >12, 000 sequences) have been reported to GenBank and are accessible on public web sites.

Large-scale clustering of cDNA-fingerprinting data

Clustering is one of the main mathematical challenges in large-scale gene expression analysis. We describe a clustering procedure based on a sequential k-means algorithm with additional refinements that is able to handle high-throughput data in the order of hundreds of thousands of data items measured on hundreds of variables. The practical motivation for our algorithm is oligonucleotide fingerprinting-a method for simultaneous determination of expression level for every active gene of a specific tissue-although the algorithm can be applied as well to other large-scale projects like EST clustering and qualitative clustering of DNA-chip data. As a pairwise similarity measure between two p-dimensional data points, x and y, we introduce mutual information that can be interpreted as the amount of information about x in y, and vice versa. We show that for our purposes this measure is superior to commonly used metric distances, for example, Euclidean distance. We also introduce a modified version of mutual information as a novel method for validating clustering results when the true clustering is known. The performance of our algorithm with respect to experimental noise is shown by extensive simulation studies. The algorithm is tested on a subset of 2029 cDNA clones coming from 15 different genes from a cDNA library derived from human dendritic cells. Furthermore, the clustering of these 2029 cDNA clones is demonstrated when the entire set of 76,032 cDNA clones is processed.

Toward the gene catalogue of sea urchin development: the construction and analysis of an unfertilized egg cDNA library highly normalized by oligonucleotide fingerprinting

We describe the use of oligonucleotide fingerprinting for the generation of a normalized cDNA library from unfertilized sea urchin eggs and report the preliminary analysis of this library, which resulted in the establishment of a partial gene catalogue of the sea urchin egg. In an analysis of 21,925 cDNA clones by hybridization with 217 oligonucleotide probes, we were able to identify 6291 clusters corresponding to different transcripts, ranging in size from 1 to 265 clones. This corresponds to an average 3.5-fold normalization of the starting library. The normalized library represents about one-third of all genes expressed in the sea urchin egg. To generate sequence information for the transcripts represented by the clusters, representative clones selected from 711 clusters were sequenced. The construction and preliminary analysis of the normalized library are the first steps in the assembly of an increasingly complete collection of maternal genes expressed in the sea urchin egg, which will provide a number of insights into the early development of this well-characterized model organism.