SpBase: the sea urchin genome database and web site

Echinobase: a resource to support the echinoderm research community

Echinobase (www.echinobase.org) is a model organism knowledgebase serving as a resource for the community that studies echinoderms, a phylum of marine invertebrates that includes sea urchins and sea stars. Echinoderms have been important experimental models for over 100 years and continue to make important contributions to environmental, evolutionary, and developmental studies, including research on developmental gene regulatory networks. As a centralized resource, Echinobase hosts genomes and collects functional genomic data, reagents, literature, and other information for the community. This third-generation site is based on the Xenbase knowledgebase design and utilizes gene-centric pages to minimize the time and effort required to access genomic information. Summary gene pages display gene symbols and names, functional data, links to the JBrowse genome browser, and orthology to other organisms and reagents, and tabs from the Summary gene page contain more detailed information concerning mRNAs, proteins, diseases, and protein-protein interactions. The gene pages also display 1:1 orthologs between the fully supported species Strongylocentrotus purpuratus (purple sea urchin), Lytechinus variegatus (green sea urchin), Patiria miniata (bat star), and Acanthaster planci (crown-of-thorns sea star). JBrowse tracks are available for visualization of functional genomic data from both fully supported species and the partially supported species Anneissia japonica (feather star), Asterias rubens (sugar star), and L. pictus (painted sea urchin). Echinobase serves a vital role by providing researchers with annotated genomes including orthology, functional genomic data aligned to the genomes, and curated reagents and data. The Echinoderm Anatomical Ontology provides a framework for standardizing developmental data across the phylum, and knowledgebase content is formatted to be findable, accessible, interoperable, and reusable by the research community.

Gene regulatory divergence amongst echinoderms underlies appearance of pigment cells in sea urchin development

Larvae of the sea urchin, Strongylocentrotus purpuratus, have pigmented migratory cells implicated in immune defense and gut patterning. The transcription factor SpGcm activates the expression of many pigment cell-specific genes, including those involved in pigment biosynthesis (SpPks1 and SpFmo3) and immune related genes (e.g. SpMif5). Despite the importance of this cell type in sea urchins, pigmented cells are absent in larvae of the sea star, Patiria miniata. In this study, we tested the premises that sea stars lack genes to synthesize echinochrome pigment, that the genes are present but are not expressed in the larvae, or rather that the homologous gene expression does not contribute to echinochrome synthesis. Our results show that orthologs of sea urchin pigment cell-specific genes (PmPks1, PmFmo3-1 and PmMifL1-2) are present in the sea star genome and expressed in the larvae. Although no cell lineage homologous to migratory sea urchin pigment cells is present, dynamic gene activation accomplishes a similar spatial and temporal expression profile. The mechanisms regulating the expression of these genes, though, is highly divergent. In sea stars, PmGcm lacks the central role in pigment gene expression since it is not expressed in PmPks1 and PmFmo3-1-positive cells, and knockdown of Gcm does not abrogate pigment gene expression. Pigment genes are instead expressed in the coelomic mesoderm early in development before later being expressed in the ectoderm. These findings were supported by in situ RNA hybridization and comparative scRNA-seq analyses. We conclude that simply the coexpression of Pks1 and Fmo3 orthologs in cells of the sea star is not sufficient to underlie the emergence of the larval pigment cell in the sea urchin.

More than a simple epithelial layer: multifunctional role of echinoderm coelomic epithelium

In echinoderms, the coelomic epithelium (CE) is reportedly the source of new circulating cells (coelomocytes) as well as the provider of molecular factors such as immunity-related molecules. However, its overall functions have been scarcely studied in detail. In this work, we used an integrated approach based on both microscopy (light and electron) and proteomic analyses to investigate the arm CE in the starfish Marthasterias glacialis during different physiological conditions (i.e., non-regenerating and/or regenerating). Our results show that CE cells share both ultrastructural and proteomic features with circulating coelomocytes (echinoderm immune cells). Additionally, microscopy and proteomic analyses indicate that CE cells are actively involved in protein synthesis and processing, and membrane trafficking processes such as phagocytosis (particularly of myocytes) and massive secretion phenomena. The latter might provide molecules (e.g., immune factors) and fluids for proper arm growth/regrowth. No stem cell marker was identified and no pre-existing stem cell was observed within the CE. Rather, during regeneration, CE cells undergo dedifferentiation and epithelial-mesenchymal transition to deliver progenitor cells for tissue replacement. Overall, our work underlines that echinoderm CE is not a "simple epithelial lining" and that instead it plays multiple functions which span from immunity-related roles as well as being a source of regeneration-competent cells for arm growth/regrowth.

EchinoDB: an update to the web-based application for genomic and transcriptomic data on echinoderms

BACKGROUND

