Developmental transcriptome of Aplysia californica

Molecular Cloning, Characterization, and Function of Insulin-Related Peptide 1 (IRP1) in the Haliotis discus hanna

Abalone is a popular mollusk in the marine aquaculture industry of China. However, existing challenges, like slow growth, individual miniaturization, and the absence of abundant abalone, have emerged as significant obstacles impeding its long-term progress in aquaculture. Studies have demonstrated that insulin-related peptide (IRP) is a crucial factor in the growth of marine organisms. However, limited studies have been conducted on IRP in abalone. This study indicated that the hdh-MIRP1 open reading frame (ORF) was composed of 456 base pairs, which encoded 151 amino acids. Based on the gene expression and immunofluorescence analyses, the cerebral ganglion of Haliotis discus hannai (H. discus hannai) was the primary site of hdh-MIRP1 mRNA expression. Moreover, hdh-MIRP1 expression was observed to be higher in the larger group than in the smaller group abalones. Only single nucleotide polymorphism (SNP) was related to their growth characteristics. However, approximately 82 proteins that may interact with hdh-MIRP1 were identified. The functional enrichment analysis of the 82 genes indicated that hdh-MIRP1 may be involved in the regulation of glucose metabolism and the process of growth. This study established a benchwork for further investigating the role of IRP in the growth of abalone.

Thyroid hormone-induced cell death in sea urchin metamorphic development

Thyroid hormones (THs) are important regulators of development, metabolism and homeostasis in metazoans. Specifically, they have been shown to regulate the metamorphic transitions of vertebrates and invertebrates alike. Indirectly developing sea urchin larvae accelerate the formation of juvenile structures in response to thyroxine (T4) treatment, while reducing their larval arm length. The mechanisms underlying larval arm reduction are unknown and we hypothesized that programmed cell death (PCD) is linked to this process. To test this hypothesis, we measured larval arm retraction in response to different THs (T4, T3, rT3, Tetrac) and assessed cell death in larvae using three different methods (TUNEL, YO-PRO-1 and caspase-3 activity) in the sea urchin Strongylocentrotus purpuratus. We also compared the extent of PCD in response to TH treatment before and after the invagination of the larval ectoderm, which marks the initiation of juvenile development in larval sea urchin species. We found that T4 treatment results in the strongest reduction of larval arms but detected a significant increase of PCD in response to T4, T3 and Tetrac in post-ingression but not pre-ingression larvae. As post-ingression larvae have initiated metamorphic development and therefore allocate resources to both larval and the juvenile structures, these results provide evidence that THs regulate larval development differentially via PCD. PCD in combination with cell proliferation likely has a key function in sea urchin development.

Regulation of gene expression during ontogeny of physiological function in the brackishwater amphipod Gammarus chevreuxi

Embryonic development is a complex process involving the co-ordinated onset and integration of multiple morphological features and physiological functions. While the molecular basis of morphological development in embryos is relatively well known for traditional model species, the molecular underpinning of the development of physiological functions is not. Here, we used global gene expression profiling to investigate the transcriptional changes associated with the development of morphological and physiological function in the amphipod crustacean Gammarus chevreuxi. We compared the transcriptomes at three timepoints during the latter half of development, characterised by different stages of the development of heart form and function: 10 days post fertilisation (dpf, Early: no heart structure visible), 15 dpf (Middle: heart present but not fully functional), and 18 dpf (Late: regular heartbeat). Gene expression profiles differed markedly between developmental stages, likely representing a change in the activity of different processes throughout the latter period of G. chevreuxi embryonic development. Differentially expressed genes belonged to one of three distinct clusters based on their expression patterns across development. One of these clusters, which included key genes relating to cardiac contractile machinery and calcium handling, displayed a pattern of sequential up-regulation throughout the developmental period studied. Further analyses of these transcripts could reveal genes that may influence the onset of a regular heartbeat. We also identified morphological and physiological processes that may occur alongside heart development, such as development of digestive caeca and the cuticle. Elucidating the mechanisms underpinning morphological and physiological development of non-model organisms will support improved understanding of conserved mechanisms, addressing the current phylogenetic gap between relatively well known model species.

Developmental series of gene expression clarifies maternal mRNA provisioning and maternal-to-zygotic transition in a reef-building coral

BACKGROUND

Maternal mRNA provisioning of oocytes regulates early embryogenesis. Maternal transcripts are degraded as zygotic genome activation (ZGA) intensifies, a phenomenon known as the maternal-to-zygotic transition (MZT). Here, we examine gene expression over nine developmental stages in the Pacific rice coral, Montipora capitata, from eggs and embryos at 1, 4, 9, 14, 22, and 36 h-post-fertilization (hpf), as well as swimming larvae (9d), and adult colonies.

RESULTS

Weighted Gene Coexpression Network Analysis revealed four expression peaks, identifying the maternal complement, two waves of the MZT, and adult expression. Gene ontology enrichment revealed maternal mRNAs are dominated by cell division, methylation, biosynthesis, metabolism, and protein/RNA processing and transport functions. The first MZT wave occurs from ~4-14 hpf and is enriched in terms related to biosynthesis, methylation, cell division, and transcription. In contrast, functional enrichment in the second MZT wave, or ZGA, from 22 hpf-9dpf, includes ion/peptide transport and cell signaling. Finally, adult expression is enriched for functions related to signaling, metabolism, and ion/peptide transport. Our proposed MZT timing is further supported by expression of enzymes involved in zygotic transcriptional repression (Kaiso) and activation (Sox2), which peak at 14 hpf and 22 hpf, respectively. Further, DNA methylation writing (DNMT3a) and removing (TET1) enzymes peak and remain stable past ~4 hpf, suggesting that methylome programming occurs before 4 hpf.

CONCLUSIONS

Our high-resolution insight into the coral maternal mRNA and MZT provides essential baseline information to understand parental carryover effects and the sensitivity of developmental success under increasing environmental stress.

Hormonal Regulation of Programmed Cell Death in Sea Urchin Metamorphosis

Programmed cell death (PCD) has been identified as a key process in the metamorphic transition of indirectly developing organisms such as frogs and insects. Many marine invertebrate species with indirect development and biphasic life cycles face the challenge of completing the metamorphic transition of the larval body into a juvenile when they settle into the benthic habitat. Some key characteristics stand out during this transition in comparison to frogs and insects: (1) the transition is often remarkably fast and (2) the larval body is largely abandoned and few structures transition into the juvenile stage. In sea urchins, a group with a drastic and fast metamorphosis, development and destruction of the larval body is regulated by endocrine signals. Here we provide a brief review of the basic regulatory mechanisms of PCD in animals. We then narrow our discussion to metamorphosis with a specific emphasis on sea urchins with indirect life histories and discuss the function of thyroid hormones and histamine in larval development, metamorphosis and settlement of the sea urchin Strongylocentrotus purpuratus. We were able to annotate the large majority of PCD related genes in the sea urchin S. purpuratus and ongoing studies on sea urchin metamorphosis will shed light on the regulatory architecture underlying this dramatic life history transition. While we find overwhelming evidence for hormonal regulation of PCD in animals, especially in the context of metamorphosis, the mechanisms in many marine invertebrate groups with indirect life histories requires more work. Hence, we propose that studies of PCD in animals requires functional studies in whole organisms rather than isolated cells. We predict that future work, targeting a broader array of organisms will not only help to reveal important new functions of PCD but provide a fundamentally new perspective on its use in a diversity of taxonomic, developmental, and ecological contexts.

Temporally balanced selection during development of larval Pacific oysters (Crassostrea gigas) inherently preserves genetic diversity within offspring

Balancing selection is one of the mechanisms which has been proposed to explain the maintenance of genetic diversity in species across generations. For species with large populations and complex life histories, however, heterogeneous selection pressures may create a scenario in which the net effects of selection are balanced across developmental stages. With replicated cultures and a pooled sequencing approach, we show that genotype-dependent mortality in larvae of the Pacific oyster (Crassostrea gigas) is largely temporally dynamic and inconsistently in favour of a single genotype or allelic variant at each locus. Overall, the patterns of genetic change we observe to be taking place are more complex than what would be expected under classical examples of additive or dominant genetic interactions. They are also not easily explained by our current understanding of the effects of genetic load. Collectively, temporally heterogeneous selection pressures across different larval developmental stages may act to maintain genetic diversity, while also inherently sheltering genetic load within oyster populations.

Transcriptome analysis provides genome annotation and expression profiles in the central nervous system of Lymnaea stagnalis at different ages

BACKGROUND

The pond snail, Lymnaea stagnalis (L. stagnalis), has served as a valuable model organism for neurobiology studies due to its simple and easily accessible central nervous system (CNS). L. stagnalis has been widely used to study neuronal networks and recently gained popularity for study of aging and neurodegenerative diseases. However, previous transcriptome studies of L. stagnalis CNS have been exclusively carried out on adult L. stagnalis only. As part of our ongoing effort studying L. stagnalis neuronal growth and connectivity at various developmental stages, we provide the first age-specific transcriptome analysis and gene annotation of young (3 months), adult (6 months), and old (18 months) L. stagnalis CNS.

RESULTS

