Transcriptomic analysis of the highly derived radial body plan of a sea urchin

Integrating Complex Life Cycles in Comparative Developmental Biology

The goal of comparative developmental biology is identifying mechanistic differences in embryonic development between different taxa and how these evolutionary changes have led to morphological and organizational differences in adult body plans. Much of this work has focused on direct-developing species in which the adult forms straight from the embryo and embryonic modifications have direct effects on the adult. However, most animal lineages are defined by indirect development, in which the embryo gives rise to a larval body plan and the adult forms by transformation of the larva. Historically, much of our understanding of complex life cycles is viewed through the lenses of ecology and zoology. In this review, we discuss the importance of establishing developmental rather than morphological or ecological criteria for defining developmental mode and explicitly considering the evolutionary implications of incorporating complex life cycles into broad developmental comparisons of embryos across metazoans.

A case of hermaphroditism in the gonochoristic sea urchin, Strongylocentrotus purpuratus, reveals key mechanisms of sex determination†

Sea urchins are usually gonochoristic, with all of their five gonads either testes or ovaries. Here, we report an unusual case of hermaphroditism in the purple sea urchin, Strongylocentrotus purpuratus. The hermaphrodite is self-fertile, and one of the gonads is an ovotestis; it is largely an ovary with a small segment containing fully mature sperm. Molecular analysis demonstrated that each gonad producedviable gametes, and we identified for the first time a somatic sex-specific marker in this phylum: Doublesex and mab-3 related transcription factor 1 (DMRT1). This finding also enabled us to analyze the somatic tissues of the hermaphrodite, and we found that the oral tissues (including gut) were out of register with the aboral tissues (including tube feet) enabling a genetic lineage analysis. Results from this study support a genetic basis of sex determination in sea urchins, the viability of hermaphroditism, and distinguish gonad determination from somatic tissue organization in the adult.

Transcriptomic analysis of sea star development through metamorphosis to the highly derived pentameral body plan with a focus on neural transcription factors

The Echinodermata is characterized by a secondarily evolved pentameral body plan. While the evolutionary origin of this body plan has been the subject of debate, the molecular mechanisms underlying its development are poorly understood. We assembled a de novo developmental transcriptome from the embryo through metamorphosis in the sea star Parvulastra exigua. We use the asteroid model as it represents the basal-type echinoderm body architecture. Global variation in gene expression distinguished the gastrula profile and showed that metamorphic and juvenile stages were more similar to each other than to the pre-metamorphic stages, pointing to the marked changes that occur during metamorphosis. Differential expression and gene ontology (GO) analyses revealed dynamic changes in gene expression throughout development and the transition to pentamery. Many GO terms enriched during late metamorphosis were related to neurogenesis and signalling. Neural transcription factor genes exhibited clusters with distinct expression patterns. A suite of these genes was up-regulated during metamorphosis (e.g. Pax6, Eya, Hey, NeuroD, FoxD, Mbx, and Otp). In situ hybridization showed expression of neural genes in the CNS and sensory structures. Our results provide a foundation to understand the metamorphic transition in echinoderms and the genes involved in development and evolution of pentamery.

Transcriptomic analysis of Nodal - and BMP- associated genes during development to the juvenile seastar in Parvulastra exigua (Asterinidae)

The molecular mechanisms underlying development of the pentameral body of adult echinoderms are poorly understood but are important to solve with respect to evolution of a unique body plan that contrasts with the bilateral body plan of other deuterostomes. As Nodal and BMP2/4 signalling is involved in axis formation in larvae and development of the echinoderm body plan, we used the developmental transcriptome generated for the asterinid seastar Parvulastra exigua to investigate the temporal expression patterns of Nodal and BMP2/4 genes from the embryo and across metamorphosis to the juvenile. For echinoderms, the Asteroidea represents the basal-type body architecture with a distinct (separated) ray structure. Parvulastra exigua has lecithotrophic development forming the juvenile soon after gastrulation providing ready access to the developing adult stage. We identified 39 genes associated with the Nodal and BMP2/4 network in the P. exigua developmental transcriptome. Clustering analysis of these genes resulted in 6 clusters with similar temporal expression patterns across development. A co-expression analysis revealed genes that have similar expression profiles as Nodal and BMP2/4. These results indicated genes that may have a regulatory relationship in patterning morphogenesis of the juvenile seastar. Developmental RNA-seq analyses of Parvulastra exigua show changes in Nodal and BMP2/4 signalling genes across the metamorphic transition. We provide the foundation for detailed analyses of this cascade in the evolution of the unusual pentameral echinoderm body and its deuterostome affinities.

