Gene structure in the sea urchin Strongylocentrotus purpuratus based on transcriptome analysis

Genomic signatures of exceptional longevity and negligible aging in the long-lived red sea urchin

The red sea urchin (Mesocentrotus franciscanus) is one of the Earth's longest-living animals, reported to live more than 100 years with indeterminate growth, life-long reproduction, and no increase in mortality rate with age. To understand the genetic underpinnings of longevity and negligible aging, we constructed a chromosome-level assembly of the red sea urchin genome and compared it to that of short-lived sea urchin species. Genome-wide syntenic alignments identified chromosome rearrangements that distinguish short- and long-lived species. Expanded gene families in long-lived species play a role in innate immunity, sensory nervous system, and genome stability. An integrated network of genes under positive selection in the red sea urchin was involved in genomic regulation, mRNA fidelity, protein homeostasis, and mitochondrial function. Our results implicated known longevity genes in sea urchin longevity but also revealed distinct molecular signatures that may promote long-term maintenance of tissue homeostasis, disease resistance, and negligible aging.

Parentage influence on gene expression under acidification revealed through single-embryo sequencing

The dissolution of anthropogenic carbon dioxide (CO2 ) in seawater has altered its carbonate chemistry in the process of ocean acidification (OA). OA affects the viability of marine species. In particular, calcifying organisms and their early planktonic larval stages are considered vulnerable. These organisms often utilize energy reserves for metabolism rather than growth and calcification as supported by bulk RNA-sequencing (RNA-seq) experiments. Yet, transcriptomic profiling of a bulk sample reflects the average gene expression of the population, neglecting the variations between individuals, which forms the basis for natural selection. Here, we used single-embryo RNA-seq on larval sea urchin Heliocidaris crassispina, which is a commercially and ecologically valuable species in East Asia, to document gene expression changes to OA at an individual and family level. Three paternal half-sibs groups were fertilized and exposed to 3 pH conditions (ambient pH 8.0, 7.7 and 7.4) for 12 h prior to sequencing and oxygen consumption assay. The resulting transcriptomic profile of all embryos can be distinguished into four clusters, with differences in gene expressions that govern biomineralization, cell differentiation and patterning, as well as metabolism. While these responses were influenced by pH conditions, the male identities also had an effect. Specifically, a regression model and goodness of fit tests indicated a significant interaction between sire and pH on the probability of embryo membership in different clusters of gene expression. The single-embryo RNA-seq approach is promising in climate stressor research because not only does it highlight potential impacts before phenotypic changes were observed, but it also highlights variations between individuals and lineages, thus enabling a better determination of evolutionary potential.

Proteomic analysis of digestive tract peptidases and lipases from the invasive gastropod Pomacea canaliculata

BACKGROUND

The invasive gastropod Pomacea canaliculata has received great attention in the last decades as a result of its negative impact on crops agriculture, yet knowledge of their digestive physiology remains incomplete, particularly the enzymatic breakdown of macromolecules such as proteins and lipids.

RESULTS

Discovery proteomics revealed aspartic peptidases, cysteine peptidases, serine peptidases, metallopeptidases and threonine peptidases, as well as acid and neutral lipases and phospholipases along the digestive tract of P. canaliculata. Peptides specific to peptidases (139) and lipases (14) were quantified by targeted mass spectrometry. Digestion begins in the mouth via diverse salivary peptidases (nine serine peptidases; seven cysteine peptidases, one aspartic peptidase and 22 metallopeptidases) and then continues in the oesophagus (crop) via three luminal metallopeptidases (Family M12) and six serine peptidases (Family S1). Downstream, the digestive gland provides a battery of enzymes composed of aspartic peptidase (one), cysteine peptidases (nine), serine peptidases (12) and metallopeptidases (24), including aminopeptidases, carboxypeptidases and dipeptidases). The coiled gut has M1 metallopeptidases that complete the digestion of small peptides. Lipid extracellular digestion is completed by triglyceride lipases.

CONCLUSION

From an integrative physiological and anatomical perspective, P. canaliculata shows an unexpected abundance and diversity of peptidases, which participate mainly in extracellular digestion. Moreover, the previously unknown occurrence of luminal lipases from the digestive gland is reported for the first time. Salivary and digestive glands were the main tissues involved in the synthesis and secretion of these enzymes, but plausibly the few luminally exclusive peptidases are secreted by ventrolateral pouches or epithelial unicellular glands. © 2022 The Authors. Pest Management Science published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Functional characterization of the MyD88 homologs in Strongylocentrotus purpuratus

Toll-like receptor signaling is an evolutionarily conserved pathway to induce the expression of immune mediators in response to encounters with pathogens. MyD88 is a central adapter connecting the intracellular domain of the receptors to downstream kinases. Here, we conducted a comprehensive assessment of the MyD88 family in an echinoderm, Strongylocentrotus purpuratus. Of five SpMyD88s only two closely related proteins, SpMyD88A and SpMyD88B, are functional in mammalian cell lines as their overexpression facilitates the activation of the downstream transcription factor NF-κB. This requires the presence of the endogenous mammalian MyD88s, and domain swapping indicated that the death domains of S. purpuratus MyD88 are unable to efficiently connect to the respective domains of the vertebrate IRAK kinases. This suggests that the interaction surfaces between the signaling mediators in this conserved signaling pathway are not as conserved as previously thought but were likely shaped and evolved by pathogenic selection over evolutionary timescales.

Different sea urchin RAG-like genes were domesticated to carry out different functions

The closely linked recombination activating genes (RAG1 and RAG2) in vertebrates encode the core of the RAG recombinase that mediates the V(D)J recombination of the immunoglobulin and T-cell receptor genes. RAG1 and RAG2 homologues (RAG1L and RAG2L) are present in multiple invertebrate phyla, including mollusks, nemerteans, cnidarians, and sea urchins. However, the function of the invertebrates' RAGL proteins is yet unknown. The sea urchins contain multiple RAGL genes that presumably originated in a common ancestral transposon. In this study, we demonstrated that two different RAG1L genes in the sea urchin Paracentrutus lividus (PlRAG1La and PlRAG1Lb) lost their mobility and, along with PlRAG2L, were fully domesticated to carry out different functions. We found that the examined echinoid RAGL homologues have distinct expression profiles in early developmental stages and in adult tissues. Moreover, the predicted structure of the proteins suggests that while PlRAG1La could maintain its endonuclease activity and create a heterotetramer with PlRAG2L, the PlRAG1Lb adopted a different function that does not include an interaction with DNA nor a collaboration with PlRAG2L. By characterizing the different RAG homologues in the echinoid lineage, we hope to increase the knowledge about the evolution of these genes and shed light on their domestication processes.

Active DNA demethylation of developmental cis-regulatory regions predates vertebrate origins

