Possible cooption of a VEGF-driven tubulogenesis program for biomineralization in echinoderms.

Biomineralization is the process by which living organisms use minerals to form hard structures that protect and support them. Biomineralization is believed to have evolved rapidly and independently in different phyla utilizing preexisting components. The mechanistic understanding of the regulatory networks that drive biomineralization and their evolution is far from clear. Sea urchin skeletogenesis is an excellent model system for studying both gene regulation and mineral uptake and deposition. The sea urchin calcite spicules are formed within a tubular cavity generated by the skeletogenic cells controlled by vascular endothelial growth factor (VEGF) signaling. The VEGF pathway is essential for biomineralization in echinoderms, while in many other phyla, across metazoans, it controls tubulogenesis and vascularization. Despite the critical role of VEGF signaling in sea urchin spiculogenesis, the downstream program it activates was largely unknown. Here we study the cellular and molecular machinery activated by the VEGF pathway during sea urchin spiculogenesis and reveal multiple parallels to the regulation of vertebrate vascularization. Human VEGF rescues sea urchin VEGF knockdown, vesicle deposition into an internal cavity plays a significant role in both systems, and sea urchin VEGF signaling activates hundreds of genes, including biomineralization and interestingly, vascularization genes. Moreover, five upstream transcription factors and three signaling genes that drive spiculogenesis are homologous to vertebrate factors that control vascularization. Overall, our findings suggest that sea urchin spiculogenesis and vertebrate vascularization diverged from a common ancestral tubulogenesis program, broadly adapted for vascularization and specifically coopted for biomineralization in the echinoderm phylum.

Spinochrome Identification and Quantification in Pacific Sea Urchin Shells, Coelomic Fluid and Eggs Using HPLC-DAD-MS.

The high-performance liquid chromatography method coupled with diode array and mass spectrometric detector (HPLC-DAD-MS) method for quinonoid pigment identification and quantification in sea urchin samples was developed and validated. The composition and quantitative ratio of the quinonoid pigments of the shells of 16 species of sea urchins, collected in the temperate (Sea of Japan) and tropical (South-China Sea) climatic zones of the Pacific Ocean over several years, were studied. The compositions of the quinonoid pigments of sea urchins Maretia planulata, Scaphechinus griseus, Laganum decagonale and Phyllacanthus imperialis were studied for the first time. A study of the composition of the quinonoid pigments of the coelomic fluid of ten species of sea urchins was conducted. The composition of quinonoid pigments of Echinarachnius parma jelly-like egg membrane, of Scaphechinus mirabilis developing embryos and pluteus, was reported for the first time. In the case of Scaphechinus mirabilis, we have shown that the compositions of pigment granules of the shell epidermis, coelomic fluid, egg membrane, developing embryos and pluteus are different, which should enable a fuller understanding of the functions of pigments at different stages of life.

The first complete mitochondrial genome of the sand dollar Sinaechinocyamus mai (Echinoidea: Clypeasteroida).

Morphologic and molecular data often lead to different hypotheses of phylogenetic relationships. Such incongruence has been found in the echinoderm class Echinoidea. In particular, the phylogenetic status of the order Clypeasteroida is not well resolved. Complete mitochondrial genomes are currently available for 29 echinoid species, but no clypeasteroid had been sequenced to date. DNA extracted from a single live individual of Sinaechinocyamus mai was sequenced with 10× Genomics technology. This first complete mitochondrial genome (mitogenome) for the order Clypeasteroida is 15,756 base pairs in length. Phylogenomic analysis based on 34 ingroup taxa belonging to nine orders of the class Echinoidea show congruence between our new genetic inference and published trees based on morphologic characters, but also includes some intriguing differences that imply the need for additional investigation.

A new ophiocistioid with soft-tissue preservation from the Silurian Herefordshire Lagerstätte, and the evolution of the holothurian body plan.