Ocean acidification induces distinct transcriptomic responses across life history stages of the sea urchin Heliocidaris erythrogramma

Ocean acidification (OA) from seawater uptake of rising carbon dioxide emissions impairs development in marine invertebrates, particularly in calcifying species. Plasticity in gene expression is thought to mediate many of these physiological effects, but how these responses change across life history stages remains unclear. The abbreviated lecithotrophic development of the sea urchin Heliocidaris erythrogramma provides a valuable opportunity to analyse gene expression responses across a wide range of life history stages, including the benthic, post-metamorphic juvenile. We measured the transcriptional response to OA in H. erythrogramma at three stages of the life cycle (embryo, larva, and juvenile) in a controlled breeding design. The results reveal a broad range of strikingly stage-specific impacts of OA on transcription, including changes in the number and identity of affected genes; the magnitude, sign, and variance of their expression response; and the developmental trajectory of expression. The impact of OA on transcription was notably modest in relation to gene expression changes during unperturbed development and much smaller than genetic contributions from parentage. The latter result suggests that natural populations may provide an extensive genetic reservoir of resilience to OA. Taken together, these results highlight the complexity of the molecular response to OA, its substantial life history stage specificity, and the importance of contextualizing the transcriptional response to pH stress in light of normal development and standing genetic variation to better understand the capacity for marine invertebrates to adapt to OA.

Transcriptomic analysis of sea cucumber (Holothuria leucospilota) coelomocytes revealed the echinoderm cytokine response during immune challenge

Background

The sea cucumber Holothuria leucospilota belongs to echinoderm, which is evolutionally the most primitive group of deuterostomes. Sea cucumber has a cavity between its digestive tract and the body wall that is filled with fluid and suspended coelomic cells similar to blood cells. The humoral immune response of the sea cucumber is based on the secretion of various immune factors from coelomocytes into the coelomic cavity. The aim of this study is to lay out a foundation for the immune mechanisms in echinoderms and their origins in chordates by using RNA-seq.

Results

Sea cucumber primary coelomocytes were isolated from healthy H. leucospilota and incubated with lipopolysaccharide (LPS, 10 μg/ml), polyinosinic-polycytidylic acid [Poly (I:C), 10 μg/ml] and heat-inactived Vibrio harveyi (10⁷ cell/ml) for 24 h, respectively. After high-throughput mRNA sequencing on an Illumina HiSeq2500, a de novo transcriptome was assembled and the Unigenes were annotated. Thirteen differentially expressed genes (DEGs) were selected randomly from our data and subsequently verified by using RT-qPCR. The results of RT-qPCR were consistent with those of the RNA-seq (R² = 0.61). The top 10 significantly enriched signaling pathways and immune-related pathways of the common and unique DEGs were screened from the transcriptome data. Twenty-one cytokine candidate DEGs were identified, which belong to 4 cytokine families, namely, BCL/CLL, EPRF1, IL-17 and TSP/TPO. Gene expression in response to LPS dose-increased treatment (0, 10, 20 and 50 μg/ml) showed that IL-17 family cytokines were significantly upregulated after 10 μg/ml LPS challenge for 24 h.

Conclusion

A de novo transcriptome was sequenced and assembled to generate the gene expression profiling across the sea cucumber coelomocytes treated with LPS, Poly (I:C) and V. harveyi. The cytokine genes identified in DEGs could be classified into 4 cytokine families, in which the expression of IL-17 family cytokines was most significantly induced after 10 μg/ml LPS challenge for 24 h. Our findings have laid the foundation not only for the research of molecular mechanisms related to the immune response in echinoderms but also for their origins in chordates, particularly in higher vertebrates.

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.

Transcriptome analysis to characterize the genes related to gonad growth and fatty acid metabolism in the sea urchin Strongylocentrotus intermedius

BACKGROUND

