Whole genome analysis of a schistosomiasis-transmitting freshwater snail

Different metazoan parasites, different transcriptomic responses, with new insights on parasitic castration by digenetic trematodes in the schistosome vector snail Biomphalaria glabrata

BACKGROUND

Gastropods of the genus Biomphalaria (Family Planorbidae) are exploited as vectors by Schistosoma mansoni, the most common causative agent of human intestinal schistosomiasis. Using improved genomic resources, overviews of how Biomphalaria responds to S. mansoni and other metazoan parasites can provide unique insights into the reproductive, immune, and other systems of invertebrate hosts, and their responses to parasite challenges.

RESULTS

Using Illumina-based RNA-Seq, we compared the responses of iM line B. glabrata at 2, 8, and 40 days post-infection (dpi) to single infections with S. mansoni, Echinostoma paraensei (both digenetic trematodes) or Daubaylia potomaca (a nematode parasite of planorbid snails). Responses were compared to unexposed time-matched control snails. We observed: (1) each parasite provoked a distinctive response with a predominance of down-regulated snail genes at all time points following exposure to either trematode, and of up-regulated genes at 8 and especially 40dpi following nematode exposure; (2) At 2 and 8dpi with either trematode, several snail genes associated with gametogenesis (particularly spermatogenesis) were down-regulated. Regarding the phenomenon of trematode-mediated parasitic castration in molluscs, we define for the first time a complement of host genes that are targeted, as early as 2dpi when trematode larvae are still small; (3) Differential gene expression of snails with trematode infection at 40dpi, when snails were shedding cercariae, was unexpectedly modest and revealed down-regulation of genes involved in the production of egg mass proteins and peptide processing; and (4) surprisingly, D. potomaca provoked up-regulation at 40dpi of many of the reproduction-related snail genes noted to be down-regulated at 2 and 8dpi following trematode infection. Happening at a time when B. glabrata began to succumb to D. potomaca, we hypothesize this response represents an unexpected form of fecundity compensation. We also document expression patterns for other Biomphalaria gene families, including fibrinogen domain-containing proteins (FReDs), C-type lectins, G-protein coupled receptors, biomphalysins, and protease and protease inhibitors.

CONCLUSIONS

Our study is relevant in identifying several genes involved in reproduction that are targeted by parasites in the vector snail B. glabrata and that might be amenable to manipulation to minimize their ability to serve as vectors of schistosomes.

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.

Accelerating snail vector genomics

The three most important genera of snails for the transmission of schistosomes are Bulinus, Biomphalaria and Oncomelania. Each of these genera, found in two distantly related families, includes species that act as the intermediate host for one of the three most widespread schistosome species infecting humans, Schistosoma haematobium, S. mansoni and S. japonicum, respectively. An important step in the fight against schistosomiasis in Asia has been taken with the publication of the article "Chromosome-level genome assembly of Oncomelania hupensis: the intermediate snail host of Schistosoma japonicum", which means that genomes for all three major genera, including species across three continents, are now available in the public domain. This includes the first genomes of African snail vectors, namely Biomphalaria sudanica, Bi. pfeifferi and Bulinus truncatus, as well as high-quality chromosome level assemblies for South American Bi. glabrata. Most importantly, the wealth of new genomic and transcriptomic data is helping to establish the specific molecular mechanisms that underly compatibility between snails and their schistosomes, which although diverse and complex, may help to identify potential targets dictating host parasite interactions that can be utilised in future transmission control strategies. This new work on Oncomelania hupensis and indeed studies on other snail vectors, which provide deep insights into the genome, will stimulate research that may well lead to new and much needed control interventions.

DNA Methylation Machinery in Gastropod Mollusks

The role of DNA methylation in mollusks is just beginning to be understood. This review synthesizes current knowledge on this potent molecular hallmark of epigenetic control in gastropods-the largest class of mollusks and ubiquitous inhabitants of diverse habitats. Their DNA methylation machinery shows a high degree of conservation in CG maintenance methylation mechanisms, driven mainly by DNMT1 homologues, and the presence of MBD2 and MBD2/3 proteins as DNA methylation readers. The mosaic-like DNA methylation landscape occurs mainly in a CG context and is primarily confined to gene bodies and housekeeping genes. DNA methylation emerges as a critical regulator of reproduction, development, and adaptation, with tissue-specific patterns being observed in gonadal structures. Its dynamics also serve as an important regulatory mechanism underlying learning and memory processes. DNA methylation can be affected by various environmental stimuli, including as pathogens and abiotic stresses, potentially impacting phenotypic variation and population diversity. Overall, the features of DNA methylation in gastropods are complex, being an essential part of their epigenome. However, comprehensive studies integrating developmental stages, tissues, and environmental conditions, functional annotation of methylated regions, and integrated genomic-epigenomic analyses are lacking. Addressing these knowledge gaps will advance our understanding of gastropod biology, ecology, and evolution.

The genome and transcriptome of the snail Biomphalaria sudanica s.l.: immune gene diversification and highly polymorphic genomic regions in an important African vector of Schistosoma mansoni

BACKGROUND

Control and elimination of schistosomiasis is an arduous task, with current strategies proving inadequate to break transmission. Exploration of genetic approaches to interrupt Schistosoma mansoni transmission, the causative agent for human intestinal schistosomiasis in sub-Saharan Africa and South America, has led to genomic research of the snail vector hosts of the genus Biomphalaria. Few complete genomic resources exist, with African Biomphalaria species being particularly underrepresented despite this being where the majority of S. mansoni infections occur. Here we generate and annotate the first genome assembly of Biomphalaria sudanica sensu lato, a species responsible for S. mansoni transmission in lake and marsh habitats of the African Rift Valley. Supported by whole-genome diversity data among five inbred lines, we describe orthologs of immune-relevant gene regions in the South American vector B. glabrata and present a bioinformatic pipeline to identify candidate novel pathogen recognition receptors (PRRs).

RESULTS

De novo genome and transcriptome assembly of inbred B. sudanica originating from the shoreline of Lake Victoria (Kisumu, Kenya) resulted in a haploid genome size of ~ 944.2 Mb (6,728 fragments, N50 = 1.067 Mb), comprising 23,598 genes (BUSCO = 93.6% complete). The B. sudanica genome contains orthologues to all described immune genes/regions tied to protection against S. mansoni in B. glabrata, including the polymorphic transmembrane clusters (PTC1 and PTC2), RADres, and other loci. The B. sudanica PTC2 candidate immune genomic region contained many PRR-like genes across a much wider genomic region than has been shown in B. glabrata, as well as a large inversion between species. High levels of intra-species nucleotide diversity were seen in PTC2, as well as in regions linked to PTC1 and RADres orthologues. Immune related and putative PRR gene families were significantly over-represented in the sub-set of B. sudanica genes determined as hyperdiverse, including high extracellular diversity in transmembrane genes, which could be under pathogen-mediated balancing selection. However, no overall expansion in immunity related genes was seen in African compared to South American lineages.

CONCLUSIONS

The B. sudanica genome and analyses presented here will facilitate future research in vector immune defense mechanisms against pathogens. This genomic/transcriptomic resource provides necessary data for the future development of molecular snail vector control/surveillance tools, facilitating schistosome transmission interruption mechanisms in Africa.

A haplotype-like, chromosome-level assembled and annotated genome of Biomphalaria glabrata, an important intermediate host of schistosomiasis and the best studied model of schistosomiasis vector snails

Schistosomiasis is one of the world's most devastating parasitic diseases, afflicting 251 million people globally. The Neotropical snail Biomphalaria glabrata is an important intermediate host of the human blood fluke Schistosoma mansoni and a predominant model for schistosomiasis research. To fully exploit this model snail for biomedical research, here we report a haplotype-like, chromosome-level assembled and annotated genome of the homozygous iM line of B. glabrata that we developed at the University of New Mexico. Using multiple sequencing platforms, including Illumina, PacBio, and Omni-C sequencing, 18 sequence contact matrices representing 18 haploid chromosomes (2n = 36) were generated (337x genome coverage), and 96.5% of the scaffold sequences were anchored to the 18 chromosomes. Protein-coding genes (n = 34,559), non-coding RNAs (n = 2,406), and repetitive elements (42.52% of the genome) were predicted for the whole genome, and detailed annotations for individual chromosomes were also provided. Using this genomic resource, we have investigated the genomic structure and organization of the Toll-like receptor (TLR) and fibrinogen-domain containing protein (FReD) genes, the two important immune-related gene families. Notably, TLR-like genes are scattered on 13 chromosomes. In contrast, almost all (39 of 40) fibrinogen-related genes (FREPs) (immunoglobulin superfamily (IgSF) + fibrinogen (FBG)) are clustered within a 5-million nucleotide region on chromosome 13, yielding insight into mechanisms involved in the diversification of FREPs. This is the first genome of schistosomiasis vector snails that has been assembled at the chromosome level, annotated, and analyzed. It serves as a valuable resource for a deeper understanding of the biology of vector snails, especially Biomphalaria snails.

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.