Using the above three age cohorts, our study generated 55-69 millions of 150 bp paired-end RNA sequencing reads using the Illumina NovaSeq 6000 platform. Of these reads, ~ 74% were successfully mapped to the reference genome of L. stagnalis. Our reference-based transcriptome assembly predicted 42,478 gene loci, of which 37,661 genes encode coding sequences (CDS) of at least 100 codons. In addition, we provide gene annotations using Blast2GO and functional annotations using Pfam for ~ 95% of these sequences, contributing to the largest number of annotated genes in L. stagnalis CNS so far. Moreover, among 242 previously cloned L. stagnalis genes, we were able to match ~ 87% of them in our transcriptome assembly, indicating a high percentage of gene coverage. The expressional differences for innexins, FMRFamide, and molluscan insulin peptide genes were validated by real-time qPCR. Lastly, our transcriptomic analyses revealed distinct, age-specific gene clusters, differentially expressed genes, and enriched pathways in young, adult, and old CNS. More specifically, our data show significant changes in expression of critical genes involved in transcription factors, metabolisms (e.g. cytochrome P450), extracellular matrix constituent, and signaling receptor and transduction (e.g. receptors for acetylcholine, N-Methyl-D-aspartic acid, and serotonin), as well as stress- and disease-related genes in young compared to either adult or old snails.

CONCLUSIONS

Together, these datasets are the largest and most updated L. stagnalis CNS transcriptomes, which will serve as a resource for future molecular studies and functional annotation of transcripts and genes in L. stagnalis.

Development and Interrogation of a Transcriptomic Resource for the Giant Triton Snail (Charonia tritonis)

Gastropod molluscs are among the most abundant species that inhabit coral reef ecosystems. Many are specialist predators, along with the giant triton snail Charonia tritonis (Linnaeus, 1758) whose diet consists of Acanthaster planci (crown-of-thorns starfish), a corallivore known to consume enormous quantities of reef-building coral. C. tritonis are considered vulnerable due to overexploitation, and a decline in their populations is believed to have contributed to recurring A. planci population outbreaks. Aquaculture is considered one approach that could help restore natural populations of C. tritonis and mitigate coral loss; however, numerous questions remain unanswered regarding their life cycle, including the molecular factors that regulate their reproduction and development. In this study, we have established a reference C. tritonis transcriptome derived from developmental stages (embryo and veliger) and adult tissues. This was used to identify genes associated with cell signalling, such as neuropeptides and G protein-coupled receptors (GPCRs), involved in endocrine and olfactory signalling. A comparison of developmental stages showed that several neuropeptide precursors are exclusively expressed in post-hatch veligers and functional analysis found that FFamide stimulated a significant (20.3%) increase in larval heart rate. GPCRs unique to veligers, and a diversity of rhodopsin-like GPCRs located within adult cephalic tentacles, all represent candidate olfactory receptors. In addition, the cytochrome P450 superfamily, which participates in the biosynthesis and degradation of steroid hormones and lipids, was also found to be expanded with at least 91 genes annotated, mostly in gill tissue. These findings further progress our understanding of C. tritonis with possible application in developing aquaculture methods.

ATP signaling in the integrative neural center of Aplysia californica

ATP and its ionotropic P2X receptors are components of the most ancient signaling system. However, little is known about the distribution and function of purinergic transmission in invertebrates. Here, we cloned, expressed, and pharmacologically characterized the P2X receptors in the sea slug Aplysia californica—a prominent neuroscience model. AcP2X receptors were successfully expressed in Xenopus oocytes and displayed activation by ATP with two-phased kinetics and Na⁺-dependence. Pharmacologically, they were different from other P2X receptors. The ATP analog, Bz-ATP, was a less effective agonist than ATP, and PPADS was a more potent inhibitor of the AcP2X receptors than the suramin. AcP2X were uniquely expressed within the cerebral F-cluster, the multifunctional integrative neurosecretory center. AcP2X receptors were also detected in the chemosensory structures and the early cleavage stages. Therefore, in molluscs, rapid ATP-dependent signaling can be implicated both in development and diverse homeostatic functions. Furthermore, this study illuminates novel cellular and systemic features of P2X-type ligand-gated ion channels for deciphering the evolution of neurotransmitters.

Global Analysis of Gene Expression Profiles Provides Novel Insights into the Development and Evolution of the Large Crustacean Eriocheir sinensis

Chinese mitten crab (Eriocheir sinensis) is an important aquaculture species in Crustacea. Functional analysis, although essential, has been hindered due to the lack of sufficient genomic or transcriptomic resources. In this study, transcriptome sequencing was conducted on 59 samples representing diverse developmental stages (fertilized eggs, zoea, megalopa, three sub-stages of larvae, juvenile crabs, and adult crabs) and different tissues (eyestalk, hepatopancreas, and muscle from juvenile crabs, and eyestalk, hepatopancreas, muscle, heart, stomach, gill, thoracic ganglia, intestine, ovary, and testis from adult crabs) of E. sinensis. A comprehensive reference transcriptome was assembled, including 19,023 protein-coding genes. Hierarchical clustering based on 128 differentially expressed cuticle-related genes revealed two distinct expression patterns during the early larval developmental stages, demonstrating the distinct roles of these genes in “crab-like” cuticle formation during metamorphosis and cuticle calcification after molting. Phylogenetic analysis of 1406 one-to-one orthologous gene families identified from seven arthropod species and Caenorhabditis elegans strongly supported the hypothesis that Malacostraca and Branchiopoda do not form a monophyletic group. Furthermore, Branchiopoda is more phylogenetically closely related to Hexapoda, and the clade of Hexapoda and Branchiopoda and the clade of Malacostraca belong to the Pancrustacea. This study offers a high-quality transcriptome resource for E. sinensis and demonstrates the evolutionary relationships of major arthropod groups. The differentially expressed genes identified in this study facilitate further investigation of the cuticle-related gene expression networks which are likely associated with “crab-like” cuticle formation during metamorphosis and cuticle calcification after molting.

The transcriptome analysis of the whole-body of the gastropod mollusk Limax flavus and screening of putative antimicrobial peptide and protein genes

The gastropod mollusk Limax flavus, one of the most widespread pests in China, is used to treat infectious diseases in traditional Chinese medicine. However, little genomic information is available for this non-model species. In this study, the whole-body transcriptome of L. flavus was sequenced using next generation sequencing technology. A total of 6.81 Gb clean reads were obtained, which were assembled into 150,766 transcripts with 132,206 annotated unigenes. Functionally classification assigned 30,542 unigenes to 56 Gene Ontology terms, 16,745 unigenes were divided into 26 euKaryotic Ortholog Groups of proteins categories, and 13,854 unigenes were assigned to 230 Kyoto Encyclopedia of Genes and Genomes pathways. Furthermore, we identified 17,251 simple sequence repeats and several kinds of antimicrobial peptide and protein (AMPs) genes. The transcriptome data of L. flavus will provide a valuable genomic resource for further studies on this species, and the AMPs identified in L. flavus will support its medical potential.

De novo transcriptome of the whole-body of the gastropod mollusk Philomycus bilineatus, a pest with medical potential in China

Philomycus bilineatus is a highly common gastropod mollusk pest in China and is also utilized to treat infectious diseases. However, no genomic resources are available for this non-model species. In the present study, the transcriptomic analysis of P. bilineatus was completed. After sequencing using the next generation sequencing technology, 9.11 Gb of clean reads were obtained, which led to the assembly and annotation of 145,523 transcripts and 125,690 unigenes. Unigenes were functionally classified using Gene Ontology (GO), euKaryotic Ortholog Groups of proteins (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG). A total of 27,554 unigenes were assigned into 55 GO terms, 13,989 unigenes were differentiated into 26 KOG categories, and 16,368 unigenes were assigned to 229 KEGG pathways. Furthermore, 16,614 simple sequence repeats (SSRs), 38 olfactory genes, and 40 antimicrobial peptide/protein genes were identified. The transcriptome profile of P. bilineatus will provide a valuable genomic resource for further study, will promote the development of new pest management strategies through interference of chemosensory communication, and will support potential medicinal uses of this species.

Transcriptional profiles of early stage red sea urchins (Mesocentrotus franciscanus) reveal differential regulation of gene expression across development

The red sea urchin, Mesocentrotus franciscanus, is an ecologically important kelp forest species that also serves as a valuable fisheries resource. In this study, we have assembled and annotated a developmental transcriptome for M. franciscanus that represents eggs and six stages of early development (8- to 16-cell, morula, hatched blastula, early gastrula, prism and early pluteus). Characterization of the transcriptome revealed distinct patterns of gene expression that corresponded to major developmental and morphological processes. In addition, the period during which maternally-controlled transcription was terminated and the zygotic genome was activated, the maternal-to-zygotic transition (MZT), was found to begin during early cleavage and persist through the hatched blastula stage, an observation that is similar to the timing of the MZT in other sea urchin species. The presented developmental transcriptome will serve as a useful resource for investigating, in both an ecological and fisheries context, how the early developmental stages of this species respond to environmental stressors.

Transcriptome analysis of air-breathing land slug, Incilaria fruhstorferi reveals functional insights into growth, immunity, and reproduction

