Low coverage sequencing of three echinoderm genomes: the brittle star Ophionereis fasciata, the sea star Patiriella regularis, and the sea cucumber Australostichopus mollis

Genomic and 16S metabarcoding data of Holothuria tubulosa Gmelin, 1791

Holothuria tubulosa Gmelin, 1791 is an edible sea cucumber species widespread in the Mediterranean Sea with ecological and increasing economic importance. Genome data of holothurian species is limited and the availability of genomic data resources is crucial in understanding their biology and adaptability mechanisms. This dataset presents the raw genome sequence data of H. tubulosa sequenced on an Illumina NextSeq 2000 platform. Genome size estimation was performed based on k-mer frequency approach. Additionally, the bacterial microbiome in the stomach and intestine of H. tubulosa collected from the Strymonian Gulf (North Aegean Sea, Greece) through 16S rRNA amplicon metabarcoding sequencing is reported. Sequencing was performed on an Illumina MiSeq platform. Analysis was conducted using the QIIME2 software package, the DADA2 algorithm and a trained classifier for taxonomy assignment. The datasets presented in this work serve as valuable resources for a comprehensive investigation of H. tubulosa at the genome level and for comparative genomics and echinoderms gut microbial studies.

The Holothuria leucospilota genome elucidates sacrificial organ expulsion and bioadhesive trap enriched with amyloid-patterned proteins

Some tropical sea cucumbers of the family Holothuriidae can efficiently repel or even fatally ensnare predators by sacrificially ejecting a bioadhesive matrix termed the Cuvierian organ (CO), so named by the French zoologist Georges Cuvier who first described it in 1831. Still, the precise mechanisms for how adhesiveness genetically arose in CO and how sea cucumbers perceive and transduce danger signals for CO expulsion during defense have remained unclear. Here, we report the first high-quality, chromosome-level genome assembly of Holothuria leucospilota, an ecologically significant sea cucumber with prototypical CO. The H. leucospilota genome reveals characteristic long-repeat signatures in CO-specific outer-layer proteins, analogous to fibrous proteins of disparate species origins, including spider spidroin and silkworm fibroin. Intriguingly, several CO-specific proteins occur with amyloid-like patterns featuring extensive intramolecular cross-β structures readily stainable by amyloid indicator dyes. Distinct proteins within the CO connective tissue and outer surface cooperate to give the expelled matrix its apparent tenacity and adhesiveness, respectively. Genomic evidence offers further hints that H. leucospilota directly transduces predator-induced mechanical pressure onto the CO surface through mediation by transient receptor potential channels, which culminates in acetylcholine-triggered CO expulsion in part or in entirety. Evolutionarily, innovative events in two distinct regions of the H. leucospilota genome have apparently spurred CO's differentiation from the respiratory tree to a lethal defensive organ against predators.

Twinkle twinkle brittle star: the draft genome of Ophioderma brevispinum (Echinodermata: Ophiuroidea) as a resource for regeneration research

Background

Echinoderms are established models in experimental and developmental biology, however genomic resources are still lacking for many species. Here, we present the draft genome of Ophioderma brevispinum, an emerging model organism in the field of regenerative biology. This new genomic resource provides a reference for experimental studies of regenerative mechanisms.

Results

We report a de novo nuclear genome assembly for the brittle star O. brevispinum and annotation facilitated by the transcriptome assembly. The final assembly is 2.68 Gb in length and contains 146,703 predicted protein-coding gene models. We also report a mitochondrial genome for this species, which is 15,831 bp in length, and contains 13 protein-coding, 22 tRNAs, and 2 rRNAs genes, respectively. In addition, 29 genes of the Notch signaling pathway are identified to illustrate the practical utility of the assembly for studies of regeneration.

Conclusions

The sequenced and annotated genome of O. brevispinum presented here provides the first such resource for an ophiuroid model species. Considering the remarkable regenerative capacity of this species, this genome will be an essential resource in future research efforts on molecular mechanisms regulating regeneration.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-022-08750-y).

Architecture and evolution of the cis-regulatory system of the echinoderm kirrelL gene

