The transcriptome of the early life history stages of the California Sea Hare Aplysia californica

APGW/AKH Precursor from Rotifer Brachionus plicatilis and the DNA Loss Model Explain Evolutionary Trends of the Neuropeptide LWamide, APGWamide, RPCH, AKH, ACP, CRZ, and GnRH Families

In the year 2002, DNA loss model (DNA-LM) postulated that neuropeptide genes to emerged through codons loss via the repair of damaged DNA from ancestral gene namely Neuropeptide Precursor Predictive (NPP), which organization correspond two or more neuropeptides precursors evolutive related. The DNA-LM was elaborated according to amino acids homology among LWamide, APGWamide, red pigment-concentrating hormone (RPCH), adipokinetic hormones (AKHs) and in silico APGW/RPCH NPPAPGW/AKH NPP were proposed. With the above principle, it was proposed the evolution of corazonin (CRZ), gonadotropin-releasing hormone (GnRH), AKH, and AKH/CRZ (ACP), but any NPP never was considered. However, the evolutive relation via DNA-LM among these neuropeptides precursors not has been established yet. Therefore, the transcriptomes from crabs Callinectes toxotes and Callinectes arcuatus were used to characterized ACP and partial CRZ precursors, respectively. BLAST alignment with APGW/RPCH NPP and APGW/AKH NPP allow identified similar NPP in the rotifer Brachionus plicatilis and other invertebrates. Moreover, three bioinformatics algorithms and manual verification were used to purify 13,778 sequences, generating a database with 719 neuropeptide precursors. Phylogenetic trees with the DNA-LM parameters showed that some ACP, CRZ, AKH2 and two NPP share nodes with GnRH from vertebrates and some of this neuropeptide had nodes in invertebrates. Whereas the phylogenetic tree with standard parameters do not showed previous node pattern. Robinson-Foulds metric corroborates the differences among phylogenetic trees. Homology relationship showed four putative orthogroups; AKH4, CRZ, and protostomes GnRH had individual group. This is the first demonstration of NPP in species and would explain the evolution neuropeptide families by the DNA-LM.

Transcriptome Dynamics of an Oyster Larval Response to a Conspecific Cue-Mediated Settlement Induction in the Pacific Oyster Crassostrea gigas

The molecular mechanisms underlying the conspecific cue-mediated larval settlement in Crassostrea gigas is not yet fully understood. In this study, we described and compared the transcriptomes of competent pediveligers (Pedi) and conspecific cue-induced postlarvae (PL). A total of 2383 candidate transcripts were identified: 740 upregulated and 1643 downregulated transcripts, after settlement. Gene Ontology analysis revealed active chitin binding, calcium ion binding, and extracellular region processes in both stages. Results showed that the differential expression trend of six candidate transcripts were consistent between the quantitative real-time PCR and transcriptome data. The differential transcript expression related to shell formation showed closely linked dynamics with a gene regulatory network that may involve the interplay of various hormone receptors, neurotransmitters, and neuropeptide receptors working together in a concerted way in the Pedi and PL stages. Our results highlight the transcriptome dynamics underlying the settlement of oysters on conspecific adult shells and demonstrate the potential use of this cue as an attractant for wild and hatchery-grown oyster larval attachment on artificial substrates. It also suggests the possible involvement of an ecdysone signal pathway that may be linked to a neuroendocrine-biomineralization crosstalk in C. gigas settlement.

Proteome of larval metamorphosis induced by epinephrine in the Fujian oyster Crassostrea angulata

Background

The Fujian oyster Crassostrea angulata is an economically important species that has typical settlement and metamorphosis stages. The development of the oyster involves complex morphological and physiological changes, the molecular mechanisms of which are as yet unclear.

Results

In this study, changes in proteins were investigated during larval settlement and metamorphosis of Crassostrea angulata using epinephrine induction. Protein abundance and identity were characterized using label-free quantitative proteomics, tandem mass spectrometry (MS/ MS), and Mascot methods. The results showed that more than 50% (764 out of 1471) of the quantified proteins were characterized as differentially expressed. Notably, more than two-thirds of the differentially expressed proteins were down-regulated in epinephrine-induced larvae. The results showed that “metabolic process” was closely related to the development of settlement and metamorphosis; 5 × 10^− 4 M epinephrine induced direct metamorphosis of larvae and was non-toxic. Calmodulin and MAPK pathways were involved in the regulation of settlement of the oyster. Expression levels of immune-related proteins increased during metamorphosis. Hepatic lectin-like proteins, cadherins, calmodulin, calreticulin, and cytoskeletal proteins were involved in metamorphosis. The nervous system may be remodeled in larval metamorphosis induced by epinephrine. Expression levels of proteins that were enriched in the epinephrine signaling pathway may reflect the developmental stage of the larvae, that may reflect whether or not larvae were directly involved in metamorphosis when the larvae were treated with epinephrine.

