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

A pancreatic exocrine-like cell regulatory circuit operating in the upper stomach of the sea urchin Strongylocentrotus purpuratus larva

Background

Digestive cells are present in all metazoans and provide the energy necessary for the whole organism. Pancreatic exocrine cells are a unique vertebrate cell type involved in extracellular digestion of a wide range of nutrients. Although the organization and regulation of this cell type is intensively studied in vertebrates, its evolutionary history is still unknown. In order to understand which are the elements that define the pancreatic exocrine phenotype, we have analyzed the expression of genes that contribute to specification and function of this cell-type in an early branching deuterostome, the sea urchin Strongylocentrotus purpuratus.

Results

We defined the spatial and temporal expression of sea urchin orthologs of pancreatic exocrine genes and described a unique population of cells clustered in the upper stomach of the sea urchin embryo where exocrine markers are co-expressed. We used a combination of perturbation analysis, drug and feeding experiments and found that in these cells of the sea urchin embryo gene expression and gene regulatory interactions resemble that of bona fide pancreatic exocrine cells. We show that the sea urchin Ptf1a, a key transcriptional activator of digestive enzymes in pancreatic exocrine cells, can substitute for its vertebrate ortholog in activating downstream genes.

Conclusions

Collectively, our study is the first to show with molecular tools that defining features of a vertebrate cell-type, the pancreatic exocrine cell, are shared by a non-vertebrate deuterostome. Our results indicate that the functional cell-type unit of the vertebrate pancreas may evolutionarily predate the emergence of the pancreas as a discrete organ. From an evolutionary perspective, these results encourage to further explore the homologs of other vertebrate cell-types in traditional or newly emerging deuterostome systems.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0686-0) contains supplementary material, which is available to authorized users.

Do echinoderm genomes measure up?

Echinoderm genome sequences are a corpus of useful information about a clade of animals that serve as research models in fields ranging from marine ecology to cell and developmental biology. Genomic information from echinoids has contributed to insights into the gene interactions that drive the developmental process at the molecular level. Such insights often rely heavily on genomic information and the kinds of questions that can be asked thus depend on the quality of the sequence information. Here we describe the history of echinoderm genomic sequence assembly and present details about the quality of the data obtained. All of the sequence information discussed here is posted on the echinoderm information web system, Echinobase.org.

Tools for sea urchin genomic analysis