DNA methylation [5-methylcytosine (5mC)] is a repressive gene-regulatory mark required for vertebrate embryogenesis. Genomic 5mC is tightly regulated through the action of DNA methyltransferases, which deposit 5mC, and ten-eleven translocation (TET) enzymes, which participate in its active removal through the formation of 5-hydroxymethylcytosine (5hmC). TET enzymes are essential for mammalian gastrulation and activation of vertebrate developmental enhancers; however, to date, a clear picture of 5hmC function, abundance, and genomic distribution in nonvertebrate lineages is lacking. By using base-resolution 5mC and 5hmC quantification during sea urchin and lancelet embryogenesis, we shed light on the roles of nonvertebrate 5hmC and TET enzymes. We find that these invertebrate deuterostomes use TET enzymes for targeted demethylation of regulatory regions associated with developmental genes and show that the complement of identified 5hmC-regulated genes is conserved to vertebrates. This work demonstrates that active 5mC removal from regulatory regions is a common feature of deuterostome embryogenesis suggestive of an unexpected deep conservation of a major gene-regulatory module.

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).

Evolutionary analyses of genes in Echinodermata offer insights towards the origin of metazoan phyla

Despite recent studies discussing the evolutionary impacts of gene duplications and losses among metazoans, the genomic basis for the evolution of phyla remains enigmatic. Here, we employ phylogenomic approaches to search for orthologous genes without known functions among echinoderms, and subsequently use them to guide the identification of their homologs across other metazoans. Our final set of 14 genes was obtained via a suite of homology prediction tools, gene expression data, gene ontology, and generating the Strongylocentrotus purpuratus phylome. The gene set was subjected to selection pressure analyses, which indicated that they are highly conserved and under negative selection. Their presence across broad taxonomic depths suggests that genes required to form a phylum are ancestral to that phylum. Therefore, rather than de novo gene genesis, we posit that evolutionary forces such as selection on existing genomic elements over large timescales may drive divergence and contribute to the emergence of phyla.

An approach to quantitate maternal transcripts localized in sea urchin egg cortex using RT-qPCR with accurate normalization

The sea urchin egg cortex is a peripheral region of eggs comprising a cell membrane and adjacent cytoplasm, which contains actin and tubulin cytoskeleton, cortical granules and some proteins required for early development. Method for isolation of cortices from sea urchin eggs and early embryos was developed in 1970s. Since then, this method has been reliable tool to study protein localization and cytoskeletal organization in cortex of unfertilized eggs and embryos during first cleavages. This study was aimed to estimate the reliability of RT-qPCR to analyze levels of maternal transcripts that are localized in egg cortex. Firstly, we selected seven potential reference genes, 28S, Cycb, Ebr1, GAPDH, Hmg1, Smtnl1 and Ubb, the transcripts of which are maternally deposited in sea urchin eggs. The candidate reference genes were ranked by five different algorithms (BestKeeper, CV, ΔCt, geNorm and NormFinder) based on calculated level of stability in both eggs as well as isolated cortices. Our results showed that gene ranking differs in total RNA and mRNA samples, though Ubb is most suitable reference gene in both cases. To validate feasibility of comparative analysis of eggs and isolated egg cortices, we selected Daglb-2 as a gene of interest, which transcripts are potentially localized in cortex according to transcriptome analysis, and observed increased level of Daglb-2 in egg cortices by RT-qPCR. This suggests that proposed RNA isolation method with subsequent quantitative RT-qPCR analysis can be used to determine cortical association of transcripts in sea urchin eggs.

Developmental gene regulatory network connections predicted by machine learning from gene expression data alone

Gene regulatory network (GRN) inference can now take advantage of powerful machine learning algorithms to complement traditional experimental methods in building gene networks. However, the dynamical nature of embryonic development-representing the time-dependent interactions between thousands of transcription factors, signaling molecules, and effector genes-is one of the most challenging arenas for GRN prediction. In this work, we show that successful GRN predictions for a developmental network from gene expression data alone can be obtained with the Priors Enriched Absent Knowledge (PEAK) network inference algorithm. PEAK is a noise-robust method that models gene expression dynamics via ordinary differential equations and selects the best network based on information-theoretic criteria coupled with the machine learning algorithm Elastic Net. We test our GRN prediction methodology using two gene expression datasets for the purple sea urchin, Stronglyocentrotus purpuratus, and cross-check our results against existing GRN models that have been constructed and validated by over 30 years of experimental results. Our results find a remarkably high degree of sensitivity in identifying known gene interactions in the network (maximum 81.58%). We also generate novel predictions for interactions that have not yet been described, which provide a resource for researchers to use to further complete the sea urchin GRN. Published ChIPseq data and spatial co-expression analysis further support a subset of the top novel predictions. We conclude that GRN predictions that match known gene interactions can be produced using gene expression data alone from developmental time series experiments.

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

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

Pattern of Repetitive Element Transcription Segregate Cell Lineages during the Embryogenesis of Sea Urchin Strongylocentrotus purpuratus

Repetitive elements (REs) occupy a significant part of eukaryotic genomes and are shown to play diverse roles in genome regulation. During embryogenesis of the sea urchin, a large number of REs are expressed, but the role of these elements in the regulation of biological processes remains unknown. The aim of this study was to identify the RE expression at different stages of embryogenesis. REs occupied 44% of genomic DNA of Strongylocentrotus purpuratus. The most prevalent among these elements were the unknown elements-in total, they contributed 78.5% of REs (35% in total genome occupancy). It was revealed that the transcription pattern of genes and REs changes significantly during gastrulation. Using the de novo transcriptome assembly, we showed that the expression of RE is independent of its copy number in the genome. We also identified copies that are expressed. Only active RE copies were used for mapping and quantification of RE expression in the single-cell RNA sequencing data. REs expression was observed in all cell lineages and they were detected as population markers. Moreover, the primary mesenchyme cell (PMC) line had the greatest diversity of REs among the markers. Our data suggest a role for RE in the organization of developmental domains during the sea urchin embryogenesis at the single-cell resolution level.

Dynamically expressed single ELAV/Hu orthologue elavl2 of bees is required for learning and memory

Changes in gene expression are a hallmark of learning and memory consolidation. Little is known about how alternative mRNA processing, particularly abundant in neuron-specific genes, contributes to these processes. Prototype RNA binding proteins of the neuronally expressed ELAV/Hu family are candidates for roles in learning and memory, but their capacity to cross-regulate and take over each other's functions complicate substantiation of such links. Honey bees Apis mellifera have only one elav/Hu family gene elavl2, that has functionally diversified by increasing alternative splicing including an evolutionary conserved microexon. RNAi knockdown demonstrates that ELAVL2 is required for learning and memory in bees. ELAVL2 is dynamically expressed with altered alternative splicing and subcellular localization in mushroom bodies, but not in other brain regions. Expression and alternative splicing of elavl2 change during memory consolidation illustrating an alternative mRNA processing program as part of a local gene expression response underlying memory consolidation.

