Transcriptome profiling of Lymnaea stagnalis (Gastropoda) for ecoimmunological research

Molecular Insights into the Marine Gastropod Olivancillaria urceus: Transcriptomic and Proteopeptidomic Approaches Reveal Polypeptides with Putative Therapeutic Potential

The marine environment is a rich source of new biotechnologies and products. Bottom trawling for shrimp species such as Xiphopenaeus kroyeri and Farfantepenaeus brasiliensis leads to the unintentional capture of non-target species, known as bycatch, which includes a variety of marine life that are often discarded without economic value. A common bycatch species on the southeast coast of Brazil is Olivancillaria urceus (O. urceus), a carnivorous gastropod that feeds mainly on bivalves. Despite its abundance, this species is still little studied, especially for biotechnological applications. Other gastropods such as Conus are known for their diverse and potent toxins, which offer great potential for pharmacological discoveries. In this study, an omics approach, including transcriptomics and proteopeptidomics, was applied to explore O. urceus at the molecular level. The transcriptome of the muscle foot/mantle led to the annotation of 19,097 genes via Gene Ontology, identifying 20 toxin-like transcripts identified considering the Gastropod class. The proteome fraction confirmed 2179 transcripts, including sequences with toxin activity, such as conotoxin precursors, Conodipine-P3, and BPTI/Kunitz domain-containing proteins. In addition, 9663 peptides of 1484 precursor proteins were detected in the peptide fraction, including 2 sequences representing neurotoxins. The identification of these sequences could lead to the discovery of new molecules with therapeutic potential.

The genome of the simultaneously hermaphroditic snail Lymnaea stagnalis reveals an evolutionary expansion of FMRFamide-like receptors

The great pond snail Lymnaea stagnalis has served as a model organism for over a century in diverse disciplines such as neurophysiology, evolution, ecotoxicology and developmental biology. To support both established uses and newly emerging research interests we have performed whole genome sequencing (avg.176 × depth), assembly and annotation of a single individual derived from an inbred line. These efforts resulted in a final assembly of 943 Mb (L50 = 257; N50 = 957,215) with a total of 22,499 predicted gene models. The mitogenome was found to be 13,834 bp long and similarly organized as in other lymnaeid species, with minor differences in location of tRNA genes. As a first step towards understanding the hermaphroditic reproductive biology of L. stagnalis, we identified molecular receptors, specifically nuclear receptors (including newly discovered 2xDNA binding domain-NRs), G protein-coupled receptors, and receptor tyrosine kinases, that may be involved in the cellular specification and maintenance of simultaneously active male and female reproductive systems. A phylogenetic analysis of one particular family of GPCRs (Rhodopsin neuropeptide FMRFamide-receptor-like genes) shows a remarkable expansion that coincides with the occurrence of simultaneous hermaphroditism in the Euthyneura gastropods. As some GPCRs and NRs also showed qualitative differences in expression in female (albumen gland) and male (prostate gland) organs, it is possible that separate regulation of male and female reproductive processes may in part have been enabled by an increased abundance of receptors in the transition from a separate-sexed state to a hermaphroditic condition. These findings will support efforts to pair receptors with their activating ligands, and more generally stimulate deeper insight into the mechanisms that underlie the modes of action of compounds involved in neuroendocrine regulation of reproduction, induced toxicity, and development in L. stagnalis, and molluscs in general.

Transcriptome-based network analysis reveals hub immune genes and pathways of hepatopancreas against LPS in Amphioctopus fangsiao

The hepatopancreas is the biggest digestive organ in Amphioctopus fangsiao (A. fangsiao), but also undertakes critical functions like detoxification and immune defense. Generally, pathogenic bacteria or endotoxin from the gut microbiota would be arrested and detoxified in the hepatopancreas, which could be accompanied by the inevitable immune responses. In recent years, studies related to cephalopods immune have been increasing, but the molecular mechanisms associated with the hepatopancreatic immunity are still unclear. In this study, lipopolysaccharide (LPS), a major component of the cell wall of Gram-negative bacteria, was used for imitating bacteria infection to stimulate the hepatopancreas of A. fangsiao. To investigate the immune process happened in A. fangsiao hepatopancreas, we performed transcriptome analysis of hepatopancreas tissue after LPS injection, and identified 2615 and 1943 differentially expressed genes (DEGs) at 6 and 24 h post-injection, respectively. GO and KEGG enrichment analysis showed that these DEGs were mainly involved in immune-related biological processes and signaling pathways, including ECM-receptor interaction signaling pathway, Phagosome signaling pathway, Lysosome signaling pathway, and JAK-STAT signaling pathways. The function relationships between these DEGs were further analyzed through protein-protein interaction (PPI) networks. It was found that Mtor, Mapk14 and Atm were the three top interacting DEGs under LPS stimulation. Finally, 15 hub genes involving multiple KEGG signaling pathways and PPI relationships were selected for qRT-PCR validation. In this study, for the first time we explored the molecular mechanisms associated with hepatopancreatic immunity in A. fangsiao using a PPI networks approach, and provided new insights for understanding hepatopancreatic immunity in A. fangsiao.

