Pyrosequencing-based transcriptomic resources in the pond snail Lymnaea stagnalis, with a focus on genes involved in molecular response to diquat-induced stress

A combined bioinformatics and LC-MS based approach for the development and benchmarking of a comprehensive database of Lymnaea CNS proteins

Applications of key technologies in biomedical research, such as qRT-PCR or LC-MS based proteomics, are generating large biological (-omics) data sets which are useful for the identification and quantification of biomarkers in any research area of interest. Genome, transcriptome and proteome databases are already available for a number of model organisms including vertebrates and invertebrates. However, there is insufficient information available for protein sequences of certain invertebrates, such as the great pond snail Lymnaea stagnalis, a model organism that has been used highly successfully in elucidating evolutionarily conserved mechanisms of memory function and dysfunction. Here we used a bioinformatics approach to designing and benchmarking a comprehensive CNS proteomics database (LymCNS-PDB) for the identification of proteins from the CNS of Lymnaea by LC-MS based proteomics. LymCNS-PDB was created by using the Trinity TransDecoder bioinformatics tool to translate amino acid sequences from mRNA transcript assemblies obtained from a published Lymnaea transcriptomics database. The blast-style MMSeq2 software was used to match all translated sequences to UniProtKB sequences for molluscan proteins, including Lymnaea and other molluscs. LymCNS-PDB contains 9,628 identified matched proteins that were benchmarked by performing LC-MS based proteomics analysis with proteins isolated from the Lymnaea CNS. MS/MS analysis using the LymCNS-PDB database led to the identification of 3,810 proteins. Only 982 proteins were identified by using a non-specific molluscan database. LymCNS-PDB provides a valuable tool that will enable us to perform quantitative proteomics analysis of protein interactomes involved in several CNS functions in Lymnaea, including learning and memory and age-related memory decline.

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.

Neuropeptide Localization in Lymnaea stagnalis: From the Central Nervous System to Subcellular Compartments

Due to the relatively small number of neurons (few tens of thousands), the well-established multipurpose model organism Lymnaea stagnalis, great pond snail, has been extensively used to study the functioning of the nervous system. Unlike the more complex brains of higher organisms, L. stagnalis has a relatively simple central nervous system (CNS) with well-defined circuits (e.g., feeding, locomotion, learning, and memory) and identified individual neurons (e.g., cerebral giant cell, CGC), which generate behavioral patterns. Accumulating information from electrophysiological experiments maps the network of neuronal connections and the neuronal circuits responsible for basic life functions. Chemical signaling between synaptic-coupled neurons is underpinned by neurotransmitters and neuropeptides. This review looks at the rapidly expanding contributions of mass spectrometry (MS) to neuropeptide discovery and identification at different granularity of CNS organization. Abundances and distributions of neuropeptides in the whole CNS, eleven interconnected ganglia, neuronal clusters, single neurons, and subcellular compartments are captured by MS imaging and single cell analysis techniques. Combining neuropeptide expression and electrophysiological data, and aided by genomic and transcriptomic information, the molecular basis of CNS-controlled biological functions is increasingly revealed.

Molecular Pathways and Pigments Underlying the Colors of the Pearl Oyster Pinctada margaritifera var. cumingii (Linnaeus 1758)

The shell color of the Mollusca has attracted naturalists and collectors for hundreds of years, while the molecular pathways regulating pigment production and the pigments themselves remain poorly described. In this study, our aim was to identify the main pigments and their molecular pathways in the pearl oyster Pinctada margaritifera—the species displaying the broadest range of colors. Three inner shell colors were investigated—red, yellow, and green. To maximize phenotypic homogeneity, a controlled population approach combined with common garden conditioning was used. Comparative analysis of transcriptomes (RNA-seq) of P. margaritifera with different shell colors revealed the central role of the heme pathway, which is involved in the production of red (uroporphyrin and derivates), yellow (bilirubin), and green (biliverdin and cobalamin forms) pigments. In addition, the Raper–Mason, and purine metabolism pathways were shown to produce yellow pigments (pheomelanin and xanthine) and the black pigment eumelanin. The presence of these pigments in pigmented shell was validated by Raman spectroscopy. This method also highlighted that all the identified pathways and pigments are expressed ubiquitously and that the dominant color of the shell is due to the preferential expression of one pathway compared with another. These pathways could likely be extrapolated to many other organisms presenting broad chromatic variation.

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 pond snail Lymnaea stagnalis