Here we release a new version of EchinoDB, EchinoDB v2.0 ( https://echinodb.uncc.edu ). EchinoDB is a database of genomic and transcriptomic data on echinoderms. The initial database consisted of groups of 749,397 orthologous and paralogous transcripts arranged in orthoclusters by sequence similarity.

RESULTS

The updated version of EchinoDB includes two new major datasets: the RNA-Seq data of the brittle star Ophioderma brevispinum and the high-quality genomic assembly data of the green sea urchin Lytechinus variegatus. In addition, we enabled keyword searches for annotated data and installed an updated version of Sequenceserver to allow Basic Local Alignment Search Tool (BLAST) searches. The data are downloadable in FASTA format. The first version of EchinoDB appeared in 2016 and was implemented in GO on a local server. The new version has been updated using R Shiny to include new features and improvements in the application. Furthermore, EchinoDB now runs entirely in the cloud for increased reliability and scaling.

CONCLUSION

EchinoDB serves a user base drawn from the fields of phylogenetics, developmental biology, genomics, physiology, neurobiology, and regeneration. As use cases, we illustrate the function of EchinoDB in retrieving components of signaling pathways involved in the tissue regeneration process of different echinoderms, including the emerging model species Ophioderma brevispinum. Moreover, we use EchinoDB to shed light on the conservation of the molecular components involved in two echinoderm-specific phenomena: spicule matrix proteins involved in the formation of stereom endoskeleton and the tensilin protein that contributes to the capacity of the connective tissues to quickly change its mechanical properties. The genes involved in the former had been previously studied in echinoids, while gene sequences involved in the latter had been previously described in holothuroids. Specifically, we ask (a) if the biomineralization-related proteins previously reported only in sea urchins are also present in other, non-echinoid, echinoderms and (b) if tensilin, the protein responsible for the control of stiffness of the mutable collagenous tissue, previously described in sea cucumbers, is conserved across the phylum.

A Chromosome-level Genome Assembly of the Highly Heterozygous Sea Urchin Echinometra sp. EZ Reveals Adaptation in the Regulatory Regions of Stress Response Genes

Echinometra is the most widespread genus of sea urchin and has been the focus of a wide range of studies in ecology, speciation, and reproduction. However, available genetic data for this genus are generally limited to a few select loci. Here, we present a chromosome-level genome assembly based on 10x Genomics, PacBio, and Hi-C sequencing for Echinometra sp. EZ from the Persian/Arabian Gulf. The genome is assembled into 210 scaffolds totaling 817.8 Mb with an N50 of 39.5 Mb. From this assembly, we determined that the E. sp. EZ genome consists of 2n = 42 chromosomes. BUSCO analysis showed that 95.3% of BUSCO genes were complete. Ab initio and transcript-informed gene modeling and annotation identified 29,405 genes, including a conserved Hox cluster. E. sp. EZ can be found in high-temperature and high-salinity environments, and we therefore compared E. sp. EZ gene families and transcription factors associated with environmental stress response ("defensome") with other echinoid species with similar high-quality genomic resources. While the number of defensome genes was broadly similar for all species, we identified strong signatures of positive selection in E. sp. EZ noncoding elements near genes involved in environmental response pathways as well as losses of transcription factors important for environmental response. These data provide key insights into the biology of E. sp. EZ as well as the diversification of Echinometra more widely and will serve as a useful tool for the community to explore questions in this taxonomic group and beyond.

Integrating experience with databases, bioinformatics, and wet lab exercises for students in an introductory genetics course

Students often have no exposure to the incredible amount of genomic and proteomic data that is freely and easily accessible online. Becoming familiar with these resources, and seeing how they could be applied to a specific research question, is a prerequisite for students to apply them to their own scientific development. Many students have to "see it" and "do it" before they "get it." This paper describes a teaching laboratory for undergraduate Genetics students that combines exploration of several publicly available databases with some simple bioinformatic exercises and "'real" live experience in a wet lab exercise. The goal is to teach skills in applying genomic data to a real scientific question. In this exercise, students identify a target protein after exploring several protein and signal transduction databases, such as the Kyoto Encyclopedia of Genes and Genomes database. They then search for the encoding RNA in a newly available sea star mature egg transcriptome database and for the DNA in an existing sea star genome database. The students design primers against specific regions or domains in their target RNA and amplify these by reverse transcription PCR against RNA purified from fresh sea star eggs. The PCR reactions are analyzed by agarose gel electrophoresis. It is hoped that the combination of the computational biology exercises with the real lab work will excite the students and stimulate them to explore this exciting new biology further.

The arm of the starfish: The far-reaching applications of Patiria miniata as a model system in evolutionary, developmental, and regenerative biology

Many species of echinoderms have long been considered model research organisms in biology. Historically, much of this research has focused on the embryology of sea urchins and the use of their extensive gene regulatory networks as a tool to understand how the genome controls cell state specification and patterning. The establishment of Patiria miniata, the bat sea star, as a research organism has allowed us to expand on the concepts explored with sea urchins, viewing these genetic networks through a comparative lens, gaining great insight into the evolutionary mechanisms that shape developmental diversity. Extensive molecular tools have been developed in P. miniata, designed to explore gene expression dynamics and build gene regulatory networks. Echinoderms also have a robust set of bioinformatic and computational resources, centered around echinobase.org, an extensive database containing multiomic, developmental, and experimental resources for researchers. In addition to comparative evolutionary development, P. miniata is a promising system in its own right for studying whole body regeneration, metamorphosis and body plan development, as well as marine disease.

A conserved node in the regulation of Vasa between an induced and an inherited program of primordial germ cell specification

Primordial germ cells (PGCs) are specified by diverse mechanisms in early development. In some animals, PGCs are specified via inheritance of maternal determinants, while in others, in a process thought to represent the ancestral mode, PGC fate is induced by cell interactions. Although the terminal factors expressed in specified germ cells are widely conserved, the mechanisms by which these factors are regulated can be widely diverse. Here we show that a post-translational mechanism of germ cell specification is conserved between two echinoderm species thought to employ divergent germ line segregation strategies. Sea urchins segregate their germ line early by an inherited mechanism. The DEAD-box RNA - helicase Vasa, a conserved germline factor, becomes enriched in the PGCs by degradation in future somatic cells by the E3-ubiquitin-ligase Gustavus (Gustafson et al., 2011). This post-translational activity occurs early in development, substantially prior to gastrulation. Here we test this process in germ cell specification of sea star embryos, which use inductive signaling mechanisms after gastrulation for PGC fate determination. We find that Vasa-GFP protein becomes restricted to the PGCs in the sea star even though the injected mRNA is present throughout the embryo. Gustavus depletion, however, results in uniform accumulation of the protein. These data demonstrate that Gustavus-mediated Vasa turnover in somatic cells is conserved between species with otherwise divergent PGC specification mechanisms. Since Gustavus was originally identified in Drosophila melanogaster to have similar functions in Vasa regulation (Kugler et al., 2010), we conclude that this node of Vasa regulation in PGC formation is ancestral and evolutionarily transposable from the ancestral, induced PGC specification program to an inherited PGC specification mechanism.

Chlorin Endogenous to the North Pacific Brittle Star Ophiura sarsii for Photodynamic Therapy Applications in Breast Cancer and Glioblastoma Models

Photodynamic therapy (PDT) represents a powerful avenue for anticancer treatment. PDT relies on the use of photosensitizers—compounds accumulating in the tumor and converted from benign to cytotoxic upon targeted photoactivation. We here describe (3S,4S)-14-Ethyl-9-(hydroxymethyl)-4,8,13,18-tetramethyl-20-oxo-3-phorbinepropanoic acid (ETPA) as a major metabolite of the North Pacific brittle stars Ophiura sarsii. As a chlorin, ETPA efficiently produces singlet oxygen upon red-light photoactivation and exerts powerful sub-micromolar phototoxicity against a panel of cancer cell lines in vitro. In a mouse model of glioblastoma, intravenous ETPA injection combined with targeted red laser irradiation induced strong necrotic ablation of the brain tumor. Along with the straightforward ETPA purification protocol and abundance of O. sarsii, these studies pave the way for the development of ETPA as a novel natural product-based photodynamic therapeutic.

Echinobase: leveraging an extant model organism database to build a knowledgebase supporting research on the genomics and biology of echinoderms

Echinobase (www.echinobase.org) is a third generation web resource supporting genomic research on echinoderms. The new version was built by cloning the mature Xenopus model organism knowledgebase, Xenbase, refactoring data ingestion pipelines and modifying the user interface to adapt to multispecies echinoderm content. This approach leveraged over 15 years of previous database and web application development to generate a new fully featured informatics resource in a single year. In addition to the software stack, Echinobase uses the private cloud and physical hosts that support Xenbase. Echinobase currently supports six echinoderm species, focused on those used for genomics, developmental biology and gene regulatory network analyses. Over 38 000 gene pages, 18 000 publications, new improved genome assemblies, JBrowse genome browser and BLAST + services are available and supported by the development of a new echinoderm anatomical ontology, uniformly applied formal gene nomenclature, and consistent orthology predictions. A novel feature of Echinobase is integrating support for multiple, disparate species. New genomes from the diverse echinoderm phylum will be added and supported as data becomes available. The common code development design of the integrated knowledgebases ensures parallel improvements as each resource evolves. This approach is widely applicable for developing new model organism informatics resources.

The tweety Gene Family: From Embryo to Disease

The tweety genes encode gated chloride channels that are found in animals, plants, and even simple eukaryotes, signifying their deep evolutionary origin. In vertebrates, the tweety gene family is highly conserved and consists of three members—ttyh1, ttyh2, and ttyh3—that are important for the regulation of cell volume. While research has elucidated potential physiological functions of ttyh1 in neural stem cell maintenance, proliferation, and filopodia formation during neural development, the roles of ttyh2 and ttyh3 are less characterized, though their expression patterns during embryonic and fetal development suggest potential roles in the development of a wide range of tissues including a role in the immune system in response to pathogen-associated molecules. Additionally, members of the tweety gene family have been implicated in various pathologies including cancers, particularly pediatric brain tumors, and neurodegenerative diseases such as Alzheimer’s and Parkinson’s disease. Here, we review the current state of research using information from published articles and open-source databases on the tweety gene family with regard to its structure, evolution, expression during development and adulthood, biochemical and cellular functions, and role in human disease. We also identify promising areas for further research to advance our understanding of this important, yet still understudied, family of genes.

FGF signalling plays similar roles in development and regeneration of the skeleton in the brittle star Amphiura filiformis

Regeneration as an adult developmental process is in many aspects similar to embryonic development. Although many studies point out similarities and differences, no large-scale, direct and functional comparative analyses between development and regeneration of a specific cell type or structure in one animal exist. Here, we use the brittle star Amphiura filiformis to characterise the role of the FGF signalling pathway during skeletal development in embryos and arm regeneration. In both processes, we find ligands expressed in ectodermal cells that flank underlying skeletal mesenchymal cells, which express the receptors. Perturbation of FGF signalling showed inhibited skeleton formation in both embryogenesis and regeneration, without affecting other key developmental processes. Differential transcriptome analysis finds mostly differentiation genes rather than transcription factors to be downregulated in both contexts. Moreover, comparative gene analysis allowed us to discover brittle star-specific differentiation genes. In conclusion, our results show that the FGF pathway is crucial for skeletogenesis in the brittle star, as in other deuterostomes, and provide evidence for the re-deployment of a developmental gene regulatory module during regeneration.

More Than One-to-Four via 2R: Evidence of an Independent Amphioxus Expansion and Two-Gene Ancestral Vertebrate State for MyoD-Related Myogenic Regulatory Factors (MRFs)

The evolutionary transition from invertebrates to vertebrates involved extensive gene duplication, but understanding precisely how such duplications contributed to this transition requires more detailed knowledge of specific cases of genes and gene families. Myogenic differentiation (MyoD) has long been recognized as a master developmental control gene and member of the MyoD family of bHLH transcription factors (myogenic regulatory factors [MRFs]) that drive myogenesis across the bilaterians. Phylogenetic reconstructions within this gene family are complicated by multiple instances of gene duplication and loss in several lineages. Following two rounds of whole-genome duplication (2R WGD) at the origin of the vertebrates, the ancestral function of MRFs is thought to have become partitioned among the daughter genes, so that MyoD and Myf5 act early in myogenic determination, whereas Myog and Myf6 are expressed later, in differentiating myoblasts. Comparing chordate MRFs, we find an independent expansion of MRFs in the invertebrate chordate amphioxus, with evidence for a parallel instance of subfunctionalization relative to that of vertebrates. Conserved synteny between chordate MRF loci supports the 2R WGD events as a major force in shaping the evolution of vertebrate MRFs. We also resolve vertebrate MRF complements and organization, finding a new type of vertebrate MRF gene in the process, which allowed us to infer an ancestral two-gene state in the vertebrates corresponding to the early- and late-acting types of MRFs. This necessitates a revision of previous conclusions about the simple one-to-four origin of vertebrate MRFs.

Na+/H+ exchangers differentially contribute to midgut fluid sodium and proton concentration in the sea urchin larva

Regulation of ionic composition and pH is a requisite of all digestive systems in the animal kingdom. Larval stages of the marine superphylum Ambulacraria, including echinoderms and hemichordates, were demonstrated to have highly alkaline conditions in their midgut with the underlying epithelial transport mechanisms being largely unknown. Using ion-selective microelectrodes, the present study demonstrated that pluteus larvae of the purple sea urchin have highly alkaline pH (pH ∼9) and low [Na+] (∼120 mmol l-1) in their midgut fluids, compared with the ionic composition of the surrounding seawater. We pharmacologically investigated the role of Na+/H+ exchangers (NHE) in intracellular pH regulation and midgut proton and sodium maintenance using the NHE inhibitor 5-(n-ethyl-n-isopropyl)amiloride (EIPA). Basolateral EIPA application decreased midgut pH while luminal application via micro-injections increased midgut [Na+], without affecting pH. Immunohistochemical analysis demonstrated a luminal localization of NHE-2 (SpSlc9a2) in the midgut epithelium. Specific knockdown of spslc9a2 using Vivo-Morpholinos led to an increase in midgut [Na+] without affecting pH. Acute acidification experiments in combination with quantitative PCR analysis and measurements of midgut pH and [Na+] identified two other NHE isoforms, Spslc9a7 and SpSlc9a8, which potentially contribute to the regulation of [Na+] and pH in midgut fluids. This work provides new insights into ion regulatory mechanisms in the midgut epithelium of sea urchin larvae. The involvement of NHEs in regulating pH and Na+ balance in midgut fluids shows conserved features of insect and vertebrate digestive systems and may contribute to the ability of sea urchin larvae to cope with changes in seawater pH.

CRISPR-Cas9 editing of non-coding genomic loci as a means of controlling gene expression in the sea urchin

We seek to manipulate gene function here through CRISPR-Cas9 editing of cis-regulatory sequences, rather than the more typical mutation of coding regions. This approach would minimize secondary effects of cellular responses to nonsense mediated decay pathways or to mutant protein products by premature stops. This strategy also allows for reducing gene activity in cases where a complete gene knockout would result in lethality, and it can be applied to the rapid identification of key regulatory sites essential for gene expression. We tested this strategy here with genes of known function as a proof of concept, and then applied it to examine the upstream genomic region of the germline gene Nanos2 in the sea urchin, Strongylocentrotus purpuratus. We first used CRISPR-Cas9 to target established genomic cis-regulatory regions of the skeletogenic cell transcription factor, Alx1, and the TGF-β signaling ligand, Nodal, which produce obvious developmental defects when altered in sea urchin embryos. Importantly, mutation of cis-activator sites (Alx1) and cis-repressor sites (Nodal) result in the predicted decreased and increased transcriptional output, respectively. Upon identification of efficient gRNAs by genomic mutations, we then used the same validated gRNAs to target a deadCas9-VP64 transcriptional activator to increase Nodal transcription directly. Finally, we paired these new methodologies with a more traditional, GFP reporter construct approach to further our understanding of the transcriptional regulation of Nanos2, a key gene required for germ cell identity in S. purpuratus. With a series of reporter assays, upstream Cas9-promoter targeted mutagenesis, coupled with qPCR and in situ RNA hybridization, we concluded that the promoter of Nanos2 drives strong mRNA expression in the sea urchin embryo, indicating that its primordial germ cell (PGC)-specific restriction may rely instead on post-transcriptional regulation. Overall, we present a proof-of-principle tool-kit of Cas9-mediated manipulations of promoter regions that should be applicable in most cells and embryos for which CRISPR-Cas9 is employed.

Rab35 regulates skeletogenesis and gastrulation by facilitating actin remodeling and vesicular trafficking

Rab35 is a small GTPase that regulates plasma membrane to early endosome vesicular trafficking and mediates actin remodeling to form actin-rich cellular structures. While the function of Rab35 in the cellular context has been examined, its role during development has not been well studied. In this study, we take advantage of the sea urchin's high fecundity, external fertilization, and transparent embryos to determine the function of Rab35 during development. We found that loss of function of Rab35 results in defects in skeletogenesis and gastrulation, which were rescued by co-injection of sea urchin Rab35. The loss of Rab35's function results in decreased endocytosis and impaired exocytosis, which may be important for skeletogenesis and gastrulation. Skeletal spicules of Rab35 knockdown embryos have reduced organized actin compared to the control, supporting the notion that Rab35 regulates actin dynamics. In addition, the skeletal and gastrulation defects induced by Rab35 knockdown were rescued by co-injection with Fascin, an actin-bundling protein, indicating that proper actin dynamics play a critical role for both skeletogenesis and gastrulation. Overall, results indicate that through its role in mediating vesicular trafficking and actin remodeling, Rab35 is an important regulator of embryonic structure formation in early development.

Systematic comparison of sea urchin and sea star developmental gene regulatory networks explains how novelty is incorporated in early development

The extensive array of morphological diversity among animal taxa represents the product of millions of years of evolution. Morphology is the output of development, therefore phenotypic evolution arises from changes to the topology of the gene regulatory networks (GRNs) that control the highly coordinated process of embryogenesis. A particular challenge in understanding the origins of animal diversity lies in determining how GRNs incorporate novelty while preserving the overall stability of the network, and hence, embryonic viability. Here we assemble a comprehensive GRN for endomesoderm specification in the sea star from zygote through gastrulation that corresponds to the GRN for sea urchin development of equivalent territories and stages. Comparison of the GRNs identifies how novelty is incorporated in early development. We show how the GRN is resilient to the introduction of a transcription factor, pmar1, the inclusion of which leads to a switch between two stable modes of Delta-Notch signaling. Signaling pathways can function in multiple modes and we propose that GRN changes that lead to switches between modes may be a common evolutionary mechanism for changes in embryogenesis. Our data additionally proposes a model in which evolutionarily conserved network motifs, or kernels, may function throughout development to stabilize these signaling transitions.

Extracellular matrix gene expression during arm regeneration in Amphiura filiformis

Extracellular matrix (ECM) plays a dynamic role during tissue development and re-growth. Body part regeneration efficiency relies also on effective ECM remodelling and deposition. Among invertebrates, echinoderms are well known for their striking regenerative abilities since they can rapidly regenerate functioning complex structures. To gather insights on the involvement of ECM during arm regeneration, the brittle star Amphiura filiformis was chosen as experimental model. Eight ECM genes were identified and cloned, and their spatio-temporal and quantitative expression patterns were analysed by means of whole mount in situ hybridisation and quantitative PCR on early and advanced regenerative stages. Our results show that almost none of the selected ECM genes are expressed at early stages of regeneration, suggesting a delay in their activation that may be responsible for the high regeneration efficiency of these animals, as described for other echinoderms and in contrast to most vertebrates. Moreover, at advanced stages, these genes are spatially and temporally differentially expressed, suggesting that the molecular regulation of ECM deposition/remodelling varies throughout the regenerative process. Phylogenetic analyses of the identified collagen-like genes reveal complex evolutionary dynamics with many rounds of duplications and losses and pinpointed their homologues in selected vertebrates. The study of other ECM genes will allow a better understanding of ECM contribution to brittle star arm regeneration.

Ocean acidification promotes broad transcriptomic responses in marine metazoans: a literature survey

For nearly a decade, the metazoan-focused research community has explored the impacts of ocean acidification (OA) on marine animals, noting that changes in ocean chemistry can impact calcification, metabolism, acid-base regulation, stress response and behavior in organisms that hold high ecological and economic value. Because OA interacts with several key physiological processes in marine organisms, transcriptomics has become a widely-used method to characterize whole organism responses on a molecular level as well as inform mechanisms that explain changes in phenotypes observed in response to OA. In the past decade, there has been a notable rise in studies that examine transcriptomic responses to OA in marine metazoans, and here we attempt to summarize key findings across these studies. We find that organisms vary dramatically in their transcriptomic responses to pH although common patterns are often observed, including shifts in acid-base ion regulation, metabolic processes, calcification and stress response mechanisms. We also see a rise in transcriptomic studies examining organismal response to OA in a multi-stressor context, often reporting synergistic effects of OA and temperature. In addition, there is an increase in studies that use transcriptomics to examine the evolutionary potential of organisms to adapt to OA conditions in the future through population and transgenerational experiments. Overall, the literature reveals complex organismal responses to OA, in which some organisms will face more dramatic consequences than others. This will have wide-reaching impacts on ocean communities and ecosystems as a whole.

Improved protocols for BAC insert DNA isolation, BAC end sequencing and FISH for construction of BAC based physical map of genes on the chromosomes

Bacterial artificial chromosome (BAC) library is an important genomic resource useful in targeted marker development, positional cloning, physical mapping and a substrate for genome sequencing for better understanding the genome organization of a species. The present manuscript elucidates the improvement in protocols for economical and efficient BAC insert DNA isolation, BAC end sequencing and FISH for physical localization on the metaphase chromosome complements. BAC clones of Clarias magur, maintained in 384-well plate format in our laboratory, were used in this study. The protocols gave consistent and efficient results. We use routinely these protocols for BAC insert DNA extraction, generating end sequence data of the clone and constructing DNA probes to hybridize on the metaphase spreads of C. magur using FISH for physical their localization.

Developmental transcriptomes of the sea star, Patiria miniata, illuminate how gene expression changes with evolutionary distance

Understanding how changes in developmental gene expression alter morphogenesis is a fundamental problem in development and evolution. A promising approach to address this problem is to compare the developmental transcriptomes between related species. The echinoderm phylum consists of several model species that have significantly contributed to the understanding of gene regulation and evolution. Particularly, the regulatory networks of the sea star, Patiria miniata (P. miniata), have been extensively studied, however developmental transcriptomes for this species were lacking. Here we generated developmental transcriptomes of P. miniata and compared these with those of two sea urchins species. We demonstrate that the conservation of gene expression depends on gene function, cell type and evolutionary distance. With increasing evolutionary distance the interspecies correlations in gene expression decreases. The reduction is more severe in the correlations between morphologically equivalent stages (diagonal elements) than in the correlation between morphologically distinct stages (off-diagonal elements). This could reflect a decrease in the morphological constraints compared to other constraints that shape gene expression at large evolutionary divergence. Within this trend, the interspecies correlations of developmental control genes maintain their diagonality at large evolutionary distance, and peak at the onset of gastrulation, supporting the hourglass model of phylotypic stage conservation.

Skeletal development in the sea urchin relies upon protein families that contain intrinsic disorder, aggregation-prone, and conserved globular interactive domains

The formation of the sea urchin spicule skeleton requires the participation of hydrogel-forming protein families that regulate mineral nucleation and nanoparticle assembly processes that give rise to the spicule. However, the structure and molecular behavior of these proteins is not well established, and thus our ability to understand this process is hampered. We embarked on a study of sea urchin spicule proteins using a combination of biophysical and bioinformatics techniques. Our biophysical findings indicate that recombinant variants of the two most studied spicule matrix proteins, SpSM50 and SpSM30B/C (S. purpuratus) have a conformational landscape that include a C-terminal random coil/intrinsically disordered MAPQG sequence coupled to a conserved, folded N-terminal C-type lectin-like (CTLL) domain, with SpSM50 > SpSM30B/C with regard to intrinsic disorder. Both proteins possess solvent-accessible unfolded MAQPG sequence regions where Asn, Gln, and Arg residues may be accessible for protein hydrogel interactions with water molecules. Our bioinformatics study included seven other spicule matrix proteins where we note similarities between these proteins and rare, unusual proteins that possess folded and unfolded traits. Moreover, spicule matrix proteins possess three types of sequences: intrinsically disordered, amyloid-like, and folded protein-protein interactive. Collectively these reactive domains would be capable of driving protein assembly and hydrogel formation. Interestingly, three types of global conformations are predicted for the nine member protein set, wherein we note variations in the arrangement of intrinsically disordered and interactive globular domains. These variations may reflect species-specific requirements for spiculogenesis. We conclude that the molecular landscape of spicule matrix protein families enables them to function as hydrogelators, nucleators, and assemblers of mineral nanoparticles.

Pervasive convergent evolution and extreme phenotypes define chaperone requirements of protein homeostasis

Maintaining protein homeostasis is an essential requirement for cell and organismal viability. An elaborate regulatory system within cells, the protein homeostasis network, safeguards that proteins are correctly folded and functional. At the heart of this regulatory system lies a class of specialized protein quality control enzymes called chaperones that are tasked with assisting proteins in their folding, avoiding aggregation and degradation. Failure and decline of protein homeostasis are directly associated with conditions of aging and aging-related neurodegeneration. However, it is not clear what tips the balance of protein homeostasis and leads to onset of aging and diseases. Here, using a comparative genomics approach we report general principles of maintaining protein homeostasis across the eukaryotic tree of life. Expanding a previous study of 16 eukaryotes to the quantitative analysis of 216 eukaryotic genomes, we find a strong correlation between the composition of eukaryotic chaperone networks and genome complexity that is distinct for different species kingdoms. Organisms with pronounced phenotypes clearly buck this trend. Northobranchius furzeri, the shortest-lived vertebrate and a widely used model for fragile protein homeostasis, is found to be chaperone limited while Heterocephalus glaber as the longest-lived rodent and thus an especially robust organism is characterized by above-average numbers of chaperones. Strikingly, the relative size of chaperone networks is found to generally correlate with longevity in Metazoa. Our results thus indicate that the balance in protein homeostasis may be a key variable in explaining organismal robustness.

Drosophila menthol sensitivity and the Precambrian origins of transient receptor potential-dependent chemosensation

Transient receptor potential (TRP) cation channels are highly conserved, polymodal sensors which respond to a wide variety of stimuli. Perhaps most notably, TRP channels serve critical functions in nociception and pain. A growing body of evidence suggests that transient receptor potential melastatin (TRPM) and transient receptor potential ankyrin (TRPA) thermal and electrophile sensitivities predate the protostome-deuterostome split (greater than 550 Ma). However, TRPM and TRPA channels are also thought to detect modified terpenes (e.g. menthol). Although terpenoids like menthol are thought to be aversive and/or harmful to insects, mechanistic sensitivity studies have been largely restricted to chordates. Furthermore, it is unknown if TRP-menthol sensing is as ancient as thermal and/or electrophile sensitivity. Combining genetic, optical, electrophysiological, behavioural and phylogenetic approaches, we tested the hypothesis that insect TRP channels play a conserved role in menthol sensing. We found that topical application of menthol to Drosophila melanogaster larvae elicits a Trpm- and TrpA1-dependent nocifensive rolling behaviour, which requires activation of Class IV nociceptor neurons. Further, in characterizing the evolution of TRP channels, we put forth the hypotheses that three previously undescribed TRPM channel clades (basal, αTRPM and βTRPM), as well as TRPs with residues critical for menthol sensing, were present in ancestral bilaterians. This article is part of the Theo Murphy meeting issue 'Evolution of mechanisms and behaviour important for pain'.

CNEr: A toolkit for exploring extreme noncoding conservation

Conserved Noncoding Elements (CNEs) are elements exhibiting extreme noncoding conservation in Metazoan genomes. They cluster around developmental genes and act as long-range enhancers, yet nothing that we know about their function explains the observed conservation levels. Clusters of CNEs coincide with topologically associating domains (TADs), indicating ancient origins and stability of TAD locations. This has suggested further hypotheses about the still elusive origin of CNEs, and has provided a comparative genomics-based method of estimating the position of TADs around developmentally regulated genes in genomes where chromatin conformation capture data is missing. To enable researchers in gene regulation and chromatin biology to start deciphering this phenomenon, we developed CNEr, a R/Bioconductor toolkit for large-scale identification of CNEs and for studying their genomic properties. We apply CNEr to two novel genome comparisons—fruit fly vs tsetse fly, and two sea urchin genomes—and report novel insights gained from their analysis. We also show how to reveal interesting characteristics of CNEs by coupling CNEr with existing Bioconductor packages. CNEr is available at Bioconductor (https://bioconductor.org/packages/CNEr/) and maintained at github (https://github.com/ge11232002/CNEr).

Asymmetric oceanographic processes mediate connectivity and population genetic structure, as revealed by RADseq, in a highly dispersive marine invertebrate (Parastichopus californicus)

Marine populations are typically characterized by weak genetic differentiation due to the potential for long-distance dispersal favouring high levels of gene flow. However, strong directional advection of water masses or retentive hydrodynamic forces can influence the degree of genetic exchange among marine populations. To determine the oceanographic drivers of genetic structure in a highly dispersive marine invertebrate, the giant California sea cucumber (Parastichopus californicus), we first tested for the presence of genetic discontinuities along the coast of North America in the northeastern Pacific Ocean. Then, we tested two hypotheses regarding spatial processes influencing population structure: (i) isolation by distance (IBD: genetic structure is explained by geographic distance) and (ii) isolation by resistance (IBR: genetic structure is driven by ocean circulation). Using RADseq, we genotyped 717 individuals from 24 sampling locations across 2,719 neutral SNPs to assess the degree of population differentiation and integrated estimates of genetic variation with inferred connectivity probabilities from a biophysical model of larval dispersal mediated by ocean currents. We identified two clusters separating north and south regions, as well as significant, albeit weak, substructure within regions (F_ST  = 0.002, p = .001). After modelling the asymmetric nature of ocean currents, we demonstrated that local oceanography (IBR) was a better predictor of genetic variation (R²  = .49) than geographic distance (IBD) (R²  = .18), and directional processes played an important role in shaping fine-scale structure. Our study contributes to the growing body of literature identifying significant population structure in marine systems and has important implications for the spatial management of P. californicus and other exploited marine species.

cis-Regulatory analysis for later phase of anterior neuroectoderm-specific foxQ2 expression in sea urchin embryos

The specification of anterior neuroectoderm is controlled by a highly conserved molecular mechanism in bilaterians. A forkhead family gene, foxQ2, is known to be one of the pivotal regulators, which is zygotically expressed in this region during embryogenesis of a broad range of bilaterians. However, what controls the expression of this essential factor has remained unclear to date. To reveal the regulatory mechanism of foxQ2, we performed cis-regulatory analysis of two foxQ2 genes, foxQ2a and foxQ2b, in a sea urchin Hemicentrotus pulcherrimus. In sea urchin embryos, foxQ2 is initially expressed in the entire animal hemisphere and subsequently shows narrower expression restricted to the anterior pole region. In this study, as a first step to understand the foxQ2 regulation, we focused on the later restricted expression and analyzed the upstream cis-regulatory sequences of foxQ2a and foxQ2b by using the constructs fused to short half-life green fluorescent protein. Based on deletion and mutation analyses of both foxQ2, we identified each of the five regulatory sequences, which were 4-9 bp long. Neither of the regulatory sequences contains any motifs for ectopic activation or spatial repression, suggesting that later mRNA localization is regulated in situ. We also suggest that the three amino acid loop extension-class homeobox gene Meis is involved in the maintenance of foxQ2b, the expression of which is dominant during embryogenesis.

Reference assembly and gene expression analysis of Apostichopus japonicus larval development

The transcriptome of the holothurian Apostichopus japonicus was sequenced at four developmental stages-blastula, gastrula, auricularia, pentactula-on an Illumina sequencer. Based on our RNA-seq data and the paired-end reads from 16 libraries obtained by other researchers earlier, we have achieved the currently most complete transcriptome assembly for A. japonicus with the best basic statistical parameters. An analysis of the obtained transcriptome has revealed 174 differentially expressed transcription factors, as well as stage-specific transcription factors that are most promising for further study. In addition, a total of 1,174,999 high-quality single nucleotide polymorphisms have been identified, including 58,932 indels. A GO enrichment analysis of contigs containing polymorphic loci shows the predominance of GO terms associated with immune response. The data obtained by us provide an additional basis for a deeper study of the mechanisms of the planktotrophic-type development in holothurians and can be used in commercial sea cucumber breeding programs.

Developmental effector gene regulation: Multiplexed strategies for functional analysis

The staggering complexity of the genome controls for developmental processes is revealed through massively parallel cis-regulatory analysis using new methods of perturbation and readout. The choice of combinations of these new methods is tailored to the system, question and resources at hand. Our focus is on issues that include the necessity or sufficiency of given cis-regulatory modules, cis-regulatory function in the normal spatial genomic context, and easily accessible high throughput and multiplexed analysis methods. In the sea urchin embryonic model, recombineered BACs offer new opportunities for consecutive modes of cis-regulatory analyses that answer these requirements, as we here demonstrate on a diverse suite of previously unstudied sea urchin effector genes expressed in skeletogenic cells. Positively active cis-regulatory modules were located in single Nanostring experiments per BAC containing the gene of interest, by application of our previously reported "barcode" tag vectors of which> 100 can be analyzed at one time. Computational analysis of DNA sequences that drive expression, based on the known skeletogenic regulatory state, then permitted effective identification of functional target site clusters. Deletion of these sub-regions from the parent BACs revealed module necessity, as simultaneous tests of the same regions in short constructs revealed sufficiency. Predicted functional inputs were then confirmed by site mutations, all generated and tested in multiplex formats. There emerged the simple conclusion that each effector gene utilizes a small subset of inputs from the skeletogenic GRN. These inputs may function to only adjust expression levels or in some cases necessary for expression. Since we know the GRN architecture upstream of the effector genes, we could then conceptually isolate and compare the wiring of the effector gene driver sub-circuits and identify the inputs whose removal abolish expression.

CRISPR/Cas9-mediated genome editing in sea urchins

The CRISPR (clustered regularly interspaced short palindromic repeat)/Cas9 (CRISPR-associated nuclease 9) technology enables rapid, targeted, and efficient changes in the genomes of various model organisms. The short guide RNAs (gRNAs) of the CRISPR/Cas9 system can be designed to recognize target DNA within coding regions for functional gene knockouts. Several studies have demonstrated that the CRISPR/Cas9 system efficiently and specifically targets sea urchin genes and results in expected mutant phenotypes. In addition to disrupting gene functions, modifications and additions to the Cas9 protein enable alternative activities targeted to specific sites within the genome. This includes a fusion of cytidine deaminase to Cas9 (Cas9-DA) for single nucleotide conversion in targeted sites. In this chapter, we describe detailed methods for the CRISPR/Cas9 application in sea urchin embryos, including gRNA design, in vitro synthesis of single guide RNA (sgRNA), and the usages of the CRISPR/Cas9 technology for gene knockout and single nucleotide editing. Methods for genotyping the resultant embryos are also provided for assessing efficiencies of gene editing.

The crowns have eyes: multiple opsins found in the eyes of the crown-of-thorns starfish Acanthaster planci

Background

Opsins are G protein-coupled receptors used for both visual and non-visual photoreception, and these proteins evolutionarily date back to the base of the bilaterians. In the current sequencing age, phylogenomic analysis has proven to be a powerful tool, facilitating the increase in knowledge about diversity within the opsin subclasses and, so far, at least nine types of opsins have been identified. Within echinoderms, opsins have been studied in Echinoidea and Ophiuroidea, which do not possess proper image forming eyes, but rather widely dispersed dermal photoreceptors. However, most species of Asteroidea, the starfish, possess true eyes and studying them will shed light on the diversity of opsin usage within echinoderms and help resolve the evolutionary history of opsins.

Results

Using high-throughput RNA sequencing, we have sequenced and analyzed the transcriptomes of different Acanthaster planci tissue samples: eyes, radial nerve, tube feet and a mixture of tissues from other organs. At least ten opsins were identified, and eight of them were found significantly differentially expressed in both eyes and radial nerve, with R-opsin being the most highly expressed in the eye.

Conclusion

This study provides new important insight into the involvement of opsins in visual and nonvisual photoreception. Of relevance, we found the first indication of an r-opsin photopigment expressed in a well-developed visual eye in a deuterostome animal. Additionally, we provided tissue specific A. planci transcriptomes that will aid in future Evo Devo studies.

Electronic supplementary material

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

Characterization of the TLR Family in Branchiostoma lanceolatum and Discovery of a Novel TLR22-Like Involved in dsRNA Recognition in Amphioxus

Toll-like receptors (TLRs) are important for raising innate immune responses in both invertebrates and vertebrates. Amphioxus belongs to an ancient chordate lineage which shares key features with vertebrates. The genomic research on TLR genes in Branchiostoma floridae and Branchiostoma belcheri reveals the expansion of TLRs in amphioxus. However, the repertoire of TLRs in Branchiostoma lanceolatum has not been studied and the functionality of amphioxus TLRs has not been reported. We have identified from transcriptomic data 30 new putative TLRs in B. lanceolatum and all of them are transcribed in adult amphioxus. Phylogenetic analysis showed that the repertoire of TLRs consists of both non-vertebrate and vertebrate-like TLRs. It also indicated a lineage-specific expansion in orthologous clusters of the vertebrate TLR11 family. We did not detect any representatives of the vertebrate TLR1, TLR3, TLR4, TLR5 and TLR7 families. To gain insight into these TLRs, we studied in depth a particular TLR highly similar to a B. belcheri gene annotated as bbtTLR1. The phylogenetic analysis of this novel BlTLR showed that it clusters with the vertebrate TLR11 family and it might be more related to TLR13 subfamily according to similar domain architecture. Transient and stable expression in HEK293 cells showed that the BlTLR localizes on the plasma membrane, but it did not respond to the most common mammalian TLR ligands. However, when the ectodomain of BlTLR is fused to the TIR domain of human TLR2, the chimeric protein could indeed induce NF-κB transactivation in response to the viral ligand Poly I:C, also indicating that in amphioxus, specific accessory proteins are needed for downstream activation. Based on the phylogenetic, subcellular localization and functional analysis, we propose that the novel BlTLR might be classified as an antiviral receptor sharing at least partly the functions performed by vertebrate TLR22. TLR22 is thought to be viral teleost-specific TLR but here we demonstrate that teleosts and amphioxus TLR22-like probably shared a common ancestor. Additional functional studies with other lancelet TLR genes will enrich our understanding of the immune response in amphioxus and will provide a unique perspective on the evolution of the immune system.

Secrets of the Sea Urchin Spicule Revealed: Protein Cooperativity Is Responsible for ACC Transformation, Intracrystalline Incorporation, and Guided Mineral Particle Assembly in Biocomposite Material Formation

The formation of the sea urchin spicule involves the stabilization and transformation of amorphous calcium carbonate (ACC) and assembly of ACC nanoparticle precursors into a mesoscale single crystal of fracture-resistant calcite. This process of particle assembly or attachment is under the control of a family of proteins known as the spicule matrix [Strongylocentrotus purpuratus (SpSM)] proteome. Recently, two members of this proteome, SpSM50 and the glycoprotein SpSM30B/C-G (in recombinant forms), were found to interact together via SpSM30B/C-G oligosaccharide-SpSM50 protein interactions to form hybrid protein hydrogels with unique physical properties. In this study, we investigate the mineralization properties of this hybrid hydrogel alongside the hydrogels formed by SpSM50 and SpSM30B/C-G individually. We find that the SpSM50 + SpSM30B/C-G hybrid hydrogel is synergistic with regard to surface modifications and intracrystalline inclusions of existing calcite crystals, the inhibition of ACC formation, and the kinetic destabilization of ACC to form a crystalline phase. Most importantly, the hybrid hydrogel phase assembles and organizes mineral particles into discrete clusters or domains within in vitro mineralization environments. Thus, the interactions of SpSM50 and SpSM30B/C-G, mediated by carbohydrate-protein binding, reflect the need for protein cooperativity for the ACC-to-crystalline transformation, intracrystalline void formation, and guided mineral particle assembly processes that are instrumental in spicule formation.

PdumBase: a transcriptome database and research tool for Platynereis dumerilii and early development of other metazoans

BACKGROUND

The marine polychaete annelid Platynereis dumerilii has recently emerged as a prominent organism for the study of development, evolution, stem cells, regeneration, marine ecology, chronobiology and neurobiology within metazoans. Its phylogenetic position within the spiralian/ lophotrochozoan clade, the comparatively high conservation of ancestral features in the Platynereis genome, and experimental access to any stage within its life cycle, make Platynereis an important model for elucidating the complex regulatory and functional molecular mechanisms governing early development, later organogenesis, and various features of its larval and adult life. High resolution RNA-seq gene expression data obtained from specific developmental stages can be used to dissect early developmental mechanisms. However, the potential for discovery of these mechanisms relies on tools to search, retrieve, and compare genome-wide information within Platynereis, and across other metazoan taxa.

RESULTS

To facilitate exploration and discovery by the broader scientific community, we have developed a web-based, searchable online research tool, PdumBase, featuring the first comprehensive transcriptome database for Platynereis dumerilii during early stages of development (2 h ~ 14 h). Our database also includes additional stages over the P. dumerilii life cycle and provides access to the expression data of 17,213 genes (31,806 transcripts) along with annotation information sourced from Swiss-Prot, Gene Ontology, KEGG pathways, Pfam domains, TmHMM, SingleP, and EggNOG orthology. Expression data for each gene includes the stage, the normalized FPKM, the raw read counts, and information that can be leveraged for statistical analyses of differential gene expression and the construction of genome-wide co-expression networks. In addition, PdumBase offers early stage transcriptome expression data from five further species as a valuable resource for investigators interested in comparing early development in different organisms. To understand conservation of Platynereis gene models and to validate gene annotation, most Platynereis gene models include a comprehensive phylogenetic analysis across 18 species representing diverse metazoan taxa.

CONCLUSIONS

PdumBase represents the first online resource for the early developmental transcriptome of Platynereis dumerilii. It serves as a research platform for discovery and exploration of gene expression during early stages, throughout the Platynereis life cycle, and enables comparison to other model organisms. PdumBase is freely available at http://pdumbase.gdcb.iastate.edu .

Gene expression across tissues, sex, and life stages in the sea urchin Tripneustes gratilla [Toxopneustidae, Odontophora, Camarodonta]

The pan-tropical sea urchin Tripneustes gratilla is an ecologically and economically important shallow water algal grazer. The aquaculture of T. gratilla has spurred growing interest in the population biology of the species, and by extension the generation of more molecular resources. To this purpose, de novo transcriptomes of T. gratilla were generated for two adults, a male and a female, as well as for a cohort of approximately 1000 plutei larvae. Gene expression profiles of three adult tissue samples were quantified and compared. These samples were of gonadal tissue, the neural ring, and pooled tube feet and pedicellariae. Levels of shared and different gene expression between sexes, as well as across functional categories of interest, including the immune system, toxins, genes involved in fertilization, and sensory genes are highlighted. Differences in expression of isoforms between the sexes and Sex determining Region Y-related High Mobility Group box groups is observed. Additionally an expansion of the tumor suppressor DMBT1 is observed in T. gratilla when compared to the annotated genome of the sea urchin Strongylocentrotus purpuratus. The draft transcriptome of T. gratilla is presented here in order to facilitate more genomic level analysis of emerging model sea urchin systems.

Echinobase: an expanding resource for echinoderm genomic information

Echinobase, a web accessible information system of diverse genomics and biological data for the echinoderm clade, grew out of SpBase, the first echinoderm genome project for sea urchin, Strongylocentrotus purpuratus. Sea urchins and their relatives are utilitarian research models in fields ranging from marine biology to developmental biology and gene regulatory systems. Echinobase is a user-friendly web interface that links an array of biological data that would otherwise have been tedious and frustrating for researchers to extract and organize. The system hosts a powerful gene search engine, genomics browser and other bioinformatics tools to investigate genomics and high throughput data. The Echinobase information system now serves genomic information for eight echinoderm species: S. purpuratus, Strongylocentrotus fransciscanus, Allocentrotus fragilis, Lytechinus variegatus, Patiria miniata, Parastichopus parvimensis and Ophiothrix spiculata, Eucidaris tribuloides. Herein lies a description of the web information system, genomics data types and content hosted by Echinobase.org. The goal of Echinobase is to connect genomic information to various experimental data and accelerate the research in field of molecular biology, developmental process, gene regulatory networks and more recently engineering biological systems0.

Database URL:http://www.echinobase.org

Glycosylation Fosters Interactions between Model Sea Urchin Spicule Matrix Proteins. Implications for Embryonic Spiculogenesis and Biomineralization

The formation of embryonic mineralized skeletal elements (spicules) in the sea urchin requires the participation of proteins, many of which may interact with one another and assist in the creation of an extracellular matrix wherein mineral formation takes place. To probe this, we created a sea urchin spicule recombinant model protein pair system wherein we tested the interactions between two major spicule proteins, SpSM50 and the glycoprotein, SpSM30B/C. Both proteins are strong hydrogelators that manipulate early and later events in mineral formation. We discovered that the anionic glycan moieties of SpSM30B/C are required for interaction with the SpSM50 protein and that these interactions are Ca(II)-independent. In addition, when these proteins form a complex, they create hybrid hydrogel particles that are physically distinct from their individual counterparts. Thus, glycan-mediated interactions play an important role in in vitro spicule protein assembly and most likely within the spicule itself.

Omics approaches to study gene regulatory networks for development in echinoderms

Gene regulatory networks (GRNs) describe the interactions for a developmental process at a given time and space. Historically, perturbation experiments represent one of the key methods for analyzing and reconstructing a GRN, and the GRN governing early development in the sea urchin embryo stands as one of the more deeply dissected so far. As technology progresses, so do the methods used to address different biological questions. Next-generation sequencing (NGS) has become a standard experimental technique for genome and transcriptome sequencing and studies of protein-DNA interactions and DNA accessibility. While several efforts have been made toward the integration of different omics approaches for the study of the regulatory genome in many animals, in a few cases, these are applied with the purpose of reconstructing and experimentally testing developmental GRNs. Here, we review emerging approaches integrating multiple NGS technologies for the prediction and validation of gene interactions within echinoderm GRNs. These approaches can be applied to both 'model' and 'non-model' organisms. Although a number of issues still need to be addressed, advances in NGS applications, such as assay for transposase-accessible chromatin sequencing, combined with the availability of embryos belonging to different species, all separated by various evolutionary distances and accessible to experimental regulatory biology, place echinoderms in an unprecedented position for the reconstruction and evolutionary comparison of developmental GRNs. We conclude that sequencing technologies and integrated omics approaches allow the examination of GRNs on a genome-wide scale only if biological perturbation and cis-regulatory analyses are experimentally accessible, as in the case of echinoderm embryos.

Discovery of novel representatives of bilaterian neuropeptide families and reconstruction of neuropeptide precursor evolution in ophiuroid echinoderms

Neuropeptides are a diverse class of intercellular signalling molecules that mediate neuronal regulation of many physiological and behavioural processes. Recent advances in genome/transcriptome sequencing are enabling identification of neuropeptide precursor proteins in species from a growing variety of animal taxa, providing new insights into the evolution of neuropeptide signalling. Here, detailed analysis of transcriptome sequence data from three brittle star species, Ophionotus victoriae, Amphiura filiformis and Ophiopsila aranea, has enabled the first comprehensive identification of neuropeptide precursors in the class Ophiuroidea of the phylum Echinodermata. Representatives of over 30 bilaterian neuropeptide precursor families were identified, some of which occur as paralogues. Furthermore, homologues of endothelin/CCHamide, eclosion hormone, neuropeptide-F/Y and nucleobinin/nesfatin were discovered here in a deuterostome/echinoderm for the first time. The majority of ophiuroid neuropeptide precursors contain a single copy of a neuropeptide, but several precursors comprise multiple copies of identical or non-identical, but structurally related, neuropeptides. Here, we performed an unprecedented investigation of the evolution of neuropeptide copy number over a period of approximately 270 Myr by analysing sequence data from over 50 ophiuroid species, with reference to a robust phylogeny. Our analysis indicates that the composition of neuropeptide ‘cocktails’ is functionally important, but with plasticity over long evolutionary time scales.

The developmental transcriptomes of two sea biscuit species with differing larval types

BACKGROUND

Larval developmental patterns are extremely varied both between and within phyla, however the genetic mechanisms leading to this diversification are poorly understood. We assembled and compared the developmental transcriptomes for two sea biscuit species (Echinodermata: Echinoidea) with differing patterns of larval development, to provide a resource for investigating the evolution of alternate life cycles. One species (Clypeaster subdepressus) develops via an obligately feeding larva which metamorphoses 3-4 weeks after fertilization; the other (Clypeaster rosaceus) develops via a rare, intermediate larval type-facultative feeding- and can develop through metamorphosis entirely based on egg provisioning in under one week.

RESULTS

Overall, the two transcriptomes are highly similar, containing largely orthologous contigs with similar functional annotation. However, we found distinct differences in gene expression patterns between the two species. Larvae from C. rosaceus, the facultative planktotroph, turned genes on at earlier stages and had less differentiation in gene expression between larval stages, whereas, C. subdepressus showed a higher degree of stage-specific gene expression.

CONCLUSION

This study is the first genetic analysis of a species with facultatively feeding larvae. Our results are consistent with known developmental differences between the larval types and raise the question of whether earlier onset of developmental genes is a key step in the evolution of a reduced larval period. By publishing a transcriptome for this rare, intermediate, larval type, this study adds developmental breadth to the current genetic resources, which will provide a valuable tool for future research on echinoderm development as well as studies on the evolution of development in general.

Translatome analysis at the egg-to-embryo transition in sea urchin

Early embryogenesis relies on the translational regulation of maternally stored mRNAs. In sea urchin, fertilization triggers a dramatic rise in translation activity, necessary for the onset of cell division. Here, the full spectrum of the mRNAs translated upon fertilization was investigated by polysome profiling and sequencing. The translatome of the early sea urchin embryo gave a complete picture of the polysomal recruitment dynamics following fertilization. Our results indicate that only a subset of maternal mRNAs were selectively recruited onto polysomes, with over-represented functional categories in the translated set. The increase in translation upon fertilization depends on the formation of translation initiation complexes following mTOR pathway activation. Surprisingly, mTOR pathway inhibition differentially affected polysomal recruitment of the newly translated mRNAs, which thus appeared either mTOR-dependent or mTOR-independent. Therefore, our data argue for an alternative to the classical cap-dependent model of translation in early development. The identification of the mRNAs translated following fertilization helped assign translational activation events to specific mRNAs. This translatome is the first step to a comprehensive analysis of the molecular mechanisms governing translation upon fertilization and the translational regulatory networks that control the egg-to-embryo transition as well as the early steps of embryogenesis.

AID/APOBEC-like cytidine deaminases are ancient innate immune mediators in invertebrates

In the course of both innate and adaptive immunity, cytidine deaminases within the activation induced cytidine deaminase (AID)/apolipoprotein B editing complex (APOBEC) family modulate immune responses by mutating specific nucleic acid sequences of hosts and pathogens. The evolutionary emergence of these mediators, however, seems to coincide precisely with the emergence of adaptive immunity in vertebrates. Here, we show a family of genes in species within two divergent invertebrate phyla-the echinoderm Strongylocentrotus purpuratus and the brachiopod Lingula anatina-that encode proteins with similarities in amino acid sequence and enzymatic activities to the vertebrate AID/APOBECs. The expression of these invertebrate factors is enriched in tissues undergoing constant, direct interactions with microbes and can be induced upon pathogen challenge. Our findings suggest that AID/APOBEC proteins, and their function in immunity, emerged far earlier than previously thought. Thus, cytidine deamination is probably an ancient innate immune mechanism that predates the protostome/deuterostome divergence.

Genome-wide use of high- and low-affinity Tbrain transcription factor binding sites during echinoderm development

Sea stars and sea urchins are model systems for interrogating the types of deep evolutionary changes that have restructured developmental gene regulatory networks (GRNs). Although cis-regulatory DNA evolution is likely the predominant mechanism of change, it was recently shown that Tbrain, a Tbox transcription factor protein, has evolved a changed preference for a low-affinity, secondary binding motif. The primary, high-affinity motif is conserved. To date, however, no genome-wide comparisons have been performed to provide an unbiased assessment of the evolution of GRNs between these taxa, and no study has attempted to determine the interplay between transcription factor binding motif evolution and GRN topology. The study here measures genome-wide binding of Tbrain orthologs by using ChIP-sequencing and associates these orthologs with putative target genes to assess global function. Targets of both factors are enriched for other regulatory genes, although nonoverlapping sets of functional enrichments in the two datasets suggest a much diverged function. The number of low-affinity binding motifs is significantly depressed in sea urchins compared with sea star, but both motif types are associated with genes from a range of functional categories. Only a small fraction (∼10%) of genes are predicted to be orthologous targets. Collectively, these data indicate that Tbr has evolved significantly different developmental roles in these echinoderms and that the targets and the binding motifs in associated cis-regulatory sequences are dispersed throughout the hierarchy of the GRN, rather than being biased toward terminal process or discrete functional blocks, which suggests extensive evolutionary tinkering.

Molecular mechanisms of fission in echinoderms: Transcriptome analysis

Echinoderms are capable of asexual reproduction by fission. An individual divides into parts due to changes in the strength of connective tissue of the body wall. The structure of connective tissue and the mechanisms of variations in its strength in echinoderms remain poorly studied. An analysis of transcriptomes of individuals during the process of fission provides a new opportunity to understand the mechanisms of connective tissue mutability. In the holothurian Cladolabes schmeltzii, we have found a rather complex organization of connective tissue. Transcripts of genes encoding a wide range of structural proteins of extracellular matrix, as well as various proteases and their inhibitors, have been discovered. All these molecules may constitute a part of the mechanism of connective tissue mutability. According to our data, the extracellular matrix of echinoderms is substantially distinguished from that of vertebrates by the lack of elastin, fibronectins, and tenascins. In case of fission, a large number of genes of transcription factors and components of different signaling pathways are expressed. Products of these genes are probably involved in regulation of asexual reproduction, connective tissue mutability, and preparation of tissues for subsequent regeneration. It has been shown that holothurian tensilins are a special group of tissue inhibitors of metalloproteinases, which has formed within the class Holothuroidea and is absent from other echinoderms. Our data can serve a basis for the further study of the mechanisms of extracellular matrix mutability, as well as the mechanisms responsible for asexual reproduction in echinoderms.

Toxigenic effects of two benthic diatoms upon grazing activity of the sea urchin: morphological, metabolomic and de novo transcriptomic analysis

Diatoms are unicellular algae playing a key role as photosynthetic organisms in the world's ocean food webs. The chemical ecology of planktonic diatoms is well documented, but few studies have reported on the effects of benthic diatoms on their consumers, also due to difficulties in the collection, quantification and massive culturing of benthic species. Here for the first time we investigate the effects of feeding on two abundantly occurring benthic diatoms, Nanofrustulum shiloi and Cylindrotheca closterium, isolated from the leaves of the seagrass Posidonia oceanica, on the sea urchin Paracentrotus lividus. Adult P. lividus were fed for one month on diets of either one of the two diatoms and on the green alga Ulva rigida, used as a feeding control. By combining morphological, metabolomic and de novo transcriptomic approaches, we demonstrate toxigenic effect on embryos generated by females fed with these benthic diatoms. Furthermore, chemical analysis reveal the presence of polyunsaturated aldehydes only for N. shiloi, and a high production of other oxylipins (cytotoxic compounds on their grazers and on cancer cell lines) for both diatoms, including some additional peaks not correlated to the canonic oxylipins commonly observed in planktonic diatoms. These findings open new perspectives in the study of diatom secondary metabolites influencing their grazers.

Global analysis of primary mesenchyme cell cis-regulatory modules by chromatin accessibility profiling

Background

The developmental gene regulatory network (GRN) that underlies skeletogenesis in sea urchins and other echinoderms is a paradigm of GRN structure, function, and evolution. This transcriptional network is deployed selectively in skeleton-forming primary mesenchyme cells (PMCs) of the early embryo. To advance our understanding of this model developmental GRN, we used genome-wide chromatin accessibility profiling to identify and characterize PMC cis-regulatory modules (CRMs).

Results

ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) analysis of purified PMCs provided a global picture of chromatin accessibility in these cells. We used both ATAC-seq and DNase-seq (DNase I hypersensitive site sequencing) to identify > 3000 sites that exhibited increased accessibility in PMCs relative to other embryonic cell lineages, and provide both computational and experimental evidence that a large fraction of these sites represent bona fide skeletogenic CRMs. Putative PMC CRMs were preferentially located near genes differentially expressed by PMCs and consensus binding sites for two key transcription factors in the PMC GRN, Alx1 and Ets1, were enriched in these CRMs. Moreover, a high proportion of candidate CRMs drove reporter gene expression specifically in PMCs in transgenic embryos. Surprisingly, we found that PMC CRMs were partially open in other embryonic lineages and exhibited hyperaccessibility as early as the 128-cell stage.

Conclusions

Our work provides a comprehensive picture of chromatin accessibility in an early embryonic cell lineage. By identifying thousands of candidate PMC CRMs, we significantly enhance the utility of the sea urchin skeletogenic network as a general model of GRN architecture and evolution. Our work also shows that differential chromatin accessibility, which has been used for the high-throughput identification of enhancers in differentiated cell types, is a powerful approach for the identification of CRMs in early embryonic cells. Lastly, we conclude that in the sea urchin embryo, CRMs that control the cell type-specific expression of effector genes are hyperaccessible several hours in advance of gene activation.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4542-z) contains supplementary material, which is available to authorized users.

