Immune repertoire in the transcriptome of Littorina littorea reveals new trends in lophotrochozoan proto-complement evolution

Littoporins: Novel actinoporin-like proteins in caenogastropod genus Littorina

In the course of searching for genes controlling the immune system in caenogastropod mollusks, we characterized and phylogenetically placed five new actinoporin-like cytolysins expressed in periwinkles of the genus Littorina. These newly discovered proteins, named littoporins (LitP), contain a central cytolysin/lectin domain and exhibit a predicted protein fold that is almost identical to the three-dimensional structures of actinoporins. Two of these proteins, LitP-1 and LitP-2, were found to be upregulated in L. littorea kidney tissues and immune cells in response to natural and experimental infection with the trematode Himasthla elongata, suggesting their potential role as perforins in the systemic anti-trematode immune response. The primary sequence divergence of littoporins is hypothesized to be attributed to the taxonomic range of cell membranes they can recognize and permeabilize.

Parvatrema spp. (Digenea, Gymnophallidae) with parthenogenetic metacercariae: diversity, distribution and host specificity in the palaearctic

There are several species of gymnophallid digeneans in the genus Parvatrema that are unique in developing metacercariae that reproduce by parthenogenesis in the second intermediate host. Transmission of these digeneans takes place in coastal ecosystems of the North Pacific and North Atlantic seas. The first intermediate hosts are bivalves, the second ones are gastropods, and the definitive hosts are migratory birds. We integrated data accumulated over 25 years of research and differentiated a complex of five closely related species. They differ in the molluscan second intermediate hosts, distribution ranges, and life cycles patterns. The type I life cycle includes two generations of parthenogenetic metacercariae, followed by development of metacercariae which are invasive for the definitive host. In the type II life cycle, the number of generations of parthenogenetic metacercariae is unlimited, and they can also produce cercariae. These cercariae emerge into the environment and can infect new individuals of the second intermediate host. We conclude that the type I life cycle is a derived option that has evolved as a better fit to transmission in the unstable conditions in the intertidal zone. Another evolutionary trend in Parvatrema is transition from inhabiting the extrapallial space of the gastropod second intermediate host to endoparasitism in its mantle and internal organs. rDNA sequence analysis highlighted that Parvatrema spp. with parthenogenetic metacercariae form a monophyletic clade and suggested the Pacific origin of the group, with two transfers to the North Atlantic and colonisation of new second intermediate host species. Apparently the group formed in the late Pliocene-Pleistocene and diversified as a result of recurrent isolation in inshore refugia during glacial periods. We argue that parthenogenetic metacercariae in Parvatrema may serve as a model for early digenean evolution, demonstrating the first steps of adopting the molluscan first intermediate host and becoming tissue parasites.

Neuronal gene expression in two generations of the marine parasitic worm, Cryptocotyle lingua

Trematodes, or flukes, undergo intricate anatomical and behavioral transformations during their life cycle, yet the functional changes in their nervous system remain poorly understood. We investigated the molecular basis of nervous system function in Cryptocotyle lingua, a species of relevance for fisheries. Transcriptomic analysis revealed a streamlined molecular toolkit with the absence of key signaling pathways and ion channels. Notably, we observed the loss of nitric oxide synthase across the Platyhelminthes. Furthermore, we identified upregulated neuronal genes in dispersal larvae, including those involved in aminergic pathways, synaptic vesicle trafficking, TRPA channels, and surprisingly nitric oxide receptors. Using neuronal markers and in situ hybridization, we hypothesized their functional relevance to larval adaptations and host-finding strategies. Additionally, employing a behavior quantification toolkit, we assessed cercaria motility, facilitating further investigations into the behavior and physiology of parasitic flatworms. This study enhances our understanding of trematode neurobiology and provides insights for targeted antiparasitic strategies.

Unveiling novel RNA viruses in trematodes parasitizing common periwinkle: Implications for host-parasite interactions