The freshwater snail Lymnaea stagnalis has a long research history, but only relatively recently has it emerged as an attractive model organism to study molecular mechanisms in the areas of developmental biology and translational medicine such as learning/memory and neurodegenerative diseases. The species has the advantage of being a hermaphrodite and can both cross- and self-mate, which greatly facilitates genetic approaches. The establishment of body-handedness, or chiromorphogenesis, is a major topic of study, since chirality is evident in the shell coiling. Chirality is maternally inherited, and only recently a gene-editing approach identified the actin-related gene Lsdia1 as the key handedness determinant. This short article reviews the natural habitat, life cycle, major research questions and interests, and experimental approaches.

Metabolic profiles of Oncomelania hupensis after molluscicidal treatment: Carbohydrate metabolism targeted and energy deficiency

Oncomelania hupensis is the intermediate host of Schistosoma japonicum, one of the Schistosoma species that can cause human schistosomiasis. Molluscicidal treatment remains the primary means to control snail. Niclosamide is the only molluscicide recommended by the World Health Organization, and it has been used throughout schistosomiasis-endemic areas in China for almost 30 years. In our previous studies on transcriptomics, morphology, and enzymology of snails after molluscicidal treatment, two effective molluscicides were used, 50% wettable powder of niclosamide ethanolamine salt (WPN) and a new molluscicide derived from niclosamide, the salt of quinoid-2', 5-dichloro-4'-nitro-salicylanilide (LDS, simplified for Liu Dai Shui Yang An). Genes involved in cell structure mintenance, inhibition of neurohumoral transmission, and energy metabolism showed significant differential expression after molluscicide treatments. Damages in the structure of liver and muscle cells were accompanied by inhibited activities of enzymes related to carbohydrate metabolism and energy supply. This study was designed to clarify the dynamic metabolic process by metabonomics, together with the previous transcriptomic and enzymological profiles, to identify potential metabolite markers and metabolism pathways that related to the toxic mechanism of the molluscicide. In total, 56 metabolites were identified for O. hupensis, and 75% of these metabolites consisted of amino acids and derivatives, organic acids, and nucleic acid components. The concentration of glucose, maltose, succinate, choline, and alanine changed significantly after molluscicide treatments. These changes in metabolites mainly occurred in the process of carbohydrate metabolism, energy metabolism, and amino acid metabolism, primarily related to glycolysis/gluconeogenesis, oxidative phosphorylation, and transamination by KEGG pathway identification. Most of the identified pathways were also related to those differentially expressed unigenes and observed enzymes from our previous studies. Inhibited aerobic respiration and oxidative phosphorylation, and energy deficiency were implied further to be the leading causes of the final death of snails after molluscicide treatments. The hypothesised mathematical model in this study identified the rational hysteresis to explain the inconsistency of responses of unigenes, enzymes, and metabolites to molluscicide treatments. This study contributes to the comprehensive understanding of the molluscicidal mechanism in the metabolic process and this could assist in improving existing molluscicide formulations or development of new molluscicides.

Functional shell matrix proteins tentatively identified by asymmetric snail shell morphology

Molluscan shell matrix proteins (SMPs) are essential in biomineralization. Here, we identify potentially important SMPs by exploiting the asymmetric shell growth in snail, Lymnaea stagnalis. Asymmetric shells require bilaterally asymmetric expression of SMP genes. We examined expression levels of 35,951 transcripts expressed in the left and right sides of mantle tissue of the pond snail, Lymnaea stagnalis. This transcriptome dataset was used to identify 207 SMPs by LC-MS/MS. 32 of the 207 SMP genes show asymmetric expression patterns, which were further verified for 4 of the 32 SMPs using quantitative PCR analysis. Among asymmetrically expressed SMPs in dextral snails, those that are more highly expressed on the left side than the right side are 3 times more abundant than those that are more highly expressed on the right than the left, suggesting potentially inhibitory roles of SMPs in shell formation. The 32 SMPs thus identified have distinctive features, such as conserved domains and low complexity regions, which may be essential in biomineralization.

The unlimited potential of the great pond snail, Lymnaea stagnalis