A comprehensive review on the dynamic role of toll-like receptors (TLRs) in frontier aquaculture research and as a promising avenue for fish disease management

Toll-like receptors (TLRs) represent a conserved group of germline-encoded pattern recognition receptors (PRRs) that recognize pathogen-associated molecular patterns (PAMPs) and play a crucial role in inducing the broadly acting innate immune response against pathogens. In recent years, the detection of 21 different TLR types in various fish species has sparked interest in exploring the potential of TLRs as targets for boosting immunity and disease resistance in fish. This comprehensive review offers the latest insights into the diverse facets of fish TLRs, highlighting their history, classification, architectural insights through 3D modelling, ligands recognition, signalling pathways, crosstalk, and expression patterns at various developmental stages. It provides an exhaustive account of the distinct TLRs induced during the invasion of specific pathogens in various fish species and delves into the disparities between fish TLRs and their mammalian counterparts, highlighting the specific contribution of TLRs to the immune response in fish. Although various facets of TLRs in some fish, shellfish, and molluscs have been described, the role of TLRs in several other aquatic organisms still remained as potential gaps. Overall, this article outlines frontier aquaculture research in advancing the knowledge of fish immune systems for the proper management of piscine maladies.

One-health approach on the future application of snails: a focus on snail-transmitted parasitic diseases

Snails are fascinating molluscs with unique morphological and physiological adaptive features to cope with various environments. They have traditionally been utilized as food and feed sources in many regions of the world. The future exploitation of alternative nutrient sources, like snails, is likely to increase further. Snails, however, also serve as an intermediate host for several zoonotic parasites. A category of parasitic infections, known as snail-transmitted parasitic diseases (STPDs), is harmful to humans and animals and is mainly driven by various trematodes, cestodes, and nematodes. The environment plays a crucial role in transmitting these parasites, as suitable habitats and conditions can facilitate their growth and proliferation in snails. In light of diverse environmental settings and biologically categorized snail species, this review evaluates the dynamics of significant STPDs of zoological importance. Additionally, possible diagnostic approaches for the prevention of STPDs are highlighted. One-health measures must be considered when employing snails as an alternative food or feed source to ensure the safety of snail-based products and prevent any adverse effects on humans, animals, and the environment.

Hemoglobin wonders: a fascinating gas transporter dive into molluscs

Hemoglobin (Hb) has been identified in at least 14 molluscan taxa so far. Research spanning over 130 years on molluscan Hbs focuses on their genes, protein structures, functions, and evolution. Molluscan Hbs are categorized into single-, two-, and multiple-domain chains, including red blood cell, gill, and extracellular Hbs, based on the number of globin domains and their respective locations. These Hbs exhibit variation in assembly, ranging from monomeric and dimeric to higher-order multimeric forms. Typically, molluscan Hbs display moderately high oxygen affinity, weak cooperativity, and varying pH sensitivity. Hb's potential role in antimicrobial pathways could augment the immune defense of bivalves, which may be a complement to their lack of adaptive immunity. The role of Hb as a respiratory protein in bivalves likely originated from the substitution of hemocyanin. Molluscan Hbs demonstrate adaptive evolution in response to environmental changes via various strategies (e.g. increasing Hb types, multimerization, and amino acid residue substitutions at key sites), enhancing or altering functional properties for habitat adaptation. Concurrently, an increase in Hb assembly diversity, coupled with a downward trend in oxygen affinity, is observed during molluscan differentiation and evolution. Hb in Protobranchia, Heteroconchia, and Pteriomorphia bivalves originated from separate ancestors, with Protobranchia inheriting a relative ancient molluscan Hb gene. In bivalves, extracellular Hbs share a common origin, while gill Hbs likely emerged from convergent evolution. In summary, research on molluscan Hbs offers valuable insights into the origins, biological variations, and adaptive evolution of animal Hbs.

Core gut microbes Cloacibacterium and Aeromonas associated with different gastropod species could be persistently transmitted across multiple generations

BACKGROUND

Studies on the gut microbiota of animals have largely focused on vertebrates. The transmission modes of commensal intestinal bacteria in mammals have been well studied. However, in gastropods, the relationship between gut microbiota and hosts is still poorly understood. To gain a better understanding of the composition of gut microbes and their transmission routes in gastropods, a large-scale and long-term experiment on the dynamics and transmission modes of gut microbiota was conducted on freshwater snails.

RESULTS

We analyzed 244 microbial samples from the digestive tracts of freshwater gastropods and identified Proteobacteria and Bacteroidetes as dominant gut microbes. Aeromonas, Cloacibacterium, and Cetobacterium were identified as core microbes in the guts, accounting for over 50% of the total sequences. Furthermore, both core bacteria Aeromonas and Cloacibacterium, were shared among 7 gastropod species and played an important role in determining the gut microbial community types of both wild and cultured gastropods. Analysis of the gut microbiota at the population level, including wild gastropods and their offspring, indicated that a proportion of gut microbes could be consistently vertically transmitted inheritance, while the majority of the gut microbes resulted from horizontal transmission. Comparing cultured snails to their wild counterparts, we observed an increasing trend in the proportion of shared microbes and a decreasing trend in the number of unique microbes among wild gastropods and their offspring reared in a cultured environment. Core gut microbes, Aeromonas and Cloacibacterium, remained persistent and dispersed from wild snails to their offspring across multiple generations. Interestingly, under cultured environments, the gut microbiota in wild gastropods could only be maintained for up to 2 generations before converging with that of cultured snails. The difference observed in gut bacterial metabolism functions was associated with this transition. Our study also demonstrated that the gut microbial compositions in gastropods are influenced by developmental stages and revealed the presence of Aeromonas and Cloacibacterium throughout the life cycle in gastropods. Based on the dynamics of core gut microbes, it may be possible to predict the health status of gastropods during their adaptation to new environments. Additionally, gut microbial metabolic functions were found to be associated with the adaptive evolution of gastropods from wild to cultured environments.

CONCLUSIONS

Our findings provide novel insights into the dynamic processes of gut microbiota colonization in gastropod mollusks and unveil the modes of microbial transmission within their guts. Video Abstract.

Chromosome-level genome assembly of the deep-sea snail Phymorhynchus buccinoides provides insights into the adaptation to the cold seep habitat

BACKGROUND

The deep-sea snail Phymorhynchus buccinoides belongs to the genus Phymorhynchus (Neogastropoda: Raphitomidae), and it is a dominant specie in the cold seep habitat. As the environment of the cold seep is characterized by darkness, hypoxia and high concentrations of toxic substances such as hydrogen sulfide (H2S), exploration of the diverse fauna living around cold seeps will help to uncover the adaptive mechanisms to this unique habitat. In the present study, a chromosome-level genome of P. buccinoides was constructed and a series of genomic and transcriptomic analyses were conducted to explore its molecular adaptation mechanisms to the cold seep environments.

RESULTS

The assembled genome size of the P. buccinoides was approximately 2.1 Gb, which is larger than most of the reported snail genomes, possibly due to the high proportion of repetitive elements. About 92.0% of the assembled base pairs of contigs were anchored to 34 pseudo-chromosomes with a scaffold N50 size of 60.0 Mb. Compared with relative specie in the shallow water, the glutamate regulative and related genes were expanded in P. buccinoides, which contributes to the acclimation to hypoxia and coldness. Besides, the relatively high mRNA expression levels of the olfactory/chemosensory genes in osphradium indicate that P. buccinoides might have evolved a highly developed and sensitive olfactory organ for its orientation and predation. Moreover, the genome and transcriptome analyses demonstrate that P. buccinoides has evolved a sulfite-tolerance mechanism by performing H2S detoxification. Many genes involved in H2S detoxification were highly expressed in ctenidium and hepatopancreas, suggesting that these tissues might be critical for H2S detoxification and sulfite tolerance.

CONCLUSIONS

In summary, our report of this chromosome-level deep-sea snail genome provides a comprehensive genomic basis for the understanding of the adaptation strategy of P. buccinoides to the extreme environment at the deep-sea cold seeps.

The genome and transcriptome of the snail Biomphalaria sudanica s.l.: Immune gene diversification and highly polymorphic genomic regions in an important African vector of Schistosoma mansoni

Background

