Deep mRNA sequencing of the Tritonia diomedea brain transcriptome provides access to gene homologues for neuronal excitability, synaptic transmission and peptidergic signalling

Tracing RNA viruses associated with Nudibranchia gastropods

Background

Nudibranchia is an under-studied taxonomic group of gastropods, including more than 3,000 species with colourful and extravagant body shapes and peculiar predatory and defensive strategies. Although symbiosis with bacteria has been reported, no data are available for the nudibranch microbiome nor regarding viruses possibly associated with these geographically widespread species.

Methods

Based on 47 available RNA sequencing datasets including more than two billion reads of 35 nudibranch species, a meta-transcriptome assembly was constructed. Taxonomic searches with DIAMOND, RNA-dependent-RNA-polymerase identification with palmscan and viral hallmark genes identification by VirSorter2 in combination with CheckV were applied to identify genuine viral genomes, which were then annotated using CAT.

Results

A total of 20 viral genomes were identified as bona fide viruses, among 552 putative viral contigs resembling both RNA viruses of the Negarnaviricota, Pisuviricota, Kitrinoviricota phyla and actively transcribing DNA viruses of the Cossaviricota and Nucleocytoviricota phyla. The 20 commonly identified viruses showed similarity with RNA viruses identified in other RNA-seq experiments and can be putatively associated with bacteria, plant and arthropod hosts by co-occurence analysis. The RNA samples having the highest viral abundances showed a heterogenous and mostly sample-specific distribution of the identified viruses, suggesting that nudibranchs possess diversified and mostly unknown viral communities.

Comparative Analysis of Neuropeptides in Homologous Interneurons and Prohormone Annotation in Nudipleuran Sea Slugs

Despite substantial research on neuronal circuits in nudipleuran gastropods, few peptides have been implicated in nudipleuran behavior. In this study, we expanded the understanding of peptides in this clade, using three species with well-studied nervous systems, Hermissenda crassicornis, Melibe leonina, and Pleurobranchaea californica. For each species, we performed sequence homology analysis of de novo transcriptome predictions to identify homologs to 34 of 36 prohormones previously characterized in the gastropods Aplysia californica and Lymnaea stagnalis. We then used single-cell mass spectrometry to characterize peptide profiles in homologous feeding interneurons: the multifunctional ventral white cell (VWC) in P. californica and the small cardioactive peptide B large buccal (SLB) cells in H. crassicornis and M. leonina. The neurons produced overlapping, but not identical, peptide profiles. The H. crassicornis SLB cells expressed peptides from homologs to the FMRFamide (FMRFa), small cardioactive peptide (SCP), LFRFamide (LFRFa), and feeding circuit activating peptides prohormones. The M. leonina SLB cells expressed peptides from homologs to the FMRFa, SCP, LFRFa, and MIP-related peptides prohormones. The VWC, previously shown to express peptides from the FMRFa and QNFLa (a homolog of A. californica pedal peptide 4) prohormones, was shown to also contain SCP peptides. Thus, each neuron expressed peptides from the FMRFa and SCP families, the H. crassicornis and M. leonina SLB cells expressed peptides from the LFRFa family, and each neuron contained peptides from a prohormone not found in the others. These data suggest each neuron performs complex co-transmission, which potentially facilitates a multifunctional role in feeding. Additionally, the unique feeding characteristics of each species may relate, in part, to differences in the peptide profiles of these neurons. These data add chemical insight to enhance our understanding of the neuronal basis of behavior in nudipleurans and other gastropods.

Transcriptomics provides a robust framework for the relationships of the major clades of cladobranch sea slugs (Mollusca, Gastropoda, Heterobranchia), but fails to resolve the position of the enigmatic genus Embletonia

Background

The soft-bodied cladobranch sea slugs represent roughly half of the biodiversity of marine nudibranch molluscs on the planet. Despite their global distribution from shallow waters to the deep sea, from tropical into polar seas, and their important role in marine ecosystems and for humans (as targets for drug discovery), the evolutionary history of cladobranch sea slugs is not yet fully understood.