A nomenclature for echinoderm genes

Echinoderm embryos and larvae are prominent experimental model systems for studying developmental mechanisms. High-quality, assembled, annotated genome sequences are now available for several echinoderm species, including representatives from most classes. The increased availability of these data necessitates the development of a nomenclature that assigns universally interpretable gene symbols to echinoderm genes to facilitate cross-species comparisons of gene functions, both within echinoderms and across other phyla. This paper describes the implementation of an improved set of echinoderm gene nomenclature guidelines that both communicates meaningful orthology information in protein-coding gene symbols and names and establishes continuity with nomenclatures developed for major vertebrate model organisms, including humans. Differences between the echinoderm gene nomenclature guidelines and vertebrate guidelines are examined and explained. This nomenclature incorporates novel solutions to allow for several types of orthologous relationships, including the single echinoderm genes with multiple vertebrate co-orthologs that result from whole-genome-duplication events. The current version of the Echinoderm Gene Nomenclature Guidelines can be found at https://www.echinobase.org/gene/static/geneNomenclature.jsp

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

Global patterns of enhancer activity during sea urchin embryogenesis assessed by eRNA profiling

We used capped analysis of gene expression with sequencing (CAGE-seq) to profile eRNA expression and enhancer activity during embryogenesis of a model echinoderm: the sea urchin, Strongylocentrotus purpuratus. We identified more than 18,000 enhancers that were active in mature oocytes and developing embryos and documented a burst of enhancer activation during cleavage and early blastula stages. We found that a large fraction (73.8%) of all enhancers active during the first 48 h of embryogenesis were hyperaccessible no later than the 128-cell stage and possibly even earlier. Most enhancers were located near gene bodies, and temporal patterns of eRNA expression tended to parallel those of nearby genes. Furthermore, enhancers near lineage-specific genes contained signatures of inputs from developmental gene regulatory networks deployed in those lineages. A large fraction (60%) of sea urchin enhancers previously shown to be active in transgenic reporter assays was associated with eRNA expression. Moreover, a large fraction (50%) of a representative subset of enhancers identified by eRNA profiling drove tissue-specific gene expression in isolation when tested by reporter assays. Our findings provide an atlas of developmental enhancers in a model sea urchin and support the utility of eRNA profiling as a tool for enhancer discovery and regulatory biology. The data generated in this study are available at Echinobase, the public database of information related to echinoderm genomics.

Guardian of the Genome: An Alternative RAG/Transib Co-Evolution Hypothesis for the Origin of V(D)J Recombination

The appearance of adaptive immunity in jawed vertebrates is termed the immunological 'Big Bang' because of the short evolutionary time over which it developed. Underlying it is the recombination activating gene (RAG)-based V(D)J recombination system, which initiates the sequence diversification of the immunoglobulins and lymphocyte antigen receptors. It was convincingly argued that the RAG1 and RAG2 genes originated from a single transposon. The current dogma postulates that the V(D)J recombination system was established by the split of a primordial vertebrate immune receptor gene into V and J segments by a RAG1/2 transposon, in parallel with the domestication of the same transposable element in a separate genomic locus as the RAG recombinase. Here, based on a new interpretation of previously published data, we propose an alternative evolutionary hypothesis suggesting that two different elements, a RAG1/2 transposase and a Transib transposon invader with RSS-like terminal inverted repeats, co-evolved to work together, resulting in a functional recombination process. This hypothesis offers an alternative understanding of the acquisition of recombinase function by RAGs and the origin of the V(D)J system.

MeDAS: a Metazoan Developmental Alternative Splicing database

Alternative splicing is widespread throughout eukaryotic genomes and greatly increases transcriptomic diversity. Many alternative isoforms have functional roles in developmental processes and are precisely temporally regulated. To facilitate the study of alternative splicing in a developmental context, we created MeDAS, a Metazoan Developmental Alternative Splicing database. MeDAS is an added-value resource that re-analyses publicly archived RNA-seq libraries to provide quantitative data on alternative splicing events as they vary across the time course of development. It has broad temporal and taxonomic scope and is intended to assist the user in identifying trends in alternative splicing throughout development. To create MeDAS, we re-analysed a curated set of 2232 Illumina polyA+ RNA-seq libraries that chart detailed time courses of embryonic and post-natal development across 18 species with a taxonomic range spanning the major metazoan lineages from Caenorhabditis elegans to human. MeDAS is freely available at https://das.chenlulab.com both as raw data tables and as an interactive browser allowing searches by species, tissue, or genomic feature (gene, transcript or exon ID and sequence). Results will provide details on alternative splicing events identified for the queried feature and can be visualised at the gene-, transcript- and exon-level as time courses of expression and inclusion levels, respectively.

Identification and functional characterization of piwi1 gene in sea cucumber, Apostichopus japonicas

The sea cucumber (Apostichopus japonicus) is an economically important mariculture species in Asia. However, the genetic breeding of sea cucumbers is difficult because the sexes cannot be identified by appearance. Therefore, studies on sex-related genes are helpful in revealing the mechanisms of sex determination and differentiation in sea cucumbers. P-element induced wimpy testis (piwi) is a germ cell marker involved in gametogenesis in vertebrates; however, the expression pattern and function during gametogenesis remain unclear in sea cucumbers. In this study, we identified a piwi homolog gene in A. japonicus (Ajpiwi1) and investigated its expression pattern, and function. Ajpiwi1 is a maternal factor and is ubiquitously expressed in adult tissues, including the ovary and testis. Ajpiwi1 expression is strong in early oocytes, spermatocytes, and spermatogonia; weak in mature oocytes; and undetected in spermatids and intra-gonadal somatic cells. The knockdown of Ajpiwi1 by RNA interference (RNAi) led to the downregulation of other conserved sex-related genes such as dmrt1, foxl2, and germ cell-less. Therefore, Ajpiwi1 might play a critical role during gametogenesis in A. japonicus. This study creates new possibilities for studying sex-related gene functions in the sea cucumber and builds a gene function research platform based on RNAi for the first time.

Phylogenetic Analyses of Glycosyl Hydrolase Family 6 Genes in Tunicates: Possible Horizontal Transfer

Horizontal gene transfer (HGT) is the movement of genetic material between different species. Although HGT is less frequent in eukaryotes than in bacteria, several instances of HGT have apparently shaped animal evolution. One well-known example is the tunicate cellulose synthase gene, CesA, in which a gene, probably transferred from bacteria, greatly impacted tunicate evolution. A Glycosyl Hydrolase Family 6 (GH6) hydrolase-like domain exists at the C-terminus of tunicate CesA, but not in cellulose synthases of other organisms. The recent discovery of another GH6 hydrolase-like gene (GH6-1) in tunicate genomes further raises the question of how tunicates acquired GH6. To examine the probable origin of these genes, we analyzed the phylogenetic relationship of GH6 proteins in tunicates and other organisms. Our analyses show that tunicate GH6s, the GH6-1 gene, and the GH6 part of the CesA gene, form two independent, monophyletic gene groups. We also compared their sequence signatures and exon splice sites. All tunicate species examined have shared splice sites in GH6-containing genes, implying ancient intron acquisitions. It is likely that the tunicate CesA and GH6-1 genes existed in the common ancestor of all extant tunicates.