HpBase: A genome database of a sea urchin, Hemicentrotus pulcherrimus

To understand the mystery of life, it is important to accumulate genomic information for various organisms because the whole genome encodes the commands for all the genes. Since the genome of Strongylocentrotus purpratus was sequenced in 2006 as the first sequenced genome in echinoderms, the genomic resources of other North American sea urchins have gradually been accumulated, but no sea urchin genomes are available in other areas, where many scientists have used the local species and reported important results. In this manuscript, we report a draft genome of the sea urchin Hemincentrotus pulcherrimus because this species has a long history as the target of developmental and cell biology in East Asia. The genome of H. pulcherrimus was assembled into 16,251 scaffold sequences with an N50 length of 143 kbp, and approximately 25,000 genes were identified in the genome. The size of the genome and the sequencing coverage were estimated to be approximately 800 Mbp and 100×, respectively. To provide these data and information of annotation, we constructed a database, HpBase (http://cell-innovation.nig.ac.jp/Hpul/). In HpBase, gene searches, genome browsing, and blast searches are available. In addition, HpBase includes the "recipes" for experiments from each lab using H. pulcherrimus. These recipes will continue to be updated according to the circumstances of individual scientists and can be powerful tools for experimental biologists and for the community. HpBase is a suitable dataset for evolutionary, developmental, and cell biologists to compare H. pulcherrimus genomic information with that of other species and to isolate gene information.