Results

To enlarge the current knowledge on the phylogenetic relationships, we generated new transcriptome data for 19 species of cladobranch sea slugs and two additional outgroup taxa (Berthella plumula and Polycera quadrilineata). We complemented our taxon sampling with previously published transcriptome data, resulting in a final data set covering 56 species from all but one accepted cladobranch superfamilies. We assembled all transcriptomes using six different assemblers, selecting those assemblies that provided the largest amount of potentially phylogenetically informative sites. Quality-driven compilation of data sets resulted in four different supermatrices: two with full coverage of genes per species (446 and 335 single-copy protein-coding genes, respectively) and two with a less stringent coverage (667 genes with 98.9% partition coverage and 1767 genes with 86% partition coverage, respectively). We used these supermatrices to infer statistically robust maximum-likelihood trees. All analyses, irrespective of the data set, indicate maximal statistical support for all major splits and phylogenetic relationships at the family level. Besides the questionable position of Noumeaella rubrofasciata, rendering the Facelinidae as polyphyletic, the only notable discordance between the inferred trees is the position of Embletonia pulchra. Extensive testing using Four-cluster Likelihood Mapping, Approximately Unbiased tests, and Quartet Scores revealed that its position is not due to any informative phylogenetic signal, but caused by confounding signal.

Conclusions

Our data matrices and the inferred trees can serve as a solid foundation for future work on the taxonomy and evolutionary history of Cladobranchia. The placement of E. pulchra, however, proves challenging, even with large data sets and various optimization strategies. Moreover, quartet mapping results show that confounding signal present in the data is sufficient to explain the inferred position of E. pulchra, again leaving its phylogenetic position as an enigma.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01944-0.

FMRF-NH2 -related neuropeptides in Biomphalaria spp., intermediate hosts for schistosomiasis: Precursor organization and immunohistochemical localization

Freshwater snails of the genus Biomphalaria serve as intermediate hosts for the digenetic trematode Schistosoma mansoni, the etiological agent for the most widespread form of intestinal schistosomiasis. As neuropeptide signaling in host snails can be altered by trematode infection, a neural transcriptomics approach was undertaken to identify peptide precursors in Biomphalaria glabrata, the major intermediate host for S. mansoni in the Western Hemisphere. Three transcripts that encode peptides belonging to the FMRF-NH2 -related peptide (FaRP) family were identified in B. glabrata. One transcript encoded a precursor polypeptide (Bgl-FaRP1; 292 amino acids) that included eight copies of the tetrapeptide FMRF-NH2 and single copies of FIRF-NH2 , FLRF-NH2 , and pQFYRI-NH2 . The second transcript encoded a precursor (Bgl-FaRP2; 347 amino acids) that comprised 14 copies of the heptapeptide GDPFLRF-NH2 and 1 copy of SKPYMRF-NH2 . The precursor encoded by the third transcript (Bgl-FaRP3; 287 amino acids) recapitulated Bgl-FaRP2 but lacked the full SKPYMRF-NH2 peptide. The three precursors shared a common signal peptide, suggesting a genomic organization described previously in gastropods. Immunohistochemical studies were performed on the nervous systems of B. glabrata and B. alexandrina, a major intermediate host for S. mansoni in Egypt. FMRF-NH2 -like immunoreactive (FMRF-NH2 -li) neurons were located in regions of the central nervous system associated with reproduction, feeding, and cardiorespiration. Antisera raised against non-FMRF-NH2 peptides present in the tetrapeptide and heptapeptide precursors labeled independent subsets of the FMRF-NH2 -li neurons. This study supports the participation of FMRF-NH2 -related neuropeptides in the regulation of vital physiological and behavioral systems that are altered by parasitism in Biomphalaria.

Ion channel profiling of the Lymnaea stagnalis ganglia via transcriptome analysis