Genomic insights of body plan transitions from bilateral to pentameral symmetry in Echinoderms

Echinoderms are an exceptional group of bilaterians that develop pentameral adult symmetry from a bilaterally symmetric larva. However, the genetic basis in evolution and development of this unique transformation remains to be clarified. Here we report newly sequenced genomes, developmental transcriptomes, and proteomes of diverse echinoderms including the green sea urchin (L. variegatus), a sea cucumber (A. japonicus), and with particular emphasis on a sister group of the earliest-diverged echinoderms, the feather star (A. japonica). We learned that the last common ancestor of echinoderms retained a well-organized Hox cluster reminiscent of the hemichordate, and had gene sets involved in endoskeleton development. Further, unlike in other animal groups, the most conserved developmental stages were not at the body plan establishing phase, and genes normally involved in bilaterality appear to function in pentameric axis development. These results enhance our understanding of the divergence of protostomes and deuterostomes almost 500 Mya.

Li et al. investigate the evolution and genetic basis of the adult pentameral body plan in echinoderms using genomic, transcriptomic, and proteomic data. They determine that the last common ancestor of echinoderms contained an organized Hox cluster and endoskeleton genes, and suggest that cooption of bilateral development genes was involved in evolution of the pentameric body plan.

Unique age-related transcriptional signature in the nervous system of the long-lived red sea urchin Mesocentrotus franciscanus

The red sea urchin, Mesocentrotus franciscanus, is one the earth’s longest-lived animals, reported to live more than 100 years with indeterminate growth, life-long reproduction and no increase in mortality rate with age. To gain insight into mechanisms associated with longevity and negligible senescence, age-related transcriptional profiles were examined in tissues of the red sea urchin. Genome-wide transcriptional profiling using RNA-Seq revealed few age-related changes in gene expression in muscle and esophagus tissue. In contrast, radial nerve showed an unexpected level of complexity with the expression of 3,370 genes significantly altered more than two-fold with age, including genes involved in nerve function, signaling, metabolism, transcriptional regulation and chromatin modification. There was an age-related upregulation in expression of genes involved in synaptogenesis, axonogenesis and neuroprotection suggesting preservation of neuronal processes with age. There was also an upregulation in expression of positive regulators and key components of the AMPK pathway, autophagy, proteasome function, and the unfolded protein response. This unique age-related gene expression profile in the red sea urchin nervous system may play a role in mitigating the detrimental effects of aging in this long-lived animal.

Maternal control of early patterning in sea urchin embryos

Sea urchin development has been studied extensively for more than a century and considered regulative since the first experimental evidence. Further investigations have repeatedly supported this standpoint by revealing the presence of inductive mechanisms that alter cell fate decisions at early cleavage stages and flexibility of development in response to environmental conditions. Some features indicate that sea urchin development is not completely regulative, but actually includes determinative events. In 16-cell embryos, mesomeres and macromeres represent multipotency, while the cell fate of most vegetal micromeres is restricted. It is known that the mature sea urchin eggs are polarized by the asymmetrical distribution of some maternal mRNAs and proteins. Spatially-distributed maternal factors are necessary for the orientation of the primary animal-vegetal axis, which is established by both maternal and zygotic mechanisms later in development. The secondary dorsal-ventral axis is conditionally specified later in development. Dorsal-ventral polarity is very liable during the early cleavages, though more recent data argue that its direction may be oriented by maternal asymmetry. In this review, we focus on the role of maternal factors in initial embryonic patterning during the first cleavages of sea urchin embryos before activation of the embryonic genome.

A Peak of H3T3 Phosphorylation Occurs in Synchrony with Mitosis in Sea Urchin Early Embryos

The sea urchin embryo provides a valuable system to analyse the molecular mechanisms orchestrating cell cycle progression and mitosis in a developmental context. However, although it is known that the regulation of histone activity by post-translational modification plays an important role during cell division, the dynamics and the impact of these modifications have not been characterised in detail in a developing embryo. Using different immuno-detection techniques, we show that the levels of Histone 3 phosphorylation at Threonine 3 oscillate in synchrony with mitosis in Sphaerechinus granularis early embryos. We present, in addition, the results of a pharmacological study aimed at analysing the role of this key histone post-translational modification during sea urchin early development.

Ion channels and signaling pathways used in the fast polyspermy block

Fertilization of an egg by multiple sperms, polyspermy, is lethal to most sexually reproducing species. To combat the entry of additional sperm into already fertilized eggs, organisms have developed various polyspermy blocks. One such barrier, the fast polyspermy block, uses a fertilization-activated depolarization of the egg membrane to electrically inhibit supernumerary sperm from entering the egg. The fast block is commonly used by eggs of oviparous animals with external fertilization. In this review, we discuss the history of the fast block discovery, as well as general features shared by all organisms that use this polyspermy block. Given the diversity of habitats of external fertilizers, the fine details of the fast block-signaling pathways differ drastically between species, including the identity of the depolarizing ions. We highlight the known molecular mediators of these signaling pathways in amphibians and echinoderms, with a fine focus on ion channels that signal these fertilization-evoked depolarizations. We also discuss the investigation for a fast polyspermy block in mammals and teleost fish, and we outline potential fast block triggers. Since the first electrical recordings made on eggs in the 1950s, the fields of developmental biology and electrophysiology have substantially matured, and yet we are only now beginning to discern the intricate molecular mechanisms regulating the fast block to polyspermy.

Enhancement of oxidative stress contributes to increased pathogenicity of the invasive pine wood nematode

Reactive oxygen species (ROS) play important roles in defence response of host plants versus pathogens. While generation and detoxification of ROS is well understood, how varied ability of different isolates of pathogens to overcome host ROS, or ROS contribution to a particular isolate's pathogenicity, remains largely unexplored. Here, we report that transcriptional regulation of the ROS pathway, in combination with the insulin pathway, increases the pathogenicity of invasive species Bursaphelenchus xylophilus. The results showed a positive correlation between fecundity and pathogenicity of different nematode isolates. The virulent isolates from introduced populations in Japan, China and Europe had significantly higher fecundity than native avirulent isolates from the USA. Increased expression of Mn-SOD and reduced expression of catalase/ GPX-5 and H₂O₂ accumulation during invasion are associated with virulent strains. Additional H₂O₂ could improve fecundity of Bu. xylophilus. Furthermore, depletion of Mn-SOD decreased fecundity and virulence of Bu. xylophilus, while the insulin pathway is significantly affected. Thus, we propose that destructive pathogenicity of Bu. xylophilus to pines is partly owing to upregulated fecundity modulated by the insulin pathway in association with the ROS pathway and further enhanced by H₂O₂ oxidative stress. These findings provide a better understanding of pathogenic mechanisms in plant-pathogen interactions and adaptive evolution of invasive species. This article is part of the theme issue 'Biotic signalling sheds light on smart pest management'.