The gene regulatory network (GRN) that underlies echinoderm skeletogenesis is a prominent model of GRN architecture and evolution. KirrelL is an essential downstream effector gene in this network and encodes an Ig-superfamily protein required for the fusion of skeletogenic cells and the formation of the skeleton. In this study, we dissected the transcriptional control region of the kirrelL gene of the purple sea urchin, Strongylocentrotus purpuratus. Using plasmid- and bacterial artificial chromosome-based transgenic reporter assays, we identified key cis-regulatory elements (CREs) and transcription factor inputs that regulate Sp-kirrelL, including direct, positive inputs from two key transcription factors in the skeletogenic GRN, Alx1 and Ets1. We next identified kirrelL cis-regulatory regions from seven other echinoderm species that together represent all classes within the phylum. By introducing these heterologous regulatory regions into developing sea urchin embryos we provide evidence of their remarkable conservation across ~500 million years of evolution. We dissected in detail the kirrelL regulatory region of the sea star, Patiria miniata, and demonstrated that it also receives direct inputs from Alx1 and Ets1. Our findings identify kirrelL as a component of the ancestral echinoderm skeletogenic GRN. They support the view that GRN subcircuits, including specific transcription factor–CRE interactions, can remain stable over vast periods of evolutionary history. Lastly, our analysis of kirrelL establishes direct linkages between a developmental GRN and an effector gene that controls a key morphogenetic cell behavior, cell–cell fusion, providing a paradigm for extending the explanatory power of GRNs.

Facilitating population genomics of non-model organisms through optimized experimental design for reduced representation sequencing

BACKGROUND

Genome-wide data are invaluable to characterize differentiation and adaptation of natural populations. Reduced representation sequencing (RRS) subsamples a genome repeatedly across many individuals. However, RRS requires careful optimization and fine-tuning to deliver high marker density while being cost-efficient. The number of genomic fragments created through restriction enzyme digestion and the sequencing library setup must match to achieve sufficient sequencing coverage per locus. Here, we present a workflow based on published information and computational and experimental procedures to investigate and streamline the applicability of RRS.

RESULTS

In an iterative process genome size estimates, restriction enzymes and size selection windows were tested and scaled in six classes of Antarctic animals (Ostracoda, Malacostraca, Bivalvia, Asteroidea, Actinopterygii, Aves). Achieving high marker density would be expensive in amphipods, the malacostracan target taxon, due to the large genome size. We propose alternative approaches such as mitogenome or target capture sequencing for this group. Pilot libraries were sequenced for all other target taxa. Ostracods, bivalves, sea stars, and fish showed overall good coverage and marker numbers for downstream population genomic analyses. In contrast, the bird test library produced low coverage and few polymorphic loci, likely due to degraded DNA.

CONCLUSIONS

Prior testing and optimization are important to identify which groups are amenable for RRS and where alternative methods may currently offer better cost-benefit ratios. The steps outlined here are easy to follow for other non-model taxa with little genomic resources, thus stimulating efficient resource use for the many pressing research questions in molecular ecology.

Draft Genome of the Sea Cucumber Holothuria glaberrima, a Model for the Study of Regeneration

Regeneration is one of the most fascinating and yet least understood biological processes. Echinoderms, one of the closest related invertebrate groups to humans, can contribute to our understanding of the genetic basis of regenerative processes. Among echinoderms, sea cucumbers have the ability to grow back most of their body parts following injury, including the intestine and nervous tissue. The cellular and molecular events underlying these abilities in sea cucumbers have been most extensively studied in the species Holothuria glaberrima. However, research into the regenerative abilities of this species has been impeded due to the lack of adequate genomic resources. Here, we report the first draft genome assembly of H. glaberrima and demonstrate its value for future genetic studies. Using only short sequencing reads, we assembled the genome into 89,105 scaffolds totaling 1.1 gigabases with an N50 of 25 kilobases. Our BUSCO assessment of the genome resulted in 894 (91.4%) complete and partial genes from 978 genes queried. We incorporated transcriptomic data from several different life history stages to annotate 51,415 genes in our final assembly. To demonstrate the usefulness of the genome, we fully annotated the melanotransferrin (Mtf) gene family, which have a potential role in the regeneration of the sea cucumber intestine. Using these same data, we extracted the mitochondrial genome, which showed high conservation to that of other holothuroids. Thus, these data will be a critical resource for ongoing studies of regeneration and other studies in sea cucumbers.

BLAST-QC: automated analysis of BLAST results

BACKGROUND