Only a limited number of animal species lend themselves to becoming model organisms in multiple biological disciplines: one of these is the great pond snail, Lymnaea stagnalis. Extensively used since the 1970s to study fundamental mechanisms in neurobiology, the value of this freshwater snail has been also recognised in fields as diverse as host-parasite interactions, ecotoxicology, evolution, genome editing and 'omics', and human disease modelling. While there is knowledge about the natural history of this species, what is currently lacking is an integration of findings from the laboratory and the field. With this in mind, this article aims to summarise the applicability of L. stagnalis and points out that this multipurpose model organism is an excellent, contemporary choice for addressing a large range of different biological questions, problems and phenomena.

Physiological and pharmacological characterization of a molluscan neuronal efflux transporter; evidence for age-related transporter impairment

Plasma membrane efflux transporters play crucial roles in the removal and release of both harmful and beneficial substances from the interior of cells and tissue types in virtually every extant species. They contribute to the clearance of a broad spectrum of exogenous and endogenous toxicants and harmful metabolites, including the reactive lipid aldehyde byproducts of lipid peroxidation that are a hallmark of cellular ageing. Here, we tested whether declining transporter functionality may contribute to functional decline in a snail model of neuronal ageing. Through measuring the removal of 5(6)-carboxyfluorescein, a known substrate for membrane efflux transporters, we provide, for the first time, physiological evidence for the existence of probenecid-, MK571- and glutathione-sensitive efflux transporters in (gastropod) neurons and demonstrate that their functionality declines with age. Our data support the idea that waning cellular detoxification capacity might be a significant factor in the escalation of (lipo-)toxicity observed in neuronal ageing.

Effects on tadpole snail gene expression after exposure to vinclozolin

The fungicide vinclozolin (Vz) is an endocrine disruptor with known anti-androgenic activity in vertebrates. However, there is a lack of information about the Vz mode of action in invertebrates, although some studies have shown that this compound can produce alterations in different species. Transcriptional activity was analyzed in the freshwater snail Physella acuta in order to elucidate putative cellular processes altered by this chemical during a response. In order to identify potential molecular biomarkers, a de novo transcriptome was generated for this species that constitutes a valuable source for future studies. This data, together with some already available data, permitted the identification of several genes related to detoxification mechanisms (Cyp2u1, Cyp3a7, Cyp4f22, GSTo1, GSTt2, and MRP1), stress response (Hsp20.4, Hsp17, Hsp16.6, and Cu,Zn-SOD), the hormonal system (Estrogen Receptor and Hsp90), apoptosis (Casp3), and copper homeostasis (ATOX1). Using quantitative Real-Time polymerase chain reaction, mRNA levels of these genes were examined in snails exposed to 20 or 200 µg/L Vz for 24 h. The results showed an overall weak response, with downregulation of Hsp20.4 and no statistically significant change for the other genes. These findings suggest that P. acuta can manage the concentrations of Vz found in the environment with no relevant activation of the pathways analyzed, although additional studies are needed for longer exposure times and including other metabolic pathways. The new genes described open the range of processes that can be studied at the molecular level in toxicity tests.

Combining the assessment of apical endpoints and gene expression in the freshwater snail Physa acuta after exposure to reclaimed water

Post-treatment wastewater reuses are diverse. Recreational and environmental restoration uses of reclaimed water (RW) can be potentially harmful to aquatic organisms. In this work the freshwater snail Physa acuta was exposed to RW (100%) and its dilution (RW 50%). A simple laboratory mixture of three emerging pollutants was used to address the complex problem of mixture toxicity of RW. Hence fortified reclaimed water (FRW), obtained by adding fluoxetine (400 μg FLX/L), perfluorooctane sulphonic acid (90 μg PFOS/L) and methylparaben (9 μg MP/L), was tested at two dilution percentages: 100% and 50%. The effects of the laboratory mixture of FLX, PFOS and MP on the test medium were also studied. Long-lasting effects, together with early molecular responses, were assessed. Fecundity (cumulative egg production) over 21 days and the hatching of produced eggs (F1) after another 21-day embryonic exposure were monitored. The gene expression of three genes was analysed after 24 h of exposure: two endocrine-related nuclear receptors (ERR and RXR) and one stress protein gene (Hsp70). This reproduction test, with additional assessments of the F1 recovered eggs' hatching success, showed that both RW and FRW significantly reduced fecundity. F1 hatching was affected only by FRW. The gene expression results showed that the RXR response was strikingly similar to the fecundity response, which suggests that this nuclear receptor is involved in the reproductive pathways of gastropods. ERR remained virtually unaltered. Hsp70 was overexpressed by the laboratory mixture in the test medium, but no effect was observed in the fortification of RW. This opposite effect and lack of response for F1 hatching produced by the laboratory mixture in the test medium highlighted the difficulty of predicting mixture effects. The experimental approach allowed us to test the effects caused by RW on P. acuta at different biological organisation levels. Thus, the combination of molecular biomarkers and ecological relevant endpoints is a good strategy to test complex mixtures like RW as it provides a framework to link mechanisms of action and whole organism effects when it is almost impossible to detect the pollutant(s) that cause toxic effects.