Sea urchin gonads of both sexes, commonly termed "roe", are highly valued seafood delicacies, and Strongylocentrotus intermedius is considered one of the tastiest sea urchins. In order to produce high-quality gonads for consumption and clarify the mechanism of gonad growth and development of the sea urchin, more genetic information, especially at the transcriptome level, is needed.

OBJECTIVE

A more thorough understanding of sea urchin gonad growth and development in both sexes could enable regulation of these processes at several stages with the aim of suppressing gametogenesis in order to produce high-quality gonads for consumption.

METHODS

The adult sea urchins S. intermedius were cultured for 3 months, and were sampled for the gonadal transcriptome analysis which has been performed on the RNAs of three male and female adults of S. intermedius in each gonad development stage.

RESULTS

Illumina sequencing raw sequence data was deposited in the NCBI Sequence Read Archive (SRA) database (PRJNA532998). It generated 560,196,356 raw reads and 548,956,944 clean reads were acquired, which were assembled into 107,850 transcripts with 44,124 genes. Comparative analysis showed the differentially expressed genes (DEGs) from 114 to 2566. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were used to determine the functional significance of these DEGs. We have selected 9 genes related to growth and 12 genes related to fatty acid biosynthesis and metabolism in sea urchin gonads.

CONCLUSION

These data for sea urchins were intended to provide markers for gonad growth and development that can be accumulated for use in aquaculture applications.

Little evidence of adaptation potential to ocean acidification in sea urchins living in "Future Ocean" conditions at a CO2 vent

Ocean acidification (OA) can be detrimental to calcifying marine organisms, with stunting of invertebrate larval development one of the most consistent responses. Effects are usually measured by short-term, within-generation exposure, an approach that does not consider the potential for adaptation. We examined the genetic response to OA of larvae of the tropical sea urchin Echinometra sp. C. raised on coral reefs that were either influenced by CO2 vents (pH ~ 7.9, future OA condition) or nonvent control reefs (pH 8.2). We assembled a high quality de novo transcriptome of Echinometra embryos (8 hr) and pluteus larvae (48 hr) and identified 68,056 SNPs. We tested for outlier SNPs and functional enrichment in embryos and larvae raised from adults from the control or vent sites. Generally, highest F ST values in embryos were observed between sites (intrinsic adaptation, most representative of the gene pool in the spawned populations). This comparison also had the highest number of outlier loci (40). In the other comparisons, classical adaptation (comparing larvae with adults from the control transplanted to either the control or vent conditions) and reverse adaptation (larvae from the vent site returned to the vent or explanted at the control), we only observed modest numbers of outlier SNPs (6-19) and only enrichment in two functional pathways. Most of the outliers detected were silent substitutions without adaptive potential. We conclude that there is little evidence of realized adaptation potential during early development, while some potential (albeit relatively low) exists in the intrinsic gene pool after more than one generation of exposure.

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.

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.

A phylogenomic resolution of the sea urchin tree of life

BACKGROUND

Echinoidea is a clade of marine animals including sea urchins, heart urchins, sand dollars and sea biscuits. Found in benthic habitats across all latitudes, echinoids are key components of marine communities such as coral reefs and kelp forests. A little over 1000 species inhabit the oceans today, a diversity that traces its roots back at least to the Permian. Although much effort has been devoted to elucidating the echinoid tree of life using a variety of morphological data, molecular attempts have relied on only a handful of genes. Both of these approaches have had limited success at resolving the deepest nodes of the tree, and their disagreement over the positions of a number of clades remains unresolved.

RESULTS

We performed de novo sequencing and assembly of 17 transcriptomes to complement available genomic resources of sea urchins and produce the first phylogenomic analysis of the clade. Multiple methods of probabilistic inference recovered identical topologies, with virtually all nodes showing maximum support. In contrast, the coalescent-based method ASTRAL-II resolved one node differently, a result apparently driven by gene tree error induced by evolutionary rate heterogeneity. Regardless of the method employed, our phylogenetic structure deviates from the currently accepted classification of echinoids, with neither Acroechinoidea (all euechinoids except echinothurioids), nor Clypeasteroida (sand dollars and sea biscuits) being monophyletic as currently defined. We show that phylogenetic signal for novel resolutions of these lineages is strong and distributed throughout the genome, and fail to recover systematic biases as drivers of our results.

CONCLUSIONS