The Basic Local Alignment Search Tool (BLAST) from NCBI is the preferred utility for sequence alignment and identification for bioinformatics and genomics research. Among researchers using NCBI's BLAST software, it is well known that analyzing the results of a large BLAST search can be tedious and time-consuming. Furthermore, with the recent discussions over the effects of parameters such as '-max_target_seqs' on the BLAST heuristic search process, the use of these search options are questionable. This leaves using a stand-alone parser as one of the only options of condensing these large datasets, and with few available for download online, the task is left to the researcher to create a specialized piece of software anytime they need to analyze BLAST results. The need for a streamlined and fast script that solves these issues and can be easily implemented into a variety of bioinformatics and genomics workflows was the initial motivation for developing this software.

RESULTS

In this study, we demonstrate the effectiveness of BLAST-QC for analysis of BLAST results and its desirability over the other available options. Applying genetic sequence data from our bioinformatic workflows, we establish BLAST_QC's superior runtime when compared to existing parsers developed with commonly used BioPerl and BioPython modules, as well as C and Java implementations of the BLAST_QC program. We discuss the 'max_target_seqs' parameter, the usage of and controversy around the use of the parameter, and offer a solution by demonstrating the ability of our software to provide the functionality this parameter was assumed to produce, as well as a variety of other parsing options. Executions of the script on example datasets are given, demonstrating the implemented functionality and providing test-cases of the program. BLAST-QC is designed to be integrated into existing software, and we establish its effectiveness as a module of workflows or other processes.

CONCLUSIONS

BLAST-QC provides the community with a simple, lightweight and portable Python script that allows for easy quality control of BLAST results while avoiding the drawbacks of other options. This includes the uncertain results of applying the -max_target_seqs parameter or relying on the cumbersome dependencies of other options like BioPerl, Java, etc. which add complexity and run time when running large data sets of sequences. BLAST-QC is ideal for use in high-throughput workflows and pipelines common in bioinformatic and genomic research, and the script has been designed for portability and easy integration into whatever type of processes the user may be running.

Identification of molecular markers for superior quantitative traits in a novel sea cucumber strain by comparative microRNA-mRNA expression profiling

To investigate the adaptability of Apostichopus japonicus (A. japonicus) strain "Anyuan No. 1" in the South China Sea, field monitoring and microRNA-mRNA integrated analyses were conducted between "Anyuan No. 1" and a regular A. japonicus population from Wendeng (Shandong Province, as a control) in the Xiapu farming area in Fujian Province, China. The results showed that "Anyuan No. 1" exhibited greater body weight increase and a higher number of papillae compared to the control during two and a half months of field monitoring. Comparative microRNA (miRNA) and mRNA transcriptome analyses identified 12 differentially expressed miRNAs (DEMs) and 165 differentially expressed genes (DEGs) in "Anyuan No. 1" compared to the control. Long-chain specific acyl-CoA dehydrogenase (ACADL), transmembrane protein 251 (TMEM251), dehydrogenase/reductase SDR family protein 7-like (Dhrs7), insulin-like growth factor-binding protein 7 (IGFBP-7), CDK5 regulatory subunit-associated protein 1 (CDK5RAP1), visual pigment-like receptor peropsin, 39S ribosomal protein, miR-10, miR-153, miR-7, and miR-3529 were identified as gene and miRNA candidates correlated with superior economic traits in "Anyuan No. 1". Collectively, "Anyuan No. 1" is suitable for large-scale cultivation extension due to its better adaptability to the South China Sea area. Furthermore, we identified "miR10-ACADL" as a potential module for further molecular marker-assisted selective breeding of A. japonicus.

Microscopic ossicle analyses and the complete mitochondrial genome sequence of Holothuria (Roweothuria) polii (Echinodermata; Holothuroidea) provide new information to support the phylogenetic positioning of this sea cucumber species