BACKGROUND

The pond snail Lymnaea stagnalis (L. stagnalis) has been widely used as a model organism in neurobiology, ecotoxicology, and parasitology due to the relative simplicity of its central nervous system (CNS). However, its usefulness is restricted by a limited availability of transcriptome data. While sequence information for the L. stagnalis CNS transcripts has been obtained from EST libraries and a de novo RNA-seq assembly, the quality of these assemblies is limited by a combination of low coverage of EST libraries, the fragmented nature of de novo assemblies, and lack of reference genome.

RESULTS

In this study, taking advantage of the recent availability of a preliminary L. stagnalis genome, we generated an RNA-seq library from the adult L. stagnalis CNS, using a combination of genome-guided and de novo assembly programs to identify 17,832 protein-coding L. stagnalis transcripts. We combined our library with existing resources to produce a transcript set with greater sequence length, completeness, and diversity than previously available ones. Using our assembly and functional domain analysis, we profiled L. stagnalis CNS transcripts encoding ion channels and ionotropic receptors, which are key proteins for CNS function, and compared their sequences to other vertebrate and invertebrate model organisms. Interestingly, L. stagnalis transcripts encoding numerous putative Ca2+ channels showed the most sequence similarity to those of Mus musculus, Danio rerio, Xenopus tropicalis, Drosophila melanogaster, and Caenorhabditis elegans, suggesting that many calcium channel-related signaling pathways may be evolutionarily conserved.

CONCLUSIONS

Our study provides the most thorough characterization to date of the L. stagnalis transcriptome and provides insights into differences between vertebrates and invertebrates in CNS transcript diversity, according to function and protein class. Furthermore, this study provides a complete characterization of the ion channels of Lymnaea stagnalis, opening new avenues for future research on fundamental neurobiological processes in this model system.

The Distribution and Possible Roles of Small Cardioactive Peptide in the Nudibranch Melibe leonina

The neuropeptide small cardioactive peptide (SCP) plays an integrative role in exciting various motor programs involved in feeding and locomotion in a number of gastropod species. In this study, immunohistochemistry, using monoclonal antibodies against SCPB, was used to localize SCPB-like-immunoreactive neurons in the central nervous system, and map their connections to various tissues, in the nudibranch, Melibe leonina. Approximately 28-36 SCPB-like-immunoreactive neurons were identified in the M. leonina brain, as well as one large neuron in each of the buccal ganglia. The neuropil of the pedal ganglia contained the most SCPB-like-immunoreactive varicosities, although only a small portion of these were due to SCPB-like-immunoreactive neurons in the same ganglion. This suggests that much of the SCPB-like immunoreactivity in the neuropil of the pedal ganglia was from neurons in other ganglia that projected through the pedal-pedal connectives or the connectives from the cerebral and pleural ganglia. We also observed extensive SCPB innervation along the length of the esophagus. Therefore, we investigated the impact of SCPB on locomotion in intact animals, as well as peristaltic contractions of the isolated esophagus. Injection of intact animals with SCPB at night led to a significant increase in crawling and swimming, compared to control animals injected with saline. Furthermore, perfusion of isolated brains with SCPB initiated expression of the swim motor program. Application of SCPB to the isolated quiescent esophagus initiated rhythmic peristaltic contractions, and this occurred in preparations both with and without the buccal ganglia being attached. All these data, taken together, suggest that SCPB could be released at night to arouse animals and enhance the expression of both feeding and swimming motor programs in M. leonina.

Multi-omics investigations within the Phylum Mollusca, Class Gastropoda: from ecological application to breakthrough phylogenomic studies