This study characterized novel RNA viruses, parasites of parasites, or hyperparasites identified during transcriptomic analyses of two trematode species, Cryptocotyle lingua and Himasthla elongata, infecting a sea snail, Littorina littorea. According to the viral genome structures and phylogenetic analysis, Cryptolin alternavirus (ClRNAV1), Cryptolin calicivirus (ClRNAV2) and Himastelon rhabdovirus 1 (HeRNAV1) were respectively classified within the families Alternaviridae, Caliciviridae and Rhabdoviridae. They replicate species-specifically in two studied phases of trematode live cycle: intramolluscan parthenogenetic rediae and free-swimming cercariae. ClRNAV1 showed significantly higher expression in C. lingua cercariae relative to rediae. HeRNAV1's similarity to rabies viruses raises questions about its potential effects on the nervous system of H. elongata. This 'trematode rabies' could enable the use of genetically modified viruses for developing new methods to control the spread and intensity of diseases caused by trematodes.

Invertebrate C1q Domain-Containing Proteins: Molecular Structure, Functional Properties and Biomedical Potential

C1q domain-containing proteins (C1qDC proteins) unexpectedly turned out to be widespread molecules among a variety of invertebrates, despite their lack of an integral complement system. Despite the wide distribution in the genomes of various invertebrates, data on the structure and properties of the isolated and characterized C1qDC proteins, which belong to the C1q/TNF superfamily, are sporadic, although they hold great practical potential for the creation of new biotechnologies. This review not only summarizes the current data on the properties of already-isolated or bioengineered C1qDC proteins but also projects further strategies for their study and biomedical application. It has been shown that further broad study of the carbohydrate specificity of the proteins can provide great opportunities, since for many of them only interactions with pathogen-associated molecular patterns (PAMPs) was evaluated and their antimicrobial, antiviral, and fungicidal activities were studied. However, data on the properties of C1qDC proteins, which researchers originally discovered as lectins and therefore studied their fine carbohydrate specificity and antitumor activity, intriguingly show the great potential of this family of proteins for the creation of targeted drug delivery systems, vaccines, and clinical assays for the differential diagnosis of cancer. The ability of invertebrate C1qDC proteins to recognize patterns of aberrant glycosylation of human cell surfaces and interact with mammalian immunoglobulins indicates the great biomedical potential of these molecules.

The primitive complement system in molluscs

The complement system is an important immune defense mechanism that plays essential roles in both innate and adaptive immunity of vertebrates. Since complement components are identified in deuterostome and even primitive protostome species, the origin and evolution of complement system in invertebrates have been of great interest. Recently, research on the complement system in mollusc immunity has been increasing due to their importance in worldwide aquaculture, and their phylogenetic position. Complement components including C3, C1q domain containing protein (C1qDCP), C-type lectin (CTL), ficolin-like, mannose-binding lectin (MBL)-associated serine proteases like (MASPL), and factor B have been identified, suggesting the existence of complement system in molluscs. The lectin pathway has been outlined in molluscs, which is initiated by CTL with CCP domain and MASPL protein to generate C3 cleavage fragments. The molluscan C1qDCP exhibits the capability to bind human IgG, indicating the existence of possible C1qDCP-mediated activation pathway in molluscs. The activation of C3 regulates the expressions of immune effectors (cytokines and antibacterial peptides), mediates the haemocyte phagocytosis, and inhibits the bacterial growth. Some MACPF domain containing proteins may replace the missing terminal pathway in molluscs. This article provides a review of complement system in molluscs, including its components, activation mechanisms and functions in the immune response of molluscs.

The localization of Toll and Imd pathway and complement system components and their response to Vibrio infection in the nemertean Lineus ruber

BACKGROUND

Innate immunity is the first line of defense against pathogens. In animals, the Toll pathway, the Imd pathway, the complement system, and lectins are well-known mechanisms involved in innate immunity. Although these pathways and systems are well understood in vertebrates and arthropods, they are understudied in other invertebrates.

RESULTS