Sea cucumbers (Holothuroidea) are ecologically important organisms for their bioturbation and alkalinization activities of the seabed. These species are extensively fished as they are considered luxury food. Sea cucumbers are also relevant for biomedical studies and the production of bioactive compounds. A few initiatives are recently evaluating sea cucumbers as novel aquaculture species. The aim of this study was to provide morphological and genetic information useful for the identification of Holothuria polii, the white spot sea cucumber (a common species of the Mediterranean Sea). We generated the complete sequence of the mitochondrial DNA (mtDNA) genome of this species and combined it with a detailed ossicle characterization of the sequenced specimen by scanning electron microscopic analysis. Ossicles (known also as sclerites) are anatomical features that can discriminate Holothuroidea species, including the closely related ones of the genus Holothuria. The complete mitochondrial genome was assembled, functionally annotated and then used to evaluate the phylogenetic relationship of H. polii against the other few Holothuroidea species for which the whole mtDNA was available. The 15,907 bp H. polii mtDNA sequence has the same gene order already reported for H. scabra, H. forskali and other species of the same class. Cox1 and 16S gene sequences were informative for species identification across the genus and could be used for the authentication of commercialized Holothuria spp. The mitochondrial genome sequence presented here provides the basis to a future analysis of the variability of H. polii populations in the Mediterranean region.

Sea cucumber genome provides insights into saponin biosynthesis and aestivation regulation

Echinoderms exhibit several fascinating evolutionary innovations that are rarely seen in the animal kingdom, but how these animals attained such features is not well understood. Here we report the sequencing and analysis of the genome and extensive transcriptomes of the sea cucumber Apostichopus japonicus, a species from a special echinoderm group with extraordinary potential for saponin synthesis, aestivation and organ regeneration. The sea cucumber does not possess a reorganized Hox cluster as previously assumed for all echinoderms, and the spatial expression of Hox7 and Hox11/13b potentially guides the embryo-to-larva axial transformation. Contrary to the typical production of lanosterol in animal cholesterol synthesis, the oxidosqualene cyclase of sea cucumber produces parkeol for saponin synthesis and has "plant-like" motifs suggestive of convergent evolution. The transcriptional factors Klf2 and Egr1 are identified as key regulators of aestivation, probably exerting their effects through a clock gene-controlled process. Intestinal hypometabolism during aestivation is driven by the DNA hypermethylation of various metabolic gene pathways, whereas the transcriptional network of intestine regeneration involves diverse signaling pathways, including Wnt, Hippo and FGF. Decoding the sea cucumber genome provides a new avenue for an in-depth understanding of the extraordinary features of sea cucumbers and other echinoderms.

Ancient Origin of the CARD-Coiled Coil/Bcl10/MALT1-Like Paracaspase Signaling Complex Indicates Unknown Critical Functions

The CARD–coiled coil (CC)/Bcl10/MALT1-like paracaspase (CBM) signaling complexes composed of a CARD–CC family member (CARD-9, -10, -11, or -14), Bcl10, and the type 1 paracaspase MALT1 (PCASP1) play a pivotal role in immunity, inflammation, and cancer. Targeting MALT1 proteolytic activity is of potential therapeutic interest. However, little is known about the evolutionary origin and the original functions of the CBM complex. Type 1 paracaspases originated before the last common ancestor of planulozoa (bilaterians and cnidarians). Notably in bilaterians, Ecdysozoa (e.g., nematodes and insects) lacks Bcl10, whereas other lineages have a Bcl10 homolog. A survey of invertebrate CARD–CC homologs revealed such homologs only in species with Bcl10, indicating an ancient common origin of the entire CBM complex. Furthermore, vertebrate-like Syk/Zap70 tyrosine kinase homologs with the ITAM-binding SH2 domain were only found in invertebrate organisms with CARD–CC/Bcl10, indicating that this pathway might be related to the original function of the CBM complex. Moreover, the type 1 paracaspase sequences from invertebrate organisms that have CARD–CC/Bcl10 are more similar to vertebrate paracaspases. Functional analysis of protein–protein interactions, NF-κB signaling, and CYLD cleavage for selected invertebrate type 1 paracaspase and Bcl10 homologs supports this scenario and indicates an ancient origin of the CARD–CC/Bcl10/paracaspase signaling complex. By contrast, many of the known MALT1-associated activities evolved fairly recently, indicating that unknown functions are at the basis of the protein conservation. As a proof-of-concept, we provide initial evidence for a CBM- and NF-κB-independent neuronal function of the Caenorhabditis elegans type 1 paracaspase malt-1. In conclusion, this study shows how evolutionary insights may point at alternative functions of MALT1.

The complex simplicity of the brittle star nervous system

Background

Brittle stars (Ophiuroidea, Echinodermata) have been increasingly used in studies of animal behavior, locomotion, regeneration, physiology, and bioluminescence. The success of these studies directly depends on good working knowledge of the ophiuroid nervous system.

