Uncovering the Cellular and Molecular Mechanisms of Synapse Formation and Functional Specificity Using Central Neurons of Lymnaea stagnalis

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.

Tissue-specific evaluation of suitable reference genes for RT-qPCR in the pond snail, Lymnaea stagnalis

Reverse transcription quantitative PCR (RT-qPCR) is a robust technique for the quantification and comparison of gene expression. To obtain reliable results with this method, one or more reference genes must be employed to normalize expression measurements among treatments or tissue samples. Candidate reference genes must be validated to ensure that they are stable prior to use in qPCR experiments. The pond snail (Lymnaea stagnalis) is a common research organism, particularly in the areas of learning and memory, and is an emerging model for the study of biological asymmetry, biomineralization, and evolution and development. However, no systematic assessment of qPCR reference genes has been performed in this animal. Therefore, the aim of our research was to identify stable reference genes to normalize gene expression data from several commonly studied tissues in L. stagnalis as well as across the entire body. We evaluated a panel of seven reference genes across six different tissues in L. stagnalis with RT-qPCR. The genes included: elongation factor 1-alpha, glyceraldehyde-3-phosphate dehydrogenase, beta-actin, beta-tubulin, ubiquitin, prenylated rab acceptor protein 1, and a voltage gated potassium channel. These genes exhibited a wide range of expression levels among tissues. The tissue-specific stability of each of the genes was consistent when measured by the standard stability assessment algorithms: geNorm, NormFinder, BestKeeper, and RefFinder. Our data indicate that the most stable reference genes vary among the tissues that we examined (central nervous system, tentacles, lips, penis, foot, mantle). Our results were generally congruent with those obtained from similar studies in other molluscs. Given that a minimum of two reference genes are recommended for data normalization, we provide suggestions for strong pairs of reference genes for single- and multi-tissue analyses of RT-qPCR data in L. stagnalis.

The importance of identified neurons in gastropod molluscs to neuroscience

Gastropod molluscs have large neurons that are uniquely identifiable across individuals and across species based on neuroanatomical and neurochemical criteria, facilitating research into neural signaling and neural circuits. Novel neuropeptides have been identified through RNA sequencing and mass spectroscopic analysis of single neurons. The roles of peptides and other signaling molecules including second messengers have been placed in the context of small circuits that control simple behaviors. Despite the stereotypy, neurons vary over time in their activity in large ensembles. Furthermore, there is both intra-species and inter-species variation in synaptic properties and gene expression. Research on gastropod identified neurons highlights the features that might be expected to be stable in more complex systems when trying to identify cell types.