A newly identified scallop MyD88 interacts with TLR and functions in innate immunity

Myeloid differentiation factor-88 (MyD88) is a key adaptor of the toll-like receptor (TLR) signaling pathway and plays a crucial role in innate immune signal transduction in animals. However, the MyD88-mediated signal transduction mechanism in shellfish has not been well studied. In this study, a new MyD88 gene, CfMyD88-2, was identified in the Zhikong scallop, Chlamys farreri. The 1779 bp long open reading frame encodes 592 amino acids. The N-terminus of CfMyD88-2 contains a conserved death domain (DD), followed by a TIR (TLR/Interleukin-1 Receptor) domain. The results of the multi-sequence comparison showed that the TIR domain sequences were highly conserved. Phylogenetic analysis revealed that CfMyD88-2 was first associated with Mizuhopecten yessoensis MyD88-4 and Argopecten irradians MyD88-4. CfMyD88-2 mRNA was expressed in all scallop tissues, as detected by qRT-PCR, and the expression level was the highest in the mantle and hepatopancreas. In addition, CfMyD88-2 mRNA expression significantly increased after pathogen-associated molecular patterns (PAMPs, such as lipopolysaccharide, peptidoglycan, or polyinosinic-polycytidylic acid) stimulation. The results of the co-immunoprecipitation experiments in HEK293T cells showed that both CfMyD88-1 and CfMyD88-2 interacted with the TLR protein of scallops, suggesting the existence of more than one functional TLR-MyD88 signaling axis in scallops. Dual luciferase reporter gene assays indicated that the overexpressed CfMyD88-2 in HEK293T cells activated interferon (IFN) α, IFN-β, IFN-γ, and NF-κB reporter genes, indicating that the protein has multiple functions. The results of the subcellular localization experiment uncovered that CfMyD88-2 was mainly localized in the cytoplasm of human cells. In summary, the novel identified CfMyD88-2 can respond to the challenge of PAMPs, participate in TLR immune signaling, and may activate downstream effector genes such as NF-κB gene. These research results will be useful in advancing the theory of innate immunity in invertebrates and provide a reference for the selection of disease-resistant scallops in the future.

Transcriptome sequencing of the endangered land snail Karaftohelix adamsi from the Island Ulleung: De novo assembly, annotation, valuation of fitness genes and SSR markers

BACKGROUND

The Bradybaenidae snail Karaftohelix adamsi is endemic to Korea, with the species tracked from Island Ulleung in North Gyeongsang Province of South Korea. K. adamsi has been classified under the Endangered Wildlife Class II species of Korea and poses a severe risk of extinction following habitat disturbances. With no available information at the DNA (genome) or mRNA (transcriptome) level for the species, conservation by utilizing informed molecular resources seems difficult.

OBJECTIVE

In this study, we used the Illumina short-read sequencing and Trinity de novo assembly to draft the reference transcriptome of K. adamsi.

RESULTS

After assembly, 13,753 unigenes were obtained of which 10,511 were annotated to public databases (a maximum of 10,165 unigenes found homologs in PANM DB). A total of 6,351, 3,535, 358, and 3,407 unigenes were ascribed to the functional categories under KOG, GO, KEGG, and IPS, respectively. The transcripts such as the HSP 70, aquaporin, TLR, and MAPK, among others, were screened as putative functional resources for adaptation. DNA transposons were found to be thickly populated in comparison to retrotransposons in the assembled unigenes. Further, 2,164 SSRs were screened with the promiscuous presence of dinucleotide repeats such as AC/GT and AG/CT.

CONCLUSION

