Derivedness Index for Estimating Degree of Phenotypic Evolution of Embryos: A Study of Comparative Transcriptomic Analyses of Chordates and Echinoderms.

Species retaining ancestral features, such as species called living fossils, are often regarded as less derived than their sister groups, but such discussions are usually based on qualitative enumeration of conserved traits. This approach creates a major barrier, especially when quantifying the degree of phenotypic evolution or degree of derivedness, since it focuses only on commonly shared traits, and newly acquired or lost traits are often overlooked. To provide a potential solution to this problem, especially for inter-species comparison of gene expression profiles, we propose a new method named "derivedness index" to quantify the degree of derivedness. In contrast to the conservation-based approach, which deals with expressions of commonly shared genes among species being compared, the derivedness index also considers those that were potentially lost or duplicated during evolution. By applying our method, we found that the gene expression profiles of penta-radial phases in echinoderm tended to be more highly derived than those of the bilateral phase. However, our results suggest that echinoderms may not have experienced much larger modifications to their developmental systems than chordates, at least at the transcriptomic level. In vertebrates, we found that the mid-embryonic and organogenesis stages were generally less derived than the earlier or later stages, indicating that the conserved phylotypic period is also less derived. We also found genes that potentially explain less derivedness, such as Hox genes. Finally, we highlight technical concerns that may influence the measured transcriptomic derivedness, such as read depth and library preparation protocols, for further improvement of our method through future studies. We anticipate that this index will serve as a quantitative guide in the search for constrained developmental phases or processes.

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

BACKGROUND: Proteomic studies of skeletal proteins have revealed large, complex mixtures of proteins occluded within the mineral. Many skeletal proteomes contain rapidly evolving proteins with repetitive domains, further complicating our understanding. In echinoderms, proteomic analysis of the skeletal proteomes of mineralized tissues of the sea urchin Strongylocentrotus purpuratus prominently featured spicule matrix proteins with repetitive sequences linked to a C-type lectin domain. A comparative study of the brittle star Ophiocoma wendtii skeletal proteome revealed an order of magnitude fewer proteins containing C-type lectin domains. A number of other proteins conserved in the skeletons of the two groups were identified. Here we report the complete skeletal proteome of the sea star Patiria miniata and compare it to that of the other echinoderm groups. RESULTS: We have identified eighty-five proteins in the P. miniata skeletal proteome. Forty-two percent of the proteins were determined to be homologous to proteins found in the S. purpuratus skeletal proteomes. An additional 34 % were from similar functional classes as proteins in the urchin proteomes. Thirteen percent of the P. miniata proteins had homologues in the O. wendtii skeletal proteome with an additional 29% showing similarity to brittle star skeletal proteins. The P. miniata skeletal proteome did not contain any proteins with C-lectin domains or with acidic repetitive regions similar to the sea urchin or brittle star spicule matrix proteins. MSP130 proteins were also not found. We did identify a number of proteins homologous between the three groups. Some of the highly conserved proteins found in echinoderm skeletons have also been identified in vertebrate skeletons. CONCLUSIONS: The presence of proteins conserved in the skeleton in three different echinoderm groups indicates these proteins are important in skeleton formation. That a number of these proteins are involved in skeleton formation in vertebrates suggests a common origin for some of the fundamental processes co-opted for skeleton formation in deuterostomes. The proteins we identify suggest transport of proteins and calcium via endosomes was co-opted to this function in a convergent fashion. Our data also indicate that modifications to the process of skeleton formation can occur through independent co-option of proteins following species divergence as well as through domain shuffling.

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