Background

Incilaria (= Meghimatium) fruhstorferi is an air-breathing land slug found in restricted habitats of Japan, Taiwan and selected provinces of South Korea (Jeju, Chuncheon, Busan, and Deokjeokdo). The species is on a decline due to depletion of forest cover, predation by natural enemies, and collection. To facilitate the conservation of the species, it is important to decide on a number of traits related to growth, immunity and reproduction addressing fitness advantage of the species.

Results

The visceral mass transcriptome of I. fruhstorferi was enabled using the Illumina HiSeq 4000 sequencing platform. According to BUSCO (Benchmarking Universal Single-Copy Orthologs) method, the transcriptome was considered complete with 91.8% of ortholog genes present (Single: 70.7%; Duplicated: 21.1%). A total of 96.79% of the raw read sequences were processed as clean reads. TransDecoder identified 197,271 contigs that contained candidate-coding regions. Of a total of 50,230 unigenes, 34,470 (68.62% of the total unigenes) annotated to homologous proteins in the Protostome database (PANM-DB). The GO term and KEGG pathway analysis indicated genes involved in metabolism, phosphatidylinositol signalling system, aminobenzoate degradation, and T-cell receptor signalling pathway. Many genes associated with molluscan innate immunity were categorized under pathogen recognition receptor, TLR signalling pathway, MyD88 dependent pathway, endogenous ligands, immune effectors, antimicrobial peptides, apoptosis, and adaptation-related. The reproduction-associated unigenes showed homology to protein fem-1, spermatogenesis-associated protein, sperm associated antigen, and testis expressed sequences, among others. In addition, we identified key growth-related genes categorized under somatotrophic axis, muscle growth, chitinases and collagens. A total of 4822 Simple Sequence Repeats (SSRs) were also identified from the unigene sequences of I. fruhstorferi.

Conclusions

This is the first available genomic information for non-model land slug, I. fruhstorferi focusing on genes related to growth, immunity, and reproduction, with additional focus on microsatellites and repeating elements. The transcriptome provides access to greater number of traits of unknown relevance in the species that could be exploited for in-depth analyses of evolutionary plasticity and making informed choices during conservation planning. This would be appropriate for understanding the dynamics of the species on a priority basis considering the ecological, health, and social benefits.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5526-3) contains supplementary material, which is available to authorized users.

Workshop report: Cryopreservation of aquatic biomedical models

The genetic resources of aquatic biomedical model organisms are the products of millions of years of evolution, decades of scientific development, and hundreds of millions of dollars of research funding investment. Genetic resources (e.g., specific alleles, transgenes, or combinations) of each model organism can be considered a form of scientific wealth that can be accumulated and exchanged, typically in the form of live animals or germplasm. Large-scale maintenance of live aquatic organisms that carry these genetic resources is inefficient, costly, and risky. In situ maintenance may be substantially enhanced and backed up by combining cryopreserved germplasm repositories and genetic information systems with live animal culture. Unfortunately, cryopreservation has not advanced much beyond the status of an exploratory research for most aquatic species, lacks widespread application, and methods for successful cryopreservation remain poorly defined. For most aquatic species biological materials other than sperm or somatic cells are not comprehensively banked to represent and preserve a broad range of genetic diversity for each species. Therefore, new approaches and standardization are needed for repository-level application to ensure reproducible recovery of cryopreserved materials. Additionally, development of new technologies is needed to address preservation of novel biological materials, such as eggs and embryos of aquatic species. To address these goals, the Office of Research Infrastructure Programs (ORIP) of the National Institutes of Health (NIH) hosted the Cryopreservation of Aquatic Biomedical Models Workshop on January 7 to 8, 2017, in conjunction with the 8th Aquatic Animal Models of Human Disease Conference in Birmingham, Alabama. The goals of the workshop were to assess the status of germplasm cryopreservation in various biomedical aquatic models and allow representatives of the scientific community to develop and prioritize a consensus of specific actionable recommendations that will move the field of cryopreservation of aquatic resources forward. This workshop included sessions devoted to new approaches for cryopreservation of aquatic species, discussion of current efforts and approaches in preservation of aquatic model germplasm, consideration of needs for standardization of methods to support reproducibility, and enhancement of repository development by establishment of scalable high-throughput technologies. The following three broad recommendations were forwarded from workshop attendees: 1: Establish a comprehensive, centralized unit ("hub") to programmatically develop training for and documentation of cryopreservation methods for aquatic model systems. This would include development of species-specific protocols and approaches, outreach programs, community development and standardization, freezing services and training of the next generation of experts in aquatic cryopreservation. 2: Provide mechanisms to support innovative technical advancements that will increase the reliability, reproducibility, simplicity, throughput, and efficiency of the cryopreservation process, including vitrification and pipelines for sperm, oocytes, eggs, embryos, larvae, stem cells, and somatic cells of all aquatic species. This recommendation encompasses basic cryopreservation knowledge and engineering technology, such as microfluidics and automated processing technologies. 3: Implement mechanisms that allow the various aquatic model stock centers to increase their planning, personnel, ability to secure genetic resources and to promote interaction within an integrated, comprehensive repository network for aquatic model species repositories.

The evolution of ependymin-related proteins

Background

Ependymins were originally defined as fish-specific secreted glycoproteins involved in central nervous system plasticity and memory formation. Subsequent research revealed that these proteins represent a fish-specific lineage of a larger ependymin-related protein family (EPDRs). EPDRs have now been identified in a number of bilaterian animals and have been implicated in diverse non-neural functions. The recent discoveries of putative EPDRs in unicellular holozoans and an expanded EPDR family with potential roles in conspecific communication in crown-of-thorns starfish suggest that the distribution and diversity of EPDRs is significantly broader than currently understood.

Results

We undertook a systematic survey to determine the distribution and evolution of EPDRs in eukaryotes. In addition to Bilateria, EPDR genes were identified in Cnidaria, Placozoa, Porifera, Choanoflagellatea, Filasterea, Apusozoa, Amoebozoa, Charophyta and Percolozoa, and tentatively in Cercozoa and the orphan group Malawimonadidae. EPDRs appear to be absent from prokaryotes and many eukaryote groups including ecdysozoans, fungi, stramenopiles, alveolates, haptistans and cryptistans. The EPDR family can be divided into two major clades and has undergone lineage-specific expansions in a number of metazoan lineages, including in poriferans, molluscs and cephalochordates. Variation in a core set of conserved residues in EPDRs reveals the presence of three distinct protein types; however, 3D modelling predicts overall protein structures to be similar.

Conclusions

Our results reveal an early eukaryotic origin of the EPDR gene family and a dynamic pattern of gene duplication and gene loss in animals. This research provides a phylogenetic framework for the analysis of the functional evolution of this gene family.

Electronic supplementary material

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

The developmental transcriptome atlas of the spoon worm Urechis unicinctus (Echiurida: Annelida)

Background

Echiurida is one of the most intriguing major subgroups of annelida because, unlike most other annelids, echiurids lack metameric body segmentation as adults. For this reason, transcriptome analyses from various developmental stages of echiurid species can be of substantial value for understanding precise expression levels and the complex regulatory networks during early and larval development.

Results

A total of 914 million raw RNA-Seq reads were produced from 14 developmental stages of Urechis unicinctus and were de novo assembled into contigs spanning 63,928,225 bp with an N50 length of 2700 bp. The resulting comprehensive transcriptome database of the early developmental stages of U. unicinctus consists of 20,305 representative functional protein-coding transcripts. Approximately 66% of unigenes were assigned to superphylum-level taxa, including Lophotrochozoa (40%). The completeness of the transcriptome assembly was assessed using benchmarking universal single-copy orthologs; 75.7% of the single-copy orthologs were presented in our transcriptome database. We observed 3 distinct patterns of global transcriptome profiles from 14 developmental stages and identified 12,705 genes that showed dynamic regulation patterns during the differentiation and maturation of U. unicinctus cells.

Conclusions

We present the first large-scale developmental transcriptome dataset of U. unicinctus and provide a general overview of the dynamics of global gene expression changes during its early developmental stages. The analysis of time-course gene expression data is a first step toward understanding the complex developmental gene regulatory networks in U. unicinctus and will furnish a valuable resource for analyzing the functions of gene repertoires in various developmental phases.

Cloning and molecular characterization of two ferritins from red abalone Haliotis rufescens and their expressions in response to bacterial challenge at juvenile and adult life stages