Disturbances in reproduction and expression of steroidogenic enzymes in aquatic invertebrates exposed to components of the herbicide Roundup

Exposure of organisms to environmental contaminants is a growing concern. We have investigated the effects of the individual active ingredients of the herbicide Roundup (glyphosate and diquat dibromide [DD]) since Roundup causes alterations in reproduction, mortality, and development in the aquatic snail Lymnaea palustris. Snails chronically treated with elevated but ecologically relevant levels of DD exhibit reduction in fecundity ( p < 0.05), while fecundity in glyphosate-treated snails is comparable to or exceeds control levels. To investigate a possible mechanism for the reproductive disturbance, we monitored levels of steroid acute regulatory (StAR) protein in whole snails and observed a correlation in StAR protein decrease with treatment with Roundup, glyphosate, or DD. We detect StAR in organs where steroid biosynthesis occurs (ovotestis, brain, kidney); StAR protein is reduced following chronic exposure to Roundup, glyphosate, or DD ( p < 0.01). Estradiol and testosterone concentrations in hemolymph were measured by enzyme-linked immunosorbent assay following 3-week exposure of snails to 3.5 mg/L glyphosate or 140 µg/L DD. Testosterone levels decrease in DD-treated groups ( p < 0.05); a trend of lower testosterone is also observed in glyphosate-treated groups ( p > 0.05). Estradiol concentration is greater than or equal to control levels in glyphosate, but decreased in DD ( p < 0.05). Because of its role in the conversion of testosterone to estradiol, we monitored abundance of aromatase and observed a reduction ( p < 0.05) in DD-treated snails (consistent with the drop in fecundity and estradiol levels) and a comparable level to control in glyphosate-treated snails (consistent with their high fecundity and estradiol levels). Although the toxicity of commercially-available Roundup to aquatic animals may have many contributing factors including its inactive surfactant, the constituent of Roundup associated with the greatest reproductive disturbances and observed developmental abnormalities of offspring is DD. This study details the analysis of particular herbicide constituents and their effect on specific targets in the reproductive pathway.

A full life-cycle bioassay with Cantareus aspersus shows reproductive effects of a glyphosate-based herbicide suggesting potential endocrine disruption

A full life-cycle (240 days) bioassay using the terrestrial snail, Cantareus aspersus, allowing exposure during embryogenesis and/or the growth and reproduction phases, was used to assess the effects of Bypass®, a glyphosate-based herbicide (GlyBH), on a range of endpoints, including parameters under endocrine control. As a positive control, a mixture (R-A) made of diquat (Reglone®) and nonylphenols (NP, Agral®), known for its endocrine disrupting effects in other organisms, was tested. At environmental concentrations, both pesticides (R-A mixture and GlyBH) enhanced growth but reduced reproduction. The R-A mixture acted mainly on the fecundity through a delay in egg-laying of approximately 20 days and a strongly reduced number of clutches. This latter dysfunction may be caused by a permanent eversion of the penis, suggesting a disrupting effect at the neuro-endocrine level, which prevented normal mating. GlyBH acted on fertility, possibly due to a decrease in the fertilization of eggs laid by adults exposed during their embryonic development. These results, associated with the absence of observed effects on gonad histology of GlyBH exposed snails, suggested that the underlying mechanisms are neuro-endocrine.

The mechanisms of nickel toxicity in aquatic environments: An adverse outcome pathway analysis