Evolutionary History of GLIS Genes Illuminates Their Roles in Cell Reprograming and Ciliogenesis

The GLIS family transcription factors, GLIS1 and GLIS3, potentiate generation of induced pluripotent stem cells (iPSCs). In contrast, another GLIS family member, GLIS2, suppresses cell reprograming. To understand how these disparate roles arose, we examined evolutionary origins and genomic organization of GLIS genes. Comprehensive phylogenetic analysis shows that GLIS1 and GLIS3 originated during vertebrate whole genome duplication, whereas GLIS2 is a sister group to the GLIS1/3 and GLI families. This result is consistent with their opposing functions in cell reprograming. Glis1 evolved faster than Glis3, losing many protein-interacting motifs. This suggests that Glis1 acquired new functions under weakened evolutionary constraints. In fact, GLIS1 induces induced pluripotent stem cells more strongly. Transcriptomic data from various animal embryos demonstrate that glis1 is maternally expressed in some tetrapods, whereas vertebrate glis3 and invertebrate glis1/3 genes are rarely expressed in oocytes, suggesting that vertebrate (or tetrapod) Glis1 acquired a new expression domain and function as a maternal factor. Furthermore, comparative genomic analysis reveals that glis1/3 is part of a bilaterian-specific gene cluster, together with rfx3, ndc1, hspb11, and lrrc42. Because known functions of these genes are related to cilia formation and function, the last common ancestor of bilaterians may have acquired this cluster by shuffling gene order to establish more sophisticated epithelial tissues involving cilia. This evolutionary study highlights the significance of GLIS1/3 for cell reprograming, development, and diseases in ciliated organs such as lung, kidney, and pancreas.

The developmental transcriptome for Lytechinus variegatus exhibits temporally punctuated gene expression changes

Embryonic development is arguably the most complex process an organism undergoes during its lifetime, and understanding this complexity is best approached with a systems-level perspective. The sea urchin has become a highly valuable model organism for understanding developmental specification, morphogenesis, and evolution. As a non-chordate deuterostome, the sea urchin occupies an important evolutionary niche between protostomes and vertebrates. Lytechinus variegatus (Lv) is an Atlantic species that has been well studied, and which has provided important insights into signal transduction, patterning, and morphogenetic changes during embryonic and larval development. The Pacific species, Strongylocentrotus purpuratus (Sp), is another well-studied sea urchin, particularly for gene regulatory networks (GRNs) and cis-regulatory analyses. A well-annotated genome and transcriptome for Sp are available, but similar resources have not been developed for Lv. Here, we provide an analysis of the Lv transcriptome at 11 timepoints during embryonic and larval development. Temporal analysis suggests that the gene regulatory networks that underlie specification are well-conserved among sea urchin species. We show that the major transitions in variation of embryonic transcription divide the developmental time series into four distinct, temporally sequential phases. Our work shows that sea urchin development occurs via sequential intervals of relatively stable gene expression states that are punctuated by abrupt transitions.

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.

The evolutionary road to invertebrate thyroid hormone signaling: Perspectives for endocrine disruption processes

Thyroid hormones (THs) are the only iodine-containing hormones that play fundamental roles in chordates and non-chordates. The chemical nature, mode of action and the synthesis of THs are well established in mammals and other vertebrates. Although thyroid-like hormones have been detected in protostomes and non-chordate deuterostomes, TH signaling is poorly understood as compared to vertebrates, particularly in protostomes. Therefore, the central objective of this article is to review TH system components and TH-induced effects in non-vertebrate chordates, non-chordate deuterostomes and protostomes based on available genomes and functional information. To accomplish this task, we integrate here the available knowledge on the THs signaling across non-vertebrate chordates, non-chordate deuterostomes and protostomes by considering studies encompassing TH system components and physiological actions of THs. We also address the possible interactions of thyroid disrupting chemicals and their effects in protostomes and non-chordate deuterostomes. Finally, the perspectives on current and future challenges are discussed.

Genome-wide identification of binding sites and gene targets of Alx1, a pivotal regulator of echinoderm skeletogenesis

Alx1 is a conserved regulator of skeletogenesis in echinoderms and evolutionary changes in Alx1 sequence and expression have played a pivotal role in modifying programs of skeletogenesis within the phylum. Alx1 regulates a large suite of effector genes that control the morphogenetic behaviors and biomineral-forming activities of skeletogenic cells. To better understand the gene regulatory control of skeletogenesis by Alx1, we used genome-wide ChIP-seq to identify Alx1-binding sites and direct gene targets. Our analysis revealed that many terminal differentiation genes receive direct transcriptional inputs from Alx1. In addition, we found that intermediate transcription factors previously shown to be downstream of Alx1 all receive direct inputs from Alx1. Thus, Alx1 appears to regulate effector genes by indirect, as well as direct, mechanisms. We tested 23 high-confidence ChIP-seq peaks using GFP reporters and identified 18 active cis-regulatory modules (CRMs); this represents a high success rate for CRM discovery. Detailed analysis of a representative CRM confirmed that a conserved, palindromic Alx1-binding site was essential for expression. Our work significantly advances our understanding of the gene regulatory circuitry that controls skeletogenesis in sea urchins and provides a framework for evolutionary studies.

Trans-Species Polymorphism in Mitochondrial Genome of Camarodont Sea Urchins

Mitochondrial (mt) genomes of the sea urchins Strongylocentrotus intermedius and Mesocentrotus nudus demonstrate the identical patterns of intraspecific length variability of the ND6 gene, consisting of 489 bp (S variant) and 498 bp (L variant), respectively. For both species, the ND6 length difference is due to the 488A>G substitution, which changes the stop codon TAG in S variant for a tryptophan codon TGG in L variant and elongates the corresponding ND6 protein by three additional amino acids, Trp-Leu-Trp. The phylogenetic analysis based on mt genomes of sea urchins and related echinoderm groups from GenBank has shown the S and L ND6 variants as shared among the camarodont sea urchins; the rest of the echinoderms demonstrate the S variant only. The data suggest that the ND6 488A>G substitution can be the first example of the trans-species polymorphism in sea urchins, persisting at least since the time of the Odontophora diversification at the Eocene/Oligocene boundary (approximately 34 million years ago), which was characterized by an abrupt climate change and significant global ocean cooling. Alternative hypotheses, including the convergent RNA editing and/or codon reassignment, are not supported by direct comparisons of the ND6 gene sequences with the corresponding transcripts using the basic local alignment search tool (BLAST) of full sea urchin transcriptomes.