To shed light on immunity in the nemertean Lineus ruber, we performed a transcriptomic survey and identified the main components of the Toll pathway (e.g., myD88, dorsal/dif/NFκB-p65), the Imd pathway (e.g., imd, relish/NFκB-p105/100), the complement system (e.g., C3, cfb), and some lectins (FreD-Cs and C-lectins). In situ hybridization showed that TLRβ1, TLRβ2, and imd are expressed in the nervous system; the complement gene C3-1 is expressed in the gut; and the lectins are expressed in the nervous system, the blood, and the gut. To reveal their potential role in defense mechanisms, we performed immune challenge experiments, in which Lineus ruber specimens were exposed to the gram-negative bacteria Vibrio diazotrophicus. Our results show the upregulation of specific components of the Toll pathway (TLRα3, TLRβ1, and TLRβ2), the complement system (C3-1), and lectins (c-lectin2 and fred-c5).

CONCLUSIONS

Therefore, similarly to what occurs in other invertebrates, our study shows that components of the Toll pathway, the complement system, and lectins are involved in the immune response in the nemertean Lineus ruber. The presence of these pathways and systems in Lineus ruber, but also in other spiralians; in ecdysozoans; and in deuterostomes suggests that these pathways and systems were involved in the immune response in the stem species of Bilateria.

Immune diversity in lophotrochozoans, with a focus on recognition and effector systems

Lophotrochozoa is one of the most species-rich but immunologically poorly explored phyla. Although lack of acquired response in a narrow sense, lophotrochozoans possess various genetic mechanisms that enhance the diversity and specificity of innate immune system. Here, we review the recent advances of comparative immunology studies in lophotrochozoans with focus on immune recognition and effector systems. Haemocytes and coelomocytes are general important yet understudied player. Comparative genomics studies suggest expansion and functional divergence of lophotrochozoan immune reorganization systems is not as "homogeneous and simple" as we thought including the large-scale expansion and molecular divergence of pattern recognition receptors (PRRs) (TLRs, RLRs, lectins, etc.) and signaling adapters (MyD88s etc.), significant domain recombination of immune receptors (RLR, NLRs, lectins, etc.), extensive somatic recombination of fibrinogenrelated proteins (FREPs) in snails. Furthermore, there are repeatedly identified molecular mechanisms that generate immune effector diversity, including high polymorphism of antimicrobial peptides and proteins (AMPs), reactive oxygen and nitrogen species (RONS) and cytokines. Finally, we argue that the next generation omics tools and the recently emerged genome editing technicism will revolutionize our understanding of innate immune system in a comparative immunology perspective.

Transcriptomic profiling of immune-associated molecules in the coelomocytes of lugworm Arenicola marina (Linnaeus, 1758)

Organization and functioning of immune system remain unevenly studied in different taxa of lophotrochozoan animals. We analyzed transcriptomic data on coelomocytes of the lugworm Arenicola marina (Linnaeus, 1758; Annelida, Polychaeta) to gain insights into the molecular mechanisms involved in polychaete immunity. Coelomocytes are specialized motile cells populating coelomic fluid of annelids, responsible for cellular defense reactions and providing humoral immune factors. The transcriptome was enriched with immune-related transcripts by challenging the cells in vitro with lipopolysaccharides of Escherichia coli and Zymosan from Saccharomyces cerevisiae. Our analysis revealed a multifaceted and complex internal defense system of the lugworm. A. marina possesses orthologs of proto-complement-like factors: six thioester-containing proteins, a complement-like receptor, and a MASP-related serine protease (MReM2). A. marina coelomocytes employ pattern-recognition receptors to detect pathogens and regulate immune responses. Among them, there are 18 Toll-like receptors and various putative lectin-like proteins with evolutionary conserved and taxa-specific domains. C-type lectins and a novel family of Gal-binding and CUB domains containing receptors were the most abundant in the transcriptome. The array of pore-forming proteins in the coelomocytes was surprisingly reduced compared to that of other invertebrate species. We characterized a set of conserved proteins metabolizing reactive oxygen species and nitric oxide and expanded the arsenal of potential antimicrobial peptides. Phenoloxidase activity in immune cells of lugworm is mediated only by laccase enzyme. The described repertoire of immune-associated molecules provides valuable candidates for further functional and comparative research on the immunity of annelids.