Current ecological risk assessment and water quality regulations for nickel (Ni) use mechanistically based, predictive tools such as biotic ligand models (BLMs). However, despite many detailed studies, the precise mechanism(s) of Ni toxicity to aquatic organisms remains elusive. This uncertainty in the mechanism(s) of action for Ni has led to concern over the use of tools like the BLM in some regulatory settings. To address this knowledge gap, the authors used an adverse outcome pathway (AOP) analysis, the first AOP for a metal, to identify multiple potential mechanisms of Ni toxicity and their interactions with freshwater aquatic organisms. The analysis considered potential mechanisms of action based on data from a wide range of organisms in aquatic and terrestrial environments on the premise that molecular initiating events for an essential metal would potentially be conserved across taxa. Through this analysis the authors identified 5 potential molecular initiating events by which Ni may exert toxicity on aquatic organisms: disruption of Ca²⁺ homeostasis, disruption of Mg²⁺ homeostasis, disruption of Fe^2+/3+ homeostasis, reactive oxygen species-induced oxidative damage, and an allergic-type response of respiratory epithelia. At the organ level of biological organization, these 5 potential molecular initiating events collapse into 3 potential pathways: reduced Ca²⁺ availability to support formation of exoskeleton, shell, and bone for growth; impaired respiration; and cytotoxicity and tumor formation. At the level of the whole organism, the organ-level responses contribute to potential reductions in growth and reproduction and/or alterations in energy metabolism, with several potential feedback loops between each of the pathways. Overall, the present AOP analysis provides a robust framework for future directed studies on the mechanisms of Ni toxicity and for developing AOPs for other metals. Environ Toxicol Chem 2017;36:1128-1137. © 2016 SETAC.

De novo transcriptome sequencing and analysis of freshwater snail (Radix balthica) to discover genes and pathways affected by exposure to oxazepam

Pharmaceuticals are increasingly found in aquatic ecosystems due to the non-efficiency of waste water treatment plants. Therefore, aquatic organisms are frequently exposed to a broad diversity of pharmaceuticals. Freshwater snail Radix balthica has been chosen as model to study the effects of oxazepam (psychotropic drug) on developmental stages ranging from trochophore to hatching. In order to provide a global insight of these effects, a transcriptome deep sequencing has been performed on exposed embryos. Eighteen libraries were sequenced, six libraries for three conditions: control, exposed to the lowest oxazepam concentration with a phenotypic effect (delayed hatching) (TA) and exposed to oxazepam concentration found in freshwater (TB). A total of 39,759,772 filtered raw reads were assembled into 56,435 contigs having a mean length of 1579.68 bp and mean depth of 378.96 reads. 44.91% of the contigs have at least one annotation. The differential expression analysis between the control condition and the two exposure conditions revealed 146 contigs differentially expressed of which 144 for TA and two for TB. 34.0% were annotated with biological function. There were four mainly impacted processes: two cellular signalling systems (Notch and JNK) and two biosynthesis pathways (Polyamine and Catecholamine pathways). This work reports a large-scale analysis of differentially transcribed genes of R. balthica exposed to oxazepam during egg development until hatching. In addition, these results enriched the de novo database of potential ecotoxicological models.

The Shell of the Invasive Bivalve Species Dreissena polymorpha: Biochemical, Elemental and Textural Investigations

The zebra mussel Dreissena polymorpha is a well-established invasive model organism. Although extensively used in environmental sciences, virtually nothing is known of the molecular process of its shell calcification. By describing the microstructure, geochemistry and biochemistry/proteomics of the shell, the present study aims at promoting this species as a model organism in biomineralization studies, in order to establish a bridge with ecotoxicology, while sketching evolutionary conclusions. The shell of D. polymorpha exhibits the classical crossed-lamellar/complex crossed lamellar combination found in several heterodont bivalves, in addition to an external thin layer, the characteristics of which differ from what was described in earlier publication. We show that the shell selectively concentrates some heavy metals, in particular uranium, which predisposes D. polymorpha to local bioremediation of this pollutant. We establish the biochemical signature of the shell matrix, demonstrating that it interacts with the in vitro precipitation of calcium carbonate and inhibits calcium carbonate crystal formation, but these two properties are not strongly expressed. This matrix, although overall weakly glycosylated, contains a set of putatively calcium-binding proteins and a set of acidic sulphated proteins. 2D-gels reveal more than fifty proteins, twenty of which we identify by MS-MS analysis. We tentatively link the shell protein profile of D. polymorpha and the peculiar recent evolution of this invasive species of Ponto-Caspian origin, which has spread all across Europe in the last three centuries.

Metabolomics to Explore Imidacloprid-Induced Toxicity in the Central Nervous System of the Freshwater Snail Lymnaea stagnalis