Ferritins are ubiquitous proteins with a pivotal role in iron storage and homeostasis, and in host defense responses during infection by pathogens in several organisms, including mollusks. In this study, we characterized two ferritin homologues in the red abalone Haliotis rufescens, a species of economic importance for Chile, USA and Mexico. Two ferritin subunits (Hrfer1 and Hrfer2) were cloned. Hrfer1 cDNA is an 807 bp clone containing a 516 bp open reading frame (ORF) that corresponds to a novel ferritin subunit in H. rufescens. Hrfer2 cDNA is an 868 bp clone containing a 516 bp ORF that corresponds to a previously reported ferritin subunit, but in this study 5'- and 3'-UTR sequences were additionally found. We detected a putative Iron Responsive Element (IRE) in the 5'-UTR sequence, suggesting a posttranscriptional regulation of Hrfer2 translation by iron. The deduced protein sequences of both cDNAs possessed the motifs and domains required in functional ferritin subunits. Expression patterns of both ferritins in different tissues, during different developmental stages, and in response to bacterial (Vibrio splendidus) exposure were examined. Both Hrfer1 and Hrfer2 are most expressed in digestive gland and gonad. Hrfer1 mRNA levels increased about 34-fold along with larval developmental process, attaining the highest level in the creeping post-larvae. Exogenous feeding is initiated at the creeping larva stage; thus, the increase of Hrfer1 may suggest and immunity-related role upon exposure to bacteria. Highest Hrfer2 expression levels were detected at trochophore stage; which may be related with early shell formation. Upon challenge with, the bacteria an early mild induction of Hrfer2 (2 h post-challenge), followed by a stronger induction of Hrfer1 at 15 h post-challenge, was observed in haemocytes from adult abalones. While maximal upregulation of both genes in the whole individual occurred at 24 h post-challenge, in juveniles. A significant increase in ferritin protein levels from 6 h to 24 h post-challenge was also detected. Our results suggest an involvement of Hrfer1 and Hrfer2, and of ferritin proteins in the immune response of H. rufescens to bacterial infection.

Description and initial characterization of metatranscriptomic nidovirus-like genomes from the proposed new family Abyssoviridae, and from a sister group to the Coronavirinae, the proposed genus Alphaletovirus

Transcriptomics has the potential to discover new RNA virus genomes by sequencing total intracellular RNA pools. In this study, we have searched publicly available transcriptomes for sequences similar to viruses of the Nidovirales order. We report two potential nidovirus genomes, a highly divergent 35.9 kb likely complete genome from the California sea hare Aplysia californica, which we assign to a nidovirus named Aplysia abyssovirus 1 (AAbV), and a coronavirus-like 22.3 kb partial genome from the ornamented pygmy frog Microhyla fissipes, which we assign to a nidovirus named Microhyla alphaletovirus 1 (MLeV). AAbV was shown to encode a functional main proteinase, and a translational readthrough signal. Phylogenetic analysis suggested that AAbV represents a new family, proposed here as Abyssoviridae. MLeV represents a sister group to the other known coronaviruses. The importance of MLeV and AAbV for understanding nidovirus evolution, and the origin of terrestrial nidoviruses are discussed.

Immune repertoire in the transcriptome of Littorina littorea reveals new trends in lophotrochozoan proto-complement evolution

The evolution of complement system in invertebrates is poorly investigated. While the repertoire of complement genes in several Ecdysozoa lineages is found substantially different from that of Deuterostomia, the composition and function of the complement in the second protostome lineage, Lophotrochozoa, remains unclear. Here we report the general description of new transcriptomic data on the common periwinkle, Littorina littorea, and trace the evolutionary trajectories of the ancestral proto-complement repertoire. The repertoire is defined as immune cascade providing the minimum set of C3-associated molecules required for C3b amplification, opsonization of the targets and their phagocytosis: thioester protein (TEP) C3, serine protease C2/factor B (Bf) and complement receptors (CR). The reference transcriptome of L. littorea was built from the dual-species RNA-seq experiment with the periwinkle and its tissue digenean parasite Himasthla elongata. Five TEPs, including the ortholog of the C3, are found expressed in the in the mollusk's inflamed tissues. The homolog of the complement receptors CR1/CR2 is also expressed, however the ortholog of Bf is not. The extensive phylogenetic analysis showed that the C3 ortholog and the complement receptors are retained in all key lophotrochozoan taxa: Mollusca, Annelida and Brachiopoda. However, the Bf ortholog was lost at least three times independently in different lineages: i) Cephalopoda, ii) a common ancestor of all Gastropoda and iii) one of the Annelida lineage, Clitellata. Both C3 and Bf molecules were retained in bivalve species, brachiopods and annelid worms from the Polychaeta lineage. Hypothetically, the function of the lost Bf in these animals can be compensated by Factor L (Lf) - the serine protease first found in L. littorea and homologous to both, the Bf and the arthropod factor C (Cf). The contrast differences in proto-complement repertoire between the sister mollusk' taxa, Bivalvia and Gastropoda (the conserved and modified sets, respectively), can underlie differences in their susceptibility to digenean infection.

Polymorphic adaptations in metazoans to establish and maintain photosymbioses

Mutualistic symbioses are common throughout the animal kingdom. Rather unusual is a form of symbiosis, photosymbiosis, where animals are symbiotic with photoautotrophic organisms. Photosymbiosis is found among sponges, cnidarians, flatworms, molluscs, ascidians and even some amphibians. Generally the animal host harbours a phototrophic partner, usually a cyanobacteria or a unicellular alga. An exception to this rule is found in some sea slugs, which only retain the chloroplasts of the algal food source and maintain them photosynthetically active in their own cytosol - a phenomenon called 'functional kleptoplasty'. Research has focused largely on the biodiversity of photosymbiotic species across a range of taxa. However, many questions with regard to the evolution of the ability to establish and maintain a photosymbiosis are still unanswered. To date, attempts to understand genome adaptations which could potentially lead to the evolution of photosymbioses have only been performed in cnidarians. This knowledge gap for other systems is mainly due to a lack of genetic information, both for non-symbiotic and symbiotic species. Considering non-photosymbiotic species is, however, important to understand the factors that make symbiotic species so unique. Herein we provide an overview of the diversity of photosymbioses across the animal kingdom and discuss potential scenarios for the evolution of this association in different lineages. We stress that the evolution of photosymbiosis is probably based on genome adaptations, which (i) lead to recognition of the symbiont to establish the symbiosis, and (ii) are needed to maintain the symbiosis. We hope to stimulate research involving sequencing the genomes of various key taxa to increase the genomic resources needed to understand the most fundamental question: how have animals evolved the ability to establish and maintain a photosymbiosis?

Transcriptome analysis of the threatened snail Ellobium chinense reveals candidate genes for adaptation and identifies SSRs for conservation genetics

Ellobium chinense (Pfeiffer, 1854) is a brackish pulmonate species that inhabits the bases of mangrove trees and is most commonly found in salt grass meadows. Threats to mangrove ecosystems due to habitat degradation and overexploitation have threatened the species with extinction. In South Korea, E. chinense has been assessed as vulnerable, but there are limited data on its population structure and distribution. The nucleotide and protein sequences for this species are not available in databases, which limits the understanding of adaptation-related traits. We sequenced an E. chinense cDNA library using the Illumina platform, and the subsequent bioinformatics analysis yielded 227,032 unigenes. Of these unigenes, 69,088 were annotated to matched protein and nucleotide sequences in databases, for an annotation rate of 30.42%. Among the predominant gene ontology terms, cellular and metabolic processes (under the biological process category), membrane and cell (under the cellular component category), and binding and catalytic activity (under the molecular function category) were noteworthy. In addition, 4850 unigenes were distributed to 15 Kyoto Encyclopaedia of Genes and Genomes based enrichment categories. Among the candidate genes related to adaptation, angiotensin I converting enzyme, adenylate cyclase activating polypeptide, and AMP-activated protein kinase were the most prominent. A total of 15,952 simple sequence repeats (SSRs) were identified in sequences of > 1 kb in length. The di- and trinucleotide repeat motifs were the most common. Among the repeat motif types, AG/CT, AC/GT, and AAC/GTT dominated. Our study provides the first comprehensive genomics dataset for E. chinense, which favors conservation programs for the restoration of the species and provides sufficient evidence for genetic variability among the wild populations.

Deviating from the Norm: Peculiarities of Aplysia cf. californica Early Cleavage Compared to Traditional Spiralian Models

Spiralia represents one of the main clades of bilaterally symmetrical metazoans (Bilateria). This group is of particular interest due to the remarkable conservation of its early developmental pattern despite of the high diversity of larval and adult body plans. Variations during embryogenesis are considered powerful tools to determine ancestral and derived characters under a phylogenetic framework. By direct observation of embryos cultured in vitro, we analyzed the early cleavage of the euopisthobranchs Aplysia cf. californica. We used tubulin immunocytochemistry to stain mitotic spindles during early cleavages, and followed each division with the aid of an autofluorescent compound inside yolk platelets, which differed from the characteristic pink-brownish pigment of the vegetal cytoplasm in zygotes and early embryos. We found that this species exhibits an unequal cleavage characterized by ooplasmic segregation, oblique inclination of mitotic spindles, and differences in size and positioning of the asters in relation to the cellular cortex. Furthermore, we detected asynchrony in cleavage timing between the two large macromeres C and D, which increases the number of cleavage rounds required to reach a particular cell stage in comparison to other spiralians. Here, we report the presence of a transient and previously undescribed U-shaped embryo in this species. The present detailed description of A. californica early development deviates considerably from stereotypical patterns described in other spiralians. Our observations demonstrate that early spiralian development can be more plastic than previously thought.

Sequencing and de novo assembly of visceral mass transcriptome of the critically endangered land snail Satsuma myomphala: Annotation and SSR discovery