The oxygen carrier M101 alleviates complement activation, which may be beneficial for donor organ preservation

During organ transplantation, ischemia/reperfusion injury and pre-formed anti-HLA antibodies are the main cause of delayed graft function and recovery through the activation of the complement system. By supplying oxygen during transplantation, M101 is suspected to avoid complement activation, however, a direct effect exerted by M101 on this pathway is unknown. This was tested by using functional assays (lymphocytotoxic crossmatch test, C3d Luminex-based assay, 50% complement hemolysis [CH50], and 50% alternative complement pathway [AP50/AH50]), and quantitative assays (C3, C3a, C4, C5, C5a, C6, C7, C8, C9 and sC5b-9). M101 interferes with the anti-HLA lymphocytotoxic crossmatch assay, and this effect is complement-dependent as M101 inhibits the classical complement pathway (CH50) in a dose-dependent and stable manner. Such inhibition was independent from a proteolytic effect (fractions C3 to C9) but related to a dose-dependent inhibition of the C3 convertase as demonstrated by exploring downstream the release of the anaphylatoxins (C3a and C5a), C3d, and sC5b-9. The C3 convertase inhibition in the presence of M101 was further demonstrated in the AP50/AH50 assay. In conclusion, the use of M101 avoids the activation of the complement pathway, which constitutes an additional advantage for this extracellular hemoglobin to preserve grafts from ischemia/reperfusion injury and preformed anti-HLA antibodies.

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.

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.

Sticky problems: extraction of nucleic acids from molluscs

Traditional molecular methods and omics-techniques across molluscan taxonomy increasingly inform biology of Mollusca. Recovery of DNA and RNA for such studies is challenged by common biological properties of the highly diverse molluscs. Molluscan biomineralization, adhesive structures and mucus involve polyphenolic proteins and mucopolysaccharides that hinder DNA extraction or copurify to inhibit enzyme-catalysed molecular procedures. DNA extraction methods that employ the detergent hexadecyltrimethylammoniumbromide (CTAB) to remove these contaminants importantly facilitate molecular-level study of molluscs. Molluscan pigments may stain DNA samples and interfere with spectrophotometry, necessitating gel electrophoresis or fluorometry for accurate quantification. RNA can reliably be extracted but the 'hidden break' in 28S rRNA of molluscs (like most protostomes) causes 18S and 28S rRNA fragments to co-migrate electrophoretically. This challenges the standard quality control based on the ratio of 18S and 28S rRNA, developed for deuterostome animals. High-AT content in molluscan rRNA prevents the effective purification of polyadenylated mRNA. Awareness of these matters aids the continuous expansion of molecular malacology, enabling work also with museum specimens and next-generation sequencing, with the latter imposing unprecedented demands on DNA quality. Alternative methods to extract nucleic acids from molluscs are available from literature and, importantly, from communications with others who study the molecular biology of molluscs. 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.

RNA-seq: The early response of the snail Physella acuta to the digenetic trematode Echinostoma paraensei