Control and elimination of schistosomiasis is an arduous task, with current strategies proving inadequate to break transmission. Exploration of genetic approaches to interrupt Schistosoma mansoni transmission, the causative agent for human intestinal schistosomiasis in sub-Saharan Africa and South America, has led to genomic research of the snail vector hosts of the genus Biomphalaria. Few complete genomic resources exist, with African Biomphalaria species being particularly underrepresented despite this being where the majority of S. mansoni infections occur. Here we generate and annotate the first genome assembly of Biomphalaria sudanica sensu lato, a species responsible for S. mansoni transmission in lake and marsh habitats of the African Rift Valley. Supported by whole-genome diversity data among five inbred lines, we describe orthologs of immune-relevant gene regions in the South American vector B. glabrata and present a bioinformatic pipeline to identify candidate novel pathogen recognition receptors (PRRs).

Results

De novo genome and transcriptome assembly of inbred B. sudanica originating from the shoreline of Lake Victoria (Kisumu, Kenya) resulted in a haploid genome size of ~944.2 Mb (6732 fragments, N50=1.067 Mb), comprising 23,598 genes (BUSCO=93.6% complete). The B. sudanica genome contains orthologues to all described immune genes/regions tied to protection against S. mansoni in B. glabrata. The B. sudanica PTC2 candidate immune genomic region contained many PRR-like genes across a much wider genomic region than has been shown in B. glabrata, as well as a large inversion between species. High levels of intra-species nucleotide diversity were seen in PTC2, as well as in regions linked to PTC1 and RADres orthologues. Immune related and putative PRR gene families were significantly over-represented in the sub-set of B. sudanica genes determined as hyperdiverse, including high extracellular diversity in transmembrane genes, which could be under pathogen-mediated balancing selection. However, no overall expansion in immunity related genes were seen in African compared to South American lineages.

Conclusions

The B. sudanica genome and analyses presented here will facilitate future research in vector immune defense mechanisms against pathogens. This genomic/transcriptomic resource provides necessary data for the future development of molecular snail vector control/surveillance tools, facilitating schistosome transmission interruption mechanisms in Africa.

Metagenomic Analysis Reveals Variations in Gut Microbiomes of the Schistosoma mansoni-Transmitting Snails Biomphalaria straminea and Biomphalaria glabrata

Biomphalaria snails play a crucial role in the transmission of the human blood fluke Schistosoma mansoni. The gut microbiota of intermediate hosts is known to influence their physiological functions, but little is known about its composition and role in Biomphalaria snails. To gain insights into the biological characteristics of these freshwater intermediate hosts, we conducted metagenomic sequencing on Biomphalaria straminea and B. glabrata to investigate variations in their gut microbiota. This study revealed that the dominant members of the gut microbiota in B. glabrata belong to the phyla Bacteroidetes and Proteobacteria, which were also found to be the top two most abundant gut bacteria in B. straminea. We identified Firmicutes, Acidovorax and Bosea as distinctive gut microbes in B. straminea, while Aeromonas, Cloacibacterium and Chryseobacterium were found to be dependent features of the B. glabrata gut microbiota. We observed significant differences in the community structures and bacterial functions of the gut microbiota between the two host species. Notably, we found a distinctive richness of antibiotic resistance genes (ARGs) associated with various classes of antibiotics, including bacitracin, chloramphenicol, tetracycline, sulfonamide, penicillin, cephalosporin_ii and cephalosporin_i, fluoroquinolone, aminoglycoside, beta-lactam, multidrug and trimethoprim, in the digestive tracts of the snails. Furthermore, this study revealed the potential correlations between snail gut microbiota and the infection rate of S. mansoni using Spearman correlation analysis. Through metagenomic analysis, our study provided new insights into the gut microbiota of Biomphalaria snails and how it is influenced by host species, thereby enhancing our understanding of variant patterns of gut microbial communities in intermediate hosts. Our findings may contribute to future studies on gastropod-microbe interactions and may provide valuable knowledge for developing snail control strategies to combat schistosomiasis in the future.

Molluscicidal and Cercaricidal Effects of Myrciaria floribunda Essential Oil Nanoemulsion

Schistosomiasis is a tropical disease transmitted in an aqueous environment by cercariae from the Schistosoma genus. This disease affects 200 million people living in risk areas around the world. The control of schistosomiasis is realized by chemotherapy, wastewater sanitation, health education, and mollusk control using molluscicidal agents. This work evaluates the effects of a nanoemulsion containing essential oil from Myrciaria floribunda leaves as a molluscicidal and cercaricidal agent against Biomphalaria glabrata mollusks and Schistosoma mansoni cercariae. The Myrciaria floribunda essential oil from leaves showed nerolidol, β-selinene, 1,8 cineol, and zonarene as major constituents. The formulation study suggested the F3 formulation as the most promising nanoemulsion with polysorbate 20 and sorbitan monooleate 80 (4:1) with 5% (w/w) essential oil as it showed a smaller droplet size of approximately 100 nm with a PDI lower than 0.3 and prominent bluish reflection. Furthermore, this nanoemulsion showed stability after 200 days under refrigeration. The Myrciaria floribunda nanoemulsion showed LC50 values of 48.11 µg/mL, 29.66 µg/mL, and 47.02 µg/mL in Biomphalaria glabrata embryos, juveniles, and adult mollusks, respectively, after 48 h and 83.88 µg/mL for Schistosoma mansoni cercariae after 2 h. In addition, a survival of 80% was observed in Danio rerio, and the in silico toxicity assay showed lower overall human toxicity potential to the major compounds in the essential oil compared to the reference molluscicide niclosamide. These results suggest that the nanoemulsion of Myrciaria floribunda leaves may be a promising alternative for schistosomiasis control.

Genome assembly of the Korean intertidal mud-creeper Batillaria attramentaria

Batillaridae is a common gastropod family that occurs abundantly in the shallow coastal zone of the intertidal mudflats of the northwest Pacific Ocean, Australasia, and North America. In this family, Batillaria attramentaria is known for its biological invasion and colonization in estuarine and intertidal zones. It can endure and adapt the harsh intertidal conditions such as frequent temperature alteration, salinity, and air exposure. Therefore, we sequenced and assembled this Korean batillariid genome to get insight into its intertidal adaptive features. Approximately 53 Gb of DNA sequences were generated, and 863 scaffolds were assembled into a draft genome of 0.715 Gb with 97.1% BUSCO completeness value. A total of 40,596 genes were predicted. We estimated that B. attramentaria and Conus consors diverged about 230 million years ago (MYA) based on the phylogenetic analysis of closely related gastropod species. This genome study sets the footstep for genomics studies among native and introduced Batillaria populations and the Batillaridae family members.

A High-Quality Chromosome-Level Genome Assembly of a Snail Cipangopaludina cathayensis (Gastropoda: Viviparidae)

Cipangopaludina cathayensis (Gastropoda: Prosobranchia; Mesogastropoda; Viviparidae) is widely distributed in the freshwater habitats of China. It is an economically important snail with high edible and medicinal value. However, the genomic resources and the reference genome of this snail are lacking. In this study, we assembled the first chromosome-level genome of C. cathayensis. The preliminary assembly genome was 1.48 Gb in size, with a contig N50 size of 93.49 Mb. The assembled sequences were anchored to nine pseudochromosomes using Hi-C data. The final genome after Hi-C correction was 1.48 Gb, with a contig N50 of 98.49 Mb and scaffold N50 of 195.21 Mb. The anchored rate of the chromosome was 99.99%. A total of 22,702 protein-coding genes were predicted. Phylogenetic analyses indicated that C. cathayensis diverged with Bellamya purificata approximately 158.10 million years ago. There were 268 expanded and 505 contracted gene families in C. cathayensis when compared with its most recent common ancestor. Five putative genes under positive selection in C. cathayensis were identified (false discovery rate <0.05). These genome data provide a valuable resource for evolutionary studies of the family Viviparidae, and for the genetic improvement of C. cathayensis.

Understanding Snail Mucus Biosynthesis and Shell Biomineralisation through Genomic Data Mining of the Reconstructed Carbohydrate and Glycan Metabolic Pathways of the Giant African Snail (Achatina fulica)