Our investigation substantially augments the molecular resources available for sea urchins, providing the first transcriptomes for many of its main lineages. Using this expanded genomic dataset, we resolve the position of several clades in agreement with early molecular analyses but in disagreement with morphological data. Our efforts settle multiple phylogenetic uncertainties, including the position of the enigmatic deep-sea echinothurioids and the identity of the sister clade to sand dollars. We offer a detailed assessment of evolutionary scenarios that could reconcile our findings with morphological evidence, opening up new lines of research into the development and evolutionary history of this ancient clade.

Getting a Head with Ptychodera flava Larval Regeneration

Severe injury to the central nervous system of chordates often results in permanent and irreversible mental and physical challenges. While some chordates are able to repair and/or regenerate portions of their nervous system, no chordate has been shown to be able to regenerate all regions of its central nervous system after catastrophic injury or amputation. Some hemichordates, on the other hand, are able to efficiently regenerate all neural structures, including their dorsal, hollow neural tube after complete ablation. Solitary hemichordates are marine acorn worms and a sister group to the echinoderms. The hemichordate Ptychodera flava progresses from a pelagic, feeding tornaria larva to a tripartite benthic worm with an anterior proboscis, a middle collar region, and a long posterior trunk. The adult worm regenerates all body parts when bisected in the trunk, but it was unknown whether the regeneration process was present in tornaria larvae. Now, we show that P. flava larvae are capable of robust regeneration after bisection through the sagittal, coronal, and axial planes. We also use antibody staining to show that the apical sensory organ regenerates a rich, serotonin-positive complex of cells within two weeks after amputation. Cells labeled with 5-ethynyl-2'-deoxyuridine confirm that regeneration is occurring through epimorphic processes as new cells are added at the cut site and throughout the regenerating tissue. This study verifies that P. flava larvae can be used for future functional studies aimed at identifying the genetic and morphological mechanisms controlling central nervous system regeneration in a stem deuterostome.

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.

Expression of genes and proteins of the pax-six-eya-dach network in the metamorphic sea urchin: Insights into development of the enigmatic echinoderm body plan and sensory structures

BACKGROUND

Photoreception-associated genes of the Pax-Six-Eya-Dach network (PSEDN) are deployed for many roles in addition to photoreception development. In this first study of PSEDN genes during development of the pentameral body in sea urchins, we investigated their spatial expression in Heliocidaris erythrogramma.

RESULTS

Expression of PSEDN genes in the hydrocoele of early (Dach, Eya, Six1/2) and/or late (Pax6, Six3/6) larvae, and the five hydrocoele lobes, the first morphological expression of pentamery, supports a role in body plan development. Pax6, Six1/2, and Six3/6 were localized to the primary and/or secondary podia and putative sensory/neuronal cells. Six1/2 and Six3/6 were expressed in the neuropil region in the terminal disc of the podia. Dach was localized to spines. Sequential up-regulation of gene expression as new podia and spines formed was evident. Rhabdomeric opsin and pax6 protein were localized to cells in the primary podia and spines.

CONCLUSIONS

Our results support roles for PSEDN genes in development of the pentameral body plan, contributing to our understanding of how the most unusual body plan in the Bilateria may have evolved. Development of sensory cells within the Pax-Six expression field is consistent with the role of these genes in sensory cell development in diverse species. Developmental Dynamics 247:239-249, 2018. © 2017 Wiley Periodicals, Inc.

Gaskell revisited: new insights into spinal autonomics necessitate a revised motor neuron nomenclature

Several concepts developed in the nineteenth century have formed the basis of much of our neuroanatomical teaching today. Not all of these were based on solid evidence nor have withstood the test of time. Recent evidence on the evolution and development of the autonomic nervous system, combined with molecular insights into the development and diversification of motor neurons, challenges some of the ideas held for over 100 years about the organization of autonomic motor outflow. This review provides an overview of the original ideas and quality of supporting data and contrasts this with a more accurate and in depth insight provided by studies using modern techniques. Several lines of data demonstrate that branchial motor neurons are a distinct motor neuron population within the vertebrate brainstem, from which parasympathetic visceral motor neurons of the brainstem evolved. The lack of an autonomic nervous system in jawless vertebrates implies that spinal visceral motor neurons evolved out of spinal somatic motor neurons. Consistent with the evolutionary origin of brainstem parasympathetic motor neurons out of branchial motor neurons and spinal sympathetic motor neurons out of spinal motor neurons is the recent revision of the organization of the autonomic nervous system into a cranial parasympathetic and a spinal sympathetic division (e.g., there is no sacral parasympathetic division). We propose a new nomenclature that takes all of these new insights into account and avoids the conceptual misunderstandings and incorrect interpretation of limited and technically inferior data inherent in the old nomenclature.