The Axial Organ and the Pharynx Are Sites of Hematopoiesis in the Sea Urchin

Background: The location of coelomocyte proliferation in adult sea urchins is unknown and speculations since the early 1800s have been based on microanatomy and tracer uptake studies. In adult sea urchins (Strongylocentrotus purpuratus) with down-regulated immune systems, coelomocyte numbers increase in response to immune challenge, and whether some or all of these cells are newly proliferated is not known. The gene regulatory network that encodes transcription factors that control hematopoiesis in embryonic and larval sea urchins has not been investigated in adults. Hence, to identify the hematopoietic tissue in adult sea urchins, cell proliferation, expression of phagocyte specific genes, and expression of genes encoding transcription factors that function in the conserved regulatory network that controls hematopoiesis in embryonic and larval sea urchins were investigated for several tissues. Results: Cell proliferation was induced in adult sea urchins either by immune challenge through injection of heat-killed Vibrio diazotrophicus or by cell depletion through aspiration of coelomic fluid. In response to either of these stimuli, newly proliferated coelomocytes constitute only about 10% of the cells in the coelomic fluid. In tissues, newly proliferated cells and cells that express SpTransformer proteins (formerly Sp185/333) that are markers for phagocytes are present in the axial organ, gonad, pharynx, esophagus, and gut with no differences among tissues. The expression level of genes encoding transcription factors that regulate hematopoiesis show that both the axial organ and the pharynx have elevated expression compared to coelomocytes, esophagus, gut, and gonad. Similarly, an RNAseq dataset shows similar results for the axial organ and pharynx, but also suggests that the axial organ may be a site for removal and recycling of cells in the coelomic cavity. Conclusions: Results presented here are consistent with previous speculations that the axial organ may be a site of coelomocyte proliferation and that it may also be a center for cellular removal and recycling. A second site, the pharynx, may also have hematopoietic activity, a tissue that has been assumed to function only as part of the intestinal tract.

Culture of and experiments with sea urchin embryo primary mesenchyme cells

Skeletogenesis in the sea urchin embryo gives rise to a pair of intricate endoskeletal spicules. Deposition of these skeletal elements in the early larva is the outcome of a morphogenetic program that begins with maternal inputs in the early zygote and results in the specification of the large micromere-primary mesenchyme cell (PMC) lineage. PMCs are of considerable interest as a model system, not only to dissect the mechanism of specific developmental processes, but also to investigate their evolution and the unrivaled level of control over the formation of a graded, mechanically robust, yet single crystalline biomineral. The ability to study gene regulatory circuits, cellular behavior, signaling pathways, and molecular players involved in biomineralization is significantly boosted by the high level of autonomy of PMCs. In fact, in the presence of horse serum, micromeres differentiate into PMCs and produce spicules in vitro, separated from the embryonic milieu. PMC culture eliminates indirect effects that can complicate the interpretation of experiments in vivo, offers superior spatiotemporal control, enables PMC-specific readouts, and is compatible with most imaging and characterization techniques. In this chapter, we provide an updated protocol, based on the pioneering work by Okazaki and Wilt, for the isolation of micromeres and subsequent culture of PMCs, as well as protocols for fixation and staining for fluorescent microscopy, preparation of cell cultures for electron microscopy, and the isolation of RNA.

Genomic resources for the study of echinoderm development and evolution

Echinoderms are important research models for a wide range of biological questions. In particular, echinoderm embryos are exemplary models for dissecting the molecular and cellular processes that drive development and testing how these processes can be modified through evolution to produce the extensive morphological diversity observed in the phylum. Modern attempts to characterize these processes depend on some level of genomic analysis; from querying annotated gene sets to functional genomics experiments to identify candidate cis-regulatory sequences. Given how essential these data have become, it is important that researchers using available datasets or performing their own genome-scale experiments understand the nature and limitations of echinoderm genomic analyses. In this chapter we highlight the current state of echinoderm genomic data and provide methodological considerations for common approaches, including analysis of transcriptome and functional genomics datasets.

Two more Posterior Hox genes and Hox cluster dispersal in echinoderms

BACKGROUND

Hox genes are key elements in patterning animal development. They are renowned for their, often, clustered organisation in the genome, with supposed mechanistic links between the organisation of the genes and their expression. The widespread distribution and comparable functions of Hox genes across the animals has led to them being a major study system for comparing the molecular bases for construction and divergence of animal morphologies. Echinoderms (including sea urchins, sea stars, sea cucumbers, feather stars and brittle stars) possess one of the most unusual body plans in the animal kingdom with pronounced pentameral symmetry in the adults. Consequently, much interest has focused on their development, evolution and the role of the Hox genes in these processes. In this context, the organisation of echinoderm Hox gene clusters is distinctive. Within the classificatory system of Duboule, echinoderms constitute one of the clearest examples of Disorganized (D) clusters (i.e. intact clusters but with a gene order or orientation rearranged relative to the ancestral state).

RESULTS

Here we describe two Hox genes (Hox11/13d and e) that have been overlooked in most previous work and have not been considered in reconstructions of echinoderm Hox complements and cluster organisation. The two genes are related to Posterior Hox genes and are present in all classes of echinoderm. Importantly, they do not reside in the Hox cluster of any species for which genomic linkage data is available.

CONCLUSION

Incorporating the two neglected Posterior Hox genes into assessments of echinoderm Hox gene complements and organisation shows that these animals in fact have Split (S) Hox clusters rather than simply Disorganized (D) clusters within the Duboule classification scheme. This then has implications for how these genes are likely regulated, with them no longer covered by any potential long-range Hox cluster-wide, or multigenic sub-cluster, regulatory mechanisms.

Phylogenetic Distribution of CMP-Neu5Ac Hydroxylase (CMAH), the Enzyme Synthetizing the Proinflammatory Human Xenoantigen Neu5Gc