Reconstructing the evolutionary assembly of animal body plans is challenging when there are large morphological gaps between extant sister taxa, as in the case of echinozoans (echinoids and holothurians). However, the inclusion of extinct taxa can help bridge these gaps. Here we describe a new species of echinozoan, Sollasina cthulhu, from the Silurian Herefordshire Lagerstätte, UK. Sollasina cthulhu belongs to the ophiocistioids, an extinct group that shares characters with both echinoids and holothurians. Using physical-optical tomography and computer reconstruction, we visualize the internal anatomy of S. cthulhu in three dimensions, revealing inner soft tissues that we interpret as the ring canal, a key part of the water vascular system that was previously unknown in fossil echinozoans. Phylogenetic analyses strongly suggest that Sollasina and other ophiocistioids represent a paraphyletic group of stem holothurians, as previously hypothesized. This allows us to reconstruct the stepwise reduction of the skeleton during the assembly of the holothurian body plan, which may have been controlled by changes in the expression of biomineralization genes.

MITF: an evolutionarily conserved transcription factor in the sea urchin Paracentrotus lividus.

Microphthalmia-associated transcription factor (MITF) is a member of MYC superfamily, associated with melanocyte cells, as it was discovered in depigmented mice. However, over the last years it was found to be involved in many cellular signaling pathways, among which oncogenesis, osteoclast differentiation, and stress response. In mammals, Mitf gene mutations can cause diverse syndromes affecting pigmentation of eyes or skin, bone defects and melanomas. As MITF protein homologs were also found in some invertebrates, we have isolated and characterized the MITF cDNAs from the sea urchin Paracentrotus lividus, referred to as Pl-Mitf. The in silico study of the secondary and tertiary structure of Pl-Mitf protein showed high conserved regions mostly lying in the DNA binding domain. To understand the degree of evolutionary conservation of MITF, a phylogenetic analysis was performed comparing the Pl-Mitf deduced protein with proteins from different animal species. Moreover, the analysis of temporal and spatial expression pattern of Pl-Mitf mRNA showed that it was expressed from the onset of gastrulation of the sea urchin embryo to the pluteus larva, specifically in primary mesenchymes cells (PMCs), the sea urchin skeletogenic cells, and in the forming archenteron, the larval gut precursor. In silico protein-protein interactions analysis was used to understand the association of MITF with other proteins. Our results put in evidence the conservation of the MITF protein among vertebrates and invertebrates and may provide new perspectives on the pathways underlying sea urchin development, even if further functional analyses are needed.

A phylogenomic resolution of the sea urchin tree of life.

BACKGROUND: Echinoidea is a clade of marine animals including sea urchins, heart urchins, sand dollars and sea biscuits. Found in benthic habitats across all latitudes, echinoids are key components of marine communities such as coral reefs and kelp forests. A little over 1000 species inhabit the oceans today, a diversity that traces its roots back at least to the Permian. Although much effort has been devoted to elucidating the echinoid tree of life using a variety of morphological data, molecular attempts have relied on only a handful of genes. Both of these approaches have had limited success at resolving the deepest nodes of the tree, and their disagreement over the positions of a number of clades remains unresolved. RESULTS: We performed de novo sequencing and assembly of 17 transcriptomes to complement available genomic resources of sea urchins and produce the first phylogenomic analysis of the clade. Multiple methods of probabilistic inference recovered identical topologies, with virtually all nodes showing maximum support. In contrast, the coalescent-based method ASTRAL-II resolved one node differently, a result apparently driven by gene tree error induced by evolutionary rate heterogeneity. Regardless of the method employed, our phylogenetic structure deviates from the currently accepted classification of echinoids, with neither Acroechinoidea (all euechinoids except echinothurioids), nor Clypeasteroida (sand dollars and sea biscuits) being monophyletic as currently defined. We show that phylogenetic signal for novel resolutions of these lineages is strong and distributed throughout the genome, and fail to recover systematic biases as drivers of our results. CONCLUSIONS: Our investigation substantially augments the molecular resources available for sea urchins, providing the first transcriptomes for many of its main lineages. Using this expanded genomic dataset, we resolve the position of several clades in agreement with early molecular analyses but in disagreement with morphological data. Our efforts settle multiple phylogenetic uncertainties, including the position of the enigmatic deep-sea echinothurioids and the identity of the sister clade to sand dollars. We offer a detailed assessment of evolutionary scenarios that could reconcile our findings with morphological evidence, opening up new lines of research into the development and evolutionary history of this ancient clade.