Results

Here, we describe the arm nervous system at different levels of organization, including the microanatomy of the radial nerve cord and peripheral nerves, ultrastructure of the neural tissue, and localization of different cell types using specific antibody markers. We standardize the nomenclature of nerves and ganglia, and provide an anatomically accurate digital 3D model of the arm nervous system as a reference for future studies. Our results helped identify several general features characteristic to the adult echinoderm nervous system, including the extensive anatomical interconnections between the ectoneural and hyponeural components, neuroepithelial organization of the central nervous system, and the supporting scaffold of the neuroepithelium formed by radial glial cells. In addition, we provide further support to the notion that the echinoderm radial glia is a complex and diverse cell population. We also tested the suitability of a range of specific cell-type markers for studies of the brittle star nervous system and established that the radial glial cells are reliably labeled with the ERG1 antibodies, whereas the best neuronal markers are acetylated tubulin, ELAV, and synaptotagmin B. The transcription factor Brn1/2/4 – a marker of neuronal progenitors – is expressed not only in neurons, but also in a subpopulation of radial glia. For the first time, we describe putative ophiuroid proprioceptors associated with the hyponeural part of the central nervous system.

Conclusions

Together, our data help establish both the general principles of neural architecture common to the phylum Echinodermata and the specific ophiuroid features.

Electronic supplementary material

The online version of this article (10.1186/s12983-017-0247-4) contains supplementary material, which is available to authorized users.

Holothurians as a Model System to Study Regeneration

Echinoderms possess an incredible regenerative capacity. Within this phylum, holothurians, better known as sea cucumbers, can regenerate most of their internal and external organs. While regeneration has been studied in several species, the most recent and extensive studies have been done in the species Holothuria glaberrima, the focus of most of our discussion. This chapter presents the model system and integrates the work that has been done to determine the major steps that take place, during regeneration of the intestinal and nervous system, from wound healing to the reestablishment of original function. We describe the cellular and molecular events associated with the regeneration processes and also describe the techniques that have been used, discuss the results, and explain the gaps in our knowledge that remain. We expect that the information provided here paves the road for new and young investigators to continue the study of the amazing potential of regeneration in members of the Echinodermata and how these studies will shed some light into the mechanisms that are common to many regenerative processes.

An Organismal Model for Gene Regulatory Networks in the Gut-Associated Immune Response

The gut epithelium is an ancient site of complex communication between the animal immune system and the microbial world. While elements of self-non-self receptors and effector mechanisms differ greatly among animal phyla, some aspects of recognition, regulation, and response are broadly conserved. A gene regulatory network (GRN) approach provides a means to investigate the nature of this conservation and divergence even as more peripheral functional details remain incompletely understood. The sea urchin embryo is an unparalleled experimental model for detangling the GRNs that govern embryonic development. By applying this theoretical framework to the free swimming, feeding larval stage of the purple sea urchin, it is possible to delineate the conserved regulatory circuitry that regulates the gut-associated immune response. This model provides a morphologically simple system in which to efficiently unravel regulatory connections that are phylogenetically relevant to immunity in vertebrates. Here, we review the organism-wide cellular and transcriptional immune response of the sea urchin larva. A large set of transcription factors and signal systems, including epithelial expression of interleukin 17 (IL17), are important mediators in the activation of the early gut-associated response. Many of these have homologs that are active in vertebrate immunity, while others are ancient in animals but absent in vertebrates or specific to echinoderms. This larval model provides a means to experimentally characterize immune function encoded in the sea urchin genome and the regulatory interconnections that control immune response and resolution across the tissues of the organism.

High-throughput targeted screening in triple-negative breast cancer cells identifies Wnt-inhibiting activities in Pacific brittle stars

Pro-proliferative oncogenic signaling is one of the hallmarks of cancer. Specific targeting of such signaling pathways is one of the main approaches to modern anti-cancer drug discovery, as opposed to more traditional search for general cytotoxic agents. Natural products, especially from marine sources, represent a largely untapped source of chemical diversity, which so far have mostly been screened for cytotoxicity. Here we present a pioneering pipeline of high-throughput screening of marine-based activities targeted against the Wnt signaling pathway, which is one of the key factors in oncogenic transformation, growth and metastasis in different cancers, including the devastating triple-negative breast cancer (TNBC) currently lacking any targeted therapies. This pipeline consisted of collection and characterization of numerous invertebrates during the SokhoBio expedition to the Kuril Basin in North Pacific, preparation of extracts from these specimen, and their screening in dedicated assays monitoring Wnt signaling in TNBC cells. This approach yielded a number of promising hits, including highly specific anti-Wnt activities targeting multiple levels within the Wnt pathway from Ophiura irrorata and other Pacific brittle stars.