Modern toxicology is seeking new testing methods to better understand toxicological effects. One of the most concerning chemicals is the neonicotinoid pesticide imidacloprid. Although imidacloprid is designed to target insects, recent studies have shown adverse effects on nontarget species. Metabolomics was applied to investigate imidacloprid-induced sublethal toxicity in the central nervous system of the freshwater snail Lymnaea stagnalis. The snails (n = 10 snails) were exposed for 10 days to increasing imidacloprid concentrations (0.1, 1, 10, and 100 μg/L). The comparison between control and exposure groups highlighted the involvement and perturbation of many biological pathways. The levels of several metabolites belonging to different metabolite classes were significantly changed by imidacloprid exposure. A change in the amino acids and nucleotide metabolites like tryptophan, proline, phenylalanine, uridine, and guanosine was found. Many fatty acids were down-regulated, and the levels of the polyamines, spermidine and putrescine, were found to be increased which is an indication of neuron cell injury. A turnover increase between choline and acetylcholine led us to hypothesize an increase in cholinergic gene expression to overcome imidacloprid binding to the nicotinic acetylcholine receptors. Metabolomics revealed imidacloprid induced metabolic changes at low and environmentally relevant concentration in a nontarget species and generated a novel mechanistic hypothesis.

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.

Next-generation proteomics: toward customized biomarkers for environmental biomonitoring

Because of their ecological representativeness, invertebrates are commonly employed as test organisms in ecotoxicological assessment; however, to date, biomarkers employed for these species were the result of a direct transposition from vertebrates, despite deep evolutionary divergence. To gain efficiency in the diagnostics of ecosystem health, specific biomarkers must be developed. In this sense, next-generation proteomics enables the specific identification of proteins involved in key physiological functions or defense mechanisms, which are responsive to ecotoxicological challenges. However, the analytical investment required restricts use in biomarker discovery. Routine biomarker validation and assays rely on more conventional mass spectrometers. Here, we describe how proteomics remains a challenge for ecotoxicological test organisms because of the lack of appropriate protein sequences databases, thus restricting the analysis on conserved and ubiquitous proteins. These limits and some strategies used to overcome them are discussed. These new tools, such as proteogenomics and targeted proteomics, should result in new biomarkers specific to relevant environmental organisms and applicable to routine ecotoxicological assessment.

Comparative survey of bacterial and archaeal communities in high arsenic shallow aquifers using 454 pyrosequencing and traditional methods

A survey of bacterial and archaeal community structure was carried out in 10 shallow tube wells in a high arsenic groundwater system located in Hetao Basin, Inner Mongolia by 16S rRNA gene based two-step nested PCR-DGGE, clone libraries and 454 pyrosequencing. 12 bacterial and 18 archaeal DGGE bands and 26-136 species-level OTUs were detected for all the samples. 299 bacterial and 283 archaeal 16S rRNA gene clones for two typical samples were identified by phylogenetic analysis. Most of the results from these different methods were consistent with the dominant bacterial populations. But the proportions of the microbial populations were mostly different and the bacterial communities in most of these samples from pyrosequencing were both more abundant and more diverse than those from the traditional methods. Even after quality filtering, pyrosequencing revealed some populations including Alishewanella, Sulfuricurvum, Arthrobacter, Sporosarcina and Algoriphagus which were not detected with traditional techniques. The most dominant bacterial populations in these samples identified as some arsenic, iron, nitrogen and sulfur reducing and oxidizing related populations including Acinetobacter, Pseudomonas, Flavobacterium, Brevundimonas, Massilia, Planococcus, and Aquabacterium and archaeal communities Nitrosophaera and Methanosaeta. Acinetobacter and Pseudomonas were distinctly abundant in most of these samples. Methanogens were found as the dominant archeal population with three methods. From the results of traditional methods, the dominant archaeal populations apparently changed from phylum Thaumarchaeota to Euryarchaeota with the arsenic concentrations increasing. But this structure dynamic change was not revealed with pyrosequencing. Our results imply that an integrated approach combining the traditional methods and next generation sequencing approaches to characterize the microbial communities in high arsenic groundwater is recommended.

Proteomic Analysis of the Reproductive Organs of the Hermaphroditic Gastropod Lymnaea stagnalis Exposed to Different Endocrine Disrupting Chemicals