The enzyme CMP-N-acetylneuraminic acid hydroxylase (CMAH) is responsible for the synthesis of N-glycolylneuraminic acid (Neu5Gc), a sialic acid present on the cell surface proteins of most deuterostomes. The CMAH gene is thought to be present in most deuterostomes, but it has been inactivated in a number of lineages, including humans. The inability of humans to synthesize Neu5Gc has had several evolutionary and biomedical implications. Remarkably, Neu5Gc is a xenoantigen for humans, and consumption of Neu5Gc-containing foods, such as red meats, may promote inflammation, arthritis, and cancer. Likewise, xenotransplantation of organs producing Neu5Gc can result in inflammation and organ rejection. Therefore, knowing what animal species contain a functional CMAH gene, and are thus capable of endogenous Neu5Gc synthesis, has potentially far-reaching implications. In addition to humans, other lineages are known, or suspected, to have lost CMAH; however, to date reports of absent and pseudogenic CMAH genes are restricted to a handful of species. Here, we analyze all available genomic data for nondeuterostomes, and 322 deuterostome genomes, to ascertain the phylogenetic distribution of CMAH. Among nondeuterostomes, we found CMAH homologs in two green algae and a few prokaryotes. Within deuterostomes, putatively functional CMAH homologs are present in 184 of the studied genomes, and a total of 31 independent gene losses/pseudogenization events were inferred. Our work produces a list of animals inferred to be free from endogenous Neu5Gc based on the absence of CMAH homologs and are thus potential candidates for human consumption, xenotransplantation research, and model organisms for investigation of human diseases.

From genome to anatomy: The architecture and evolution of the skeletogenic gene regulatory network of sea urchins and other echinoderms

The skeletogenic gene regulatory network (GRN) of sea urchins and other echinoderms is one of the most intensively studied transcriptional networks in any developing organism. As such, it serves as a preeminent model of GRN architecture and evolution. This review summarizes our current understanding of this developmental network. We describe in detail the most comprehensive model of the skeletogenic GRN, one developed for the euechinoid sea urchin Strongylocentrotus purpuratus, including its initial deployment by maternal inputs, its elaboration and stabilization through regulatory gene interactions, and its control of downstream effector genes that directly drive skeletal morphogenesis. We highlight recent comparative studies that have leveraged the euechinoid GRN model to examine the evolution of skeletogenic programs in diverse echinoderms, studies that have revealed both conserved and divergent features of skeletogenesis within the phylum. Last, we summarize the major insights that have emerged from analysis of the structure and evolution of the echinoderm skeletogenic GRN and identify key, unresolved questions as a guide for future work.

Comparative Studies of Gene Expression Kinetics: Methodologies and Insights on Development and Evolution

Across the animal kingdom, embryos of closely related species show high morphological similarity despite genetic and environmental distances. Deciphering the molecular mechanisms that underlie morphological conservation and those that support embryonic adaptation are keys to understand developmental robustness and evolution. Comparative studies of developmental gene regulatory networks can track the genetic changes that lead to evolutionary novelties. However, these studies are limited to a relatively small set of genes and demand extensive experimental efforts. An alternative approach enabled by next-generation sequencing, is to compare the expression kinetic of large sets of genes between different species. The advantages of these comparisons are that they can be done relatively easily, for any species and they provide information of all expressed genes. The challenge in these experiments is to compare the kinetic profiles of thousands of genes between species that develop in different rates. Here we review recent comparative studies that tackled the challenges of accurate staging and large-scale analyses using different computational approaches. These studies reveal how correct temporal scaling exposes the striking conservation of developmental gene expression between morphologically similar species. Different clustering approaches are used to address various comparative questions and identify the conservation and divergence of large gene sets. We discuss the unexpected contribution of housekeeping genes to the interspecies correlations and how this contribution distorts the hourglass pattern generated by developmental genes. Overall, we demonstrate how comparative studies of gene expression kinetics can provide novel insights into the developmental constraints and plasticity that shape animal body plans.

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

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.

SpTransformer proteins from the purple sea urchin opsonize bacteria, augment phagocytosis, and retard bacterial growth

The purple sea urchin, Strongylocentrotus purpuratus, has a complex and robust immune system that is mediated by a number of multi-gene families including the SpTransformer (SpTrf) gene family (formerly Sp185/333). In response to immune challenge from bacteria and various pathogen-associated molecular patterns, the SpTrf genes are up-regulated in sea urchin phagocytes and express a diverse array of SpTrf proteins. We show here that SpTrf proteins from coelomocytes and isolated by nickel affinity (cNi-SpTrf) bind to Gram-positive and Gram-negative bacteria and to Baker's yeast, Saccharomyces cerevisiae, with saturable kinetics and specificity. cNi-SpTrf opsonization of the marine bacteria, Vibrio diazotrophicus, augments phagocytosis, however, opsonization by the recombinant protein, rSpTrf-E1, does not. Binding by cNi-SpTrf proteins retards growth rates significantly for several species of bacteria. SpTrf proteins, previously thought to be strictly membrane-associated, are secreted from phagocytes in short term cultures and bind V. diazotrophicus that are located both outside of and within phagocytes. Our results demonstrate anti-microbial activities of native SpTrf proteins and suggest variable functions among different SpTrf isoforms. Multiple isoforms may act synergistically to detect a wide array of pathogens and provide flexible and efficient host immunity.

A SLC4 family bicarbonate transporter is critical for intracellular pH regulation and biomineralization in sea urchin embryos

Efficient pH regulation is a fundamental requisite of all calcifying systems in animals and plants but with the underlying pH regulatory mechanisms remaining largely unknown. Using the sea urchin larva, this work identified the SLC4 HCO₃^- transporter family member SpSlc4a10 to be critically involved in the formation of an elaborate calcitic endoskeleton. SpSlc4a10 is specifically expressed by calcifying primary mesenchyme cells with peak expression during de novo formation of the skeleton. Knock-down of SpSlc4a10 led to pH regulatory defects accompanied by decreased calcification rates and skeleton deformations. Reductions in seawater pH, resembling ocean acidification scenarios, led to an increase in SpSlc4a10 expression suggesting a compensatory mechanism in place to maintain calcification rates. We propose a first pH regulatory and HCO₃^- concentrating mechanism that is fundamentally linked to the biological precipitation of CaCO₃. This knowledge will help understanding biomineralization strategies in animals and their interaction with a changing environment.

Many marine organisms such as mussels, sea urchins or corals, have skeletons and shells, which – due to their beautiful colors and shapes – are often desirable collector pieces. These structures are made from calcium and carbonate ions that react to form calcium carbonate crystals in a process known as biomineralization.

In sea urchin larvae, for example, the skeleton is built by so-called primary mesenchyme cells, which work similar to the bone forming cells in mammals. These mesenchyme cells use calcium from the sea water, which travels to the site where the shell starts to form. About half of the carbonate comes from carbon dioxide that the animals make as they breathe, but it is not known how the other half gets to the site of biomineralization.

Producing a skeleton generates acid, and marine animals need to be able to regulate their pH levels, as too acidic environments can dissolve the calcium carbonate and threatening to destroy the developing shell. How cells accumulate enough carbonate to make their shells, and how they cope with acidity is still poorly understood.

Here, Hu et al. address this problem by studying purple sea urchin larvae, revealing that they use ion transporters to gather bicarbonate from seawater. These structures are part of a group of bicarbonate transporters known as the ‘SLC4 transporter family’, which sit across the membrane of the mesenchyme cells and move the bicarbonate ions along. As the sea urchin larvae develop, the levels of the transporter protein start to rise in mesenchyme cells, peaking around the time they are producing the skeleton.