Satsuma myomphala is critically endangered through loss of natural habitats, predation by natural enemies, and indiscriminate collection. It is a protected species in Korea but lacks genomic resources for an understanding of varied functional processes attributable to evolutionary success under natural habitats. For assessing the genetic information of S. myomphala, we performed for the first time, de novo transcriptome sequencing and functional annotation of expressed sequences using Illumina Next-Generation Sequencing (NGS) platform and bioinformatics analysis. We identified 103,774 unigenes of which 37,959, 12,890, and 17,699 were annotated in the PANM (Protostome DB), Unigene, and COG (Clusters of Orthologous Groups) databases, respectively. In addition, 14,451 unigenes were predicted under Gene Ontology functional categories, with 4581 assigned to a single category. Furthermore, 3369 sequences with 646 having Enzyme Commission (EC) numbers were mapped to 122 pathways in the Kyoto Encyclopedia of Genes and Genomes Pathway database. The prominent protein domains included the Zinc finger (C2H2-like), Reverse Transcriptase, Thioredoxin-like fold, and RNA recognition motif domain. Many unigenes with homology to immunity, defense, and reproduction-related genes were screened in the transcriptome. We also detected 3120 putative simple sequence repeats (SSRs) encompassing dinucleotide to hexanucleotide repeat motifs from >1kb unigene sequences. A list of PCR primers of SSR loci have been identified to study the genetic polymorphisms. The transcriptome data represents a valuable resource for further investigations on the species genome structure and biology. The unigenes information and microsatellites would provide an indispensable tool for conservation of the species in natural and adaptive environments.

Transcriptome sequencing and de novo characterization of Korean endemic land snail, Koreanohadra kurodana for functional transcripts and SSR markers