The transcriptome-guided discovery of molecular resources in K. adamsi will not only serve as a basis for functional genomics studies but also provide sustainable tools to be utilized for the protection of the species in the wild. Moreover, the development of polymorphic SSRs is valuable for the identification of species from newer habitats and cross-species genotyping.

Yolk proteins of the schistosomiasis vector snail Biomphalaria glabrata revealed by multi-omics analysis

Vitellogenesis is the most important process in animal reproduction, in which yolk proteins play a vital role. Among multiple yolk protein precursors, vitellogenin (Vtg) is a well-known major yolk protein (MYP) in most oviparous animals. However, the nature of MYP in the freshwater gastropod snail Biomphalaria glabrata remains elusive. In the current study, we applied bioinformatics, tissue-specific transcriptomics, ovotestis-targeted proteomics, and phylogenetics to investigate the large lipid transfer protein (LLTP) superfamily and ferritin-like family in B. glabrata. Four members of LLTP superfamily (BgVtg1, BgVtg2, BgApo1, and BgApo2), one yolk ferritin (Bg yolk ferritin), and four soma ferritins (Bg ferritin 1, 2, 3, and 4) were identified in B. glabrata genome. The proteomic analysis demonstrated that, among the putative yolk proteins, BgVtg1 was the yolk protein appearing in the highest amount in the ovotestis, followed by Bg yolk ferritin. RNAseq profile showed that the leading synthesis sites of BgVtg1 and Bg yolk ferritin are in the ovotestis (presumably follicle cells) and digestive gland, respectively. Phylogenetic analysis indicated that BgVtg1 is well clustered with Vtgs of other vertebrates and invertebrates. We conclude that, vitellogenin (BgVtg1), not yolk ferritin (Bg yolk ferritin), is the major yolk protein precursor in the schistosomiasis vector snail B. glabrata.

Expression levels of serine proteases, their homologs, and prophenoloxidase in the Eriocheir sinensis with hepatopancreatic necrosis syndrome (HPNS) and their expression regulation by Runt

The occurrence of hepatopancreatic necrosis syndrome (HPNS) has seriously affected the sustainable development of Chinese mitten crab (Eriocheir sinensis) farming industry. Limited studies have focused on the immune responses in crabs with HPNS. Serine proteases (SPs) and SP homologues (SPHs) play important roles in the innate immunity of crustaceans. This study investigated the effects of HPNS on the expression levels of genes related to prophenoloxidase (proPO) activation system, and the relationship between Runt transcription factor and the transcriptions of these genes. Eight SPs and five SPHs (SPH1-4, Mas) were identified from E. sinensis. SPs contain a catalytic triad of "HDS", while SPHs lack a catalytic residue. SPs and SPHs all contain a conservative Tryp_SPc domain. Evolutionary analysis showed that EsSPs, EsSPHs, EsPO, and EsRunt were clustered with SPs, SPHs, POs, and Runts of other arthropods, respectively. In crabs with HPNS, the expression levels of six SPs (1, 3, 4, 6, 7, and 8), five SPHs, and PO were significantly upregulated in the hepatopancreas. The knockdown of EsRunt could evidently decrease the expression levels of four SPs (3, 4, 5 and 8), five SPHs (SPH1-4, Mas), and PO. Therefore, the occurrence of HPNS activates the proPO system. Furthermore, the expression levels of partial genes related to proPO system were regulated by Runt. The activation of innate immune system may be a strategy for crabs with HPNS to improve immunity and fight diseases. Our study provides a new understanding of the relationship between HPNS and innate immunity.

Molecular identification of a novel antimicrobial peptide in giant Triton snail Charonia tritonis: mRNA profiles for tissues and its potential antibacterial activity

Antimicrobial peptides (AMPs) play an important role in innate immunity against microorganisms. AMPs is an effective antibacterial agent, and the chances of causing pathogens to develop is very low. However, there is little information about AMPs in the giant Triton snail Charonia tritonis. In this research, an antimicrobial peptide gene (termed Ct-20534) was identified in C. tritonis. The open reading frame of Ct-20534 is 381 bp in size and it encodes a basic peptide precursor containing 126 amino acids. Ct-20534 gene was found to be expressed in all five tissues examined by real-time fluorescence quantitative PCR (qPCR), but the highest expression was found in the proboscis. This is the first report that antibacterial peptides have been found in C. tritonis, and it has been proved that Ct-20534 has antibacterial activity against Gram-positive bacteria and Gram-negative bacteria, among which the activity of Staphylococcus aureus is most significantly inhibited, this suggests that the newly discovered antimicrobial peptides in C. tritonis may play an important role in the immune system and bacterial resistance of C. tritonis. This study presents the discovery of a newly identified antibacterial peptide from C. tritonis, with its structural properties fully characterized and potent antibacterial activity confirmed. The results provide essential fundamental data for the development of preventive and therapeutic measures against aquatic animal diseases, which in turn can promote the sustainable and stable growth of the aquaculture industry and create economic benefits. Additionally, this research lays the foundation for future development of novel anti-infective drugs.