The giant African snail (Order Stylommatophora: Family Achatinidae), Achatina fulica (Bowdich, 1822), is the most significant and invasive land snail pest. The ecological adaptability of this snail involves high growth rate, reproductive capacity, and shell and mucus production, driven by several biochemical processes and metabolism. The available genomic information for A. fulica provides excellent opportunities to hinder the underlying processes of adaptation, mainly carbohydrate and glycan metabolic pathways toward the shell and mucus formation. The authors analysed the 1.78 Gb draft genomic contigs of A. fulica to identify enzyme-coding genes and reconstruct biochemical pathways related to the carbohydrate and glycan metabolism using a designed bioinformatic workflow. Three hundred and seventy-seven enzymes involved in the carbohydrate and glycan metabolic pathways were identified based on the KEGG pathway reference in combination with protein sequence comparison, structural analysis, and manual curation. Fourteen complete pathways of carbohydrate metabolism and seven complete pathways of glycan metabolism supported the nutrient acquisition and production of the mucus proteoglycans. Increased copy numbers of amylases, cellulases, and chitinases highlighted the snail advantage in food consumption and fast growth rate. The ascorbate biosynthesis pathway identified from the carbohydrate metabolic pathways of A. fulica was involved in the shell biomineralisation process in association with the collagen protein network, carbonic anhydrases, tyrosinases, and several ion transporters. Thus, our bioinformatic workflow was able to reconstruct carbohydrate metabolism, mucus biosynthesis, and shell biomineralisation pathways from the A. fulica genome and transcriptome data. These findings could reveal several evolutionary advantages of the A. fulica snail, and will benefit the discovery of valuable enzymes for industrial and medical applications.

A genome sequence for Biomphalaria pfeifferi, the major vector snail for the human-infecting parasite Schistosoma mansoni

BACKGROUND

Biomphalaria pfeifferi is the world's most widely distributed and commonly implicated vector snail species for the causative agent of human intestinal schistosomiasis, Schistosoma mansoni. In efforts to control S. mansoni transmission, chemotherapy alone has proven insufficient. New approaches to snail control offer a way forward, and possible genetic manipulations of snail vectors will require new tools. Towards this end, we here offer a diverse set of genomic resources for the important African schistosome vector, B. pfeifferi.

METHODOLOGY/PRINCIPAL FINDINGS

Based largely on PacBio High-Fidelity long reads, we report a genome assembly size of 772 Mb for B. pfeifferi (Kenya), smaller in size than known genomes of other planorbid schistosome vectors. In a total of 505 scaffolds (N50 = 3.2Mb), 430 were assigned to 18 large linkage groups inferred to represent the 18 known chromosomes, based on whole genome comparisons with Biomphalaria glabrata. The annotated B. pfeifferi genome reveals a divergence time of 3.01 million years with B. glabrata, a South American species believed to be similar to the progenitors of B. pfeifferi which undertook a trans-Atlantic colonization < five million years ago.

CONCLUSIONS/SIGNIFICANCE

The genome for this preferentially self-crossing species is less heterozygous than related species known to be preferential out-crossers; its smaller genome relative to congeners may similarly reflect its preference for selfing. Expansions of gene families with immune relevance are noted, including the FReD gene family which is far more similar in its composition to B. glabrata than to Bulinus truncatus, a vector for Schistosoma haematobium. Provision of this annotated genome will help better understand the dependencies of trematodes on snails, enable broader comparative insights regarding factors contributing to susceptibility/ resistance of snails to schistosome infections, and provide an invaluable resource with respect to identifying and manipulating snail genes as potential targets for more specific snail control programs.

Using freshwater snail Biomphalaria glabrata (Say, 1818) as a biological model for ecotoxicology studies: a systematic review

Over time, a growing increase in human pollutants in the aquatic environment has been observed. The global presence of residues in water bodies reinforces the need to develop improved methods to detect them and evaluate their ecotoxicological effects in aquatic environments. Thus, this study aimed to present the main assays using Biomphalaria glabrata as a biological model for ecotoxicological studies. We performed a systematic literature review with data published up to June 2022 on the Web of Science, SCOPUS, Science Direct, PubMed, and SciELO databases. Thirty studies were selected for this review after screening. Biomphalaria glabrata has been studied as an ecotoxicological model for different substances through toxicity, embryotoxicity, cytotoxicity, genotoxicity, and bioaccumulation assays. Studies evaluating the impact of B. glabrata exposure to several substances have reported effects on their offspring, as well as toxicity and behavioral and reproductive effects. This review presents various assays using B. glabrata as a biological model for ecotoxicological studies. The use of a representative species of ecosystems from tropical regions is a necessary tool for tropical environmental monitoring. It was observed that the freshwater snail B. glabrata was effective for the evaluation of the ecotoxicity of several types of chemical substances, but further studies are needed to standardize the model.

Adsorption and detoxification of glyphosate and aminomethylphosphonic acid by montmorillonite clays

The co-occurrence of mixtures of glyphosate (GLP) and aminomethylphosphonic acid (AMPA) in contaminated water, soil, sediment, and plants is a cause for concern due to potential threats to the ecosystem and human health. Major routes of exposure include contact with contaminated water and soil and through consumption of crops containing GLP and AMPA residues. Calcium montmorillonite (CM) and acid-processed montmorillonite (APM) clays were investigated for their ability to tightly sorb and detoxify GLP and AMPA mixtures. In vitro adsorption and desorption isotherms and thermodynamic analysis indicated saturable Langmuir binding of both chemicals with high capacities, affinities, enthalpies, and free energies of sorption and low desorption rates. In silico computational modeling indicated that both GLP and AMPA can be readily absorbed onto clay surfaces through electrostatic interactions and hydrogen bonding. The safety and efficacy of the clays were confirmed using well-established living organisms, including an aquatic cnidarian (Hydra vulgaris), a soil nematode (Caenorhabditis elegans), and a floating plant (Lemna minor). Low levels of clay inclusion (0.05% and 0.2%) in the culture medium resulted in increased growth and protection against chemical mixtures based on multiple endpoints. Results indicated that montmorillonite clays may be used to bind mixtures of GLP and AMPA in water, soil, and plants.

Transposon-derived transcription factors across metazoans

Transposable elements (TE) could serve as sources of new transcription factors (TFs) in plants and some other model species, but such evidence is lacking for most animal lineages. Here, we discovered multiple independent co-options of TEs to generate 788 TFs across Metazoa, including all early-branching animal lineages. Six of ten superfamilies of DNA transposon-derived conserved TF families (ZBED, CENPB, FHY3, HTH-Psq, THAP, and FLYWCH) were identified across nine phyla encompassing the entire metazoan phylogeny. The most extensive convergent domestication of potentially TE-derived TFs occurred in the hydroid polyps, polychaete worms, cephalopods, oysters, and sea slugs. Phylogenetic reconstructions showed species-specific clustering and lineage-specific expansion; none of the identified TE-derived TFs revealed homologs in their closest neighbors. Together, our study established a framework for categorizing TE-derived TFs and informing the origins of novel genes across phyla.

High-Density Genetic Linkage Map of the Southern Blue-ringed Octopus (Octopodidae: Hapalochlaena maculosa)

Genetic linkage maps provide a useful resource for non-model genomes and can aid in genome reassembly to form more contiguous pseudo-chromosomes. We present the first linkage map of any cephalopod, H. maculosa, composed of 47 linkage groups (LG). A total of 2166 single nucleotide polymorphisms and 2455 presence–absence variant loci were utilised by Lep-Map3 in linkage map construction. The map length spans 2016.62 cM with an average marker distance of 0.85 cM. Integration of the recent H. maculosa genome allowed 1151 scaffolds comprising 34% of the total genomic sequence to be orientated and/or placed using 1278 markers across all 47 LG. The linkage map generated provides a new perspective on HOX gene distribution in octopods. In the H. maculosa linkage map three (SCR, LOX4 and POST1) of six identified HOX genes (HOX1/LAB, SCR, LOX2, LOX4, LOX5, POST1) were located within the same LG (LG 9). The generation of a linkage map for H. maculosa has provided a valuable resource for understanding the evolution of cephalopod genomes and will provide a base for future work.

Identification of Schistosoma mansoni miracidia attractant candidates in infected Biomphalaria glabrata using behaviour-guided comparative proteomics

Schistosomiasis, caused by infection with Schistosoma digenetic trematodes, is one of the deadliest neglected tropical diseases in the world. The Schistosoma lifecycle involves the miracidial infection of an intermediate freshwater snail host, such as Biomphalaria glabrata. Dispersing snail host-derived Schistosoma miracidia attractants has been considered a method of minimising intermediate host infections and, by extension, human schistosomiasis. The attractiveness of B. glabrata to miracidia is known to be reduced following infection; however, the relationship between duration of infection and attractiveness is unclear. Excretory-secretory proteins (ESPs) most abundant in attractive snail conditioned water (SCW) are key candidates to function as miracidia attractants. This study analysed SCW from B. glabrata that were naïve (uninfected) and at different time-points post-miracidia exposure (PME; 16h, 1-week, 2-weeks and 3-weeks PME) to identify candidate ESPs mediating Schistosoma mansoni miracidia behaviour change, including aggregation and chemoklinokinesis behaviour (random motion, including slowdown and increased turning rate and magnitude). Miracidia behaviour change was only observed post-addition of naïve and 3W-PME SCW, with other treatments inducing significantly weaker behaviour changes. Therefore, ESPs were considered attractant candidates if they were shared between naïve and 3W-PME SCW (or exclusive to the former), contained a predicted N-terminal signal peptide and displayed low identity (<50%) to known proteins outside of the Biomphalaria genus. Using these criteria, a total of 6 ESP attractant candidates were identified, including acetylcholine binding protein-like proteins and uncharacterised proteins. Tissue-specific RNA-seq analysis of the genes encoding these 6 ESPs indicated relatively high gene expression within various B. glabrata tissues, including the foot, mantle and kidney. Acetylcholine binding protein-like proteins were highly promising due to their high abundance in naïve and 3W-PME SCW, high specificity to B. glabrata and high expression in the ovotestis, from which attractants have been previously identified. In summary, this study used proteomics, guided by behavioural assays, to identify miracidia attractant candidates that should be further investigated as potential biocontrols to disrupt miracidia infection and minimise schistosomiasis.