To analyze the response of the snail Physella acuta to Echinostoma paraensei, a compatible digenetic trematode, Illumina RNA-seq data were collected from snails with early infection (5 snails at 2 days post-exposure [DPE]) and established infection (4 snails, 8 DPE), and 7 control (unexposed) snails. A reference transcriptome (325,563 transcripts, including 98% of eukaryotic universal single-copy orthologs; BUSCO) and a draft P. acuta genome (employing available genomic Illumina reads; 799,945 scaffolds, includes 88% BUSCO genes) were assembled to guide RNA-seq analyses. Parasite exposure of P. acuta led to 10,195 differentially expressed (DE) genes at 2 DPE and 8,876 DE genes at 8 DPE with only 18% of up-regulated and 22% of down-regulated sequences shared between these time points. Gene ontology (GO) analysis yielded functional annotation of only 1.2% of DE genes but did not indicate major changes in biological activities of P. acuta between 2 and 8 DPE. Increased insights were achieved by analysis of expression profiles of 460 immune-relevant DE transcripts, identified by BLAST and InterProScan. Physella acuta has expanded gene families that encode immune-relevant domains, including CD109/TEP, GTPase IMAP, Limulus agglutination factor (dermatopontin), FReD (≥82 sequences with fibrinogen-related domains), and transcripts that combine C-type lectin (C-LECT) and C1q domains, novel among metazoa. Notably, P. acuta expressed sequences from these immune gene families at all time points, but the assemblages of unique transcripts from particular immune gene families differed between 2 and 8 DPE. The shift in profiles of DE immune genes, from early exposure to parasite establishment, suggests that compatible P. acuta initially respond to infection but switch to express immune genes that likely are less effective against E. paraensei but counter other types of (opportunistic) pathogens and parasites. We propose that the latter expression profile is part of an extended phenotype of E. paraensei, imposed upon P. acuta through parasite manipulation of the host, following successful parasite establishment in the snail after 2 DPE.

Transcriptomic analysis of Manila clam Ruditapes philippinarum under lipopolysaccharide challenge provides molecular insights into immune response

The Manila clam, Ruditapes philippinarum, is an economically important shellfish in marine aquaculture. A better understanding of the immune system in R. philippinarum will provide the basis for the development of strategies to mitigate the impact of infectious diseases affecting this species but can also be of relevance for other bivalves of commercial interest. In this study, the transcriptional response of the Manila clam under lipopolysaccharide (LPS) challenge was characterized using RNA sequencing. The transcriptomes of LPS challenged group of clams (LH1, LH2 and LH3), and the PBS control group (CH1, CH2 and CH3), were sequenced with the Illumina HiSeq platform. Compared with the unigene expression profile of the control group, 223 unigenes were up-regulated and 389 unigenes were down-regulated in the LPS challenged group. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that signal transduction, defense response, and immune-related pathways such as Chemokine signaling pathway, Complement and coagulation cascades, NOD-like receptor signaling pathway, and Inflammatory mediator regulation of TRP channels in sensory system were the most highly enriched pathways among the genes that were differentially expressed under LPS challenge. This study present understanding of the molecular basis underpinning response to LPS challenge and provides useful information for future work on the molecular mechanism of pathogen resistance and immunity in Manila clam.

Toll-like signaling pathway in the transcriptome of Littorina littorea

The transcriptome of the caenogastropod mollusk Littorina littorea was scanned for the presence of sequences encoding Toll-like receptors (TLRs) and corresponding proteins involved in downstream TLR signaling pathway. In the transcriptomic snapshots of hemocytes and kidney tissues, 45 complete TLRs encoded by 35 genes were identified. Out of the 59 non-TLR molecules involved in a canonical TLR signaling pathway, 35 genes were classified as homologous and could be placed within the TLR-mediated MyD88-and MAPK-dependent circuitries. No reference vertebrate adapters TIRAP, TRIF and TRAM were identified in the transcriptome. The results of RNA-seq experiments with an immune challenge (rediae of the digenean Himasthla elongata) indicate that four TLRs (LlTLR1, 3, 5 and 8) and a set of upregulated genes involved in signal transduction (LlMyd88, LlTNFα, LlCASP8, LlFADD, LlNFKBIA (IkBα), LlIRAK1, LlSTAT1, LlMAPK14 (P38), LlMAP2K1 (MEK1/2), LlIRF3 and LlIRF5) may participate in the anti-digenean immune response of L. littorea.

Transcriptomic resources for evolutionary studies in flat periwinkles and related species