BACKGROUND: While formation of mineralized tissue is characteristic of many animal taxa, the proteins that interact with mineral are diverse and appear in many cases to be of independent origin. Extracellular matrix proteins involved in mineralization do share some common features. They tend to be disordered, secreted proteins with repetitive, low complexity. The genes encoding these proteins are often duplicated and undergo concerted evolution, further diversifying the repetitive domains. This makes it difficult to identify mineralization genes and the proteins they encode using bioinformatics techniques. Here we describe the use of proteomics to identify mineralization genes in an ophiuroid echinoderm, Ophiocoma wendtii (O. wendtii). RESULTS: We have isolated the occluded proteins within the mineralized tissue of the brittle star Ophiocoma wendtii. The proteins were analyzed both unfractionated and separated on SDS-PAGE gels. Each slice was analyzed using mass spectroscopy and the amino acid sequence of the most prevalent peptides was obtained. This was compared to both an embryonic transcriptome from the gastrula stage when skeleton is being formed and a tube foot (an adult mineralized tissue) transcriptome. Thirty eight proteins were identified which matched known proteins or protein domains in the NCBI databases. These include C-type lectins, ECM proteins, Kazal-type protease inhibitors, matrix metalloproteases as well as more common cellular proteins. Many of these are similar to those found in the sea urchin Strongylocentrotus purpuratus (S. purpuratus) skeleton. We did not, however, identify clear homologs to the sea urchin spicule matrix proteins, and the number of C-type lectin containing genes was much reduced compared to sea urchins. Also notably absent was MSP-130. CONCLUSIONS: Our results show an overall conservation of the types of proteins found in the mineralized tissues of two divergent groups of echinoderms, as well as in mineralized tissues in general. However, the extensive gene duplication and concerted evolution seen in the spicule matrix proteins found in the sea urchin skeleton was not observed in the brittle star.

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.

Tracing misalignment and unequal crossing over during evolution of the repetitive region of the SM50 gene in sea urchins.

The location of misalignment and unequal crossover involved in concerted evolution of tandemly repetitive sequences is difficult to document owing to the homogeneity of sequences that are subject to this process. However, the repetitive domain of the SM50 gene in sea urchins contains variation, within the gene itself, between alleles, and between species that has allowed us to determine where misalignment and unequal crossing over occurred during evolution of this gene. We have therefore analyzed the SM50 repeat regions in a variety of species to determine where recent changes in repeat numbers have occurred, and from this have deduced the mechanisms that lead to these changes. We next tried to determine whether recent misalignment and unequal crossover has produced allelic variation in current populations of sea urchins. We found SM50 alleles within three species that have different numbers of repeats. This marks the first reported documentation of allelic variation in the number of repeats in the SM50 gene. We also show how a single unequal crossover event could have produced the allelic variation. We have found that substitutions and small deletions in the sequences within the repeats can substantially affect how misalignment occurs, resulting in different patterns of repeats after concerted evolution.

Sea urchin metalloproteases: a genomic survey of the BMP-1/tolloid-like, MMP and ADAM families.

Analysis of the Strongylocentrotus purpuratus genome has revealed approximately 240 metalloprotease genes, and they represent all 23 families expressed in vertebrates. EST/cDNA sequencing and microarray analysis show that nearly 70% are represented in embryo RNA. Among them are many metalloproteases with demonstrated developmental roles in other systems-BMP-1/TLD (tolloid) (astacins), MMPs (matrix metalloproteases) and the ADAMs (disintegrin/metalloproteases). The developmental functions of these kinds of metalloproteases include modifying the extracellular matrix, regulating signaling pathways or modulating cellular adhesive properties. The unexpectedly large number of BMP-1/TLD-like protease genes (23) results primarily from expansion of a set encoding an unusual domain conserved in structure and primary sequence only in nematode astacins. Such proteases may have interesting developmental functions because the expression patterns of several are highly regulated along the primary axis at times when cell differentiation and morphogenesis begin. The size of the sea urchin MMP family and the clustered arrangement of many of its members are similar to vertebrates, but phylogenetic analyses suggest that different ancestral genes were independently amplified in sea urchins and vertebrates. One expansion appears to be genes encoding MMPs that have putative transmembrane domains and may be membrane-tethered (MT). Interestingly, the genes encoding TIMPs, inhibitors of MMPs, have also been amplified and the 10 genes are tandemly arranged in a single cluster. In contrast, there are fewer ADAM and ADAMTS genes in sea urchins, but they represent all but one of the chordate-specific groups. The genome sequence now opens the door to experimental manipulations designed to understand how modulation of the extracellular environment affects development.

Overexpression and mechanistic characterization of blastula protease 10, a metalloprotease involved in sea urchin embryogenesis and development.