Two novel mollusk short-form ApeC-containing proteins act as pattern recognition proteins for peptidoglycan

The Apextrin C-terminal (ApeC) domain is a new protein domain largely specific to aquatic invertebrates. In amphioxus, a short-form ApeC-containing protein (ACP) family is capable of binding peptidoglycan (PGN) and agglutinating bacteria via its ApeC domain. However, the functions of ApeC in other phyla remain unknown. Here we examined 130 ACPs from gastropods and bivalves, the first and second biggest mollusk classes. They were classified into nine groups based on their phylogenetics and architectures, including three groups of short-form ACPs, one group of apextrins and two groups of ACPs of complex architectures. No groups have orthologs in other phyla and only four groups have members in both gastropods and bivalves, suggesting that mollusk ACPs are highly diversified. We selected one bivalve ACP (CgACP1; from the oyster Crossostrea gigas) and one gastropod ACP (BgACP1; from the snail Biomphalaria glabrata) for functional experiments. Both are highly-expressed, secreted short-form ACPs and hence comparable to the amphioxus ACPs previously reported. We found that recombinant CgACP1 and BgACP1 bound with yeasts and several bacteria with different affinities. They also agglutinated these microbes, but showed no inhibiting or killing effects. Further analyses show that both ACPs had high affinities to the Lys-type PGN from S. aureus but weak or no affinities to the DAP-type PGN from Bacillus subtilis. Both recombinant ACPs displayed weak or no affinities to other microbial cell wall components, including lipopolysaccharide (LPS), lipoteichoic acid (LTA), zymosan A, chitin, chitosan and cellulose, as well as to several PGN moieties, including muramyl dipeptide (MDP), N-acetylglucosamine (GlcNAc) and N-acetylmuramic acid (MurNAc). Besides, CgACP1 had the highest expression in the gill and could be greatly up-regulated quickly after bacterial challenge. This is reminiscent of the amphioxus ACP1/2 which serve as essential mucus lectins in the gill. Taken together, the current findings from mollusk and amphioxus ACPs suggest several basic common traits for the ApeC domains, including the high affinity to Lys-type PGN, the bacterial binding and agglutinating capacity, and the role as mucus proteins to protect the mucosal surface.

Evaluation of the developmental toxicity of solvents, metals, drugs, and industrial chemicals using a freshwater snail (Biomphalaria glabrata) assay

A freshwater snail assay was employed to assess the embryotoxicity of solvents including acetone, methanol, ethanol, isopropanol, dimethyl-sulfoxide, glycerin, metals/metalloids including mercuric chloride (HgCl₂), cadmium chloride (CdCl₂,), antimony salts Sb⁺³ and Sb⁺⁵, drugs including colchicine, hydroxyurea, cyclophosphamide, an industrial chemical sodium azide (SA), an anionic surfactant dodecyl sodium sulfate-(DSS), H₂O₂ and sodium chloride (NaCl). The assay consists of exposing Biomphalaria glabrata egg masses (EM) to the substances for 96-hr and following up embryo/snail development for lethality, abnormal morphology (teratogenicity), and day of hatching up to day 10 or 14 after spawning. Based upon concentration-response relationships, LC_50%s (embryolethality), EC_50%s (teratogenicity) and IC_50%s (hatching retardation) and 95%CIs were determined for tested chemicals. The LOECs indicated that HgCl₂ (37 nM) and CdCl₂ (140 nM) are potent embryotoxic agents in snails. Teratogenic indices (TI = LC₅₀/EC₅₀) for almost all tested chemicals were lower than or close to unity suggesting that these compounds were not teratogenic in this assay. The snail assay may be adequately performed in a cost-effective standardized protocol which enables testing a number of environmental chemicals over a broad concentration range. The snail assay needs to undergo further validation to be recognized for an internationally harmonized hazard identification in ecotoxicity risk assessment.

Large-scale and small-scale population genetic structure of the medically important gastropod species Bulinus truncatus (Gastropoda, Heterobranchia)

Background

Gastropod snails remain strongly understudied, despite their important role in transmitting parasitic diseases. Knowledge of their distribution and population dynamics increases our understanding of the processes driving disease transmission. We report the first study to use high-throughput sequencing (HTS) to elucidate the population genetic structure of the hermaphroditic snail Bulinus truncatus (Gastropoda, Heterobranchia) on a regional (17–150 km) and inter-regional (1000–5400 km) scale. This snail species acts as an intermediate host of Schistosoma haematobium and Schistosoma bovis, which cause human and animal schistosomiasis respectively.

Methods

Bulinus truncatus snails were collected in Senegal, Cameroon, Egypt and France and identified through DNA barcoding. A single-end genotyping-by-sequencing (GBS) library, comprising 87 snail specimens from the respective countries, was built and sequenced on an Illumina HiSeq 2000 platform. Reads were mapped against S. bovis and S. haematobium reference genomes to identify schistosome infections, and single nucleotide polymorphisms (SNPs) were scored using the Stacks pipeline. These SNPs were used to estimate genetic diversity, assess population structure and construct phylogenetic trees of B. truncatus.

Results

A total of 10,750 SNPs were scored and used in downstream analyses. The phylogenetic analysis identified five clades, each consisting of snails from a single country but with two distinct clades within Senegal. Genetic diversity was low in all populations, reflecting high selfing rates, but varied between locations due to habitat variability. Significant genetic differentiation and isolation by distance patterns were observed at both spatial scales, indicating that gene flow is not strong enough to counteract the effects of population bottlenecks, high selfing rates and genetic drift. Remarkably, the population genetic differentiation on a regional scale (i.e. within Senegal) was as large as that between populations on an inter-regional scale. The blind GBS technique was able to pick up parasite DNA in snail tissue, demonstrating the potential of HTS techniques to further elucidate the role of snail species in parasite transmission.

Conclusions

HTS techniques offer a valuable toolbox to further investigate the population genetic patterns of intermediate schistosome host snails and the role of snail species in parasite transmission.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05445-x.

The divergent genome of Scorpion Group 1 (SG1) intracellular bacterium from the venom glands of Vaejovis smithi (Scorpiones: Vaejovidae)

Scorpions were among the first animals on land around 430 million years ago. Like many arachnids, scorpions have evolved complex venoms used to paralyze their prey and for self-defense. Here we sequenced and analyzed the metagenomic DNA from venom glands from Vaejovis smithi scorpions. A metagenome-assembled genome (MAG) of 624,025 bp was obtained corresponding to the previously reported Scorpion Group 1 (SG1). The SG1 genome from venom glands had a low GC content (25.8%) characteristic of reduced genomes, many hypothetical genes and genes from the reported minimal set of bacterial genes. Phylogenomic reconstructions placed the uncultured SG1 distant from other reported bacteria constituting a taxonomic novelty. By PCR we detected SG1 in all tested venom glands from 30 independent individuals. Microscopically, we observed SG1 inside epithelial cells from the venom glands using FISH and its presence in scorpion embryos suggested that SG1 is transferred from mother to offspring.

Compatibility between snails and schistosomes: insights from new genetic resources, comparative genomics, and genetic mapping

The freshwater snail Biomphalaria glabrata is an important intermediate host of the parasite Schistosoma mansoni that causes human intestinal schistosomiasis. To better understand vector snail biology and help advance innovative snail control strategies, we have developed a new snail model consisting of two homozygous B. glabrata lines (iM line and iBS90) with sharply contrasting schistosome-resistance phenotypes. We produced and compared high-quality genome sequences for iM line and iBS90 which were assembled from 255 (N50 = 22.7 Mb) and 346 (N50 = 19.4 Mb) scaffolds, respectively. Using F2 offspring bred from the two lines and the newly generated iM line genome, we constructed 18 linkage groups (representing the 18 haploid chromosomes) covering 96% of the genome and identified three new QTLs (quantitative trait loci), two involved in snail resistance/susceptibility and one relating to body pigmentation. This study provides excellent genomic resources for unveiling complex vector snail biology, reveals genomic difference between resistant and susceptible lines, and offers novel insights into genetic mechanism of the compatibility between snail and schistosome.