The flat periwinkles, Littorina fabalis and L. obtusata, comprise two sister gastropod species that have an enormous potential to elucidate the mechanisms involved in ecological speciation in the marine realm. However, the molecular resources currently available for these species are still scarce. In order to circumvent this limitation, we used RNA-seq data to characterize the transcriptome of four individuals from each species sampled in different locations across the Iberian Peninsula. Four de novo transcriptome assemblies were generated, as well as a pseudo-reference using the L. saxatilis reference transcriptome as backbone. After transcripts' annotation, variant calling resulted in the identification of 19,072 to 45,340 putatively species-diagnostic SNPs. The discriminatory power of a subset of these SNPs was validated by implementing an independent genotyping assay to characterize reference populations, resulting in an accurate classification of individuals into each species and in the identification of hybrids between the two. These data comprise valuable genomic resources for a wide range of evolutionary and conservation studies in flat periwinkles and related taxa.

Biomphalaria glabrata immunity: Post-genome advances

The freshwater snail, Biomphalaria glabrata, is an important intermediate host in the life cycle for the human parasite Schistosoma mansoni, the causative agent of schistosomiasis. Current treatment and prevention strategies have not led to a significant decrease in disease transmission. However, the genome of B. glabrata was recently sequenced to provide additional resources to further our understanding of snail biology. This review presents an overview of recently published, post-genome studies related to the topic of snail immunity. Many of these reports expand on findings originated from the genome characterization. These novel studies include a complementary gene linkage map, analysis of the genome of the B. glabrata embryonic (Bge) cell line, as well as transcriptomic and proteomic studies looking at snail-parasite interactions and innate immune memory responses towards schistosomes. Also included are biochemical investigations on snail pheromones, neuropeptides, and attractants, as well as studies investigating the frontiers of molluscan epigenetics and cell signaling were also included. Findings support the current hypotheses on snail-parasite strain compatibility, and that snail host resistance to schistosome infection is dependent not only on genetics and expression, but on the ability to form multimeric molecular complexes in a timely and tissue-specific manner. The relevance of cell immunity is reinforced, while the importance of humoral factors, especially for secondary infections, is supported. Overall, these studies reflect an improved understanding on the diversity, specificity, and complexity of molluscan immune systems.

A New Assessment of Thioester-Containing Proteins Diversity of the Freshwater Snail Biomphalaria glabrata

Thioester-containing proteins (TEPs) superfamily is known to play important innate immune functions in a wide range of animal phyla. TEPs are involved in recognition, and in the direct or mediated killing of several invading organisms or pathogens. While several TEPs have been identified in many invertebrates, only one TEP (named BgTEP) has been previously characterized in the freshwater snail, Biomphalaria glabrata. As the presence of a single member of that family is particularly intriguing, transcriptomic data and the recently published genome were used to explore the presence of other BgTEP related genes in B. glabrata. Ten other TEP members have been reported and classified into different subfamilies: Three complement-like factors (BgC3-1 to BgC3-3), one α-2-macroblobulin (BgA2M), two macroglobulin complement-related proteins (BgMCR1, BgMCR2), one CD109 (BgCD109), and three insect TEP (BgTEP2 to BgTEP4) in addition to the previously characterized BgTEP that we renamed BgTEP1. This is the first report on such a level of TEP diversity and of the presence of macroglobulin complement-related proteins (MCR) in mollusks. Gene structure analysis revealed alternative splicing in the highly variable region of three members (BgA2M, BgCD109, and BgTEP2) with a particularly unexpected diversity for BgTEP2. Finally, different gene expression profiles tend to indicate specific functions for such novel family members.

Immune response of a caenogastropod host: A case study of Littorina littorea and its digenean parasites

The caenogastropod mollusk Littorina littorea is a promising experimental model for comparative studies on host/parasite immune conflict. Several different digenean parasites use L. littorea as the first intermediate host, overcoming snail immune reactions by a wide range of tactics that are radically different among different digenean species and at different developmental parasite stages. The immune system of L. littorea is rather effective against digenean Himasthla elongata invasion, and even successfully established parasite induces a chronic host immune reaction, present at a low but stable level, that may be involved in the selection of derived parasitic clones in long lived self-sustaining infrapopulations (SSI) of rediae. An anti-digenean response in L. littorea is not systemic (non-generalized) yet tissue specific, mostly reliant on cellular rather than humoral reactions. The repertoire of immune pattern-recognizing receptors in the common periwinkle comprises diverse secreted and membrane-attached lectin molecules, as the main drivers of snail immune discrimination of digenean parasites. Comparative studies suggest that the characteristic vulnerability to digenean parasitism of L. littorea, and gastropods in general, is in part due the overall organization of immunity relative to other classes of molluscs, e.g. the immune strategy of bivalves seems to rely on less specific cellular reactions and a more generalized systemic humoral immunity. This difference may arise from the molecular features of the selective retention of their taxon-specific complement-like molecular complexes, which diverged in common ancestors of Bivalvia and Gastropoda.