Nodal and BMP expression during the transition to pentamery in the sea urchin Heliocidaris erythrogramma: insights into patterning the enigmatic echinoderm body plan

BACKGROUND

The molecular mechanisms underlying the development of the unusual echinoderm pentameral body plan and their likeness to mechanisms underlying the development of the bilateral plans of other deuterostomes are of interest in tracing body plan evolution. In this first study of the spatial expression of genes associated with Nodal and BMP2/4 signalling during the transition to pentamery in sea urchins, we investigate Heliocidaris erythrogramma, a species that provides access to the developing adult rudiment within days of fertilization.

RESULTS

BMP2/4, and the putative downstream genes, Six1/2, Eya, Tbx2/3 and Msx were expressed in the earliest morphological manifestation of pentamery during development, the five hydrocoele lobes. The formation of the vestibular ectoderm, the specialized region overlying the left coelom that forms adult ectoderm, involved the expression of putative Nodal target genes Chordin, Gsc and BMP2/4 and putative BMP2/4 target genes Dlx, Msx and Tbx. The expression of Nodal, Lefty and Pitx2 in the right ectoderm, and Pitx2 in the right coelom, was as previously observed in other sea urchins.

CONCLUSION

That genes associated with Nodal and BMP2/4 signalling are expressed in the hydrocoele lobes, indicates that they have a role in the developmental transition to pentamery, contributing to our understanding of how the most unusual body plan in the Bilateria may have evolved. We suggest that the Nodal and BMP2/4 signalling cascades might have been duplicated or split during the evolution to pentamery.

Autonomy and integration in complex parasite life cycles

Complex life cycles are common in free-living and parasitic organisms alike. The adaptive decoupling hypothesis postulates that separate life cycle stages have a degree of developmental and genetic autonomy, allowing them to be independently optimized for dissimilar, competing tasks. That is, complex life cycles evolved to facilitate functional specialization. Here, I review the connections between the different stages in parasite life cycles. I first examine evolutionary connections between life stages, such as the genetic coupling of parasite performance in consecutive hosts, the interspecific correlations between traits expressed in different hosts, and the developmental and functional obstacles to stage loss. Then, I evaluate how environmental factors link life stages through carryover effects, where stressful larval conditions impact parasites even after transmission to a new host. There is evidence for both autonomy and integration across stages, so the relevant question becomes how integrated are parasite life cycles and through what mechanisms? By highlighting how genetics, development, selection and the environment can lead to interdependencies among successive life stages, I wish to promote a holistic approach to studying complex life cycle parasites and emphasize that what happens in one stage is potentially highly relevant for later stages.

De Novo Adult Transcriptomes of Two European Brittle Stars: Spotlight on Opsin-Based Photoreception

Next generation sequencing (NGS) technology allows to obtain a deeper and more complete view of transcriptomes. For non-model or emerging model marine organisms, NGS technologies offer a great opportunity for rapid access to genetic information. In this study, paired-end Illumina HiSeqTM technology has been employed to analyse transcriptomes from the arm tissues of two European brittle star species, Amphiura filiformis and Ophiopsila aranea. About 48 million Illumina reads were generated and 136,387 total unigenes were predicted from A. filiformis arm tissues. For O. aranea arm tissues, about 47 million reads were generated and 123,324 total unigenes were obtained. Twenty-four percent of the total unigenes from A. filiformis show significant matches with sequences present in reference online databases, whereas, for O. aranea, this percentage amounts to 23%. In both species, around 50% of the predicted annotated unigenes were significantly similar to transcripts from the purple sea urchin, the closest species to date that has undergone complete genome sequencing and annotation. GO, COG and KEGG analyses were performed on predicted brittle star unigenes. We focused our analyses on the phototransduction actors involved in light perception. Firstly, two new echinoderm opsins were identified in O. aranea: one rhabdomeric opsin (homologous to vertebrate melanopsin) and one RGR opsin. The RGR-opsin is supposed to be involved in retinal regeneration while the r-opsin is suspected to play a role in visual-like behaviour. Secondly, potential phototransduction actors were identified in both transcriptomes using the fly (rhabdomeric) and mammal (ciliary) classical phototransduction pathways as references. Finally, the sensitivity of O.aranea to monochromatic light was investigated to complement data available for A. filiformis. The presence of microlens-like structures at the surface of dorsal arm plate of O. aranea could potentially explain phototactic behaviour differences between the two species. The results confirm (i) the ability of these brittle stars to perceive light using opsin-based photoreception, (ii) suggest the co-occurrence of both rhabdomeric and ciliary photoreceptors, and (iii) emphasise the complexity of light perception in this echinoderm class.