Common aquatic pollutants modify hemocyte immune responses in Biomphalaria glabrata

Disruptions to reproductive health in wildlife species inhabiting polluted environments is often found to occur alongside compromised immunity. However, research on impacts of aquatic pollution on freshwater mollusc immune responses is limited despite their importance as vectors of disease (Schistosomiasis) in humans, cattle and wild mammals. We developed an in vitro ‘tool-kit’ of well-characterized quantitative immune tests using Biomphalaria glabrata hemocytes. We exposed hemocytes to environmentally-relevant concentrations of common aquatic pollutants (17β-estradiol, Bisphenol-A and p,p’-DDE) and measured key innate immune responses including motility, phagocytosis and encapsulation. Additionally, we tested an extract of a typical domestic tertiary treated effluent as representative of a ‘real-world’ mixture of chemicals. Encapsulation responses were stimulated by p,p’-DDE at low doses but were suppressed at higher doses. Concentrations of BPA (above 200 ng/L) and p,p’-DDE (above 500 ng/L) significantly inhibited phagocytosis compared to controls, whilst hemocyte motility was reduced by all test chemicals and the effluent extract in a dose-dependent manner. All responses occurred at chemical concentrations considered to be below the cytotoxic thresholds of hemocytes. This is the first time a suite of in vitro tests has been developed specifically in B. glabrata with the purpose of investigating the impacts of chemical pollutants and an effluent extract on immunity. Our findings indicate that common aquatic pollutants alter innate immune responses in B. glabrata, suggesting that pollutants may be a critical, yet overlooked, factor impacting disease by modulating the dynamics of parasite transmission between molluscs and humans.

Biomphalaria embryotoxicity test (BET): 60 years of research crossing boundaries for developing standard protocols

Snail's embryotoxicity test is a suitable approach for toxicity assay of traditional and emerging pollutants, environmental risk assessment, as well as screening and development of new molluscicides. Among the snail species, Biomphalaria spp. has been indicated as a promising model system for developing standardized test protocols for assessing the chemical toxicity using early developmental stages. Thus, the current study aimed to review the data available in the scientific literature concerning the experimental approach, type of chemicals and the response of multiple biomarkers (survival, hatching rate, development delays, morphological and behavior changes) in snail embryos applied in toxicity tests. Revised data showed that the use of Biomphalaria embryos to assess chemical toxicity began in 1962. Snail's embryotoxicity test was applied mainly for analyzing the toxicity and development of new molluscicides, while its use in ecotoxicological studies is emerging. Biomphalaria glabrata was the main species analyzed. Embryos exposed to chemicals showed bioaccumulation, mortality, hatching inhibition, development delays, and morphological malformations, which were classified into four categories (hydropic, shell, cephalic and unspecified malformations). Besides, research gaps and recommendations for future research are indicated. Overall, the results showed that the Biomphalaria embryotoxicity test (BET) is a suitable tool for toxicity and health risk assessment.

Molecular characterization of thioester-containing proteins in Biomphalaria glabrata and their differential gene expression upon Schistosoma mansoni exposure

Schistosomiasis is a disease caused by trematode parasites of the genus Schistosoma that affects approximately 200 million people worldwide. Schistosomiasis has been a persistent problem in endemic areas as there is no vaccine available, currently used anti-helmintic medications do not prevent reinfection, and most concerning, drug resistance has been documented in laboratory and field isolates. Thus, alternative approaches to curtail this human disease are warranted. Understanding the immunobiology of the obligate intermediate host of these parasites, which include the freshwater snail Biomphalaria glabrata, may facilitate the development of novel methods to stop or reduce transmission to humans. Molecules from the thioester-containing protein (TEP) superfamily have been shown to be involved in immunological functions in many animals including corals and humans. In this study we identified, characterized, and compared TEP transcripts and their expression upon S. mansoni exposure in resistant and susceptible strains of B. glabrata snails. Results showed the expression of 11 unique TEPs in B. glabrata snails. These transcripts present high sequence identity at the nucleotide and putative amino acid levels between susceptible and resistant strains. Further analysis revealed differences in several TEPs’ constitutive expression levels between resistant and susceptible snail strains, with C3-1, C3-3, and CD109 having higher constitutive expression levels in the resistant (BS90) strain, whereas C3-2 and TEP-1 showed higher constitutive expression levels in the susceptible (NMRI) strain. Furthermore, TEP-specific response to S. mansoni miracidia exposure reiterated their differential expression, with resistant snails upregulating the expression of both TEP-4 and TEP-3 at 2 h and 48 h post-exposure, respectively. Further understanding the diverse TEP genes and their functions in invertebrate animal vectors will not only expand our knowledge in regard to this ancient family of immune proteins, but also offer the opportunity to identify novel molecular targets that could aid in the efforts to develop control methods to reduce schistosomiasis transmission.

A Dissenters’ View on AppleSnail Immunobiology

We stand as dissenters against the acceptance of scientific knowledge that has not been built on empirical data. With this in mind, this review synthesizes selected aspects of the immunobiology of gastropods and of apple snails (Ampullariidae) in particular, from morphological to molecular and “omics” studies. Our trip went through more than two centuries of history and was guided by an evo-devo hypothesis: that the gastropod immune system originally developed in the mesenchymal connective tissue of the reno-pericardial complex, and that in that tissue some cells differentiated into hematopoietically committed progenitor cells that integrate constitutive hemocyte aggregations in the reno-pericardial territory, whether concentrated in the pericardium or the kidney in a species-specific manner. However, some of them may be freed from those aggregations, circulate in the blood, and form distant contingent aggregations anywhere in the body, but always in response to intruders (i.e., pathogens or any other immune challenge). After that, we reviewed the incipient immunology of the Ampullariidae by critically revising the findings in Pomacea canaliculata and Marisa cornuarietis, the only ampullariid species that have been studied in this respect, and we attempted to identify the effectors and the processes in which they are involved. Particularly for P. canaliculata, which is by far the most studied species, we ask which hemocytes are involved, in which tissues or organs are integrated, and what cellular reactions to intruders this species has in common with other animals. Furthermore, we wondered what humoral factors could also integrate its internal defense system. Among the cellular defenses, we give an outstanding position to the generation of hemocyte nodules, which seems to be an important process for these snails, serving the isolation and elimination of intruders. Finally, we discuss hematopoiesis in apple snails. There have been contrasting views about some of these aspects, but we envision a hematopoietic system centered in the constitutive hemocyte islets in the ampullariid kidney.

Analysis of rhodopsin G protein-coupled receptor orthologs reveals semiochemical peptides for parasite (Schistosoma mansoni) and host (Biomphalaria glabrata) interplay

Schistosomiasis is a medically significant disease caused by helminth parasites of the genus Schistosoma. The schistosome life cycle requires chemically mediated interactions with an intermediate (aquatic snail) and definitive (human) host. Blocking parasite development within the snail stage requires improved understanding of the interactions between the snail host and the Schistosoma water-borne free-living form (miracidium). Innovations in snail genomics and aquatic chemical communication provide an ideal opportunity to explore snail-parasite coevolution at the molecular level. Rhodopsin G protein-coupled receptors (GPCRs) are of particular interest in studying how trematode parasites navigate towards their snail hosts. The potential role of GPCRs in parasites makes them candidate targets for new antihelminthics that disrupt the intermediate host life-cycle stages, thus preventing subsequent human infections. A genomic-bioinformatic approach was used to identify GPCR orthologs between the snail Biomphalaria glabrata and miracidia of its obligate parasite Schistosoma mansoni. We show that 8 S. mansoni rhodopsin GPCRs expressed within the miracidial stage share overall amino acid similarity with 8 different B. glabrata rhodopsin GPCRs, particularly within transmembrane domains, suggesting conserved structural features. These GPCRs include an orphan peptide receptor as well as several with strong sequence homologies with rhabdomeric opsin receptors, a serotonin receptor, a sulfakinin (SK) receptor, an allatostatin-A (buccalin) receptor and an FMRFamide receptor. Buccalin and FMRFa peptides were identified in water conditioned by B. glabrata, and we show synthetic buccalin and FMRFa can stimulate significant rates of change of direction and turn-back responses in S. mansoni miracidia. Ortholog GPCRs were identified in S. mansoni miracidia and B. glabrata. These GPCRs may detect similar ligands, including snail-derived odorants that could facilitate miracidial host finding. These results lay the foundation for future research elucidating the mechanisms by which GPCRs mediate host finding which can lead to the potential development of novel anti-schistosome interventions.