Transcriptome-wide identification and characterization of the Sox gene family and microsatellites for Corbicula fluminea

The Asian clam, Corbicula fluminea, is a commonly consumed small freshwater bivalve in East Asia. However, available genetic information of this clam is still limited. In this study, the transcriptome of female C. fluminea was sequenced using the Illumina HiSeq 2500 platform. A total of 89,563 unigenes were assembled with an average length of 859 bp, and 36.7% of them were successfully annotated. Six members of Sox gene family namely SoxB1, SoxB2, SoxC, SoxD, SoxE and SoxF were identified. Based on these genes, the divergence time of C. fluminea was estimated to be around  476 million years ago. Furthermore, a total of 3,117 microsatellites were detected with a distribution density of 1:12,960 bp. Fifty of these microsatellites were randomly selected for validation, and 45 of them were successfully amplified with 31 polymorphic ones. The data obtained in this study will provide useful information for future genetic and genomic studies in C. fluminea.

Extensive Tandem Duplication Events Drive the Expansion of the C1q-Domain-Containing Gene Family in Bivalves

C1q-domain-containing (C1qDC) proteins are rapidly emerging as key players in the innate immune response of bivalve mollusks. Growing experimental evidence suggests that these highly abundant secretory proteins are involved in the recognition of microbe-associated molecular patterns, serving as lectin-like molecules in the bivalve proto-complement system. While a large amount of functional data concerning the binding specificity of the globular head C1q domain and on the regulation of these molecules in response to infection are quickly accumulating, the genetic mechanisms that have led to the extraordinary lineage-specific expansion of the C1qDC gene family in bivalves are still largely unknown. The analysis of the chromosome-scale genome assembly of the Eastern oyster Crassostrea virginica revealed that the 476 oyster C1qDC genes, far from being uniformly distributed along the genome, are located in large clusters of tandemly duplicated paralogs, mostly found on chromosomes 7 and 8. Our observations point out that the evolutionary process behind the development of a large arsenal of C1qDC lectin-like molecules in marine bivalves is still ongoing and likely based on an unequal crossing over.

A C1qDC (CgC1qDC-6) with a collagen-like domain mediates hemocyte phagocytosis and migration in oysters

Most of the bivalve C1q domain containing proteins (C1qDCs) are either only composed of the globular head domain, or contain an N-terminal coiled-coil domain, presumed to cover a role in oligomerization. On the other hand, collagen regions, widespread in vertebrate C1qDCs, are very uncommon in bivalves. In the present study, a C1qDC with a collagen-like domain (designated CgC1qDC-6) was identified from the Pacific oyster Crassostrea gigas and its possible involvement in immune responses was also characterized. The coding sequence of CgC1qDC-6 was of 756 bp, encoding a peptide of 251 amino acids with an N-terminal signal peptide, a central collagen-like domain, and a C-terminal ghC1q domain. CgC1qDC-6 was clustered with the C1qDCs from several mollusks in the phylogenetic tree. CgC1qDC-6 was detected at both mRNA and protein levels in all tested tissues including hepatopancreas, gonad, gill, mantle, adductor muscle, and hemocytes. The recombinant CgC1qDC-6 protein (rCgC1qDC-6) exhibited binding activity to various pathogen-associated molecular patterns (PAMPs) including LPS, PGN, mannose and Poly I:C, and microorganisms including Gram-negative bacteria (Escherichia coli and Vibrio splendidus), Gram-positive bacteria (Micrococcus luteus and Staphylococcus aureus), and fungus (Pichia pastoris). The phagocytic rates of oyster hemocytes towards V. splendidus pre-incubation with rCgC1qDC-6 were significantly enhanced (p < 0.05). In the chemotaxis assay, rCgC1qDC-6 could mediate the migration of oyster hemocytes in a dose-dependent manner, which exhibited a positive chemotactic effect at low concentration (<10 nM). These results collectively indicated that CgC1qDC-6 could serve as a pattern recognition receptor and mediate the hemocyte phagocytosis and migration to eliminate the invading pathogens.