Conclusion

The study provides insight into proteins that function in energy metabolism, immune responses, settlement and metamorphosis, and shell formation in C. angulata. The results contribute valuable information for further research on larval settlement and metamorphosis.

Graphical abstract

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.

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.

Proteomics Studies on the three Larval Stages of Development and Metamorphosis of Babylonia areolata

The ivory shell, Babylonia areolata, is a commercially important aquaculture species in the southeast coast of mainland China. The middle veliger stage, later veliger stage, and juvenile stage are distinct larval stages in B. areolata development. In this study, we used label-free quantification proteomics analysis of the three developmental stages of B. areolata. We identified a total of 5,583 proteins, of which 1,419 proteins expression level showed significant differential expression. The results of gene ontology enrichment analysis showed that the number of proteins involved in metabolic and cellular processes were the most abundant. Those proteins mostly had functions such as binding, catalytic activity and transporter activity. The results of Kyoto Encyclopedia of Genes and Genomes enrichment analysis showed that the number of proteins involved in the ribosome, carbon metabolism, and lysosome pathways were the most abundant, indicating that protein synthesis and the immune response were active during the three stages of development. This is the first study to use proteomics and real-time PCR to study the early developmental stages of B. areolata, which could provide relevant data on gastropod development. Our results provide insights into the novel aspects of protein function in shell formation, body torsion, changes in feeding habits, attachment and metamorphosis, immune-related activities in B. areolata larvae.

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.

Proteomic analysis of trochophore and veliger larvae development in the small abalone Haliotis diversicolor

Background

Haliotis diversicolor is commercially important species. The trochophore and veliger are distinct larval stages in gastropod development. Their development involves complex morphological and physiological changes. We studied protein changes during the embryonic development of H. diversicolor using two dimensional electrophoresis (2-DE) and label-free methods, tandem mass spectrometry (MS/ MS), and Mascot for protein identification.

Results

A total of 150 2-DE gel spots were identified. Protein spots showed upregulation of 15 proteins and downregulation of 28 proteins as H. diversicolor developed from trochophore to veliger larvae. Trochophore and veliger larvae were compared using a label-free quantitative proteomic approach. A total of 526 proteins were identified from both samples, and 104 proteins were differentially expressed (> 1.5 fold). Compared with trochophore larvae, veliger larvae had 55 proteins upregulated and 49 proteins downregulated. These differentially expressed proteins were involved in shell formation, energy metabolism, cellular and stress response processes, protein synthesis and folding, cell cycle, and cell fate determination. Compared with the 5 protein (fructose-bisphosphate aldolase, 14–3-3ε, profilin, actin-depolymerizing factor (ADF)/cofilin) and calreticulin) expression patterns, the mRNA expression exhibited similar patterns except gene of fructose-bisphosphate aldolase.

Conclusion

Our results provide insight into novel aspects of protein function in shell formation, torsion, and nervous system development, and muscle system differentiation in H. diversicolor larvae. “Quality control” proteins were identified to be involved in abalone larval development.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4203-7) contains supplementary material, which is available to authorized users.

Analysis of de novo sequencing and transcriptome assembly and lignocellulolytic enzymes gene expression of Coriolopsis gallica HTC

White-rot basidiomycete Coriolopsis gallica HTC is one of the main biodegraders of poplar. In our previous study, we have shown the strong capacity of C. gallica HTC to degrade lignocellulose. In this study, equal amounts of total RNA fromC. Gallica HTC cultures grown in different conditions were pooled together. Illumina paired-end RNA sequencing was performed, and 13.2 million 90-bp paired-end reads were generated. We chose the Merged Assembly of Oases data-set for the following blast searches and gene ontology analyses. The reads were assembled de novo into 28,034 transcripts (≥ 100 bp) using combined assembly strategy MAO. The transcripts were annotated using Blast2GO. In all, 18,810 transcripts (≥100 bp) achieved BLASTX hits, of which, 7048 transcripts had GO term and 2074 had ECs. The expression level of 11 lignocellulolytic enzyme genes from the assembled C. gallica HTC transcriptome were detected by real-time quantitative polymerase chain reaction. The results showed that expression levels of these genes were affected by carbon source and nitrogen source at the level of transcription. The current abundant transcriptome data allowed the identification of many new transcripts in C. gallica HTC. Data provided here represent the most comprehensive and integrated genomic resources for cloning and identifying genes of interest from C. gallica HTC. Characterization of C. gallica HTC transcriptome provides an effective tool to understand mechanisms underlying cellular and molecular functions of C. gallica HTC.