Comparative Developmental Transcriptomics Reveals Rewiring of a Highly Conserved Gene Regulatory Network during a Major Life History Switch in the Sea Urchin Genus Heliocidaris

The ecologically significant shift in developmental strategy from planktotrophic (feeding) to lecithotrophic (nonfeeding) development in the sea urchin genus Heliocidaris is one of the most comprehensively studied life history transitions in any animal. Although the evolution of lecithotrophy involved substantial changes to larval development and morphology, it is not known to what extent changes in gene expression underlie the developmental differences between species, nor do we understand how these changes evolved within the context of the well-defined gene regulatory network (GRN) underlying sea urchin development. To address these questions, we used RNA-seq to measure expression dynamics across development in three species: the lecithotroph Heliocidaris erythrogramma, the closely related planktotroph H. tuberculata, and an outgroup planktotroph Lytechinus variegatus. Using well-established statistical methods, we developed a novel framework for identifying, quantifying, and polarizing evolutionary changes in gene expression profiles across the transcriptome and within the GRN. We found that major changes in gene expression profiles were more numerous during the evolution of lecithotrophy than during the persistence of planktotrophy, and that genes with derived expression profiles in the lecithotroph displayed specific characteristics as a group that are consistent with the dramatically altered developmental program in this species. Compared to the transcriptome, changes in gene expression profiles within the GRN were even more pronounced in the lecithotroph. We found evidence for conservation and likely divergence of particular GRN regulatory interactions in the lecithotroph, as well as significant changes in the expression of genes with known roles in larval skeletogenesis. We further use coexpression analysis to identify genes of unknown function that may contribute to both conserved and derived developmental traits between species. Collectively, our results indicate that distinct evolutionary processes operate on gene expression during periods of life history conservation and periods of life history divergence, and that this contrast is even more pronounced within the GRN than across the transcriptome as a whole.

De novo transcriptome of the European brittle star Amphiura filiformis pluteus larvae

BACKGROUND

In non-classical model species, Next Generation Sequencing increases the ability to analyze the expression of transcripts/genes. In this study, paired-end Illumina HiSeq sequencing technology has been employed to describe a larval transcriptome generated from 64 h post-fertilization pluteus larvae of the brittle star Amphiura filiformis. We focused our analysis on the detection of actors involved in the opsin based light perception, respectively the opsins and the phototransduction actors.

METHODS & RESULTS

In this research, about 47 million high quality reads were generated and 86,572 total unigenes were predicted after de novo assembly. Of all the larval unigenes, 18% show significant matches with reference online databases. 46% of annotated larval unigenes were significantly similar to transcripts from the purple sea urchin. COG, GO and KEGG analyses were performed on predicted unigenes. Regarding the opsin-based photoreception process, even if possible actors of ciliary and rhabdomeric phototransduction cascades were detected, no ciliary or rhabdomeric opsin was identified in these larvae. Additionally, partial non-visual RGR (retinal G protein coupled receptor) opsin mRNAs were identified,possibly indicating the presence of visual cycle reaction in early pluteus larvae. The eye morphogene Pax 6 was also identified in the pluteus transcriptome.

CONCLUSIONS

Contrary to sea-urchin larvae, brittle star larvae appear to be characterized by an absence of visual-like opsins. These RNA-seq data also provide a useful resource for the echinoderm research community and researchers with an interest in larval biology.