Transcriptome analysis of response mechanism to ammonia stress in Asian clam (Corbicula fluminea)

Corbicula fluminea is highly sensitive to ammonia, and its response mechanism to ammonia stress is unclear. In this study, C. fluminea was exposed to different levels of ammonia (control group, 10 mg/L, and 25 mg/L) for 24 h and 48 h. A comparative analysis of transcriptome sequencing (RNA-seq) of C. fluminea digestive gland showed that the expression of 6742 genes (11.54%) was significantly affected by ammonia stress. The TLR, NF-κB, FOXO, and apoptotic signaling pathways were involved in the regulation. The differential expression of 14 genes was confirmed by real-time PCR. In summary, the response mechanism of C. fluminea digestive gland under ammonia stress may be different from that of oxidative stress in marine vertebrates. Also, the NMDAR-mediated pathway may not be the main mechanism in the response to ammonia stress in C. fluminea. The present study is a preliminary study for further investigation into ammonia toxicity in shellfish.

Derivatives of the lectin complement pathway in Lophotrochozoa

A plethora of non-overlapping immune molecular mechanisms in metazoans is the most puzzling issue in comparative immunobiology. No valid evolutionary retrospective on these mechanisms has been developed. In this study, we aimed to reveal the origin and evolution of the immune complement-like system in Lophotrochozoa. For this, we analyzed publicly available transcriptomes of prebilaterian and lophotrochozoan species, mapping lineage-specific molecular events on the phylogenetic tree. We found that there were no orthologs of mannose-binding lectin (MBL) and ficolins (FCN) in Lophotrochozoa but C1q-like proteins (C1qL), bearing both a collagen domain and a globular C1q domain, were omnipresent in them. This suggests that among all complement-like activators the C1qL-specific domain architecture was an evolutionarily first. Two novel protostomian MASP-Related Molecules, MReM1 and MReM2, might hypothetically compensate for the loss of a prebilaterian MASP-orthologous gene and act in complex with C1qL and C1qDC as a "proto-activator" of an ancient "proto-complement". We proposed a new model of the complement evolution predicting that numerous lineage-specific complement-like systems should have evolved from a stem "antique" molecular complex. First evolved in the common ancestor of coelomic animals, the "antique" humoral complex consisted of a TEP molecule, the common ancestor of TEP-associated proteases (C2/Bf/Сf/Lf), the common ancestor of MASP-like proteases (MASP/C1r/C1s, MReM1/MReM2) and multimeric recognition proteins (C1q-, MBL- and FCN-homologs). Further evolutionary specialization and expansion of the complex was independent and lineage-specific, examples being the mammalian complement system and the Apogastropoda complement-like complex. The latter includes an impressive array of multimeric recognition proteins, the variable immunoglobulin and lectin domain containing molecules (VIgL), homologous to C1q, MBL, FCN and other lectins. Four novel polymorphic subfamilies of VIgLs were found to be expressed in Apogastropoda: C1q-related proteins (QREP), zona pellucida-related proteins (ZREP), Scavenger Receptor Cys-Rich-related proteins (SREP) and HPA-lectin related proteins (HREP). The transcriptional response of fibrinogen-related proteins of VIgL family (LlFREP), LlQREP and LlSREP to infestation of common periwinkle, Littorina littorea, with digenean parasite Himasthla elongata correlates with that of LlMReM1, supporting the model suggested in this study.