The transcriptomic analysis of Planorbarius corneus hemocytes (Gastropoda) naturally infected with Bilharziella polonica (Schistosomatidae)

The study of molluscan innate immunity is essential for understanding the evolution of the immune system. An advance in the knowledge of their immune system can be achieved by increasing the number of model species. Our study focuses on the immunity of Planorbarius corneus, a pulmonate snail widely distributed in Eurasia. These snails are intermediate hosts of many trematodes, including Bilharziella polonica (Schistosomatidae). In this paper we obtained and analyzed transcriptomes of hemocytes of uninfected snails Planorbarius corneus and snails naturally infected with Bilharziella polonica. The transcriptomes were found to contain transcripts encoding all major groups of immune factors previously described for other gastropods. Pathogen-recognition molecules were the most diverse group of immune factors. Comparison of the transcriptomes of the infected and the uninfected molluscs showed that the expression of some genes changed during infection. Our results extend the knowledge of immune responses of pulmonate snails to trematode invasion and promote P. corneus as a new model for the study of molluscan defence reactions.

Cellular and humoral immune response between snail hosts and their parasites

In invertebrates, the innate immune system protects against a wide range of microbiological infections. Several immunological processes are involved in the interactive immune response between snails and their parasites, including phagocytosis, nitric oxide synthesis, phenol oxidase activity, lysozymes, and lectin formation. The immunological responses connected to the interaction between snails and parasites are discussed in detail in the current research. Understanding the nature of these interactive reactions will enable scientists to explore approaches to eliminate and cure parasitic infections.

BJ, Silvio B, Fabio T. Identification and characterization of the kynurenine pathway in the pond snail Lymnaea stagnalis

Dysregulation of the kynurenine pathway (KP) is implicated in many human diseases and disorders, from immunological, metabolic, neurodegenerative, and neuropsychiatric conditions to cancer, and represents an appealing target for new therapeutic approaches. In this intricate scenario, invertebrates, like Lymnaea stagnalis (LS), provide a flexible tool to unravel the complexity of the KP. Starting from the available LS genome and transcriptome, we identified putative transcripts of all KP enzymes containing an ORF; each predicted protein possessed a high degree of sequence conservation to known orthologues of other invertebrate and vertebrate model organisms. Sequences were confirmed by qualitative PCR and sequencing. At the same time, the qRT-PCR analysis revealed that Lym IDO-like, Lym TDO-like, Lym AFMID-like, Lym KMO-like, Lym AADAT-like, Lym KYAT I/III-like, Lym KYNU-like, Lym HAAO-like, and Lym ACMSD-like showed widespread tissue expression. Then, tryptophan, kynurenine, kynurenic acid, anthranilic acid, 3-hydroxy-kynurenine, xanthurenic acid, picolinic acid, and quinolinic acid were identified in the hemolymph of LS by UHPLC-Q exactive mass spectrometer. Our study provides the most thorough characterization to date of the KP in an invertebrate model, supporting the value of LS for future functional studies of this pathway at the cellular, synaptic, and behavioral levels.

In search of the Aplysia immunome: an in silico study

The immune repertoires of mollusks beyond commercially important organisms such as the pacific oyster Crassostrea gigas or vectors for human pathogens like the bloodfluke planorb Biomphalaria glabrata are understudied. Despite being an important model for neural aging and the role of inflammation in neuropathic pain, the immune repertoire of Aplysia californica is poorly understood. Recent discovery of a neurotropic nidovirus in Aplysia has highlighted the need for a better understanding of the Aplysia immunome. To address this gap in the literature, the Aplysia reference genome was mined using InterProScan and OrthoFinder for putative immune genes. The Aplysia genome encodes orthologs of all critical components of the classical Toll-like receptor (TLR) signaling pathway. The presence of many more TLRs and TLR associated adapters than known from vertebrates suggest yet uncharacterized, novel TLR associated signaling pathways. Aplysia also retains many nucleotide receptors and antiviral effectors known to play a key role in viral defense in vertebrates. However, the absence of key antiviral signaling adapters MAVS and STING in the Aplysia genome suggests divergence from vertebrates and bivalves in these pathways. The resulting immune gene set of this in silico study provides a basis for interpretation of future immune studies in this important model organism.