Blastula protease 10 (BP10) is a metalloenzyme involved in sea urchin embryogenesis, which has been assigned to the astacin family of zinc-dependent endopeptidases. It shows greatest homology with the mammalian tolloid-like genes and contains conserved structural motifs consistent with astacin, tolloid, and bone morphogenetic protein 1. Astacin, a crustacean digestive enzyme, has been proposed to carry out hydrolysis via a metal-centered mechanism that involves a metal-coordinated "tyrosine switch." It has not been determined if the more structurally complex members of this family involved in eukaryotic development share this mechanism. The recombinant BP10 has been overexpressed in Escherichia coli, its metalloenzyme nature has been confirmed, and its catalytic properties have been characterized through kinetic studies. BP10 shows significant hydrolysis toward gelatin both in its native zinc-containing form and copper derivative. The copper derivative of BP10 shows a remarkable 960% rate acceleration toward the hydrolysis of the synthetic substrate N-benzoyl-arginine-p-nitroanilide when compared with the zinc form. The enzyme also shows calcium-dependent activation. These are the first thorough mechanistic studies reported on BP10 as a representative of the more structurally complex members of astacin-type enzymes in deuterostomes, which can add supporting data to corroborate the metal-centered mechanism proposed for astacin and the role of the coordinated Tyr. We have demonstrated the first mechanistic study of a tolloid-related metalloenzyme involved in sea urchin embryogenesis.

Identification of cis-regulatory elements involved in transcriptional regulation of the sea urchin SpFoxB gene.

The SpFoxB gene is transiently expressed first in the mesoderm, then in the endoderm and oral ectoderm during sea urchin gastrulation. Perturbations of a number of proteins involved in endomesoderm specification have been shown to alter the mRNA levels of SpFoxB, but the cis-regulatory elements required for expression of SpFoxB have not been examined. In order to investigate this, we have screened the SpFoxB gene for sequences that can drive its expression. Both positive and negative cis-regulatory elements were found to be present. An enhancer was found that contains four GATA sites and four YY1 sites clustered within 210 base pairs (bp), as well as three lef/tcf binding sites. Electrophoretic mobility shifts indicate that the lef/tcf sites bind a complex of proteins that include beta-catenin in early cleavage, but not during subsequent stages of development. The GATA and YY1 sites bind nuclear proteins prior to SpFoxB transcription, and this binding diminishes coincident with cessation of transcription. Deletion of the GATA/YY1 sites causes a significant decrease in transcription. The DNA binding site of the SpFoxB protein has been determined, and Fox binding sites are found within the 5' UTR of SpFoxB.

Development of calcareous skeletal elements in invertebrates.

Most metazoans require skeletal support systems. While the formation of bones and teeth in vertebrates has been well studied, endo- and exoskeleton development of non-vertebrates, especially calcification during terminal differentiation, has been neglected. Biomineralization of skeletons in invertebrates presents interesting research opportunities. We undertake here to survey some of the better understood examples of skeletal development in selected invertebrates. The differentiation of the skeletal spicules of euechinoid larvae and other non-vertebrate deuterostomes, the shells of molluscs, and the calcification of crustacean carapaces are surveyed. The diversity of these different kinds of animals and our present limited understanding make it difficult to identify unifying themes, but there certainly are unifying questions: How is the mineral precursor secreted? What is the nature of the interaction of mineral with the matrix proteins of the skeleton? Is there any conservation of protein domains in matrix proteins found in skeletal elements from different phyla? Are there common strategies in the development of organs that form mineralized structures?

Injection of myo-inositol reverses the effects of lithium on sea urchin blastomeres.

Lithium is known to cause sea urchin blastomeres destined to give rise to epithelium rather than to differentiate into gut or skeleton. While it has been proposed that lithium alters development by interfering with the inositol-tris phosphate-protein kinase C (IP3 -PKC) signaling pathway, the mechanism of action of lithium in sea urchins has remained elusive. Here we describe experiments that examine the hypothesis that lithium exerts its effect on sea urchin embryos via the IP3 -PKC pathway. We make use of methods developed to isolate epithelial precursor cells from the animal hemisphere of cleavage 16-cell stage embryos. Pairs of cells were isolated and one of each pair was injected with either myo-inositol or its inactive isomer, epi-inositol. Rhodamine dextran was co-injected as a lineage tracer to follow the fate of injected cells. We demonstrate that injected myo-inositol, but not epi-inositol, can reverse the effects of lithium on sea urchin blastomeres. This is direct evidence that lithium affects the IP3 -PKC pathway in sea urchins, and that this pathway plays an important role in cell fate determination.