Differential Gene Expression during Larval Metamorphic Development in the Pearl Oyster, Pinctada fucata, Based on Transcriptome Analysis

P. fucata experiences a series of transformations in appearance, from swimming larvae to sessile juveniles, during which significant changes in gene expression likely occur. Thus, P. fucata could be an ideal model in which to study the molecular mechanisms of larval metamorphosis during development in invertebrates. To study the molecular driving force behind metamorphic development in larvae of P. fucata, transcriptomes of five larval stages (trochophore, D-shape, umbonal, eyespots, and spats) were sequenced using an Illumina HiSeq™ 2000 system and assembled and characterized with the transcripts of six tissues. As a result, a total of 174,126 unique transcripts were assembled and 60,999 were annotated. The number of unigenes varied among the five larval stages. Expression profiles were distinctly different between trochophore, D-shape, umbonal, eyespots, and spats larvae. As a result, 29 expression trends were sorted, of which eight were significant. Among others, 80 development-related, differentially expressed unigenes (DEGs) were identified, of which the majority were homeobox-containing genes. Most DEGs occurred among trochophore, D-shaped, and UES (umbonal, eyespots, and spats) larvae as verified by qPCR. Principal component analysis (PCA) also revealed significant differences in expression among trochophore, D-shaped, and UES larvae with ten transcripts identified but no matching annotations.

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.

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.

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

Deciphering life history transcriptomes in different environments

We compared whole transcriptome variation in six pre-adult stages and seven adult female ages in two populations of cactophilic Drosophila mojavensis reared on two host plants to understand how differences in gene expression influence standing life history variation. We used singular value decomposition (SVD) to identify dominant trajectories of life cycle gene expression variation, performed pairwise comparisons of stage and age differences in gene expression across the life cycle, identified when genes exhibited maximum levels of life cycle gene expression, and assessed population and host cactus effects on gene expression. Life cycle SVD analysis returned four significant components of transcriptional variation, revealing functional enrichment of genes responsible for growth, metabolic function, sensory perception, neural function, translation and ageing. Host cactus effects on female gene expression revealed population- and stage-specific differences, including significant host plant effects on larval metabolism and development, as well as adult neurotransmitter binding and courtship behaviour gene expression levels. In 3- to 6-day-old virgin females, significant upregulation of genes associated with meiosis and oogenesis was accompanied by downregulation of genes associated with somatic maintenance, evidence for a life history trade-off. The transcriptome of D. mojavensis reared in natural environments throughout its life cycle revealed core developmental transitions and genome-wide influences on life history variation in natural populations.

Assembly and annotation of a non-model gastropod (Nerita melanotragus) transcriptome: a comparison of de novo assemblers

Background

The sequencing, de novo assembly and annotation of transcriptome datasets generated with next generation sequencing (NGS) has enabled biologists to answer genomic questions in non-model species with unprecedented ease. Reliable and accurate de novo assembly and annotation of transcriptomes, however, is a critically important step for transcriptome assemblies generated from short read sequences. Typical benchmarks for assembly and annotation reliability have been performed with model species. To address the reliability and accuracy of de novo transcriptome assembly in non-model species, we generated an RNAseq dataset for an intertidal gastropod mollusc species, Nerita melanotragus, and compared the assembly produced by four different de novo transcriptome assemblers; Velvet, Oases, Geneious and Trinity, for a number of quality metrics and redundancy.

Results

Transcriptome sequencing on the Ion Torrent PGM™ produced 1,883,624 raw reads with a mean length of 133 base pairs (bp). Both the Trinity and Oases de novo assemblers produced the best assemblies based on all quality metrics including fewer contigs, increased N50 and average contig length and contigs of greater length. Overall the BLAST and annotation success of our assemblies was not high with only 15-19% of contigs assigned a putative function.

Conclusions

We believe that any improvement in annotation success of gastropod species will require more gastropod genome sequences, but in particular an increase in mollusc protein sequences in public databases. Overall, this paper demonstrates that reliable and accurate de novo transcriptome assemblies can be generated from short read sequencers with the right assembly algorithms.

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.

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.