Dynamic evolutionary history of spiralian-specific TALE homeobox genes in mollusks

Homeobox genes play essential roles in the early development of many animals. Although the repertoire of most homeobox genes, including three amino acid loop extension (TALE)-type homeobox genes, is conserved in animals, spiralian-TALE (SPILE) genes are a notable exception. In this study, SPILE genes were extracted from the genomic data of 22 mollusk species and classified into four clades (-A/C, -B, -D, and -E) to determine which SPILE genes exhibit dynamic repertoire changes. While SPILE-D and -E duplications were rarely observed, SPILE-B duplication was observed in the bivalve lineage and SPILE-A/C duplication was observed in multiple clades. Conversely, most or all SPILE genes were lost in cephalopods and in some gastropod lineages. SPILE gene expression patterns were also analyzed in multiple mollusk species using publicly available RNA-seq data. The majority of SPILE genes examined, particularly those in the A/C- and B-clades, were specifically expressed during early development, suggesting that most SPILE genes exert specific roles in early development. This comprehensive cataloging and characterization revealed a dynamic evolutionary history, including SPILE-A/C and -B gene duplications and the loss of SPILE genes in several lineages. Furthermore, this study provides a useful resource for studying the molecular mechanism of spiralian early development and the evolution of young and lineage-specific transcription factors.

Metagenome-Assembled Genomes Reveal Mechanisms of Carbohydrate and Nitrogen Metabolism of Schistosomiasis-Transmitting Vector Biomphalaria Glabrata

Biomphalaria glabrata transmits schistosomiasis mansoni which poses considerable risks to hundreds of thousands of people worldwide, and is widely used as a model organism for studies on the snail-schistosome relationship. Gut microbiota plays important roles in multiple aspects of host including development, metabolism, immunity, and even behavior; however, detailed information on the complete diversity and functional profiles of B. glabrata gut microbiota is still limited. This study is the first to reveal the gut microbiome of B. glabrata based on metagenome-assembled genome (MAG). A total of 28 gut samples spanning diet and age were sequenced and 84 individual microbial genomes with ≥ 70% completeness and ≤ 5% contamination were constructed. Bacteroidota and Proteobacteria were the dominant bacteria in the freshwater snail, unlike terrestrial organisms harboring many species of Firmicutes and Bacteroidota. The microbial consortia in B. glabrata helped in the digestion of complex polysaccharide such as starch, hemicellulose, and chitin for energy supply, and protected the snail from food poisoning and nitrate toxicity. Both microbial community and metabolism of B. glabrata were significantly altered by diet. The polysaccharide-degrading bacterium Chryseobacterium was enriched in the gut of snails fed with high-digestibility protein and high polysaccharide diet (HPHP). Notably, B. glabrata as a mobile repository can escalate biosafety issues regarding transmission of various pathogens such as Acinetobacter nosocomialis and Vibrio parahaemolyticus as well as multiple antibiotic resistance genes in the environment and to other organisms. IMPORTANCE The spread of aquatic gastropod Biomphalaria glabrata, an intermediate host of Schistosoma mansoni, exacerbates the burden of schistosomiasis disease worldwide. This study provides insights into the importance of microbiome for basic biological activities of freshwater snails, and offers a valuable microbial genome resource to fill the gap in the analysis of the snail-microbiota-parasite relationship. The results of this study clarified the reasons for the high adaptability of B. glabrata to diverse environments, and further illustrated the role of B. glabrata in accumulation of antibiotic resistance in the environment and spread of various pathogens. These findings have important implications for further exploration of the control of snail dissemination and schistosomiasis from a microbial perspective.

How many single-copy orthologous genes from whole genomes reveal deep gastropod relationships?

The Gastropoda contains 80% of existing mollusks and is the most diverse animal class second only to the Insecta. However, the deep phylogeny of gastropods has been controversial for a long time. Especially the position of Patellogastropoda is a major uncertainty. Morphology and some mitochondria studies concluded that Patellogastropoda is likely to be sister to all other gastropods (Orthogastropoda hypothesis), while transcriptomic and other mitogenomic studies indicated that Patellogastropoda and Vetigastropoda are sister taxa (Psilogastropoda). With the release of high-quality genomes, orthologous genes can be better identified and serve as powerful candidates for phylogenetic analysis. The question is, given the current limitations on the taxon sampling side, how many markers are needed to provide robust results. Here, we identified single-copy orthologous genes (SOGs) from 14 gastropods species with whole genomes available which cover five main gastropod subclasses. We generated different datasets from 395 to 1610 SOGs by allowing species missing in different levels. We constructed gene trees of each SOG, and inferred species trees from different collections of gene trees. We found as the number of SOGs increased, the inferred topology changed from Patellogastropoda being sister to all other gastropods to Patellogastropoda being sister to Vetigastropoda + Neomphalina (Psilogastropoda s.l.), with considerable support. Our study thus rejects the Orthogastropoda concept showing that the selection of the representative species and use of sufficient informative sites greatly influence the analysis of deep gastropod phylogeny.

Mitochondrial genome of Bulinus truncatus (Gastropoda: Lymnaeoidea): Implications for snail systematics and schistosome epidemiology

Many freshwater snails of the genus Bulinus act as intermediate hosts in the life-cycles of schistosomes in Africa and adjacent regions. Currently, 37 species of Bulinus representing four groups are recognised. The mitochondrial cytochrome c oxidase subunit 1 (cox1) gene has shown utility for identifying and differentiating Bulinus species and groups, but taxonomic relationships based on genetic data are not entirely consistent with those inferred using morphological and biological features. To underpin future systematic studies of members of the genus, we characterised here the mitochondrial genome of Bulinus truncatus (from a defined laboratory strain) using a combined second- and third-generation sequencing and informatics approach, enabling taxonomic comparisons with other planorbid snails for which mitochondrial (mt) genomes were available. Analyses showed consistency in gene order and length among mitochondrial genomes of representative planorbid snails, with the lowest and highest nucleotide diversities being in the cytochrome c oxidase and nicotinamide dehydrogenase subunit genes, respectively. This first mt genome for a representative of the genus Bulinus should provide a useful resource for future investigations of the systematics, population genetics, epidemiology and/or ecology of Bulinus and related snails. The sequencing and informatic workflow employed here should find broad applicability to a range of other snail intermediate hosts of parasitic trematodes.

An Overview of Transcriptional Responses of Schistosome-Susceptible (M line) or -Resistant (BS-90) Biomphalaria glabrata Exposed or Not to Schistosoma mansoni Infection

Background

We seek to provide a comprehensive overview of transcriptomics responses of immune-related features of the gastropod Biomphalaria glabrata (Bg) following exposure to Schistosoma mansoni (Sm), a trematode causing human schistosomiasis. Responses of schistosome-susceptible (M line, or SUS) and -resistant (BS-90, or RES) Bg strains were characterized following exposure to Sm for 0.5, 2, 8 or 40 days post-exposure (dpe).

Methods

RNA-Seq and differential expression analysis were undertaken on 56 snails from 14 groups. We considered 7 response categories: 1) constitutive resistance factors; 2) constitutive susceptibility factors; 3) generalized stress responses; 4) induced resistance factors; 5) resistance factors suppressed in SUS snails; 6) suppressed/manipulated factors in SUS snails; and 7) tolerance responses in SUS snails. We also undertook a gene co-expression network analysis. Results from prior studies identifying schistosome resistance/susceptibility factors were examined relative to our findings.

Results

A total of 792 million paired-end reads representing 91.2% of the estimated 31,985 genes in the Bg genome were detected and results for the 7 categories compiled and highlighted. For both RES and SUS snails, a single most supported network of genes with highly correlated expression was found.

Conclusions

1) Several constitutive differences in gene expression between SUS and RES snails were noted, the majority over-represented in RES; 2) There was little indication of a generalized stress response shared by SUS and RES snails at 0.5 or 2 dpe; 3) RES snails mounted a strong, multi-faceted response by 0.5 dpe that carried over to 2 dpe; 4) The most notable SUS responses were at 40 dpe, in snails shedding cercariae, when numerous features were either strongly down-regulated indicative of physiological distress or parasite manipulation, or up-regulated, suggestive of tolerance or survival-promoting effects; 5) Of 55 genes previously identified in genome wide mapping studies, 29 (52.7%) were responsive to Sm, as were many familiar resistance-associated genes (41.0%) identified by other means; 6) Both network analysis and remarkably specific patterns of expression of lectins and G protein-coupled receptors in categories 4, 6 and 7 were indicative of orchestrated responses of different suites of genes in SUS or RES snails following exposure to Sm.