The Korean endemic land snail Koreanohadra kurodana (Gastropoda: Bradybaenidae) found in humid areas of broadleaf forests and shrubs have been considered vulnerable as the number of individuals are declining in recent years. The species is poorly characterized at the genomic level that limits the understanding of functions at the molecular and genetics level. In the present study, we performed de novo transcriptome sequencing to produce a comprehensive transcript dataset of visceral mass tissue of K. kurodana by the Illumina paired-end sequencing technology. Over 234 million quality reads were assembled to a total of 315,924 contigs and 191,071 unigenes, with an average and N50 length of 585.6 and 715 bp and 678 and 927 bp, respectively. Overall, 36.32 % of the unigenes found matches to known protein/nucleotide sequences in the public databases. The direction of the unigenes to functional categories was determined using COG, GO, KEGG, and InterProScan protein domain search. The GO analysis search resulted in 22,967 unigenes (12.02 %) being categorized into 40 functional groups. The KEGG annotation revealed that metabolism pathway genes were enriched. The most prominent protein motifs include the zinc finger, ribonuclease H, reverse transcriptase, and ankyrin repeat domains. The simple sequence repeats (SSRs) identified from >1 kb length of unigenes show a dominancy of dinucleotide repeat motifs followed with tri- and tetranucleotide motifs. A number of unigenes were putatively assessed to belong to adaptation and defense mechanisms including heat shock proteins 70, Toll-like receptor 4, AMP-activated protein kinase, aquaporin-2, etc. Our data provide a rich source for the identification and functional characterization of new genes and candidate polymorphic SSR markers in K. kurodana. The availability of transcriptome information ( http://bioinfo.sch.ac.kr/submission/ ) would promote the utilization of the resources for phylogenetics study and genetic diversity assessment.

Transcriptomic Analysis of the Endangered Neritid Species Clithon retropictus: De Novo Assembly, Functional Annotation, and Marker Discovery

An aquatic gastropod belonging to the family Neritidae, Clithon retropictus is listed as an endangered class II species in South Korea. The lack of information on its genomic background limits the ability to obtain functional data resources and inhibits informed conservation planning for this species. In the present study, the transcriptomic sequencing and de novo assembly of C. retropictus generated a total of 241,696,750 high-quality reads. These assembled to 282,838 unigenes with mean and N50 lengths of 736.9 and 1201 base pairs, respectively. Of these, 125,616 unigenes were subjected to annotation analysis with known proteins in Protostome DB, COG, GO, and KEGG protein databases (BLASTX; E ≤ 0.00001) and with known nucleotides in the Unigene database (BLASTN; E ≤ 0.00001). The GO analysis indicated that cellular process, cell, and catalytic activity are the predominant GO terms in the biological process, cellular component, and molecular function categories, respectively. In addition, 2093 unigenes were distributed in 107 different KEGG pathways. Furthermore, 49,280 simple sequence repeats were identified in the unigenes (>1 kilobase sequences). This is the first report on the identification of transcriptomic and microsatellite resources for C. retropictus, which opens up the possibility of exploring traits related to the adaptation and acclimatization of this species.

Developmental Transcriptome Analysis and Identification of Genes Involved in Larval Metamorphosis of the Razor Clam, Sinonovacula constricta

The razor clam Sinonovacula constricta is an important commercial species. The deficiency of developmental transcriptomic data is becoming the bottleneck of further researches on the mechanisms underlying settlement and metamorphosis in early development. In this study, de novo transcriptome sequencing was performed for S. constricta at different early developmental stages by using Illumina HiSeq 2000 paired-end (PE) sequencing technology. A total of 112,209,077 PE clean reads were generated. De novo assembly generated 249,795 contigs with an average length of 585 bp. Gene annotation resulted in the identification of 22,870 unigene hits against the NCBI database. Eight unique sequences related to metamorphosis were identified and analyzed using real-time PCR. The razor clam reference transcriptome would provide useful information on early developmental and metamorphosis mechanisms and could be used in the genetic breeding of shellfish.

De novo Transcriptome Generation and Annotation for Two Korean Endemic Land Snails, Aegista chejuensis and Aegista quelpartensis, Using Illumina Paired-End Sequencing Technology

Aegista chejuensis and Aegista quelpartensis (Family-Bradybaenidae) are endemic to Korea, and are considered vulnerable due to declines in their population. The limited genetic resources for these species restricts the ability to prioritize conservation efforts. We sequenced the transcriptomes of these species using Illumina paired-end technology. Approximately 257 and 240 million reads were obtained and assembled into 198,531 and 230,497 unigenes for A. chejuensis and A. quelpartensis, respectively. The average and N50 unigene lengths were 735.4 and 1073 bp, respectively, for A. chejuensis, and 705.6 and 1001 bp, respectively, for A. quelpartensis. In total, 68,484 (34.5%) and 77,745 (33.73%) unigenes for A. chejuensis and A. quelpartensis, respectively, were annotated to databases. Gene Ontology terms were assigned to 23,778 (11.98%) and 26,396 (11.45) unigenes, for A. chejuensis and A. quelpartensis, respectively, while 5050 and 5838 unigenes were mapped to 117 and 124 pathways in the Kyoto Encyclopedia of Genes and Genomes database. In addition, we identified and annotated 9542 and 10,395 putative simple sequence repeats (SSRs) in unigenes from A. chejuensis and A. quelpartensis, respectively. We designed a list of PCR primers flanking the putative SSR regions. These microsatellites may be utilized for future phylogenetics and conservation initiatives.

Transcriptome and Gene Expression Analysis of Cylas formicarius (Coleoptera: Brentidae) During Different Development Stages

The sweet potato weevil, Cylas formicarius (F.) (Coleoptera: Brentidae), is an important pest of sweet potato worldwide. However, there is limited knowledge on the molecular mechanisms underlying growth and differentiation of C. formicarius The transcriptomes of the eggs, second instar larvae, third instar larvae (L3), pupae, females, and males of C. formicarius were sequenced using Illumina sequencing technology for obtaining global insights into developing transcriptome characteristics and elucidating the relative functional genes. A total of 54,255,544 high-quality reads were produced, trimmed, and de novo assembled into 115,281 contigs. 61,686 unigenes were obtained, with an average length of 1,009 nt. Among these unigenes, 17,348 were annotated into 59 Gene Ontology (GO) terms and 12,660 were assigned to 25 Cluster of Orthologous Groups classes, whereas 24,796 unigenes were mapped to 258 pathways. Differentially expressed unigenes between various developmental stages of C. formicarius were detected. Higher numbers of differentially expressed genes (DEGs) were recorded in the eggs versus L3 and eggs versus male samples (2,141 and 2,058 unigenes, respectively) than the others. Genes preferentially expressed in each stage were also identified. GO and pathway-based enrichment analysis were used to further investigate the functions of the DEGs. In addition, the expression profiles of ten DEGs were validated by quantitative real-time PCR. The transcriptome profiles presented in this study and these DEGs detected by comparative analysis of different developed stages of C. formicarius will facilitate the understanding of the molecular mechanism of various living process and will contribute to further genome-wide research.

Structure, phylogeny, and expression of the frizzled-related gene family in the lophotrochozoan annelid Platynereis dumerilii

Background

Wnt signaling pathways are highly conserved signal transduction pathways important for axis formation, cell fate specification, and organogenesis throughout metazoan development. Within the various Wnt pathways, the frizzled transmembrane receptors (Fzs) and secreted frizzled-related proteins (sFRPs) play central roles in receiving and antagonizing Wnt signals, respectively. Despite their importance, very little is known about the frizzled-related gene family (fzs & sfrps) in lophotrochozoans, especially during early stages of spiralian development. Here we ascertain the frizzled-related gene complement in six lophotrochozoan species, and determine their spatial and temporal expression pattern during early embryogenesis and larval stages of the marine annelid Platynereis dumerilii.

Results

Phylogenetic analyses confirm conserved homologs for four frizzled receptors (Fz1/2/7, Fz4, Fz5/8, Fz9/10) and sFRP1/2/5 in five of six lophotrochozoan species. The sfrp3/4 gene is conserved in one, divergent in two, and evidently lost in three lophotrochozoan species. Three novel fz-related genes (fzCRD1-3) are unique to Platynereis. Transcriptional profiling and in situ hybridization identified high maternal expression of fz1/2/7, expression of fz9/10 and fz1/2/7 within animal and dorsal cell lineages after the 32-cell stage, localization of fz5/8, sfrp1/2/5, and fzCRD-1 to animal-pole cell lineages after the 80-cell stage, and no expression for fz4, sfrp3/4, and fzCRD-2, and -3 in early Platynereis embryos. In later larval stages, all frizzled-related genes are expressed in distinct patterns preferentially in the anterior hemisphere and less in the developing trunk.

Conclusions

Lophotrochozoans have retained a generally conserved ancestral bilaterian frizzled-related gene complement (four Fzs and two sFRPs). Maternal expression of fz1/2/7, and animal lineage-specific expression of fz5/8 and sfrp1/2/5 in early embryos of Platynereis suggest evolutionary conserved roles of these genes to perform Wnt pathway functions during early cleavage stages, and the early establishment of a Wnt inhibitory center at the animal pole, respectively. Numerous frizzled receptor-expressing cells and embryonic territories were identified that might indicate competence to receive Wnt signals during annelid development. An anterior bias for frizzled-related gene expression in embryos and larvae might point to a polarity of Wnt patterning systems along the anterior–posterior axis of this annelid.

Electronic supplementary material

The online version of this article (doi:10.1186/s13227-015-0032-4) contains supplementary material, which is available to authorized users.

Identification of genes related to learning and memory in the brain transcriptome of the mollusc, Hermissenda crassicornis

The sea slug Hermissenda crassicornis (Mollusca, Gastropoda, Nudibranchia) has been studied extensively in associative learning paradigms. However, lack of genetic information previously hindered molecular-level investigations. Here, the Hermissenda brain transcriptome was sequenced and assembled de novo, producing 165,743 total transcripts. Orthologs of 95 genes implicated in learning were identified. These included genes for a serotonin receptor and a GABA-B receptor subunit that had not been previously described in molluscs, as well as an adenylyl cyclase gene not previously described in gastropods. This study illustrates the Hermissenda transcriptome's potential as an important genetic tool in future learning and memory research.

Transcriptome Changes during the Life Cycle of the Red Sponge, Mycale phyllophila (Porifera, Demospongiae, Poecilosclerida)

Sponges are an ancient metazoan group with broad ecological, evolutionary, and biotechnological importance. As in other marine invertebrates with a biphasic life cycle, the developing sponge undergoes a significant morphological, physiological, and ecological transformation during settlement and metamorphosis. In this study, we compare new transcriptome datasets for three life cycle stages of the red sponge (Mycale phyllophila) to test whether gene expression (as in the model poriferan, Amphimedon queenslandica) also varies more after settlement and metamorphosis. In contrast to A. queenslandica, we find that the transcriptome of M. phyllophila changes more during the earlier pre-competent larva/post-larva transition that spans these defining events. We also find that this transition is marked by a greater frequency of significantly up-regulated Gene Ontology terms including those for morphogenesis, differentiation, and development and that the transcriptomes of its pre-competent larvae and adult are distinct. The life cycle transcriptome variation between M. phyllophila and A. queenslandica may be due to their long separate evolutionary histories and corresponding differences in developmental rates and timing. This study now calls for new transcriptome datasets of M. phyllophila and other sponges, which will allow for tests of the generality of our life cycle expression differences and for the greater exploitation of poriferans in both basic and applied research.

Redefining metamorphosis in spiny lobsters: molecular analysis of the phyllosoma to puerulus transition in Sagmariasus verreauxi

The molecular understanding of crustacean metamorphosis is hindered by small sized individuals and inability to accurately define molt stages. We used the spiny lobster Sagmariasus verreauxi where the large, transparent larvae enable accurate tracing of the transition from a leaf-shaped phyllosoma to an intermediate larval-juvenile phase (puerulus). Transcriptomic analysis of larvae at well-defined stages prior to, during, and following this transition show that the phyllosoma-puerulus metamorphic transition is accompanied by vast transcriptomic changes exceeding 25% of the transcriptome. Notably, genes previously identified as regulating metamorphosis in other crustaceans do not fluctuate during this transition but in the later, morphologically-subtle puerulus-juvenile transition, indicating that the dramatic phyllosoma-puerulus morphological shift relies on a different, yet to be identified metamorphic mechanism. We examined the change in expression of domains and gene families, with focus on several key genes. Our research implies that the separation in molecular triggering systems between the phyllosoma-puerulus and puerulus-juvenile transitions might have enabled the extension of the oceanic phase in spiny lobsters. Study of similar transitions, where metamorphosis is uncoupled from the transition into the benthic juvenile form, in other commercially important crustacean groups might show common features to point on the evolutionary advantage of this two staged regulation.

Transcriptome-wide analysis of the response of the thecosome pteropod Clio pyramidata to short-term CO2 exposure

Thecosome pteropods, a group of calcifying holoplanktonic mollusks, have recently become a research focus due to their potential sensitivity to increased levels of anthropogenic dissolved CO2 in seawater and the accompanying ocean acidification. Some populations, however, already experience high CO2 in their natural distribution during diel vertical migrations. To achieve a better understanding of the mechanisms of pteropod calcification and physiological response to this sort of short duration CO2 exposure, we characterized the gene complement of Clio pyramidata, a cosmopolitan diel migratory thecosome, and investigated its transcriptomic response to experimentally manipulated CO2 conditions. Individuals were sampled from the Northwest Atlantic in the fall of 2011 and exposed to ambient conditions (~380ppm) and elevated CO2 (~800ppm, similar to levels experienced during a diel vertical migration) for ~10h. Following this exposure the respiration rate of the individuals was measured. We then performed RNA-seq analysis, assembled the C. pyramidata transcriptome de novo, annotated the genes, and assessed the differential gene expression patterns in response to exposure to elevated CO2. Within the transcriptome, we identified homologs of genes with known roles in biomineralization in other mollusks, including perlucin, calmodulin, dermatopontin, calponin, and chitin synthases. Respiration rate was not affected by short-term exposure to CO2. Gene expression varied greatly among individuals, and comparison between treatments indicated that C. pyramidata down-regulated a small number of genes associated with aerobic metabolism and up-regulated genes that may be associated with biomineralization, particularly collagens and C-type lectins. These results provide initial insight into the effects of short term CO2 exposure on these important planktonic open-ocean calcifiers, pairing respiration rate and the gene expression level of response, and reveal candidate genes for future ecophysiological, biomaterial and phylogenetic studies.

Lectin-like molecules in transcriptome of Littorina littorea hemocytes

The common periwinkle Littorina littorea was introduced in the list of models for comparative immunobiology as a representative of phylogenetically important taxon Caenogastropoda. Using Illumina sequencing technology, we de novo assembled the transcriptome of Littorina littorea hemocytes from 182 million mRNA-Seq pair-end 100 bp reads into a total of 15,526 contigs clustered in 4472 unigenes. The transcriptome profile was analyzed for presence of carbohydrate-binding molecules in a variety of architectural contexts. Hemocytes' repertoire of lectin-like proteins bearing conserved carbohydrate-recognition domains (CRDs) is highly diversified, including 11 of 15 lectin families earlier described in animals, as well as the novel members of lectin family found for the first time in mollusc species. The new molluscan lineage-specific domain combinations were confirmed by cloning and sequencing, including the fuco-lectin related molecules (FLReMs) composed of N-terminal region with no sequence homology to any known protein, a middle Fucolectin Tachylectin-4 Pentaxrin (FTP) domain, and a C-terminal epidermal growth factor (EGF) repeat region. The repertoire of lectin-like molecules is discussed in terms of their potential participation in the receptor phase of immune response. In total, immune-associated functions may be attributed to 70 transcripts belonging to 6 lectin families. These lectin-like genes show low overlap between species of invertebrates, suggesting relatively rapid evolution of immune-associated genes in the group. The repertoire provides valuable candidates for further characterization of the gene functions in mollusc immunity.

Comparative analysis of early ontogeny in Bursatella leachii and Aplysia californica

Opisthobranch molluscs exhibit fascinating body plans associated with the evolution of shell loss in multiple lineages. Sea hares in particular are interesting because Aplysia californica is a well-studied model organism that offers a large suite of genetic tools. Bursatella leachii is a related tropical sea hare that lacks a shell as an adult and therefore lends itself to comparative analysis with A. californica. We have established an enhanced culturing procedure for B. leachii in husbandry that enabled the study of shell formation and loss in this lineage with respect to A. californica life staging.

A detailed staging scheme for late larval development in Strongylocentrotus purpuratus focused on readily-visible juvenile structures within the rudiment

BACKGROUND

The purple sea urchin, Strongylocentrotus purpuratus, has long been the focus of developmental and ecological studies, and its recently-sequenced genome has spawned a diversity of functional genomics approaches. S. purpuratus has an indirect developmental mode with a pluteus larva that transforms after 1-3 months in the plankton into a juvenile urchin. Compared to insects and frogs, mechanisms underlying the correspondingly dramatic metamorphosis in sea urchins remain poorly understood. In order to take advantage of modern techniques to further our understanding of juvenile morphogenesis, organ formation, metamorphosis and the evolution of the pentameral sea urchin body plan, it is critical to assess developmental progression and rate during the late larval phase. This requires a staging scheme that describes developmental landmarks that can quickly and consistently be used to identify the stage of individual living larvae, and can be tracked during the final two weeks of larval development, as the juvenile is forming.

RESULTS

Notable structures that are easily observable in developing urchin larvae are the developing spines, test and tube feet within the juvenile rudiment that constitute much of the oral portion of the adult body plan. Here we present a detailed staging scheme of rudiment development in the purple urchin using soft structures of the rudiment and the primordia of these juvenile skeletal elements. We provide evidence that this scheme is robust and applicable across a range of temperature and feeding regimes.

CONCLUSIONS

Our proposed staging scheme provides both a useful method to study late larval development in the purple urchin, and a framework for developing similar staging schemes across echinoderms. Such efforts will have a high impact on evolutionary developmental studies and larval ecology, and facilitate research on this important deuterostome group.

Regeneration of Aplysia bag cell neurons is synergistically enhanced by substrate-bound hemolymph proteins and laminin

We have investigated Aplysia hemolymph as a source of endogenous factors to promote regeneration of bag cell neurons. We describe a novel synergistic effect between substrate-bound hemolymph proteins and laminin. This combination increased outgrowth and branching relative to either laminin or hemolymph alone. Notably, the addition of hemolymph to laminin substrates accelerated growth cone migration rate over ten-fold. Our results indicate that the active factor is either a high molecular weight protein or protein complex and is not the respiratory protein hemocyanin. Substrate-bound factor(s) from central nervous system-conditioned media also had a synergistic effect with laminin, suggesting a possible cooperation between humoral proteins and nervous system extracellular matrix. Further molecular characterization of active factors and their cellular targets is warranted on account of the magnitude of the effects reported here and their potential relevance for nervous system repair.

Reproduction-related genes in the pearl oyster genome

Molluscan reproduction has been a target of biological research because of the various reproductive strategies that have evolved in this phylum. It has also been studied for the development of fisheries technologies, particularly aquaculture. Although fundamental processes of reproduction in other phyla, such as vertebrates and arthropods, have been well studied, information on the molecular mechanisms of molluscan reproduction remains limited. The recently released draft genome of the pearl oyster Pinctada fucata provides a novel and powerful platform for obtaining structural information on the genes and proteins involved in bivalve reproduction. In the present study, we analyzed the pearl oyster draft genome to screen reproduction-related genes. Analysis was mainly conducted for genes reported from other molluscs for encoding orthologs of reproduction-related proteins in other phyla. The gene search in the P. fucata gene models (version 1.1) and genome assembly (version 1.0) were performed using Genome Browser and BLAST software. The obtained gene models were then BLASTP searched against a public database to confirm the best-hit sequences. As a result, more than 40 gene models were identified with high accuracy to encode reproduction-related genes reported for P. fucata and other molluscs. These include vasa, nanos, doublesex- and mab-3-related transcription factor, 5-hydroxytryptamine (5-HT) receptors, vitellogenin, estrogen receptor, and others. The set of reproduction-related genes of P. fucata identified in the present study constitute a new tool for research on bivalve reproduction at the molecular level.

A genome-wide analysis of annexins from parasitic organisms and their vectors

In this study, we conduct an in-depth analysis of annexin proteins from a diverse range of invertebrate taxa, including the major groups that contain the parasites and vector organisms that are harmful to humans and domestic animals. Using structure-based amino acid sequence alignments and phylogenetic analyses, we present a classification for this protein group and assign names to sequences with ambiguous annotations in public databases. Our analyses reveal six distinct annexin clades, and the mapping of genes encoding annexins to the genome of the human blood fluke Schistosoma mansoni supports the hypothesis of gene duplication as a major evolutionary event in annexin genesis. This study illuminates annexin diversity from a novel perspective using contemporary phylogenetic hypotheses of eukaryote evolution, and will aid the consolidation of annexin protein identities in public databases and provide a foundation for future functional analysis and characterisation of these proteins in parasites of socioeconomic importance.

Characterization of GdFFD, a D-amino acid-containing neuropeptide that functions as an extrinsic modulator of the Aplysia feeding circuit

During eukaryotic translation, peptides/proteins are created using L-amino acids. However, a D-amino acid-containing peptide (DAACP) can be produced through post-translational modification via an isomerase enzyme. General approaches to identify novel DAACPs and investigate their function, particularly in specific neural circuits, are lacking. This is primarily due to the difficulty in characterizing this modification and due to the limited information on neural circuits in most species. We describe a multipronged approach to overcome these limitations using the sea slug Aplysia californica. Based on bioinformatics and homology to known DAACPs in the land snail Achatina fulica, we targeted two predicted peptides in Aplysia, GFFD, similar to achatin-I (GdFAD versus GFAD, where dF stands for D-phenylalanine), and YAEFLa, identical to fulyal (YdAEFLa versus YAEFLa), using stereoselective analytical methods, i.e. MALDI MS fragmentation analysis and LC-MS/MS. Although YAEFLa in Aplysia was detected only in an all L-form, we found that both GFFD and GdFFD were present in the Aplysia CNS. In situ hybridization and immunolabeling of GFFD/GdFFD-positive neurons and fibers suggested that GFFD/GdFFD might act as an extrinsic modulator of the feeding circuit. Consistent with this hypothesis, we found that GdFFD induced robust activity in the feeding circuit and elicited egestive motor patterns. In contrast, the peptide consisting of all L-amino acids, GFFD, was not bioactive. Our data indicate that the modification of an L-amino acid-containing neuropeptide to a DAACP is essential for peptide bioactivity in a motor circuit, and thus it provides a functional significance to this modification.

Analysis of a deep transcriptome from the mantle tissue of Patella vulgata Linnaeus (Mollusca: Gastropoda: Patellidae) reveals candidate biomineralising genes

The gastropod Patella vulgata is abundant on rocky shores in Northern Europe and a significant grazer of intertidal algae. Here we report the application of Illumina sequencing to develop a transcriptome from the adult mantle tissue of P. vulgata. We obtained 47,237,104 paired-end reads of 51 bp, trialled de novo assembly methods and settled on the additive multiple K method followed by redundancy removal as resulting in the most comprehensive assembly. This yielded 29,489 contigs of at least 500 bp in length. We then used three methods to search for candidate genes relevant to biomineralisation: searches via BLAST and Hidden Markov Models for homologues of biomineralising genes from other molluscs, searches for predicted proteins containing tandem repeats and searches for secreted proteins that lacked a transmembrane domain. From the results of these searches we selected 15 contigs for verification by RT-PCR, of which 14 were successfully amplified and cloned. These included homologues of Pif-177/BSMP, Perlustrin, SPARC, AP24, Follistatin-like and Carbonic anhydrase, as well as three containing extensive G-X-Y repeats as found in nacrein. We selected two for further verification by in situ hybridisation, demonstrating expression in the larval shell field. We conclude that de novo assembly of Illumina data offers a cheap and rapid route to a predicted transcriptome that can be used as a resource for further biological study.

Suppressive subtractive hybridisation transcriptomics provides a novel insight into the functional role of the hypobranchial gland in a marine mollusc

The hypobranchial gland present in gastropods is an organ whose function is not clearly understood. Involved in mucus production, within members of the family Muricidae it is also the source of the ancient dye Tyrian purple and its bioactive precursors. To gain further insights into hypobranchial gland biology, suppressive subtractive hybridisation was performed on hypobranchial gland and mantle tissue from the marine snail Dicathais orbita creating a differentially expressed cDNA library. 437 clones were randomly sequenced, analysed and annotated and 110 sequences had their functions putatively identified. Importantly this approach identified a putative gene involved in Tyrian purple biosynthesis, an arylsulphatase gene. Confirmation of the upregulation of arylsulphatase in the hypobranchial gland compared to the mantle was demonstrated using quantitative real-time PCR. Other genes identified as playing an important role in the hypobranchial gland were those involved in mucus protein synthesis, choline ester regulation, protein and energy production. This study confirms that the hypobranchial gland is involved in the production of mucus secretion and also identifies it as a site of chemical interaction and biosynthesis. This study lays the foundation for a better understanding of the enzymatic production of Tyrian purple precursors within the gland.

Post-embryonic transcriptomes of the prawn Macrobrachium rosenbergii: multigenic succession through metamorphosis

Like many metazoans, the freshwater prawn Macrobrachium rosenbergii begins its post-embryonic life with a set of morphologically distinct planktonic larval stages, followed by a benthic post-larval stage during which the maturing organism differs from the larvae both ecologically and physiologically. Understanding of the molecular basis underlying morphogenesis in crustaceans is limited to the observation that methyl farnesoate, the non-epoxidated form of the insect juvenile hormone, acts as the active crustacean juvenoid. Molt steroids were also linked to morphogenesis and several other molecular pathways, such as Hedgehog and Wnt, are known to underlie morphogenesis in all metazoans examined and, as such, are thought to do the same in crustaceans. Using next generation sequencing, we deep-sequenced the transcriptomes of several larval and post-larval stages. De novo assembly, followed by bioinformatics analysis, revealed that many novel transcripts are over-expressed in either larvae- or post-larvae-stage prawn, shedding light on the molecular basis underlying M. rosenbergii metamorphosis. Fast larval molting rates and periodic morphological changes were reflected in over-expression of transcripts annotated to the cell cycle, DNA replication and morphogenic pathways (i.e., Hedgehog and Wnt). Further characterization of transcripts assigned to morphogenic pathways by real-time RT-PCR reconfirmed their over-expression in larvae, albeit with a more complex expression pattern when examined in the individual developmental stages. The expression level of an orthologue of cytochrome P450, 15A1, known to epoxidize methyl farnesoate in insects, was increased in the late larval and early post-larval stages, in accordance with the role of methyl farnesoate in crustacean metamorphosis. This study exemplifies the applicability of a high-throughput sequencing approach for studying complex traits, including metamorphosis, providing new insight into this unexplored area of crustacean research.

Single-neuron transcriptome and methylome sequencing for epigenomic analysis of aging

Enormous heterogeneity in transcription and signaling is the feature that slows down progress in our understanding of the mechanisms of normal aging and age-related diseases. This is critical for neurobiology of aging where the enormous diversity of neuronal populations presents a significant challenge in experimental design. Here, we introduce Aplysia as a model for genomic analysis of aging at the single-cell level and provide protocols for integrated transcriptome and methylome profiling of individually identified neurons during the aging process. These single-cell RNA-seq and DNA methylation assays (methyl-capture/methyl enrichment) are compatible with all major next generation sequencing platforms (we used Roche/454 and SOLiD/Life Technologies as illustrative examples) and can be used to integrate an epigenetic signature with transcriptional output. The described sequencing library construction protocol provides both quantitative and directional information from transcriptional profiling of individual cells. Our results also confirm that different copies of DNA in polyploid Aplysia neurons behave similarly with respect to their DNA methylation.

Single-cell semiconductor sequencing

RNA-seq or transcriptome analysis of individual cells and small-cell populations is essential for virtually any biomedical field. It is especially critical for developmental, aging, and cancer biology as well as neuroscience where the enormous heterogeneity of cells present a significant methodological and conceptual challenge. Here we present two methods that allow for fast and cost-efficient transcriptome sequencing from ultra-small amounts of tissue or even from individual cells using semiconductor sequencing technology (Ion Torrent, Life Technologies). The first method is a reduced representation sequencing which maximizes capture of RNAs and preserves transcripts' directionality. The second, a template-switch protocol, is designed for small mammalian neurons. Both protocols, from cell/tissue isolation to final sequence data, take up to 4 days. The efficiency of these protocols has been validated with single hippocampal neurons and various invertebrate tissues including individually identified neurons within a simpler memory-forming circuit of Aplysia californica and early (1-, 2-, 4-, 8-cells) embryonic and developmental stages from basal metazoans.

Pyrosequencing of Haliotis diversicolor transcriptomes: insights into early developmental molluscan gene expression

Background

The abalone Haliotis diversicolor is a good model for study of the settlement and metamorphosis, which are widespread marine ecological phenomena. However, information on the global gene backgrounds and gene expression profiles for the early development of abalones is lacking.

Methodology/Principal Findings

In this study, eight non-normalized and multiplex barcode-labeled transcriptomes were sequenced using a 454 GS system to cover the early developmental stages of the abalone H. diversicolor. The assembly generated 35,415 unigenes, of which 7,566 were assigned GO terms. A global gene expression profile containing 636 scaffolds/contigs was constructed and was proven reliable using qPCR evaluation. It indicated that there may be existing dramatic phase transitions. Bioprocesses were proposed, including the ‘lock system’ in mature eggs, the collagen shells of the trochophore larvae and the development of chambered extracellular matrix (ECM) structures within the earliest postlarvae.

Conclusion

This study globally details the first 454 sequencing data for larval stages of H. diversicolor. A basic analysis of the larval transcriptomes and cluster of the gene expression profile indicates that each stage possesses a batch of specific genes that are indispensable during embryonic development, especially during the two-cell, trochophore and early postlarval stages. These data will provide a fundamental resource for future physiological works on abalones, revealing the mechanisms of settlement and metamorphosis at the molecular level.

Urotensin II in invertebrates: from structure to function in Aplysia californica

Neuropeptides are ancient signaling molecules that are involved in many aspects of organism homeostasis and function. Urotensin II (UII), a peptide with a range of hormonal functions, previously has been reported exclusively in vertebrates. Here, we provide the first direct evidence that UII-like peptides are also present in an invertebrate, specifically, the marine mollusk Aplysia californica. The presence of UII in the central nervous system (CNS) of Aplysia implies a more ancient gene lineage than vertebrates. Using representational difference analysis, we identified an mRNA of a protein precursor that encodes a predicted neuropeptide, we named Aplysia urotensin II (apUII), with a sequence and structural similarity to vertebrate UII. With in-situ hybridization and immunohistochemistry, we mapped the expression of apUII mRNA and its prohormone in the CNS and localized apUII-like immunoreactivity to buccal sensory neurons and cerebral A-cluster neurons. Mass spectrometry performed on individual isolated neurons, and tandem mass spectrometry on fractionated peptide extracts, allowed us to define the posttranslational processing of the apUII neuropeptide precursor and confirm the highly conserved cyclic nature of the mature neuropeptide apUII. Electrophysiological analysis of the central effects of a synthetic apUII suggests it plays a role in satiety and/or aversive signaling in feeding behaviors. Finding the homologue of vertebrate UII in the numerically small CNS of an invertebrate animal model is important for gaining insights into the molecular mechanisms and pathways mediating the bioactivity of UII in the higher metazoan.

Distinct expression patterns of glycoprotein hormone subunits in the lophotrochozoan Aplysia: implications for the evolution of neuroendocrine systems in animals

Glycoprotein hormones (GPHs) comprise a group of signaling molecules critical for major metabolic and reproductive functions. In vertebrates they include chorionic gonadotropin, LH, FSH, and TSH. The active hormones are characterized by heterodimerization between a common α and hormone-specific β subunit, which activate leucine-rich repeat-containing G protein coupled receptors. To date, genes referred to as GPHα2 and GPHβ5 have been the only glycoprotein hormone subunits identified in invertebrates, suggesting that other GPHα and GPHβ subunits diversified during vertebrate evolution. Still the functions of GPHα2 and GPHβ5 remain largely unknown for both vertebrates and invertebrates. To further understand the evolution and putative function of these subunits, we cloned and analyzed phylogenetically two glycoprotein subunits, AcaGPHα and AcaGPHβ, from the sea hare Aplysia californica. Model based three-dimensional predictions of AcaGPHβ confirm the presence of a complete cysteine knot, two hairpin loops, and a long loop. As in the human GPHβ5 subunit the seatbelt structure is absent in AcaGPHβ. We also found that AcaGPHα and AcaGPHβ subunits are expressed in larval stages of Aplysia, and we present a detailed expression map of the subunits in the adult central nervous system using in situ hybridizations. Both subunits are expressed in subpopulations of pleural and buccal mechanosensory neurons, suggesting a neuronal modulatory function of these subunits in Aplysia. Furthermore it supports the model of a relatively diffuse neuroendocrine-like system in molluscs, where specific primary sensory neurons release peptides extrasynaptically (paracrine secretion). This is in contrast to vertebrates and insects, in which releasing and stimulating factor from centralized sensory regions of the central nervous system ultimately regulate hormone release in peripheral glands.

Sequencing and de novo analysis of Crassostrea angulata (Fujian oyster) from 8 different developing phases using 454 GSFlx

Research on the mechanism for early development of shellfish, such as body plan, shell formation, settlement and metamorphosis is currently an active research field. However, studies were still limited and not deep enough because of the lack of genomic resources such as genome or transcriptome sequences. In the present research, de novo transcriptome sequencing was performed for Crassostrea angulata, the most economically important cultured oyster species in China, at eight early developmental stages using the 454 sequencing technology. A total of 555,215 reads were produced with an average length of 309 nucleotides that were then assembled into 10,462 contigs. As determined by GO annotation and KEGG pathway mapping, functional annotation of the unigenes recovered diverse biological functions and processes. Six unique sequences related to settlement, metamorphosis and growth were subsequently analyzed by real-time PCR. Given the lack of whole genome information for oysters, transcriptome and de novo analysis of C. angulata from the eight different developing phases will provide important and useful information on early development mechanism and help genetic breeding of shellfish.