Identification and expression of multiple CYP1-like and CYP3-like genes in the bivalve mollusk Mytilus edulis

Various sequencing projects over the last several years have aided the discovery of previously uncharacterized invertebrate sequences, including new cytochrome P450 genes (CYPs). Here we present data on the identification and characterization of two CYP1-like and three CYP3-like genes from the bivalve mollusk Mytilus edulis, and assess their potential as biomarkers based on their responses to several known vertebrate aryl hydrocarbon receptor (AHR) agonists. Quantitative real-time PCR was used to measure CYP transcript levels in digestive gland, labial palps, adductor muscle, gill, foot, and different regions of the mantle. Levels of both CYP1-like genes were highest in digestive gland, whereas labial palps had the highest expression levels of the three CYP3-like genes followed by digestive gland and outer margin of the mantle. Mussels were exposed by injection to the AHR agonists, β-naphthoflavone (BNF; 25 μg g(-1)), 3,3',4,4',5-polychlorinated biphenyl (PCB126; 2 μg g(-1)), or 6-formylindolo[3,2-b]carbazole (FICZ; 0.1 μg g(-1)), or to Aroclor 1254 (a mixture of PCBs; 50 μg g(-1)) for 24 h, followed by CYP expression analysis. There was no statistically significant change in expression of either of the CYP1-like genes after exposure to the various AHR agonists. The CYP3-like-1 gene was significantly up-regulated by BNF in gill tissues and the CYP3-like-2 gene was up-regulated in digestive gland by PCB126 and in gill tissue by BNF. These results suggest that distinct mechanisms of CYP gene activation could be present in M. edulis, although the importance of the CYP1-like and CYP3-like genes for xenobiotic and endogenous lipids biotransformation requires additional investigation.

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.

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.

De novo sequencing and transcriptome analysis of the central nervous system of mollusc Lymnaea stagnalis by deep RNA sequencing

The pond snail Lymnaea stagnalis is among several mollusc species that have been well investigated due to the simplicity of their nervous systems and large identifiable neurons. Nonetheless, despite the continued attention given to the physiological characteristics of its nervous system, the genetic information of the Lymnaea central nervous system (CNS) has not yet been fully explored. The absence of genetic information is a large disadvantage for transcriptome sequencing because it makes transcriptome assembly difficult. We here performed transcriptome sequencing for Lymnaea CNS using an Illumina Genome Analyzer IIx platform and obtained 81.9 M of 100 base pair (bp) single end reads. For de novo assembly, five programs were used: ABySS, Velvet, OASES, Trinity and Rnnotator. Based on a comparison of the assemblies, we chose the Rnnotator dataset for the following blast searches and gene ontology analyses. The present dataset, 116,355 contigs of Lymnaea transcriptome shotgun assembly (TSA), contained longer sequences and was much larger compared to the previously reported Lymnaea expression sequence tag (EST) established by classical Sanger sequencing. The TSA sequences were subjected to blast analyses against several protein databases and Aplysia EST data. The results demonstrated that about 20,000 sequences had significant similarity to the reported sequences using a cutoff value of 1e-6, and showed the lack of molluscan sequences in the public databases. The richness of the present TSA data allowed us to identify a large number of new transcripts in Lymnaea and molluscan species.

Developmental transcriptome of Aplysia californica

Genome-wide transcriptional changes in development provide important insight into mechanisms underlying growth, differentiation, and patterning. However, such large-scale developmental studies have been limited to a few representatives of Ecdysozoans and Chordates. Here, we characterize transcriptomes of embryonic, larval, and metamorphic development in the marine mollusc Aplysia californica and reveal novel molecular components associated with life history transitions. Specifically, we identify more than 20 signal peptides, putative hormones, and transcription factors in association with early development and metamorphic stages-many of which seem to be evolutionarily conserved elements of signal transduction pathways. We also characterize genes related to biomineralization-a critical process of molluscan development. In summary, our experiment provides the first large-scale survey of gene expression in mollusc development, and complements previous studies on the regulatory mechanisms underlying body plan patterning and the formation of larval and juvenile structures. This study serves as a resource for further functional annotation of transcripts and genes in Aplysia, specifically and molluscs in general. A comparison of the Aplysia developmental transcriptome with similar studies in the zebra fish Danio rerio, the fruit fly Drosophila melanogaster, the nematode Caenorhabditis elegans, and other studies on molluscs suggests an overall highly divergent pattern of gene regulatory mechanisms that are likely a consequence of the different developmental modes of these organisms.