Gastropods are the largest and most diverse class of mollusc and include species that are well studied within the areas of taxonomy, aquaculture, biomineralization, ecology, microbiome and health. Gastropod research has been expanding since the mid-2000s, largely due to large-scale data integration from next-generation sequencing and mass spectrometry in which transcripts, proteins and metabolites can be readily explored systematically. Correspondingly, the huge data added a great deal of complexity for data organization, visualization and interpretation. Here, we reviewed the recent advances involving gastropod omics ('gastropodomics') research from hundreds of publications and online genomics databases. By summarizing the current publicly available data, we present an insight for the design of useful data integrating tools and strategies for comparative omics studies in the future. Additionally, we discuss the future of omics applications in aquaculture, natural pharmaceutical biodiscovery and pest management, as well as to monitor the impact of environmental stressors.

GABA as a Neurotransmitter in Gastropod Molluscs

The neurotransmitter gamma-aminobutyric acid (GABA) is widely distributed in the mammalian central nervous system, where it acts as a major mediator of synaptic inhibition. GABA also serves as a neurotransmitter in a range of invertebrate phyla, including arthropods, echinoderms, annelids, nematodes, and platyhelminthes. This article reviews evidence supporting the neurotransmitter role of GABA in gastropod molluscs, with an emphasis on its presence in identified neurons and well-characterized neural circuits. The collective findings indicate that GABAergic signaling participates in the selection and specification of motor programs, as well as the bilateral coordination of motor circuits. While relatively few in number, GABAergic neurons can influence neural circuits via inhibitory, excitatory, and modulatory synaptic actions. GABA's colocalization with peptidergic and classical neurotransmitters can broaden its integrative capacity. The functional properties of GABAergic neurons in simpler gastropod systems may provide insight into the role of this neurotransmitter phenotype in more complex brains.

Single neuron serotonin receptor subtype gene expression correlates with behaviour within and across three molluscan species

The marine mollusc, Pleurobranchaea californica varies daily in whether it swims and this correlates with whether serotonin (5-HT) enhances the strength of synapses made by the swim central pattern generator neuron, A1/C2. Another species, Tritonia diomedea, reliably swims and does not vary in serotonergic neuromodulation. A third species, Hermissenda crassicornis, never produces this behaviour and lacks the neuromodulation. We found that expression of particular 5-HT receptor subtype (5-HTR) genes in single neurons correlates with swimming. Orthologues to seven 5-HTR genes were identified from whole-brain transcriptomes. We isolated individual A1/C2 neurons and sequenced their RNA or measured 5-HTR gene expression using absolute quantitative PCR. A1/C2 neurons isolated from Pleurobranchaea that produced a swim motor pattern just prior to isolation expressed 5-HT2a and 5-HT7 receptor genes, as did all Tritonia samples. These subtypes were absent from A1/C2 isolated from Pleurobranchaea that did not swim on that day and from Hermissenda A1/C2 neurons. Expression of other receptors was not correlated with swimming. This suggests that these 5-HTRs may mediate the modulation of A1/C2 synaptic strength and play an important role in swimming. Furthermore, it suggests that regulation of receptor expression could underlie daily changes in behaviour as well as evolution of behaviour.

cAMP, Ca2+, pHi, and NO Regulate H-like Cation Channels That Underlie Feeding and Locomotion in the Predatory Sea Slug Pleurobranchaea californica

A systems approach to regulation of neuronal excitation in the mollusc Pleurobranchaea has described novel interactions of cyclic AMP-gated cation current (INa,cAMP), Ca2+, pHi, and NO. INa,cAMP appears in many neurons of feeding and locomotor neuronal networks. It is likely one of the family of hyperpolarization-activated, cyclic-nucleotide-gated currents (h-current) of vertebrate and invertebrate pacemaker networks. There are two isoforms. Ca2+ regulates both voltage dependence and depolarization-sensitive inactivation in both isoforms. The Type 1 INa,cAMP of the feeding network is enhanced by intracellular acidification. A direct dependence of INa,cAMP on cAMP allows the current to be used as a reporter on cAMP concentrations in the cell, and from there to the intrinsic activities of the synthetic adenyl cyclase and the degradative phosphodiesterase. Type 2 INa,cAMP of the locomotor system is activated by serotonergic inputs, while Type 1 of the feeding network is thought to be regulated peptidergically. NO synthase activity is high in the CNS, where it differs from standard neuronal NO synthase in not being Ca2+ sensitive. NO acidifies pHi, potentiating Type 1, and may act to open proton channels. A cGMP pathway does not mediate NO effects as in other systems. Rather, nitrosylation likely mediates its actions. An integrated model of the action of cAMP, Ca2+, pHi, and NO in the feeding network postulates that NO regulates proton conductance to cause neuronal excitation in the cell body on the one hand, and relief of activity-induced hyperacidification in fine dendritic processes on the other.