Many studies have reported perturbations of mollusc reproduction following exposure to low concentrations (ng/L range) of endocrine disrupting chemicals (EDCs). However, the mechanisms of action of these molecules on molluscs are still poorly understood. Investigation of the modifications of protein expression in organisms exposed to chemicals using proteomic methods can provide a broader and more comprehensive understanding of adverse impacts of pollution on organisms than conventional biochemical biomarkers (e.g., heat-shock proteins, metallothioneins, GST, EROD). In this study we have investigated the impacts of four chemicals, which exhibit different endocrine disrupting properties in vertebrates, on the proteome of the hermaphroditic freshwater pulmonate gastropod Lymnaea stagnalis after 21 days of exposure. Testosterone, tributyltin, chlordecone and cyproterone acetate were chosen as tested compounds as they can induce adverse effects on the reproduction of this snail. The 2D-DIGE method was used to identify proteins whose expression was affected by these compounds. In addition to modifying the expression of proteins involved in the structure and function of the cytoskeleton, chemicals had impacts on the expression of proteins involved in the reproduction of L. stagnalis. Exposure to 19.2 µg/L of chlordecone increased the abundance of ovipostatin, a peptide transmitted during mating through seminal fluid, which reduces oviposition in this species. The expression of yolk ferritin, the vitellogenin equivalent in L. stagnalis, was reduced after exposure to 94.2 ng Sn/L of tributyltin. The identification of yolk ferritin and the modification of its expression in snails exposed to chemicals were refined using western blot analysis. Our results showed that the tested compounds influenced the abundance of yolk ferritin in the reproductive organs. Alteration in proteins involved in reproductive pathways (e.g., ovipostatin and yolk ferritin) could constitute relevant evidence of interaction of EDCs with reproductive pathways that are under the control of the endocrine system of L. stagnalis.

afterParty: turning raw transcriptomes into permanent resources

Background

Next-generation DNA sequencing technologies have made it possible to generate transcriptome data for novel organisms quickly and cheaply, to the extent that the effort required to annotate and publish a new transcriptome is greater than the effort required to sequence it. Often, following publication, details of the annotation effort are only available in summary form, hindering subsequent exploitation of the data. To promote best-practice in annotation and to ensure that data remain accessible, we have written afterParty, a web application that allows users to assemble, annotate and publish novel transcriptomes using only a web browser.

Results

afterParty is a robust web application that implements best-practice transcriptome assembly, annotation, browsing, searching, and visualization. Users can turn a collection of reads (from Roche 454 chemistry) or assembled contigs (from any sequencing chemistry, including Illumina Solexa RNA-Seq) into a searchable, browsable transcriptome resource and quickly make it publicly available. Contigs are functionally annotated based on similarity to known sequences and protein domains. Once assembled and annotated, transcriptomes derived from multiple species or libraries can be compared and searched. afterParty datasets can either be created using the existing afterParty server, or using local instances that can be built easily using a virtual machine. afterParty includes powerful visualization tools for transcriptome dataset exploration and uses a flexible annotation architecture which will allow additional types of annotation to be added in the future.

Conclusions

afterParty's main use case scenario is one in which a working biologist has generated a large volume of transcribed sequence data and wishes to turn it into a useful resource that has some durability. By reducing the effort, bioinformatics skills, and computational resources needed to annotate and publish a transcriptome, afterParty will facilitate the annotation and sharing of sequence data that would otherwise remain unavailable. A typical metazoan transcriptome containing several tens of thousands of contigs can be annotated in a few minutes of interactive time and a few days of computational time.

Analysis of the Elodea nuttallii transcriptome in response to mercury and cadmium pollution: development of sensitive tools for rapid ecotoxicological testing

Toxic metals polluting aquatic ecosystems are taken up by inhabitants and accumulate in the food web, affecting species at all trophic levels. It is therefore important to have good tools to assess the level of risk represented by toxic metals in the environment. Macrophytes are potential organisms for the identification of metal-responsive biomarkers but are still underrepresented in ecotoxicology. In the present study, we used next-generation sequencing to investigate the transcriptomic response of Elodea nuttallii exposed to enhanced concentrations of Hg and Cd. We de novo assembled more than 60 000 contigs, of which we found 170 to be regulated dose-dependently by Hg and 212 by Cd. Functional analysis showed that these genes were notably related to energy and metal homeostasis. Expression analysis using nCounter of a subset of genes showed that the gene expression pattern was able to assess toxic metal exposure in complex environmental samples and was more sensitive than other end points (e.g., bioaccumulation, photosynthesis, etc.). In conclusion, we demonstrate the feasibility of using gene expression signatures for the assessment of environmental contamination, using an organism without previous genetic information. This is of interest to ecotoxicology in a wider sense given the possibility to develop specific and sensitive bioassays.