Draft genome of the sea cucumber Apostichopus japonicus and genetic polymorphism among color variants

The Japanese sea cucumber (Apostichopus japonicus Selenka 1867) is an economically important species as a source of seafood and ingredient in traditional medicine. It is mainly found off the coasts of northeast Asia. Recently, substantial exploitation and widespread biotic diseases in A. japonicus have generated increasing conservation concern. However, the genomic knowledge base and resources available for researchers to use in managing this natural resource and to establish genetically based breeding systems for sea cucumber aquaculture are still in a nascent stage. A total of 312 Gb of raw sequences were generated using the Illumina HiSeq 2000 platform and assembled to a final size of 0.66 Gb, which is about 80.5% of the estimated genome size (0.82 Gb). We observed nucleotide-level heterozygosity within the assembled genome to be 0.986%. The resulting draft genome assembly comprising 132 607 scaffolds with an N50 value of 10.5 kb contains a total of 21 771 predicted protein-coding genes. We identified 6.6-14.5 million heterozygous single nucleotide polymorphisms in the assembled genome of the three natural color variants (green, red, and black), resulting in an estimated nucleotide diversity of 0.00146. We report the first draft genome of A. japonicus and provide a general overview of the genetic variation in the three major color variants of A. japonicus. These data will help provide a comprehensive view of the genetic, physiological, and evolutionary relationships among color variants in A. japonicus, and will be invaluable resources for sea cucumber genomic research.

Echinoderm development and evolution in the post-genomic era

The highly recognizable animals within the phylum Echinodermata encompass an enormous disparity of adult and larval body plans. The extensive knowledge of sea urchin development has culminated in the description of the exquisitely detailed gene regulatory network (GRN) that governs the specification of various embryonic territories. This information provides a unique opportunity for comparative studies in other echinoderm taxa to understand the evolution and developmental mechanisms underlying body plan change. This review focuses on recent work that has utilized new genomic resources and systems-level experiments to address questions of evolutionary developmental biology. In particular, we synthesize the results of several recent studies from various echinoderm classes that have explored the development and evolution of the larval skeleton, which is a major feature that distinguishes the two predominant larval subtypes within the Phylum. We specifically examine the ways in which GRNs can evolve, either through cis regulatory and/or protein-level changes in transcription factors. We also examine recent work comparing evolution across shorter time scales that occur within and between species of sea urchin, and highlight the kinds of questions that can be addressed by these comparisons. The advent of new genomic and transcriptomic datasets in additional species from all classes of echinoderm will continue to empower the use of these taxa for evolutionary developmental studies.

Low coverage sequencing of three echinoderm genomes: the brittle star Ophionereis fasciata, the sea star Patiriella regularis, and the sea cucumber Australostichopus mollis

BACKGROUND

There are five major extant groups of Echinodermata: Crinoidea (feather stars and sea lillies), Ophiuroidea (brittle stars and basket stars), Asteroidea (sea stars), Echinoidea (sea urchins, sea biscuits, and sand dollars), and Holothuroidea (sea cucumbers). These animals are known for their pentaradial symmetry as adults, unique water vascular system, mutable collagenous tissues, and endoskeletons of high magnesium calcite. To our knowledge, the only echinoderm species with a genome sequence available to date is Strongylocentrotus pupuratus (Echinoidea). The availability of additional echinoderm genome sequences is crucial for understanding the biology of these animals.

FINDINGS

Here we present assembled draft genomes of the brittle star Ophionereis fasciata, the sea star Patiriella regularis, and the sea cucumber Australostichopus mollis from Illumina sequence data with coverages of 12.5x, 22.5x, and 21.4x, respectively.

CONCLUSIONS

These data provide a resource for mining gene superfamilies, identifying non-coding RNAs, confirming gene losses, and designing experimental constructs. They will be important comparative resources for future genomic studies in echinoderms.