Stopping the production of the transporter hindered the larvae from building normal skeletons and also made their cells more acidic. It turns out that bicarbonate is more than a skeleton ingredient – it also helps to buffer the acid made in the process. Bicarbonate ions can combine with acidic molecules to form water and carbon dioxide. Bicarbonate pumped in from the sea neutralises the acidic molecules made during calcium carbonate formation, which helps to stabilize pH levels.

When the acidity of the water was experimentally increased, it prompted the sea urchins to produce more of the SLC4 transporters, revealing that they may have another role to play. Their acid-neutralizing capability helped the animals to cope with changes in their environment. Taking on more bicarbonate could therefore help to compensate for rising acidity, allowing skeleton production to carry on as normal.

This last finding is important in the context of ocean acidification. As the amount of carbon dioxide in the atmosphere increases, more of the gas dissolves in the sea. The chemical reactions that follow make the water more acidic and decreases the pH levels of the sea. Understanding how animals make their skeletons and shells, and manage acid, could reveal how they will cope as the environment changes in the future.

Engineering Pichia pastoris for Efficient Production of a Novel Bifunctional Strongylocentrotus purpuratus Invertebrate-Type Lysozyme

Lysozymes are known as ubiquitously distributed immune effectors with hydrolytic activity against peptidoglycan, the major bacterial cell wall polymer, to trigger cell lysis. In the present study, the full-length cDNA sequence of a novel sea urchin Strongylocentrotus purpuratus invertebrate-type lysozyme (sp-iLys) was synthesized according to the codon usage bias of Pichia pastoris and was cloned into a constitutive expression plasmid pPIC9K. The resulting plasmid, pPIC9K-sp-iLys, was integrated into the genome of P. pastoris strain GS115. The bioactive recombinant sp-iLys was successfully secreted into the culture broth by positive transformants. The highest lytic activity of 960 U/mL of culture supernatant was reached in fed-batch fermentation. Using chitin affinity chromatography and gel-filtration chromatography, recombinant sp-iLys was produced with a yield of 94.5 mg/L and purity of > 99%. Recombinant sp-iLys reached its peak lytic activity of 8560 U/mg at pH 6.0 and 30 °C and showed antimicrobial activities against Gram-negative bacteria (Vibrio vulnificus, Vibrio parahemolyticus, and Aeromonas hydrophila) and Gram-positive bacteria (Staphylococcus aureus and Bacillus subtilis). In addition, recombinant sp-iLys displayed isopeptidase activity which reached the peak at pH 7.5 and 37 °C with the presence of 0.05 M Na⁺. In conclusion, this report describes the heterologous expression of recombinant sp-iLys in P. pastoris on a preparative-scale, which possesses lytic activity and isopeptidase activity. This suggests that sp-iLys might play an important role in the innate immunity of S. purpuratus.

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.

Developmental transcriptomics of the brittle star Amphiura filiformis reveals gene regulatory network rewiring in echinoderm larval skeleton evolution

Background

Amongst the echinoderms the class Ophiuroidea is of particular interest for its phylogenetic position, ecological importance and developmental and regenerative biology. However, compared to other echinoderms, notably echinoids (sea urchins), relatively little is known about developmental changes in gene expression in ophiuroids. To address this issue, we have generated and assembled a large RNAseq data set of four key stages of development in the brittle star Amphiura filiformis and a de novo reference transcriptome of comparable quality to that of a model echinoderm—the sea urchin Strongylocentrotus purpuratus. Furthermore, we provide access to the new data via a web interface: http://www.echinonet.eu/shiny/Amphiura_filiformis/.

Results

We have identified highly conserved genes associated with the development of a biomineralised skeleton. We also identify important class-specific characters, including the independent duplication of the msp130 class of genes in different echinoderm classes and the unique occurrence of spicule matrix (sm) genes in echinoids. Using a new quantification pipeline for our de novo transcriptome, validated with other methodologies, we find major differences between brittle stars and sea urchins in the temporal expression of many transcription factor genes. This divergence in developmental regulatory states is more evident in early stages of development when cell specification begins, rather than when cells initiate differentiation.

Conclusions

Our findings indicate that there has been a high degree of gene regulatory network rewiring and clade-specific gene duplication, supporting the hypothesis of a convergent evolution of larval skeleton development in echinoderms.

Electronic supplementary material

The online version of this article (10.1186/s13059-018-1402-8) contains supplementary material, which is available to authorized users.

EchinoBase: Tools for Echinoderm Genome Analyses

The echinoderms are a phylum of invertebrate deuterostome animals that constitute important research models for a number of biological disciplines. EchinoBase ( www.echinobase.org ) is an echinoderm-specific genome database and web information system that provides a platform for the interrogation and exploration of echinoderm genomic data. This chapter outlines the datasets available on EchinoBase; from assembled genomes and genome annotations, to spatial and quantitative expression data, as well as functional genomics datasets. We also highlight the bioinformatic tools available on the website to facilitate rapid inquiries using these data (genome browsers, precompiled BLAST databases, etc.), and suggest optimized strategies for performing these inquiries. We conclude with a perspective on how one could integrate various genomic resources to predict putative noncoding regulatory regions. The available datasets and analyses they permit provide the basic components required for developing an understanding of how echinoderm genomes are regulated, especially during early development, and provides a platform for comparative genomic inquiries among species in this phylum.

Evolution of ligands, receptors and metabolizing enzymes of thyroid signaling

Thyroid hormones (THs) play important roles in vertebrates such as the control of the metabolism, development and seasonality. Given the pleiotropic effects of thyroid disorders (developmental delay, mood disorder, tachycardia, etc), THs signaling is highly investigated, specially using mammalian models. In addition, the critical role of TH in controlling frog metamorphosis has led to the use of Xenopus as another prominent model to study THs action. Nevertheless, animals regarded as non-model species can also improve our understanding of THs signaling. For instance, studies in amphioxus highlighted the role of Triac as a bona fide thyroid hormone receptor (TR) ligand. In this review, we discuss our current understanding of the THs signaling in the different taxa forming the metazoans (multicellular animals) group. We mainly focus on three actors of the THs signaling: the ligand, the receptor and the deiodinases, enzymes playing a critical role in THs metabolism. By doing so, we also pinpoint many key questions that remain unanswered. How can THs accelerate metamorphosis in tunicates and echinoderms while their TRs have not been yet demonstrated as functional THs receptors in these species? Do THs have a biological effect in insects and cnidarians even though they do not have any TR? What is the basic function of THs in invertebrate protostomia? These questions can appear disconnected from pharmacological issues and human applications, but the investigation of THs signaling at the metazoans scale can greatly improve our understanding of this major endocrinological pathway.