Sequences of Circadian Clock Proteins in the Nudibranch Molluscs Hermissenda crassicornis, Melibe leonina, and Tritonia diomedea

While much is known about the genes and proteins that make up the circadian clocks in vertebrates and several arthropod species, much less is known about the clock genes in many other invertebrates, including nudibranchs. The goal of this project was to identify the RNA and protein products of putative clock genes in the central nervous system of three nudibranchs, Hermissenda crassicornis, Melibe leonina, and Tritonia diomedea. Using previously published transcriptomes (Hermissenda and Tritonia) and a new transcriptome (Melibe), we identified nudibranch orthologs for the products of five canonical clock genes: brain and muscle aryl hydrocarbon receptor nuclear translocator like protein 1, circadian locomotor output cycles kaput, non-photoreceptive cryptochrome, period, and timeless. Additionally, orthologous sequences for the products of five related genes-aryl hydrocarbon receptor nuclear translocator like, photoreceptive cryptochrome, cryptochrome DASH, 6-4 photolyase, and timeout-were determined. Phylogenetic analyses confirmed that the nudibranch proteins were most closely related to known orthologs in related invertebrates, such as oysters and annelids. In general, the nudibranch clock proteins shared greater sequence similarity with Mus musculus orthologs than Drosophila melanogaster orthologs, which is consistent with the closer phylogenetic relationships recovered between lophotrochozoan and vertebrate orthologs. The suite of clock-related genes in nudibranchs includes both photoreceptive and non-photoreceptive cryptochromes, as well as timeout and possibly timeless. Therefore, the nudibranch clock may resemble the one exhibited in mammals, or possibly even in non-drosopholid insects and oysters. The latter would be evidence supporting this as the ancestral clock for bilaterians.

De novo assembly, annotation, and characterization of the whole brain transcriptome of male and female Syrian hamsters

Hamsters are an ideal animal model for a variety of biomedical research areas such as cancer, virology, circadian rhythms, and behavioural neuroscience. The use of hamsters has declined, however, most likely due to the dearth of genetic tools available for these animals. Our laboratory uses hamsters to study acute social stress, and we are beginning to investigate the genetic mechanisms subserving defeat-induced behavioural change. We have been limited, however, by the lack of genetic resources available for hamsters. In this study, we sequenced the brain transcriptome of male and female Syrian hamsters to generate the necessary resources to continue our research. We completed a de novo assembly and after assembly optimization, there were 113,329 transcripts representing 14,530 unique genes. This study is the first to characterize transcript expression in both female and male hamster brains and offers invaluable information to promote understanding of a host of important biomedical research questions for which hamsters are an excellent model.

Comparative transcriptomics enlarges the toolkit of known developmental genes in mollusks

BACKGROUND

Mollusks display a striking morphological disparity, including, among others, worm-like animals (the aplacophorans), snails and slugs, bivalves, and cephalopods. This phenotypic diversity renders them ideal for studies into animal evolution. Despite being one of the most species-rich phyla, molecular and in silico studies concerning specific key developmental gene families are still scarce, thus hampering deeper insights into the molecular machinery that governs the development and evolution of the various molluscan class-level taxa.

RESULTS