Genome of the ramshorn snail Biomphalaria straminea-an obligate intermediate host of schistosomiasis

BACKGROUND

Schistosomiasis, or bilharzia, is a parasitic disease caused by trematode flatworms of the genus Schistosoma. Infection by Schistosoma mansoni in humans results when cercariae emerge into water from freshwater snails in the genus Biomphalaria and seek out and penetrate human skin. The snail Biomphalaria straminea is native to South America and is now also present in Central America and China, and represents a potential vector host for spreading schistosomiasis. To date, genomic information for the genus is restricted to the neotropical species Biomphalaria glabrata. This limits understanding of the biology and management of other schistosomiasis vectors, such as B. straminea.

FINDINGS

Using a combination of Illumina short-read, 10X Genomics linked-read, and Hi-C sequencing data, our 1.005 Gb B. straminea genome assembly is of high contiguity, with a scaffold N50 of 25.3 Mb. Transcriptomes from adults were also obtained. Developmental homeobox genes, hormonal genes, and stress-response genes were identified, and repeat content was annotated (40.68% of genomic content). Comparisons with other mollusc genomes (including Gastropoda, Bivalvia, and Cephalopoda) revealed syntenic conservation, patterns of homeobox gene linkage indicative of evolutionary changes to gene clusters, expansion of heat shock protein genes, and the presence of sesquiterpenoid and cholesterol metabolic pathway genes in Gastropoda. In addition, hormone treatment together with RT-qPCR assay reveal a sesquiterpenoid hormone responsive system in B. straminea, illustrating that this renowned insect hormonal system is also present in the lophotrochozoan lineage.

CONCLUSION

This study provides the first genome assembly for the snail B. straminea and offers an unprecedented opportunity to address a variety of phenomena related to snail vectors of schistosomiasis, as well as evolutionary and genomics questions related to molluscs more widely.

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.

Immunological Diversity Is a Cornerstone of Organismal Defense and Allorecognition across Metazoa

The ongoing arms race between hosts and microbes has fueled the evolution of novel strategies for diversifying the molecules involved in immune responses. Characterization of immune systems from an ever-broadening phylogenetic range of organisms reveals that there are many mechanisms by which this diversity can be generated and maintained. Diversification strategies operate at the level of populations, genomes, genes, and even individual transcripts. Lineage-specific innovations have been cataloged within the immune systems of both invertebrates and vertebrates. Furthermore, somatic diversification of immune receptor genes has now been described in jawless vertebrates and some invertebrate species. In addition to pathogen detection, immunological diversity plays important roles in several distinct allorecognition systems. In this Brief Review, we highlight some of the evolutionary innovations employed by a variety of metazoan species to generate the molecular diversity required to detect a vast array of molecules in the context of both immune response and self/nonself-recognition.

Not all Is SET for Methylation: Evolution of Eukaryotic Protein Methyltransferases

Dynamic posttranslational modifications to canonical histones that constitute the nucleosome (H2A, H2B, H3, and H4) control all aspects of enzymatic transactions with DNA. Histone methylation has been studied heavily for the past 20 years, and our mechanistic understanding of the control and function of individual methylation events on specific histone arginine and lysine residues has been greatly improved over the past decade, driven by excellent new tools and methods. Here, we will summarize what is known about the distribution and some of the functions of protein methyltransferases from all major eukaryotic supergroups. The main conclusion is that protein, and specifically histone, methylation is an ancient process. Many taxa in all supergroups have lost some subfamilies of both protein arginine methyltransferases (PRMT) and the heavily studied SET domain lysine methyltransferases (KMT). Over time, novel subfamilies, especially of SET domain proteins, arose. We use the interactions between H3K27 and H3K36 methylation as one example for the complex circuitry of histone modifications that make up the "histone code," and we discuss one recent example (Paramecium Ezl1) for how extant enzymes that may resemble more ancient SET domain KMTs are able to modify two lysine residues that have divergent functions in plants, fungi, and animals. Complexity of SET domain KMT function in the well-studied plant and animal lineages arose not only by gene duplication but also acquisition of novel DNA- and histone-binding domains in certain subfamilies.

Bioinformatics Approaches for Analyzing Multigene Families Encoding Immune Receptors

Genome sequences are quickly becoming available from a variety of organisms, providing researchers with an abundance of previously inaccessible information and an important source of insight into immune mechanisms. There are a variety of methods to accurately characterize genes from new genome sequences, but immune receptors pose special challenges for these techniques. Immune receptors, particularly those that directly recognize pathogens, often diverge rapidly among species and are commonly found in large, complex multigene families. Because of these characteristics, immune receptors tend to be overlooked or misannotated in large-scale genomic surveys. We describe here a strategy to characterize homologs of immune receptors and to identify putative receptors from newly assembled genome or transcriptome sequences. The description of these protocols is aimed at a typical immunologist and does not rely on substantial a priori knowledge of bioinformatics. The approach is based on using low-stringency sequence searches to identify divergent homologs. For receptors with multiple domains, the intersection of low-stringency searches can be used to identify divergent receptor sequences with high confidence. For multigene families, these predictions can be refined using sequence conservation among gene family paralogs. Assembled genome sequences serve as a critical foundation for subsequent functional characterization and remove long-standing barriers in understanding the evolution of immune recognition systems.

Identification of peptidoglycan recognition proteins in hemocytes and kidney of common periwinkle Littorinalittorea

Peptidoglycan Recognition Proteins (PGRPs) are a diverse group of proteins involved in innate immunity. In particular, PGRPs have been shown to participate in immune pattern recognition in various mollusks. However, they have not been described in Caenogastropoda, a large molluscan group comprising sea, freshwater and land snails. In this study, four short PGRPs with molecular weights ranging from 21 to 34 kDa and their isoforms were identified and structurally characterized in the kidney and hemocytic transcriptomes of a caenogastropod mollusk Littorina littorea. All of them (LlPGRP1-4) are secretory, possess a signal peptide and a characteristic N-terminal N-acetylmuramoyl-l-alanine amidase (Ami) domain with conserved Zn²⁺ binding- and amidase catalytic sites. The shortest proteins, LlPGRP1 and LlPGRP2, have no additional conserved motifs on the N-terminus. In longer and most abundantly expressed LlPGRP3 and LlPGRP4 the Ami-domain is combined with an N-terminal SH3-domain and a cysteine-rich motif, respectively. Expression analysis showed that LlPGRPs of the common periwinkle were uninvolved in the immune response to infection with trematode Himasthla elongata though they might act in antibacterial defense.

Mollusc N-glycosylation: Structures, Functions and Perspectives

Molluscs display a sophisticated N-glycan pattern on their proteins, which is, in terms of involved structural features, even more diverse than that of vertebrates. This review summarises the current knowledge of mollusc N-glycan structures, with a focus on the functional aspects of the corresponding glycoproteins. Furthermore, the potential of mollusc-derived biomolecules for medical applications is addressed, emphasising the importance of mollusc research.

Identification and diversity of fibrinogen-related protein (FREP) gene family in Haliotis discus hannai, H. rufescens, and H. laevigata and their responses to Vibrio parahemolyticus infection

Fibrinogen-related proteins (FREPs) are distributed universally in vertebrates and invertebrates. These proteins contain fibrinogen-like (FBG) domains in their C-terminal region and involve in immune responses and other aspects of physiology in invertebrates. In this study, 54 proteins that contain FBG domains or a fibrinogen_c domain were identified in Haliotis discus hannai. Comparatively, 88 and 63 FREPs were identified from the genomes of H. rufescens and H. laevigata. Most FREPs of abalones had a conserved motif containing a bound calcium ion site and a second conserved motif containing a polymerization pocket site. By sequence analysis, 394 SNPs and 11 Indels were identified in 20 FREP genes of the whole genome of H. discus hannai; 992 SNPs and 42 Indels were found in 64 FREPs of H. rufescens, and 192 SNPs and 12 Indels were found in 21 FREPs of H. laevigata. Among these SNPs, 92 missense mutation sites were identified in 26 FREP genes of H. rufescens, and 12 were identified in 8 FREP genes of H. laevigata. Due to the poor genomic integrity, annotations of the SNPs or Indels in H. discus hannai did not yield missense mutant sites. FREP genes with polymorphisms were ubiquitously expressed in all the tested tissues; however, the expression is lowest in the hemolymph. In response to Vibrio parahemolyticus infection, expression of FREP genes was significantly upregulated at different exposure times in gills, hepatopancreas, and hemolymph in H. discus hannai. Overall, this study documented the FREP genes of abalones and shed light on the role of FREPs in the innate immune system of these aquaculture species for the prevention and control of diseases.