Effects of Nodal inhibition on development of temnopleurid sea urchins.

Adult rudiment formation in some temnopleurids begins with the formation of a cell mass that is pinched off the left ectoderm in early larval development. The cell mass forms the adult rudiment with the left coelomic pouch of the mesodermal region. However, details of the mechanisms to establish position of the cell mass are still unknown. We analyzed the inhibiting effect of Nodal, a factor for morphogenesis of the oral region and right side, for location of the cell mass, in four temnopleurids. Pulse inhibition, at least 5 min inhibition, during coelomic pouch formation allowed a cell mass to form on both sides, whereas treatments after that period did not. These results indicate that Nodal signaling controls the oral-aboral axis before gastrulation and then affects the position of the cell mass and adult rudiment up to coelomic pouch formation. They also indicate that the position of the adult rudiment under Nodal signaling pathways is conserved in temnopleurids, as adult rudiment formation is dependent on the cell mass.

Structure and function of echinoderm telomerase RNA.

Telomerase is a ribonucleoprotein (RNP) enzyme that requires an integral telomerase RNA (TR) subunit, in addition to the catalytic telomerase reverse transcriptase (TERT), for enzymatic function. The secondary structures of TRs from the three major groups of species, ciliates, fungi, and vertebrates, have been studied extensively and demonstrate dramatic diversity. Herein, we report the first comprehensive secondary structure of TR from echinoderms-marine invertebrates closely related to vertebrates-determined by phylogenetic comparative analysis of 16 TR sequences from three separate echinoderm classes. Similar to vertebrate TR, echinoderm TR contains the highly conserved template/pseudoknot and H/ACA domains. However, echinoderm TR lacks the ancestral CR4/5 structural domain found throughout vertebrate and fungal TRs. Instead, echinoderm TR contains a distinct simple helical region, termed eCR4/5, that is functionally equivalent to the CR4/5 domain. The urchin and brittle star eCR4/5 domains bind specifically to their respective TERT proteins and stimulate telomerase activity. Distinct from vertebrate telomerase, the echinoderm TR template/pseudoknot domain with the TERT protein is sufficient to reconstitute significant telomerase activity. This gain-of-function of the echinoderm template/pseudoknot domain for conferring telomerase activity presumably facilitated the rapid structural evolution of the eCR4/5 domain throughout the echinoderm lineage. Additionally, echinoderm TR utilizes the template-adjacent P1.1 helix as a physical template boundary element to prevent nontelomeric DNA synthesis, a mechanism used by ciliate and fungal TRs. Thus, the chimeric and eccentric structural features of echinoderm TR provide unparalleled insights into the rapid evolution of telomerase RNP structure and function.

Cartilage acidic protein 1, a new member of the beta-propeller protein family with amyloid propensity.

Cartilage acidic protein1 (CRTAC1) is an extracellular matrix protein of chondrogenic tissue in humans and its presence in bacteria indicate it is of ancient origin. Structural modeling of piscine CRTAC1 reveals it belongs to the large family of beta-propeller proteins that in mammals have been associated with diseases, including amyloid diseases such as Alzheimer's. In order to characterize the structure/function evolution of this new member of the beta-propeller family we exploited the unique characteristics of piscine duplicate genes Crtac1a and Crtac1b and compared their structural and biochemical modifications with human recombinant CRTAC1. We demonstrate that CRTAC1 has a beta-propeller structure that has been conserved during evolution and easily forms high molecular weight thermo-stable aggregates. We reveal for the first time the propensity of CRTAC1 to form amyloid-like structures, and hypothesize that the aggregating property of CRTAC1 may be related to its disease-association. We further contribute to the general understating of CRTAC1's and beta-propeller family evolution and function. Proteins 2017; 85:242-255. © 2016 Wiley Periodicals, Inc.