Next-generation sequencing was used to retrieve transcriptomes of representatives of seven out of the eight recent class-level taxa of mollusks. Similarity searches, phylogenetic inferences, and a detailed manual curation were used to identify and confirm the orthology of numerous molluscan Hox and ParaHox genes, which resulted in a comprehensive catalog that highlights the evolution of these genes in Mollusca and other metazoans. The identification of a specific molluscan motif in the Hox paralog group 5 and a lophotrochozoan ParaHox motif in the Gsx gene is described. Functional analyses using KEGG and GO tools enabled a detailed description of key developmental genes expressed in important pathways such as Hedgehog, Wnt, and Notch during development of the respective species. The KEGG analysis revealed Wnt8, Wnt11, and Wnt16 as Wnt genes hitherto not reported for mollusks, thereby enlarging the known Wnt complement of the phylum. In addition, novel Hedgehog (Hh)-related genes were identified in the gastropod Lottia cf. kogamogai, demonstrating a more complex gene content in this species than in other mollusks.

CONCLUSIONS

The use of de novo transcriptome assembly and well-designed in silico protocols proved to be a robust approach for surveying and mining large sequence data in a wide range of non-model mollusks. The data presented herein constitute only a small fraction of the information retrieved from the analysed molluscan transcriptomes, which can be promptly employed in the identification of novel genes and gene families, phylogenetic inferences, and other studies using molecular tools. As such, our study provides an important framework for understanding some of the underlying molecular mechanisms involved in molluscan body plan diversification and hints towards functions of key developmental genes in molluscan morphogenesis.

Deep sequencing of transcriptomes from the nervous systems of two decapod crustaceans to characterize genes important for neural circuit function and modulation

BACKGROUND

Crustaceans have been studied extensively as model systems for nervous system function from single neuron properties to behavior. However, lack of molecular sequence information and tools have slowed the adoption of these physiological systems as molecular model systems. In this study, we sequenced and performed de novo assembly for the nervous system transcriptomes of two decapod crustaceans: the Jonah crab (Cancer borealis) and the American lobster (Homarus americanus).

RESULTS

Forty-two thousand, seven hundred sixty-six and sixty thousand, two hundred seventy-three contigs were assembled from C. borealis and H. americanus respectively, representing 9,489 and 11,061 unique coding sequences. From these transcripts, genes associated with neural function were identified and manually curated to produce a characterization of multiple gene families important for nervous system function. This included genes for 34 distinct ion channel types, 17 biogenic amine and 5 GABA receptors, 28 major transmitter receptor subtypes including glutamate and acetylcholine receptors, and 6 gap junction proteins - the Innexins.

CONCLUSION

With this resource, crustacean model systems are better poised for incorporation of modern genomic and molecular biology technologies to further enhance the interrogation of fundamentals of nervous system function.

De novo transcriptome assembly of the marine gastropod Reishia clavigera for supporting toxic mechanism studies

The intertidal whelk Reishia clavigera is commonly used as a biomonitor of chemical contamination in the marine environment along Western Pacific region, and as a model for mechanistic studies of organotin-mediated imposex development. However, limited genomic resources of R. clavigera have restricted its role for the investigation of molecular mechanisms of such endocrine disruptions. This study, therefore, aimed to establish tissue-specific transcriptomes of the digestive gland, gonad, head ganglia, penis and the remaining body part of the male and female R. clavigera. By combining the results, a global transcriptome was obtained. A total of 578,134,720 high-quality filtered reads were obtained using Illumina sequencing. The R. clavigera transcriptome comprised of 38,466 transcripts and 32,798 unigenes with predicted open reading frames. The average length of transcripts was 1,709bp with N50 of 2,236bp. Based on sequence similarity searches against public databases, 28,657 transcripts and 24,403 unigenes had at least one BLAST hit. There were 17,530 transcripts and 14,897 unigenes annotated with at least one Gene Ontology (GO) term. Moreover, 5,776 transcripts and 5,137 unigenes were associated with 333 Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathways. The numbers of unigenes were similar among the five target tissues and between sexes, but tissue-specific expression profiles were revealed by multivariate analyses. Based on the functional annotation, putative steroid hormone-associated unigenes were identified. In particular, we highlighted the presence of steroid hormone receptor homologues that could be the targets for mechanistic studies of the organotin-mediated imposex development in marine gastropods. This newly generated transcriptome assembly of R. clavigera provides a valuable molecular resource for ecotoxicological and environmental genomic studies.