Toll like receptors and their evolution in the lymnaeid freshwater snail species Radix auricularia and Lymnaea stagnalis, key intermediate hosts for zoonotic trematodes

One of the major evolutionarily conserved pathways in innate immunity of invertebrates is the toll-like receptor (TLR) pathway. However, little is known of the TLR protein family in gastropod molluscs despite their role in the transmission of human diseases, especially the common lymnaeid freshwater snail species Radix auricularia and Lymnaea stagnalis, key intermediate hosts of zoonotic trematodes. Using comparative genomics and gene prediction approaches utilising the freshwater snail Biomphalaria glabrata genome as a reference ten putative TLR proteins were identified in both R. auricularia and L. stagnalis. Phylogenetic analyses revealed that unlike other molluscs the lymnaeid species also possessed class 1 TLRs, previously thought to be unique to B. glabrata. Gene duplication events were also seen across the TLR classes in the lymnaeids with several of the genes appearing to exist as potential tandem elements in R. auricularia. Each predicted TLR was shown to possess the typical the leucine-rich repeat extracellular and TIR intracellular domains and both single cysteine clusters and multiple cysteine clusters TLRs were identified in both lymnaeid species. Principle component analyses of 3D models of the predicted TLRs showed that class 1 and 5 proteins did not cluster based on similarity of structure, suggested to be potential adaptation to a range of pathogens. This study provides the first detailed account of TLRs in lymnaeids and affords a platform for further research into the role of these proteins into susceptibility and compatibility of these snails with trematodes and their role in transmission.

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.

Examining adaptive evolution of immune activity: opportunities provided by gastropods in the age of 'omics'

Parasites threaten all free-living organisms, including molluscs. Understanding the evolution of immune defence traits in natural host populations is crucial for predicting their long-term performance under continuous infection risk. Adaptive trait evolution requires that traits are subject to selection (i.e. contribute to organismal fitness) and that they are heritable. Despite broad interest in the evolutionary ecology of immune activity in animals, the understanding of selection on and evolutionary potential of immune defence traits is far from comprehensive. For instance, empirical observations are only rarely in line with theoretical predictions of immune activity being subject to stabilizing selection. This discrepancy may be because ecoimmunological studies can typically cover only a fraction of the complexity of an animal immune system. Similarly, molecular immunology/immunogenetics studies provide a mechanistic understanding of immunity, but neglect variation that arises from natural genetic differences among individuals and from environmental conditions. Here, we review the current literature on natural selection on and evolutionary potential of immune traits in animals, signal how merging ecological immunology and genomics will strengthen evolutionary ecological research on immunity, and indicate research opportunities for molluscan gastropods for which well-established ecological understanding and/or 'immune-omics' resources are already available. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

New Insights Into Biomphalysin Gene Family Diversification in the Vector Snail Biomphalaria glabrata

Aerolysins initially characterized as virulence factors in bacteria are increasingly found in massive genome and transcriptome sequencing data from metazoans. Horizontal gene transfer has been demonstrated as the main way of aerolysin-related toxins acquisition in metazoans. However, only few studies have focused on their potential biological functions in such organisms. Herein, we present an extensive characterization of a multigene family encoding aerolysins - named biomphalysin - in Biomphalaria glabrata snail, the intermediate host of the trematode Schistosoma mansoni. Our results highlight that duplication and domestication of an acquired bacterial toxin gene in the snail genome result in the acquisition of a novel and diversified toxin family. Twenty-three biomphalysin genes were identified. All are expressed and exhibited a tissue-specific expression pattern. An in silico structural analysis was performed to highlight the central role played by two distinct domains i) a large lobe involved in the lytic function of these snail toxins which constrained their evolution and ii) a small lobe which is structurally variable between biomphalysin toxins and that matched to various functional domains involved in moiety recognition of targets cells. A functional approach suggests that the repertoire of biomphalysins that bind to pathogens, depends on the type of pathogen encountered. These results underline a neo-and sub-functionalization of the biomphalysin toxins, which have the potential to increase the range of effectors in the snail’s immune arsenal.