EchinoDB, an application for comparative transcriptomics of deeply-sampled clades of echinoderms.

BACKGROUND: One of our goals for the echinoderm tree of life project (http://echinotol.org) is to identify orthologs suitable for phylogenetic analysis from next-generation transcriptome data. The current dataset is the largest assembled for echinoderm phylogeny and transcriptomics. We used RNA-Seq to profile adult tissues from 42 echinoderm specimens from 24 orders and 37 families. In order to achieve sampling members of clades that span key evolutionary divergence, many of our exemplars were collected from deep and polar seas. DESCRIPTION: A small fraction of the transcriptome data we produced is being used for phylogenetic reconstruction. Thus to make a larger dataset available to researchers with a wide variety of interests, we made a web-based application, EchinoDB (http://echinodb.uncc.edu). EchinoDB is a repository of orthologous transcripts from echinoderms that is searchable via keywords and sequence similarity. CONCLUSIONS: From transcripts we identified 749,397 clusters of orthologous loci. We have developed the information technology to manage and search the loci their annotations with respect to the Sea Urchin (Strongylocentrotus purpuratus) genome. Several users have already taken advantage of these data for spin-off projects in developmental biology, gene family studies, and neuroscience. We hope others will search EchinoDB to discover datasets relevant to a variety of additional questions in comparative biology.

Ancestral state reconstruction by comparative analysis of a GRN kernel operating in echinoderms.

Diverse sampling of organisms across the five major classes in the phylum Echinodermata is beginning to reveal much about the structure and function of gene regulatory networks (GRNs) in development and evolution. Sea urchins are the most studied clade within this phylum, and recent work suggests there has been dramatic rewiring at the top of the skeletogenic GRN along the lineage leading to extant members of the euechinoid sea urchins. Such rewiring likely accounts for some of the observed developmental differences between the two major subclasses of sea urchins-cidaroids and euechinoids. To address effects of topmost rewiring on downstream GRN events, we cloned four downstream regulatory genes within the skeletogenic GRN and surveyed their spatiotemporal expression patterns in the cidaroid Eucidaris tribuloides. We performed phylogenetic analyses with homologs from other non-vertebrate deuterostomes and characterized their spatiotemporal expression by quantitative polymerase chain reaction (qPCR) and whole-mount in situ hybridization (WMISH). Our data suggest the erg-hex-tgif subcircuit, a putative GRN kernel, exhibits a mesoderm-specific expression pattern early in Eucidaris development that is directly downstream of the initial mesodermal GRN circuitry. Comparative analysis of the expression of this subcircuit in four echinoderm taxa allowed robust ancestral state reconstruction, supporting hypotheses that its ancestral function was to stabilize the mesodermal regulatory state and that it has been co-opted and deployed as a unit in mesodermal subdomains in distantly diverged echinoderms. Importantly, our study supports the notion that GRN kernels exhibit structural and functional modularity, locking down and stabilizing clade-specific, embryonic regulatory states.

Morphological and genetic analyses reveal a cryptic species complex in the echinoid Echinocardium cordatum and rule out a stabilizing selection explanation.

Preliminary analyses revealed the presence of at least five mitochondrial clades within the widespread sea urchin Echinocardium cordatum (Spatangoida). In this study, we analyzed the genetic (two mitochondrial and two nuclear sequence loci) and morphological characteristics (20 indices) from worldwide samples of this taxon to establish the species limits, morphological diversity and differentiation. Co-occurring spatangoid species were also analyzed with mitochondrial DNA. The nuclear sequences confirm that mitochondrial lineages correspond to true genetic entities and reveal that two clades (named A and B1) hybridize in their sympatry area, although a more closely related pair of clades (B1 and B2), whose distributions widely overlap, does not display hybridization. The morphology of all E. cordatum clade pairs was significantly differentiated, but no morphological diagnostic character was evidenced. By contrast, other spatangoid species pairs that diverged more recently than the E. cordatum clades display clear diagnostic characters. Morphological diversity thus appears responsible for the absence of diagnostic characters, ruling out stabilizing selection, a classical explanation for cryptic species. Alternative classical explanations are (i) environmental plasticity or (ii) a high diversity of genes determining morphology, maintained by varying environmental conditions. We suggest a new hypothesis that the observed morphological diversity is selectively neutral and reflects high effective population sizes in the E. cordatum complex. It is supported by the higher abundance of this taxon compared with other taxa, a trend for the genetic and morphological diversity to be correlated in Europe, and the higher genetic and morphological diversities found in clades of E cordatum (except B1) than in other spatangoid samples in Europe. However, the Pacific clades do not confirm these trends.

Understanding trophic relationships among Caribbean sea urchins.

The species Echinometra lucunter, Echinometra viridis, Lytechinus variegatus, Tripneustes ventricosus, and Diadema antillarum are the most common sea urchins of littoral habitats in the Caribbean. T. ventricosus and L. variegatus are associated with seagrass beds, while the other three species usually inhabit hardground substrates. Food preferences of these species are well documented and they are commonly accepted as being primarily herbivorous-omnivorous; nevertheless, few of them have previously been characterized isotopically. We used this approach for assessing the isotopic characterization of five echinoids. We established the trophic position of two groups of co-occurring species and quantified the contribution of food resources in the diet of Echinometra lucunter, considered the most common sea urchin in the Caribbean region. The species T. ventricosus and D. antillarum showed the highest values of δ15N. Sea urchins exhibited similar values of δ13C varying from -11.6 ± 0.63 to -10.4 ± 0.99%. The echinoid E. lucunter displayed the lowest values of carbon, from -15.40 ± 0.76%. Significant differences among species were found for δ15N and δ13C. Seaweed communities exhibited no differences among sites for overall δ15N (F= 1.300, df= 3, p= 0.301), but we found spatial differences for δ13C (F= 7.410, df= 3, p= 0.001). The ellipse-based metrics of niche width analysis found that the hardground biotope species (D. antillarum, E. lucunter, and E. viridis) did not overlap each other. Similar results were obtained for the co-occurring species of the seagrass biotope; however, the distance between these species was closer than that of the hardground biotope species. The Bayesian mixing models run for E. lucunter at all four localities found differences in food resources contribution. The algae D. menstrualis, C. crassa and B. triquetrum dominated in CGD; whereas C. nitens, Gracilaria spp., and D. caribaea represented the main contributor algae to the diet of E. lucunter at LQY. In Culebra Island, no dominance of any particular algae was detected in TMD, where six of the eight species exhibited a similar contribution. Similarities in δ15N between D. antillarum and T. ventricosus may hint towards a similar trophic level for these species, although T. ventricosus is widely accepted as an omnivore, while D. antillarum is considered a generalist herbivore. The lack of overlap among species in the two biotopes seems to indicate a resource partitioning strategy to avoid niche competition among co-occurring species.

Evolution of gamete attraction molecules: evidence for purifying selection in speract and its receptor, in the pantropical sea urchin Diadema.

Many free-spawning marine invertebrates, such as sea urchins, lack any courtship or assortative mating behavior. Mate recognition in such cases occur at the gametic level, and molecules present on the sperm and egg are major determinants of species-specific fertilization. These molecules must also coevolve in relation to each other in order to preserve functional integrity. When sea urchins release their gametes in seawater, diffusible molecules from the egg, termed sperm-activating peptides, activate and attract the sperm to swim toward the egg, initiating a series of interactions between the gametes. Although the compositions and diversity of such sperm-activating peptides have been characterized in a variety of sea urchins, little is known about the evolution of their genes. Here we characterize the genes encoding the sperm-activating peptide of the egg (speract) and its receptor on the sperm, and examine their evolutionary dynamics in the sea urchin genus Diadema, in the interest of determining whether they are involved in reproductive isolation between the species. We found evidence of purifying selection on several codon sites in both molecules and of selectively neutral evolution in others. The diffusible speract peptide that activates sperm is invariant across species, indicating that Diadema egg peptides do not discriminate between con- and hetero-specific sperm at this stage of the process. Speract and its receptor do not contribute to reproductive isolation in Diadema.

Polysaccharide Constituents of Three Types of Sea Urchin Shells and Their Anti-Inflammatory Activities.

As a source of potent anti-inflammatory traditional medicines, the quantitative chromatographic fingerprints of sea urchin shell polysaccharides were well established via pre-column derivatization high performance liquid chromatography (HPLC) analysis. Based on the quantitative results, the content of fucose and glucose could be used as preliminary distinguishing indicators among three sea urchin shell species. Besides, the anti-inflammatory activities of the polysaccharides from sea urchin shells and their gonads were also determined. The gonad polysaccharide of Anthocidaris crassispina showed the most potent anti-inflammatory activity among all samples tested.

Genomic organization of Hox and ParaHox clusters in the echinoderm, Acanthaster planci.

The organization of echinoderm Hox clusters is of interest due to the role that Hox genes play in deuterostome development and body plan organization, and the unique gene order of the Hox complex in the sea urchin Strongylocentrotus purpuratus, which has been linked to the unique development of the axial region. Here, it has been reported that the Hox and ParaHox clusters of Acanthaster planci, a corallivorous starfish found in the Pacific and Indian oceans, generally resembles the chordate and hemichordate clusters. The A. planci Hox cluster shared with sea urchins the loss of one of the medial Hox genes, even-skipped (Evx) at the anterior of the cluster, as well as organization of the posterior Hox genes.

The transcriptome of the NZ endemic sea urchin Kina (Evechinus chloroticus).

BACKGROUND: Sea urchins are studied as model organisms for developmental and systems biology and also produce highly valued food products. Evechinus chloroticus (Kina) is a sea urchin species that is indigenous to New Zealand. It is the type member of the Evechinus genus based on its morphological characteristics. Previous research has focused on identifying physical factors affecting commercial roe quality of E. chloroticus, but there is almost no genetic information available for E. chloroticus. E. chloroticus is the only species in its genus and has yet to be subject to molecular phylogenetic analysis. RESULTS: In this study we performed a de novo transcriptome assembly of Illumina sequencing data. A total of 123 million 100 base length paired-end reads were generated using RNA-Seq libraries from a range of E. chloroticus tissues from two individuals obtained from Fiordland, New Zealand. The assembly resulted in a set of 75,002 transcripts with an accepted read coverage and length, of which 24,655 transcripts could be functionally annotated using protein similarity. Transcripts could be further annotated with Gene Ontology, KEGG Orthology and InterPro terms. With this sequence data we could perform the first phylogenetic analysis of E. chloroticus to other species of its family using multiple genes. When sequences for the mitochondrial nitrogen dehydrogenase genes were compared, E. chloroticus remained outside of a family level clade, which indicated E. chloroticus is indeed a genetically distinct genus within its family. CONCLUSIONS: This study has produced a large set of E. chloroticus transcripts/proteins along with functional annotations, vastly increasing the amount of genomic data available for this species. This provides a resource for current and future studies on E. chloroticus, either to increase its commercial value, or its use as a model organism. The phylogenetic results provide a basis for further analysis of relationships between E. chloroticus, its family members, and its evolutionary history.

Discovery of urchin-associated densoviruses (family Parvoviridae) in coastal waters of the Big Island, Hawaii.

Echinoderms are important constituents of marine ecosystems, where they may influence the recruitment success of benthic flora and fauna, and are important consumers of detritus and plant materials. There are currently no described viruses of echinoderms. We used a viral metagenomic approach to examine viral consortia within three urchins - Colobocentrotus atratus, Tripneustes gratilla and Echinometra mathaei - which are common constituents of reef communities in the Hawaiian archipelago. Metagenomic libraries revealed the presence of bacteriophages and densoviruses (family Parvoviridae) in tissues of all three urchins. Densoviruses are known typically to infect terrestrial and aquatic arthropods. Urchin-associated densoviruses were detected by quantitative PCR in all tissues tested, and were also detected in filtered suspended matter (>0.2 µm) from plankton and in sediments at several locations near to where the urchins were collected for metagenomic analysis. To the best of our knowledge, this is the first report of echinoderm-associated viruses, which extends the known host range of parvoviruses.

Effects of delayed metamorphosis on larval survival, metamorphosis, and juvenile performance of four closely related species of tropical sea urchins (genus Echinometra).

We report here, the effects of extended competency on larval survival, metamorphosis, and postlarval juvenile growth of four closely related species of tropical sea urchins, Echinometra sp. A (Ea), E. mathaei (Em), Echinometra sp. C (Ec), and E. oblonga (Eo). Planktotrophic larvae of all four species fed on cultured phytoplankton (Chaetoceros gracilis) attained metamorphic competence within 22-24 days after fertilization. Competent larvae were forced to delay metamorphosis for up to 5 months by preventing them from settling in culture bottles with continuous stirring on a set of 10 rpm rotating rollers and larval survival per monthly intervals was recorded. Larval survival was highest at 24 days, when competence was attained (0 delayed period), and there were no significant differences among the four species. Larvae that had experienced a prolonged delay had reduced survival rate, metamorphosis success, and juvenile survival, but among older larvae, Em had the highest success followed by Ea, Eo, and Ec. Juveniles from larvae of all four species that metamorphosed soon after becoming competent tended to have higher growth rates (test diameter and length of spines) than juveniles from larvae that metamorphosed after a prolonged period of competence with progressively slower growth the longer the prolonged period. Despite the adverse effects of delaying metamorphosis on growth parameters, competent larvae of all four species were able to survive up to 5 months and after metamorphosis grew into 1-month-old juveniles in lab condition. Overall, delayed larvae of Em showed significantly higher larval survival, metamorphosis, and juvenile survival than Ea and Eo, while Ec showed the lowest values in these performances. Em has the most widespread distribution of these species ranging from Africa to Hawaii, while Ec probably has the most restricted distribution. Consequently, differences in distribution may be related to differences in the ability to delay metamorphosis.

Phylogenomics of strongylocentrotid sea urchins.

BACKGROUND: Strongylocentrotid sea urchins have a long tradition as model organisms for studying many fundamental processes in biology including fertilization, embryology, development and genome regulation but the phylogenetic relationships of the group remain largely unresolved. Although the differing isolating mechanisms of vicariance and rapidly evolving gamete recognition proteins have been proposed, a stable and robust phylogeny is unavailable. RESULTS: We used a phylogenomic approach with mitochondrial and nuclear genes taking advantage of the whole-genome sequencing of nine species in the group to establish a stable (i.e. concordance in tree topology among multiple lies of evidence) and robust (i.e. high nodal support) phylogenetic hypothesis for the family Strongylocentrotidae. We generated eight draft mitochondrial genome assemblies and obtained 13 complete mitochondrial genes for each species. Consistent with previous studies, mitochondrial sequences failed to provide a reliable phylogeny. In contrast, we obtained a very well-supported phylogeny from 2301 nuclear genes without evidence of positive Darwinian selection both from the majority of most-likely gene trees and the concatenated fourfold degenerate sites: ((P. depressus, (M. nudus, M. franciscanus), (H. pulcherrimus, (S. purpuratus, (S. fragilis, (S. pallidus, (S. droebachiensis, S. intermedius)). This phylogeny was consistent with a single invasion of deep-water environments followed by a holarctic expansion by Strongylocentrotus. Divergence times for each species estimated with reference to the divergence times between the two major clades of the group suggest a correspondence in the timing with the opening of the Bering Strait and the invasion of the holarctic regions. CONCLUSIONS: Nuclear genome data contains phylogenetic signal informative for understanding the evolutionary history of this group. However, mitochondrial genome data does not. Vicariance can explain major patterns observed in the phylogeny. Other isolating mechanisms are appropriate to explore in this system to help explain divergence patterns not well supported by vicariance, such as the effects of rapidly evolving gamete recognition proteins on isolating populations. Our findings of a stable and robust phylogeny, with the increase in mitochondrial and nuclear comparative genomic data, provide a system in which we can enhance our understanding of molecular evolution and adaptation in this group of sea urchins.