De novo assembly and annotation of the European abalone Haliotis tuberculata transcriptome

The European abalone Haliotis tuberculata is a delicacy and consequently a commercially valuable gastropod species. Aquaculture production and wild populations are subjected to multiple climate-associated stressors and anthropogenic pressures, including rising sea-surface temperatures, ocean acidification and an emerging pathogenic Vibrio infection. Transcript expression data provides a valuable resource for understanding abalone responses to variation in the biotic and abiotic environment. To generate an extensive transcriptome, we performed next-generation sequencing of RNA on larvae exposed to temperature and pH variation and on haemolymph of adults from two wild populations after experimental infection with Vibrio harveyi. We obtained more than 1.5 billion raw paired-end reads, which were assembled into 328,519 contigs. Filtration and clustering produced a transcriptome of 41,099 transcripts, of which 10,626 (25.85%) were annotated with Blast hits, and 7380 of these were annotated with Gene Ontology (GO) terms in Blast2Go. A differential expression analysis comparing all samples from the two life stages identified 5690 and 10,759 transcripts with significantly higher expression in larvae and adult haemolymph respectively. This is the greatest sequencing effort yet in the Haliotis genus, and provides the first high-throughput transcriptomic resource for H. tuberculata.

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.

Sequencing, De Novo Assembly, and Annotation of the Transcriptome of the Endangered Freshwater Pearl Bivalve, Cristaria plicata, Provides Novel Insights into Functional Genes and Marker Discovery

BACKGROUND

The freshwater mussel Cristaria plicata (Bivalvia: Eulamellibranchia: Unionidae), is an economically important species in molluscan aquaculture due to its use in pearl farming. The species have been listed as endangered in South Korea due to the loss of natural habitats caused by anthropogenic activities. The decreasing population and a lack of genomic information on the species is concerning for environmentalists and conservationists. In this study, we conducted a de novo transcriptome sequencing and annotation analysis of C. plicata using Illumina HiSeq 2500 next-generation sequencing (NGS) technology, the Trinity assembler, and bioinformatics databases to prepare a sustainable resource for the identification of candidate genes involved in immunity, defense, and reproduction.

RESULTS

The C. plicata transcriptome analysis included a total of 286,152,584 raw reads and 281,322,837 clean reads. The de novo assembly identified a total of 453,931 contigs and 374,794 non-redundant unigenes with average lengths of 731.2 and 737.1 bp, respectively. Furthermore, 100% coverage of C. plicata mitochondrial genes within two unigenes supported the quality of the assembler. In total, 84,274 unigenes showed homology to entries in at least one database, and 23,246 unigenes were allocated to one or more Gene Ontology (GO) terms. The most prominent GO biological process, cellular component, and molecular function categories (level 2) were cellular process, membrane, and binding, respectively. A total of 4,776 unigenes were mapped to 123 biological pathways in the KEGG database. Based on the GO terms and KEGG annotation, the unigenes were suggested to be involved in immunity, stress responses, sex-determination, and reproduction. A total of 17,251 cDNA simple sequence repeats (cSSRs) were identified from 61,141 unigenes (size of >1 kb) with the most abundant being dinucleotide repeats.

CONCLUSIONS

This dataset represents the first transcriptome analysis of the endangered mollusc, C. plicata. The transcriptome provides a comprehensive sequence resource for the conservation of genetic information in this species and enrichment of the genetic database. The development of molecular markers will assist in the genetic improvement of C. plicata.

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.