Impact of the redox-cycling herbicide diquat on transcript expression and antioxidant enzymatic activities of the freshwater snail Lymnaea stagnalis

The presence of pesticides in the environment results in potential unwanted effects on non-target species. Freshwater organisms inhabiting water bodies adjacent to agricultural areas, such as ditches, ponds and marshes, are good models to test such effects as various pesticides may reach these habitats through several ways, including aerial drift, run-off, and drainage. Diquat is a non-selective herbicide used for crop protection or for weed control in such water bodies. In this study, we investigated the effects of diquat on a widely spread aquatic invertebrate, the holarctic freshwater snail Lymnaea stagnalis. Due to the known redox-cycling properties of diquat, we studied transcript expression and enzymatic activities relative to oxidative and general stress in the haemolymph and gonado-digestive complex (GDC). As diquat is not persistent, snails were exposed for short times (5, 24, and 48 h) to ecologically relevant concentrations (22.2, 44.4, and 222.2 μg l(-1)) of diquat dibromide. RT-qPCR was used to quantify the transcription of genes encoding catalase (cat), a cytosolic superoxide dismutase (Cu/Zn-sod), a selenium-dependent glutathione peroxidase (gpx), a glutathione reductase (gred), the retinoid X receptor (rxr), two heat shock proteins (hsp40 and hsp70), cortactin (cor) and the two ribosomal genes r18S and r28s. Enzymatic activities of SOD, Gpx, Gred and glutathione S-transferase (GST) were investigated in the GDC using spectrophoto/fluorometric methods. Opposite trends were obtained in the haemolymph depending on the herbicide concentration. At the lowest concentration, effects were mainly observed after 24 h of exposure, with over-transcription of cor, hsp40, rxr, and sod, whereas higher concentrations down-regulated the expression of most of the studied transcripts, especially after 48 h of exposure. In the GDC, earlier responses were observed and the fold-change magnitude was generally much higher: transcription of all target genes increased significantly (or non-significantly for cat) after 5 h of exposure, and went back to control levels afterwards, suggesting the onset of an early response to oxidative stress associated to the unbalance of reactive oxygen species (ROS) in hepatocytes. Although increases obtained for Gred and SOD activities were globally consistent with their respective transcript expressions, up-regulation of transcription was not always correlated with increase of enzymatic activity, indicating that diquat might affect steps downstream of transcription. However, constitutive levels of enzymatic activities were at least maintained. In conclusion, diquat was shown to affect expression of the whole set of studied transcripts, reflecting their suitability as markers of early response to oxidative stress in L. stagnalis.

Characterization of transcriptomes from sexual and asexual lineages of a New Zealand snail (Potamopyrgus antipodarum)

Understanding the evolution and maintenance of sexual reproduction is one of the central challenges of evolutionary biology, yet we know very little about how sex influences molecular evolution. The New Zealand freshwater snail Potamopyrgus antipodarum is ideally suited to address this knowledge gap because obligately sexual individuals often coexist with multiple independently derived obligately asexual lineages. This unusual situation allows direct comparisons both between sexual and asexual P. antipodarum and across populations that differ in the relative frequency of sexual individuals. As such, P. antipodarum has received a great deal of attention as a model system for the maintenance of sex in nature and is also used as a model for environmental toxicology and biological invasions. Molecular genetic resources for P. antipodarum will thus be useful to investigators in a variety of biological fields. We used 454 sequencing of cDNA libraries to generate transcriptomes from two sexual and two asexual P. antipodarum lineages. A de novo assembly of 116.7 Mb of sequence reads produced 41 396 contigs, and sequence similarity-based Gene Ontology annotations were obtained for 3740 contigs. We detected 408 315 SNP loci and 7315 microsatellite loci, which together represent the first genome-scale resource available for P. antipodarum. Raw 454 read sequences, contig sequences, annotation data and polymorphism data are publicly available in a searchable online database and for download at http://www.biology.uiowa.edu/neiman/transcriptome.php.