The Sea Urchin Genome Project Web site, SpBase ( http://SpBase.org ), in association with a suite of publicly available sequence comparison tools provides a platform from which to analyze genes and genomic sequences of sea urchin. This information system is specifically designed to support laboratory bench studies in cell and molecular biology. In particular these tools and datasets have supported the description of the gene regulatory networks of the purple sea urchin S. purpuratus. This chapter details methods to undertake in the first steps to find genes and noncoding regulatory sequences for further analysis.

Myogenesis in the sea urchin embryo: the molecular fingerprint of the myoblast precursors

Background

In sea urchin larvae the circumesophageal fibers form a prominent muscle system of mesodermal origin. Although the morphology and later development of this muscle system has been well-described, little is known about the molecular signature of these cells or their precise origin in the early embryo. As an invertebrate deuterostome that is more closely related to the vertebrates than other commonly used model systems in myogenesis, the sea urchin fills an important phylogenetic gap and provides a unique perspective on the evolution of muscle cell development.

Results

Here, we present a comprehensive description of the development of the sea urchin larval circumesophageal muscle lineage beginning with its mesodermal origin using high-resolution localization of the expression of several myogenic transcriptional regulators and differentiation genes. A few myoblasts are bilaterally distributed at the oral vegetal side of the tip of the archenteron and first appear at the late gastrula stage. The expression of the differentiation genes Myosin Heavy Chain, Tropomyosin I and II, as well as the regulatory genes MyoD2, FoxF, FoxC, FoxL1, Myocardin, Twist, and Tbx6 uniquely identify these cells. Interestingly, evolutionarily conserved myogenic factors such as Mef2, MyoR and Six1/2 are not expressed in sea urchin myoblasts but are found in other mesodermal domains of the tip of the archenteron. The regulatory states of these domains were characterized in detail. Moreover, using a combinatorial analysis of gene expression we followed the development of the FoxF/FoxC positive cells from the onset of expression to the end of gastrulation. Our data allowed us to build a complete map of the Non-Skeletogenic Mesoderm at the very early gastrula stage, in which specific molecular signatures identify the precursors of different cell types. Among them, a small group of cells within the FoxY domain, which also express FoxC and SoxE, have been identified as plausible myoblast precursors. Together, these data support a very early gastrula stage segregation of the myogenic lineage.

Conclusions

From this analysis, we are able to precisely define the regulatory and differentiation signatures of the circumesophageal muscle in the sea urchin embryo. Our findings have important implications in understanding the evolution of development of the muscle cell lineage at the molecular level. The data presented here suggest a high level of conservation of the myogenic specification mechanisms across wide phylogenetic distances, but also reveal clear cases of gene cooption.

Sequencing and analysis of the gastrula transcriptome of the brittle star Ophiocoma wendtii

BACKGROUND

The gastrula stage represents the point in development at which the three primary germ layers diverge. At this point the gene regulatory networks that specify the germ layers are established and the genes that define the differentiated states of the tissues have begun to be activated. These networks have been well-characterized in sea urchins, but not in other echinoderms. Embryos of the brittle star Ophiocoma wendtii share a number of developmental features with sea urchin embryos, including the ingression of mesenchyme cells that give rise to an embryonic skeleton. Notable differences are that no micromeres are formed during cleavage divisions and no pigment cells are formed during development to the pluteus larval stage. More subtle changes in timing of developmental events also occur. To explore the molecular basis for the similarities and differences between these two echinoderms, we have sequenced and characterized the gastrula transcriptome of O. wendtii.

METHODS

Development of Ophiocoma wendtii embryos was characterized and RNA was isolated from the gastrula stage. A transcriptome data base was generated from this RNA and was analyzed using a variety of methods to identify transcripts expressed and to compare those transcripts to those expressed at the gastrula stage in other organisms.

RESULTS

Using existing databases, we identified brittle star transcripts that correspond to 3,385 genes, including 1,863 genes shared with the sea urchin Strongylocentrotus purpuratus gastrula transcriptome. We characterized the functional classes of genes present in the transcriptome and compared them to those found in this sea urchin. We then examined those members of the germ-layer specific gene regulatory networks (GRNs) of S. purpuratus that are expressed in the O. wendtii gastrula. Our results indicate that there is a shared 'genetic toolkit' central to the echinoderm gastrula, a key stage in embryonic development, though there are also differences that reflect changes in developmental processes.

CONCLUSIONS

The brittle star expresses genes representing all functional classes at the gastrula stage. Brittle stars and sea urchins have comparable numbers of each class of genes and share many of the genes expressed at gastrulation. Examination of the brittle star genes in which sea urchin orthologs are utilized in germ layer specification reveals a relatively higher level of conservation of key regulatory components compared to the overall transcriptome. We also identify genes that were either lost or whose temporal expression has diverged from that of sea urchins.

Molecular profiling of marine fauna: integration of omics with environmental assessment of the world's oceans

Many species that contribute to the commercial and ecological richness of our marine ecosystems are harbingers of environmental change. The ability of organisms to rapidly detect and respond to changes in the surrounding environment represents the foundation for application of molecular profiling technologies towards marine sentinel species in an attempt to identify signature profiles that may reside within the transcriptome, proteome, or metabolome and that are indicative of a particular environmental exposure event. The current review highlights recent examples of the biological information obtained for marine sentinel teleosts, mammals, and invertebrates. While in its infancy, such basal information can provide a systems biology framework in the detection and evaluation of environmental chemical contaminant effects on marine fauna. Repeated evaluation across different seasons and local marine environs will lead to discrimination between signature profiles representing normal variation within the complex milieu of environmental factors that trigger biological response in a given sentinel species and permit a greater understanding of normal versus anthropogenic-associated modulation of biological pathways, which prove detrimental to marine fauna. It is anticipated that incorporation of contaminant-specific molecular signatures into current risk assessment paradigms will lead to enhanced wildlife management strategies that minimize the impacts of our industrialized society on marine ecosystems.

In silico characterization of the neural alpha tubulin gene promoter of the sea urchin embryo Paracentrotus lividus by phylogenetic footprinting

During Paracentrotus lividus sea urchin embryo development one alpha and one beta tubulin genes are expressed specifically in the neural cells and they are early end output of the gene regulatory network that specifies the neural commitment. In this paper we have used a comparative genomics approach to identify conserved regulatory elements in the P. lividus neural alpha tubulin gene. To this purpose, we have first isolated a genomic clone containing the entire gene plus 4.5 Kb of 5' upstream sequences. Then, we have shown by gene transfer experiments that its non-coding region drives the spatio-temporal gene expression corresponding substantially to that of the endogenous gene. In addition, we have identified by genome and EST sequence analysis the S. purpuratus alpha tubulin orthologous gene and we propose a revised annotation of some tubulin family members. Moreover, by computational techniques we delineate at least three putative regulatory regions located both in the upstream region and in the first intron containing putative binding sites for Forkhead and Nkx transcription factor families.

Assembly of 500,000 inter-specific catfish expressed sequence tags and large scale gene-associated marker development for whole genome association studies

Twelve cDNA libraries from two species of catfish have been sequenced, resulting in the generation of nearly 500,000 ESTs.

Background

Through the Community Sequencing Program, a catfish EST sequencing project was carried out through a collaboration between the catfish research community and the Department of Energy's Joint Genome Institute. Prior to this project, only a limited EST resource from catfish was available for the purpose of SNP identification.

Results

A total of 438,321 quality ESTs were generated from 8 channel catfish (Ictalurus punctatus) and 4 blue catfish (Ictalurus furcatus) libraries, bringing the number of catfish ESTs to nearly 500,000. Assembly of all catfish ESTs resulted in 45,306 contigs and 66,272 singletons. Over 35% of the unique sequences had significant similarities to known genes, allowing the identification of 14,776 unique genes in catfish. Over 300,000 putative SNPs have been identified, of which approximately 48,000 are high-quality SNPs identified from contigs with at least four sequences and the minor allele presence of at least two sequences in the contig. The EST resource should be valuable for identification of microsatellites, genome annotation, large-scale expression analysis, and comparative genome analysis.

Conclusions

This project generated a large EST resource for catfish that captured the majority of the catfish transcriptome. The parallel analysis of ESTs from two closely related Ictalurid catfishes should also provide powerful means for the evaluation of ancient and recent gene duplications, and for the development of high-density microarrays in catfish. The inter- and intra-specific SNPs identified from all catfish EST dataset assembly will greatly benefit the catfish introgression breeding program and whole genome association studies.

MytiBase: a knowledgebase of mussel (M. galloprovincialis) transcribed sequences

BACKGROUND

Although bivalves are among the most studied marine organisms due to their ecological role, economic importance and use in pollution biomonitoring, very little information is available on the genome sequences of mussels. This study reports the functional analysis of a large-scale Expressed Sequence Tag (EST) sequencing from different tissues of Mytilus galloprovincialis (the Mediterranean mussel) challenged with toxic pollutants, temperature and potentially pathogenic bacteria.

RESULTS

We have constructed and sequenced seventeen cDNA libraries from different Mediterranean mussel tissues: gills, digestive gland, foot, anterior and posterior adductor muscle, mantle and haemocytes. A total of 24,939 clones were sequenced from these libraries generating 18,788 high-quality ESTs which were assembled into 2,446 overlapping clusters and 4,666 singletons resulting in a total of 7,112 non-redundant sequences. In particular, a high-quality normalized cDNA library (Nor01) was constructed as determined by the high rate of gene discovery (65.6%). Bioinformatic screening of the non-redundant M. galloprovincialis sequences identified 159 microsatellite-containing ESTs. Clusters, consensuses, related similarities and gene ontology searches have been organized in a dedicated, searchable database http://mussel.cribi.unipd.it.

CONCLUSION

We defined the first species-specific catalogue of M. galloprovincialis ESTs including 7,112 unique transcribed sequences. Putative microsatellite markers were identified. This annotated catalogue represents a valuable platform for expression studies, marker validation and genetic linkage analysis for investigations in the biology of Mediterranean mussels.

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.

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.

Ca(2+) signaling occurs via second messenger release from intraorganelle synthesis sites

Cyclic ADP-ribose is an important Ca(2+)-mobilizing cytosolic messenger synthesized from beta-NAD(+) by ADP-ribosyl cyclases (ARCs). However, the focus upon ectocellular mammalian ARCs (CD38 and CD157) has led to confusion as to how extracellular enzymes generate intracellular messengers in response to stimuli. We have cloned and characterized three ARCs in the sea urchin egg and found that endogenous ARCbeta and ARCgamma are intracellular and located within the lumen of acidic, exocytotic vesicles, where they are optimally active. Intraorganelle ARCs are shielded from cytosolic substrate and targets by the organelle membrane, but this barrier is circumvented by nucleotide transport. We show that a beta-NAD(+) transporter provides ARC substrate that is converted luminally to cADPR, which, in turn, is shuttled out to the cytosol via a separate cADPR transporter. Moreover, nucleotide transport is integral to ARC activity physiologically because three transport inhibitors all inhibited the fertilization-induced Ca(2+) wave that is dependent upon cADPR. This represents a novel signaling mechanism whereby an extracellular stimulus increases the concentration of a second messenger by promoting messenger transport from intraorganelle synthesis sites to the cytosol.

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.

Distinct profiles of expressed sequence tags during intestinal regeneration in the sea cucumber Holothuria glaberrima

Repair and regeneration are key processes for tissue maintenance, and their disruption may lead to disease states. Little is known about the molecular mechanisms that underline the repair and regeneration of the digestive tract. The sea cucumber Holothuria glaberrima represents an excellent model to dissect and characterize the molecular events during intestinal regeneration. To study the gene expression profile, cDNA libraries were constructed from normal, 3-day, and 7-day regenerating intestines of H. glaberrima. Clones were randomly sequenced and queried against the nonredundant protein database at the National Center for Biotechnology Information. RT-PCR analyses were made of several genes to determine their expression profile during intestinal regeneration. A total of 5,173 sequences from three cDNA libraries were obtained. About 46.2, 35.6, and 26.2% of the sequences for the normal, 3-days, and 7-days cDNA libraries, respectively, shared significant similarity with known sequences in the protein database of GenBank but only present 10% of similarity among them. Analysis of the libraries in terms of functional processes, protein domains, and most common sequences suggests that a differential expression profile is taking place during the regeneration process. Further examination of the expressed sequence tag dataset revealed that 12 putative genes are differentially expressed at significant level (R > 6). Experimental validation by RT-PCR analysis reveals that at least three genes (unknown C-4677-1, melanotransferrin, and centaurin) present a differential expression during regeneration. These findings strongly suggest that the gene expression profile varies among regeneration stages and provide evidence for the existence of differential gene expression.

Analysis of the immune-inducible genes of Plutella xylostella using expressed sequence tags and cDNA microarray

In the present study, the complex gene expression responses of Plutella xylostella to microbial challenges and injury were surveyed using a newly constructed expressed sequence tag (EST) clone collection and cDNA microarray analysis. A total of 1132 P. xylostella ESTs were cloned, annotated and categorized by their putative functions; these included proteases, protease inhibitors, recognition molecules and anti-microbial peptides. GeneOntology revealed that 4% of the P. xylostella ESTs corresponded to immunity-related genes potentially involved in innate immunity. We then used microarray analysis to identify 44 genes that were differentially expressed with at least a two-fold expression difference in P. xylostella before and after pathogen challenge. Together, our EST categorization and microarray profiling analyses allowed us to identify 70 genes that should be considered candidate immune response genes, providing important new insights into the molecular events that occur during the innate immune response in P. xylostella.

A genomic view of the sea urchin nervous system

The sequencing of the Strongylocentrotus purpuratus genome provides a unique opportunity to investigate the function and evolution of neural genes. The neurobiology of sea urchins is of particular interest because they have a close phylogenetic relationship with chordates, yet a distinctive pentaradiate body plan and unusual neural organization. Orthologues of transcription factors that regulate neurogenesis in other animals have been identified and several are expressed in neurogenic domains before gastrulation indicating that they may operate near the top of a conserved neural gene regulatory network. A family of genes encoding voltage-gated ion channels is present but, surprisingly, genes encoding gap junction proteins (connexins and pannexins) appear to be absent. Genes required for synapse formation and function have been identified and genes for synthesis and transport of neurotransmitters are present. There is a large family of G-protein-coupled receptors, including 874 rhodopsin-type receptors, 28 metabotropic glutamate-like receptors and a remarkably expanded group of 161 secretin receptor-like proteins. Absence of cannabinoid, lysophospholipid and melanocortin receptors indicates that this group may be unique to chordates. There are at least 37 putative G-protein-coupled peptide receptors and precursors for several neuropeptides and peptide hormones have been identified, including SALMFamides, NGFFFamide, a vasotocin-like peptide, glycoprotein hormones and insulin/insulin-like growth factors. Identification of a neurotrophin-like gene and Trk receptor in sea urchin indicates that this neural signaling system is not unique to chordates. Several hundred chemoreceptor genes have been predicted using several approaches, a number similar to that for other animals. Intriguingly, genes encoding homologues of rhodopsin, Pax6 and several other key mammalian retinal transcription factors are expressed in tube feet, suggesting tube feet function as photosensory organs. Analysis of the sea urchin genome presents a unique perspective on the evolutionary history of deuterostome nervous systems and reveals new approaches to investigate the development and neurobiology of sea urchins.

Analysis of cytoskeletal and motility proteins in the sea urchin genome assembly

The sea urchin embryo is a classical model system for studying the role of the cytoskeleton in such events as fertilization, mitosis, cleavage, cell migration and gastrulation. We have conducted an analysis of gene models derived from the Strongylocentrotus purpuratus genome assembly and have gathered strong evidence for the existence of multiple gene families encoding cytoskeletal proteins and their regulators in sea urchin. While many cytoskeletal genes have been cloned from sea urchin with sequences already existing in public databases, genome analysis reveals a significantly higher degree of diversity within certain gene families. Furthermore, genes are described corresponding to homologs of cytoskeletal proteins not previously documented in sea urchins. To illustrate the varying degree of sequence diversity that exists within cytoskeletal gene families, we conducted an analysis of genes encoding actins, specific actin-binding proteins, myosins, tubulins, kinesins, dyneins, specific microtubule-associated proteins, and intermediate filaments. We conducted ontological analysis of select genes to better understand the relatedness of urchin cytoskeletal genes to those of other deuterostomes. We analyzed developmental expression (EST) data to confirm the existence of select gene models and to understand their differential expression during various stages of early development.

Sp-Smad2/3 mediates patterning of neurogenic ectoderm by nodal in the sea urchin embryo

Nodal functions in axis and tissue specification during embryogenesis. In sea urchin embryos, Nodal is crucial for specification of oral ectoderm and is thought to pattern neurogenesis in the animal plate. To determine if Nodal functions directly in suppressing neuron differentiation we have prepared mutant forms of Sp-Smad2/3. Expressing an activated form produces embryos similar to embryos overexpressing Nodal, but with fewer neurons. In chimeras in which Nodal is suppressed, cells expressing activated Sp-Smad2/3 form oral ectoderm, but not neurons. In embryos with vegetal signaling blocked, neurons do not form if activated Smad2/3 is co-expressed. Expression of dominant negative mutants produces embryos identical to those resulting from blocking Nodal expression. In chimeras overexpressing Nodal, cells expressing dominant negative Sp-Smad2/3 form aboral ectoderm and give rise to neurons. In permanent blastula chimeras dominant negative Sp-Smad2/3 is able to suppress the effects of Nodal permitting neuron differentiation. In these chimeras Nodal expression in one half suppresses neural differentiation across the interface. Anti-phospho-Smad3 reveals that the cells adjacent to cells expressing Nodal have nuclear immunoreactivity. We conclude Sp-Smad2/3 is a component of the Nodal signaling pathway in sea urchins and that Nodal diffuses short distances to suppress neural differentiation.

Molecular diversity and evolution of the Siglec family of cell-surface lectins

Sialic acids are a family of acidic sugars with a 9-carbon backbone, prominently expressed in animals of deuterostome lineage. Siglecs are the largest family of vertebrate endogenous receptors that recognize glycoconjugates containing sialic acids. Although a few Siglecs are well-conserved throughout vertebrate evolution and show similar binding preference regardless of the species of origin, most others, particularly the CD33-related subfamily of Siglecs, show marked inter-species differences in repertoire, sequence, and binding preference. The diversification of CD33-related Siglecs may be driven by direct competition against pathogens, and/or by necessity to catch up with the changing landscape of endogenous glycans, which may in turn be changing to escape exploitation by other pathogens.

Intermediary metabolism in sea urchin: the first inferences from the genome sequence

The genome sequence of the purple sea urchin Strongylocentrotus purpuratus recently became available. We report the results of functional annotation and initial analysis of more than 2300 proteins predicted to be involved in metabolite transport and enzymatic conversion in sea urchin. The comparison of various reconstructed biosynthetic and catabolic pathways in sea urchin to those known in other genomes suggests the overall similarity of the sea urchin metabolism to that of the vertebrates, with relatively small but non-trivial differences from both vertebrates and protostomes. There are several examples of two parallel, non-orthologous solutions for the same molecular function in sea urchin, in contrast with the other completely sequenced metazoans that tend to contain just one version of the same function. There are also genes that appear to be close phylogenetic neighbors of plant or bacterial homologs, as opposed to homologs in other Metazoa. The evolutionary and functional significance of these variations is discussed.

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

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

Expression of Nox genes in rat organs, mouse oocytes, and sea urchin eggs

Degenerate primers were designed to isolate new homologs of Nox family genes in rat organs and sea urchin eggs. The primers were capable of amplifying Nox1, Nox2, Nox3, Nox4, Duox1 and Duox2 but not Nox5, and failed to isolate novel homologs in rat. However, a novel homolog (named as Nox-U1) was identified in sea urchin eggs. In the most conserved region (amino acid 336--417 in human Nox2) Nox-U1 has the highest identity with Nox2, which appears to be abundant in mouse oocytes. However, phylogenetic analysis of the entire sequence has revealed that Nox-U1 is closer to Nox4 or Nox5 than Nox2 or Nox3. Histidine residues assumed to be responsible for heme ligation, motifs for FAD- and NADPH-binding, and two asparagine-linked glycosylation sites are conserved.

Expression of Spgatae, the Strongylocentrotus purpuratus ortholog of vertebrate GATA4/5/6 factors

Spgatae is the sea urchin ortholog of the vertebrate gata4/5/6 genes, as confirmed by phylogenetic analysis. The accumulation of Spgatae transcripts during embryonic development and the spatial pattern of expression are reported here. Expression was first detected in the 15 h blastula. The number of Spgatae RNA molecules increases steadily during blastula stages, with expression peaking during gastrulation. After gastrulation is complete, the level of expression decreases until the end of embryogenesis. Whole mount in situ hybridization showed that Spgatae transcripts were first detected in a ring of prospective mesoderm cells in the vegetal plate. Spgatae expression then expands to include the entire vegetal plate at the mesenchyme blastula stage. During gastrulation Spgatae is expressed at the blastopore, and at prism stage strongly in the hindgut and midgut but not foregut, and also in mesoderm cells at the tip of the archenteron. Towards the end of embryogenesis, expression in the hindgut decreases. The terminal pattern of expression is in midgut plus coelomic pouches.

Evolution of nuclearly encoded mitochondrial genes in Metazoa

All Metazoan nuclear genomes underwent a continuous process of both complete and partial genetic material gain and loss. The forces modulating these events are also subject to the strict interaction between nuclear and mitochondrial (mt) genome. In this context we investigate the evolution of nuclear genes encoding proteins which target the mitochondrion, with a particular attention to genes involved in oxidative phosphorylation (OXPHOS), one of the most ancient and conserved functions. To examine thoroughly the evolutionary strategies that preserve OXPHOS and coordinate the two cellular genomes, a comparative analysis has been carried out for 78 OXPHOS gene families in several Metazoa (insects, tunicates, fishes and mammals). We demonstrate that the duplication rate of OXPHOS genes increases passing from invertebrates to vertebrates consistently with the total increase in genome size, but all species are prone to negatively select OXPHOS duplicates compared to the general trend of nuclear gene families. These results are consistent with the 'balance hypothesis' and, at least in insects, the expression of duplicate genes is low and strongly testis-biased.

A first generation physical map of the medaka genome in BACs essential for positional cloning and clone-by-clone based genomic sequencing

In order to realize the full potential of the medaka as a model system for developmental biology and genetics, characterized genomic resources need to be established, culminating in the sequence of the medaka genome. To facilitate the map-based cloning of genes underlying induced mutations and to provide templates for clone-based genomic sequencing, we have created a first-generation physical map of the medaka genome in bacterial artificial chromosome (BAC) clones. In particular, we exploited the synteny to the closely related genome of the pufferfish, Takifugu rubripes, by marker content mapping. As a first step, we clustered 103,144 public medaka EST sequences to obtain a set of 21,121 non-redundant sequence entities. Avoiding oversampling of gene-dense regions, 11,254 of EST clusters were successfully matched against the draft sequence of the fugu genome, and 2363 genes were selected for the BAC map project. We designed 35mer oligonucleotide probes from the selected genes and hybridized them against 64,500 BAC clones of strains Cab and Hd-rR, representing 14-fold coverage of the medaka genome. Our data set is further supplemented with 437 results generated from PCR-amplified inserts of medaka cDNA clones and BAC end-fragment markers. Our current, edited, first generation medaka BAC map consists of 902 map segments that cover about 74% of the medaka genome. The map contains 2721 markers. Of these, 2534 are from expressed sequences, equivalent to a non-redundant set of 2328 loci. The 934 markers (724 different) are anchored to the medaka genetic map. Thus, genetic map assignments provide immediate access to underlying clones and contigs, simplifying molecular access to candidate gene regions and their characterization.

Development of a citrus genome-wide EST collection and cDNA microarray as resources for genomic studies

A functional genomics project has been initiated to approach the molecular characterization of the main biological and agronomical traits of citrus. As a key part of this project, a citrus EST collection has been generated from 25 cDNA libraries covering different tissues, developmental stages and stress conditions. The collection includes a total of 22,635 high-quality ESTs, grouped in 11,836 putative unigenes, which represent at least one third of the estimated number of genes in the citrus genome. Functional annotation of unigenes which have Arabidopsis orthologues (68% of all unigenes) revealed gene representation in every major functional category, suggesting that a genome-wide EST collection was obtained. A Citrus clementina Hort. ex Tan. cv. Clemenules genomic library, that will contribute to further characterization of relevant genes, has also been constructed. To initiate the analysis of citrus transcriptome, we have developed a cDNA microarray containing 12,672 probes corresponding to 6875 putative unigenes of the collection. Technical characterization of the microarray showed high intra- and inter-array reproducibility, as well as a good range of sensitivity. We have also validated gene expression data achieved with this microarray through an independent technique such as RNA gel blot analysis.

Animal phylogenomics: multiple interspecific genome comparisons

The utility of DNA sequence information for phylogenetics and phylogeography is now well known. Rather than attempt to summarize studies addressing this well-demonstrated utility, this chapter focuses on fundamental approaches and techniques that implement the collection of DNA sequence data for comparative phylogenetic purposes in a genomic context (phylogenomics). Whole genome sequencing approaches have changed the way we think about phylogenetics and have opened the way for new perspectives on "old" phylogenetics concerns. Some of these concerns are which gene regions to use and how much sequence information is needed for robust phylogenetic inference. Whole genome sequences of a few animal model organisms have gone a long way to implement approaches to better understand these important phylogenetic concerns. This chapter also addresses how genomics has made it more important for a clear understanding of orthology of gene regions in comparative biology. Finally, genome-enabled technologies that are affecting comparative biology are also discussed.

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.

GOblet: a platform for Gene Ontology annotation of anonymous sequence data

GOblet is a comprehensive web server application providing the annotation of anonymous sequence data with Gene Ontology (GO) terms. It uses a variety of different protein databases (human, murines, invertebrates, plants, sp-trembl) and their respective GO mappings. The user selects the appropriate database and alignment threshold and thereafter submits single or multiple nucleotide or protein sequences. Results are shown in different ways, e.g. as survey statistics for the main GO categories for all sequences or as detailed results for each single sequence that has been submitted. In its newest version, GOblet allows the batch submission of sequences and provides an improved display of results with the aid of Java applets. All output data, together with the Java applet, are packed to a downloadable archive for local installation and analysis. GOblet can be accessed freely at http://goblet.molgen.mpg.de.

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.

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.